Commits

Takahiro Tsutsumi authored 4df743585f0
Added to collect the beam stats for serial

casatests/stakeholder/test_stk_alma_pipeline_imaging.py

Modified
943 943 deconvolver='hogbom', restoration=True, pbcor=True, \
944 944 weighting='briggsbwtaper', robust=0.5, npixels=0, niter=20000, \
945 945 threshold='0.354Jy', interactive=0, usemask='auto'
946 946 '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
947 947 lownoisethreshold=2.0, negativethreshold=0.0, \
948 948 minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
949 949 minpercentchange=1.0, fastnoise=False, restart=True, \
950 950 calcres=False, calcpsf=False, savemodel='none', \
951 951 parallel=True, verbose=True)
952 952
953 - # retrieve per-channel beam statistics (only in parallel)
953 + # retrieve per-channel beam statistics
954 954 bmin_dict, bmaj_dict, pa_dict = \
955 955 self.cube_beam_stats(image=img+'.image')
956 956
957 957 tclean(vis=self.msfile, imagename=file_name+'1', spw=['0'], \
958 958 field='1', imsize=[80, 80], cell=['1.1arcsec'], \
959 959 phasecenter='ICRS 00:45:54.3836 -073.15.29.413', stokes='I', \
960 960 antenna=['0,1,2,3,4,5,6,7,8'], scan=['8,12,16'], \
961 961 intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
962 962 specmode='cube', nchan=508, start='220.2526743594GHz', \
963 963 width='0.2441741MHz', outframe='LSRK', \
987 987 deconvolver='hogbom', restoration=True, pbcor=True, \
988 988 weighting='briggsbwtaper', robust=0.5, npixels=0, niter=20000, \
989 989 threshold='0.354Jy', interactive=0, usemask='auto'
990 990 '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
991 991 lownoisethreshold=2.0, negativethreshold=0.0, \
992 992 minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
993 993 minpercentchange=1.0, fastnoise=False, restart=True, \
994 994 calcres=False, calcpsf=False, savemodel='none', \
995 995 restoringbeam='common', parallel=False, verbose=True)
996 996
997 + # retrieve per-channel beam statistics
998 + bmin_dict, bmaj_dict, pa_dict = \
999 + self.cube_beam_stats(image=img+'.psf')
1000 +
997 1001
998 1002 report0 = th.checkall(imgexist = self.image_list(img, 'standard'))
999 1003
1000 1004 # .image report(test_standard_cube_briggsbwtaper)
1001 1005 im_stats_dict = self.image_stats(image=img+'.image', fit_region = \
1002 1006 'ellipse[[11.47881897deg, -73.25881015deg], [9.0414arcsec, 8.4854arcsec], 90.00000000deg]')
1003 1007
1004 1008 # test_standard_cube_briggsbwtaper.exp_im_stats
1005 1009 exp_im_stats = {'com_bmaj': [False, 8.900085299026527],
1006 1010 'com_bmin': [False, 6.193896355268416],
1239 1243 # im_stats, mask_stats, pb_stats, psf_stats,\
1240 1244 # model_stats, resid_stats, sumwt_stats]
1241 1245 savedict['im_stats_dict']=im_stats_dict
1242 1246 savedict['mask_stats_dict']=mask_stats_mod_dict
1243 1247 savedict['pb_stats_dict']=pb_stats_dict
1244 1248 savedict['psf_stats_dict']=psf_stats_dict
1245 1249 savedict['model_stats_dict']=model_stats_dict
1246 1250 savedict['resid_stats_dict']=resid_stats_dict
1247 1251 savedict['sumwt_stats_dict']=sumwt_stats_dict
1248 1252
1249 - if self.parallel:
1250 - savedict['bmin_dict']=bmin_dict
1251 - savedict['bmaj_dict']=bmaj_dict
1252 - savedict['pa_dict']=pa_dict
1253 + #if self.parallel:
1254 + # for serial get it from *.psf
1255 + savedict['bmin_dict']=bmin_dict
1256 + savedict['bmaj_dict']=bmaj_dict
1257 + savedict['pa_dict']=pa_dict
1253 1258
1254 1259 self.save_dict_to_file(test_name, savedict, test_name+'_cas13317mod_stats')
1255 1260
1256 1261
1257 1262 self.assertTrue(th.check_final(pstr = report), \
1258 1263 msg = failed)
1259 1264
1260 1265 # End of test_standard_cube_briggsbwtaper
1261 1266 #-------------------------------------------------#
1262 1267
4072 4077 restoration=True, pbcor=True, weighting='briggsbwtaper', robust=0.5,\
4073 4078 npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4074 4079 interactive=0, usemask='auto-multithresh', \
4075 4080 sidelobethreshold=1.25, noisethreshold=5.0, \
4076 4081 lownoisethreshold=2.0, negativethreshold=0.0, \
4077 4082 minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4078 4083 minpercentchange=1.0, fastnoise=False, restart=True, \
4079 4084 savemodel='none', calcres=False, calcpsf=False, \
4080 4085 restoringbeam='common', parallel=False, verbose=True)
4081 4086
4087 + # retrieve per-channel beam statistics
4088 + bmin_dict, bmaj_dict, pa_dict = \
4089 + self.cube_beam_stats(img+'.psf')
4090 +
4082 4091 report0 = th.checkall(imgexist = self.image_list(img, 'mosaic'))
4083 4092
4084 4093 # .image report (test_mosaic_cube_briggsbwtaper)
4085 4094 im_stats_dict = self.image_stats(img+'.image', fit_region = \
4086 4095 'ellipse[[11.48661818deg, -73.26292371deg], [8.2211arcsec, 7.4698arcsec], 90.00000000deg]', field_regions = \
4087 4096 ['circle[[00:45:54.383559, -73.15.29.41306], 22.45arcsec]',
4088 4097 'circle[[00:45:49.435664, -73.15.35.13742], 22.45arcsec]',
4089 4098 'circle[[00:45:53.057440, -73.15.50.79016], 22.45arcsec]',
4090 4099 'circle[[00:45:50.762696, -73.15.13.76177], 22.45arcsec]',
4091 4100 'circle[[00:45:58.006248, -73.15.45.06040], 22.45arcsec]',
4364 4373 # model_stats, resid_stats, sumwt_stats, + wt_stats (mosaic)
4365 4374 savedict['im_stats_dict']=im_stats_dict
4366 4375 savedict['mask_stats_dict']=mask_stats_mod_dict
4367 4376 savedict['pb_stats_dict']=pb_stats_mod_dict
4368 4377 savedict['psf_stats_dict']=psf_stats_dict
4369 4378 savedict['model_stats_dict']=model_stats_dict
4370 4379 savedict['resid_stats_dict']=resid_stats_dict
4371 4380 savedict['sumwt_stats_dict']=sumwt_stats_dict
4372 4381 savedict['wt_stats_dict']=wt_stats_dict
4373 4382
4374 - if self.parallel:
4375 - savedict['bmin_dict']=bmin_dict
4376 - savedict['bmaj_dict']=bmaj_dict
4377 - savedict['pa_dict']=pa_dict
4383 + #if self.parallel:
4384 + savedict['bmin_dict']=bmin_dict
4385 + savedict['bmaj_dict']=bmaj_dict
4386 + savedict['pa_dict']=pa_dict
4378 4387
4379 4388 self.save_dict_to_file(test_name,savedict, test_name+'_cas13317_stats')
4380 4389
4381 4390 self.assertTrue(th.check_final(pstr = report), \
4382 4391 msg = failed)
4383 4392
4384 4393 # End of test_mosaic_cube_briggsbwtaper
4385 4394 #-------------------------------------------------#
4386 4395 @stats_dict(test_dict)
4387 4396 # @unittest.skip("")
5784 5793 # model_stats, resid_stats, sumwt_stats, + wt_stats (mosaic)
5785 5794 savedict['im_stats_dict']=im_stats_dict
5786 5795 savedict['mask_stats_dict']=mask_stats_mod_dict
5787 5796 savedict['pb_stats_dict']=pb_stats_dict
5788 5797 savedict['psf_stats_dict']=psf_stats_dict
5789 5798 savedict['model_stats_dict']=model_stats_dict
5790 5799 savedict['resid_stats_dict']=resid_stats_dict
5791 5800 savedict['sumwt_stats_dict']=sumwt_stats_dict
5792 5801 savedict['wt_stats_dict']=wt_stats_dict
5793 5802
5794 - #if self.parallel:
5795 - # savedict['bmin_dict']=bmin_dict
5796 - # savedict['bmaj_dict']=bmaj_dict
5797 - # savedict['pa_dict']=pa_dict
5798 -
5799 5803 self.save_dict_to_file(test_name,savedict, test_name+'_cas13317_stats')
5800 5804
5801 5805 self.assertTrue(th.check_final(pstr = report), \
5802 5806 msg = failed)
5803 5807
5804 5808 # End of test_mosaic_cube_eph_briggsbwtaper
5805 5809 #-------------------------------------------------#
5806 5810 @stats_dict(test_dict)
5807 5811 # @unittest.skip("")
5808 5812 def test_mosaic_mfs_eph(self):

Everything looks good. We'll let you know here if there's anything you should know about.

Add shortcut