Simplify test_wcsgeom_read_write a bit
@@ -88,9 +88,8 @@ def test_wcsgeom_read_write(tmpdir, npix, binsz, coordsys, proj, skydir, axes):
geom0 = WcsGeom.create(npix=npix, binsz=binsz,
proj=proj, coordsys=coordsys, axes=axes)
- shape = (np.max(geom0.npix[0]), np.max(geom0.npix[1]))
hdu_bands = geom0.make_bands_hdu(hdu='BANDS')
- hdu_prim = fits.PrimaryHDU(np.zeros(shape).T)
+ hdu_prim = fits.PrimaryHDU()
hdu_prim.header.update(geom0.make_header())
filename = str(tmpdir / 'wcsgeom.fits')
