Commits

Takahiro Tsutsumi authored e6046b15508 
Removed the commented tclean commands
f1bf919f2a6
master

casa5/gcwrap/python/scripts/tests/test_stk_alma_pipeline_imaging.py

Modified
596 596 gridder='standard', mosweight=False, \
597 597 deconvolver='hogbom', restoringbeam='common', restoration=True, pbcor=True, \
598 598 weighting='briggs', robust=0.5, npixels=0, niter=20000, \
599 599 threshold='0.354Jy', interactive=0, usemask='auto'
600 600 '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
601 601 lownoisethreshold=2.0, negativethreshold=0.0, \
602 602 minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
603 603 minpercentchange=1.0, fastnoise=False, restart=True, \
604 604 calcres=False, calcpsf=False, savemodel='none', \
605 605 parallel=True, verbose=True)
606 -​
607 - # retrieve per-channel beam statistics (only in parallel)
608 - #bmin_dict, bmaj_dict, pa_dict = \
609 - # self.cube_beam_stats(image=img+'.image')
610 -​
611 - #tclean(vis=self.msfile, imagename=file_name+'1', spw=['0'], \
612 - # field='1', imsize=[80, 80], cell=['1.1arcsec'], \
613 - # phasecenter='ICRS 00:45:54.3836 -073.15.29.413', stokes='I', \
614 - # antenna=['0,1,2,3,4,5,6,7,8'], scan=['8,12,16'], \
615 - # intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
616 - # specmode='cube', nchan=508, start='220.2526743594GHz', \
617 - # width='0.2441741MHz', outframe='LSRK', \
618 - # perchanweightdensity=False, gridder='standard', \
619 - # mosweight=False, usepointing=False, \
620 - # pblimit=0.2, deconvolver='hogbom', restoration=True, \
621 - # restoringbeam='common', pbcor=True, weighting='briggs', \
622 - # robust=0.5, npixels=0, niter=0, threshold='0.354Jy', \
623 - # nsigma=0.0, interactive=0, usemask='auto-multithresh', \
624 - # sidelobethreshold=1.25, noisethreshold=5.0, \
625 - # lownoisethreshold=2.0, negativethreshold=0.0, \
626 - # minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
627 - # minpercentchange=1.0, fastnoise=False, restart=True, \
628 - # calcres=False, calcpsf=False, savemodel='none', \
629 - # parallel=False, verbose=True)
630 -​
631 606 else:
632 607 tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
633 608 spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
634 609 scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
635 610 datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
636 611 '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
637 612 nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
638 613 outframe='LSRK', perchanweightdensity=False, \
639 614 usepointing=False, pblimit=0.2, nsigma=0.0, \
640 615 gridder='standard', mosweight=False, \
965 940 gridder='standard', mosweight=False, \
966 941 deconvolver='hogbom', restoration=True, restoringbeam='common', pbcor=True, \
967 942 weighting='briggs', robust=0.5, npixels=0, niter=20000, \
968 943 threshold='0.354Jy', interactive=0, usemask='auto'
969 944 '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
970 945 lownoisethreshold=2.0, negativethreshold=0.0, \
971 946 minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
972 947 minpercentchange=1.0, fastnoise=False, restart=True, \
973 948 calcres=False, calcpsf=False, savemodel='none', \
974 949 parallel=True, verbose=True)
975 -​
976 - # retrieve per-channel beam statistics
977 - #bmin_dict, bmaj_dict, pa_dict = \
978 - # self.cube_beam_stats(image=img+'.image')
979 -​
980 - #tclean(vis=self.msfile, imagename=file_name+'1', spw=['0'], \
981 - # field='1', imsize=[80, 80], cell=['1.1arcsec'], \
982 - # phasecenter='ICRS 00:45:54.3836 -073.15.29.413', stokes='I', \
983 - # antenna=['0,1,2,3,4,5,6,7,8'], scan=['8,12,16'], \
984 - # intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
985 - # specmode='cube', nchan=508, start='220.2526743594GHz', \
986 - # width='0.2441741MHz', outframe='LSRK', \
987 - # perchanweightdensity=True, gridder='standard', \
988 - # mosweight=False, usepointing=False, \
989 - # pblimit=0.2, deconvolver='hogbom', restoration=True, \
990 - # restoringbeam='common', pbcor=True, weighting='briggs', \
991 - # robust=0.5, npixels=0, niter=0, threshold='0.354Jy', \
992 - # nsigma=0.0, interactive=0, usemask='auto-multithresh', \
993 - # sidelobethreshold=1.25, noisethreshold=5.0, \
994 - # lownoisethreshold=2.0, negativethreshold=0.0, \
995 - # minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
996 - # minpercentchange=1.0, fastnoise=False, restart=True, \
997 - # calcres=False, calcpsf=False, savemodel='none', \
998 - # parallel=False, verbose=True)
999 -​
1000 950 else:
1001 951 tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
1002 952 spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
1003 953 scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
1004 954 datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
1005 955 '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
1006 956 nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
1007 957 outframe='LSRK', perchanweightdensity=True, \
1008 958 usepointing=False, pblimit=0.2, nsigma=0.0, \
1009 959 gridder='standard', mosweight=False, \
1337 1287 gridder='standard', mosweight=False, \
1338 1288 deconvolver='hogbom', restoration=True, restoringbeam='common', pbcor=True, \
1339 1289 weighting='briggsbwtaper', robust=0.5, npixels=0, niter=20000, \
1340 1290 threshold='0.354Jy', interactive=0, usemask='auto'
1341 1291 '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
1342 1292 lownoisethreshold=2.0, negativethreshold=0.0, \
1343 1293 minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
1344 1294 minpercentchange=1.0, fastnoise=False, restart=True, \
1345 1295 calcres=False, calcpsf=False, savemodel='none', \
1346 1296 parallel=True, verbose=True)
1347 -​
1348 - # retrieve per-channel beam statistics
1349 - #bmin_dict, bmaj_dict, pa_dict = \
1350 - # self.cube_beam_stats(image=img+'.image')
1351 -​
1352 - #tclean(vis=self.msfile, imagename=file_name+'1', spw=['0'], \
1353 - # field='1', imsize=[80, 80], cell=['1.1arcsec'], \
1354 - # phasecenter='ICRS 00:45:54.3836 -073.15.29.413', stokes='I', \
1355 - # antenna=['0,1,2,3,4,5,6,7,8'], scan=['8,12,16'], \
1356 - # intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
1357 - # specmode='cube', nchan=508, start='220.2526743594GHz', \
1358 - # width='0.2441741MHz', outframe='LSRK', \
1359 - # perchanweightdensity=True, gridder='standard', \
1360 - # mosweight=False, usepointing=False, \
1361 - # pblimit=0.2, deconvolver='hogbom', restoration=True, \
1362 - # restoringbeam='common', pbcor=True, weighting='briggsbwtaper', \
1363 - # robust=0.5, npixels=0, niter=0, threshold='0.354Jy', \
1364 - # nsigma=0.0, interactive=0, usemask='auto-multithresh', \
1365 - # sidelobethreshold=1.25, noisethreshold=5.0, \
1366 - # lownoisethreshold=2.0, negativethreshold=0.0, \
1367 - # minbeamfrac=0.08, growiterations=75, dogrowprune=True, \
1368 - # minpercentchange=1.0, fastnoise=False, restart=True, \
1369 - # calcres=False, calcpsf=False, savemodel='none', \
1370 - # parallel=False, verbose=True)
1371 -​
1372 1297 else:
1373 1298 tclean(vis=self.msfile, imagename=file_name+'1', field='1', \
1374 1299 spw=['0'], imsize=[80, 80], antenna=['0,1,2,3,4,5,6,7,8'], \
1375 1300 scan=['8,12,16'], intent='OBSERVE_TARGET#ON_SOURCE', \
1376 1301 datacolumn='data', cell=['1.1arcsec'], phasecenter='ICRS '
1377 1302 '00:45:54.3836 -073.15.29.413', stokes='I', specmode='cube', \
1378 1303 nchan=508, start='220.2526743594GHz', width='0.2441741MHz',\
1379 1304 outframe='LSRK', perchanweightdensity=True, \
1380 1305 usepointing=False, pblimit=0.2, nsigma=0.0, \
1381 1306 gridder='standard', mosweight=False, \
4390 4315 restoration=True, restoringbeam='common', pbcor=True, \
4391 4316 weighting='briggs', robust=0.5,\
4392 4317 npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4393 4318 interactive=0, usemask='auto-multithresh', \
4394 4319 sidelobethreshold=1.25, noisethreshold=5.0, \
4395 4320 lownoisethreshold=2.0, negativethreshold=0.0, \
4396 4321 minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4397 4322 minpercentchange=1.0, fastnoise=False, restart=True, \
4398 4323 savemodel='none', calcres=False, calcpsf=False, \
4399 4324 parallel=True, verbose=True)
4400 -​
4401 - # retrieve per-channel beam statistics
4402 - #bmin_dict, bmaj_dict, pa_dict = \
4403 - # self.cube_beam_stats(img+'.image')
4404 -​
4405 - #tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4406 - # antenna=['0,1,2,3,4,5,6,7,8'], scan=['8,12,16'], \
4407 - # intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4408 - # imagename=file_name+'1', imsize=[108, 108], \
4409 - # cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4410 - # ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4411 - # start='220.2526743594GHz', width='0.2441741MHz', \
4412 - # outframe='LSRK', perchanweightdensity=False, \
4413 - # gridder='mosaic', mosweight=True, \
4414 - # usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4415 - # restoration=True, restoringbeam='common', pbcor=True, \
4416 - # weighting='briggs', robust=0.5, npixels=0, niter=0, \
4417 - # threshold='0.354Jy', nsigma=0.0, interactive=0, usemask='auto'
4418 - # '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
4419 - # lownoisethreshold=2.0, negativethreshold=0.0, \
4420 - # minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4421 - # minpercentchange=1.0, fastnoise=False, restart=True, \
4422 - # savemodel='none', calcres=False, calcpsf=False,
4423 - # parallel=False, verbose=True)
4424 -​
4425 4325 else:
4426 4326 tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4427 4327 antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4428 4328 intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4429 4329 imagename=file_name+'1', imsize=[108, 108], \
4430 4330 cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4431 4331 ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4432 4332 start='220.2526743594GHz', width='0.2441741MHz', \
4433 4333 outframe='LSRK', perchanweightdensity=False, \
4434 4334 gridder='mosaic', mosweight=True, \
4755 4655 tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4756 4656 antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4757 4657 intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4758 4658 imagename=file_name+'1', imsize=[108, 108], \
4759 4659 cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4760 4660 ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4761 4661 start='220.2526743594GHz', width='0.2441741MHz', \
4762 4662 outframe='LSRK', perchanweightdensity=True, \
4763 4663 gridder='mosaic', mosweight=True, \
4764 4664 usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4765 - restoration=True, pbcor=True, weighting='briggs', robust=0.5,\
4665 + restoration=True, restoringbeam='common', pbcor=True, \
4666 + weighting='briggs', robust=0.5,\
4766 4667 npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4767 4668 interactive=0, usemask='auto-multithresh', \
4768 4669 sidelobethreshold=1.25, noisethreshold=5.0, \
4769 4670 lownoisethreshold=2.0, negativethreshold=0.0, \
4770 4671 minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4771 4672 minpercentchange=1.0, fastnoise=False, restart=True, \
4772 4673 savemodel='none', calcres=False, calcpsf=False, \
4773 4674 parallel=True, verbose=True)
4774 -​
4775 - # retrieve per-channel beam statistics
4776 - bmin_dict, bmaj_dict, pa_dict = \
4777 - self.cube_beam_stats(img+'.image')
4778 -​
4779 - tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4780 - antenna=['0,1,2,3,4,5,6,7,8'], scan=['8,12,16'], \
4781 - intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4782 - imagename=file_name+'1', imsize=[108, 108], \
4783 - cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4784 - ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4785 - start='220.2526743594GHz', width='0.2441741MHz', \
4786 - outframe='LSRK', perchanweightdensity=True, \
4787 - gridder='mosaic', mosweight=True, \
4788 - usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4789 - restoration=True, restoringbeam='common', pbcor=True, \
4790 - weighting='briggs', robust=0.5, npixels=0, niter=0, \
4791 - threshold='0.354Jy', nsigma=0.0, interactive=0, usemask='auto'
4792 - '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
4793 - lownoisethreshold=2.0, negativethreshold=0.0, \
4794 - minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4795 - minpercentchange=1.0, fastnoise=False, restart=True, \
4796 - savemodel='none', calcres=False, calcpsf=False,
4797 - parallel=False, verbose=True)
4798 -​
4799 4675 else:
4800 4676 tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
4801 4677 antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
4802 4678 intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
4803 4679 imagename=file_name+'1', imsize=[108, 108], \
4804 4680 cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
4805 4681 ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
4806 4682 start='220.2526743594GHz', width='0.2441741MHz', \
4807 4683 outframe='LSRK', perchanweightdensity=True, \
4808 4684 gridder='mosaic', mosweight=True, \
4809 4685 usepointing=False, pblimit=0.2, deconvolver='hogbom', \
4810 4686 restoration=True, pbcor=True, weighting='briggs', robust=0.5,\
4811 4687 npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
4812 4688 interactive=0, usemask='auto-multithresh', \
4813 4689 sidelobethreshold=1.25, noisethreshold=5.0, \
4814 4690 lownoisethreshold=2.0, negativethreshold=0.0, \
4815 4691 minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
4816 4692 minpercentchange=1.0, fastnoise=False, restart=True, \
4817 4693 savemodel='none', calcres=False, calcpsf=False, \
4818 4694 restoringbeam='common', parallel=False, verbose=True)
4819 4695 ​
4820 - # retrieve per-channel beam statistics
4821 - bmin_dict, bmaj_dict, pa_dict = \
4822 - self.cube_beam_stats(img+'.psf')
4696 + # retrieve per-channel beam statistics
4697 + bmin_dict, bmaj_dict, pa_dict = \
4698 + self.cube_beam_stats(img+'.psf')
4823 4699 ​
4824 4700 report0 = th.checkall(imgexist = self.image_list(img, 'mosaic'))
4825 4701 ​
4826 4702 # .image report (test_mosaic_cube_pcwdT)
4827 4703 im_stats_dict = self.image_stats(img+'.image', fit_region = \
4828 4704 'ellipse[[11.48661818deg, -73.26292371deg], [8.2211arcsec, 7.4698arcsec], 90.00000000deg]', field_regions = \
4829 4705 ['circle[[00:45:54.383559, -73.15.29.41306], 22.45arcsec]',
4830 4706 'circle[[00:45:49.435664, -73.15.35.13742], 22.45arcsec]',
4831 4707 'circle[[00:45:53.057440, -73.15.50.79016], 22.45arcsec]',
4832 4708 'circle[[00:45:50.762696, -73.15.13.76177], 22.45arcsec]',
5170 5046 restoration=True, restoringbeam='common', \
5171 5047 pbcor=True, weighting='briggsbwtaper', robust=0.5,\
5172 5048 npixels=0, niter=20000, threshold='0.354Jy', nsigma=0.0, \
5173 5049 interactive=0, usemask='auto-multithresh', \
5174 5050 sidelobethreshold=1.25, noisethreshold=5.0, \
5175 5051 lownoisethreshold=2.0, negativethreshold=0.0, \
5176 5052 minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
5177 5053 minpercentchange=1.0, fastnoise=False, restart=True, \
5178 5054 savemodel='none', calcres=False, calcpsf=False, \
5179 5055 parallel=True, verbose=True)
5180 -​
5181 - # retrieve per-channel beam statistics
5182 - #bmin_dict, bmaj_dict, pa_dict = \
5183 - # self.cube_beam_stats(img+'.image')
5184 -​
5185 - #tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
5186 - # antenna=['0,1,2,3,4,5,6,7,8'], scan=['8,12,16'], \
5187 - # intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
5188 - # imagename=file_name+'1', imsize=[108, 108], \
5189 - # cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
5190 - # ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
5191 - # start='220.2526743594GHz', width='0.2441741MHz', \
5192 - # outframe='LSRK', perchanweightdensity=True, \
5193 - # gridder='mosaic', mosweight=True, \
5194 - # usepointing=False, pblimit=0.2, deconvolver='hogbom', \
5195 - # restoration=True, restoringbeam='common', pbcor=True, \
5196 - # weighting='briggsbwtaper', robust=0.5, npixels=0, niter=0, \
5197 - # threshold='0.354Jy', nsigma=0.0, interactive=0, usemask='auto'
5198 - # '-multithresh', sidelobethreshold=1.25, noisethreshold=5.0, \
5199 - # lownoisethreshold=2.0, negativethreshold=0.0, \
5200 - # minbeamfrac=0.1, growiterations=75, dogrowprune=True, \
5201 - # minpercentchange=1.0, fastnoise=False, restart=True, \
5202 - # savemodel='none', calcres=False, calcpsf=False,
5203 - # parallel=False, verbose=True)
5204 -​
5205 5056 else:
5206 5057 tclean(vis=self.msfile, field='SMIDGE_NWCloud', spw=['0'], \
5207 5058 antenna=['0,1,2,3,4,5,6,7,8'],scan=['8,12,16'], \
5208 5059 intent='OBSERVE_TARGET#ON_SOURCE', datacolumn='data', \
5209 5060 imagename=file_name+'1', imsize=[108, 108], \
5210 5061 cell=['1.1arcsec'], phasecenter='ICRS 00:45:54.3836'
5211 5062 ' -073.15.29.413', stokes='I', specmode='cube', nchan=508, \
5212 5063 start='220.2526743594GHz', width='0.2441741MHz', \
5213 5064 outframe='LSRK', perchanweightdensity=True, \
5214 5065 gridder='mosaic', mosweight=True, \

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

Add shortcut