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

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

1,，行列坐标转换地理坐标

Xp = dGeoTrans [0] +Xpixel*dGeoTrans [1]+Yline*dGeoTrans [2];

Yp = dGeoTrans [3] + Xpixel*dGeoTrans [4] + YlineL*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;