Closed arevill2 closed 3 years ago
@arevill2, this looks great. I have two primary points of feedback.
get_centroids_XY
, I recommend adding a cell where we create a single rectangular data container (say for Edinburgh) and then query the ('bands', '1') field to show how one gets data from it. Then, it would be good to show querying some spatial fields like 'x', 'y', and 'radius'. This is also a good place to show the 'pixels' unit conversion. Then, I would show the projection image of the Edinburgh rectangle here instead of after the timing. This should help set the stage for then showing the timing to show that none of this takes much more time than the rasterio way.Minor point: in the timing cell, you can replace for example point_center[0:2,0]
with point_center[:,0]
. Just a little cleaner.
Updates from Dan's feedback.