Extracting Topographic Elevations using Python

[@wasumanawko4825]

I have no words to appreciate you. all of your video tutorials are very helpful for many researchers. Your presentation is also smart. please keep it up.

[@yagitgadhok5938]

Hi, spyder is showing a error of invalid projection:no database context specified can you please help me resolve it.

[@nitingoyal5515]

can you please tell, how to extract coordinates from image into csv file

[@adampax]

This was very helpful, thanks for taking the time to make such a clear tutorial :)

[@saradindudey5011]

Please help me in that, I really need this heatmap thing cause there are no such tutorial on that topic

[@dahmanigis4121]

Hi, when i run this # reading DEM dem = rasterio.open('DEM.tif') row, col = dem.index(longitude, latitude) dem_data = dem.read(1) print('The elvation of' + name + ':' + str(dem_data[row,col]) + 'meters') I have error: IndexError Traceback (most recent call last) Cell In [22], line 5 3 row, col = dem.index(longitude, latitude) 4 dem_data = dem.read(1) ----> 5 print('The elvation of' + name + ':' + str(dem_data[row,col]) + 'meters') IndexError: index -977727 is out of bounds for axis 0 with size 5715

[@rrex7948]

Good stuff, but you skipped over what is going on with dem.read(1)

[@waltersamuel3671]

Love your videos !! It has really helped with getting work done for my research .. Cheers

[@chuxinliu]

How do I find DEM file for, let's say, Florida?

[@geodeltalabs]

You can use USGS EarthExplorer. There is a tutorial that I have done showing how to download SRTM dem for any location :)

[@chuxinliu]

@@geodeltalabs Thank you! It is very helpful!

[@anujregmi4582]

is there a way to extract the entire area elevation instead of just points

[@geodeltalabs]

How do you intend to define the 'entire area elevation'? as long as you have a sufficient coverage of points to cover each pixel, that would basically give you all the different elevation values.

[@xudongfan6962]

Hi, thanks for your video! Could you explain how should we choose the unit for longitude and latitude? I put longitude and latitude in the 'feet' unit in the 15 line. But the row and col out of the boundary. Should we convert the lon and lat into a specific unit? If i type 'dem.crs', it outputs CRS.from_epsg(26917)

[@rcssiddharth5986]

Yes, both the CRS of the Shapefile and DEM raster have to be in same units. Moreover, we could reproject geodataframe simply using gdf = gdf.to_crs(3857). Then, use rioxarray to reproject (which is simpler than rasterio). import rioxarray as rxr -> dem_raster = rxr.open_rasterio('filepath_to_raster') -> dem_reproject = dem_raster.rio.reproject('EPSG:3857') -> dem_reproject.to_raster('output_file_path'). Now, you can use the above code to extraction the elevation perfectly.

[@saradindudey5011]

Amazing one , I need a help can you make a video how to create a heatmap using python using shape points and save it in tiff raster format? It will be really helpful. Scenario: suppose you have shape points of wildfires with “brightness“ column and you want to create heatmap of that column using that shapefile points. And you want to save it in tiff

[@MdAthick]

Sir, I am a PhD student. I want you python for processing geospatial data, could you please suggest me some books, materials and links for Python for geospatial analysis.

[@geodeltalabs]

Hi, there are good resources by Packt Publishing. www.packtpub.com/ I recommend you check some out there.

[@MdAthick]

@@geodeltalabs Thanks a lot for your suggestion. Could you please tell me the keywords to search in packtpup ?

[@yaarobmohammad1020]

Thank you so much