【遙感影像】Python GDAL 畫素與座標對應
阿新 • • 發佈:2018-11-30
轉: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)座標 fori 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 版權宣告:本文為博主原創文章,轉載請附上博文連結!