python 读取geotiff_python与地理空间分析（一）

一. 概述

在气象数据分析中，地理空间要素是一个必须考虑的关键特征项，也是重要的影响因素。例如气温会随着海拔的升高而降低，地形的坡向朝向也会影响风速的分布，此外，典型的地形会形成特定的气候条件，也是数据挖掘中可以利用的区域划分标准。数据分析中，地理空间分析往往能提供有效的信息，辅助进行决策。随着航空遥感行业的发展，积累的卫星数据也成为了数据挖掘的重要数据来源。
地理空间分析有好多软件可以支持，包括Arcgis，QGIS等软件平台，本系列文章将会着重分享python在地理空间分析的应用。主要包括地理空间数据的介绍，常用的python包，对矢量数据的处理，对栅格数据的处理，以及常用的算法和示例。
地理空间数据包括几十种文件格式和数据库结构，而且还在不断更新和迭代，无法一一列举。本文将讨论一些常用的地理空间数据，对地理空间分析的对象做一个大概的了解。
地理空间数据最重要的组成部分：

•地理位置 标示数据空间位置的信息，例如经纬度•主题信息 主题信息涵盖的范围广泛，可以是卫星的光谱数据，也可以是对应地理位置的统计数据如人口，面积，长度等•此外，空间索引一般也是地理空间数据的共性，空间索引创建了一个向导，能够让软件无需扫描数据集中的每一行记录快速定位查询结果。

二. 矢量数据格式

矢量数据目前是最常用的地理空间格式，因为它是存储空间信息最有效的一种方式。矢量数据只存储几何图元，包括点、线和多边形。可供选择的矢量格式数目很饿的，开源矢量库OGR支持超过86种矢量格式，FME支持超过188种矢量格式。下面将简单介绍击中常用的数据格式：

Shapefile文件

使用最普遍的地理空间格式是Esri的Shapefile文件，可以通过arcgis软件或者第三方转换工具包将任意格式的矢量数据转换成Shapefile文件格式，例如OGR库，python中解析Shapefile文件的模块大都是基于OGR库的。 Shapefile文件格式是有多个文件组成，最重要的文件包括3种，.shp、.shx、.dbf文件。•.shp 用于存储要素几何的主文件，其中包含几何图形(读取的对象)•.shx 形状索引文件，可以加快访问速度•.dbf 数据库文件，包括几何要素的属性信息•其他 .prj 以WKT格式存储的地图投影信息 Shapefile文件是一个整体，重命名文件和编辑文件时，需要相关文件保持一致，否则很容易导致数据错误无效。

CAD文件

CAD是工程制图常用的文件格式，一般是通过AutoCAD软件创建的，格式后缀一般为DXF或DWG虽然，在工程制图中也会应用到一些地理空间分析的算法，但是一般的地理空间分析工具对其支持非常少，在实际应用中也很少见。

标签和标记格式

当用arcgis或者QGIS等软件打开数据的时候，会产生相应xml结尾的文件&#x