c++ gdal定义投影

在 C++ 中使用 GDAL 定义投影，可以按照以下步骤进行：

1. 引入 GDAL 头文件：

#include "gdal_priv.h"

2. 创建一个 GDAL 数据集对象：

GDALDataset *poDataset = (GDALDataset *) GDALOpen("your_file.tif", GA_ReadOnly);

3. 获取数据集的投影信息：

const char* pszProjection = poDataset->GetProjectionRef();

4. 根据获取到的投影信息创建一个投影对象：

OGRSpatialReference oSRS;
oSRS.importFromWkt(&pszProjection);

5. 定义新的投影信息：

const char* pszNewProjection = "PROJCS[\"WGS 84 / UTM zone 51N\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",0],PARAMETER[\"central_meridian\",117],PARAMETER[\"scale_factor\",0.9996],PARAMETER[\"false_easting\",500000],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH],AUTHORITY[\"EPSG\",\"32651\"]]";

6. 将新的投影信息设置到数据集中：

poDataset->SetProjection(pszNewProjection);

7. 释放内存：

GDALClose((GDALDatasetH)poDataset);

以上就是在 C++ 中使用 GDAL 定义投影的基本步骤。