Merge pull request #10695 from rouault/stacit_single_feature

STACIT: support top-level Feature in addition to FeatureCollection
OSGeo · Sep 7, 2024 · f65642c · f65642c
2 parents 072db35 + 267a190
commit f65642c
Show file tree

Hide file tree

Showing 2 changed files with 28 additions and 3 deletions.
diff --git a/autotest/gdrivers/stacit.py b/autotest/gdrivers/stacit.py
@@ -361,3 +361,20 @@ def test_stacit_force_opening_no_match():
 
     drv = gdal.IdentifyDriverEx("data/byte.tif", allowed_drivers=["STACIT"])
     assert drv is None
+
+
+###############################################################################
+# Test opening a top-level Feature
+
+
+def test_stacit_single_feature(tmp_vsimem):
+
+    j = json.loads(open("data/stacit/test.json", "rb").read())
+    j = j["features"][0]
+
+    filename = str(tmp_vsimem / "feature.json")
+    with gdaltest.tempfile(filename, json.dumps(j)):
+        ds = gdal.Open(filename)
+        assert ds is not None
+        assert ds.RasterXSize == 20
+        assert ds.GetRasterBand(1).Checksum() == 4672
diff --git a/frmts/stacit/stacitdataset.cpp b/frmts/stacit/stacitdataset.cpp
@@ -852,11 +852,19 @@ bool STACITDataset::Open(GDALOpenInfo *poOpenInfo)
                 return false;
         }
         const auto oRoot = oDoc.GetRoot();
-        const auto oFeatures = oRoot.GetArray("features");
+        auto oFeatures = oRoot.GetArray("features");
         if (!oFeatures.IsValid())
         {
-            CPLError(CE_Failure, CPLE_AppDefined, "Missing features");
-            return false;
+            if (oRoot.GetString("type") == "Feature")
+            {
+                oFeatures = CPLJSONArray();
+                oFeatures.Add(oRoot);
+            }
+            else
+            {
+                CPLError(CE_Failure, CPLE_AppDefined, "Missing features");
+                return false;
+            }
         }
         for (const auto &oFeature : oFeatures)
         {