forked from akleroy/phangs_imaging_scripts
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathprocess_singledish.py
77 lines (55 loc) · 3.07 KB
/
process_singledish.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
# run in working_dir = '../singledish/data/'
os.system('rm -rf ALMA_TP.M74.v0p2.image')
importfits(fitsimage='ALMA_TP.M74.v0p2.blanked.fits',
imagename='ALMA_TP.M74.v0p2.image')
imreframe(imagename='ALMA_TP.M74.v0p2.image',
restfreq='230.53800GHz')
os.system('rm -rf ALMA_TP.M74.v0p2_rebin.image')
imrebin(imagename='ALMA_TP.M74.v0p2.image',
outfile='ALMA_TP.M74.v0p2_rebin.image',
factor=[1,1,7,1])
exportfits(imagename='ALMA_TP.M74.v0p2_rebin.image',
fitsimage='ALMA_TP.M74.v0p2_rebin.fits',
velocity=True, dropdeg=True, overwrite=True)
# import/export to strip degenerate axes
os.system('rm -rf ALMA_TP.NGC_1087.CO21.v0p2.gildas.image')
importfits(fitsimage='ALMA_TP.NGC_1087.CO21.v0p2.gildas.blanked.fits',
imagename='ALMA_TP.NGC_1087.CO21.v0p2.gildas.image')
exportfits(imagename='ALMA_TP.NGC_1087.CO21.v0p2.gildas.image',
fitsimage='ALMA_TP.NGC_1087.CO21.v0p2.casa.fits',
velocity=True, dropdeg=True, overwrite=True)
# import/export to strip degenerate axes, delete "BLANK" keyword
importfits(fitsimage='ALMA_TP.NGC_1300.CO21.v0p2.image.VLSRK.gildas.blanked.fits',
imagename='ALMA_TP.NGC_1300.CO21.v0p2.image.VLSRK.gildas.image'
, overwrite=True)
exportfits(imagename='ALMA_TP.NGC_1300.CO21.v0p2.image.VLSRK.gildas.image',
fitsimage='ALMA_TP.NGC_1300.CO21.v0p2.image.VLSRK.casa.fits',
velocity=True, dropdeg=True, overwrite=True)
# import/export to strip degenerate axes, delete "BLANK" keyword
importfits(fitsimage='ALMA_TP.NGC_1385.CO21.v0p2.image.VLSRK.gildas.blanked.fits',
imagename='ALMA_TP.NGC_1385.CO21.v0p2.image.VLSRK.gildas.image'
, overwrite=True)
exportfits(imagename='ALMA_TP.NGC_1385.CO21.v0p2.image.VLSRK.gildas.image',
fitsimage='ALMA_TP.NGC_1385.CO21.v0p2.image.VLSRK.casa.fits',
velocity=True, dropdeg=True, overwrite=True)
# import/export to strip degenerate axes, delete "BLANK" keyword
importfits(fitsimage='ALMA_TP.NGC_1433.CO21.v0p2.gildas.NuclearFix.blanked.fits',
imagename='ALMA_TP.NGC_1433.CO21.v0p2.gildas.NuclearFix.image'
, overwrite=True)
exportfits(imagename='ALMA_TP.NGC_1433.CO21.v0p2.gildas.NuclearFix.image',
fitsimage='ALMA_TP.NGC_1433.CO21.v0p2.casa.fits',
velocity=True, dropdeg=True, overwrite=True)
# import/export to strip degenerate axes, delete "BLANK" keyword
importfits(fitsimage='ALMA_TP.NGC_1566.CO21.v0p2.gildas.blanked.fits',
imagename='ALMA_TP.NGC_1566.CO21.v0p2.gildas.image',
overwrite=True)
exportfits(imagename='ALMA_TP.NGC_1566.CO21.v0p2.gildas.image',
fitsimage='ALMA_TP.NGC_1566.CO21.v0p2.casa.fits',
velocity=True, dropdeg=True, overwrite=True)
# redo export
importfits(fitsimage='ALMA_TP.NGC1365.CO21.v0p2.image.fits',
imagename='ALMA_TP.NGC1365.CO21.v0p2.image',
overwrite=True)
exportfits(imagename='ALMA_TP.NGC1365.CO21.v0p2.image',
fitsimage='ALMA_TP.NGC1365.CO21.v0p2.image.casa.fits',
velocity=True, dropdeg=True, overwrite=True)