Python GDAL在GIS数据处理中的应用.docx

GIS的数据可以有多存在方式，可以将数据以某种方式存成TXT，EXCEL等，GIS数据的表现形式往往是衣服很漂亮的地图，而数据也是GIS的核心， GIS的数据往往是很多的，我们在做数据处理的时候往往希望能够批量处理……而数据处理从某种意义上来说就是按照一定的规则和方法对其进行加工，通过这样的操作使得这些数据成为符合我们项目需求的成果。数据处理中经常遇到的就是数据转换，GIS的数据格式太多了……Python作为一门很强大的开源的动态语言，以其灵活型和间接性征服了很多人，GDAL将空间数据抽象，用一个抽象的模型来表达空间数据的实体，同时它提供了很多操作数据的函数，OGR作为GDAL的分支提供了很多对矢量数据的操作。这些心得是我在学习Python和GDAL的时候所做的一些记录，希望这些能给大家学习上或工作中带来帮助，这些心得分对矢量数据和栅格数据的处理做了简单的介绍，也通过两个完整的例子，说明了GDAL的强大。语言是一中工具，用于将人的大脑中的想法，告诉给计算机的工具。一门语言要学精通是很难的，我常常在和朋友交流的时候说，语言，了解了它的特点，语法，能简单的写出几个像样的函数，至于更深入的，用的时候借助GOOGLE等搜索引擎。 GDAL作为空间数据描述处理的工具，对它提供的函数，我只是用到了其中的冰山一角，我也和很多人一样，是初学者，也需要不断的努力，这篇心得中的两个例子分别对矢量数据的处理和栅格数据的处理做了介绍，这两个例子并不是多么的深奥，只是为了体验，就像我前面说的，发现这种数据处理的规则，然后借助这样的工具，帮助我们取得学习上或者工作上的胜利…… 在接触一个新的东西的时候，我喜欢从这新东西的一些名词和结构入手。记得以前一位老师说过，一个新的东西，肯定会有自己的专业术语，然后就是体系结构等……，我们就从这个开始吧GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。OGR是GDAL项目的一个分支，功能与GDAL类似，只不过它提供对矢量数据的支持。有很多著名的GIS类产品都使用了GDAL/OGR库，包括ESRI的ArcGIS，Google Earth和跨平台的GRASS GIS系统。利用GDAL/OGR库，可以使基于Linux的地理空间数据管理系统提供对矢量和栅格文件数据的支持。GDALGDAL（Geospatial Data Abstraction Library）提供对多种栅格数据的支持，包括Arc/Info ASCII Grid(asc)，GeoTiff (tiff)，Erdas Imagine Images(img)，ASCII DEM(dem) 等格式。GDAL使用抽象数据模型(Abstract Data Model)来解析它所支持的数据格式，抽象数据模型包括数据集(Dataset)，坐标系统，仿射地理坐标转换(Affine Geo Transform)， 大地控制点(GCPs)， 元数据(Metadata)，栅格波段(Raster Band)，颜色表(Color Table)，子数据集域(Subdatasets Domain)，图像结构域(Image_Structure Domain)，XML域(XML:Domains)。详细的结构描述请访问/gdal_datamodel.htmlGDAL的核心类结构设计如图所示： 其中的类说明如下： GDALMajorObject类：带有元数据的对象。 GDALDdataset类：通常是从一个栅格文件中提取的相关联的栅格波段集合和这些波段的元数据；GDALDdataset也负责所有栅格波段的地理坐标转换(georeferencing transform)和坐标系定义。? GDALDriver类：文件格式驱动类，GDAL会为每一个所支持的文件格式创建一个该类的实体，来管理该文件格式。 GDALDriverManager类：文件格式驱动管理类，用来管理GDALDriver类。OGROGR提供对矢量数据格式的读写支持，它所支持的文件格式包括：ESRI Shapefiles， S-57， SDTS， PostGIS，Oracle Spatial， Mapinfo mid/mif ， Mapinfo TAB。更详细的请访问/ogr/ogr_formats.htmlOGR包括如下几部分：Geometry：类Geometry (包括OGRGeometry等类)封装了OpenGIS的矢量数据模型，并提供了一些几何操作，WKB(Well Knows Binary)和WKT(Well Known Text)格式之间的相互转换，以及空