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 work on Landsat 7 slc-off images and I would like to use the gdal_fillnodata.py command to fill-in the missing pixels. The problem that I have is that the gdal_fillnodata.py uses the "inverse distance weighting" (IDW) interpolation technique and I want to use the kriging geostatistical technique.

Do you have any idea how can I do it?

Thank you!

share|improve this question
add comment

1 Answer

AFAIK, you can't. There are no other algorithms defined in GDALFillNodata():

http://gdal.org/gdal__alg_8h.html#a0a079afef61968f224b159d48423d1c6

If the note mentioned in the docs applies to you:

This algorithm is generally suitable for interpolating missing regions of fairly continuously varying rasters (such as elevation models for instance). It is also suitable for filling small holes and cracks in more irregularly varying images (like airphotos). It is generally not so great for interpolating a raster from sparse point data - see the algorithms defined in gdal_grid.h for that case.

you are probably fine with GDALFillNodata().

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.