(e)提供数据转换方面的信息。 ' b5 N2 l/ B/ i! s- A) j m/ p
2、空间元数据的分类 ( z- F8 V4 m: x( H
1)高层元数据(数据集系列Metadata),描述整个数据集的元数据,包括数据集区域采样原则,数据库的有效期,数据的时间跨度、分辨率以及方法等。是用户用于概括性查询数据集的主要内容。
3 k5 I O; r; n! S+ r7 R: h0 ^- [* v2)中层元数据(数据集Metadata),既可以作为数据集系列Metadata的组成部分,也可以作为后面数据集属性以及要素等内容的父Metadata数据集系列。全面反映数据集的内容 。
4 u, { [; Q9 k1 H# v! Q; {5 e- m3)底层元数据(要素、属性的类型和实例Metadata),包括最近更新日期,位置纲量,存在问题标识(如数据的丢失原因),数据处理过程等。是元数据体系中详细描述现实世界的重要部分。
% Y( p* ^4 G; d6 I+ R1 R, I( g3、空间元数据的内容 " n( f5 I7 u. \1 \, x: Q. E
对空间元数据所要描述的一般内容进行层次化和范式化,指定出可供参考与遵循的空间元数据标准的内容框架。 9 u1 L' s! U" \" M9 o; h3 A
第一层是目录层,主要用于对数据集信息进行宏观描述,适合在数字地球的国家级空间信息交换中心或区域以及全球范围内管理和查询空间信息时使用。
2 l. f/ x) |8 l4 {& @( t第二层是空间元数据标准的主体,由八个基本内容部分和四个引用部分组成。 ( R4 K0 ?0 i+ L" A. N. y" u
4、元数据的获取
: |( p: n7 H# n+ F7 t, ?+ {, m数据收集前,得到的是根据要建设的数据库的内容而设计的元数据,包括数据类型、数据覆盖范围、使用仪器说明、数据变量表示、数据收集方法、数据时间、数据潜在利用等。 - g, X1 T& [/ s0 [; [' l) o
数据收集中,随数据的形成同步产生的元数据,例如在测量海洋要素数据时,测点的水平和垂直位置、深度、温度等是同时得到的。
( I; U; D( T0 y3 P# f数据收集后,根据需要产生的,包括数据处理过程描述、数据的利用情况、数据质量评估、数据集大小、数据存放路径等 。 6 B5 O7 {( C5 F, N& o$ z
六、空间数据的互操作 + a7 L% Y) ?$ F' I8 f
1、互操作含义
& {8 G1 r6 E" T, ^% w( D 指异构环境下两个或两个以上的实体,尽管它们实现的语言、执行的环境和基于的模型不同,但它们可以互相通信和协作,以完成某一特定任务,这些实体包括程序、对象、系统运行环境等。 " }/ m$ _$ n9 S
互操作地理信息处理,是指数字系统的这些能力: / k' s3 T# y9 c) ~$ f, E
1)自由地交换所有关于地球的信息,即所有关于地表上的、空中的、地球表面以下的对象的信息。
2 [- {* i, S7 S7 e* U2)通过网络协作运行能够操作这些信息的软件。 ! A6 {# ]. f2 d. c. n
概括为自由交换地理空间信息及协作运行空间信息处理的软件。 / L+ O# H4 T5 F& y
2、GIS互操作类型
" i. E2 \+ e2 T+ U. X1)软件的互操作,强调软件功能块间的相互调用;
8 V Y0 e& h* O1 L9 \; t- ~+ a2)数据的互操作,强调数据集之间相互透明的访问;
5 o- u4 L4 g4 S" }" G0 V% P% U9 e3)语义湖操作,强调信息的共享,在一定语义约束下(对地理现象共同的理解下)的互操作。 ; J: y1 F; q* O8 n
3、GIS 互操作问题
7 b( `$ k# I8 {' K/ v% v2 {0 e5 c目前,所建立的GIS均被认为是信息孤岛,不同系统之间存在互操作问题,因为: 2 R/ I) U8 J5 m( q6 @
1)没有统一的标准,各自采用不同的数据格式、数据存储和数据处理方法; # o2 l2 }0 b2 z( c$ l
2)系统的开发均建立在具体、相互独立和封闭的平台,且不同应用部门对地理现象有不同的理解,导致对地理信息有不同的定义,使得不同应用系统之间在共同协作时无法进行信息交流和数据共享。 ( J. K- H, P. C2 r6 n: G$ D
4、急需实现异构GIS间的互操作 ; w8 q9 M& V2 w: b* t
1)解决基础数据的共享问题的需要; V6 j3 V: y6 S$ q
2)GIS应用趋向多学科综合和集成化;
7 G e3 r' E4 Y8 B3)GIS走向社会化的需要;
. i! ~2 E7 p2 y7 [( V4)也是Internet GIS发展的需要。
0 C3 ^# H) {0 C) z- l9 q! ~" }0 U5、GIS互操作现状
) c; g5 x! g* p7 I( j. P目前,主要有两种方法初步实现互操作: & n/ p0 Q' I; B
1)OPEN GIS规范,通过规定统一的系统设计和开发软件工具的框架,OGC(Open GIS Consortium)OPEN GIS 协会为实现GIS间的互操作制定了OPEN GIS规范。 , }% n* ]8 J$ X- }
2)构件(组件)技术,构件(组件)技术也是实现互操作的可行方法。程序设计中的组件技术,可以在许多不同平台下使用,受之启发,可将GIS某功能包装成独立的组件,使之可以在不同的系统环境下调用。这样可实现系统功能的相互调用。 : {& I4 ^0 S7 u! D
七、Open GIS规范
) `+ m/ V) }9 _1 D* K) B: h( V0 n1、 含义:
* f% b7 l4 b$ s+ Y/ @! h& h4 n OGIS,也叫开放式地理数据交换规程,它是由开放地理信息系统协会(Open GIS Consortium)制定的一系列开放标准和接口。Open GIS规范是OGC规范的最高层次,是利用软件统一地表示地理数据和地理处理的规范系统。
& j( w# j2 P) p1 k5 ~1 E# H1 E2、目的: ( E3 _5 M) Q Q0 N; }
在传统GIS软件与高带宽的异构地学处理环境中架起一座桥梁,具体通过信息基础设施,把地理空间数据和地理处理资源集成到主流的计算机技术中,促使可互操作的商业地理信息处理软件的广泛应用。 % q7 e8 E+ ?( U
3、特点:
/ @( d0 O3 @1 a2 J1)是一种统一的规范,使用户和开发者能进行互操作;
; O) c' m( E- q1 p2)能克服烦琐的批处理及导入、导出障碍,在分布操作系统异构数据库环境下获取数据及数据处理功能资源; - \, `2 s, H3 D( T5 i/ z# I
3)由于Open GIS独立于具体平台,它只能是抽象层的概念描述,而不是具体的实现。
. o& \- H$ J( U8 T1 V) x& M" P. v4、 多数据格式是多源空间数据集成的瓶颈,是OpenGIS出现的基础 ' @+ o3 L% e; C8 U* ]! ]( R: m
1) 多语义性
* Y: P1 Z7 d" g 由于地理系统的研究对象的多种类特点决定了地理信息的多语义性。一个GIS研究的决不会是一个孤立的地理语义,但不同系统解决问题的侧重点也有所不同,因而会存在语义分异问题。
+ ?3 U: m7 z& @0 }/ n2 a2) 多时空性和多尺度
% P$ r# Y" q- Q8 {一个GIS系统中的数据源既有同一时间不同空间的数据系列;也有同一空间不同时间序列的数据。还会根据系统需要而采用不同尺度对地理空间进行表达,不同的观察尺度具有不同的比例尺和不同的精度。
- @1 x" n) a! _9 ]: r5 M& d. u/ s3) 获取手段多源性
1 H- X" D$ e2 }, x Y- D: [% C 获取地理空间的数据的方法有多种多样,包括来自现有系统、图表、遥感手段、GPS手段、统计调查、实地勘测等。 . r) P9 r S+ @4 w0 b9 C
4) 存储格式多源性 ! U" U# |; S( V# c7 y: y
图形数据又可以分为栅格格式和矢量格式两类。传统的GIS一般将属性数据放在关系数据库中,而将图形数据存放在专门的图形文件中。不同的GIS软件采取不同的文件存储格式。
$ R7 n3 }/ d5 F# @6 q5、OpenGIS规范的作用
3 \$ |& z! D8 ]3 a% W+ f$ |通过OpenGIS规范把商业部门、集成部门、用户、研究人员、数据提供商等连接到一起,通过必要的软件工具和通信技术,为各种用户提供对地理信息的共享和互操作。
% h7 _* Y/ A" B8 U& S6、 互操作地理信息的工作方式(OGIS框架) $ ~8 I3 v) b w1 ` e. g/ F
如何实现OpenGIS规范,OpenGIS规范并没有提出具体的标准实施模式,其框架主要由三部分组成
) f* |4 G' P+ a' `1 \9 v+ d0 L( L1)开放的地理数据模型(Open Geodata Model,OGM) ; l0 f" l& m$ j2 R" t0 y+ O& L
包含认可的类型和结构集合(将地理现实抽象为实体(特征)和现象(层)),通过这一集合,可表示任何地理模型。
, |; T u% y2 Q( V' Q! c2)OGIS服务模型(Open Service Model,OSM)
. L: r* i% K/ k6 [ O! w 定义地学数据服务的对象模型,由一组相互可操作的软件构件集组成,为对特征的访问提供对象管理、获取、操作、交换等服务设施。 1 c5 h) q9 X" U/ h
3)信息群模型(Information Communities Model) 4 V' f( ^6 y% W5 j- d- Z
信息群指共享数据的用户群,可以是数据提供者、使用者。不同用户对数据理解不同,引起语义上交流障碍。
; g4 b# H4 i" X 信息群模型,主要任务是解决具有统一的OGM(开放地理数据模型)及语义描述机制的一个信息部门内部以及不同OGM及语义描述的信息部门之间的数据共享问题。采用的主要方法是语义转换,使具有不同特征类定义以及语义模式的信息用户群之间实现语义的互操作。
2 @6 b0 E$ j ]# D. t& Y§4-6 空间数据标准目录& W4 }& o: o+ @, t4 C9 K* d
第五章 空间数据的处理 6 k7 O- a+ d+ e4 v
§5-1 坐标变换 6 G! a2 h8 f4 x+ D8 w7 |3 G' m
一、图幅数据的坐标变换
k1 c$ x# R! Z1、比例尺变换:乘系数
8 z% o, K' o, c2、变形误差改正:
3 G" M' L4 e7 I& a# P# [ 通过控制点利用高次变换、二次变换和仿射变换加以改正 & t. r: c( ~2 }. t6 ]8 a
3、坐标旋转和平移 - }# k9 e+ b- |' q3 e R
即数字化坐标变换,利用仿射变换改正。
5 N2 @6 j" F3 }$ P* u& l' f4、投影变换:
5 d! g) K( i( S& `# V三种方法。
& x1 m' }. L. m+ Z7 o( q二、几何纠正 1 N: N' A% w' |9 _
1、高次变换 ' K6 N. q' G5 [6 q% R
其中A、B代表二次以上高次项之和。上式是高次曲线方程,符合上式的变换称为高次变换。式中有12个未知数,所以在进行高次变换时,需要有6对以上控制点的坐标和理论值,才能求出待定系数。
8 {' A. q) x0 f# a) y2、二次变换 & Z, K' _$ a# @ N+ w# Y' z9 j
当不考虑高次变换方程中的A和B时,则变成二次曲线方程,称为二次变换。二次变换适用于原图有非线性变形的情况,至少需要5对控制点的坐标及其理论值,才能解算待定系数。
& T z% H, k3 V; p+ Z8 y2 ^" f3、仿射变换
! M0 a" v5 }; q0 N实质是两坐标系间的旋转变换。 ) b ?1 ?. f2 ]# ~$ s. X
设图纸变形引起x,y两个方向比例尺不同,当x,y比例尺相同时,为相似变换。 8 X9 h$ B, W: r" f, b4 h
特性:?
- Y$ E" Y' B: [) v- n•直线变换后仍为直线;? + N% Y4 f9 q( @
•平行线变换后仍为平行线;?
" l) [& M1 o3 ~2 @•不同方向上的长度比发生变化。?
' G! c! _* y m/ B 求解上式中的6个未知数,需不在一直线上的3对已知控制点,由于误差,需多余观测,所以,用于图幅定向至少需要四对控制点。
4 F8 C! w6 p6 h' H三、地图投影变换 ) `% e# V( b. S! P
假定原图点的坐标为x,y(称为旧坐标),新图点的坐标为X,Y(称为新坐标),则由旧坐标变换为新坐标的基本方程式为:
& L/ j: Q5 C4 B; A2 r1、解析变换法? * F( u5 J+ W( m4 S
1)反解变换法(又称间接变换法) 5 K! G ^) G6 t: Z$ H/ n
2)正解变换法(又称直接变换法)
2 z" R" I) o' `: k, s p: g2、数值变换法 6 q& `: [0 A4 U# {( }# [1 Y! t/ W4 g
利用若干同名数字化点(对同一点在两种投影中均已知其坐标的点),采用插值法、有限差分法或多项式逼近的方法,即用数值变换法来建立两投影间的变换关系式。
" w/ q9 V1 h4 n2 f5 B3、数值解析变换法
, a) ^6 ^9 V$ x2 ]3 Z7 G当已知新投影的公式,但不知原投影的公式时,可先通过数值变换求出原投影点的地理坐标φ,λ,然后代入新投影公式中,求出新投影点的坐标。即:? & D8 r' a0 [3 b5 }$ Y7 X# j; _
§5-2 图形编辑
- v. n0 _* q8 {; F& s7 Y& M图形编辑又叫数据编辑、数字化编辑,是指对地图资料数字化后的数据进行编辑加工,其主要的目的是在改正数据差错的同时,相应地改正数字化资料的图形。 . c3 V3 d- b5 y6 f W9 E
图形编辑是一交互处理过程, GIS具备的图形编辑功能的要求是: * N7 R2 A7 f% d1 ~
1)具有友好的人机界面,即操作灵活、易于理解、响应迅速等; 6 J. l3 a' Y$ w4 f
2)具有对几何数据和属性编码的修改功能,如点、线、面的增加、删除、修改等;
% X; w' O: I3 v& R' Y- y! A( a3)具有分层显示和窗口操作功能,便于用户的使用。
4 f: p. m7 w+ t6 R一、编辑操作
/ `9 }0 Q3 m1 {" ?' E" P$ Z. X1、结点的编辑
7 H3 B- H2 l9 p, g5 C7 }1)结点吻合(Snap)
+ N. |( n) o# _/ N0 `. i8 ?9 _ 或称结点匹配、结点咬合,结点附和。 ; }3 g+ z1 _9 N- D' Q+ M. b
方法: % k V! j# [) Y% Z- g
A、 结点移动,用鼠标将其它两点移到另一点;
$ C) U6 R( d4 r) JB、 鼠标拉框,用鼠标拉一个矩形,落入该矩形内的结点坐标通过求它们的中间坐标匹配成一致; + t9 W! f5 o: U
C、 求交点,求两条线的交点或其延长线的交点,作为吻合的结点;
9 \6 L ~) X- R5 S; _: ND、自动匹配,给定一个吻合容差,或称为咬合距,在图形数字化时或之后,将容差范围内的结点自动吻合成一点。
+ [$ R* B3 M/ H9 P* g一般,若结点容差设置合理,大多数结点能够吻合在一起,但有些情况还需要使用前三种方法进行人工编辑。
: e' v1 R# _, a: E2 Z$ A2)结点与线的吻合 ; l9 P, i5 K. U0 G7 C
在数字化过程中,常遇到一个结点与一个线状目标的中间相交。由于测量或数字化误差,它不可能完全交于线目标上,需要进行编辑,称为结点与线的吻合。
2 j& P1 C m# f- M编辑的方法:
# z* P4 Z6 w! y% h& ]4 Y7 FA、 结点移动,将结点移动到线目标上。 ( p( N& U5 C& F/ x
B、 使用线段求交;
' X& O8 q: B" [+ ~# p HC、 自动编辑,在给定容差内,自动求交并吻合在一起。
# D; d0 `( n! Q6 R( }9 w3)需要考虑两种情况 6 `2 {* G9 g- t* q$ I( w2 d
A、 要求坐标一致,而不建立拓扑关系;如 高架桥(不需打断,直接移动)
! k; f# Q: R: [ {' qB、 不仅坐标一致,且要建立之间的空间关联关系;如 道路交*口(需要打断)
7 U3 n U$ z9 T1 R3 q/ ?$ {4)清除假结点(伪结点)
3 T: a, U& @5 ]3 D: b由仅有两个线目标相关联的结点成为假结点。
! w) ~' a: M; S1 ?8 D0 \有些系统要将这种假结点清除掉(如ARC/INFO),即将目标A 和B合并成一条,使它们之间不存在结点; $ p7 _0 J+ N ^. V, F# y
但有些系统并不要求清除假结点,如Geostar,因为它们并不影响空间查询、分析和制图。
" o) W8 W% ~1 N3 ?" x; t2、图形编辑
0 B. Y5 ?& v4 _' L* Y& r! f包括用鼠标增加或删除一个点、线、面实体,移动、旋转一个点、线、面实体。
* m \: k% o# [' d1 a- K1)删除和增加一个顶点 ; k s" u2 e# O& `3 p8 d3 n6 ]
删除顶点,在数据库中不用整体删除与目标有关的数据,只是在原来存储的位置重写一次坐标,拓扑关系不变。
6 x3 a1 z6 L& ?; P 增加顶点,则操作和处理都要复杂。不能在原来的存储位置上重写,需要给一个新的目标标识号,在新位置上重写,而将原来的目标删除,此时需要做一系列处理,调整空间拓扑关系。
* d y6 W; L7 t7 B& o2 t! A- @" ?2)移动一个顶点
7 P. L" t' A: B! p r 移动顶点只涉及某个点的坐标,不涉及拓扑关系的维护,较简单。 ' w& P- ]) b2 B/ H7 R
3)删除一段弧段 : C0 F- D3 D" C1 J0 ^3 U
复杂,先要把原来的弧段打断,存储上原来的弧段实际被删除,拓扑关系需要调整和变化. . D2 t! G: n, r) q5 J/ o2 i/ V! a
3、数据检查与清理 * j) T5 R' _, J* Z
数据检查指拓扑关系的检查,结点是否匹配,是否存在悬挂弧段,多边形是否封闭,是否有假结点。
& w' o8 t, W, ]5 u& J9 B! Y* A/ e9 J要求系统能将有错误或不正确的拓扑关系的点、线和面用不同的颜色和符号表示出来,以便于人工检查和修改。
- D; [$ {; d: s/ w7 p0 M& c7 V数据清理则是用自动的方法清除空间数据的错误.
1 M' l9 k$ f" j% \; |( d1 r 例如给定一个结点吻合的容差使该容差范围内的结点自动吻合在一起,并建立拓扑关系。给定悬挂弧段容差,将小于该容差的短弧自动删除。在Arc/info中用Data Clean 命令,在Geostar中选择整体结点匹配菜单。 + r- s# t8 a; k2 ^3 l+ M* h W
4、撤消与恢复编辑 ; ^0 a: p0 H5 M, u8 o
Undo,Redo功能是必要的。但功能的实现是困难的。当撤消编辑,即恢复目标,要恢复目标的标识和坐标、拓扑关系。这一处理过程相当复杂. 2 h( j+ W2 A4 ^/ w% c
因此,有些GIS不在图形编辑时实时建立和维护拓扑关系,如Arc/Info等,而在图形编辑之后,发Clean 或Build命令重新建立拓扑关系。这样,在每次进行任何一次编辑,都要重新Clean 或Build,对用户不便。 2 R9 o3 E; p% x: h- L8 l
二、关键算法 ( C$ t c+ ]7 Y
1、点的捕捉 - o; I% n: v1 N
设光标点为S(x,y), 1 h6 r, F- ~7 _6 ]- Y# H' X$ A6 o
某一点状要素的坐标为A(X,Y)
x) E1 d4 K" i2 v- f可设一捕捉半径D(通常为3~5个象素,这主要由屏幕的分辩率和屏幕的尺寸决定)。 & l; [$ z1 h; h8 C8 ]4 W# m/ V+ {5 P
若S和A的距离d小于D则认为捕捉成功,即认为找到的点是A,否则失败,继续搜索其它点。
2 B$ W( _) b6 T0 \$ E乘方运算影响了搜索的速度,因此,把距离d的计算改为: . f0 `" O- d+ h/ k6 i
捕捉范围由圆改为矩形,这可大大加快搜索速度。 9 c9 L. s" b/ ?, g7 h" a. X' \4 E
2、线的捕捉 # {/ m' t& a3 J
设光标点坐标为S(x,y),D为捕捉半径,线的坐标为(x1,y1),(x2,y2),…(xn,yn)。通过计算S到该线的每个直线段的距离d。.
9 c3 z8 O$ r: ^: `; P5 [1 t6 A若min(d1,d2,…dn-1)<D,则认为光标S捕捉到了该条线,否则为未捕捉到。
" ]' A3 x4 R+ U* _$ X加快线捕捉的速度的方法:
/ n3 U$ U- Q2 X, L8 C+ A1)在实际的捕捉中,可每计算一个距离di就进行一次比较,若di<D,则捕捉成功,不需再进行下面直线段到点S的距离计算了。
5 [( T4 ^. `2 N6 w2)把不可能被光标捕捉到的线,用简单算法去除。 0 L, ^- |5 ?$ F; n' e3 G
3)对于线段也采用类似的方法处理。
" v6 {7 ^7 F- t' w# G/ f+ b8 C& [% D4 K4)简化距离公式: 4 G- z0 \4 g# k( m
点S(x,y)到直线段(x1,y1),(x2,y2)的距离d的计算公式为: 8 c/ U. w n& F$ Q& N5 G3 I
& T' B* i( p7 n; [) U
3、面的捕捉
. F- k% p" r) X [2 N/ o实际上就是判断光标点S(x,y)是否在多边形内,若在多边形内则说明捕捉到。 3 l' Q% U0 j! T6 b1 T
判断点是否在多边形内的算法主要有垂线法或转角法。 ( F! p; W2 \( Z% p* O+ Y0 t4 _
垂线法的基本思想是从光标点引垂线(实际上可以是任意方向的射线),计算与多边形的交点个数。 ; f& S# I4 e) a
若交点个数为奇数则说明该点在多边形内;若交点个数为偶数,则该点在多边形外。 9 }( \( ` @( B, N5 J
加快速度的方法:
; U8 }7 T$ L: J8 C1)找出该多边形的外接矩形,若光标点落在该矩形中,才有可能捕捉到该面,否则放弃对该多边形的进一步计算和判断。
+ o) X9 Y: w6 t2)对不可能有交点的线段应通过简单的坐标比较迅速去除。
$ w/ M4 F4 A1 m7 O1 y; d3)运用计算交点的技巧。
n, ~& `5 p t) ~, o* Q6 i8 `" y4、图形编辑的数据组织—空间索引 : f0 [; s! Y9 `+ N8 g( F
为加速检索,需要分层建索引,主要方法有格网索引和四*树索引。 6 r8 y# m- x7 a" _
1)格网索引 s2 {7 w# [$ i
a、每个要素在一个或多个网格中 0 _' L: p2 k7 t( l/ D2 q
b、每个网格可含多个要素 , k' t$ y6 V; J9 h2 n
c、要素不真正被网格分割 , " ?, Q- |* q1 n8 `
2)四*树索引
# q" n4 J3 f H( Y线性四*树和层次四*树都可以用来进行空间索引。
* l0 ^+ [/ B, f# |9 BA、线性四*树,先采用Morton或Peano码,再根据空间对象覆盖的范围进行四*树分割。 ! ?- e( `+ ]$ _& T8 p% [
B、层次四*树,需要记录中间结点和父结点与子结点之间的指针,若某个地物覆盖了哪个中间结点,还要记录该空间对象的标识。 ) Q1 u% s3 y' w$ S
§5-3 拓扑关系的自动建立
& Y$ ~7 O, A, C6 x& ?, B- a- _0 @1、在图形采集和编辑中实时建立
- r* |6 c& x; Y: _(见龚建雅的书)
# X) [% G9 k3 F% k, N; N' j二、多边形拓扑关系自动建立
/ p+ F$ @4 `' j1、链的组织
- l2 i+ N3 p, A; n6 x# w( Z) s1)找出在链的中间相交的情况,自动切成新链;
. j8 `$ x0 o% v1 x* F7 g2)把链按一定顺序存储,并把链按顺序编号。
P- S" c8 P' i# C1 D; @2 y* m2、结点匹配
& K4 Q: V! Z& n, M1) 把一定限差内的链的端点作为一个结点,其坐标值取多个端点的平均值。
+ B6 f6 P( r9 T( {, a2)对结点顺序编号。
, J1 V. F0 J( N+ k$ X R3、检查多边形是否闭合
4 R& z. M7 ~& \: [5 @$ l通过判断一条链的端点是否有与之匹配的端点来进行.
* D9 H4 D! d- z多边形不闭合的原因:
0 ], ^) H3 j* c+ V+ Z( V1)由于结点匹配限差的问题,造成应匹配的端点未匹配;
6 O/ U& ?# w4 d) ]0 Z2)由于数字化误差较大,或数字化错误,这些可以通过图形编辑或重新确定匹配限差来确定。 4 D0 e; ^/ z1 }$ c; ~
3)还可能这条链本身就是悬挂链,不需参加多边形拓扑,这种情况下可以作一标记,使之不参加下一阶段拓扑建立多边形的工作。 ! c4 E2 b9 b9 a2 S
4、建立多边形 , ^8 k* Y+ g- u, s% P3 }
1)概念
! O! K$ w6 w+ p- A% e1 V9 E. B4 ta、顺时针方向构多边形:指多边形是在链的右侧。
+ J5 _2 Z0 a! ~" n& w- mb、最*右边的链:指从链的一个端点出发,在这条链的方向上最右边的第一条链,实质上它也是左边最近链。a的最右边的链为d
8 ]) v4 r' t1 {' g! Rc、多边形面积的计算 ' ^- {8 G2 ]8 U! o3 F4 D; {4 E4 t
2)建立多边形的基本过程 ' h( ~' o2 |+ k
1° 顺序取一个结点为起始结点,取完为止;取过该结点的任一条链作为起始链。
8 A. ^$ e: u9 _* A2° 取这条链的另一结点,找这个结点上,*这条链最右边的链,作为下一条链。 : N8 B' p" q- h3 ]; f) _
3° 是否回到起点:是,已形成一多边形,记录之,并转4°;否,转2°。
& M; K+ V* Y/ o/ P! Z4°取起始点上开始的,刚才所形成多边形的最后一条边作为新的起始链,转2°;若这条链已用过两次,即已成为两个多边形的边,则转1°。
. O$ S7 C7 n4 v/ P; _5、岛的判断 0 ]0 X+ c8 I9 F/ b! s, k" P
找出多边形互相包含的情况. 0 r1 Z6 S# I& ~( [* F
1°、计算所有多边形的面积。
\7 K. R, q1 P, W2°、分别对面积为正的多边形和面积为负的多边形排序。 3 n: ?& G, D4 L4 K: q
3°、从面积为正的多边形中,顺序取每个多边形,取完为止。若负面积多边形个数为0,则结束。 2 R4 _4 B2 o1 g+ I0 ?! `% J
4°、找出该多边形所包含的所有面积为负的多边形,并把这些面积为负的多边形加入到包含它们的多边形中,转3°。
D9 G5 I, s$ c, ~4 J8 n正面积多边形包含的负面积多边形是关键. 7 ]) R# c% x" J
1°、找出所有比该正面积多边形面积小的负面积多边形。 ) ]9 F% N0 c+ F$ I8 D. \/ Z7 C
2°、用外接矩形法去掉不可能包含的多边形。即负面积多边形的外接矩形不和该正面积多边形的外接矩形相交或被包含时,则不可能为该正面积多边形包含。
5 W5 A' L: n8 {7 @# n) y3°、取负面积多边形上的一点,看是否在正面积多边形内,若在内,则被包含;若在外,则不被包含。 4 o5 |; t' E; y, r( U b4 M- R! v
6、确定多边形的属性
0 v; L3 ~. K4 A1 i& Q" O多边形以内点标识。内点与多边形匹配后,内点的属性常赋于多边形. % Y! f- J- d, t2 D2 g% \
) Z6 ^: i% b5 T' B7 P3 u, b
§5-4 图形的裁剪、合并和图幅接边
% s4 j% A: \. f) [9 b一、图形的裁剪--开窗处理
% v; u' S1 A7 R0 g \- Y- w1、方式:
; u0 D' V& g6 C3 N! T, B" }/ W 正窗:提取窗口内的数据。 开负窗:提取窗口外的数据子集。
7 }5 I$ P9 a' a _1 u/ p8 x( U 矩形窗和多边形窗。
) [7 z6 ~$ X. c$ x# h# N& W2、算法:
/ E l. c4 d' Y6 g1 }( {* J0 i 包括点、线、面的窗口裁剪---计算机图形学。
1 R7 H+ E* ~2 a+ d4 s! s 而不规则多边形开窗------相当于多边形叠置处理。
7 c* X9 {1 }2 d# C+ S7 K/ _二、图形合并---数据文件合并 . t9 X7 a/ Q( t( o7 k; f+ r' S
一幅图内的多层数据合并在一起;
% A& {8 R- }" c9 i( \ 或将相邻的多幅图的同一层数据合并.
4 U1 ]4 t y9 g/ R P 涉及到空间拓扑关系的重建。对于多边形,由于同一个目标在两幅图内已形成独立的多边形,合并时,需去除公共边界,属性合并,具体算法,删去共同线段。 实际处理过程是先删除两个多边形,解除空间关系后,删除公共边,再重建拓扑。 * e5 w5 d# y7 j5 P9 @1 h; f
三、图幅接边—形成无缝数据库 ( \9 y" P- R' n6 }; _
几何裂缝:指由数据文件边界分开的一个地物的两部分不能精确地衔接。--几何接边
$ G. Z# e& J! V4 h逻辑裂缝:同一地物地物编码不同或具有不同的属性信息,如公路的宽度,等高线高程等。---逻辑接边 % w. K- d4 p( s; h* G! W% t
1、识别或提取相邻图幅。--要求图幅编号合理
. V5 l J }: ~& E2、几何接边 ! z' D+ ?4 O4 a o1 s
3、逻辑接边
7 i$ T! N4 G% p5 g: \9 \1)检查同一地物在相邻图幅的地物编码和属性值是否一致,不一致,进行人工编辑。
' R- @* ^2 P% z- y0 B- E2)将同一地物在相邻图幅的空间数据在逻辑上连在一起。
6 i; v; z3 u" J: J9 ]/ v§5-5 空间插值
) H' W! i! t p5 |3 T空间插值:
% ` k1 i. Y1 V. s4 U( Q6 i$ k内插:在已观测点的区域内估算未观测点的数据的过程;
. F& t! [+ s3 ~0 j外推:在已观测点的区域外估算未观测点的数据的过程.—--预测。
3 Y( @4 N8 P9 }/ X" y8 _一、边界内插
, c6 z- A/ | i6 w 首先假定任何重要的变化都发生在区域的边界上,边界内的变化则是均匀的、同质的。
3 B" i: ?, h1 G! A% e( s 边界内插的方法之一是泰森多边形法。 % Q6 P' Z8 p: i5 l
泰森多边形法的基本原理是,未知点的最佳值由最邻近的观测值产生。
, H( r6 o. o8 M6 \0 B4 K B二、趋势面分析
& w# w0 K# _2 _4 h: f是一种多项式回归分析技术。多项式回归的基本思想是用多项式表示线或面,按最小二乘法原理对数据点进行拟合,拟合时假定数据点的空间坐标X、Y为独立变量,而表示特征值的Z坐标为因变量。 $ J5 |# V4 t: T# W2 _" S3 Y
1、当数据为一维时,
# h0 g# ?4 M5 M6 ~1)线性回归: % T; _+ x/ F/ W# H
2、数据是二维的 4 F* m1 N e9 g" C5 q( s
二元二次或高次多项式 5 Z# Y5 C' X: L; e
三、局部内插 8 }9 A7 J. L* l
利用局部范围内的已知采样点的数据内插出未知点的数据。 ( _& f; j& g+ ^8 L
1、线性内插
1 N7 S0 H; u9 ^7 Y; @: g将内插点周围的3个数据点的数据值带入多项式,即可解算出系数a0、a1、a2 。 ! K3 }4 ]- [$ Z) w7 E
2、双线性多项式内插
% o, P. {8 |& n! B; n将内插点周围的4个数据点的数据值带入多项式,即可解算出系数a0、a1、a2、a3 。
1 k/ @7 m" v$ ?$ n& T' x3、双三次多项式(样条函数)内插 : }- M4 P, R4 g3 \; h
是一种分段函数,每次只用少量的数据点,故内插速度很快;样条函数通过所有的数据点,故可用于精确的内插;可用于平滑处理。
3 I/ S K1 c& P0 x1 e0 N, i 双三次多项式内插的多项式函数为: 6 C& s: P! Y# n# t/ G1 ^1 V2 B9 W
3 E Q+ C g+ v
四、移动平均法
6 Z% e6 d, z& q" v+ P8 C* l在局部范围(或称窗口)内计算个数据点的平均值
8 m C# J$ J# F* k' ^. Q: R二维平面的移动平均法也可用相同的公式,但位置Xi应被坐标矢量Xi代替。
7 k( A5 y2 s& ~" P O 窗口的大小对内插的结果有决定性的影响。 4 ]6 [6 R9 @: @! ?
小窗口将增强近距离数据的影响; 大窗口将增强远距离数据的影响,减小近距离数据的影响。
+ {+ t& `( \3 D' {当观测点的相互位置越近,其数据的相似性越强;当观测点的相互位置越远,其数据的相似性越低。
- I7 q# ?9 H, E6 \加权移动平均法:λi是采样点i对应的权值
. S4 F- q# ]" L& y/ {. d, E8 _" D! P1 i加权平均内插的结果随使用的函数及其参数、采样点的分布、窗口的大小等的不同而变化。通 , a B0 A6 Z2 {0 u) F3 O! L
常使用的采样点数为6—8点。对于不规则分布的采样点需要不断地改变窗口的大小、形状和方向,以获取一定数量的采样点。
: M: h! f) h6 \' R§5-6 数据压缩与光滑 1 V, S u- z; q c
一、数据压缩
+ h V6 T$ ^8 a+ S' Y7 X矢量数据压缩
2 s5 R, u. B7 L# N/ m" A4 ]栅格数据压缩
5 Z; c3 B8 c$ A- K( M3 Z1、 Douglas—Peucker
6 ? a' e5 ]$ o* X5 F5 O压缩效果好,但必须在对整条曲线数字化完成后才能进行,且计算量较大;
# x( k6 L, w5 n' Y+ F6 q2、垂距法
5 h, `( k8 r/ ~$ x6 H每次顺序取曲线上的三个点,计算中间点与其它两点连线的垂线距离d,并与限差D比较。若d<D,则中间点去掉;若d≥D,则中间点保留。然后顺序取下三个点继续处理,直到这条线结束。
! {9 c: F. E* d; S7 y压缩算法好,可在数字化时实时处理,每次判断下一个数字化的点,且计算量较小; " G6 T8 P3 R. e# T( W
3、光栏法 , ]7 r% v* H. i
定义一个扇形区域,通过判断曲线上的点在扇形外还是在扇形内,确定保留还是舍去。 6 g0 A6 v1 w3 i T# \- v8 c% [
算法简单,速度快,但有时会将曲线的弯曲极值点p值去掉而失真。
( Q) \5 _. i, b6 e, V# |6 P; g二、曲线光滑(拟合) % N: m! C: }+ M% \' w6 A
是假象曲线为一组离散点,寻找形式较简单、性能良好的曲线解析式。 4 [9 X0 T8 {* F T9 I
插值方式:曲线通过给定的离散点。如拉格朗日插值,三次样条曲线 ' C8 n* i, }3 V5 [) r0 A+ h3 |* K
逼近方式:曲线尽量逼近给定离散点。如贝塞尔和B样条曲线。
( u: U4 I& g2 Q5 i9 R3 s) t一、矢量向栅格转换 . C1 ~" R! _) t" Z0 @+ C+ A% O
点:简单的坐标变换
, m) s% C: O6 @+ {; y3 b5 |线:线的栅格化
. F4 h" e$ M/ X2 j* [- \+ l G面:线的栅格化 +面填充 / w6 ]3 a* R; X& ]
(一)线的栅格化
$ q* { H4 p$ A& I1、DDA法(数字微分分析法)? 5 ~7 S; W# e# q; I% ?7 C
2、Bresenham算法 . h( W3 ^. q! p0 A& a
(二)面(多边形)的填充方法
; Z5 a4 [3 |9 E1 N1、内部点扩散法(种子扩散法)? * H3 L" I$ `3 |3 w
2、扫描法?
3 E; R# A; g4 c2 ?; G% Y& c3、边填充算法 # o; C. x2 i0 f5 a
二、栅格向矢量转换 $ O8 @' |6 y. F2 K
从栅格单元转换为几何图形的过程为矢量化; 4 |: c" g! O* p" W% e3 S/ S u% x' v# ]
(一)要求(矢量化过程应保持):
, h4 U2 F ]* b& `% w; Y1) 栅->矢转换为拓扑转换,即保持实体原有的连通性、邻接性等;
# |) `! Q' \4 ^- O2) 转换实体保持正确的外形。
4 h* e4 }4 |8 b) G) H0 G, x* A(二)方法
6 J- C4 S- X# Z* e0 |方法一,实际应用中大多数采用人工矢量化法,如扫描矢量化,该法工作量大,成为GIS数据输入、更新的瓶颈问题之一。 ( j7 e! C m/ H/ v+ i" V
方法二,程序转化转换(全自动或半自动)
- ~7 U6 D# l/ X8 o过程为:
5 r) d0 v$ n2 Q& Y0 }1、边界提取 4 k- e6 P* P; y
2、二值化 # y( Y( _7 C8 Y4 z" A( Q
3、二值图像的预处理 ( u6 T$ z0 q* I) _* J! P" U! E
4、细化:1)剥皮法 2)骨架法 ) _# ^ G3 Y% C) E: F
5、跟踪 6、拓扑化