代碼具有詳細(xì)說明注釋

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

創(chuàng)作挑戰(zhàn)賽新人創(chuàng)作獎(jiǎng)勵(lì)來咯，堅(jiān)持創(chuàng)作打卡瓜分現(xiàn)金大獎(jiǎng)

總結(jié)

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

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。