added calculate zonal statistics