CI fails, testing nhanced

xcube-dev · Jan 15, 2025 · c72acb4 · c72acb4
1 parent 987c30c
commit c72acb4
Showing 1 changed file with 42 additions and 21 deletions.
diff --git a/test/core/resampling/test_spatial.py b/test/core/resampling/test_spatial.py
@@ -89,58 +89,79 @@ def test_rectify_and_upscale_dataset(self):
     def test_reproject_utm(self):
         source_ds = self.new_5x5_dataset_regular_utm()
 
-        # similar resolution
+        # test projected CRS similar resolution
         target_gm = GridMapping.regular(
-            size=(5, 5), xy_min=(9.9886, 53.5499), xy_res=0.0006, crs=CRS_WGS84
+            size=(5, 5), xy_min=(4320080, 3382480), xy_res=80, crs="epsg:3035"
         )
         target_ds = resample_in_space(source_ds, target_gm=target_gm)
         np.testing.assert_almost_equal(
             target_ds.band_1.values,
             np.array(
                 [
-                    [7, 8, 8, 8, 9],
-                    [12, 13, 13, 13, 14],
-                    [12, 13, 13, 13, 14],
-                    [17, 18, 18, 18, 19],
-                    [22, 23, 23, 23, 24],
+                    [1, 1, 2, 3, 4],
+                    [6, 6, 7, 8, 9],
+                    [11, 12, 12, 13, 14],
+                    [16, 17, 17, 18, 19],
+                    [21, 17, 17, 18, 19],
                 ],
                 dtype=target_ds.band_1.dtype,
             ),
         )
 
-        # 1/2 resolution, test if subset calculation works as expected
+        # test projected CRS finer resolution
+        # test if subset calculation works as expected
         target_gm = GridMapping.regular(
-            size=(5, 5), xy_min=(9.9886, 53.5499), xy_res=0.0003, crs=CRS_WGS84
+            size=(5, 5), xy_min=(4320080, 3382480), xy_res=20, crs="epsg:3035"
         )
         target_ds = resample_in_space(source_ds, target_gm=target_gm)
         np.testing.assert_almost_equal(
             target_ds.band_1.values,
             np.array(
                 [
-                    [12, 12, 12, 13, 13],
-                    [17, 17, 17, 18, 18],
-                    [17, 17, 17, 18, 18],
-                    [22, 17, 17, 18, 18],
-                    [22, 22, 22, 23, 23],
+                    [15, 16, 16, 16, 16],
+                    [15, 16, 16, 16, 16],
+                    [15, 16, 16, 16, 16],
+                    [20, 21, 21, 21, 21],
+                    [20, 21, 21, 21, 21],
                 ],
                 dtype=target_ds.band_1.dtype,
             ),
         )
 
-        # test projected CRS
+        # test geographic CRS with similar resolution
         target_gm = GridMapping.regular(
-            size=(5, 5), xy_min=(4320080, 3382480), xy_res=80, crs="epsg:3035"
+            size=(5, 5), xy_min=(9.9886, 53.5499), xy_res=0.0006, crs=CRS_WGS84
         )
         target_ds = resample_in_space(source_ds, target_gm=target_gm)
         np.testing.assert_almost_equal(
             target_ds.band_1.values,
             np.array(
                 [
-                    [1, 1, 2, 3, 4],
-                    [6, 6, 7, 8, 9],
-                    [11, 12, 12, 13, 14],
-                    [16, 17, 17, 18, 19],
-                    [21, 17, 17, 18, 19],
+                    [7, 8, 8, 8, 9],
+                    [12, 13, 13, 13, 14],
+                    [12, 13, 13, 13, 14],
+                    [17, 18, 18, 18, 19],
+                    [22, 23, 23, 23, 24],
+                ],
+                dtype=target_ds.band_1.dtype,
+            ),
+        )
+
+        # test geographic CRS with 1/2 resolution
+        # test if subset calculation works as expected
+        target_gm = GridMapping.regular(
+            size=(5, 5), xy_min=(9.9886, 53.5499), xy_res=0.0003, crs=CRS_WGS84
+        )
+        target_ds = resample_in_space(source_ds, target_gm=target_gm)
+        np.testing.assert_almost_equal(
+            target_ds.band_1.values,
+            np.array(
+                [
+                    [12, 12, 12, 13, 13],
+                    [17, 17, 17, 18, 18],
+                    [17, 17, 17, 18, 18],
+                    [22, 17, 17, 18, 18],
+                    [22, 22, 22, 23, 23],
                 ],
                 dtype=target_ds.band_1.dtype,
             ),