GDAL源码剖析与开发指南(txt+pdf+epub+mobi电子书下载)


发布时间:2020-07-10 00:01:28

点击下载

作者:李民录

出版社:信息技术第一出版分社

格式: AZW3, DOCX, EPUB, MOBI, PDF, TXT

GDAL源码剖析与开发指南

GDAL源码剖析与开发指南试读:

前言

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

由于GDAL库支持很多数据格式,目前几乎所有的GIS和RS软件底层都使用GDAL来读写空间数据。

目前国内外关于GDAL的图书非常少,比较权威的只有GDAL官网的说明文档,导致用户缺少实用的参考资料来学习GDAL。针对这一现状,本书对GDAL库进行了全面的介绍,配以详尽的实例说明,来帮助读者尽快掌握和领会GDAL强大的功能。

关于本书

本书共9章和一个附录,分为三部分:第一部分为GDAL基础说明,对应第1章到第5章;第二部分为GDAL高级使用说明,对应第6章到第8章;第三部分为GDAL库提供的工具使用说明,对应第9章。

第1章对GDAL进行了概述,主要介绍GDAL的概念、特点、支持的数据格式等,对GDAL库的编译进行了说明(Windows平台下),同时介绍如何使用SWIG工具编译生成C#、Java、Python的调用接口以及如何使用Doxygen工具生成GDAL库的帮助文档等,最后分别使用四种编程语言(C++、C#、Python和Java)编写了一个简单的、使用GDAL读取图像信息的例子。

第2章对GDAL库中的空间参考类OGRSpatialReference和空间坐标转换类OGRCoordinateTransformation进行说明,并且分别使用四种编程语言(C++、C#、Python和Java)编写了一个简单的坐标转换的例子。

第3章对GDAL库中的OGR库进行说明,包括OGR库中对矢量数据读取和写入的相关类使用说明、OGR库支持的SQL查询说明,最后用四种编程语言(C++、C#、Python和Java)分别编写了一个简单的使用OGR读写矢量数据的例子。

第4章对GDAL库中读写栅格相关的类进行说明,包括GDAL库中对于栅格数据各个部分所对应的类,着重阐述栅格图像中的读写接口 RasterIO,此外对图像的金字塔、统计值信息也进行了说明,最后分别使用四种编程语言(C++、C#、Python和Java)编写了一个简单的读写栅格数据的例子。

第5章对 GDAL 库中支持的常用数据格式进行说明,也包括对栅格数据的创建选项和GDAL库中通用的配置项的说明。

第6章主要对GDAL库目前不支持的数据格式进行扩展,使GDAL可以支持新的数据格式,分为栅格数据扩展和矢量数据扩展两部分。

第7章主要对GDAL库中提供的算法进行说明,基本上涉及到遥感图像处理中的各个流程。每个算法都有完整的示例代码和算法原理说明。

第8章是对GDAL库中所使用的通用库CPL进行说明,涉及的内容主要有:GDAL库中常量和常用宏定义、错误处理相关函数、文件系统相关、文件读写、字符串及其转换、ZIP文件读写、HTTP网络文件访问和XML文件读写等。

第9章对GDAL库自带的命令行工具和Python脚本工具进行说明,这些命令行工具和脚本工具基本上可以完成大部分的栅格数据和矢量数据处理。

关于读者

本书适合地理信息系统(GIS)和遥感(RS)等相关的专业人员阅读,尤其适合开发GIS或RS相关的系统的人员,如果会C/C++的话更好不过了。由于GDAL库是一个专业性比较强的库,需要使用者对空间数据有一定的了解,否则对于GDAL库中空间坐标相关的内容会感到难以理解。

本书中大部分的示例代码都是使用 C/C++语言编写而成的,相比其他的语言,唯一的差别可能就在语法和GDAL库的部分接口上,其实大部分的思路都是一致的。

致谢

本书起源于作者在CSDN网站的博客http://blog.csdn.net/liminlu0314中关于GDAL库的一系列文章,之后在人民邮电出版社陈冀康编辑的联系下才有了出版的想法。

本书在编写过程中得到了很多同学和朋友的帮助和指点,提出了宝贵的意见,他们是王海、刘巍、张艮中、姜言、朱桂海、解云飞、吴炜、HONY、ミ西风゛、narut、小土豆~~、未来、古城、Don、Sunny、图比特-么玉顺、蔚蓝♪libert、流浪的关公、★如履薄冰★等;最后要特别感谢女友孙姝娟对所有的初稿进行了审阅和修订工作,并且在编写过程中给予我鼓励和支持。

父母的养育才有了我的今天,本书的出版也是父母用心教育的结果,养育之恩,没齿难忘。

借此书出版之际,谨向为此书出版做出贡献的同学和朋友们表示衷心的感谢。

如果在阅读过程中有任何疑问,读者可以在作者的CSDN博客中留言或者直接发邮件与作者联系。作者的电子邮箱是liminlu0314@163.com,博客地址为http://blog.csdn.net/liminlu0314。

由于本人水平有限,书中错误和不妥之处在所难免,恳请读者不吝批评指正。

CSDN网站对作者的专访

1.介绍自己学习和工作经历,曾开发过什么产品?之前的工作经验对你现在的工作有何帮助?

2008年7月毕业于中国矿业大学地理信息系统专业。

2008年~2009年年底任职于北京东方泰坦科技股份有限公司。

2010年~2010年10月任职于北京数字空间科技有限公司。

2010年10月至今任职于二十一世纪空间技术股份有限公司。

在这三家公司都是做遥感图像处理相关的算法和软件研发。

之前的工作经验我觉得是我现在的基石,我工作四年多都是在做遥感图像处理的相关工作,很少涉及到其他的行业。这点对我来说既是优点也是缺点。优点就是你可以对遥感图像处理的流程及其原理有比较深入的理解,缺点就是对于其他行业的知识面有些窄。一个是纵向深展,一个是横向拓展。

2.你何时开始接触计算机?是什么原因促使你走上编程的道路?最先使用的语言是哪个?

第一次接触计算机是在2001年,是去同学家接触到的,第一眼觉得很好奇,这就是传说中的电脑。之后高考填报志愿失误,调配到地理信息系统(GIS)专业,这个专业在当时(2004年)算是一个新兴的专业,简单的说就是用计算机来处理和管理地理信息的一种技术。这也就促使我走上了编程的道路。

最先使用的语言应该算是VB,当时大一学校开设的公共基础课就有VB编程。在之后学习专业课中就有C++,此后写程序就全部使用C++语言。在2010年开始接触了Python语言,这是一门解释性语言。之后在写一些简单的工具的时候使用Python,由于它不需要编译,只要安装一个 Python 的运行环境就可以,相比 C++来说,写一些小的工具还是很方便的。

3.你目前工作主要以 C++为主,从长期的使用经验来看,你是如何看待 C++语言的未来,以及C++与C#的关系?

个人觉得C++在很长的时间中还是主流的开发语言,尤其是对于算法要求较高的行业,比如我从事的遥感图像处理。随着卫星的分辨率越来越高,数据量也越来越大,只靠提升硬件来增加处理能力只能提升部分,所以还是需要从算法的角度来进行优化。从语言的执行效率来讲,C++比C#要高很多。但是C#对于界面的开发有着很大的优势,这点对于C++来说确实是一个短板。

所以,我觉得 C++和 C#的关系应该是互存互惠的一个关系,各司其职,C++开发一些对算法要求比较高的系统,而C#应该是朝着快速的这个方向来发展。

4.你怎么看待面向对象技术的发展前景?

面向对象技术是软件设计的一个主流方向,从以前的面向过程到现在的面向对象,使人们处理问题能够更符合人的思维逻辑。同时面向对象的方法可以有效的提升开发效率以及后续的维护和二次开发。

5.C++入门难,你怎么看待这个问题?

C++入门难,个人感觉不存在这个问题。现在很多人觉得C++入门难,我觉得是因为现在网络上很多文章都在说C++难学,这样就无意识的给这些想学习C++的人灌输了一个概念,那就是C++难学。如果事先没有这个概念,那么学习C++和学习其他的语言是一样的。

对于 C++难学可能还有另外一个方面的原因,就是 C++本身只是提供了一套语法规则,没有像 C#和 Java 之类的其他语言那么多自带的库,比如正则表达式等。如果用 C++实现类似的功能确实很复杂,但是比较幸运的是,目前有很多开源的 C++库提供了这些功能,比如被称为C++的准标准库的Boost库,里面提供了一百多个涵盖了各个方面的库。

任何语言的学习都离不开一个过程,那就是亲自写代码,调试程序,在调试过程中理解语法的意义和概念。

此外,对于学习C++,建议大家经常去看一些科技网站,关注一些C++的开源库,这样很多时候就避免了重复发明轮子的这个过程。而且很多的开源库都有很多商业软件在使用,在效率和可靠性上都是比较好的。所以多了解一些C++的开源库对自己的工作肯定有很大的帮助。

6.从博客中看到,你正在写关于GDAL源码剖析的系列文章。你最初的想法是什么?透露下,未来一个月你将写哪些内容?

当初我学习GDAL的时候(2008年左右),关于GDAL的学习资料很少,除了官方的API帮助文档之外,就很少了。慢慢的,发现积累的东西越来越多,但是没有记录,忘记的也快,为了以后用到的时候方便查阅,于是就开始写关于GDAL的博客。之后创建了一个GDAL的交流群,不到一月的时间,群就满了。很多人在我的博客下面留言,说是帮助他解决了问题,或者问我一些类似的问题。所以我觉得有必要把我学习到的关于GDAL的一些知识跟大家分享一下,独乐乐不如众乐乐,于是就整理了一下,形成现在的GDAL源码剖析的系列的文章。

接下来就是继续更新《GDAL源码剖析》这个系列的内容,之前主要写的都是关于GDAL的使用方面的内容,接下来想把GDAL的内部机制和实现原理进行一个说明,真正对得起这个系列的名字,好让广大的GDAL使用者能够知其然并且知其所以然。

7.GDAL在读取不同的数据时,会发生什么常见问题?有何方法能快速的解决问题?

GDAL在读取不同的数据时,最常见的问题可能有下面三个:

1、GDAL关于中文路径的问题,以及矢量数据属性中文乱码的问题;

2、对于C#语言调用会出现程序崩溃的问题;

3、关于RasterIO读写图像的问题。

这三个问题是在群里提问最多的,对于第一个问题,主要是由于GDAL在1.8之后的版本,增加了对于多字节的处理,默认的路径编码改成UTF-8的编码导致的。由于C++中文版默认的编码是GBK或者其他的,一般有两种方式解决,第一将中文路径编码专为UTF-8,第二设置GDAL的默认编码不是UTF-8即可解决。第二个问题最主要的原因有两个:一是GDAL的dll没有复制到程序的目录;二是由于编译GDAL的平台版本和C#的平台不一致导致,比如GDAL编译的是32系统的,C#在编译的时候必须选择X86平台,否则就会出现这个问题。第三个问题主要是对于RasterIO这个函数的参数没有理解清楚导致的,尤其是对于最后的三个参数。

关于这三个问题,我分别写了几篇博客,对其进行了说明。

8.现在你从事的工作主要以C++和GDAL库为主,两者搭配的优势和弊端分别是什么?

GDAL 库本身就是用 C/C++语言编写,所以对于 C++的使用,本身没有语言的差异,使用起来也非常的方便。通过最基本的包含头文件,引用库文件就可以使用。

C++和 GDAL 搭配的弊端主要还是 C++语法的一些限制吧,比如读取数据的时候内存的申请和释放等。对于这个问题只能从使用人员本身来进行避免。

同时GDAL使用SWIG来提供其他语言的API接口,比如C#、Java、Python、PHP等。所以在使用GDAL的时候,对于搭配的语言来说,根据个人的习惯来决定。

9.分享下你学习C++和GDAL的一些经验,有何好书籍或方法可以推荐给初学者?

学习C++的经验就是自己亲手写代码、调试,这样才能真正提高自己的能力,光看书是不行的,那就成了纸上谈兵,到真正用的时候就发现写不出来了。当然了,书还是必须要看的。我觉得学习任何语言都脱离不了这个规则吧。关于C++的书籍,《C++ Primer》是一本好书,不管是学习还是用来当作参考书都是很不错的选择。另外还有一本就是《深入理解计算机系统》,这本书对于学习C/C++的人来说是一本非常好的书,尤其是想理解C/C++程序是如何映射到系统上,以及程序是如何执行的,读者能够更好地理解程序的行为为什么是这样的,以及效率低下是如何造成的。这本书通过程序员的视角,让读者可以深深地体会到学习计算机系统的内部工作原理会对他们今后的工作有进一步的帮助。

对于学习GDAL的方法,第一看GDAL的官方文档。我曾经把GDAL的官方文档全部看完过,当时还准备翻译成中文的文档,后来由于种种原因没有完成。通过官方文档和官方提供的入门Demo,逐个理解每个类,每个函数的意义。或者先熟悉一下每个类里面的接口,有个印象,等用到的时候可以大致知道,然后再去文档里面找说明。

此外,学习 GDAL 首先要对遥感图像数据的基本信息有一个了解,然后才能更好的了解GDAL 库,比如遥感图像同普通的图像多了空间信息,这些空间信息主要就是投影信息和分辨率。有了这些概念,就可以很容易的和GDAL的API结合起来,比如GDALDataset类中的GetProjectionRef这个函数就是获取的遥感图像的投影信息。

由于GDAL是在2007年前后才开始广泛的应用,所以目前还没有出版的参考书来进行参考,只能去参考它的官方帮助文档。此外我整理了一个GDAL专栏,基本上也可以作为中文的一个参考文档。

10.你在学习或工作中,是怎么接触到CSDN?CSDN对于你的工作或学习有什么影响,起到过什么帮助?有什么故事可以分享?

第一次接触 CSDN 是在大学 C++的课堂上,老师推荐有问题可以去 CSDN,之后遇到问题,首先就是先看看CSDN有没有别人遇到同样或者此类的问题(基本上80%都可以找到解决问题的方式)。

很感谢CSDN提供的这样一个平台,可以让大家集合起来,尤其是对于CSDN的论坛和博客,这两个板块个人非常喜欢。在论坛可以得到别人的帮助和帮助别人,博客可以和大家分享自己的一些学习经验等。

11.你对CSDN有什么建议,以及你对CSDN的未来有什么期待?

对于 CSDN,希望可以在发布博客的时候提供客户端的工具,之前可以使用 Live Writer来发布,但是之后又不行了。对于有时候发布一些图片比较多的博客使用在线的编辑工具确实不太方便,希望能够提供一个客户端的工具。

最后祝CSDN的发展越来越好!第1章GDAL简介1.1什么是GDAL

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

OGR(OGR Simple Features Library)是GDAL 项目的一个分支,功能与GDAL 类似,只是它提供对矢量数据的读写支持。同时它实现了一个对空间参考信息进行处理的类,用来对空间数据的空间信息进行处理。

GDAL 最初是由Frank Warmerdam 于1998 年开始开发的,在GDAL1.3.2 版本之后,正式由开源空间信息基金会(OpenSourceGeospatialFoundation,简称OSGeo)下的GDAL/OGR项目管理委员会对其进行维护。

很多著名的GIS类产品都使用了GDAL/OGR库,包括ESRI的ArcGIS系列、Erdas2011、Google Earth和跨平台的GRASS GIS、Quantunm GIS系统等(更多使用GDAL 的软件和产品,参考GDAL网址http://trac.osgeo.org/gdal/wiki/SoftwareUsingGdal)。1.2GDAL特点

GDAL 提供对多种栅格数据的支持,包括Arc/Info ASCII Grid(asc),GeoTiff (tiff),Erdas Imagine Images(img),ASCII DEM(dem) 等格式。

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

GDAL 库还提供了一系列算法接口,比如矢量栅格化、栅格矢量化、图像校正算法以及DEM相关的算法接口等;同时对这些算法提供了一系列可以运行的程序或者脚本,方便人们使用。1.3GDAL支持的数据格式

GDAL支持的数据格式分为两种:栅格数据格式和矢量数据格式。GDAL在每次版本更新的时候,都会加入对新的格式的支持。本书使用 GDAL1.10 版本中支持的栅格数据格式,如表1-1所示,支持的矢量数据格式如表1-2所示。表1-1 GDAL支持的栅格数据格式(续表)(续表)(续表)(续表)(续表)(续表)(续表)

注意

文件的最大存储大小不仅由文件的格式决定,而且和操作系统有关系。表1-2 GDAL支持的矢量数据格式(续表)(续表)(续表)(续表)1.4GDAL源码下载

GDAL 的源码下载有两种方式。第一种是直接下载打包好的压缩包,下载地址为:http://trac.osgeo.org/gdal/wiki/DownloadSource。另外,GDAL的各个历史版本都可以在这个地址进行下载:http://download.osgeo.org/gdal/。这种下载方式很简单,此处不作说明。下载完成,使用压缩软件解压即可。

第二种下载方式是使用 SVN 源代码管理工具直接从 GDAL 的源码服务器进行下载。下载地址是:http://svn.osgeo.org/gdal/trunk。从 SVN 下载的好处是可以直接看到代码的提交者对代码做的修改信息等。需要注意的是,从这里下载的源代码不是发布的版本,而是开发版本,可能会有些bug。下面以TortoiseSVN为例对使用SVN下载做一个简单的说明,使用SVN下载之前要确保已经安装了SVN的客户端软件。

1)在电脑中新建一个文件夹用来存放GDAL的源代码。

2)在该文件夹点击右键,在弹出菜单中选择“SVN check out…”,弹出如图1-1所示的对话框。

3)对话框中的“URL of repository”中填入“http://svn.osgeo.org/gdal/trunk”,在“Checkout directory”中确定迁出的文件夹,其他选项默认,点击“OK”按钮,即开始下载GDAL的源代码。图1-1 SVN checkout 对话框1.5GDAL源码目录

下载的GDAL源代码压缩包目录如图1-2所示,使用SVN迁出的GDAL源代码目录如图1-3所示。图1-2 GDAL 源码压缩包文件目录结构

从图1-2和图1-3所示可以看出,不管用什么方式获取GDAL的源代码,它的目录结构都是一样的。下面我们就针对目录结构中的每个文件夹和文件作一个简单的说明(按照字母顺序)。图1-3 GDAL SVN源码文件目录结构1.5.1 文件夹说明

对于图1-3中所示的文件夹,详细说明如下。

alg:该文件夹存放的是GDAL库中提供的一些算法的源代码,这些算法包括但不限于:DEM生成等高线算法,图像纠正算法(几何纠正、TPS纠正、RPC纠正等),栅格矢量化算法,矢量栅格化算法,格网计算算法,PCT和RGB互转算法,分类图的小碎斑块去除算法等。

apps:该文件夹中存放的是GDAL库中提供的一些命令行工具集的源代码,这些工具集的介绍可以参考http://gdal.org/gdal_utilities.html。后面会对这些工具做一个简单的说明,其中有些工具非常有用,比如gdalinfo,可以用来查看图像的元数据信息等。

bridge:该文件夹中存放的是用来连接GDAL抽象类的定义与GDAL自己的结构体定义及其实现的源代码。后面涉及GDAL的实现原理时我们会对该文件夹作一个比较详细的介绍。

data:该文件夹中存放的是GDAL库中需要用到的一些“配置文件”(此处叫配置文件可能不太准确),这些文件主要有ESRI的投影文件,ESPG的投影文件,PCI的投影和椭球体文件,AutoCAD 的DXF 格式的头文件等。在使用GDAL 库的过程中可以在程序中使用函数CPLSetConfigOption("GDAL_DATA","C:\GDAL\data");来设置该文件夹的目录,如果没有设置GDAL会自动从环境变量中查找,环境变量的名字叫GDAL_DATA,该变量的值就是data文件夹的路径;如果还是没有找到,那么GDAL可能会提示错误,比如在写入AtuoCAD的dxf格式的时候如果不设置GDAL_DATA,那么就会提示创建不成功,后面遇到时我们会再进行说明。

doc:该文件夹存放的是用来生成 GDAL 帮助文档的一些 dox 文件。dox 文件是使用doxygen工具来进行生成的,后面我们会对doxygen工具作一个简单的介绍,并且说明在自己的工程中怎样使用doxygen生成程序的开发帮助文档。

frmts:这个文件夹可以说是GDAL代码中内容最多的一个文件夹了,每次更新GDAL的版本后这个文件夹中都会多出几个文件夹,同时在GDAL支持的文件格式中也会多出来几个新的文件格式。没错,这个文件夹存放的就是GDAL针对不同的特定图像格式解析的源代码,我们可以举几个简单的例子,比如bmp文件夹就是解析BMP图像的,hfa文件夹是用来解析Erdas的img图像格式,还有pcidsk文件夹是读取PCI的pix格式的文件等。所以这个文件夹存放的是解析各个文件格式的源代码。

gcore:这个文件夹是 GDAL 的灵魂所在,主要存放 GDAL 抽象类的数据集,图像、波段的读写接口等都是在这里面实现的。如果想知道GDAL的抽象类是怎么对图像格式进行抽象的,读者可以看看这个文件夹里面的代码。

html:这个文件夹主要用来存放GDAL生成的帮助文档,使用doxygen工具生成的GDAL帮助文档会出现在这个文件夹中。如果使用压缩包的话,该文件夹里面应该是空的。后面我们会和doc文件夹一起进行详细介绍。

m4:该文件夹存放的是以m4 为后缀名的文件。m4 文件叫Macro Processor Library,是编译基础中最核心的文件,这个文件主要是用autoconf来产生configure配置文件,继而自动生成Makefile文件。这个文件夹在Windows平台下没用,在此就不作细述了。

man:该文件夹是用来生成Linux或者其他平台下的帮助文件,Windows平台下也没用。

ogr:GDAL和OGR曾经是两个库,GDAL库负责读写栅格数据,OGR库负责读写矢量数据,然而GDAL的算法库中经常会用到矢量数据的读取,将这两个库分开会造成不便,现在这两个库已经整合到一起,OGR库就是GDAL库的一个子集。其实OGR库还是可以单独编译出来的。ogr文件夹就是存放OGR库源代码的文件夹。后面我们会对这个文件夹里的内容进行详细介绍。

port:该文件夹中存放的是CPL库的东西,CPL是一个GDAL的底层支持库,CPL库中定义了一些字符串的操作、文件处理、网页请求、数据库连接、哈希表、字符加密文件压缩等基础的函数。比如GDAL中所有的导出函数符号CPL_DLL就是在这个CPL文件夹中定义的;还有frmts文件夹中,打开文件、打开数据库、打开网络路径以及字符串的处理等都是用的port库。

swig:该文件夹主要是存放swig 的脚本。SWIG 全称叫Simplified Wrapper and InterfaceGenerator,官方网站是 http://www.swig.org/。SWIG 的作用就是可以将 C/C++写的库封装为Python、C#、Java、Perl和Ruby等其他语言的访问接口。GDAL的C#版本就是使用swig来编译出来实现的。

vb6:这个文件夹用来将GDAL编译成一个VB6的模块,对于里面具体的文件说明以及如何编译参考文件夹中的readme.txt,这里就不作具体说明了。

wince:顾名思义,这个文件夹中的文件就是用来编译Windows CE 平台下的GDAL 库的,具体编译请参考其中的说明文档。1.5.2 文件说明

aclocal.m4:同1.5.1小节介绍的m4文件夹。

autogen.sh:Linux平台下的shell文件,通过调用autoconf来产生configure配置文件。

COMMITERS:该文件中的内容是GDAL开发人员的信息,包括姓名、联系邮箱以及各自负责开发的模块说明等。

config.guess、config.sub、configure、configure.in:这四个文件在Windows下没有具体作用,在此略过。

Doxyfile:Doxyfile就是前面doc文件夹说明中提到的doxygen的工程文件,用来生成帮助文档,后面我们在介绍doxygen的使用时会对该文件进行说明。

GDALmake.opt.in:这个文件是 Linux 平台下的 GDAL 库编译配置文件,在后面的nmake.opt中我们会介绍其功能。

gdalnightlysvn.sh:Linux平台下调用svn获取GDAL源代码的一个shell脚本。

GNUmakefile:GNU的make文件。

HOWTO-RELEASE:GDAL发布版本的一些说明。

install-sh:Linux平台下GDAL的安装shell脚本。

LICENSE.TXT:GDAL的许可说明文件。

ltmain.sh:Linux平台下libtool的shell脚本,Windows下没用到。

makefile.vc:GDAL的编译文件,用来将源代码编译成dll文件,后面我们会在GDAL编译中作进一步说明。

makegdal_gen.bat:一个用来生成VS的工程文件的批处理文件,后面在GDAL编译中我们会对该文件的使用方式进行说明。

makegdal10.sln:文件夹中所有的 sln 文件都是 VS 的项目文件,文件名后的数字代表的是VS的版本号。

makegdal10.vcxproj、makegdal71.vcproj:VS的工程文件,该文件可以由makegdal_gen.bat文件自动生成,后面我们会详细介绍。

mkbindist.sh、mkgdaldist.sh、mktestdist.sh:三个shell脚本文件,Windows下没用,略过。

NEWS:介绍GDAL的新增功能以及修复的bug记录等。

nmake.opt:GDAL编译选项配置文件,在编译GDAL中,可以指定GDAL绑定的其他库等都在这里进行设置。在后面的GDAL编译中我们会详细说明。

nmake-wince.opt:编译wince版本的编译选项配置文件。

PROVENANCE.TXT:GDAL目录说明文件,如果上面说明的不够清楚,可以参考这个文件。

submake.bat:一个编译的批处理文件。

svnkeywords.sh:SVN的又一个shell脚本。

VERSION:GDAL版本信息。1.6GDAL编译

GDAL 的编译分为三个小节,第一节是常规编译,使用微软的Visual Studio 2008进行编译,如果用微软的其他编译器,编译过程与此类似。需要注意的是,GDAL已经不再支持VisualC++6.0编译器。第二小节主要说明怎样修改GDAL的配置文件,以使GDAL支持更多的数据格式以及算法等。第三小节对编译GDAL的64位版本进行说明。1.6.1 常规编译

1.使用Visual Studio 2008 IDE 编译

在 GDAL 的源代码目录中,有几个以 sln 为后缀的文件,比如 makegdal10.sln、makegdal80.sln、makegdal71.sln、makegdal90.sln 等。这些是 VisualStudio 的工程文件,文件名中的数字表示VS的版本号:71表示VS2003、80表示VS2005、90表示VS2008、10表示VS2010。我们根据自己电脑安装的VS版本,打开对应的文件,如图1-4所示(使用VS2008SP1版本,打开makegdal90.sln文件)。图1-4 GDAL 源码压缩包文件目录结构

然后我们在左侧解决方案点击右键,在弹出菜单中点击“生成”或者“重新生成”命令, GDAL 就会开始编译。输出窗口中会有编译过程中的一些信息提示,直到出现生成成功等信息,就表示GDAL已经完成编译,这时在GDAL的源代码目录中会出现gdal.lib、gdal_i.lib、gdal18.dll等文件。如果没有修改GDAL中的nmake.opt文件的话,那么C盘中会出现一个名为“C:\warmerda\bld”的文件夹,里面包含三个文件夹,分别是bin、data和html。bin文件夹中存放的是编译出来的GDAL的可执行程序,包括GDAL提供的十几个工具集;data文件夹就是1.5.1小节中介绍的data文件夹;html文件夹中存放的是各种数据格式的说明文档。

2.使用Visual Studio 2008 命令行编译

使用cmd命令行编译,首先在开始菜单中找到“所有程序→Microsoft Visual Studio 2008→Visual Studio Tools→Visual Studio 2008 命令提示”,点击“Visual Studio 2008 命令提示”弹出如图1-5所示的界面。图1-5 Visual Studio 2008命令提示

使用cd命令,切换到GDAL的源代码目录,如图1-6所示。图1-6 切换至GDAL 源码目录

切换至GDAL的源代码目录后,依次敲入下面的命令行后回车,等待编译结束即可。

nmake -f makefile.vc

nmake -f makefile.vc install

nmake -f makefile.vc devinstall

同时还有其他的命令,如下所示。

nmake -f makefile.vc clean

nmake -f makefile.vcMSVC_VER=1400 clean

nmake -f makefile.vcMSVC_VER=1400 DEBUG = 1

以上六行的命令含义依次如下。

编译GDAL库;

编译GDAL库并安装(这里安装的意思是将生成的dll、exe等文件复制到C:\warmerda\bld目录);

编译GDAL库并安装开发者模式(安装的意思同上,开发者模式意思是将开发用的include文件夹中的头文件和lib文件一同复制到C:\warmerda\bld目录,此时会在C:\warmerda\bld目录中增加两个文件夹,分别是include和lib,分别存放GDAL的头文件和lib文件,用于调用GDAL库);

清理GDAL库,同时会删除编译GDAL库所生成的临时文件,作用相当于VS环境中的清理命令;

作用与第四行相同,但是添加了一个“MSVC_VER = 1400”,表示使用VS2005 编译;

编译GDAL库的debug模式,可以用来调试GDAL源码。1.6.2 自定义编译

GDAL 的强大之处不单单在于可以读取栅格和矢量数据,同时还表现在以下三个方面:第一,可以进行矢量图形之间的一些常用操作,比如求交、求并、缓冲区等;第二,可以进行投影和坐标转换;第三,可以使用配置文件来对GDAL库支持的数据格式进行扩展。

如果使用GDAL默认的编译方式,那么上述的前两个功能您将不能使用,因为GDAL这两大功能是基于另外的两个开源库GEOS(Geometry Engine, Open Source)库和PROJ4 库来实现的。下面对这两大库分别做一个简单的说明,讲述如何修改编译文件,让GDAL能够拥有这两大功能。

1.集成GEOS

简单来说,GEOS提供了OGC规范中简单几何要素对象操作的C++语言的实现。更详细的说明可以在GEOS的官网http://geos.osgeo.org中找到。

GEOS 库提供了常用的空间操作,比如判断两个几何对象相互关系的函数 Contains()、Intersects()、Touches()和Crosses()等,计算两个几何对象的空间操作的函数Difference()、Union()以及Buffer()等。

我们假设已经有编译好的GEOS库,首先使用记事本或者其他的文本编辑器打开GDAL源代码目录下的nmake.opt 文件,找到“# Uncomment for GEOS support”这句(1.10版本大概在490行左右),如下面三行代码所示。

#GEOS_DIR = C:/warmerda/geos

#GEOS_CFLAGS = -I$(GEOS_DIR)/capi -I$(GEOS_DIR)/source/headers -DHAVE_GEOS

#GEOS_LIB = $(GEOS_DIR)/source/geos_c_i.lib

我们将它们修改为以下代码。

GEOS_DIR = F:\Work\3rdPart\geos-3.2.2

GEOS_CFLAGS = -I$(GEOS_DIR)/capi -I$(GEOS_DIR)/source/headers -DHAVE_GEOS

GEOS_LIB = $(GEOS_DIR)/source/geos_c_i.lib

代码中F:\Work\3rdPart\geos-3.2.2是GEOS存放的主目录,后面两行设置的是GEOS的头文件目录和lib文件路径,设置好后保存即可。对比结果如图1-7所示。图1-7 GEOS配置修改

保存完nmake.opt之后,我们按照1.6.1小节中的编译方式进行编译即可。编译后的GDAL将会支持图形之间的操作,因为函数主要是在OGR库中,后面会在OGR库中进行详细说明。编译完成后,我们需将geos_c.dll文件复制到gdal18.dll的同级目录下,否则会提示找不到geos_c.dll文件。

2.集成PROJ.4

PROJ.4是一套开源的坐标投影转换库,它可以完成在两套不同制图投影系统之间的转换,或者不同的椭球体或大地基准面之间的转换。GDAL中用到的坐标转换、投影转换、几何纠正、正射纠正等算法,都需要 PROJ.4 库的支持。PROJ.4 的配置同 GEOS 库的配置方法相同,在nmake.opt文件中找到PROJ.4库的位置(1.10版本大概在429行左右),如下面的三行代码所示。

#PROJ_FLAGS = -DPROJ_STATIC

#PROJ_INCLUDE = -Id:\projects\proj.4\src

#PROJ_LIBRARY = d:\projects\proj.4\src\proj_i.lib

我们将它们修改为以下代码。

#PROJ_FLAGS = -DPROJ_STATIC

PROJ_INCLUDE = -IF:\Work\3rdPart\proj-4.7.0\src

PROJ_LIBRARY = F:\Work\3rdPart\proj-4.7.0\src\proj_i.lib

代码中第一行表示是否使用静态链接的方式,第二行表示PROJ.4库存放的路径,第三行为PROJ.4库的lib文件所在路径,修改后保存即可,对比结果如图1-8所示。图1-8 PROJ库配置修改

同GEOS库一样,保存完nmake.opt之后,我们按照第一步中的编译方式进行编译即可。编译后,我们同样将proj.dll文件复制到gdal110.dll的同级目录下,否则与坐标转换相关的所有功能都不能使用。

3.集成HDF数据读取

通过上面GEOS和PROJ库的介绍,我们对于GDAL的配置文件nmake.opt已经有了一个初步的了解。下面使用GDAL读取HDF数据也是同样,我们先在HDF的官方网上下载编译好的HDF4和HDF5的库,然后修改配置文件进行编译。HDF官网提供了32位的库和64位的库,这里我们以HDF4.2.6和HDF5-1.8.7两个库为例介绍32位库的编译,对于64位库的编译后面会进行专门介绍。

我们将下载好的 HDF4.2.6 和 HDF5-1.8.7 两个库解压,然后在 nmake.opt 文件中,找到HDF的位置(1.10版本大概在330行左右),代码如下。

# Uncomment the following and update to enable NCSA HDF Release 4 support.

#HDF4_PLUGIN = NO

#HDF4_DIR =D:\warmerda\HDF41r5

#HDF4_LIB =/LIBPATH:$(HDF4_DIR)\lib Ws2_32.lib

# Uncomment the following and update to enable NCSA HDF Release 5 support.

#HDF5_PLUGIN = NO

#HDF5_DIR =c:\warmerda\supportlibs\hdf5\5-164-win

#HDF5_LIB =$(HDF5_DIR)\dll\hdf5dll.lib

我们将它们修改为下面的代码。

# Uncomment the following and update to enable NCSA HDF Release 4 support.

HDF4_PLUGIN = NO

HDF4_DIR =F:\Work\3rdPart\HDF4.2.6_win_x86

HDF4_LIB = $(HDF4_DIR)\dll\hd426m.lib $(HDF4_DIR)\dll\hm426m.lib \

$(HDF4_DIR)\lib\hd426.lib $(HDF4_DIR)\lib\hm426.lib Ws2_32.lib

# Uncomment the following and update to enable NCSA HDF Release 5 support.

HDF5_PLUGIN = NO

HDF5_DIR =F:\Work\3rdPart\HDF5-1.8.7_win_x86

HDF5_LIB =$(HDF5_DIR)\dll\hdf5dll.lib

代码修改前后的对比如图1-9所示。图1-9 HDF库配置修改

保存,然后编译 GDAL 即可,同时我们将 HDF 库中的 dll 文件夹下的 dll 文件复制到gdal110.dll的同级目录下。

按照同样的方式,我们可以通过修改配置文件使GDAL支持NetCDF、Oracle等格式。1.6.3 编译其他方面

1.makegdal_gen.bat使用

下面我们介绍怎样使用makegdal_gen.bat 文件生成Visual Studio 的工程文件。

首先我们打开 cmd 命令行窗口,使用 cd 命令切换到 GDAL 源代码目录,然后输入makegdal_gen.bat后回车,会得到该工具的一个简单实用帮助,如图1-10所示。图1-10 makegdal_gen.bat使用说明

该工具是一个带有命令行参数的批处理工具(在后面对GDAL工具集的介绍中我们会对带有参数的命令行程序以及编写带有命令行的程序作一个详细说明),通过图1-10 可以看出其基本语法如下。

makegdal_gen 7.10 > makegdal71.vcproj

makegdal_gen 8.00 > makegdal80.vcproj

通过上面的示例我们可以看出,该工具的命令行参数分别是:VS的版本号(具体版本号参考1.6.1 小节)、一个大于号“>”、输出的VS的工程名字。如果要使用该命令行生成一个VS2008版本的工程文件,我们可以输入如下命令,回车即可。

makegdal_gen 9.00 > makegdal90.vcproj

2.编译64位系统下的GDAL

GDAL的64位系统的编译,和32位系统的编译基本相同,首先我们在VS的工程中打开配置管理器,然后在活动解决方案平台的下拉列表中选择“新建”,然后弹出“新建解决方案平台”对话框,选择新平台为x64(需要在安装VS时安装64位的编译环境),然后点击“确定”按钮,最后在VS中选择X64进行编译即可,如图1-11所示。图1-11 GDAL 编译64位的工程配置

使用cmd 命令行编译基本同32 位的一样,区别在于在开始菜单选择的不是“Visual Studio 2008 命令提示”,而是“Visual Studio 2008 x64 兼容工具命令提示”。

在编译开始之前,我们还需要打开nmake.opt 文件,找到131 行处的“#WIN64 = YES”,将前面的“#”去掉,保存,然后开始编译。如果就这样编译通过的话,那么恭喜;如果不能顺利编译过去,那么我们需要按照下面的步骤进行设置。

1)我们在GDAL目录下的nmake.opt文件中找到SYM_PREFIX的定义,应该在438行左右;将“SYM_PREFIX = _”改为“SYM_PREFIX = ”,也就是将最后的下划线去掉。

2)在GDAL目录下的makefile.vc文件中,找到46行左右的代码,如下所示。

BASE_INCLUDE = /INCLUDE:_GDALSimpleImageWarp@36 \

/INCLUDE:_GDALReprojectImage@48 \

/INCLUDE:_GDALComputeMedianCutPCT@32 \

/INCLUDE:_GDALDitherRGB2PCT@28 \

/INCLUDE:_OCTNewCoordinateTransformation@8 $(VB6_SAFEARRAYSYM)

我们需要将它修改为如下代码。

BASE_INCLUDE = /INCLUDE:$(SYM_PREFIX)GDALSimpleImageWarp \

/INCLUDE:$(SYM_PREFIX)GDALReprojectImage \

/INCLUDE:$(SYM_PREFIX)GDALComputeMedianCutPCT \

/INCLUDE:$(SYM_PREFIX)GDALDitherRGB2PCT \

/INCLUDE:$(SYM_PREFIX)OCTNewCoordinateTransformation $(VB6_SAFEARRAYSYM)

修改完之后保存,然后重新编译生成GDAL即可。1.7SWIG编译1.7.1 SWIG简介

SWIG全称是Simplified Wrapper and Interface Generator,官方网站是http://www.swig.org/。SWIG是一个开发工具,能够将使用C或者C++编写的软件与其他各种高级编程语言进行嵌入联接。SWIG能应用于各种不同类型的语言,包括常用脚本编译语言例如Perl、PHP、Python、Tcl、Ruby、PHP 等,也包括非脚本编译语言,如C#、Common Lisp (CLISP、Allegro CL、CFFI、UFFI)、Java、Modula-3、OCAML、R等。SWIG最典型的应用就是转换C/C++接口,产生高级语言需要的代码来调用 C/C++代码。SWIG 是一个免费开源软件,可以在商业或非商业中自由使用、分发、修改。

下载安装Swig的时候我们注意要下载Swigwin文件而不是下载源代码。以swigwin-2.0.4为例,我们下载swigwin-2.0.4.zip文件后解压,将解压目录也添加到环境变量Path中去,否则后面需要设置swig的所在目录。检验swig是否成功设置到环境变量Path中的最简单的办法就是在运行中输入swig后回车,如果提示windows找不到swig,那么说明没有设置成功;如果出现一个黑屏一闪而过,那么说明已经设置成功了。

下面三小节需要在安装完SWIG之后才能进行编译。需要注意的是,编译下面三种语言的GDAL之前,我们首先要编译C++版本的,就是1.6节的内容。1.7.2 编译C#版本GDAL

首先,我们打开 nmake.opt 文件,找到 SWIG = swig.exe 这一句,假如没有将 swig 的目录添加到环境变量中,那么将这句后面的 swig.exe 修改为 swig.exe 的全路径,如F:\Work\3rdPart\swigwin-2.0.4\swig.exe;如果设置了环境变量,那么就不需要进行修改了。

然后,我们按照1.6节中的使用cmd命令编译GDAL的方式来进行编译,打开“Visual Studio2008 命令提示”并定位到 GDAL 源代码目录,然后依次执行如下三行命令,其实就是编译C++版本的GDAL。

nmake /f makefile.vc

nmake /f makefile.vc install

nmake /f makefile.vc devinstall

执行以上三句之后,接下来我们使用 cd 命令,进入 swig\csharp 文件夹中并使用 nmake工具编译,命令如下所示。

cd swig\csharp

nmake /f makefile.vc

nmake /f makefile.vc install

执行完以上两句命令后,系统会在csharp文件夹下生成8个dll文件,并将这8个dll文件复制到GDAL输出目录下的csharp文件夹中。

这8个dll分别如下。

gdal_csharp.dll

gdal_wrap.dll

gdalconst_csharp.dll

gdalconst_warp.dll

ogr_csharp.dll

ogr_wrap.dll

osr_csharp.dll

osr_wrap.dll

在C#工程中引用GDAL的时候,我们只需要把上面8个dll中后面以_csharp.dll结尾的四个dll加入到C#工程中即可。需要注意的是,在使用这几个dll的时候我们需要将gdal110.dll以及它依赖的其他dll都复制到同一个目录中,才能使C#程序正常运行。1.7.3 编译Python 版本GDAL

在编译Python版本之前,我们要确保自己的电脑中已经安装了Python。如果我们安装了ArcGIS软件,那么就不用再安装Python了,因为在安装ArcGIS的时候必须要安装Python。将Python的bin目录添加到环境变量Path中去,然后我们就可以开始编译Python版本的GDAL了。

打开nmake.opt文件(1.10版本大概76行左右),将PYDIR后面的路径修改为本机Python的路径。

PYDIR ="C:\Python26"

修改完之后,保存nmake.opt文件。接下来使用cd命令将命令行定位到swig\python目录后,依次输入如下命令。

python setup.py build

python setup.py install

执行完上述命令后,我们会在python的site-packages目录看到多了GDAL和OGR的文件以及一个osgeo的文件夹。我们将编译出来的gdal110.dll以及它依赖的所有DLL都复制到site-packages文件夹中的osgeo 文件夹中,最后就可以在Python 中用imort gdal 来使用GDAL了。测试Python版本的gdal是否安装成功,可以使用如下Python代码。

from osgeo import gdal

from osgeo.gdalconst import *

dataset = gdal.Open("F:\Work\Data\envi.img",GA_ReadOnly)

dataset.GetDriver().ShortName

#'ENVI' #图像格式

dataset.GetProjectionRef()

#'PROJCS["unnamed",GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHERO ID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4269"] ],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",30],PARAMET ER["central_meridian",-84.16666666666667],PARAMETER["scale_factor",0.9999],PAR AMETER["false_easting",2296587.926509186],PARAMETER["false_northing",0],UNIT["Foot (International)",0.3048]]' #图像投影信息1.7.4 编译Java 版本GDAL

在编译Java版本的GDAL时,除SWIG之外,我们还需要JRE和Ant这两个工具,没有的话需要先下载安装,其中 ANT 只需直接解压,不用安装。接下来我们用记事本打开nmake.opt,找到第86行左右位置(1.10版本),将文件中的JAVA_HOME和ANT_HOME修改为本机的路径,只要修改到主文件夹即可(就是bin文件夹所在的文件夹);如果文件夹路径有空格的话,把路径用双引号引起来,修改后的内容与以下代码类似。

# SWIG Java settings

!IFNDEF JAVA_HOME

JAVA_HOME = "C:\Program Files\Java\jdk1.6.0_26"

!ENDIF

!IFNDEF ANT_HOME

ANT_HOME = F:\Work\GDAL\gdal_need\apache-ant-1.8.4

!ENDIF

JAVADOC = $(JAVA_HOME)\bin\javadoc

JAVAC = $(JAVA_HOME)\bin\javac

JAVA = $(JAVA_HOME)\bin\java

试读结束[说明:试读内容隐藏了图片]

下载完整电子书


相关推荐

最新文章


© 2020 txtepub下载