# Check that downsampled data covers the same lat/lon extents as original data, but is truly taking a subset of the locations _, orig_lons, orig_lats = sample_raster_data(tif_bytes_out, interval=1) # Plot the raster data and sampled points plt.scatter(orig_lons, orig_lats, color='blue', alpha=1, s=1) plt.scatter(longitudes, latitudes, color='red', s=1, alpha=1) plt.gca().set_xlim(left=orig_lons.min(), right=orig_lons.min()+0.002) plt.gca().set_ylim(bottom=orig_lats.min(), top=orig_lats.min()+0.08) plt.show() plt.figure() plt.scatter(orig_lons, orig_lats, color='blue', s=.001) plt.scatter(longitudes, latitudes, color='red', s=.001) plt.show()