fitsimages = glob.glob(dir_data+"data/"+galname+"*.fits") for i in range(len(fitsimages)): myim2.eazy_importfits(fitsimages[i],defaultaxes=False) # find imported images image_co32 = glob.glob(dir_data+"data/"+galname+"*co32*image*")[0] ### step 2/10: Kelvin to Jansky conversion print("### step 2/10: Kelvin to Jansky conversion") bunit = imhead(image_co32,"list")["bunit"] synsbeam21 = imhead(image_co32,"list")["beammajor"]["value"] if bunit == "K": # CO(2-1) conversion if bunit = K myim2.easy_K2Jy(image_co32,synsbeam21,230.53800) image_co32 = image_co32 + ".jy" else: print("# skip easy_K2Jy for the CO(2-1) data") ### step 3/10: imregrid print("### step 3/10: imregrid") pbimage = glob.glob(dir_data+"data/"+galname+"*pb.image")[0] myim2.easy_imregrid(pbimage,image_co32,False) # pbimage pbimage = pbimage + ".regrid"
print("### step 1/10: importfits") fitsimages = glob.glob(dir_data + "data/" + galname + "*.fits") for i in range(len(fitsimages)): myim2.eazy_importfits(fitsimages[i], defaultaxes=False) # find imported images image_co21 = glob.glob(dir_data + "data/" + galname + "*image*")[0] ### step 2/10: Kelvin to Jansky conversion print("### step 2/10: Kelvin to Jansky conversion") bunit = imhead(image_co21, "list")["bunit"] synsbeam21 = imhead(image_co21, "list")["beammajor"]["value"] if bunit == "K": # CO(2-1) conversion if bunit = K myim2.easy_K2Jy(image_co21, synsbeam21, restfreq) image_co21 = image_co21 + ".jy" else: print("# skip easy_K2Jy for the CO(2-1) data") ### step 3/10: imregrid print("### step 3/10: imregrid") #pbimage = glob.glob(dir_data+"data/"+galname+"*pb.image")[0] #myim2.easy_imregrid(pbimage,image_co21,False) # pbimage #pbimage = pbimage + ".regrid" ### step 4/10: imsmooth print("### step 4/10: imsmooth") beam_mask = beam_size * increment_mask # beam size for the masking
fitsimages = glob.glob(dir_data + "data/" + galname + "*.fits") for i in range(len(fitsimages)): myim2.eazy_importfits(fitsimages[i]) # find imported images image_co10 = glob.glob(dir_data + "data/" + galname + "*co10*image*")[0] image_co21 = glob.glob(dir_data + "data/" + galname + "*co21*image*")[0] ### step 2/10: Kelvin to Jansky conversion print("### step 2/10: Kelvin to Jansky conversion") bunit = imhead(image_co10, "list")["bunit"] synsbeam10 = imhead(image_co10, "list")["beammajor"]["value"] if bunit == "K": # CO(1-0) conversion if bunit = K myim2.easy_K2Jy(image_co10, synsbeam10, 115.27120) image_co10 = image_co10 + ".jy" else: print("# skip easy_K2Jy for the CO(1-0) data") bunit = imhead(image_co21, "list")["bunit"] synsbeam21 = imhead(image_co21, "list")["beammajor"]["value"] if bunit == "K": # CO(2-1) conversion if bunit = K myim2.easy_K2Jy(image_co21, synsbeam21, 230.53800) image_co21 = image_co21 + ".jy" else: print("# skip easy_K2Jy for the CO(2-1) data") ### step 3/10: imregrid print("### step 3/10: imregrid")
print("### step 1/10: importfits") fitsimages = glob.glob(dir_data + "images/" + galname + "*.fits") for i in range(len(fitsimages)): myim2.eazy_importfits(fitsimages[i]) # find imported images image_co32 = glob.glob(dir_data + "images/" + galname + "*.image")[0] ### step 2/10: Kelvin to Jansky conversion print("### step 2/10: Kelvin to Jansky conversion") bunit = imhead(image_co32, "list")["bunit"] synsbeam10 = imhead(image_co32, "list")["beammajor"]["value"] if bunit == "K": # CO(1-0) conversion if bunit = K myim2.easy_K2Jy(image_co32, synsbeam10, 115.27120) image_co32 = image_co32 + ".jy" else: print("# skip easy_K2Jy for the CO(3-2) data") ### step 4/10: imsmooth print("### step 4/10: imsmooth") beam_mask = beam_size * increment_mask # beam size for the masking myim2.easy_imsmooth(image_co32, beam_mask, False) # co32 ### mv the cubes to the working directory os.system("mkdir " + dir_data + galname) os.system("mv "+dir_data+"images/"+galname+"*32*smooth "\ +dir_data+galname+"/"+galname+"_co32_"+suffix+".cube")
print("### step 1/10: importfits") fitsimages = glob.glob(dir_data + "data/" + galname + "*.fits") for i in range(len(fitsimages)): myim2.eazy_importfits(fitsimages[i], defaultaxes=False) # find imported images image_co10 = glob.glob(dir_data + "data/" + galname + "*co10*image*")[0] ### step 2/10: Kelvin to Jansky conversion print("### step 2/10: Kelvin to Jansky conversion") bunit = imhead(image_co10, "list")["bunit"] synsbeam10 = imhead(image_co10, "list")["beammajor"]["value"] if bunit == "K": # CO(2-1) conversion if bunit = K myim2.easy_K2Jy(image_co10, synsbeam10, 1900.53690000 / (1 + 4.69)) image_co10 = image_co10 + ".jy" else: print("# skip easy_K2Jy for the CO(2-1) data") ### step 3/10: imregrid print("### step 3/10: imregrid") #pbimage = glob.glob(dir_data+"data/"+galname+"*pb.image")[0] #myim2.easy_imregrid(pbimage,image_co10,False) # pbimage #pbimage = pbimage + ".regrid" ### step 4/10: imsmooth print("### step 4/10: imsmooth") beam_mask = beam_size * increment_mask # beam size for the masking