Add info to use rasterstats to calculate zonal statistics