gdalwrap tif GDALWarp

本章来源：

顺序重新调整,增加GDAL算法介绍，本章是进阶章节，阅读前请务必先熟悉 GDAL数据模型 与 栅格数据读写 章节。

GDAL Warp API（在文件 gdalwarper.h 中定义）是一个高效的进行图像变换的接口。主要由几何变换函数（ GDALTransformerFunc ），多种图像重采样方式，掩码操作选项等组成。这个接口可以对很大的图像进行处理。

使用步骤如下：

下面先介绍 GDALWarpOptions 和 GDALTransformerFunc ，然后对完整流程进行解析。

GDALWarpOptions 是对图像变换进行设置的类，成员见下：

GDALWarpOptions结构体中包含了很多参数来对变换进行设置。下面对一些比较重要的进行列举说明：

这个是一个字符串列表，用来设置图像变换过程中的一些选项，样式为 NAME=VALUE , 可以使用 CSLSetNameValue() 设置。现在支持的一些键值对设置为：

通常在计算一个数据的部分区域时，变换样区在输出区域每个边界上转换21个点到源文件中，然后在源图中计算一个足够大的窗口。这样做是为了使转换更加高效，但某些情况下，这种做法计算的窗口过小，甚至丢失大量的区域。特别是对于非线性或者是翻转的转换。例如对极地点附近的区域，从WGS84转到Polar Stereographic投影转换，或者某些根本无法进行投影转换的数据。对于这种情况，GDAL提供以下键值对设置来防止窗口计算错误：

设置GDALWarpOperation在处理图像中使用的最大内存数。单位为比特，默认值比较保守（比较小），可以根据自己的内存大小来进行调整这个值。增加这个值可以帮助提高程序的运行效率，但是需要注意内存的大小。这个大小、GDAL的缓存大小，还有你的应用程序以及系统所需要的内存加起来要小于系统的内存，否则会导致错误。一般来说，比如一个内存为256MB的系统，这个值最少设置为64MB比较合理。注意，这个值不包括GDAL读取数据使用的内存。

重采样的算法选择：

计算时用的数据类型，GDT_Unknown 为程序自动判断

源数据集

目标数据集

所需转换的波段数，0为所有波段

int数组，源图像中所需处理的波段，从1开始

int数组，输出的波段

源数据中的Alpha波段，0为禁止

目标数据中的Alpha波段，0为禁止

源数据中 nodata 的实部，NULL表示没有

源数据中 nodata 的虚部，NULL表示没有

目标数据中 nodata 的实部，NULL表示没有

目标数据中 nodata 的虚部，NULL表示没有

进度条显示函数，NULL表示没有

传给pfnProgress的回调参数

几何变换函数

几何变换函数参数

裁切线，OGRPolygonH 类型

切割线宽容度，默认为零

GDALTransformerFunc 是一个函数签名，用来计算空间点几何变换，也是 GDALWarp 的核心部分，定义在GDAL算法文件中。 GDAL 中可用的算法有：

与 GDALTransformerFunc 对应，默认有几个函数生成几何变换函数参数：

这部分看起来比较复杂，实际使用时，比较简便。只需要在 GDALWarpOptions->pfnTransformer 中指定算法名称， GDALWarpOptions->pTransformerArg 中指定对应的算法参数函数即可，如下代码所示：

每个Transformer的所需参数不同，请参见 GDAL Algorithms C API ，近似计算有些特殊，在使用Transformer时要传递几何变换函数，参考 性能优化 部分最后一条。

下面介绍使用 GDALWarp 提供的接口，实现几何变换的完整流程：

几何变换后输出文件的范围和投影仿射参数都将发生变化，输出的文件大小可能是已知，也可能是未知的。因此创建输出文件方式也有两种：

重采样创建文件示例：

投影变换创建文件示例：

选项设置参见 GDALWarpOptions选项介绍 ，设置一般需要设置输入输出数据集，需要操作的波段，几何变换函数和参数，进度条，优化参数这几项。示例如下：

执行较为简单，新建 GDALWarpOperation ,调用 Initialize 函数，然后调用 ChunkAndWarpMulti 或者 ChunkAndWarpImage ，代码一般如下：

完成后，需要释放变换参数、删除数据集。代码如下：

完整的重采样代码如下：

下面几点可以在使用变换API的时候提高处理效率。