Avoids empty geometries for parts of Points

Author: bethard

src/normit/geo/ops.py

@@ -63,6 +63,10 @@ def of(self,
         return result.intersection(sector)
 
     def part_of(self, geometry: shapely.geometry.base.BaseGeometry):
+        # points cannot be divided into parts
+        if isinstance(geometry, shapely.Point | shapely.MultiPoint):
+            return geometry
+        # calculate diameter and centroid
         d = shapely.minimum_bounding_radius(geometry) * 2
         c = geometry.centroid
         # create point at azimuth 0 (North) and rotate (negative = clockwise)

test/test_geo_ops.py

@@ -31,6 +31,10 @@ def test_inner_regions(georeader: GeoJsonDirReader):
     assert NorthEast.part_of(india).intersection(gorakhpur).area > 0
     assert SouthWest.part_of(india).intersection(gorakhpur).area == 0
 
+    # SouthWest part of a point is the point
+    kintbury = georeader.read(21517801)
+    assert SouthWest.part_of(kintbury) == kintbury
+
 
 def test_outer_regions(georeader: GeoJsonDirReader):
     az = georeader.read(162018)