Re-organize coreg.py

.pre-commit-config.yaml

@@ -14,7 +14,7 @@ repos:
           hooks:
                 - id: codespell
                   args: [
-                    '--ignore-words-list', 'nd,alos',
+                    '--ignore-words-list', 'nd,alos,inout',
                     '--ignore-regex', '\bhist\b',
                     '--'
                   ]

docs/source/code/coregistration_plot_nuth_kaab.py

@@ -11,11 +11,11 @@
 inlier_mask = ~outlines_1990.create_mask(dem_2009)
 
 nuth_kaab = xdem.coreg.NuthKaab()
-nuth_kaab.fit(dem_2009.data, dem_1990.data, transform=dem_2009.transform, inlier_mask=inlier_mask)
-dem_coreg = nuth_kaab.apply(dem_1990.data, transform=dem_1990.transform)
+nuth_kaab.fit(dem_2009, dem_1990, inlier_mask=inlier_mask)
+dem_coreg = nuth_kaab.apply(dem_1990)
 
 ddem_pre = (dem_2009.data - dem_1990.data).filled(np.nan).squeeze()
-ddem_post = (dem_2009.data - dem_coreg).filled(np.nan).squeeze()
+ddem_post = (dem_2009.data - dem_coreg.data).filled(np.nan).squeeze()
 
 nmad_pre = xdem.spatialstats.nmad(ddem_pre[inlier_mask.squeeze()])
 nmad_post = xdem.spatialstats.nmad(ddem_post[inlier_mask.squeeze()])

tests/test_examples.py

@@ -55,7 +55,8 @@ def test_array_content(self, rst_and_truevals: tuple[Raster, NDArrayf]) -> None:
 
         assert values == pytest.approx(truevals)
 
-    @pytest.mark.parametrize("rst_and_truenodata", [(ref_dem, 0), (tba_dem, 0), (ddem, 2316)])  # type: ignore
+    # Note: Following PR #329, no gaps on DEM edges after coregistration
+    @pytest.mark.parametrize("rst_and_truenodata", [(ref_dem, 0), (tba_dem, 0), (ddem, 0)])  # type: ignore
     def test_array_nodata(self, rst_and_truenodata: tuple[Raster, int]) -> None:
         """Let's also check that the data arrays have always the same number of not finite values"""
 

tests/test_spatialstats.py

@@ -390,9 +390,9 @@ def test_sample_multirange_variogram_default(self) -> None:
 
         # Check the variogram output is consistent for a random state
         df = xdem.spatialstats.sample_empirical_variogram(values=self.diff, subsample=10, random_state=42)
-        assert df["exp"][15] == pytest.approx(23.517837524414062, abs=1e-3)
+        assert df["exp"][15] == pytest.approx(5.087792205810548, abs=1e-3)
         assert df["lags"][15] == pytest.approx(5120)
-        assert df["count"][15] == 2
+        assert df["count"][15] == 5
         # With a single run, no error can be estimated
         assert all(np.isnan(df.err_exp.values))
 
@@ -1180,7 +1180,7 @@ def test_patches_method_loop_quadrant(self) -> None:
         assert all(df.columns == ["nmad", "nb_indep_patches", "exact_areas", "areas"])
 
         # Check the sampling is fixed for a random state
-        assert df["nmad"][0] == pytest.approx(1.8663623135417342, abs=1e-3)
+        assert df["nmad"][0] == pytest.approx(1.8392861049330378, abs=1e-3)
         assert df["nb_indep_patches"][0] == 100
         assert df["exact_areas"][0] == pytest.approx(df["areas"][0], rel=0.2)