Automating Raster Data Extraction using Python

[@pacomermela6497]

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!!!

[@christodoulossazeides8988]

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

[@geodeltalabs]

Glad I could help

[@ibrahimalabdulsalam1005]

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

[@niampahamidou8931]

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

[@johnowusukonduah2305]

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-cy1fo]

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

[@MrABN143]

Excellent brother 💐

[@sanjayawanninayake3819]

Good one...!!!👍

[@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!!

[@geodeltalabs]

Thanks. Will include them wherever relevant in the upcoming videos :)

[@nilpipra]

@@geodeltalabs : is there still a chance to upload this two files (US temp raster and station shape) used in this exercise?. I am doing something similar to that and it will be a great help if I have the exact file to practices,. Thanks in advance!

[@syedrafsanali6522]

Thank you very much. It was very helpful.

[@geodeltalabs]

Glad it helped!

[@syedrafsanali6522]

@@geodeltalabs Could you please check your email? I have sent something in detail. It would be a great help for me. Thanks again.

[@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.

[@abdallahabdelmajeed8120]

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?

[@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.

[@geodeltalabs]

Hi, you may follow this video and you will be able to get an insight on how to iterate through multiple raster files kzbin.info/www/bejne/labEh6WJmrF5o6M

[@joshritman2480]

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

[@not4weak]

Put the readio statement outside for loop

[@Jobsvisapro7030]

Sir can you please upload complete tutorial list about PyQgis

[@mrcactus1984]

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

[@AayushiSaini-l3g]

can you provide sample data that you have used??

[@moumitaghorai2040]

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?

[@MrABN143]

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

[@AbuZafor-yy5qh]

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

[@mmukarrammuneer8715]

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

[@geodeltalabs]

Absolutely

[@henryjayr]

can you show tis in arcMap?

[@geodeltalabs]

In arcmap this can be done using one of the inbuilt tools. Search for the tool 'Extract values to point', and it basically does the same thing

[@anitat9727]

thanks

[@geodeltalabs]

You're welcome!

[@manasranjanpanda9859]

Indeed it's an excellent tutorial but I think if we do have the shapefile of locations (includes coordinates) then just by using the extract raster values to points tool in ArcGIS, we can easily get values from n number of points that we need.

[@matthese9986]

But you have to pay for it

[@manasranjanpanda9859]

@@matthese9986 yeah I agree 😊

[@isabelahorta3063]

Hi! Thx for the tutorial! I got an error in the "print" line: --------------------------------------------------------------------------- IndexError Traceback (most recent call last) in () 11 rowIndex, colIndex = dem.index(longitude, latitude) 12 ---> 13 print(str(code) + ': ' + str(dem_data[rowIndex, colIndex])) IndexError: index 752266 is out of bounds for axis 0 with size 17336 Do you know how to solve it?

[@isabelahorta3063]

Hiiii! I already fixed it! I reprojected the shape and the raster to the same projection system. Thx.

[@geodeltalabs]

Great! glad you were able to fix it! :)