Automating Raster Data Extraction using Python

In this tutorial, you will learn how you can use a simple python script to extract data from a raster.

Extracting pixel values at point locations can easily be done using a dedicated GIS software package such as QGIS or ArcGIS, if you are doing it just for a limited number of points. However, if you were to do the same for a large number of points, this exercise can get quite tedious, repetitive and extremely challenging.

In this tutorial, we are discussing a way of doing this in just a matter of few seconds using a small Python script, regardless of how large the dataset is. To demonstrate this, we are making use of a point shapefile, which indicates a number of selected cities across the United States. Along with that, we are using a raster dataset which shows the average temperature variation across the country at a selected day in November 2021.

The two major Python libraries that are being incorporated in this tutorial are GeoPandas and Rasterio. With the help of these two libraries, we manage to extract the point coordinates in the form of Latitude and Longitude values at each point location (city), and extract the corresponding temperature value (pixel value) from the raster for each of the aforementioned points. Finally, we write the information to a Pandas DataFrame and export it as a separate CSV file using the Pandas library.

#spatialanalysis #geopandas #rasterio

[pacomermela]

Excellent. It would be interesting to show how to extract pixel values ​​from an ESRI Shapefile of polygons to create a dataframe. It would also be interesting to learn how to assign a value to each sample according to a field in the attribute table. Thanks for your videos!!!

[christodoulossazeides]

i was searching for 3 days to do that. Every one suggested a different packages as (gdal, pygmt...) and i was not able to extract the values. Your simple way solve my problem. thank you very match

[ibrahimalabdulsalam]

I should always start your videos by clicking the LIKE button.

[niampahamidou]

Thanks a lot for this video. You're the best!

[johnowusukonduah]

These video came as a live saver just when I needed it most. Thank you very much. But I want to ask if this same approach works similarly for multiband raster data? If its not so, please can you do a video to cover that? This will really help me a lot

[sonukumar-cyfo]

Thank you very much sir for your kind support. Your all material is very interesting. May I request you to make one video on object based classification

[nilpipra]

wish you could provide the sample files (raster and the points) used in this tutorials. We all learn a lot from you. You are the best!!

[syedrafsanali]

Thank you very much. It was very helpful.

[PhatPham-User]

Thanks so much for your great tutorial. I followed your step and I can extract the tiff file already but I still have a concern that I have many files from the year 1981 - up to now (the step is one day). I would like to get all data day by day but still cannot solve this problem, could you has any idea or could you make another tutorial to do that if possible. Thanks again for your kind.

[AsifMarazi]

Great video... really helped me get closer to my goal. However, now I have multiple raster files, I tried to loop over the raster files but failed to write the extracted values in a better way.

[abdallahabdelmajeed]

Thank you very much but how can we do the same for polygons not point and in time series data?

[SiddhZarathustra]

Great tutorial! Can you please share the source for download link or the temperature raster?

[joshritman]

Great content. Thanks. Where can we download the raster data (avg_temp_nov2021.tif) from?

[mrcactus]

hello, how could i extract values from a certain point from diferent rasters to table? thank you

[moumitaghorai]

I am struggling with installing rasterio. I can see it in the site-packages folder. But cannot import it, getting an error. Do you have any suggestions on this?

[MrABN]

Can we use same process in other operations in raster? Like solar panel monitor, thermal sensor raster?

[AbuZafor-yyqh]

AttributeError: 'numpy.float32' object has no attribute 'is_integer'
at the last segment this error is come.What is the solution ?

[notweak]

Put the readio statement outside for loop

[Jobsvisapro]

Sir can you please upload complete tutorial list about PyQgis

[mmukarrammuneer]

Can I done same task in Jupyter of NetCdf (Raster) file with any shaple file???