python gdal 重采样_用GDAL对影像重采样的一些要点

1、RasterIO只能使用最临近插值法

做影像金字塔时，开始用的GDALDataSet的RasterIO通过控制行列宽度来抽层，代码写起来很简洁，奇怪一直没见设置重采样参数的地方，急于实现功能，没有怎么深究，基本功能全都做好后发现分层质量很差。文档里边没找到资料介绍这些，查了很久在Maillist上发现RasterIO() level resampling should all be "nearest neighbour", and so

not introduce any new values.看来不换方法是不行了。

参考资源：http://lists.maptools.org/pipermail/gdal-dev/2005-May/005646.html

On 5/4/05, Sheykhet, Rostic wrote:

>Hello,

>

>I am trying to perform RasterIO on a band that references a color

table.

>When the source and destination image dimensions are different, I get

>back a band that has indexes into blank/invalid color table entries.

>I assume this is because of the decimation procedure. Is there a way

>to modify this behavior?

Rostic,

RasterIO() level resampling should all be "nearest