Take the 2-minute tour ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

I'm working with on a DEM (Geotiff) in Python GDAL, and I'd like to find out the spikes or pits in my terrain.
I have no idea how to proceed, does anyone have good algorithm to help me out?
Thanks,
eo

share|improve this question
1  
define 'spikes' and 'pits', and you'll find where to begin. –  nickves May 15 '13 at 10:11
add comment

1 Answer

GDAL doesn't have a ready-made utility for this that I know of. You would probably have to script something.

However, for filling sinks, GRASS has the r.fill.dir and r.terraflow functions. SAGA and some other FOSS raster GIS packages also have sink filling functions too.

Spikes pose a different problem and you'll need to consider a median filter and/or Gausian blur. See this thread for more discussion on spikes.

share|improve this answer
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.