Wechat: yu389741| Email: gisdqy@163.com

Shop:https://www.giserdqy.com/shop

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)\PostgreSQL\9.2\lib下,找到libpq.lib文件,这个文件就是要在编译GDAL/OGR时加入的那个文件。为了方便起见,可以将其拷贝的浅层目录下,在这里,我将其拷贝到C盘根目录下。同时C:\Program Files (x86)\PostgreSQL\9.2目录下的include文件夹拷贝到C盘根目录下,这个文件夹也要添加到GDAL编译文件中。

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

       # PostGIS Libraries
       #PG_INC_DIR = n:\pkg\libpq_win32\include
       #PG_LIB = n:\pkg\libpq_win32\lib\libpqdll.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)\PostgreSQL\9.2\bin\postgisgui目录下的libpq.dll和libintl.dll文件拷贝到GDAL编译生成的目录,默认是C:\warmerda\bld\bin 下面。现在通过命令行到C:\warmerda\bld\bin 目录下运行 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");

输出结果为:

转载自:https://blog.csdn.net/rybgis/article/details/12858105