轉：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 
版權宣告：本文為博主原創文章，轉載請附上博文連結！