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

http://www.gdal.org/ogr/ogr_apitut.html

又：

GDAL介绍、安装及使用 

1.GDAL介绍

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

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

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

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

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

目前ogr能够支持的数据格式包括：
Arc/Info Binary Coverage、DWG、ESRI PersonalGeoDatabase、ArcSDE、ESRI
Shapefile、GML、GRASS、Mapinfo File、Microstation DGN、ODBC、OracleSpatial和
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

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

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

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

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

行nmake /f makefile.vcinstall，然后是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:\>test1.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++)
       {    
           CStringfilename=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");
           else if (suffix == "hdr")
               driver =GetGDALDriverManager()->GetDriverByName("EHdr");
            
    
           int w=tmpDataset->GetRasterXSize();
           inth=tmpDataset->GetRasterYSize();            
           GDALDataset*maxNDV=driver->Create(maxNDVfile,w,h,1,GDT_Float32,NULL);  

     
           int xOff=0;
           int yOff=0;
           int width=w;
           int height=1;//一行一行地读取
           vector<float*> bufs;
           for(int i=0;i<datasets.size();i++)
           {
               float *buf=new float[width*height];
               bufs.push_back(buf);
               buf=NULL;
           }        
           float *newbuf=new float[width*height];
           GDALRasterBand *newpoband=maxNDV->GetRasterBand(1);        
           for(int j=0;j<h;j++)
           {                
               for(int i=0;i<datasets.size();i++)
               {
                   GDALDataset *tmpDt=datasets[i];
                   GDALRasterBand *poband=tmpDt->GetRasterBand(1);                
                   float *buf=bufs[i];
                   xOff=0;
                   yOff=j;
                   poband->RasterIO

(GF_Read,xOff,yOff,width,height,buf,width,height,GDT_Float32,0,0);
                   buf=NULL;
                   tmpDt=NULL;//在后面再释放,因为很多个指针指向同一一段内存，不能随

便释放内存哦
               }            
               //////////
               //在这可以开始运算咯
               for(int k=0;k<w;k++)//对于每个像元
               {
                   //下面这段代码以后要支持均值，极差，方差
                   for(int kk=0;kk<bufs.size();kk++)
                   {
                       float *buf=bufs[kk];
                       if(kk==0)//初始化
                       {
                           newbuf[k]=buf[k];
                       }
                       if(buf[k]>newbuf[k])//最大的,稍微改下就是最小值，
                       {
                           newbuf[k]=buf[k];
                       }                    
                   }
               }
               //////////写入数据
               newpoband->RasterIO

(GF_Write,xOff,yOff,width,height,newbuf,width,height,GDT_Float32,0,0);
           }
           delete newbuf;
           newbuf=NULL;
           for( i=0;i<bufs.size();i++)
           {
               if(bufs[i]!=NULL)
                   deletebufs[i];            
           }
           bufs.resize(0);
           delete maxNDV;
           maxNDV=NULL;
           delete driver;
           driver=NULL;
       }
       tmpDataset=NULL;
       //释放掉内存
       for( i=0;i<datasets.size();i++)
       {
           if(datasets[i]!=NULL)
           {
               delete datasets[i];
               datasets[i]=NULL;
           }
       }
       datasets.resize(0);
}