OGR教程 http://www.gdal.org…

GDAL介绍、安装及使用 

1.GDAL介绍

GDAL是一个转换各种栅格数据格式的库。包括读取、写入、转换、处理各种栅格数据格式（有些

特定的格式对一些操作如写入等不支持）。它使用了一个单一的抽象数据模型就支持了大多数的

栅格数据（GIS对栅格，矢量，3D数据模型的抽象能力实在令人叹服）。当然除了栅格操作，这

个库还同时包括了操作矢量数据的另一个有名的库ogr（转换矢量GIS数据），这样这个库就同时

具备了操作栅格和矢量数据的能力，

目前ogr能够支持的数据格式包括：
Arc/Info Binary Coverage、DWG、ESRI Personal
GeoDatabase、ArcSDE、ESRI
Shapefile、GML、GRASS、Mapinfo File、Microstation DGN、ODBC、Oracle
Spatial和
PostgreSQL等。应该说，这就基本包括了我们平常用到的所有矢量型GIS文件格式了。
Gdal支持的栅格数据格式参阅http://www.gdal.org/formats_list.html

2.GDAL安装

（1）下载gdal的安装文件，http://download.osgeo.org/gdal/gdal142.zip，解压到某目录下

，如C:\gdalsrc下。
这里我们假定VC6的安装在默认目录C:\Program Files\Microsoft Visual Studio8下。
（2）启动cmd，即打开控制台窗口。进入VC6的安装目录下，如cd C:\Program Files\Microsoft

Visual
Studio8\VC\bin\，在此目录下有个文件VCVARS32.BAT，执行一下这个文件，然后重新回

到C:\gdalsrc下。运行命令nmake /f
makefile.vc。编译完成后，用记事本打开文件

C:\gdalsrc\nmake.opt，根据自己的情况修改GDAL_HOME = 这一行，这个指的是最终GDAL的安装

目录，比如说我们安装在C:\GDAL，那么这一行就是GDAL_HOME = “C:\GDAL”，在C:\gdalsrc下执

行nmake /f makefile.vc
install，然后是nmake /f makefile.vc devinstall，然后我们需要

的东西就安装到了C:\GDAL下。

3.GDAL使用

初次使用visual studio
对工具不熟悉，有些步骤可以配置的

（1）在VS2005中新建win32控制台程序
testGDALconsole，（向导中附加选项不能选为空项目）将C:\GDAL\bin\gdal14.dll拷贝到testGDALconsole目录下的debug目录中。（否则运行时会提示找不到gdal14.dll）

（2）在工程的Library files中和Include files中分别添加GDAL的LIB文件目录和头文件目录，

这步我不会，就直接将gdal_priv.h拷贝到testGDALconsole.cpp所在目录下

（3）添加一个文件gdal_i.lib。使用
#pragma comment (lib,”..\\lib\\gdal_i.lib”)（这步没做）
（4）程序文件

样本程序如下：

#include “stdafx.h”

#include “gdal_priv.h”

#include
<iostream>

using namespace std;

int main(int argc, char*
argv[])

{

 

      
if(argc!=2)

      
{

             
return 0;

      
}

 

      
GDALDataset *poDataset;

 

   
GDALAllRegister();

 

    poDataset =
(GDALDataset *) GDALOpen( argv[1], GA_ReadOnly );

    if(
poDataset != NULL )

    {

       
cout<<“RasterXSize:”<<poDataset->GetRasterXSize()<<endl;

             
cout<<“RasterYSize:”<<poDataset->GetRasterYSize()<<endl;

             
cout<<“RasterCount:”<<poDataset->GetRasterCount()<<endl;

    }

 

      
return 0;

}

 

假定生成的程序名为test.exe，假定有一个文件名为1.jpg的文件在C:\下，在命令行下test.exe

c:\1.jpg，程序将打印出此图片的横纵坐标和raster count，如：

 

C:\>test
1.JPG

RasterXSize:800

RasterYSize:600

RasterCount:3

练习总结：

1）使用 #pragma comment
(lib,”..\\lib\\gdal_i.lib”)来引入外部的库文件

或者在Visual studio里面设置

project－>
properties-> Linker->
Input里面的Additional Dependencies里面加lib 
          
－> properties->
General-> Additional Libary Directories加Lib目录

还没有解决，未知原因，以后熟悉了在看

下面附上一个一个GDAL的读写数据例子
引自 http://blog.csdn.net/tangnf/archive/2008/10/26/3152538.aspx

//

#include “stdafx.h”
#include “fangshibo.h”

#include
<vector>

#include<math.h>
///gdal头文件
#include “..\\include\\gdal.h”
#include “..\\include\\gdal_priv.h”
#include “..\\include\\ogr_srs_api.h”
#include “..\\include\\cpl_string.h”
#include “..\\include\\cpl_conv.h”
#pragma comment (lib,”..\\lib\\gdal_i.lib”)
/////////////////////////////////////////////////////////////////////////////

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////

// The one and only application object

CWinApp theApp;
using namespace std;

///////////////////////
void MaxNDVI(vector<CString>
files,CString maxNDVfile)
{    

//   
GDALAllRegister();
   
vector<GDALDataset *>
datasets;    

       
//打开文件
       
for(int i=0;i<files.size();i++)
       
{    

           
CString
filename=files[i];            

           
GDALDataset *tmpDataset = (GDALDataset *) GDALOpen( files[i],

GA_ReadOnly);
           
if(tmpDataset!=NULL)
           
{                

               
datasets.push_back(tmpDataset);                

           
}
           
else
           
{
               
fprintf( stderr, “GDALOpen failed – %d\n%s\n”,
CPLGetLastErrorNo(),

CPLGetLastErrorMsg() );
           
}    

           
tmpDataset=NULL;//以后再释放 
       
}
       
////读取数据咯，按行来读取
       
//GDALDir
       
if(datasets.size()==0)return;
       
GDALDataset *tmpDataset=datasets[0];
       
if(tmpDataset!=NULL)
       
{    

           
GDALDriver * driver=NULL;
           
int index = maxNDVfile.ReverseFind(‘.’);
           
if (index < 0)
               
return ;
           
if (index == maxNDVfile.GetLength()-1)
               
return ;
           
CString suffix =
maxNDVfile.Right(maxNDVfile.GetLength()-1-index);      

     

           
suffix.MakeLower();
           
if (suffix == “bmp”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“BMP”);

           
else if (suffix == “jpg”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“JPEG”);

           
else if (suffix == “tif”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“GTiff”);

           
else if (suffix == “img”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“HFA”);

           
else if (suffix == “bt”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“BT”);
           
else if (suffix == “ecw”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“ECW”);

           
else if (suffix == “fits”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“FITS”);

           
else if (suffix == “gif”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“GIF”);

           
else if (suffix == “hdf”)
               
driver =
GetGDALDriverManager()->GetDriverByName(“HDF4”);