gdal中shapefile座標度轉換到柵格公尺

2021-08-16 15:25:52 字數 999 閱讀 1734

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,要是你會改目錄,找這個東西...