Question

我正在使用GDAL库。目前，我可以接受左上角和右上角，并从原始图像中切出图像。我现在要做的是，接受两个WKT点并转换为X，Y坐标来做同样的事情。我只是想知道如果我知道GeoTransform及其使用的坐标系统（WGS84）是否可以这样做？

Answer 1

我之前也遇到了这个问题，这里有一个很好的方法来进行坐标转换。

GDALDataset :: GetProjectionRef（）返回的坐标系描述仿射所暗示的地理参考坐标 GDALDataset :: GetGeoTransform（）返回的地理配准变换。

我们可以将其与OGRCoordinateTransformation一起用于为我们进行转换。

基本上代码看起来像这样：

// Load up some dataset.
dataset = (GDALDataset *) GDALOpen( mapfile, GA_ReadOnly );

// Define Geographic coordinate system - set it to WGS84.
OGRSpatialReference *poSRS_Geog = new OGRSpatialReference();
poSRS_Geog->importFromEPSG( 4326 ); // WGS84

// Define Projected coordinate system - set to the GeoTransform.
const char *sProj = dataset->GetProjectionRef();
OGRSpatialReference *poSRS_Proj = new OGRSpatialReference( sProj );

// Set up the coordinate transform (geographic-to-projected).
OGRCoordinateTransformation *poCT_Geog2Proj;
poCT_Geog2Proj = OGRCreateCoordinateTransformation( poSRS_Geog, poSRS_Proj );

// Now everything is set up and we set transforming coordinates!
// Pass Lon/Lat coordinates to the Transform function:
double x = lon;
double y = lat;
poCT_Geog2Proj->Transform( 1, &x, &y );

// Now x and y variables will contain the X/Y pixel coordinates.

这就是你可以在经度/纬度和像素坐标之间进行转换的方法。请注意，您可以将数组与Transform()一起使用，并将多个坐标转换为一起。第一个参数是要转换的坐标对的数量，第二个和第三个参数是指向x和y的指针。我只是在这里改造了一对。

注意设置逆变换同样容易：

// Set up the coordinate transform (projected-to-geographic).
OGRCoordinateTransformation *poCT_Proj2Geog;
poCT_Proj2Geog = OGRCreateCoordinateTransformation( poSRS_Proj, poSRS_Geog );

Answer 2

我使用图像的仿射变换来计算一些样本纬度/经度。我遇到的唯一问题是如果图像是朝向北方，则在计算Lat / Lon时需要将geoTransform [2]和geTransform [4]清零。
x = (int)Math.Abs(Math.Round((Latitude - geotransform[0]) / geotransform[1])); y = (int)Math.Abs(Math.Round((Longitude - geotransform[3]) / geotransform[5]));

如果您想暴力破解它，您可以执行以下操作（我这样做但它有效，但这只是伪代码）：

//Get the Pixel for the length and width, this portion is for the full image pixelXSize = AbsoluteValue((latitudeAt(Zero)-(latitudeAt(Length)-1))/imageLength); pixelYSize = AbsoluteValue((longitudeAt(Zero)-(LongitudeAt(Width)-1))/imageWidth);

//Calculate the x,y conversion for the points you want to calculate x = AbsoluteValue((latitudeToConvert-latitudeAt(Zero))/pixelXSize); y = AbsoluteValue((longitudeToConvert-longitudteAt(Zero))/pixelYSize);

这个答案可能会偏离一两个像素。如果您正在使用地理转换，则North Facing变量可能会使返回的答案陷入混乱。到目前为止，我只针对朝北的图像进行了测试。

Answer 3

我使用这种方法：

void transformCoordinatesEPSG(OGRGeometry &geometry,int from, int to) {
    OGRSpatialReference srcSpatialReference;
    OGRErr error = srcSpatialReference.importFromEPSG(from);

    #ifdef __OGRTRANSFORMDEBUG
        qDebug() << "Import EPSG  " << from << "return " << error;
    #endif

    OGRSpatialReference dstSpatialReference;
    error = error | dstSpatialReference.importFromEPSG(to);

    #ifdef __OGRTRANSFORMDEBUG
    qDebug() << "Import EPSG  " << to << "return " << error;
    #endif

    OGRCoordinateTransformation* coordTrans = OGRCreateCoordinateTransformation(&srcSpatialReference, &dstSpatialReference);
    geometry.transform(coordTrans);
}

对于lat / long来说必须是4326。

如果你有GeoTiff，是否可以使用GeoTransform将Lat / Lon点转换为X，Y？

3 个答案: