Image difference.pro
From RHESSI Wiki
(Difference between revisions)
Agopie (Talk | contribs)
(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-> ...)
Newer edit →
(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-> ...)
Newer edit →
Revision as of 17:39, 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