Image difference.pro

From RHESSI Wiki

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