OGR结果体系

在网上看到这篇文章写得非常全面,而且是中文的,就转过来了,学习一下
其中有部分内容是重复的,不过不影响,非常好的资料

主要转自http://blog.csdn.net/liminlu0314?viewmode=contents

一、简介

1.GDAL

简单地说,GDAL是一个操作各种栅格地理数据格式的库。包括读取、写入、转换、处理各种栅格数据格式(有些特定的格式对一些操作如写入等不支持)。它使用了一个单一的抽象数据模型就支持了大多数的栅格数据(GIS对栅格,矢量,3D数据模型的抽象能力实在令人叹服)。当然除了栅格操作,这个库还同时包括了操作矢量数据的另一个有名的库ogr(ogr这个库另外介绍),这样这个库就同时具备了操作栅格和矢量数据的能力。

最重要的是这个库是跨平台的,开源的!如今这个库对各种数据格式的支持相当强大,几乎所有的GIS、RS所要处理的数据格式类型都包括在内。它的主页最后那些使用了它作为底层数据处理的软件列表中,不仅有GRASSQuantum
GIS (QGIS)
,更有Google Earth

学习这个库是十分有必要的。该库在C/C++/python/ruby/VB/java/C#(这个暂时不完全支持)下都能用,而且都一样用。

 

GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。

 

2.Ogr

OGR是GDAL项目的一个分支,功能与GDAL类似,只不过它提供对矢量数据的支持。   有很多著名的GIS类产品都使用了GDAL/OGR库,包括ESRI的ArgGIS 9.2,Google Earth和跨平台的GRASS GIS系统。  

OGR提供对矢量数据格式的读写支持,它所支持的文件格式包括:ESRI Shapefiles, S-57, SDTS, PostGIS,Oracle Spatial, Mapinfo mid/mif , Mapinfo TAB。

1)OGR体系结构

OGR包括如下几部分:

Geometry:类Geometry (包括OGRGeometry等类)封装了OpenGIS的矢量数据模型,并提供了一些几何操作,WKB(Well Knows Binary)和WKT(Well Known Text)格式之间的相互转换,以及空间参考系统(投影)。

Spatial Reference:类OGRSpatialReference封装了投影和基准面的定义。

Feature:类OGRFeature封装了一个完整feature的定义,一个完整的feature包括一个geometry和geometry的一系列属性。

Feature Definition:类OGRFeatureDefn里面封装了feature的属性,类型、名称及其默认的空间参考系统等。一个OGRFeatureDefn对象通常与一个层(layer)对应。

Layer:类OGRLayer是一个抽象基类,表示数据源类OGRDataSource里面的一层要素(feature)。

Data Source:类OGRDataSource是一个抽象基类,表示含有OGRLayer对象的一个文件或一个数据库。

Drivers:类OGRSFDriver对应于每一个所支持的矢量文件格式。类OGRSFDriver由类OGRSFDriverRegistrar来注册和管理。

 

2) OGR的Geometry模型

 

OGR的Geometry模型是建立在OpenGIS的简单要素数据模型之上的。如下图所示:

 

 

图-OGR的Geometry模型关系图

 

OpenGIS的简单要素数据模型,其关系图如下所示:

 

图-OpenGIS的简单要素数据模型

 

由上面两图的对比,可以清楚的看到,OGR的Geometry模型是严格遵循OpenGIS的简单要素数据规范的。OGR的Geometry模型不仅在继承体系上与OpenGIS的简单要素数据模型一致,在函数接口上也向其靠拢,从基本的获取Geometry对象信息的方法如Dimension( )、GeometryType ( )、SRID ( )、Envelope( )、AsText( )、Boundary( )等到判定空间未知关系的方法如Equals(anotherGeometry:Geometry)、Disjoint(anotherGeometry:Geometry)、Intersects(anotherGeometry:Geometry)、Touches(anotherGeometry:Geometry)等都是符合其标准的。

3)OGR的API使用范例

以下C++示例代码展示了使用OGR的API来读取矢量数据。

 

  1. #include “ogrsf_frmts.h”  
  2. int main()  
  3. {      
  4. //注册所有的文件格式驱动      
  5. OGRRegisterAll();      
  6. //打开point.shp文件      
  7. OGRDataSource *poDS = OGRSFDriverRegistrar::Open( “point.shp”, FALSE );      
  8. //获取点层      
  9. OGRLayer *poLayer = poDS->GetLayerByName( “point” );      
  10. OGRFeature *poFeature;      
  11. //重置该层,确保从层的开始读取数据      
  12. poLayer->ResetReading();      
  13. while ( (poFeature = poLayer->GetNextFeature()) != NULL )      
  14. {          
  15. //获取该层的属性信息          
  16. OGRFeatureDefn *poFDefn = poLayer->GetLayerDefn();          
  17. int iField;          
  18. for ( iField = 0; iField < poFDefn->GetFieldCount(); iField++ )   
  19. {       
  20. //获取某一个字段的信息              
  21. OGRFieldDefn *poFieldDefn = poFDefn->GetFieldDefn( iField );            if( poFieldDefn->GetType() == OFTInteger )                  
  22. printf( “%d,”, poFeature->GetFieldAsInteger( iField ) );          
  23. }          
  24. OGRGeometry *poGeometry;          
  25. //获取feature          
  26. poGeometry = poFeature->GetGeometryRef();   
  27. //用wkbFlatten宏把wkbPoint25D类型转换为wkbPoint类型          
  28. if( poGeometry != NULL             && wkbFlatten(poGeometry->getGeometryType()) == wkbPoint )     
  29. {              
  30. OGRPoint *poPoint = (OGRPoint *) poGeometry;              
  31. printf( “%.3f,%3.f\n”, poPoint->getX(), poPoint->getY() );          
  32. }          
  33. else              
  34. printf( “no point geometry\n” );          
  35. //销毁feature          
  36. OGRFeature::DestroyFeature( poFeature );      
  37. }      
  38. //销毁数据源,以便关闭矢量文件      
  39. OGRDataSource::DestroyDataSource( poDS );  
  40. }  


利用GDAL/OGR库,可以提供对矢量和栅格文件数据的支持。

转载自:https://blog.csdn.net/L_Andy/article/details/50434524

You may also like...