GDAL 地理坐标与图像行列坐标相互转换公式

通过GetGeoTransform(dGeoTrans)获取仿射信息

 

1、行列坐标转换地理坐标

Xp = dGeoTrans [0] +Xpixel*dGeoTrans [1]+Yline*dGeoTrans [2];
Yp = dGeoTrans [3] + Xpixel*dGeoTrans [4] + Yline*dGeoTrans[5];

 

2、地理坐标转换行列坐标
double dTemp = dGeoTrans[1] * dGeoTrans[5] - dGeoTrans[2] *dGeoTrans[4];
Xpixel= (dGeoTrans[5] * (Xp - dGeoTrans[0]) -dGeoTrans[2] * (Yp - dGeoTrans[3])) / dTemp + 0.5;
Yline = (dGeoTrans[1] * (Yp - dGeoTrans[3]) -dGeoTrans[4] * (Xp - dGeoTrans[0])) / dTemp + 0.5;

参考原文：https://blog.csdn.net/zyf19930610/article/details/43487473