使用 GDAL/OGR C++ 将 GeoTiff 的 XY 点重新投影到纬度/经度

Question

我正在使用 GDAL/OGR API (C++) 将 XY 点从给定的 Tiff 重新投影到纬度/经度坐标。

我现在可以实现它，但是，我得到的坐标与我试图得到的实际坐标不对应。

我意识到该方法不是读取输入数据集，因此，没有“起点”。

事实上，使用命令行中的gdaltranslate 工具可以提供正确的坐标，因为我将 Tiff 作为参数发送。不幸的是，我不能使用这个命令。

以下是我的代码（到目前为止）和gdaltranslate 的示例，作为输出：

C++ 方法：

void reproject_coords(string map_biomass) {
        GDALDataset *dataset;
        GDALAllRegister();
        string ds = map_biomass;
        dataset = (GDALDataset *) GDALOpen(ds.c_str(), GA_ReadOnly);

        OGRSpatialReference source(dataset->GetMetadata()), target;
        OGRCoordinateTransformation *poCT;

        target.importFromEPSG(4326);
        //source.AutoIdentifyEPSG();
        //source.SetWellKnownGeogCS("32643");
        double x = 1794, y = 6644;

        //source.dumpReadable();
        poCT = OGRCreateCoordinateTransformation(&source, &target);

        if( poCT == NULL || !poCT->Transform( 1, &x, &y ) )
            printf( "Transformation failed.\n" );
        else
            printf( "(%f,%f)\n",
                    x, y );
}

输入：x = 1794 y = 6644
输出：(70.527326,0.059926)

命令行：

gdaltransform Karnataka_biomass.tif -t_srs EPSG:4326 -s_srs EPSG:32643

输入：x = 1794 y = 6644 输出：75.7461474090405 12.4826172522718

问题是是否有任何方法可以“读取”C++ 方法中的 Tiff 以获得真正的纬度/经度坐标。我需要吗？

Answer 1

您以错误的方式创建源 空间参考。您应该使用方法 GetProjectionRef() 而不是 GetMetaData()：

OGRSpatialReference source(dataset->GetProjectionRef()), target;

这应该在您的 GeoTIFF 中使用相同的 SRS 初始化源空间参考。

如果您只需要转换点，则无需使用 GeoTIFF。您只需初始化两个 OGRSpatialReference 对象并转换一个 OGRPoint。

更新 #1

这是一个例子：

GDALAllRegister();

OGRSpatialReference source, target;

source.importFromEPSG(32643);
target.importFromEPSG(4326);

OGRPoint p;
p.setX(1794);
p.setY(6644);
p.assignSpatialReference(&source);

p.transformTo(&target);

// transformed coordinates
cout << p.getX() << " | " << p.getY();

然而，请注意你的结果是正确的，你会得到与上面的代码相同的结果。当您查看转换后的 GeoTIFF 时看到不同值的原因是因为 gdaltransform 可能需要对其进行旋转/重新缩放，并且由于图像始终是矩形的，因此它可能会选择不同的左上角来完成此操作。

更新 #2

OP 澄清说他不想投射单个点，而是投射整个图像。此任务可以使用 GDALWarp API（教程here）来完成。

作为替代方案，仍然可以使用 gdaltransform 实用程序，通过 system() 调用来调用它。