GDAL/OGR 快速入门

　　本教程只涉及命令行操作。对于数据的可视化，请使用 OSGeo-Live 所提供的各款桌面 GIS 系统，例如 Quantum GIS (QGIS) 。

　　本教程分为 GDAL（栅格数据处理）和 OGR（矢量数据处理）两个部分。第一部分为 GDAL 。

　　涉及的内容有：

　　OSGeo Live 的演示文件位于 /usr/local/share/data 。这里使用 Natural Earth data 来进行演示。请先将其拷贝至 home 目录。即：

以下位置将出现一个 NaturalEarth 栅格文件和一个 tfw 文件：

Tip

建议使用 QGIS 等程序浏览数据。

　　首先，应当了解所需的数据驱动。使用 –formats 开关可以显示 gdal_translate 所支持的格式。

　　–format 开关等够显示有关数据驱动的许多信息，包括文件创建操作和允许的格式等。

　　使用 gdal_translate 完成转换操作。默认输出是 GeoTIFF：

　　使用 -of 开关控制输出格式，使用 -co 控制文件创建参数（此处为 jpg 文件质量）：

　　另外，-ot 用于控制输出数据类别（此处为转换至 16 位的 tif）

　　使用 gdalinfo 验证输出数据的格式信息。

　　使用 -outsize 可以控制输出文件的尺寸。

　　使用 gdalinfo 验证尺寸信息。

　　使用 -scale 可以控制输入/输出范围和比例。使用 -mm 开关可以显示像元值极限。

　　使用 -srcwin 可以通过定位参数（xoff yoff xsize ysize）将影象切分成两部分。使用 -projwin 可定义四角地理坐标（ulx uly lrx lry）。

　　可以建立 shp 文件显示栅格切片外框，作为数据索引。对于每幅栅格影象，将生成一个多边形显示其边界，并包含其路径。

　　可以使用 QGIS 和 ogrinfo（相见后述）查看输出的 shp 文件。

　　这里假设 HYP_50M_SR_W.tif 已被恰当地的创建，具备适当的边界。有前述可见，原始影象没有地理坐标信息，故这里假设使用的是 WGS84 地理坐标。

　　使用 gdalwarp 进行投影变换。这里将影象重采样到莫卡脱投影：

　　使用 gdalinfo 验证或显示查看。

　　这里再将影象重采样到正交投影：

　　请注意，地球两极被部分裁切了，这是因为非常靠近两极的边缘是难以重投影的，gdalwarp 放弃了这部分的数据。通过强制 gdalwarp 读取两极数据可以改善这一问题。相见栅格处理教程中的有关内容：http://trac.osgeo.org/gdal/wiki/UserDocs/RasterProcTutorial 。

　　gdal_merge.py 是用于拼接影象的脚本。这里将 east.tif 和 west.tif 拼合成一副影象：

　　拼接也可通过 gdalwarp 实现。它功能更强，但处理速度较慢：

Tip

可以先用 QGIS 等浏览矢量数据。

　　使用 -so 获取具体文件的信息：

　　若不使用以上参数，输出为综述及各个数据集的信息。

　　将结果递交给 grep 可过滤结果，例如根据 COUNTRY 字段：

　　数据可转换至其它格式。格式信息开关为 –formats 。

　　ogr2ogr 用于转换矢量数据的文件。–formats 开关可用于显示各个格式的读写支持。

　　将 countries 文件转换至 GML：

　　你可以尝试完成以下操作：

　　完成了以上最初的尝试后，以下资源将帮助你进一步学习 GDAL/OGR：

　GDAL 官方网站

　OGR 官方网站：

　GDAL 教程：