Switch to area averaging when binning hips

Author: natelust

python/lsst/pipe/tasks/fits2hips/_high_order_hips.py

@@ -198,7 +198,7 @@ def run(self, input_images: Iterable[tuple[NDArray, SkyWcs, Box2I]], healpix_id)
             self.log.info("Generating tiles for hxp level %d", hips_level)
             if zoom:
                 size = 4096 // zoom
-                binned_array = cv2.resize(output_array_hpx, (size, size), interpolation=cv2.INTER_LANCZOS4)
+                binned_array = cv2.resize(output_array_hpx, (size, size), interpolation=cv2.INTER_AREA)
             else:
                 binned_array = output_array_hpx
             # always create blocks of 512x512 as that is native shift order 9 size

python/lsst/pipe/tasks/fits2hips/_low_order_hips.py

@@ -176,7 +176,7 @@ def run(self, hpx_container: Iterable[tuple[DeferredDatasetHandle, int]]) -> Str
                 size_counter += 1
 
                 # resample the image to a smaller grid and store it for the next order
-                zoomed = cv2.resize(hpx_next_array, (256, 256), interpolation=cv2.INTER_LANCZOS4)
+                zoomed = cv2.resize(hpx_next_array, (256, 256), interpolation=cv2.INTER_AREA)
 
                 hpx_next_container.append((zoomed, hpx_next_id))
             hpx_container = hpx_next_container