casa / casa6 / 72d21ff7a17

Commits

Andrew McNichols authored 72d21ff7a17 02 Mar 2021

Fix pb math in task_simalma (CAS-13369)

The .image file produced by clean in 5.4 had primary beam correction applied. Since CAS-12502 and the change to using imtclean instead of imclean in simanalyze (called by simalma) the .image file is now sky-flat and does not need an extra scaling by the primary beam before it is feathered with the SD image.

Ensure the task_simalma code follows the algorithm:
* create SD image
* create IF flat-noise image
* regrid SD to IF
* output = feather(SD*IF.pb, IF(flat-noise) )/ IF.pb

casa5/gcwrap/python/scripts/task_simalma.py

Modified

-1537
+#                    ia.open(scaledimg)
-1538
+#                    ia.setrestoringbeam(beam=beam_tp)
-1539
+#                    ia.setbrightnessunit(bunit_tp)
-1540
+#                    ia.close()
-1541
-1542
+                msg("impbcor('"+regridimg+"', '"+pbcov+"', outfile='"+scaledimg+"',mode='multiply')",priority="info")
-1543
+                if not dryrun:
-1544
+                    impbcor(regridimg, pbcov, outfile=scaledimg,mode='multiply')
-1545
-1546
+                # de-pbcor the INT image
--
+                highimage = fileroot+"/"+imagename_int+".pbscaled"
-+
+                # not needed now that imtclean is used and .image from tclean is flat sky by default, see CAS-13370
-+
+                #highimage = fileroot+"/"+imagename_int+".pbscaled"
-1549
+                #immath(imagename=[highimage0, pbcov],
-1550
+                #       expr='IM1/IM0',outfile=highimage)
--
+                msg(" ",priority="info")
--
+                msg("Multiply interferometric image by sensitivity map (un-pbcor):",priority="info")
-+
+#                msg(" ",priority="info")
-+
+#                msg("Multiply interferometric image by sensitivity map (un-pbcor):",priority="info")
-1553
+#                msg("immath(imagename=['"+highimage0+"','"+pbcov+"'],expr='IM1*IM0',outfile='"+highimage+"')",priority="info")
-1554
+#                msg("Restore interferometric beam and brightness unit:",priority="info")
-1555
+#                msg("ia.open('"+highimage0+"')",priority="info")
-1556
+#                msg("beam_int = ia.restoringbeam()",priority="info")
-1557
+#                msg("bunit_int = ia.brightnessunit()",priority="info")
-1558
+#                msg("ia.close()",priority="info")
-1559
+#                msg("ia.open('"+highimage+"')",priority="info")
-1560
+#                msg("ia.setrestoringbeam(beam=beam_int)",priority="info")
-1561
+#                msg("ia.setbrightnessunit(bunit_int)",priority="info")
-1562
+#                msg("ia.close()",priority="info")
-1566
+#                           expr='IM1*IM0',outfile=highimage)
-1567
+#                    # restore INT beam and brightness unit
-1568
+#                    ia.open(highimage0)
-1569
+#                    beam_int = ia.restoringbeam()
-1570
+#                    bunit_int = ia.brightnessunit()
-1571
+#                    ia.close()
-1572
+#                    ia.open(highimage)
-1573
+#                    ia.setrestoringbeam(beam=beam_int)
-1574
+#                    ia.setbrightnessunit(bunit_int)
-1575
+#                    ia.close()
--
--
+                msg("impbcor('"+highimage0+"', '"+pbcov+"', outfile='"+highimage+"',mode='multiply')",priority="info")
--
+                if not dryrun:
--
+                    impbcor(highimage0, pbcov, outfile=highimage,mode='multiply')
-+
+#
-+
+#                msg("impbcor('"+highimage0+"', '"+pbcov+"', outfile='"+highimage+"',mode='multiply')",priority="info")
-+
+#                if not dryrun:
-+
+#                    impbcor(highimage0, pbcov, outfile=highimage,mode='multiply')
-1580
-1581
-1582
-1583
-1584
-1585
+                # Feathering
-1586
+                task_param = {}
-1587
+                task_param['imagename'] = outimage0
--
+                task_param['highres'] = highimage
-+
+                task_param['highres'] = highimage0
-1589
+                task_param['lowres'] = lowimage
-1590
-1591
+                msg(" ",priority="info")
-1592
+                msg(get_taskstr('feather', task_param), priority="info")
-1593
+                try:
-1594
+                    if not is_CASA6:
-1595
+                        saveinputs('feather',
-1596
+                                   fileroot+"/"+project+".feather.last",
-1597
+                                   myparams=task_param)
-1598
+                    else:

casa6

Commits

Jira Issues

Add shortcut