rasterstats

import os from rasterstats import zonal_stats, point_query # This is a simplified example. In a real scenario, you'd load actual geospatial files. # For demonstration, we'll simulate inputs or refer to a common pattern. # Assuming 'polygons.shp' and 'elevation.tif' exist in a 'data' directory for this example. # You would replace these paths with your actual data file paths. # Example for zonal_stats # stats = zonal_stats(os.environ.get('VECTOR_PATH', 'data/polygons.shp'), # os.environ.get('RASTER_PATH', 'data/elevation.tif'), # stats=['min', 'max', 'mean', 'median']) # print('Zonal Stats for first polygon:', stats[0]) # Example for point_query # point_geom = {'type': 'Point', 'coordinates': (245309.0, 1000064.0)} # Example coordinate # value = point_query(point_geom, os.environ.get('RASTER_PATH', 'data/elevation.tif')) # print('Point Query value:', value) print("To run quickstart, ensure you have 'polygons.shp' and 'elevation.tif' files.") print("Uncomment the examples in the code to execute them.") print("For full quickstart, refer to the official documentation.")