I am trying to reclassify a continuous raster(Int32) using GDAL Python/NumPy, by giving values greater than a number (in this case 60 minutes) 0, and less than/equal to the number 1. It's working for the most part, but some pixels are not changing to 1 in spite of their value being equal to the number. In other cases pixels equal to the number are 1.
I think the culprit line is:
cum_cost_array_bin = np.where(cum_cost_array_cont < (num_hours*60), 1, 0)
I have also tried:
cum_cost_array_bin = np.where(cum_cost_array_cont <= (num_hours*60), 1, 0)
The image shows the output, with the blue pixels <= 60 minutes, and the red > 60 in the continuous. I've reclassified the raster using ArcMap, and the yellow polygon represents that output, which is correct. Any ideas as to why NumPy is not changing the values based on the conditional?