working on point extraction

gregory-halverson · Nov 5, 2024 · 1bfc631 · 1bfc631
1 parent b2a35e5
commit 1bfc631
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/rasters/raster.py b/rasters/raster.py
@@ -970,6 +970,9 @@ def geolocation(self) -> Raster:
         return self.contain(geometry=self.geometry.geolocation)
 
     def to_point(self, point: Point):
+        if not self.geometry.intersects(point):
+            return np.nan
+
         index = self.geometry.index_point(point)
         value = self.array[index]