Correct faulty input geometries

wri · Oct 12, 2021 · 5b4a436 · 5b4a436
1 parent 5f45510
commit 5b4a436
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/src/main/scala/org/globalforestwatch/features/SpatialFeatureDF.scala b/src/main/scala/org/globalforestwatch/features/SpatialFeatureDF.scala
@@ -48,9 +48,11 @@ object SpatialFeatureDF {
       FeatureDF(input, featureObj, filters, spark, delimiter)
     val emptyPolygonWKB = "0106000020E610000000000000"
 
+    // ST_PrecisionReduce may create invalid geometry if it contains a "sliver" that is below the precision threshold
+    // ST_Buffer(0) fixes these invalid geometries
     featureDF
       .selectExpr(
-        s"ST_PrecisionReduce(ST_GeomFromWKB(${wkbField}), 11) AS polyshape",
+        s"ST_Buffer(ST_PrecisionReduce(ST_GeomFromWKB(${wkbField}), 11), 0) AS polyshape",
         s"struct(${featureObj.featureIdExpr}) as featureId"
       )
       .where(s"${wkbField} != '${emptyPolygonWKB}'")

diff --git a/version.sbt b/version.sbt
@@ -1 +1 @@
-version in ThisBuild := "1.7.6"
+version in ThisBuild := "1.7.7"
Original file line number	Diff line number	Diff line change
		@@ -1 +1 @@
		version in ThisBuild := "1.7.6"
		version in ThisBuild := "1.7.7"