Image difference.pro
From RHESSI Wiki
(Difference between revisions)
(New page: pro image_difference filename = 'hsi_evlist_025_050keV.fits' img_obj = hsi_image(ev_filename = filename) img_obj->set,im_time_interval='2-Sept-2002 '+[' 00:07:36', '00:07:56'] img_obj-> ...) |
|||
Line 2: | Line 2: | ||
filename = 'hsi_evlist_025_050keV.fits' | filename = 'hsi_evlist_025_050keV.fits' | ||
+ | |||
img_obj = hsi_image(ev_filename = filename) | img_obj = hsi_image(ev_filename = filename) | ||
+ | |||
img_obj->set,im_time_interval='2-Sept-2002 '+[' 00:07:36', '00:07:56'] | img_obj->set,im_time_interval='2-Sept-2002 '+[' 00:07:36', '00:07:56'] | ||
img_obj-> set, det_index_mask=[1,1,1,1,1,1,1,1,1] | img_obj-> set, det_index_mask=[1,1,1,1,1,1,1,1,1] | ||
+ | |||
img_obj-> set, image_dim= [64, 64] | img_obj-> set, image_dim= [64, 64] | ||
+ | |||
img_obj-> set, pixel_size= [.5, .5] | img_obj-> set, pixel_size= [.5, .5] | ||
+ | |||
img_obj-> set, xyoffset = [0.0, 0.0] | img_obj-> set, xyoffset = [0.0, 0.0] | ||
img_obj-> set, image_algorithm= 'Clean' | img_obj-> set, image_algorithm= 'Clean' | ||
+ | |||
img_obj-> set, clean_niter = 500 | img_obj-> set, clean_niter = 500 | ||
+ | |||
img_obj-> set, clean_beam_width_factor = 2. | img_obj-> set, clean_beam_width_factor = 2. | ||
+ | |||
image = img_obj -> getdata() | image = img_obj -> getdata() | ||
Line 20: | Line 28: | ||
point_map = make_map(p.maps) | point_map = make_map(p.maps) | ||
+ | |||
point_map.dx = .25 | point_map.dx = .25 | ||
+ | |||
point_map.dy = .25 | point_map.dy = .25 | ||
+ | |||
point_src = point_map.data | point_src = point_map.data | ||
+ | |||
point_src = congrid(point_src,64,64) | point_src = congrid(point_src,64,64) | ||
dif_map = image - point_src | dif_map = image - point_src | ||
- | |||
- | |||
end | end |
Latest revision as of 17:41, 18 August 2010
pro image_difference
filename = 'hsi_evlist_025_050keV.fits'
img_obj = hsi_image(ev_filename = filename)
img_obj->set,im_time_interval='2-Sept-2002 '+[' 00:07:36', '00:07:56']
img_obj-> set, det_index_mask=[1,1,1,1,1,1,1,1,1]
img_obj-> set, image_dim= [64, 64]
img_obj-> set, pixel_size= [.5, .5]
img_obj-> set, xyoffset = [0.0, 0.0]
img_obj-> set, image_algorithm= 'Clean'
img_obj-> set, clean_niter = 500
img_obj-> set, clean_beam_width_factor = 2.
image = img_obj -> getdata()
img_obj -> plotman, plotman = plotman
restore, 'mc_maps_point.sav ', /verbose
point_map = make_map(p.maps)
point_map.dx = .25
point_map.dy = .25
point_src = point_map.data
point_src = congrid(point_src,64,64)
dif_map = image - point_src
end