Merge pull request #213 from yyu1/optimization

Do not create array with ones, and then multiply value by nodata. In…
perrygeo · Feb 29, 2020 · 853a0c3 · 853a0c3
2 parents f967c1d + a913471
commit 853a0c3
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/rasterstats/io.py b/src/rasterstats/io.py
@@ -183,7 +183,8 @@ def boundless_array(arr, window, nodata, masked=False):
         window_shape = (wr_stop - wr_start, wc_stop - wc_start)
 
     # create an array of nodata values
-    out = np.ones(shape=window_shape) * nodata
+    out = np.empty(shape=window_shape)
+    out[:] = nodata
 
     # Fill with data where overlapping
     nr_start = olr_start - wr_start