Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

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?

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

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

Your Answer


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.