GDAL处理地理图像坐标计算

本文讲解如何使用GDAL处理地理图像时，通过使用行列号计算和转换成tiff图像的地理坐标：

tif中坐标计算的方法如下，其中Col表示该坐标点处图像的列号，ROW表示该坐标点处图像的行号， 比如图像左上角Col为0，ROW为0，图像右下角Col为图像宽度，ROW为图像高度。  Xgeo = GT(0) + Col*GT(1) + Row*GT(2)

 Ygeo = GT(3) + Col*GT(4) + Row*GT(5)

一、文中用到的系数说明：

Col表示该坐标点处图像的列号 ROW表示该坐标点处图像的行号 x:该点X坐标 y:该点Y坐标 upleft:图像左上角坐标 upright:图像右上角坐标 downleft:图像左下角坐标 downright:图像右下角坐标

二、公式推导

则根据计算公式： GT(0)=左上点的Xgeo=upleft(x) GT(3)=坐上点的Ygeo=upleft(y) upRighg(x)=GT(0)+srcWidth*GT(1)+0*GT(2) =>GT(1)=(upRight(x)-GT(0))/srcWidth; downLeft(x)=GT(0)+0*GT(1)+srcHeight*GT(2) =>GT(2)=(downLeft(x)-GT(0))/srcHeight upRight(y)=GT(3)+srcWidth*GT(4)+0*GT(5) =>GT(4)=(upRight(y)-GT(3))/srcWidth; downLeft(y)=GT(3)+0*GT(4)+srcHeight*GT(5) =>GT(5)=(downLeft(y)-GT(3))/srcHeight;

(三）系数验证

验证：右下角(假设图像为平行四边形） downright=[downleft(x)+（upright(x)-upleft(x)）, upright(y)-(upleft(y)-downleft(y))]

（1）证明右下角X方向坐标成立：downright(x)=downleft(x)+（upright(x)-upleft(x)

downright(x)=GT(0) + Col*GT(1) + Row*GT(2) =upleft(x)+srcWidth*(upRight(x)-GT(0))/srcWidth+srcHeight*(downLeft(x)-GT(0))/srcHeight =upleft(x)+upRight(x)-GT(0)+downLeft(x)-GT(0) =upleft(x)+upRight(x)-upleft(x)+downLeft(x)-upleft(x) =upRight(x)+downLeft(x)-upLeft(x) 证明成立（系数GT(0)，GT(1)，GT(2)正确）

（2）证明右下角Y方向坐标成立：

downright(y)=upright(y)-upleft(y）+downLeft(y) downright(y)=GT(3) + Col*GT(4) + Row*GT(5) =upleft(y)+srcWidth*GT(4)+srcHeight*GT(5) =upleft(y)+srcWidht*(upRight(y)-GT(3))/srcWidth+srcHeight*(downLeft(y)-GT(3))/srcHeight =upleft(y)+upRight(y)-GT(3)+downLeft(y)-GT(3) =upleft(y)+upRight(y)-upleft(y)+downLeft(y)-upleft(y) =upRight(y)+downLeft(y)-upleft(y)

证明成立（系数GT(3)，GT(4)，GT(5)正确）

(四)结果图