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...