Running
thumb_shape, thumb_wcs = enmap.thumbnail_geometry(r=cutout_size_deg*utils.degree, res=cutout_resolution_deg*utils.degree)
with (cutout_size_deg, cutout_resolution_deg) = (10, 0.5), I get a shape of (43, 41) and a wcs of tan:{cdelt:[-0.5,0.5],crval:[0,0],crpix:[21.00,22.00]}.
Whereas actually getting the shape from a thumbnail, given some imap,
thumb_base = reproject.thumbnails(imap, coords=np.deg2rad([0,0]).T, res=cutout_resolution_deg*utils.degree, r=cutout_size_deg*utils.degree)
thumb_shape, thumb_wcs = thumb_base.shape, thumb_base.wcs
gives me a thumb_shape of (41, 41) and wcs of car:{cdelt:[-0.5,0.5],crval:[0,0],crpix:[21.00,21.00]}. Is there a reason these are unequal and that the first case returns a non-square array?
I am using pixell 0.28.4 and python 3.10.2.