Python YSO.c2dGBS_v2示例

示例#1

文件： CephBPLOT.py 项目： DJRumble/PHD

        "> /dev/null",
    )
    os.system(cmd)

    ndf2fits.ndf2fits(mask)

    Fig1.show_contour(maskFITS, levels=(0.5, 1000), linewidth=3, colors=("black"))
    Fig2.show_contour(maskFITS, levels=(0.5, 1000), linewidth=3, colors=("black"))
    Fig3.show_contour(maskFITS, levels=(0.5, 1000), linewidth=3, colors=("black"))

    cmd = "rm %s %s %s" % (thresh, mask, magic)
    os.system(cmd)


# Plot YSOs
YSO.c2dGBS_v2(Fig1)
YSO.c2dGBS_v2(Fig2)
YSO.c2dGBS_v2(Fig3)

# PLOT OB stars - primary stars in cyan
Fig1.show_markers(RA_OB, DEC_OB, edgecolor="k", facecolor="none", marker="*", s=250, alpha=1)
Fig2.show_markers(RA_OB, DEC_OB, edgecolor="k", facecolor="none", marker="*", s=250, alpha=1)
Fig3.show_markers(RA_OB, DEC_OB, edgecolor="k", facecolor="none", marker="*", s=250, alpha=1)

Fig1.show_markers(RA_OBx, DEC_OBx, edgecolor="k", facecolor="cyan", marker="*", s=250, alpha=1)
Fig2.show_markers(RA_OBx, DEC_OBx, edgecolor="k", facecolor="cyan", marker="*", s=250, alpha=1)
Fig3.show_markers(RA_OBx, DEC_OBx, edgecolor="k", facecolor="cyan", marker="*", s=250, alpha=1)

# Recenter
Fig1.recenter(pos1[0], pos1[1], pos1[2], pos1[3])
Fig2.recenter(pos2[0], pos2[1], pos2[2], pos2[3])

示例#2

文件： OphPLOT_temp.py 项目： DJRumble/PHD

RA_OB = data1[:,0]	
DEC_OB = data1[:,1]

data2 = np.loadtxt('OBstars.txt',dtype='string')
RA_OBx = map(float, data2[:,2])	
DEC_OBx = map(float, data2[:,3])

Fig1 = aplpy.FITSFigure(temp)
Fig1.show_colorscale(vmax=40.0,vmin=10.0, stretch='linear',cmap='CMRmap')

#PLOT contours
#SCUBA-2
Fig1.show_contour(s850, levels=(sig850,1000), linewidth=3, colors=('k'))

#Plot YSOs
YSO.c2dGBS_v2(Fig1)

#PLOT OB stars - primary stars in cyan
Fig1.show_markers(RA_OB,DEC_OB, edgecolor='k', facecolor='none', marker='*',s=250, alpha=1)
Fig1.show_markers(RA_OBx,DEC_OBx, edgecolor='k', facecolor='cyan', marker='*',s=250, alpha=1)

#Recenter
Fig1.recenter(pos1[0],pos1[1],pos1[2],pos1[3])

#Label
Fig1.add_label(0.65, 0.49, 'A', relative=True)
Fig1.add_label(0.22, 0.42, 'B', relative=True)
Fig1.add_label(0.45, 0.26, 'C', relative=True)
Fig1.add_label(0.4, 0.15, 'E', relative=True)
Fig1.add_label(0.25, 0.08, 'F', relative=True)