Corrected map_pixel

Corrects map_pixel, while making union and align work correctly
ozak · May 28, 2020 · c5416d8 · c5416d8
1 parent 90c443b
commit c5416d8
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/georasters/georasters.py b/georasters/georasters.py
@@ -99,8 +99,8 @@ def map_pixel(point_x, point_y, cellx, celly, xmin, ymax):
     '''
     point_x = np.asarray(point_x)
     point_y = np.asarray(point_y)
-    col = np.floor((point_x - xmin) / cellx).astype(int)
-    row = np.floor((point_y - ymax) / celly).astype(int)
+    col = np.round((point_x - xmin) / cellx).astype(int)
+    row = np.round((point_y - ymax) / celly).astype(int)
     return row, col
 
 def map_pixel_inv(row, col, cellx, celly, xmin, ymax):