Calculating Zonal Statistics of a Raster using Python

This tutorial explains the basics of calculating the zonal statistics such as average, max, min values corresponding to a defined spatial area of a raster using a polygon vector layer using Python, with the aid of GeoPandas, rasterio and rasterstats libraries. Additionally, you will learn how to plot rasters and vectors on the same layer, how to use matplotlib subplots to plot maps and histograms together in the same canvas, as well as to extract the summary of the data into a Pandas DataFrame using the Pandas library.
Download exercise data - drive.google.com/drive/folder...
⭐NOTE:
‘Kite’ is an excellent tool that can help you boost your productivity when it comes to programming with its free Artificial Intelligence (AI) powered coding assistant. Kite is a plugin which provides smart completions and documentation as you’re typing-in your code, and its compatible with most of the well-known IDEs such as PyCharm, Spyder and VS Code. I absolutely love it, and recommend you guys to try it out as well.
www.kite.com/get-kite/?...

[@olssonfilip]

Just wanted to thank you for these vids, wish I had them during my time at Uni but now as a GIS professional I’m feeling even more inspired by your work! Looking forward to the next one : )

[@geodeltalabs]

Glad you like them!

[@kaustavdas5600]

@geodeltalabs: In the same code, if I want to extract the 95th percentile value, what modification in the code do I need to do ?

[@moisengwa9904]

Very Good job GeoDelta Labs. I don't have enough words to thank you. Please please kindly post the part 2 tutorial Calculating Zonal Statistics of a Raster using Python. Looking forward to that:

[@muhammadanjum524]

Dear Sir , Thank you soooooooooooooo much . You are sooo nice and kind . You are sooo expert in GIS . I hope you will continue your help .Its really very very helpful for me . I actually want to do that one . Sir, you solved my problem . Once again THANK YOU soo much .

[@abramswee]

Thank you very much for sharing this with us! much appreciated.

[@rafaeltiara9092]

Great tutorial. Thank you!

[@mulugetaworkuayele6412]

I am too late to learn from this interesting and very useful tutorial. Thank you so much!

[@geodeltalabs]

Better late than never :)

[@amitblizer4567]

Thank you. This was very helpful for me in my thesis🙏

[@jeffz7310]

I didn't need to go to Uni because of these videos, thx

[@yaarobmohammad1020]

You are the best

[@mohammadbelalalzoubi8102]

Thank you. It was a great video.

[@geodeltalabs]

You are welcome :)

[@shaileshshrestha9725]

good job bro...more geopandas pls.. one thing you have a good voice

[@DileepGuptagoodman]

Very nice video . I like sir . 👍👍👍

[@geodeltalabs]

Thanks!

[@gmjakirullahnooruddin1734]

Dear @GeoDelta Labs, Thanks for this great tutorial. Very much appreciable. Waiting for part two with multiple raster data. I have a couple questions. Do I have to project both my vector and raster dataset? If so how could I project multiple raster datasets at once? and the second question is while plotting/showing the raster data how could I show also raster scale right-side on my raster plot? Thanks in advance.

[@hugonicolaubarbosadegusmao7824]

The tutorial was very good, I learned a lot! A question, where did you download this data? I wanted to download to Brazil.

[@aakashbhatt3303]

Hi @Geo Data Labs, Great Tutorial thanks for this And please can you make tutorials on the liabrary raterio

[@australoitaliano]

@GeoDelta Labs, thank you very much for this fantastic tutorial. Just one question: how can I save the "mean" dataset in the shapefile? After I have performed the zonal statistic, I have open the shapefile with QGis but the "mean" column was not present in the attribute table. Thanks

[@muhammadanjum524]

Dear Sir, I am also trying to use QGIs and also python , but I could not my results . Can we discuss Please ?

[@raquelzafrir2752]

Hi GeoDelta thanks for the video. I have a question. Did you check the values where your raster has no data or zero values? I am getting negative values and I do not know how to correct them. Thanks!

[@umairrasool9580]

Very great and clear work you have done. i have a question please, i am working on my raster dataset for prediction like ANN, RF and CNN, i have converted the rasters into numeric and then train and test the data and got very good accuracy. Now i need to convert my test data into raster again as final prediction map but i don't know how to do this, please guide me thanks.

[@sonalisharma3421]

Hello Sir, can you please tell me how to perform zonal statistics on each band of raster stack. Thank you very much for making such useful videos.

[@manuvideos3979]

Hello, Thanks for this great tutorial. I have a question, while calculating the average, are taking in account the nodata pixels? If yes, is there a way to exclude nodata pixels while calculating average? Thanks in advance for your answer.

[@yaarobmohammad1020]

Can you please make more videos How to extract and compute volumes from raster How to make hydrogeological analysis And so on

[@geodeltalabs]

Hi, do you mean to just to do a one time computation of the volume from a raster? If you do not have a specific reason to use scripting to do that, I believe you can do this kind of processes quite easily using qgis

[@yaarobmohammad1020]

GeoDelta Labs I think that coding always better than do it with software ; because u can use you’re code always with different situations and faster Thank you so much

[@kaustavdas5600]

Thanks a lot for this! ...just one question, if I want to extract the 95th percentile at any of the districts, wherever the raster is present. What modification in the code do I need to do. Anyone's help is really appreciated

[@DileepGuptagoodman]

Please make the video for mapping of CSV data in geotiff. CSV have lat lon and values

[@JesusRomero-zw8yw]

Thanks for the explanation, but if I want to extract the areal average of each day (that is, the average area of ​​a polygon) from a netcdf file, for example precipitation or temperature, maybe some idea

[@geodeltalabs]

Hi, for this I suggest you follow the tutorials that I have done on handling netcdf files using python. You might be able to modify the code in such a way that you would be able to extract the values into a separate set based on their geographical position, and perform the statistical operations on those extracted values.

[@shakhawathossainshohel5159]

Thanks a lot for the video, it was great. But I need to learn zonal statistics for multiple netcdf file like daily TRMM. Would you please help me?

[@geodeltalabs]

Better try to convert to raster first.

[@rajuahmmad9923]

I have converted APHRODITE precipitation data from netcdf to tif file. After that, i followed your tutorial and calculated mean value for each subbasin within my watershed shapefile, But i got NaN value in three of my subbasin. Could you please point out what might be the possible cause of this problem. Thanks in advance.

[@slimieboy6264]

having the same issue with a tif converted from a netcdf. idk what is wrong tbh

[@malaymondal9953]

From which site we can get raster weather data, like wind speed, rainfall, temp etc.

[@geodeltalabs]

You can try Aphrodite's portal for rainfall and temp in netCDF format. Windspeed, not so sure.

[@andresfiallos7112]

Hello GeoDeltaLabs, thanks a lot for this great video. I have a question and was hoping someone here can help me. When I run this line "avg_rf = rasterstats.zonal_stats(districts, rainfall_array, affine = affine, stats = ['mean'], geojson_out = True)" I get the following error: "MemoryError: Unable to allocate 3.95 TiB for an array with shape (620193, 874319) and data type float64" I have been trying to resolve this issue but I havent been able to do so. I would really appreciate the help!

[@eddisonjose]

Hola muchas gracias, eso se puede hacer en pycharm? Soy neófito en programación, pero que opinas del programa pycharm? Saludos

[@geodeltalabs]

Hello Eddison, pycharm is a good IDE too, however I have always worked with Spyder IDE, as the interface is simple, but sufficient enough to get the job done. You can indeed use pycharm to do all what we are doing here

[@eddisonjose]

@@geodeltalabs hola muchas gracias por su respuesta y amabilidad

[@eddisonjose]

GeoDelta Labs Hola una pregunta mas: Que opinas de hacer estos procesos en jupyter notebook? es factible? se pueden hacer allí?🙂

[@jitutiwari7813]

Hi, first of all Thank you for this Tutorial. I am a GIS and remote Sensing student and also want to Learn python for geospatial task can you please suggest me how could i learn python. or which books or tutorial can help me. waiting for your response. and please can you tell me how to calculate zonal on bigTiff.

[@geodeltalabs]

You are welcome :)

[@sonukumar-cy1fo]

Sir, I got this error :MemoryError: Unable to allocate 21.4 TiB for an array with shape (2028294, 1449644) and data type float64. please help me

[@daanidev]

Hi, Can you please tell me how you prepare tif files? i have tif files and shape file but i am unable to map correctly just as you did in this video. Can you help me with it? Thanks

[@geodeltalabs]

Hi, did you check whether both the raster and the vector files are having the same coordinate reference?

[@daanidev]

@@geodeltalabs Hi my shape file and raster file has same crs can you please help me to show shape file over raster in python ?

[@daanidev]

@@geodeltalabs thanks issue solved (Y)

[@reshmaa8818]

How to install rasterio and rasterstats?? Can you provide your video link below

[@muhammadanjum524]

Dear Sir, this is really really an excellent video for zonal statistics . But I followed your method but could not get results . I need to discuss . Sir, Please help me .This is my very very humble request . Please . I need to discuss with you . Please let me know how and where can i discuss with you ? Please my most humble request . Please reply positively as soon as possible . Please

[@geodeltalabs]

Hi, kindly send your query to our email. Thank you.

[@muhammadanjum524]

@@geodeltalabs Sir, I have raster , I did calculate average monthly values by using QGIS raster calculator and now I want to calculate Zonal statistics . How can i send you my raster and shape file . Please help me . Sir, its my humble request that how can i send you ? Please

[@wongamasiza7877]

I wish I could do all this in R. ..I can't install some of these packages as I am not administrator. I have to log calls to our service desk every time I need to install a package. When I try installing them on my laptop I get errors that I must provide a path to GDAL.. Last time I had to use anything related to GDAL it gave me hell. Its easier to install all these packages in R... my goodness, I have to migrate to python once. I am going to die