function sxt2eit, sxtfits, eitfits ;+ ; Name: sxt2eit ; ; Purpose: align SXT image to EIT ; ; Calling Sequence: ; sxtaligned=sxt2eit(sxtfits, eitfits) ; ; History: proto-version, 23-jan-1996 ; ; Restrictions: ; only full sun EIT and SXT , simple FITS for now ;- if not data_chk(eitfits,/string) or not data_chk(sxtfits,/string) then $ message,"Must supply EIT and SXT FITs file names if not file_exist(eitfits) then message,"can't find: " + eitfits if not file_exist(sxtfits) then message,"can't find: " + sxtfits sxtimg=readfits(sxtfits,shead) ; read SXT image ; ------------- get EIT (SOHO) pointing info ---------------------- ; att_info=get_sc_att(etime) ; not yet.. eit_getlimb, eitfits, ex, ey, er, image=eitimg, header=ehead ;---- todo - fits header summary - extract time, size, pointing info --- stime=fmt_tim(fxpar(shead,'date-obs') + ' ' + fxpar(shead,'time-obs')) etime=fmt_tim(fxpar(ehead,'date-obs') + ' ' + fxpar(ehead,'time-obs')) ; assume SXT sun-center in FITS file (otherwise, call sxt_cen) spoint=[fxpar(shead,'crpix1'),fxpar(shead,'crpix2'),fxpar(shead,'crota')] ssize =[fxpar(shead,'naxis1'),fxpar(shead,'naxis2')] spoint=spoint * (1024./ssize(0)) esize =[fxpar(ehead,'naxis1'),fxpar(ehead,'naxis2')] ;*********** ; kludges for full EIT only if esize(0) ne 1024 or esize(1) ne 1024 then begin message,/info,"Sorry - only full frame EIT images for now... return,eitimg endif sxtimg=congrid(sxtimg,1024,1024) ; zoom to EIT ssize=ssize*(1024./ssize(0)) ;*********** ; -------------------------------------------------------------------- simg=coal_image(sxtimg,[ex, ey, er],spoint,ssize) ; translate SXT simg=rot(simg,-(fxpar(shead,'CROTA')),gt_pix_size()/eit_pixsize()) ; roll and demagnify return, simg end
Last revised: - Wed May 9 21:45:23 2007- F. Auchère