Use geomad instead of hdstats.

Author: SpacemanPaul

docker/requirements.txt

@@ -1,11 +1,10 @@
 
---extra-index-url https://packages.dea.ga.gov.au/
 datacube-ows>=1.9
 datacube[performance,s3]>=1.9.5
 eodatasets3>1.9
-hdstats==0.1.8.post1
+geomad
 numexpr>=2.11
-odc-algo>=1.0.1
+odc-algo>=1.1.0
 odc-apps-cloud>=0.2.2
 # For testing
 odc-apps-dc-tools>=0.2.12

tests/test_gm_ls.py

@@ -82,7 +82,7 @@ def dataset():
 
 
 def test_native_transform(dataset):
-    _ = pytest.importorskip("hdstats")
+    _ = pytest.importorskip("geomad")
 
     dataset = dataset.copy()
     stats_gmls = StatsGMLS(nodata_classes=(0,))
@@ -104,7 +104,7 @@ def test_native_transform(dataset):
 
 
 def test_result_bands_to_match_inputs(dataset):
-    _ = pytest.importorskip("hdstats")
+    _ = pytest.importorskip("geomad")
     mask_filters = {
         "cloud": [("closing", 2), ("dilation", 1)],
         "shadow": [("closing", 0), ("dilation", 5)],
@@ -130,7 +130,7 @@ def test_result_bands_to_match_inputs(dataset):
 
 
 def test_result_aux_bands_to_match_inputs(dataset):
-    _ = pytest.importorskip("hdstats")
+    _ = pytest.importorskip("geomad")
     mask_filters = {
         "cloud": [("closing", 2), ("dilation", 1)],
         "shadow": [("closing", 0), ("dilation", 5)],
@@ -165,7 +165,7 @@ def test_result_aux_bands_to_match_inputs(dataset):
 
 
 def test_resampling(dataset):
-    _ = pytest.importorskip("hdstats")
+    _ = pytest.importorskip("geomad")
     mask_filters = {
         "cloud": [("closing", 2), ("dilation", 1)],
         "shadow": [("closing", 0), ("dilation", 5)],

tests/test_gm_ls_bitmask.py

@@ -109,7 +109,7 @@ def test_fuser(dataset):
 
 
 def test_reduce(dataset):
-    _ = pytest.importorskip("hdstats")
+    _ = pytest.importorskip("geomad")
     gm = StatsGMLSBitmask(
         bands=["band_red"], offset=-0.2, scale=0.00975, output_scale=100
     )
@@ -140,7 +140,7 @@ def test_reduce(dataset):
 
 
 def test_reduce_with_filters(dataset):
-    _ = pytest.importorskip("hdstats")
+    _ = pytest.importorskip("geomad")
     mask_filters = [("closing", 2), ("dilation", 1)]
     gm = StatsGMLSBitmask(
         bands=["band_red"],
@@ -173,7 +173,7 @@ def test_reduce_with_filters(dataset):
 
 
 def test_aux_result_bands_to_match_inputs(dataset):
-    _ = pytest.importorskip("hdstats")
+    _ = pytest.importorskip("geomad")
     aux_names = dict(smad="SMAD", emad="EMAD", bcmad="BCMAD", count="COUNT")
     gm = StatsGMLSBitmask(
         bands=["band_red"],