Hello all,

I'm trying to devise an efficient way to fill voids in SRTM, using
neighboring tiles when there is edge nodata.  In the end, I want to have
tiles matching the input, just lacking any voids.

My initial thought is to write a simple gdal python script to determine
which tiles have edge nodata and only process those.  I'll then write a
script to:

 - determine which tiles are neighbors
 - merge those into a single temp dem
 - set the region to the desired tile
 - r.fillnull the temp file, creating an output with the correct extents

This should work, though I'm thinking there must be a more efficient way.
Anyone have suggestions?


grass-user mailing list

Reply via email to