Blind Image Analysis
From RHESSI Wiki
Images Reconstructed by Brian Dennis
Images were reconstructed using Clean and pixon for the first four test cases using the following instructions from Richard Schwartz:
Here is the link to the download with the files needed to run your blind image reconstructions http://hesperia.gsfc.nasa.gov/~richard/blindtest.zip Download it and unpack it into a directory. There are a few minor changes in files that will go online next week. I didn't want to chance breaking everything over the weekend by putting this all in SSW.
Get into sswidl in the directory where you unpacked these files. Use an empty one or not as you please.
From the command line
This builds your image object, obj, and loads in the simulated calib_eventlist (stacked) of your choice. The file is set to use image 0 right now but if you look in the file you can see how to load in the other three as you choose. Once you use the oc->newdata method you can't go back to acquiring data normally in that image object without doing this oc->newdata, /clear
Note that this is the calib_eventlist object, OC, one of the constituents of the image object, OBJ
Also, using the seed causes the poisson deviates to be used for the simulation. If you set seed to a fixed value, use a number gt 100, then you'll get the same values every time. Calling it without using seed will send in the expected values for the counts without any poisson simulation.
You can use this OBJ for all of the images without regeneration, and for any and all analysis routines.
To use Image I
These next two lines will cause the current image analysis routine to use the new data set
obj-> set, det_index_mask= [0B, 1B, 1B, 1B, 1B, 1B, 1B, 1B, 1B] obj-> set, det_index_mask= [1B, 1B, 1B, 1B, 1B, 1B, 1B, 1B, 1B]
I will make some changes next week so send in your requests. I will also make some new files of calib_eventlists from other images of my choosing. These were designed by John Brown. That is a big hint.
Richard Schwartz firstname.lastname@example.org
Test Image #1 (i = 0)
These pixon and Clean images show that there are two main sources with the north-east source appearing loop shaped. The pixon image shows an extended elliptical source surrounding the two main compact sources. The Pixon Modulation Profiles show that there are acceptable values of the C-statistic for all detectors except possibly for detector #8 where C = 2.06.
Test Image #2 (i = 1)
This test image appears as a line source, perhaps a series of up to four compact sources, with a separate weaker source to the south-west. Pixon shows an additional extended almost circular source surrounding these compact sources. The Pixon Modulation Profiles show surprisingly large values for the C-statistic with numbers over 4.0 for all detectors except detector #7 where C = 3.45.
Test Image #3 (i = 2)
This Clean image showing just the Clean components suggests that there is significant fine structure in this source but the details probably can't be trusted. The details do not shown up with Pixon, partly because Pixon cannot be used with pixels smaller than 1 arcsecond square.