代碼具有詳細說明注釋

所需庫 #gdalconst中的常量都加了前綴，力圖與其他的module沖突最小 from osgeo import gdal,gdalconst from osgeo import gdal_array as ga # 用于引入一個模塊的同時為該模塊取一個別名 from osgeo.gdalconst import GA_ReadOnly# 開始對柵格的操作# GDAL所有操作都需要先注冊格式# 一次性注冊所有的數據驅動，但是只能讀不能寫：gdal.AllRegister()gdal.AllRegister()#打開數據集，并傳遞數據集的名稱和所需的訪問權限（GA_ReadOnly或GA_Update）#img柵格絕對路徑dr = gdal.Open(img, GA_ReadOnly)if dr is None:print('打開柵格文件失敗')sys.exit(1)print("打開柵格后的數據")# 讀取圖像y方向上的像素個數rows = dr.RasterYSize# 讀取圖像x方向上的像素個數cols = dr.RasterXSize# 波段數bands = dr.RasterCount#存儲著柵格數據集的地理坐標信息transform = dr.GetGeoTransform()#影像左上角橫坐標xOrigin = transform[0]#影像左上角縱坐標yOrigin = transform[3]#遙感圖像的水平空間分辨率或者東西方向上的像素分辨率pixelWidth = transform[1]#遙感圖像的垂直空間分辨率或者南北方向上的像素分辨率pixelHeight = transform[5]#通常geoTransform[5] 與 geoTransform[1]相等

創作挑戰賽新人創作獎勵來咯，堅持創作打卡瓜分現金大獎

總結

以上是生活随笔為你收集整理的利用python读取栅格数据的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。