; This is a simple script which computes DEMs from a set of AIA image files. ; After it is finished (which may take a while depending on the resolution of the images), ; The DEMs are stored in memory in a structure called dem_struc (see readme and the comments ; for firdem.pro for a description of the structure's contents and how to work with it). .compile ../estimate_quantile.pro .compile ../firdem.pro .compile aia_firdem_wrapper.pro ; This search string should be such that findfile finds the set of 6 AIA images to invert: if(n_elements(search_string) eq 0) then search_string = "../../aia_test_images/aia*.fits" filenames=findfile(search_string) indices_used = [0,1,2,3,4,5] aia_prep,filenames,indices_used,indout,datout,/uncomp_delete print,filenames(indices_used) channels='A'+string(format='(%"%d")',indout.wavelnth) ; Computing DEMs for a full-resolution AIA image takes a while (~ 1 hour). For testing, ; we rebin the data to 1024x1024. Comment out this line for full resolution: datout=rebin(datout,1024,1024,n_elements(channels)) aia_firdem_wrapper,channels,indout.exptime,datout,dem_struc