def aplpy_centered_gaussian(): import aplpy import pylab as pl fig = pl.figure(1) fig.clf() F = aplpy.FITSFigure(data.centered_gaussian_hdu(), subplot=(2,3,1), figure=fig) F.show_grayscale() F._ax1.set_title("Input") hcr,ast,rep,mon = test_centered_gaussian() titles = ['hcongrid','wcsalign/starlink','pyreproject','montage'] figs = [] for ii,im in enumerate((hcr,ast,rep,mon)): figs.append(aplpy.FITSFigure(im,subplot=(2,3,ii+2), figure=fig)) figs[-1].show_grayscale() figs[-1].show_contour(data.centered_gaussian_hdu()) figs[-1]._ax1.set_title(titles[ii]) F.save('centered_gaussian_test.png')
def test_centered_gaussian(): hdu_in = data.centered_gaussian_hdu() hdr_out = headers.header_radec2 return all_tools(hdu_in,hdr_out)