遥感影像Python GDAL 像素与坐标对应

Posted 2021-01-23 再出发

转：https://blog.csdn.net/theonegis/article/details/50805520

https://blog.csdn.net/wsp_1138886114/article/details/82778409

https://blog.csdn.net/Amrser/article/details/77595603?locationNum=3&fps=1

 

Xgeo = GT(0) + Xpixel*GT(1) + Yline*GT(2)

Ygeo = GT(3) + Xpixel*GT(4) + Yline*GT(5)

 

from osgeo import gdal

gdal.AllRegister()

filePath = ‘/home/theone/Data/GreatKhingan/DEM/Slope_GreatKhingan_500m.tif‘
dataset = gdal.Open(filePath)

adfGeoTransform = dataset.GetGeoTransform()

# 左上角地理坐标
print(adfGeoTransform[0])
print(adfGeoTransform[3])

nXSize = dataset.RasterXSize #列数
nYSize = dataset.RasterYSize #行数

arrSlope = [] # 用于存储每个像素的（X，Y）坐标
for i in range(nYSize):
    row = []
    for j in range(nXSize):
        px = adfGeoTransform[0] + i * adfGeoTransform[1] + j * adfGeoTransform[2]
        py = adfGeoTransform[3] + i * adfGeoTransform[4] + j * adfGeoTransform[5]
        col = [px, py]
        row.append(col)
    arrSlope.append(row)

print(len(arrSlope))
--------------------- 
作者：TheOneGIS 
来源：CSDN 
原文：https://blog.csdn.net/theonegis/article/details/50805520 
版权声明：本文为博主原创文章，转载请附上博文链接！