接上面，可以使用exe工具实现同样的功能。

这个工具详情如下：

gdalwarp [--help-general] [--formats]
    [-s_srs srs_def] [-t_srs srs_def] [-ct string] [-to "NAME=VALUE"]* [-novshiftgrid]
    [-order n | -tps | -rpc | -geoloc] [-et err_threshold]
    [-refine_gcps tolerance [minimum_gcps]]
    [-te xmin ymin xmax ymax] [-te_srs srs_def]
    [-tr xres yres] [-tap] [-ts width height]
    [-ovr level|AUTO|AUTO-n|NONE] [-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16]
    [-srcnodata "value [value...]"] [-dstnodata "value [value...]"]
    [-srcalpha|-nosrcalpha] [-dstalpha]
    [-r resampling_method] [-wm memory_in_mb] [-multi] [-q]
    [-cutline datasource] [-cl layer] [-cwhere expression]
    [-csql statement] [-cblend dist_in_pixels] [-crop_to_cutline]
    [-if format]* [-of format] [-co "NAME=VALUE"]* [-overwrite]
    [-nomd] [-cvmd meta_conflict_value] [-setci] [-oo NAME=VALUE]*
    [-doo NAME=VALUE]*
    srcfile* dstfile

参数的具体说明见如下链接：

https://www.osgeo.cn/gdal/programs/gdalwarp.html#gdalwarp

数据情况和上篇一样。

这里直接上例子：

gdalwarp.exe -t_srs '+proj=aea +lat_1=15 +lat_2=65 +lat_0=30 +lon_0=95 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs'E:\testdata\CZIOrigin.tiff O:\CZIProject.tif

空间参考采用的是proj.4的形式，当然也可以使用EPSG格式。

和这个exe工具对应的有代码中的调用，可以使用GDALWarp()函数，这个函数在gdal_utils.h头文件下。

代码：

#include <iostream>
#include <gdal_priv.h>
#include <gdal_utils.h>

int main()
{
	GDALAllRegister();

	GDALDatasetH srcDs = GDALOpen("E:\\testdata\\CZIOrigin.tiff", GA_ReadOnly);

	char** papszArgv = NULL;
	papszArgv = CSLAddString(papszArgv, "-s_srs");
	papszArgv = CSLAddString(papszArgv, "+proj=longlat +datum=WGS84 +no_defs");
	papszArgv = CSLAddString(papszArgv, "-t_srs");
	papszArgv = CSLAddString(papszArgv, "+proj=aea +lat_1=15 +lat_2=65 +lat_0=30 +lon_0=95 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs");

	GDALWarpAppOptions* warpAppOptions = GDALWarpAppOptionsNew(papszArgv, NULL);

	GDALDatasetH projectData = GDALWarp("O:\\CZIReprojext.tif", NULL, 1, &srcDs, warpAppOptions, NULL);

	GDALClose(srcDs);
    CSLDestroy(papszArgv);
	GDALWarpAppOptionsFree(warpAppOptions);
	GDALClose(projectData);

	system("pause");
	return 0;
}

这两个的效果和之前是一样的。