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) ax.set_aspect(1) ax.axis([x1, x2, y1, y2]) for zoomregion in zoomregions_order: ZR = zoomregions[zoomregion] if ZR['inregion'] != regionname: print( f"{zoomregion} is not in {regionname}, it's in {ZR['inregion']}" ) continue
dec.set_ticks(exclude_overlapping=True) im = ax.imshow(hdu.data.squeeze(), transform=ax.get_transform(wcs.WCS(hdu.header).celestial), vmin=2300, vmax=55000, cmap=pl.cm.gray_r, origin='lower', norm=asinh_norm.AsinhNorm()) tr_fk5 = ax.get_transform("fk5") (x1, y1), (x2, y2) = (206, 174), (348, 333) ax.axis([x1, x2, y1, y2]) coredots = plotcores(ax, alpha=0.5, transform=tr_fk5, markersize=markersize, zorder=50) cb = pl.colorbar(mappable=im) cb.set_label("$S_{450 \mu m}$ [MJy sr$^{-1}$]") fig3.savefig(paths.fpath("cores_on_SCUBA_feathered.png"), bbox_inches='tight') im = ax.imshow(hdu.data.squeeze(), transform=ax.get_transform(wcs.WCS(hdu.header).celestial), vmin=2300, vmax=25000, cmap=pl.cm.gray_r, origin='lower', norm=asinh_norm.AsinhNorm()) cb.set_clim(im.get_clim()) cb.draw_all()