ra = ax.coords['ra'] ra.set_major_formatter('hh:mm:ss.s') dec = ax.coords['dec'] ra.set_axislabel("RA (J2000)", fontsize=pl.rcParams['axes.labelsize']) dec.set_axislabel("Dec (J2000)", fontsize=pl.rcParams['axes.labelsize'], minpad=0) ra.ticklabels.set_fontsize(tick_fontsize) ra.set_ticks(exclude_overlapping=True) dec.ticklabels.set_fontsize(tick_fontsize) dec.set_ticks(exclude_overlapping=True) scalebarpos = coordinates.SkyCoord("17:47:27.5", "-28:26:00.0", unit=(u.h, u.deg), frame='fk5') make_scalebar(ax, scalebarpos, length=(2.0*u.pc / distance).to(u.arcsec, u.dimensionless_angles()), color='k', label='2 pc', text_offset=1.0*u.arcsec, ) fig.savefig(paths.fpath("HC3N_contours_on_SCUBA_column.pdf"), bbox_inches='tight') core_phot_tbl = Table.read(paths.tpath("continuum_photometry.ipac"), format='ascii.ipac') cores = coordinates.SkyCoord(core_phot_tbl['RA'], core_phot_tbl['Dec'], frame='fk5') markersize = 6 tr_fk5 = ax.get_transform("fk5") coredots, = ax.plot(cores.ra, cores.dec, '.', color='lime', transform=tr_fk5, markersize=markersize, zorder=50) ax.axis([195,340,180,320])
interpolation='nearest', origin='lower', norm=asinh_norm.AsinhNorm()) tr_fk5 = ax.get_transform("fk5") #(x1,y1),(x2,y2) = (1200,434),(2142,1743) # wrong (x1,y1),(x2,y2) = tr_fk5.transform_point([bottomleft.ra.deg, bottomleft.dec.deg]),tr_fk5.transform_point([topright.ra.deg, topright.dec.deg]) (x1, y1), (x2, y2) = (toplevel_wcs.wcs_world2pix( [[bottomleft.ra.deg, bottomleft.dec.deg]], 0)[0], toplevel_wcs.wcs_world2pix( [[topright.ra.deg, topright.dec.deg]], 0)[0]) make_scalebar( ax, scalebarpos, length=(0.5 * u.pc / distance).to(u.arcsec, u.dimensionless_angles()), color='k', label='0.5 pc', text_offset=1.0 * u.arcsec, ) if regionname in ('full', 'fullN'): markersize = 1 coredots = plotcores(ax, alpha=1, transform=ax.get_transform('fk5'), markerfacecolor='none', markersize=markersize, zorder=50)
# anno = ax0.annotate('S', xy=(sgrb2s_x, sgrb2s_y), # xytext=(sgrb2s_txt_x, sgrb2s_txt_y), # color=textcolor, # arrowprops=dict(edgecolor=textcolor, # arrowstyle='-', linewidth=1)) scalebarpos = coordinates.SkyCoord("17:47:27", "-28:25:30", unit=(u.h, u.deg), frame='fk5') scalebarlen = (2 * u.pc / distance).to(u.arcsec, u.dimensionless_angles()) make_scalebar(ax0, scalebarpos, length=scalebarlen, color=textcolor, label='50" 2 pc', text_offset=1.0 * u.arcsec, fontsize=8) ra.set_ticks(size=-3, color='w') dec.set_ticks(size=-3, color='w') # cbar_ax = fig.add_axes( # [0.7, 0.11, 0.01, 0.51]) # cb = plt.colorbar(im, cax=cbar_ax, orientation='vertical') # cbar_ax.yaxis.set_ticks_position("right") # cbar_ax.yaxis.set_label_position('right') # cbar_ax.tick_params(axis='both', which='both', direction='in') # cb.set_label("$S$ [mJy beam$^{-1}$]") # cb.ax.zorder = 1e8
fig = pl.figure(1) fig.clf() ax = pl.subplot(1, 3, 1, projection=wcs.WCS(fh_s[0].header)) pl.imshow(sgrb2_hii_T_cutout.data, vmin=-1e-1, vmax=4e-0, cmap='viridis', origin='lower', interpolation='nearest') pl.title("Sgr B2 HII T") hide_labels(ax) make_scalebar( ax, left_side=coordinates.SkyCoord('266d54m39s', '-28d27m48.5s', unit=(u.deg, u.deg), frame='fk5'), length=(0.1 * u.pc / d_sgr).to(u.arcsec, u.dimensionless_angles()), label='0.1 pc', fontsize=8, ) pl.subplot(1, 3, 2) pl.imshow(orion_cutout.data * (d_ori / d_sgr).decompose()**2, vmin=-1e-1, vmax=4e-0, cmap='viridis', origin='lower', interpolation='nearest') pl.title("M42") pl.gca().get_xaxis().set_ticks([]) pl.gca().get_yaxis().set_ticks([])
vmin=vmin, vmax=vmax, cmap=cmaps[regionname], interpolation='nearest', origin='lower', norm=norms[regionname]) tr_fk5 = ax.get_transform("fk5") #(x1,y1),(x2,y2) = (1200,434),(2142,1743) # wrong (x1,y1),(x2,y2) = tr_fk5.transform_point([bottomleft.ra.deg, bottomleft.dec.deg]),tr_fk5.transform_point([topright.ra.deg, topright.dec.deg]) (x1,y1),(x2,y2) = (mywcs.wcs_world2pix([[bottomleft.ra.deg, bottomleft.dec.deg]],0)[0], mywcs.wcs_world2pix([[topright.ra.deg, topright.dec.deg]],0)[0] ) make_scalebar(ax, scalebarpos[regionname], length=(scalebarlength[regionname] / distance).to(u.arcsec, u.dimensionless_angles()), color='k', label=scalebarlabel[regionname], text_offset=1.0*u.arcmin, ) ax.set_aspect(1) ax.axis([x1,x2,y1,y2]) for zoomregion in zoomregions_order: ZR = zoomregions[zoomregion] if ZR['inregion'] != regionname: continue