gdal 建立GeoTiff柵格資料

2022-09-17 12:09:10 字數 1275 閱讀 1924

1

//定義轉換引數

2private

readonly

double d_transform = ;3//

定義座標系統

4private

readonly

string d_project = "

geogcs[\"wgs 84\",datum[\"wgs_1984\",spheroid[\"wgs 84\",6378137,298.257223563,authority[\"epsg\",\"7030\"]],authority[\"epsg\",\"6326\"]],primem[\"greenwich\",0],unit[\"degree\",0.0174532925199433],authority[\"epsg\",\"4326\"]]";

56//gdal註冊

7gdal.allregister();

8 driver = gdal.getdriverbyname("

gtiff

");

9 dataset ds = driver.create(filepath, xsize, ysize, 1, datatype.gdt_float32, new

string

);10 ds.setgeotransform(d_transform); //

影像轉換引數

11 ds.setprojection(d_project); //

投影12

//壓縮方式,設定後重新讀取發現並沒有這個metadata引數,具體情況未知

13 ds.setmetadataitem("

compression

", "

lzw", "

image_structure");

14//

初始化陣列

15 ds.getrasterband(1).writeraster(0, 0, xsize, ysize, buffer, xsize, ysize, 0, 0

);16 ds.getrasterband(1

).flushcache();

17ds.flushcache();

18ds.dispose();

1920

//其中buffer為float型別陣列,長度為[xsize*ysize];

重新獲取可以通過

string metadata = ds.getmetadata("

image_structure

");

GDAL向量轉柵格

gdal版本1.8.0.首先嘗試使用gdal rasterize小工具,但是不支援新生成輸出柵格檔案。官方說1.8版本後已經支援。但仍然無效。柵格的元資料賦值和仿射變換以後需要進一步研究。實現 include stdafx.h include include gdal priv.h include ...

GDAL 處理柵格資料

gdal 是柵格資料處理模組 ogr 是向量資料處理模組 gdal資料驅動,與ogr資料驅動類似,需要先建立某一型別的資料驅動,再建立響應的柵格資料集。一次性註冊所有的資料驅動,但是只能讀不能寫 gdal.allregister 單獨註冊某一型別的資料驅動,這樣的話可以讀也可以寫,可以新建資料集 d...

GDAL建立投影

gdal中,使用的是 wkt 串來表示投影,具體的表示內容可以參考鏈結,下面用例子簡單的介紹一下,後面表示注釋 projcs wgs 84 utm zone 52n 投影名稱 geogcs wgs 84 地理座標系統名 datum wgs 1984 水平基準面 spheroid wgs 84 637...