根据仿射变换参数实现从地理坐标转换为图像坐标的一个思考

作者：朱金灿

来源：http://blog.csdn.net/clever101

          六个参数分别是：
    geos[0]  top left x 左上角geox坐标
    geos[1]  w-e pixel resolution 东西方向像素分辨率
    geos[2]  rotation, 0 if image is "north up" 旋转角度，正北向上时为0
    geos[3]  top left y 左上角geoy坐标
    geos[4]  rotation, 0 if image is "north up" 旋转角度，正北向上时为0
    geos[5]  n-s pixel resolution 南北向像素分辨率
    x/y为图像的x/y坐标，geox/geoy为对应的投影坐标

          一般而言，从仿射变换参数实现从地理坐标转换为图像坐标（行列号坐标）遵循以下公式：

Col = GT(0) + Xgeo*GT(1)+ Ygeo*GT(2)

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

 

        上面公式中的GT(0)为仿射变换参数的第一个参数，如此类推，Xgeo为地理坐标中x方向的值，Ygeo为地理坐标中y方向的值，Col为列号，Row为行号。

 

         实际上这样转换是有误差的，因为比如0行0列的像素实际上其中心点坐标为0.5行和0.5列，如下图：

 

         以中心点代表指定位置的像素显然更为精确，在计算行列坐标时，需要将这0.5的误差给补齐，因此更精确的转换公式应该是：

Col = GT(0) + Xgeo*GT(1) + Ygeo*GT(2) + 0.5

Row = GT(3) + Xgeo*GT(4) + Ygeo*GT(5)+0.5