search

esciencecenter-digital-skills / geospatial-python

Introduction to Geospatial Raster and Vector Data with Python
https://esciencecenter-digital-skills.github.io/geospatial-python/
Other
3 stars 0 forks source link

Ep 10 zonal stats #20

Closed cpranav93 closed 2 years ago

cpranav93 commented 2 years ago

This episode deals with calculating zonal statistics on raster data using vector defined zones. We use the cropped_field.shp vector data file to defined our crop zones to calculate the NDVI statistics on these zones. To use vector fields, they need to be rasterized apriori which requires a conversion of CRS and a cropping to the correct domain size. The exercise deals with calculating zonal statistics on the NDVI classified data from the previous episode.

The zonal stats are currently output for each crop code (gewascode) with no title to the zone. The lesson can be extended to extract, simplify, and add zone titles to the statistics if needed - I thought the lesson was already long and possibly complicated so I have opted to leave that out for now.

This PR addresses issue #6

cpranav93 commented 2 years ago

Hey Ou,

I have addressed your comments and added/modified a few lines. Give it a quick look through and let me know your thoughts?

cpranav93 commented 2 years ago

Feel free to leave further comments @SarahAlidoost and @rbavery on the content of this episode.