OGR读取PostgreSQL/PostGIS格式数据

在OGR官网中可以看到OGR库支持PostgreSQL/PostGIS数据库的读写，但是编译的时候需要加入PostgreSQL client library (libpq)，下面就介绍一下如何编译GDAL/OGR使其可以支持PostgreSQL/PostGIS数据库的访问。

首先，安装PostgreSQL/PostGIS。安装步骤网上有好多，也并不复杂，不需要特殊设置，具体步骤可参考http://wenku.baidu.com/view/ac2b40ab284ac850ad0242f3.html。

安装完成后到安装目录C:Program Files (x86)PostgreSQL9.2lib下，找到libpq.lib文件，这个文件就是要在编译GDAL/OGR时加入的那个文件。为了方便起见，可以将其拷贝的浅层目录下，在这里，我将其拷贝到C盘根目录下。同时C:Program Files (x86)PostgreSQL9.2目录下的include文件夹拷贝到C盘根目录下，这个文件夹也要添加到GDAL编译文件中。

现在，打开GDAL源码的解决方案，到Make Files文件夹中中的nmake.opt文件，定位到364行，将下面三行

# PostGIS Libraries
#PG_INC_DIR = n:pkglibpq_win32include
#PG_LIB = n:pkglibpq_win32liblibpqdll.lib wsock32.lib

改为

# PostGIS Libraries
PG_INC_DIR = C:include
PG_LIB = C:libpq.lib

由于我已经将include文件夹和libpq.lib拷贝到C盘根目录下，所以PG_INC_DIR = C:include，PG_LIB = C:libpq.lib。如果你将这两个文件放到了其他目录下，那就将此处改为你存放的那个目录。

修改完成后，编译GDAL。编译成功后，将C:Program Files (x86)PostgreSQL9.2binpostgisgui目录下的libpq.dll和libintl.dll文件拷贝到GDAL编译生成的目录，默认是C:warmerdabldbin 下面。现在通过命令行到C:warmerdabldbin 目录下运行 ogrinfo.exe --formats，出现下图结果：

可以看到，OGR已经支持PostgreSQL的读写。

下面介绍在程序里面如何读取数据库中的表。新建C++ 控制台程序，按GDAL程序的搭建步骤，在项目属性里设置好包含目录和库目录，并将编译生成的gdal110.dll和libpq.dll、libintl.dll三个文件拷贝到项目的Debug文件夹下。PostgreSQL数据源的连接字符串如下：

PG:dbname=postgis host=localhost port=5432 user=postgres password=renyibin

下面的代码是连接数据库，并读取数据库中名为layer1的表，显示layer1中的要素个数。

OGRRegisterAll();
const char *filePath = "PG:dbname=postgis host=localhost port=5432 user=postgres password=renyibin";
const char *driveName = "PostgreSQL";
const char *pLayerName = "layer1";
OGRSFDriver *pDriver = NULL;
OGRLayer *pLayer = NULL;
OGRDataSource *pODS = NULL; 
pDriver = OGRSFDriverRegistrar::GetRegistrar()->GetDriverByName(driveName);
if (pDriver == NULL)
{
    cout<<"Can not get PostgreSQL Driver"<<endl;
    exit(1);
}
pODS = pDriver->Open(filePath,0);
if (pODS == NULL)
{
     cout<<"Can not open "<<filePath<<endl;
     exit(1);
}
pLayer = pODS->GetLayerByName(pLayerName);
if (pLayer == NULL)
{
    cout<<"Can not get layer "<<pLayerName<<endl;
    exit(1);
}
cout<<"feature count of layer1 is "<<pLayer->GetFeatureCount(1)<<endl;
OGRDataSource::DestroyDataSource(pODS);
system("pause");

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!