def T(fig,H70,s21): #poster image - 850 + H70/21cm contours + CO #deffine image in APLPY Fig = aplpy.FITSFigure(fig) YSO.w40_radio(Fig) OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt') ra = OBmarkers[:,0] dec = OBmarkers[:,1] Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=100, alpha=1) #cset4 = Fig.show_contour(H70,levels=(300,1200,4800,12000),colors=('blue')) #Herschel 70um data. cset1 = Fig.show_contour(s21, levels=(0.01,0.03,0.05,0.1), linewidth=4, colors=('red')) #Units: Jy/Beam, Resolution: 45'' #Figure properties Fig.show_grayscale(vmax=0.027,vmin=0., stretch='linear',invert='TRUE') Fig.tick_labels.set_font(size='medium') Fig.add_colorbar() #Fig.colorbar.set_axis_label_text('$850\,(K)$') Fig.recenter(277.8925,-2.13200, height=0.51 ,width=0.42) Fig.ticks.show() Fig.ticks.set_color('black') return
def K(before,after,freefree): big = plt.figure(figsize=(18, 12)) #deffine image in APLPY Before = aplpy.FITSFigure(before,figure=big,subplot=[0.1,0.1,0.4,0.8]) After = aplpy.FITSFigure(after,figure=big,subplot=[0.5,0.1,0.4,0.8]) YSO.w40_radio(Before) YSO.w40_radio(After) YSO.w40(Before) YSO.w40(After) #Figure properties Before.show_grayscale(vmax=0.032,vmin=0.0,invert=True,stretch='sqrt') #After.show_grayscale(vmax=0.032,vmin=0.0,invert=True,stretch='sqrt') After.show_grayscale(vmax=0.022,vmin=0.001375, stretch='sqrt',invert='TRUE') #Deffine contours cset = Before.show_contour(freefree, levels=(0.001,0.005,0.01), linewidth=(6,4,2), colors=('r')) cset2 = After.show_contour(freefree, levels=(0.001,0.005,0.01), linewidth=(6,4,2), colors=('r')) #Before.recenter(277.8925,-2.08200, height=0.1 ,width=0.15) #After.recenter(277.8925,-2.08200, height=0.1 ,width=0.15) Before.recenter(277.8525,-2.10200, height=0.03 ,width=0.035) After.recenter(277.8525,-2.10200, height=0.03 ,width=0.035) Before.tick_labels.set_font(size='medium') After.tick_labels.set_font(size='medium') Before.add_label(0.1, 0.95, r'SCUBA-2 450um',size='large', relative=True) After.add_label(0.1, 0.95, r'VLA 3.6cm',size='large', relative=True) Before.ticks.show() Before.ticks.set_color('black') After.axis_labels.hide_y() After.tick_labels.hide_y() After.ticks.show() After.ticks.set_color('black') #After.add_colorbar() #After.colorbar.set_axis_label_text('850$\mu m$ (Jy/pixel)') big.canvas.draw() return
def S(fig,H70,HII1,HII2):#,mask): #W40 temp map with Herschel Contours #deffine image in APLPY Fig = aplpy.FITSFigure(fig) #Deffine Markers #YSO.GBS(Fig) #YSO.w40(Fig) YSO.w40_radio(Fig) OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt') ra = OBmarkers[:,0] dec = OBmarkers[:,1] Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=100, alpha=1) #Deffine contours cset1 = Fig.show_contour(HII1, levels=(0.05,0.075,0.1,0.2), linewidth=4, colors=('blue')) #Units: Jy/Beam, Resolution: 45'' #cset2 = Fig.show_contour(HII2, levels=(0.001,0.005,0.0075,0.01), linewidth=2, colors=('magenta')) #Units: Jy/Beam, Resolution: 45'' #cset3 = Fig.show_contour(COmask, levels=1, linewidth=2, colors=('black')) #Mask of VanDerWeil2012 CO #cset4 = Fig.show_contour(H70,levels=(300,1200,4800,12000),colors=('k')) #Herschel 70um data. #cset3 = Fig.show_contour(flux, levels=(0.0027,0.0081,0.027), linewidth=2, colors=('black','black','magenta')) #cset5 = Fig.show_contour(YSO,levels=(15,45,75),colors=('magenta')) #Figure properties Fig.show_colorscale(vmax=30.0,vmin=10, stretch='sqrt') Fig.tick_labels.set_font(size='medium') Fig.add_colorbar() Fig.colorbar.set_axis_label_text('$T_{d}\,(K)$') Fig.recenter(277.8925,-2.13200, height=0.51 ,width=0.42) Fig.ticks.show() Fig.ticks.set_color('black') return
def I(fig,red,blue): #deffine image in APLPY Fig = aplpy.FITSFigure(fig) #Deffine contours cset_red = Fig.show_contour(red, levels=(5,10,15,20,25,30,35), linewidth=3, colors=('r')) cset_red2 = Fig.show_contour(red, levels=(30,35,100), linewidth=3, filled=('TRUE'), colors=('r')) cset_blue = Fig.show_contour(blue, levels=(5,10,15,20,25,30,35), linewidth=3, colors=('b')) cset_blue2 = Fig.show_contour(blue, levels=(30,35,100), linewidth=3, filled=('TRUE'), colors=('b')) cset = Fig.show_contour(fig, levels=(0.0027,0.0081,0.027,0.081), linewidth=3, colors=('black')) #Figure properties Fig.show_grayscale(vmax=0.027,vmin=0., stretch='linear',invert='TRUE') #Fig.add_label(0.07, 0.95, r'$850\mu m$',size='large', relative=True) #Deffine Markers OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt') ra = OBmarkers[:,0] dec = OBmarkers[:,1] Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=25, alpha=1) YSO.w40_radio(Fig) Fig.tick_labels.set_font(size='medium') Fig.add_colorbar() #Fig.colorbar.set_location('top') Fig.colorbar.set_axis_label_text('850$\mu m$ (Jy/pixel)') Fig.recenter(277.86083,-2.0885,height=0.1 ,width=0.1) Fig.ticks.show() Fig.ticks.set_color('white') return
def N(fig,SCUBA2): Fig = aplpy.FITSFigure(fig) YSO.w40(Fig) YSO.w40_radio(Fig) #3.6cm Fig.show_grayscale(vmax=0.022,vmin=0.001375, stretch='sqrt',invert='TRUE') #21cm #Fig.show_colorscale(vmax=0.025,vmin=0.0015, stretch='linear') cset3 = Fig.show_contour(SCUBA2, levels=(0.0027,0.0081,0.027), linewidth=2, colors=('black','black','magenta')) Fig.add_colorbar() Fig.colorbar.set_axis_label_text('3.6 cm (Jy/Beam)') Fig.ticks.show() return
def U(small,large): #deffine image in APLPY Fig = aplpy.FITSFigure(large) cset1 = Fig.show_contour(small, levels=(0.001,0.005,0.0075,0.01), linewidth=4, colors=('red')) #Units: Jy/Beam, Resolution: 45'' cset2 = Fig.show_contour(large, levels=(0.001,0.005,0.0075,0.01), linewidth=4, colors=('lime')) #Units: Jy/Beam, Resolution: 45'' OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt') ra = OBmarkers[:,0] dec = OBmarkers[:,1] Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=30, alpha=1) YSO.w40_radio(Fig) Fig.ticks.set_color('black') return
def Q(fig,SCUBA2): #VLA 3.6cm map Fig = aplpy.FITSFigure(fig) YSO.w40_radio(Fig) OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt') ra = OBmarkers[:,0] dec = OBmarkers[:,1] Fig.show_markers(ra,dec, edgecolor='k', facecolor='yellow', marker='*',s=200, alpha=1) cset850 = Fig.show_contour(SCUBA2, levels=(0.0025,0.0075,0.025),colors='k') Fig.show_grayscale(vmax=0.022,vmin=0.001375, stretch='sqrt',invert='TRUE') Fig.add_colorbar() Fig.recenter(277.8525,-2.09200, height=0.06 ,width=0.07) Fig.colorbar.set_axis_label_text('$S_{3.6\,cm}$ (Jy/beam)') Fig.add_beam(9.97*u.arcsecond,9.97*u.arcsecond,24.42*u.degree,color='black') Fig.ticks.show() Fig.ticks.set_color('black') return
def L(fig,HII1):#,mask): #W40 with CO observstion areas #deffine image in APLPY Fig = aplpy.FITSFigure(fig) #Deffine Markers #YSO.GBS(Fig) #YSO.oth(Fig) YSO.w40_radio(Fig) #YSO.w40(Fig) #Deffine contour cset1 = Fig.show_contour(HII1, levels=(0.01,0.05,0.1), linewidth=4, colors=('blue')) #Units: Jy/Beam, Resolution: 45'' #cset2 = Fig.show_contour(HII2, levels=(0.0015,0.015,0.0175,0.02), linewidth=2, colors=('cyan')) #Units: Jy/Beam, Resolution: 45'' cset3 = Fig.show_contour(fig, levels=(0.0027,0.0081,0.027), linewidth=2, colors=('black','black','magenta')) #Figure properties Fig.show_grayscale(vmax=0.054,vmin=0.0,invert=True) Fig.add_label(0.07, 0.95, r'$850\mu m$',size='large', relative=True) Fig.tick_labels.set_font(size='medium') Fig.add_colorbar() Fig.ticks.show() Fig.ticks.set_color('black') Fig.recenter(277.8925,-2.13200, height=0.51 ,width=0.42) return
def P(s450,s850,s450ff,s850ff): big = plt.figure(figsize=(16, 12)) #deffine image in APLPY S450 = aplpy.FITSFigure(s450,figure=big,subplot=[0.1,0.1,0.4,0.8]) S850 = aplpy.FITSFigure(s850,figure=big,subplot=[0.5,0.1,0.4,0.8]) #Figure properties S450.show_grayscale(vmax=0.185,vmin=0., stretch='linear',invert='TRUE') S850.show_grayscale(vmax=0.027,vmin=0., stretch='linear',invert='TRUE') #Deffine contours cset = S450.show_contour(s450ff, levels=(0.0104,0.0173,0.0519,100), colors=('r','k'),filled='TRUE') cset2 = S850.show_contour(s850ff, levels=(0.0015,0.0025,0.0075,100), colors=('r','k'),filled='TRUE') cset450 = S450.show_contour(s450, levels=(0.0104,0.0173,0.0519,0.1110), colors=('k')) cset850 = S850.show_contour(s850, levels=(0.0015,0.0025,0.0075,0.0161), colors=('k')) #recenter S450.recenter(277.8525,-2.10200, height=0.025 ,width=0.05) S850.recenter(277.8525,-2.10200, height=0.025 ,width=0.05) YSO.w40_radio(S450) #YSO.w40_radio(S850) OBmarkers = np.loadtxt('/data/damian/data/aquila/OBobjects.txt') ra = OBmarkers[:,0] dec = OBmarkers[:,1] S850.show_markers(ra,dec, layer='W40a', edgecolor='k', facecolor='yellow', marker='*',s=150, alpha=1) S450.show_markers(ra,dec, layer='W40b', edgecolor='k', facecolor='yellow', marker='*',s=150, alpha=1) ra = 277.849937 dec = -2.091536 S450.show_markers(ra,dec, edgecolor='k', facecolor='blue', marker='*',s=100, alpha=1) S850.show_markers(ra,dec, edgecolor='k', facecolor='blue', marker='*',s=100, alpha=1) #850 set up S850.axis_labels.hide_y() S850.tick_labels.hide_y() S850.ticks.show() S850.ticks.set_color('black') S850.add_colorbar() S850.colorbar.set_font(size='small') S850.colorbar.set_location('top') S850.colorbar.set_axis_label_text('850$\mu m$ (Jy/pixel)') S850.add_beam(14.5*u.arcsecond,14.5*u.arcsecond,0,color='black') S850.tick_labels.set_font(size='small') #450 set up S450.add_colorbar() S450.colorbar.set_location('top') S450.colorbar.set_axis_label_text('450$\mu m$ (Jy/pixel)') S450.colorbar.set_font(size='small') S450.add_beam(9.8*u.arcsecond,9.8*u.arcsecond,0,color='black') S450.ticks.show() S450.ticks.set_color('black') S450.tick_labels.set_font(size='small') big.canvas.draw() return