from osgeo import gdal
from gdalconst import
*import os, sys, time
start_time = time.time(
)os.chdir(r'e:\')
#獲取所有註冊類
gdal.allregister(
)#開啟柵格
ds = gdal.open(
'aster.img'
, ga_readonly)
if ds is
none
:print
('could not open image'
) sys.exit(1)
#獲取行列數和波段數
rows = ds.rasterysize
cols = ds.rasterxsize
bands = ds.rastercount
#獲取原點經緯度和像元長寬
transform = ds.getgeotransform(
)xorigin = transform[0]
yorigin = transform[3]
pixelwidth = transform[1]
pixelheight = transform[5]
#讀取所有資料
data =
for j in
range
(bands)
: band = ds.getrasterband(j+1)
data1 = band.readasarray(0,
0,cols,rows)
#print(data)
#讀取每行資料
rowsdata =
for j in
range
(bands)
: band = ds.getrasterband(j+1)
for i in
range
(rows)
: data = band.readasarray(
0, i, cols, i)
#讀取每列資料
colsdata =
for j in
range
(bands)
: band = ds.getrasterband(j+1)
for i in
range
(rows)
: data = band.readasarray(i,
0,i,rows)
PythonGIS,GDAL 讀取柵格資料
新手上路，把我在學習pythongis中的過程記錄，包括出現的問題以及解決辦法記錄下來，從而幫助更多的新手。我比較喜歡在geany上操作，以下 也是在geany上操作的。from osgeo import gdal 匯入gdal包，這步驟也可驗證你gdal包是否裝好，環境設定是否正確。import ...
python gdal開發環境安裝
安裝環境是windows10，python是最新版本3.6，gdal版本是2.x 選擇 add python to path一路next 直到完成。按快捷鍵 win r，輸入cmd回車，進入cmd控制台，輸入 py 3.6看是否安裝成功 有下面提示，則成功。重新進入cmd控制台，輸入 py 3.6 ...
python gdal讀寫shp檔案
環境配置見 安裝qgis，直接使用qgis帶的python直譯器，就不用了安裝各種環境和配置了，很方便。見 注釋寫的非常詳細。coding utf 8 try from osgeo import gdal from osgeo import ogr from osgeo import osr exc...