Blind Image Analysis

From RHESSI Wiki

(Difference between revisions)
Jump to: navigation, search
(Test Image #3 (i = 2))
(Test Image #3 (i = 2))
Line 44: Line 44:
== Test Image #3 (i = 2)==
== Test Image #3 (i = 2)==
-
[[Image:im_ti2_cleanuniformBWF10_pixon.jpg|border|right|text-top|500px|Image shows a Clean image in color with the pixon image overlaid as white contours. The Clean image was made with uniform weighting,a gain of 0.01, 1000 iterations, and a beam width factor of 10. The latter factor allows the individual Clean components to be shown. The Pixon image was made with 1 arcsecond pixels using all nine detectors.]]
+
[[Image:im_ti2_cleanuniformBWF10_pixon.jpg|border|right|text-top|500px|Caption]]
 +
 
 +
Image shows a Clean image in color with the pixon image overlaid as white contours. The Clean image was made with uniform weighting,a gain of 0.01, 1000 iterations, and a beam width factor of 10. The latter factor allows the individual Clean components to be shown. The Pixon image was made with 1 arcsecond pixels using all nine detectors.
----
----

Revision as of 22:02, 22 October 2009

Contents

Images Reconstructed by Brian Dennis

Images were reconstructed for the first four test cases using the following instructions from Richard Schwartz:


Instructions

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

.run blind_image_analysis_prep

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

oc->setnewdata, pcbe[*,I],seed=seed

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 richard.schwartz@nasa.gov

Test Image #1 (i = 0)

Test Image #2 (i = 1)

Test Image #3 (i = 2)

Caption

Image shows a Clean image in color with the pixon image overlaid as white contours. The Clean image was made with uniform weighting,a gain of 0.01, 1000 iterations, and a beam width factor of 10. The latter factor allows the individual Clean components to be shown. The Pixon image was made with 1 arcsecond pixels using all nine detectors.


Test Image #4 (i = 3)

Im ti3 pixon cleanun iformBWF2.jpg

Pixon image with Clean image overlaid as white contours. Clean image was made with uniform weighting, a gain of 0.05, 500 iterations,a Clean Beam Width Factor of 2.0, and 64x64 0.5 arcsecond pixels.

Personal tools
Namespaces
Variants
Actions
Navigation
Toolbox