+galname+"*_co32_"+suffix+".cube")[0] thres_co32 = rms_co32 * increment_mask * thres_masking outmask_co32=cube_co32.replace(".cube",".mask") myim2.createmask(cube_co32,thres_co32,outmask_co32) beamarea = myim2.beam_area(image_co32,increment_mask) myim2.remove_smallmask(outmask_co32,beamarea,pixelmin) ### mv to working directory os.system("rm -rf "+cube_co32) os.system("mv "+image_co32+" "+cube_co32) print("### step 9/10: immoments") myim2.moment_maps(cube_co32,chans,outmask_co32,rms_co32*thres_mom) ### pbmask print("### step 10/10: pb mask at " + str(pbcut)) mask_pb = dir_data+galname+"/"+galname+"_pb_"+suffix+".mask" peak = imhead(pbimage,mode="list")["datamax"] myim2.createmask(pbimage,peak*pbcut,mask_pb) os.chdir(dir_data+galname) ia.open(glob.glob("*pb*mask")[0]) ia.calcmask(mask=glob.glob("*pb*mask")[0]+">0", name=glob.glob("*pb*mask")[0]+".2")
cube_co21 = glob.glob(dir_data+galname+"/"\ +galname+"_"+suffix+".cube")[0] thres_co21 = rms_co21 * increment_mask * thres_masking outmask_co21 = cube_co21.replace(".cube", ".mask") myim2.createmask(cube_co21, thres_co21, outmask_co21) beamarea = myim2.beam_area(image_co21, increment_mask) myim2.remove_smallmask(outmask_co21, beamarea, pixelmin) ### mv to working directory os.system("rm -rf " + cube_co21) os.system("mv " + image_co21 + " " + cube_co21) print("### step 9/10: immoments") myim2.moment_maps(cube_co21, chans, outmask_co21, rms_co21 * thres_mom) ### pbmask print("### step 10/10: pb mask at " + str(pbcut)) #mask_pb = dir_data+galname+"/"+galname+"_pb_"+suffix+".mask" #peak = imhead(pbimage,mode="list")["datamax"] #myim2.createmask(pbimage,peak*pbcut,mask_pb) #os.chdir(dir_data+galname) #ia.open(glob.glob("*pb*mask")[0]) #ia.calcmask(mask=glob.glob("*pb*mask")[0]+">0", # name=glob.glob("*pb*mask")[0]+".2") #ia.close() #os.chdir("/Users/saito/data/mycasaimaging_scripts/scripts_image_phangs2")
myim2.remove_smallmask(mask_combine, beamarea, pixelmin) """ ### step 8/10: imsmooth print("### step 8/10: imsmooth") myim2.easy_imsmooth(image_co32,beam_size,False) # co32 myim2.easy_imsmooth(image_co21,beam_size,False) # co21 """ ### mv to working directory os.system("rm -rf " + cube_co32) os.system("rm -rf " + cube_co21) os.system("cp -r " + dir_data + "data/" + galname + "*co32* " + cube_co32) print("### step 9/10: immoments") myim2.moment_maps(cube_co32, chans, mask_combine, rms_co32 * thres_mom) ### pbmask print("### step 10/10: pb mask at " + str(pbcut)) images_moment = glob.glob(cube_co32 + ".moment*") images_moment.extend(glob.glob(cube_co21 + ".moment*")) for i in range(len(images_moment)): os.system("rm -rf " + outfile) imagename = images_moment[i] outfile = images_moment[i].replace(".ICRS", "").replace(".cube.", "") imregrid(imagename=imagename, template="J2000", output=outfile) os.system("rm -rf " + imagename) # exportfits os.system("rm -rf " + outfile + ".fits")
myim2.createmask(cube_ci10, thres_ci10, outmask_ci10) beamarea = myim2.beam_area(image_ci10, increment_mask) myim2.remove_smallmask(outmask_ci10, beamarea, pixelmin) ### step 8/10: imsmooth print("### step 8/10: imsmooth") myim2.easy_imsmooth(image_ci10, beam_size, True) # ci10 ### mv to working directory os.system("rm -rf " + cube_ci10) os.system("mv " + dir_data + "fits/" + galname + "*10*smooth " + cube_ci10) print("### step 9/10: immoments") myim2.moment_maps(cube_ci10, chans, outmask_ci10, rms_ci10 * thres_mom) ### pbmask print("### step 10/10: pb mask at " + str(pbcut)) pbimage = glob.glob(dir_data + "fits/" + galname + "*.pb")[0] mask_pb = dir_data + galname + "/" + galname + "_ci10_" + suffix + "_pb.mask" peak = imhead(pbimage, mode="list")["datamax"] myim2.createmask(pbimage, peak * pbcut, mask_pb) images_moment = glob.glob(cube_ci10 + ".moment*") for i in range(len(images_moment)): outfile = images_moment[i].replace(".cube", "") os.system("rm -rf " + outfile) immath(imagename=[images_moment[i], mask_pb], mode="evalexpr",
myim2.remove_smallmask(mask_combine, beamarea, pixelmin) ### step 8/10: imsmooth print("### step 8/10: imsmooth") myim2.easy_imsmooth(image_co10, beam_size, False) # co10 myim2.easy_imsmooth(image_co21, beam_size, False) # co21 ### mv to working directory os.system("rm -rf " + cube_co10) os.system("rm -rf " + cube_co21) os.system("mv " + dir_data + "data/" + galname + "*10*smooth " + cube_co10) os.system("mv " + dir_data + "data/" + galname + "*21*smooth " + cube_co21) print("### step 9/10: immoments") myim2.moment_maps(cube_co10, chans, mask_combine, rms_co10 * thres_mom) myim2.moment_maps(cube_co21, chans, mask_combine, rms_co21 * thres_mom) ### pbmask print("### step 10/10: pb mask at " + str(pbcut)) mask_pb = dir_data + galname + "/" + galname + "_pb_" + suffix + ".mask" peak = imhead(pbimage, mode="list")["datamax"] myim2.createmask(pbimage, peak * pbcut, mask_pb) images_moment = glob.glob(cube_co10 + ".moment*") images_moment.extend(glob.glob(cube_co21 + ".moment*")) for i in range(len(images_moment)): outfile = images_moment[i].replace(".cube", "") os.system("rm -rf " + outfile) immath(imagename=[images_moment[i], mask_pb],