py#gdal GetTransform()详解

1）GetTransform()返回6个参数，分别为：

(左上角x坐标[0]， 水平分辨率[1]，旋转参数[2]， 左上角y坐标[3]，旋转参数[4]，竖直分辨率[5])

遥感影像如果没有发生旋转，即上北下南，则旋转参数[2]与旋转参数[4]均为0

返回的参数对应的是影像当前坐标系下的坐标值，比如WGS-84影像返回左上角x坐标[0]=经度(分为百分制)，左上角y坐标[3]=纬度；水平分辨率[1]与竖直分辨率[5]一般相等，为像元的尺寸，单位与坐标系相同。

2）影像行列号和坐标值转换关系

如果某像元行列号为row，col，则坐标(x,y)为：

x = geoTransform[0] + col * geoTransform[1] + row * geoTransform[2]
y = geoTransform[3] + col * geoTransform[4] + row * geoTransform[5]

3）已知经纬度坐标，求取影像行列号

经纬度=(lon,lat)  #经度lon对应x，纬度lat对应y，旋转参数为0

通过坐标计算行列号，即求取该坐标与左上角初始坐标的像元相对位置

行列号(row,col)=

col=int((x-geoTransform[0])/geoTransform[1])

row=int((y-geoTransform[3])/geoTransform[5])

使用int取整，因为0.9对应的是第一个像元，即编号0像元。

注意：通过坐标x计算得到的是影像的列号（col num），通过坐标y计算得到的是影像的行号（row num）

标签：GetTransform,影像,py,geoTransform,坐标,左上角,col,gdal,row
来源： https://www.cnblogs.com/misakayier/p/16101478.html