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()