当前位置：首页 > 代码相关 > 正文内容

GDAL使用（四）

admin5年前 (2021-02-12)代码相关4367

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

这个工具详情如下：

```
<pre class="prism-highlight prism-language-powershell">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:\\estdataCZIOrigin.tiff O:\\CZIProject.tif

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

![image.png](https://lovedm.club/zb_users/upload/2021/02/202102121613109526206041.png)

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

代码：

```
<pre class="prism-highlight prism-language-cpp">#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;
}
```

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