标签:GetTransform 影像 py geoTransform 坐标 左上角 col gdal row
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
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。