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)
