pro eit_sub_movie,header,movie,tstart,tstop,ref_helio=ref_helio,ref_date=ref_date,$ outsize=outsize,wave=wave,lz=lz,object=object,arcsec=arcsec,list=list,$ phot_norm=phot_norm ;+ ; NAME: ; EIT_SUB_MOVIE ; ; PURPOSE: ; Create a movie of a subfield image from a sequence of full ; field iages, tracking the region of interest with solar ; rotation. ; ; CATEGORY: ; FITS processing ; ; CALLING SEQUENCE: ; EIT_SUB_MOVIE outheader,outmovie,start time,stop time [,/arcsec] ; [,ref_helio=ref_helio] [,ref_date=ref_date] [,wave=wave] ; [,outsize=outsize] [,object=object] [,list=list] ; ; INPUTS: ; tstart I = start time ; tstop I = stop time ; ; OPTIONAL INPUT KEYWORD PARAMETERS: ; ref_helio I = heliocentric Center of FOV (lon,lat) ; arcsec I = set if ref_helio is in arcsec coords (lon-lat) instead of ; helio lon-lat ; ref_date I = reference time/date of center ; outsize I = number of pixels for output image (x,y) ; wave I = wavelength of images (default=195) ; object I = object type (e.g. 'full fov') default = all images ; list I = input list instead of time range ; phot_norm I = set for using phot_norm keyword in eit_prep ; ; OUTPUTS: ; header O = output header (structure) ; movie O = output processed subfield 3D data array ; ; COMMON BLOCKS: none. ; ; SIDE EFFECTS: ; ; RESTRICTIONS: ; ; PROCEDURE: ; Read in catalog data. Read in selected images and headers. ; Force solar pointing over short period of images to be constant, ; this eliminates the "bad" limb fits. Extract selected FOV, tracking ; solar rotation. ; ; MODIFICATION HISTORY: ; Written by: J. Newmark Date. September 1997 ; 24-Nov-1997 J. Newmark - Account for CRPIX1,2 offset in FITS header ; 07-Apr-1998 J. Newmark - update ; 30-jul-1999 J. Newmark - add phot_norm keyword ;- if not keyword_set(wave) then wave = 195 if not keyword_set(object) then object = 0 if n_elements(list) eq 0 then $ list = eit_catrd(timerange=[tstart,tstop],wave=wave,object=object) read_eit,eit_file2path(list),idx1,data1,outsize=[1024,1024] ; ; use single image x,y center since limb fits vary ; subtract 1 since going from FITS to IDL. Add back in later. ; crpix1 = median(idx1.crpix1 - 1) idx1.crpix1 = crpix1 crpix2 = median(idx1.crpix2 - 1) idx1.crpix2 = crpix2 ; ; extract chosen field of view ssw_track_fov,idx1,temporary(data1),idx2,data2,ref_helio=ref_helio,$ ref_date=ref_date,outsize=outsize,arcsec=arcsec sz = size(data2) for i = 1,sz(3)-1 do fake_missing_blocks,data2,i,/before,block=16 fake_missing_blocks,data2,0,/after,block=16 eit_prep,idx2,header,movie,data=temporary(data2),/normalize,phot_norm=phot_norm idx2.crpix1 = idx2.crpix1 + 1 idx2.crpix2 = idx2.crpix2 + 1 end
Last revised: - Wed May 9 21:45:03 2007- F. Auchère