地理信息系统空间度量算法

空间查询
简单的点线面相互关系拓扑查询包括： 简单的点线面相互关系拓扑查询包括： 面面查询:如与某个多边形相邻的多边形有哪些； 面面查询 如与某个多边形相邻的多边形有哪些； 如与某个多边形相邻的多边形有哪些 面线查询:如某个多边形内包含哪些线； 面线查询 如某个多边形内包含哪些线； 如某个多边形内包含哪些线 面点查询:如某个多边形内有哪些点状地物； 面点查询 如某个多边形内有哪些点状地物； 如某个多边形内有哪些点状地物 线面查询:如某条线经过的多边形有哪些； 线面查询 如某条线经过的多边形有哪些； 如某条线经过的多边形有哪些 线线查询:如与某条河流相连的支流有哪些； 线线查询 如与某条河流相连的支流有哪些； 如与某条河流相连的支流有哪些 线点查询:如某条道路上有哪些桥梁， 线点查询 如某条道路上有哪些桥梁，某条输电线上有哪些变电 如某条道路上有哪些桥梁 站；点面查询:如某个点落在那个多边形内； 点面查询 如某个点落在那个多边形内； 如某个点落在那个多边形内 点线查询:如某个结点由哪些线相交而成； 点线查询 如某个结点由哪些线相交而成； 如某个结点由哪些线相交而成
空间信息量算
2．面积 ． 多边形的面积是一个重要指标。 多边形边界可以分 多边形的面积是一个重要指标 。 解为上下两半， 解为上下两半 ， 其面积就是上半边界下的积分值与下 半边界下的积分值之差。 半边界下的积分值之差 。 设面状物体的轮廓边界由一 …， )表示 个点的序列P 个点的序列P1 (x1 , y1), P2 ( x2, y2 ), …，Pn (xn, yn)表示 其面积为: ，其面积为
v0 v2
v 三维三角形示例
x0 y0 1 (x1− x0) (x2 − x0) 2s = = x1 y1 1 ( y1− y0) ( y2 − y0) x2 y2 1 = (x1− x0)( y2 − y0) − (x2 − x0)( y1− y0)
v0
v1 二维三角形示例
若v0、v1、v2是逆时针排列，面积是正数；反之为负。 v0、v1、v2是逆时针排列，面积是正数；反之为负。
空间查询
1、基于空间关系查询 、 空间实体间存在多种空间关系，包括拓扑、距离、 空间实体间存在多种空间关系，包括拓扑、距离、方 位等。如查找满足下列条件的城市： 位等。如查找满足下列条件的城市： 在京沪线的东部；距离京沪线不超过 公里 公里; 在京沪线的东部；距离京沪线不超过50公里 城市人口大于100万; 城市区域面积 万 城市区域面积5000平方公里 平方公里. 城市人口大于 平方公里
空间信息量算
几何量算
1．长度 ． 线状物体的长度是最基本的形态参数之一， 线状物体的长度是最基本的形态参数之一， 在矢量 数据格式下，线由点组成， 数据格式下，线由点组成，线状物体表示为一个坐标串 (Xi, Yi)，而线长度可由两点间直线距离相加得到。则线 ，而线长度可由两点间直线距离相加得到。 状物体长度的计算公式为： 状物体长度的计算公式为：
Ai = ( yi+1 + yi )(xi − xi+1 ) / 2 2 2 Xi Ai = (xi+1 + xi+1xi + xi )( yi+1 − yi ) / 6 Yi Ai = ( yi2+1 + yi+1 yi + yi2 )(xi − xi+1 ) / 6
按梯形计算重心位置
空间信息查询与量算
查询和定位空间对象， 查询和定位空间对象，并对空间对象进行量 算是GIS的基本功能之一，它是GIS进行高层次 的基本功能之一，它是 算是 的基本功能之一 进行高层次 分析的基础。 分析的基础。在GIS中，为进行高层次分析，往 中 为进行高层次分析， 往需要查询定位空间对象，并用一些简单的量测 往需要查询定位空间对象， 值对地理分布或现象进行描述，如长度、面积、 值对地理分布或现象进行描述，如长度、面积、 距离等。实际上， 距离等。实际上，空间分析首先始于空间查询和 量算，它是空间分析的定量基础。 量算，它是空间分析的定量基础。
空间信息量算
2．面积 ． 三角形的面积公式有： 三角形的面积公式有：
1 s = absin θ 2
1 s = bh 2
b2 s= 2(cot θ + cot ϑ)
s = s(s − a)(s − b)(s − c)
1 s= 4a2b2 − (a2 + b2 + c2 )2 4
1 s = (a + b + c) 2
S1
X
S1=(x2-x1)(y1+y2)/2+ (x3-x2)(y2+y3)/2 + (x4-x3)(y3+y4)/2 + (x5-x4)(y4+y5)/2
空间信息量算
3．弯曲度 ． 弯曲度是描述线状物体弯曲程度的一个重 要参数， 要参数 ， 它定义为曲线长度与曲线的两个端 点之间长度的比值， 点之间长度的比值，即 :
1 xi S= ∑ 2 i=1 xi+1
n
yi yi+1
Y
2．面积 ．
1 xi S= ∑ 2 i=1 xi+1
n
yi yi+1
o
S2
S S1
X
S=S2-S1
Y
(X5,y5) (X1,y1) (X2,y2) (X4,y4) (X3,y3)
1 n xi S= ∑ 2 i=1 xi+1
yi yi+1
o
w=
观测 的路 径长 度 起 点到 终点的 直线 距离
空间信息量算
4. 重心量算
分两种情况： 分两种情况： 1）面的重心。可以理解为多边形内的平衡点， 块被悬挂起来的平衡点。 块被悬挂起来的平衡点。
空间信息量算
面状目标重心可以通过计算梯形重 心的平均值而得到。 心的平均值而得到。将多边形的各个 顶点投影到x 轴上， 顶点投影到 x 轴上 ， 就得到一系列梯 形（如图），所有梯形重心的联合就 如图） 确定了整个多边形的重心。 确定了整个多边形的重心。
空间分析的主要方法
( 1) 基于地图的空间图形分析， GIS中的缓冲区 叠加分析、 基于地图的空间图形分析，如GIS中的缓冲区、叠加分析、数 中的缓冲区、
字高程模型,数字地面模型等； 字高程模型,数字地面模型等； ( 2) 空间动力学分析，有城市扩张模型(驱动力等) 空间动力学分析，有城市扩张模型(驱动力等)、空间价格竞争
空间查询
图形和属性的互查是最常用的查询，主要有两类： 图形和属性的互查是最常用的查询，主要有两类： 1、按属性信息的要求来查询定位空间位置，称为“属性 、按属性信息的要求来查询定位空间位置，称为“ 查图形” 如在中国行政区划图上查询人口大于4000万且 查图形”。如在中国行政区划图上查询人口大于 万且 城市人口大于1000万的省有哪些 称为 万的省有哪些?称为 查询. 城市人口大于 万的省有哪些 称为SQL查询 查询 2、根据对象的空间位置查询有关的属性信息，称为“图 、根据对象的空间位置查询有关的属性信息，称为“ 形查属性” 如一般的GIS软件都提供一个“INFO”工具 软件都提供一个“ 形查属性”。如一般的 软件都提供一个 工具 让用户利用鼠标，用点选、画线、矩形、 ，让用户利用鼠标，用点选、画线、矩形、圆、不规则多 边形等工具选中地物，并显示所查询对象的属性列表， 边形等工具选中地物，并显示所查询对象的属性列表，可 进行有关统计分析。 进行有关统计分析。
i i
其中， 为离散目标物 为离散目标物， 为该目标物权重。 为其坐标。 其中，i为离散目标物，Wi为该目标物权重。Xi与Yi为其坐标。
空间信息量算
5. 形状量算 当把城市作为单个面状目标看待时， 当把城市作为单个面状目标看待时 ， 可以直接使用面状目标的形状系数， 可以直接使用面状目标的形状系数 ， 如形 状率、圆形率、紧凑度等, 状率 、圆形率 、紧凑度等,这些指标计算较 简单，但只反映一个抽象的形状； 简单，但只反映一个抽象的形状； 当把城市作为面状目标的集合看待时 可以使用放射状指数、 ， 可以使用放射状指数 、 标准面积指数等 形状系数， 这些指标计算较复杂， 形状系数 ， 这些指标计算较复杂 ， 但反映 了城市内部的具体联系。 在多数指标中， 了城市内部的具体联系 。 在多数指标中 ， 都以圆形作为城市的标准形状。 都以圆形作为城市的标准形状。
模型(区位优势) 空间择位模型（中心地等） 模型(区位优势)、空间择位模型（中心地等）； ( 3) 空间信息分析，是指根据数据或统计方法建立的模型， 空间信息分析，是指根据数据或统计方法建立的模型，如空间
聚类、空间自相关、回归模型等。实际上， 聚类、空间自相关、回归模型等。实际上，同一种空间分析方法 和模型可以归属于不同的类型。例如， 和模型可以归属于不同的类型。例如，中心地属于空间信息分析 模型，但中心地的形成又追求服务距离最短的动力学机制， 模型，但中心地的形成又追求服务距离最短的动力学机制，可归 属于空间动力学模型。 属于空间动力学模型。
空间信息量算
2）面状分布离散目标的重心 ）
可理解为其分布中心。 可理解为其分布中心 。 其重心计算方法是取离散 目标的加权平均中心， 目标的加权平均中心，它是离散目标保持均匀分布的 平衡点。 计算公式为： 平衡点。 计算公式为：
XG
∑W X = ∑W
i i i i
i
,
YG
∑W Y = ∑W
i i i
空间信息量算
1）形状比(FORM RATIO) 形状比(FORM 形状比= / 形状比=A/L2 其中，A为区域面积，L为区域最长轴的长度。 为区域面积， 为区域最长轴的长度 为区域最长轴的长度。 其中， 为区域面积 该指标能反映城市的带状特征， 该指标能反映城市的带状特征，城市的带状特 征越明显则形状比越小。显然， 征越明显则形状比越小。显然，如果城市为狭 长带状分布，其长轴两端的联系是不便捷的。 长带状分布，其长轴两端的联系是不便捷的。