GDAL編譯與功能介紹

GDAL編譯與功能介紹

王健

2011-11-16GDAL(GeospatialDataAbstractionLibrary)是一個開源柵格空間數(shù)據(jù)轉(zhuǎn)換庫。它利用抽象數(shù)據(jù)模型來表達(dá)所支持的各種文件格式。并有一系列命令行工具來進(jìn)行數(shù)據(jù)轉(zhuǎn)換和處理。使用GDAL開發(fā)的軟件：ESRIArcGIS9.2+GoogleEarthERDASERViewerOpenEV

MapServerFWTools

IDRISI…….GDAL的編譯GDAL為在C/C++語言的基礎(chǔ)上編寫的類庫，相應(yīng)的，編譯的版本為C/C++版本。使用SWIG能使GDAL編譯為以下高級語言版本：C#JAVAPythonRubyVb6GDAL的基本應(yīng)用在使用GDAL的函數(shù)之前要先進(jìn)行驅(qū)動注冊GDALAllRegister(

)

GDAL常用函數(shù)：打開一個影像的數(shù)據(jù)集GDALOpen(const

char*

pszFilename,GDALAccess

eAccess)eAccess=GA_Update||GA_ReadOnly

GDALDataset*poDataset=(GDALDataset*)GDALOpen(Filename,GA_ReadOnly);doubleadfGeoTransform[6];poDataset->GetGeoTransform(adfGeoTransform)；adfGeoTransform[0]/*左上角x值*/adfGeoTransform[1]/*X軸分辨率*/adfGeoTransform[2]/*偏轉(zhuǎn)，指北則為0*/adfGeoTransform[3]/*左上角y值*/adfGeoTransform[4]/*偏轉(zhuǎn)，指北則為0*/adfGeoTransform[5]/*Y軸分辨率*/

獲取柵格影像的數(shù)據(jù)類型GDALDataTypetype=GDALGetRasterDataType(GDALGetRasterBand(hdataset,1));獲取影像的行列值。Intx=poDataset->GetRasterXSize()；Inty=poDataset->GetRasterYSize();獲取影像的投影參數(shù)。Constchar*pszproj=poDataset->GetProjectionRef();創(chuàng)建新的文件GDALDriver*hDriver=GetGDALDriverManager()->GetDriverByName("GTIFF");GDALDataset*poDataset2=hDriver->Create(savename,x,y,1,type,NULL);GDALDataset*GDALDriver::Create

(constchar*

pszFilename,int

nXSize,int

nYSize,int

nBands,GDALDataType

eType,char**

papszOptions

)pszFilename為儲存影像路徑nXSize為影像的行值nYSize為影像的列值nBands為波段數(shù)eType為影像數(shù)據(jù)類型papszOptions為一些特殊的控制參數(shù)，通常為NULL讀取影像數(shù)據(jù)int*Dstdn=newint[x*y];memset(Dstdn,0,sizeof(int)*x*y);poDataset->RasterIO(GF_Read,0,0,x,y,Dstdn,x,y,type,1,0,0,0,0);GDALDataset::RasterIO

(GDALRWFlag

eRWFlag, //打開影像方式，GF_Read||GF_Write

int

nXOff, //X軸偏移值int

nYOff, //Y軸偏移值int

nXSize //X軸讀取大小

int

nYSize, //Y軸讀取大小void*

pData, //影像讀取至緩存int

nBufXSize, //儲存影像X軸縮放大小int

nBufYSize, //儲存影像Y軸縮放大小GDALDataType

eBufType,//影像數(shù)據(jù)格式int

nBandCount, //波段數(shù)int*

panBandMap, //0int

nPixelSpace, //0int

nLineSpace, //0int

nBandSpace

//0

)#include<iostream>#include<stdlib.h>#include"gdal/gdal_priv.h"intmain(){GDALDataset*poDataset;doubleadfGeoTransform[6];float*win;char*slopeBuf; floatdx,dy,key,slopePct;

inti，j,n,containsNull;floatscale=0.5;//水平與垂直方向的變化因子 constchar*pszFilename="D:\\honda.jpg";constchar*pszSlopeFilename="D:\\slope-honda-2.bmp";GDALAllRegister();

poDataset=(GDALDataset*)GDALOpen(pszFilename,GA_ReadOnly);//打開文件

if(poDataset==NULL){

std::cout<<"Couldn'topendataset"<<pszFilename<<"/n";}

GDALRasterBand*poBand;

poBand=poDataset->GetRasterBand(1);

poDataset->GetGeoTransform(adfGeoTransform);constdoublecellsizeY=adfGeoTransform[5];//Y軸分辨率

constdoublecellsizeX=adfGeoTransform[1];//X軸分辨率

constfloatnullValue=(float)poBand->GetNoDataValue();//獲取缺省值

constint

nXSize=poBand->GetXSize();//列數(shù)

constint

nYSize=poBand->GetYSize();//行數(shù)

slopeBuf=(char*)CPLMalloc(sizeof(char)*nXSize);win=(float*)CPLMalloc(sizeof(float)*9);

GDALDriver*poDriver;

poDriver=GetGDALDriverManager()->GetDriverByName("BMP");//文件輸出類型為TIF

GDALDataset*poSlopeDS;GDALRasterBand*poSlopeBand;poSlopeDS=poDriver->Create(pszSlopeFilename,nXSize,nYSize,1,GDT_Byte,NULL);poSlopeBand=poSlopeDS->GetRasterBand(1);for(i=0;i<nYSize;i++){for(j=0;j<nXSize;j++){

containsNull=0;

if(i==0||j==0||i==nYSize-1||j==nXSize-1)//邊界操作{

slopeBuf[j]=nullValue;continue;}

poBand->RasterIO(GF_Read,j-1,i-1,3,3,win,3,3,GDT_Float32,0,0);//讀取3X3的矩陣

for(n=0;n<=8;n++){

if(win[n]==nullValue)//檢查矩陣中是否有NULL值{

containsNull=1;break;}}if(containsNull==1){

slopeBuf[j]=nullValue;continue;}else{

dx=((win[0]+win[3]+win[3]+win[6])- // SOBEL算子為3X3的矩陣，

(win[2]+win[5]+win[5]+win[8])); // 如下所示為下標(biāo)

// 0 1 2

dy=((win[6]+win[7]+win[7]+win[8])- // 3 4 5(win[0]+win[1]+win[1]+win[2])); // 6 7 8 key=((dx/(8*cellsizeX*scale))*(dx/(8*cellsizeX*scale)))+ ((dy/(8*cellsizeY*scale))*(dy/(8*cellsizeY*scale)));

slopePct=100*sqrt(key);

if(slopePct>830)