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