python gdal讀取(ogr)shapefile的座標表示為度,柵格資料的座標為公尺。
以下**展示讀取shapefile檔案的度座標並轉換為依賴柵格資料投影的座標公尺。
from osgeo import gdal,ogr
defworld2pixel
(padftransform, x, y):
pixel = padftransform[0] + x*padftransform[1] + y*padftransform[2]
line = padftransform[3] + x*padftransform[4] + y*padftransform[5]
return (pixel, line)
dataset = gdal.open("e:/rs_data/caijian1214/res_data.tif")
driver = ogr.getdriverbyname("esri shapefile")
datasource = driver.open(u"e:\rs_data\caijian1214\lastoutput_clip.shp")
layer = datasource.getlayer(0)
minx, maxx, miny, maxy = layer.getextent()
print("原邊界(座標系度):",minx,maxx,maxy,miny)
geotrans = dataset.getgeotransform()
ulx, uly = world2pixel(geotrans, minx, maxy)
lrx, lry = world2pixel(geotrans, maxx, miny)
print("新邊界:(座標系公尺)",ulx,uly,lrx,lry)
輸出結果如下:
pycharm中配置GDAL庫
table of contents 2 安裝 3 在pycharm中新增gdal庫 然後將其複製到python安裝路徑下的scripts資料夾裡面 開啟cmd對話方塊,定位到gdal 2.4.1 cp37 cp37m win amd64.whl的存放路徑 然後輸入pip install gdal 2...
在C 中初次使用GDAL
如果要處理柵格資料,特別是格式轉換,地理空間資料抽象庫 gdal 是非常好的。我想使用gdal進行生物多樣性建模專案,所以我看了一下gdal的c 繫結。按照如下步驟你就可以在vs 2010以及vs 2012使用gdal了,所以您需要做的是 設定系統環境變數。在path系統變數中新增包含c progr...
VC 2010中GDAL的安裝
準備檔案 vs2010 必須滴 gdal的原始碼 我用的這個版本的 1 編譯gdal,獲取include和lib和dll 2 啟動cmd,進入vc10的安裝目錄下,預設為c program files microsoft visual studio 10.0 vc bin,要是你會改目錄,找這個東西...