1

I'm trying to get slope data from the SRTM (Shuttle Radar Topography Mission) elevation dataset which provides many .hgt files. I've parsed the hgt files to get meters above sea level. When I'm calculating slope in degrees, my data tends to be quite noisy and seems a bit off.

Ex. the following slope maps (in degrees) are from a very flat area and mountainous area. The first one is really noisy and the second one seems uncalibrated.

Flat: https://i.stack.imgur.com/HGNJ8.png

Mountainous: https://i.stack.imgur.com/K02Td.png

I was told that this was because I'm calculating slope directly from the hgt files, which are in WGS84, and that I need to convert to a meter-based CRS before generating a slope map. I don't really know what CRS to use however, or how to convert between them.

I looked at some of the code on this website for the documentation for geopandas: https://geopandas.org/projections.html#re-projecting, and it seems like the re-projection example code could be helpful, but I don't know how to feed in hgt files to geopandas or convert between CRSs.

Shrey Joshi
  • 1,066
  • 8
  • 25
  • Noisy: it is normal: if you read the documentation, you see that it fails with water, and with snow in mountain (and a lot more). You see "volcano craters" instead of peaks. In any case, I would calculate with original file, and just convert the result in your output CRS. How do you do your slope calculations? with geotiff tools? There you can also specify the CRS – Giacomo Catenazzi Jan 08 '21 at 10:25
  • I'm doing slope calculations with a library called richdem in python: https://richdem.readthedocs.io/en/latest/ – Shrey Joshi Jan 09 '21 at 02:12

0 Answers0