obj = hsi_image() obj-> set, det_index_mask= [0, 1, 1, 1, 1, 1, 1, 0, 0] obj-> set, im_time_interval= '6-Jan-2004 '+['06:22:20', '06:23:00'] obj->set, im_energy_binning=[25.,35.] obj-> set, image_algorithm= 'VIS FWDFIT' obj-> set, xyoffset=[-970,73] obj-> set, use_phz_stacker= 1 ;obj-> set,PHZ_radius=10. obj-> set,image_dim=[64,64] obj-> set,pixel_size=[1.0,1.0] obj-> set, vf_multi=1 obj-> set,phz_n_roll_bins_control=[20,38,20,20,12,12,12,12,20] ;sets the number of roll bins per detector ; useful controls when calculating the visibilities obj-> set, vis_edit=1, vis_combine=1 ; remove outliers and combine conjugates our_fit=obj ->get(/vf_srcout) ;gets visibilities parameters fit_err=obj ->get(/vf_sigma) ;gets visibilities parameter errors ; you can set some initial parameters for you visibility structure our_fit.srctype='ellipse' our_fit.srcflux=[11,9] our_fit.srcfwhm=[7,6] our_fit.eccen=[0.7,0.7] our_fit.srcpa=[175,-178] our_fit.srcx=[907,904] our_fit.srcy=[248,272] obj -> set, vf_srcin=our_fit ; sets out fit params obj-> plot obj->plotman