Remote SensingGeographical Information System
: b' T5 d; o: ^" g(RS & GIS )
# P/ @3 H) J% w& `遥感及地理信息系统* d; h( p; P4 [2 u3 y. {% m
(讲稿)
1 ^' e4 c! l- F! M 第一章GIS 概述
1 K7 j. f% U& [2 Z' a. F- F, X/ G§1.1 地理信息系统的产生与发展
$ P( p- a3 e& o& O( b6 A3 i" Q一、GIS产生的背景:' O4 _4 Q1 S9 w
地理学——地图绘制(手工、计算机)——世界上第一个地理信息系统:5 J5 T- P+ h4 l4 z" l' J
CGIS(1963年)
+ G2 S2 R8 L' x; q3 S4 i二、GIS的发展5 S; F6 |& z) C3 e! [
1. 计算机技术的发展
! B5 D8 O/ w7 A: h7 F9 v9 }2. 空间技术的发展& W+ ~6 b* ~5 }8 D- e' \
卫星定位技术
5 y0 m3 q. A' t. e$ u卫星定位的特点 1.全球地面无缝覆盖,全天候定位。2.定位精度高,实时定位速度快。3.观测时间短。4.提供三维坐标。5.使用简便。' f2 T1 F( E0 z
GPS系统:美国于1978-1993年发射,共24颗卫星。其中21颗工作卫星,3颗备用卫星。GPS卫星向用户发送用于导航定位的调制波,它含有:载波(L1和L2)、测距码(C/A码和P码)和数据码(导航电文)。C/A码伪距,精度约为20米左右(民用)P码伪距,精度约为2米左右(军用)
+ k. I% J$ g9 S0 ZGLONASS系统:苏联于20世纪70年代开始建设,并于1993年启用由24颗卫星组成。目前在轨14颗卫星;预计于2010年左右使卫星数量达到“满员”状态。民用信号定位精度:30米。
" Q2 X# ?% U- D5 ?& t! Y7 S“伽利略”系统:30颗卫星组成欧空局和欧盟合作于1999年启动。2005年12月28日,首颗试验卫星Glove-A成功发射;预计2008年底前全部发射入轨。民用信号精度:1米。+ L1 Q8 T2 m5 K% F8 C5 A+ w5 ~
“北斗”系统:7 P8 e3 J# B( J7 [0 d1 U
“北斗1代”:2000-2003年发射3颗静止轨道试验导航卫星,组成了“北斗”区域导航系统(2颗工作卫星、1颗备用卫星)。2 i# l+ Y' K) n6 X. B1 X
北斗2代:“北斗”全球卫星导航系统。空间段由5颗静止轨道卫星和30颗非静止轨道卫星组成。定位精度为10米。
- \+ K" ~8 V& p5 p5 N5 o5 z卫星遥感技术: o7 f7 g. ^7 T+ x Y% W
卫星遥感的特点:能够快速、大范围获取地面影像和地表专题数据( H2 L: y- h3 G
?SRTM计划——航天飞机雷达地形测绘任务(2000年)
0 p% M: v; Y& ^ T" r3 g(Shuttle Radar Topography Mission)4 B" V U( Y! s" r: ]
覆盖范围在北纬60度至南纬56度之间。10天内,有99.968%被一次覆盖,94.59%两次,49.25%三次,24.1%四次。
9 _" t& v0 e+ V& F+ XDEM高程数据间隔为1弧秒,约30米,所包含的信息内容相当于1∶5万地形图
2 f) c6 a2 R* U2 I& rSRTM由三部分组成:主雷达天线、桅杆、机外雷达天线(Outboard radar antenna)
5 E5 f: a9 Q! _+ g x, B) M8 p* v现代GIS的相关技术:8 K: p& E& r( \; p
地理学、遥感技术、测量学、数学和统计学、制图技术、计算机科学、专家系统、计算机图形学、计算机辅助设计、数据库技术、软件工程7 Y7 u0 R" i, u0 V4 f( S1 s
§1.2 GIS的相关概念
! v" f! V' L) F! d9 v% X6 U* b/ a: p2 G1、数据:是客观事物的属性、数量、位置及其相互关系等的抽象表示,如数字、文字、符号、4 |) S0 ?1 E1 j [( A
图形等。
" {+ Z$ W0 v6 f2、信息:是向人们或机器提供关于现实世界各种事实的知识,是数据、消息中所包含的意义。它不随载体的物理形式的变化而改变。/ F7 S( q: Q% ~- I& @% Q! Z$ i
数据与信息的关系:3 W; E' j8 ]3 L( Z0 {
数据是信息的载体?信息是数据的内涵?数据本身并没有意义,数据只有对实体行为产生影响时才能成为信息?数据只有通过解释、解算才能成为信息--使用者的知识背景。
! p$ T; c& b0 l4 {' `1 c9 Z4 X信息的特点:0 ~& Z' @* u; s6 p8 b* @7 Z7 P6 R
(1)信息的客观性:任何信息都是与客观事物紧密相关的-信息正确性与精确度的保证。
% B1 g' I' T3 q3 N: b7 O& y(2)信息的适用性:不同的信息适用于与之相关的领域。
2 g4 @' ]# N* m; A$ C1 y(3)信息的传输性:信息在传输时,其原始意义并不改变。+ M0 F4 K$ y: ?* c9 k& G
(4)信息的共享性:信息的使用对其本身无损失。 W0 d& J k2 b- z
3、系统:具有特定功能的、相互有机联系的若干要素构成的一个整体。
0 W3 s* U; P* o4、信息系统:具有处理、管理和分析数据功能的系统,它能为人们的决策提供正确、有用的信息。" k8 |8 S$ ]/ G+ j0 e5 `2 L" c
5、地图:是将地理环境诸要素按照一定的数学法则,运用符号系统并经过制图综合缩绘于平面上的图形。a 由数学决定的结构b 特定的符号系统c 现象表示的取舍和概括
: v3 |1 b9 S* H+ U ?: Q. v/ R/ t. S6、地理信息:是指表示地理环境诸要素的数量、质量、分布特征及其相互联系和变化规律的数字、文字、图象和图形等的总称。0 K9 w; U1 A) b8 }
地理信息的特征:
# E* j- i- A, w P5 G, K(1)定位特征:通过地理坐标实现空间位置的识别) v- \2 X' Q. q+ Q- J" k1 M
(2)多维结构的特征:除地理位置外,还有多个专题和属性4 u" @5 V2 B6 u# w$ N7 i2 v
(3)时序特征:随时间的变化而变化,需及时采集和更新?
2 o1 I/ j6 K5 K, ]3 K- V! F* k地图与地理信息的关系:
) l c% i/ [6 D9 O, Y$ K; V地图是地理信息的载体?地图是地理信息的传统数据源?地图是GIS的查询与分析结果的表示方法' V; D/ L' W/ B
§1.3 地理信息系统的概念
$ P2 n* |6 w( J j7 w0 e地理信息系统是一门学科:是研究采集、管理、分析地理数据和输出地理信息的理论和方法的新兴、交叉学科。
6 _: O* i) H$ j" c" T0 z地理信息系统是一种决策支持系统:在计算机软、硬件系统的支持下,运用系统工程和信息科学的理论,科学管理和综合分析具有空间内涵的地理数据,以提供管理、决策、规划等所需地理信息的技术系统。
( I. E' ]/ F; X7 y; p5 N+ U) nGeographical Information System1 f8 E& t1 U9 C8 v" I
Geographic Information System* _& ~7 ? j- C; x1 B
Geo-Information System
" O" G) n: K; T0 f$ y; o' C地理信息系统中“地理”的概念并非指地理学,而是广义地指地理坐标参照系统中的坐标数据、属性数据以及以此为基础而演绎出来的知识。3 r! j( A9 T/ g
?地理信息系统的分类——按功能分类:
" _; C& ]" g8 y1 A5 f( Z# E1. 专题地理信息系统(Thematic GIS):是具有有限目标和专业特点的地理信息系统。为特2 c; J( q) K2 Y0 \- h2 Y" G
定的专门的目的服务,如水资源管理信息系统、矿产资源信息系统、农作物估产信息系
6 u/ q- A1 S# S( |* }) f/ f7 [. ?- r统、草场资源管理信息系统、水土流失信息系统、环境管理信息系统等。
. f: b- P8 ^& D3 W" Q( |2. 区域地理信息系统(Regional GIS):主要以区域综合研究和全面信息服务为目标。如国
" l! v9 r: C6 @3 T. A1 h 家级、地区级、市级或县级等。1 G& q* O: B' `6 h
3. 地理信息系统工具(GIS Tools):是一组具有图形图像数字化、存储管理、查询检索、" i# I- ?! K6 T9 _
分析运算和输出等地理信息系统基本功能的软件包。( ]' p) E: G6 g! \2 ^# _) a: |
GIS与其它系统的区别:
7 b- l0 b& V L0 e: y2 KGIS与DBMS(数据库管理系统)的区别' |) T1 p% P% [; R
1 . DBMS缺乏空间实体定义能力* [" [4 T6 C7 _) J9 p# r* s; H3 z6 H
2 . DBMS缺乏空间关系查询能力* g; a$ A" @7 k) ~4 A" H6 v- I8 X3 R
3 . GIS是能对空间数据进行分析的DBMS,GIS必须包含DBMS。& [4 A6 F- B9 Z7 p
GIS与MIS(管理信息系统)的区别
1 O2 ^9 P( H ^+ T1. GIS要对图形数据和属性数据库共同管理、分析和应用,GIS的软硬件设备要复杂、系统功能要强;" r5 I( P3 G! M6 L: b4 f
2. MIS则只有属性数据库的管理,即使存贮了图形,也是以文件形式管理,图形要素不能分解、查询、没有拓扑关系。, E ~; a# q, x1 Y1 G
GIS与其它系统的区别
8 P9 m. D" w$ Z" F# dGIS与地图数据库的区别:
# t) ]) H; g6 d地图数据库仅仅是将数字地图有组织地存放起来,不注重分析和查询,不可能去综合图形数据和属性数据进行深层次的空间分析,提供辅助决策的信息,它只是GIS的一个数据源。8 k X/ N& u9 z" F" a
GIS与CAD系统的区别:
" T/ S: X- [; `; a, g1 . 二者虽然都有参考坐标系统,都能描述图形,但CAD系统主要处理几何图形,属性库功能弱,更缺乏分析和判断能力。2. CAD的坐标无地理意义,不能进行地理坐标转换。
5 U+ e- x8 Q; R1 T v& h" O§1.4 GIS的发展历史4 Y: s9 K, x' D+ {* _: W* j
一、国外GIS的发展
9 w- }' A2 o! c+ A9 j1.起步阶段(1960s):注重空间数据的地学处理。相关组织机构纷纷成立。1963年,加拿大测量学家R. T. Tomlinson首先提出GIS这一术语,建立加拿大地理信息系统(CGIS);1969年,ESRI (环境系统研究所)建立;1969年,Integraph公司建立。
* k/ u5 {$ j/ N$ M2.发展阶段(1970s):注重空间地理信息的管理,受到政府部门、商业公司和大学的普遍重视。
- x" g0 f; M2 J3.推广应用阶段(1980s):注重空间决策支持分析。应用领域迅速扩大,与RS结合,参与解决全球性问题。1981年,ESRI ARC/INFO GIS发布;1985年,GPS成为可运行系统;1986年,MapInfo 建立;1986年,SPOT卫星首次发射;1987年,地理信息系统的国际杂志出版;1989年,Ingegraph 发布MGE;1989年全球有报价的GIS软件达70多个。/ Y4 `4 j# X0 _7 E( T5 Y m, u/ b
4.用户时代(1990至今): 注重GIS社会应用与服务-商业决策和政府管理GIS技术迅猛发展-GIS产业化。控件式GIS成为GISTools的发展方向;WebGIS蓬勃发展;三维GIS崭露头角。
- b) i, x2 ]% k( N* E0 q5.国外主流GIS软件:ARC/INFO(ArcView、ArcObject、ArcIMS)-ARCGISGENAMAP(澳大利亚)MGE(Modular GIS Environment)MapInfo(MapinfoProserver、MapX、MapXtreme、SpatialWare)ERDAS
8 R$ C$ `& ]9 C6 } ~二、我国GIS的发展# V" p1 X( ^/ o9 i0 T m. p
1.准备阶段(1970后期): 舆论准备、提出倡议、组建队伍、组织个别实验研究。1980年中科院遥感所成立第一个GIS研究室$ C, A% L+ m+ _7 u
2.试验阶段(1981-1985): 理论探索和区域性研究;制定国家地理信息系统规范。1:100万国土基础信息系统和全国土地信息系统1 :400万全国资源和环境信息系统
- N, y: |6 k+ \. w# k3.初步发展阶段(1986至20世纪末): 地理信息系统的研究被列入我国“七五”攻关课1 ?" h. ^( R/ w& R+ f* m4 H
题,并且作为一个全国性的研究领域,已逐步和国民经济建设相结合,并取得了重要进展和实际应用效益。: e) h$ ]3 m8 U5 E' o8 I l
4.快速发展阶段( 20世纪末到现在): 理论日趋成熟,应用日益广泛,三维GIS、WEBGIS 走向应用,GIS市场开始形成。5 x( Z& m' o6 f1 Y) M7 T) s- W
5.国产主流GIS软件: GeoStar(武测)MapGIS(地大)CityStar(北大)
8 t' Z: S3 E0 B, a K6 L( J% {三、地理信息系统的发展趋势
% I3 X* c- F/ Z4 K$ f(1)GIS与遥感和全球定位系统进一步结合,构成地理学日趋完善的技术体系-3S技术;
& O/ e$ o* `. G& c# T1 x1 ~7 `9 \(2)空间数据结构与数据管理的研究更加深入;* i! O# j& e! V7 f- g
(3)GIS应用模型开发日趋加强;7 z/ U4 f! T" L" ]
(4)GIS智能化; Q7 E/ M- q) ~% T, Y9 y5 E% x
(5)GIS网络化-WEBGIS;" C- a o4 Q5 V: s8 ^+ g7 t
(6)三维GIS7 g& ~! U' Q3 S" g/ R7 J
第二章遥感的基本原理及其应用& \7 t1 m* O( G
§2.1 遥感概述
6 O0 {. N; ]1 Q' Z$ i一、遥感的概念
% O# t3 M% U; u8 v20世纪60年代发展起来的对地观测综合性技术。遥感一词的英文为Remote Sensing,意思是遥远的感知。60年代由美国人EvelynPruitt(美国海军研究局)提出。6 b. e2 n# ?- J9 U3 W& ^
广义理解:遥感泛指一切无接触的远距离探测,包括对电磁场、力场、机械波(声波、地震波)等的探测。实际工作中,重力、磁力、声波、地震波等的探测被划为物探(物理探测)的范畴。因而,只有电磁波探测属于遥感的范畴。
0 y9 v& `7 ~5 R+ K4 ^狭义理解:遥感是指从不同高度的平台(Platform)上,使用各种传感器(Sensor),接收来自地球表层的各种电磁波信息,并对这些信息进行加工处理,从而对不同的地物及其特性进行远距离探测和识别的综合技术。
& \2 o9 G9 H7 N y A7 [; x" r二、遥感技术系统
: o( Q. U1 H+ Q6 N* p8 Y根据电磁波理论,不与目标物接触,从远处用探测仪器接收来自目标物的电磁波信息,通过对信息的处理和分析研究,确定目标物的属性及目标物相互间的关系的综合技术系统,称为遥感技术系统。
- b3 t8 ~4 j& ?' @; m$ q/ _遥感技术系统包括:# v* o" O" u2 z
1. 遥感信息源- z+ P3 h) H7 g
2. 空间信息获取- o9 Y" F+ r" f. w
3. 遥感数据传输与接收
~: q ^* S* [. f9 Q, y# O4. 遥感信息处理、提取与分析
- r, X! D- e# T/ m1 S2 o4 u5. 遥感信息应用 T5 ]3 ]& H# q' m& e) B X
6. 遥感技术系统及工作流程2 K i. D3 j8 Q* p
7. 遥感技术系统组成
- B4 u; ~( K Y4 E) L- J4 ]遥感信息源:任何目标物都具有发射、反射和吸收电磁波的性质,这是遥感的信息源。目标物与电磁波的相互作用.构成了目标物的电磁波特性,它是遥感探测的依据。& B8 b6 }' U. ]. l9 ?5 F. R0 r! y
空间信息获取:地物空间信息主要由搭载在遥感平台上的遥感器来获取。接收、记录目标物电磁波特征的仪器,称为遥感器。
+ F* d, _, S' l3 J/ L W- b: l遥感数据传输与接收:传感器接收到目标地物的电磁波信息,被记录在数字磁介质或胶片上。胶片是由人或回收舱送至地面回收,而数字记录信息则可通过卫星上的微波天线传输给地面卫星接收站。遥感技术系统组成
1 {2 {, V: ^8 F6 I" l6 J遥感信息处理、提取与分析:地面站对接收到的数字信息进行一系列的处理:信息恢复辐射校正卫星姿态校正投影变换通用数据格式转换地面站或用户还可根据需要进行精校正处理和专题信息处理、分类等。遥感技术系统组成2 Q% N2 P: n: N2 Y$ N2 C0 Y$ e
遥感信息应用:遥感获取信息的目的是应用。这项工作由各专业人员按不同的应用目的进行。遥感技术系统组成
4 z; q; @. k6 c5 p7 |9 r6 i, l; A3 h三、遥感的类型 ]& t1 _* ]" s
1.按遥感平台分
, l! {+ |/ k% n* M& k地面遥感:遥感器设置在地面平台上,如车载、船载、固定或活动高架平台等;# {( m4 M. T5 n
航空遥感:遥感器设置于航空器上,主要是飞机、气球、飞艇等;
, q( a2 e- v' u" Q- X8 ?3 F) [" h; ?/ _- H航天遥感:遥感器设置于环地球的航天器上,如人造地球卫星、航天飞机、空间站等;
y0 l* h* {* \0 C9 \2.按传感器的探测波段分
" n+ x3 J& p' J1 z4 f/ j紫外遥感——探测波段在0.05一0.38μm之间;8 I8 T3 T; B- K2 M7 c2 e+ O
可见光遥感—探测波段在0.38一0.76μm之间;
( H4 J- u( ? h ^红外遥感——探测波段在0.76一1000μm之间;* v; O/ I* _, n+ Z) \7 Q( e& x
微波遥感——探测波段在1mm一1m之间;
. E a3 M1 u2 f多波段遥感—指探测波段在可见光波段和红外波段范围内,再分成若干窄波段来探测目标。! G3 w" ~) G, P# N
3.按工作方式分
$ w' ~( b( H5 y0 B# d+ w主动遥感和被动遥感:主动遥感由探测器主动发射一定电磁波能量并接收目标的后向散射值量;被动遥感的传感器不向目标发射电磁波,仅被动接收目标物的自身发射和对自然辐射源的反射能量。
; G* m% z' `4 p9 ]9 _成像遥感与非成像遥感:成像遥感——传感器接收的目标电磁辐射信号可转换成(数字或模拟)图像。非成像遥感——传感器接收的目标电磁辐射信号不能形成图像。
# l2 e. } r2 _; Z$ m1 T- V- b4.按遥感的应用领域分8 |; p- m2 s- Z7 @
从大的研究领域:可分为外层空间遥感、大气层遥感、陆地遥感、海洋遥感等;
( U! |3 d4 F# C. d. o- U从具体应用领域:可分为资源遥感、环境遥感、农业遥感、林业遥感、渔业遥感、地质遥感、气象遥感、水文遥感、城市遥感、工程遥感及灾害遥感、军事遥感等,还可以划分为更细的研究对象进行各种专题应用。
3 M$ k4 L3 z: _- l3 P四、遥感的特点0 o8 W% Y3 C* Z) A
1.大面积同步观测:遥感观测不受地形阻隔等限制,探测到的地面范围大,容易发现地球上一些重要目标物空间分布的宏观规律。
. \9 P4 h' t9 c2.时效性强:遥感探测可以在短时间内对向一地区进行重复探测,发现地球上许多事物的动态变化。这对于研究地球上不同周期的动态变化非常重要
1 @5 w9 N* W# T0 m2 ], k3.数据的综合性与可比性好:遥感获得的地物电磁波特性数据综合地反映了地球上的自然、人文信息。如地球资源卫星LANDSAT所获得的地物电磁波特性均可以较综合地反映地质、地貌、土壤、植被、水文等特征。7 w3 L2 M2 Z! K- P* R
4.较高的经济与社会效益:与传统的方法相比,可以大大地节省人力、物力、财力和时间,具有很高的经济效益和社会效益。
6 h9 j* Y4 o) \+ B( {* _: F5.一定的局限性:所利用的电磁波波段有限、已经被利用的电磁波谱段对许多地物的某些特征还不能难确反映、地面调查和验证尚不可缺少。5 L& z0 e0 E x) Q
五、遥感的发展简史& ^! S2 d& G. \
1. 无纪录的地面遥感阶段(1608-1838年): 1608年,汉斯·李波尔赛制造了世界第一架望远镜,1609年伽利略制作了放大倍数3倍的科学望远镜,从而为观测远距离目标开辟了先河。6 x3 f6 P% K! Z* l4 g# z6 {) Y
2. 有记录的地面遥感阶段(1839-1857年): 对探测目标的记录与成像始于摄影技术的发明,并与望远镜相结合发展为远距离摄影。1839年,达盖尔(Daguame)发表了他和尼普斯(Niepce)拍摄的照,第一次成功地把拍摄到事物形象地记录在胶片上。6 n* _, [* O6 Z6 ^" U0 P0 m
3. 航空摄影遥感阶段(1858-1956年): 1858年,G·F·陶纳乔用系留气球拍摄了法国巴黎的“鸟瞰”像片。1903年,W·莱特和O.莱特(WilbourWright & OrvilkeWright)发明了飞机。1909年,W·莱特在意大利的森托塞尔上空用飞机进行了空中摄影。1 c* v& H. p6 c9 j
4. 航天遥感阶段(1957-):
! Y. k2 C( e* X: X) X7 b6 s1957年10月4日,苏联第一颗人造地球卫星的发射成功。5 |* Y F/ y8 ^( g" l
1960年,美国发射了TIROS-1和NOAA-1太阳同步卫星,从航天器上对地球进行长期观测。
3 U9 f* o" f+ h9 b7 [- Z0 v8 m. M' f1972年ERTS-1发射(后改名为Landsat-1),装有MSS传感器,分辨率79米。2 d0 p/ Q) M6 u4 W6 g" N% u. q
1982年Landsat-4发射,装有TM传感器,分辨率提高到30米。3 _1 Y; f8 n7 _" W# z; r
1986年法国发射SPOT-1,装有PAN和XS遥感器,分辨率提高到10米。
$ w* e' i6 m* H: y1988年9月7日中国发射的第一颗“风云1号”气象卫星。! r" x6 m9 A' ]% G9 C$ [# u: w8 w, x3 }
1999年4月15日Landsat-7发射,携带了增强型主题成像传感器(ETM+)。
% c' M3 t9 O) R2 U/ R1 @& d 1999年9月24日美国发射IKNOS,空间分辨率提高到1米。
8 z+ R7 b0 s& J; d) O! h1999年10月14日中国成功发射资源卫星1号。
5 b$ |3 J! Q4 Y" _: s( A2001年10月发射的快鸟卫星,可获取0.61米分辨率的黑白影像以及2.44米分辨率的多光谱(4波段)影像。
2 K; h f( ^% @" P, s$ c2002年其5月发射Spot-5,其性能作了重大改进,可获取2.5米分辨率的全色影像。# x( y) K" c: I4 U
六、遥感技术的发展趋势
( N1 ~% L( P! H- R8 H) T) w, H1. 多分辨率传感器并存一方面,高分辨率的图象(如1m分辨率的IKONOS)可以会有较快的发展,这是由于小卫星技术日趋成熟,发射成本越来越低,它在军事侦察与城市测绘方面有独突的作用。另一方面,低空间分辨(如1km 分辨率的NOAA 系列)将会得到较快的发展,因为实践表明这种尺度较适合于全球变化的研究。
* y' A( k5 c( ?) f' \ M9 V) v3 Q2.高光谱遥感遥感器波段宽度窄化,针对性更强,可以突出特定地物反射峰值波长的微小差异;同时,成像光谱仪等的应用,提高了地物光谱分辨力,有利于区别各类物质在不同波段的光谱响应特性。: l' A6 V% c. i% K
3. 遥感分析技术从“定性”向“定量”转变定量遥感成为遥感应用发展的热点。
2 k) ]! X: w$ U. Q7 p' s4. 建立适用于遥感图像自动解译的专家系统逐步实现遥感图像专题信息提取自动化。+ e4 f& Z9 S& [% j$ c/ r& e$ C
5. RS与GPS、GIS 的有机结合构成集成系统
9 Y0 `7 U# t% L0 y§2.2 遥感平台及遥感器
# `- p4 h* v% [ N' I& V一、遥感平台 a5 h% m! E+ z. k
遥感平台:——遥感中搭载遥感器的工具。
9 n& T% N W7 q. Q6 A1.遥感平台的种类
, D2 u. U. M3 @+ Z6 p地面平台——地面遥感+ f6 n3 R7 j) b4 Z- Q! z& y
航空平台——航空遥感
8 C! c, f% e2 A K9 `1 z7 |5 U7 `航天平台——航天遥感3 `# m% s5 E. A# S, \
按遥感平台的高度和运行方式分为:
* F2 A( L6 B' t! n% a; b% U(1)地面平台:高度在0~100m范围内,包括三脚架、车、船、遥感塔等。在平台上放置地物波谱仪、辐射计、分光光度计等,可以测定各类地物的波谱特性。对地观测研究中应用较少。主要目的:对地物进行波谱测量,近距离摄影测量,大气辐射纠正等。
) ^! Q2 ?; ?& }7 H B5 f. V' ~(2)航空平台:飞机高度在百米~10几km,包括低、中、高空飞机,以及飞艇、气球等。航空平台历史悠久,主要目的是航空摄影。低空:n+1 ,甚至可以达到3n 。B、代表性:训练样本应该反映该类地物的亮度特征;且需注意同类地物空间分布的不连续性(如;可以分别在2块草地中选取训练样本)。C、亮度分布:从亮度特征角度而言,对于同一类地物具有不同亮度特征情况,通常都要选取(同物异谱)。D、纯度的把握:一般来说,在已知类别的中间部位较高;如果过纯,则方差过小,导致许多此类地物的像元被拒分;如果不纯,则方差过大,导致不是此类地物的像元被分到该类别中,导致出现“错分”的情况;把握纯度的方法:读像元的亮度值,确定该类地物的亮度分布区间;同时,注意其它地物亮度值与其近似部分的像元亮度分布区间,依“错分损失最小”的原
. u' q) Y4 r7 I- u' i0 c 则,确定“分类界限”。
m$ q( {7 {; g, n& m非监督分类与监督分类的结合:
1 b1 w0 T6 ?; O8 J9 k( m! }通常非监督分类做为监督分类的前导;如果监督分类的最终结果要8类,那么做非监督分类时致少要12类。以此保证分类的纯度,因为方差小,错分的可能性也小。
) O) \/ s) k7 p. ~1 P+ ?3 }§2.7 遥感图像地学分析与应用; w# C4 o# z, `8 D/ N3 C
地学分析基础:
& `' s/ v5 E" {9 V: Z: D. F8 V9 ^遥感影像反映的是某一区域特定地理环境中的综合信息,它综合地反映了地球系统各要素的相关性。
! R9 E3 ^8 k. |4 T/ ~7 t遥感地学相关分析,指的是充分认识地物间的相关性,在遥感图像上寻找目标识别的相关因子,建立直接解译标志和间接解译标志,推断和识别目标本身。
l4 g5 v0 F. z9 w$ Q; J地学分析前的准备:
4 W1 y( v9 X* g ~ F+ Q1 R) i4 M地学(地理、地质等)背景知识
4 ~, L) r8 L' [+ m- Y& N遥感影像的尺度特征 {% i- r8 b$ ^7 z% Q
混合像元# F1 J/ O7 }5 P" e/ _' J
光谱分辨率
! c$ Z8 D8 c- C时间分辨率(获取时间)
. V v2 I- a' J$ E; e一、水体信息提取
& a$ f1 f& [7 j o' A5 C$ X* i+ h, c通过对遥感影像的分析,获得水体的分布、泥沙、有机质等状况和水深、水温等要素的信息,从而对一个地区的水资源和水环境等作出评价。
5 W; u5 Q+ c8 v5 n" R6 F% j1 ]* z ?水体界限的确定(水体的分布)9 X4 b, s" U: g& F3 P
水体悬浮物质的确定$ n/ n8 |/ j( N2 n- h+ {
水温的探测' R, P5 B+ S+ n! ?. t
5 ?+ ?8 }+ O; M+ ?水体污染的探测
V4 J+ l+ {, X# H水深的探测. w6 B+ R0 t# O3 g8 g+ x* |
1.水体的光谱特征 B/ s9 D+ L4 e( v o% g5 |9 M7 i: B
遥感器所接收到的辐射包括水面反射光、悬浮物反
5 v% |5 t* N- V5 S y& B: H射光、水底反射光和天空散射光。由于不同水体的水面性
0 b: Z$ Z n" S6 S质、水体中悬浮物的性质和含量、水深和水底特性等不同,4 d6 @6 b, |* ]* z7 w" Z
从而形成传感器上接收到的反射光谱特征存在差异,为遥" l3 o8 i1 ?' E" M& c
感探测水体提供了基础。
& a8 c& H0 c- `9 P- `2.水体界线的确定# a3 X! d% C9 L! M4 L9 A
在可见光范围内,水体的反射率总体上比较低,不
# l+ y4 G1 S- w9 m N超过10%,一般为4%~5%,并随着波长的增大逐渐降
; l" o! A! Z% w( @; q低,到0.6μm处约2%~3%,过了0.75μm,水体几乎
3 o) O5 v0 _9 w# ?% u" u) _成为全吸收体。因此在近红外影像上,清澈的水体呈黑色,4 Q, F! F" S( w% R. ~" H$ \) ]) F @2 S
即对近红外光全吸收。为区分水陆界线,确定地面上有无水体覆盖,应选择近红外波段的影像。
- q d6 I3 H7 R$ B4 i3.水体悬浮物质的确定
2 e) s& P$ m4 W B) W2 g ?. R遥感能够探测的水中悬浮物质主要有两种:一种是无机的泥沙;一种是有机的叶绿素。
+ `) F& ~( y8 r2 ~) H1) 泥沙的确定
% b! M/ u4 D! w5 D7 o2) 叶绿素的确定; J+ P: E( `# U; ?4 W
? 混浊水体的反射波谱曲线整体高于清水,随着悬浮泥沙浓度的增加,差别加大;, y5 z; Y, ~7 y0 A0 _
? 波长较短的可见光,如蓝光和绿光对水体的穿透能力较强,能反映出水面下一定深度的泥
# M, Q" E) x. t* w 沙分布状况;. P# q; ?/ S$ o' W- i! n! f- I7 F
? 随着悬浮泥沙浓度的加大,可见光对水体的透射能力减弱,反射能力加强。因此,在近岸的浅水区,其波浪和水流对水底泥沙的扰动作用比较强烈,使水体混浊,故遥感影像上色调较浅。而深水处扰动作用较弱,水体较清,遥感影像上水体色调较深。
* U9 C2 a( b& m7 B# Q相对水深的探测
) o" e J+ @- p2 t二、地质地貌信息提取
7 I2 u0 M- z7 C, _; p地质遥感实际上就是从遥感影像上提取地质专题信息的技术过程,地质遥感的应用目的就是为地质制图、矿产资源的探查、工程地质和水文地质调查服务。主要包括:
+ H" k1 g8 o8 S% F, P# Z8 L◆岩性的识别/ Z ^5 Y' d7 N
◆地质构造的识别
$ r( v) N$ l3 R. v1. 岩性的识别:. e* m7 L5 z( C
所谓岩性就是根据岩石的矿物成分、结构构造来判断岩石的种类。不同的岩石,由于其矿物成分、结构构造的不同,具有不同的电磁波辐射特性,在遥感影像上表现为不同的影像特征,岩性的识别就是根据各种岩石在遥感影像上特有的影像特征来判断岩石种类的过程。TM5,TM7为区分岩石性质最好的波段,各种岩石的光谱差异最明显。
" v; J9 e$ {! ^4 S V1 v$ ]; R三大岩类的影像特征:3 y) Y" C1 X# z" U
沉积岩的条带状影像特征——沉积岩最大的特点是具有成层性,具有明显的层理构造。胶结良好的沉积岩,出露充分时,可在较大范围内呈条带状延伸。在高分辨率影像上可以显示出岩层的走向和倾向。
2 X6 r+ Q; J d0 S岩浆岩影像——岩浆岩因其生成的条件与沉积岩差别很大,因此,在岩浆岩发育的地区则常常见到节理(即岩石的破裂面),而基本上看不到层理;, g3 B+ k* A4 S0 \! s
变质岩——特征构造是具有叶理(岩石中片状结构或组织的通称),和沉积岩的层理不同,它是变质作用的形成的。其中的矿物都呈定向排列。
! ?' K, E4 l- r$ B8 `2 a2. 地质构造的识别! p/ d' l" V3 Z( n6 i1 d
构造解译是遥感图像地质解译的重点,也是遥感图像用于地质上效果最好的领域。由于遥感图像视域辽阔,概括性强,用它研究规模较大的构造以及一些隐伏构造和深部构造,效果常常优于地面观察,能将工作区内各种构造形态解译出来。包括岩层产状的判读、褶皱、断层、断裂构造的判断。
2 A, D- I: }0 A地形起伏区水平岩层—封闭曲线、云朵状环带。
+ u0 s2 v* {5 [- J8 k0 \岩层三角面:某一岩层同一层面地表露头线上山脊点与相邻两沟谷点用直线连接所成的假想三角形平面(随地形不同的变态:梯形、熨斗形、半圆形、半月形)。
' q* D/ m4 ~8 J/ Q0 ^岩层三角面——遥感图像上判断和量测岩层产状的最佳标志。" W9 L/ l6 e7 J* i2 D
断层的突出影像特征就是线性,在线形构造的两侧是岩性不同的岩石,通常岩石突然消失的地带就是断层。
2 Y! z* [8 ~; g- s% S9 q 第三章GIS的组成
# M% `/ b) R! o4 C0 o- W§3.1 GIS 的计算机硬件系统0 F' k1 l% R6 N) g2 I j# A
计算机硬件是计算机系统中的实际物理装置的总称,可以是电子的、电的、磁的、机械的、光的元件或装置,是GIS的物理外壳,系统的规模、精度、速度、功能、形式、使用方法甚至软件都与硬件有极大的关系,受硬件指标的支持或制约.GIS由于其任务的复杂性和特殊性,必须由计算机设备支持.GIS硬件配置一般包括四个部分:# K, n) Q* q4 {% G- o; x
1、计算机主机;% \! U/ X' R! k) g4 V0 b( i4 ?
2、数据输入设备:数字化仪、图像扫描仪、手写笔、光笔、键盘、通讯端口等;
3 g& _4 k7 U! s# B# _# R: W$ @3、数据存贮设备:光盘刻录机、磁带机、光盘塔、活动硬盘、磁盘阵列等;
) @/ l4 e+ \; D1 P" ^% ^' L# n; M4、数据输出设备:笔式绘图仪、喷墨绘图仪(打印机)、激光打印机等。: C. C* N/ M0 V5 f9 s1 ?8 |
/ R+ E" p3 y; i6 O* n% V. X §3.2 GIS的计算机软件系统及模块划分) h1 i1 h2 I; i( ^! i1 Z6 F) O
一、软件的种类:
* p6 P* M$ q, V: h1 N* c- `* A7 X计算机系统软件:操作系统、汇编程序、编译程序库、程序使用手册、程序说明。$ H& ~/ O T6 ]5 ~
GIS软件和其他支持软件:通用的GIS软件包、数据库管理系统(DBMS)、图形图像处理系统 应用分析软件:系统开发人员或用户根据地理专题或区域分析模型编制的用于某种特定应用任务的程序,是系统功能的扩充与延伸。
. E4 P; ^$ S2 \7 E二、GIS软件系统的模块划分! ~1 e4 R) O! C# |& ?5 p" N
GIS的软件系统根据其功能的特点,可以分为以下几种软件模块:
5 w4 s3 g& ?* V6 {+ A' T1. 数据输入与编辑模块3 O* R4 f6 ^1 I) Z
2. 数据存储与管理模块6 h* {" F* M; ^! c5 U
3. 数据变换与分析模块" C6 d6 G# k& ~
4. 数据输出与显示模块2 T7 b' M/ a! Q$ J/ N
5. 用户接口模块: K6 r* K- t1 p$ s4 B1 L# n
各功能模块之间的关系
+ c- g2 B( _* B" }+ R( s( C+ M6 s
§3.3 GIS的地理数学基础9 c) |5 V5 F6 c
一、地球的形状及相关概念9 {* K. k4 x% u, J/ {
1.大地水准面:不考虑潮汐、风浪及大气压的影响,与平均海水面相重合,并延伸到大陆下面处处与铅垂线(重力方向)相垂直的水准面,称为大地水准面。
7 ?) L% x( q9 I0 d* S特点:由于地壳质量分布的不均匀,引起重力方向的变化,导致处处和重力方向成正交的大地水准面成为一个不规则的、高低起伏的曲面。该曲面无法用数学公式表达。
% B0 H9 Z! Q) [- l6 [0 [6 j2 C0 F% ?, V; o6 Q0 t; y, K4 U2 O! B
2.大地球体:闭合的大地水准面所包围的形体,叫大地球体。
6 S5 Z; j5 X! v6 m& Y3.地球椭球体:大地水准面形状虽然十分复杂,但从整体来看,起伏是微小的。大地球体是一个很接近于绕自转轴(短轴)旋转的椭球体。在大地测量和制图中用旋转椭球来代替大地球体,这个旋转球体称为地球椭球体。
8 l# ~( F' _6 X, [" p4.参考椭球体:对于国家测图或区域制图,往往采用其大小及定位定向最接近本国或本地区的地球椭球体。这样的地球椭球体称为参考椭球体。
; ^; |3 c* g" o7 S特点:参考椭球体在大小及定位定向上都不与地球椭球体一致。 定位不一致:地心不重合;定向不一致:地轴不重合或不平行 5. 国家高程基准:0 \- J. R' w* m- ]4 V) }
(1)高程基准面:地面点高程或海底点水深的统一起算面,称为高程基准面。大地水准
( {$ b$ q9 s/ Y9 g* e" F+ J面通常作为高程基准面。对于一个国家来说,一般根据一个验潮站的潮汐资料所求得的平均海水面作为国家高程基准面。7 [! A) Q) W7 y
(2)水准原点:传递高程的起算点。用精密水准测量方法,将水准原点与验潮站的水准! ?3 J8 B$ T) H" M# e2 k9 S7 I* ?, W
标尺联测,以高程基准面为零,推求水准原点的高程。
$ X- N/ S: Y9 @1 }(3)我国高程基准面使用青岛验潮站多年潮汐资料推求的平均海水面,作为中国高程基
" Q0 ?5 t1 `1 G1 S准面。
/ ?" O( U- T$ [青岛验潮站:位于青岛港内 我国水准原点:位于青岛观象山 (4)我国高程基准面的变更:
3 T8 {- }1 K4 R?1956年黄海高程系统:1950年-1956年7年间的潮汐资料推求的平均海水面,作为中; Z8 `& {+ r$ F, `' V
国高程基准面。水准原点高程:72.289米?
; r2 ]+ m2 `- q* c: v1985国家高程基准:1952年-1979年(28年)间19年的潮汐资料推求的平均海水面,) h0 w1 g2 M. v6 |+ H( M
作为新的中国高程基准面。水准原点高程:72.260米1988年1月1日启用,同时“1956年黄海高程系统”废止。
- E A4 Q# ~: S(5)我国其他海水水深起算面:
# z" }+ Z2 x3 A7 f?理论深度基准面--海图深度基准面 国家高程基准-理论最低潮面≈2.4m ?短周期平均海平面* b8 `* `3 D7 J8 k
二、地球的大小
5 `) S, f% O2 b! h# ~$ I; O+ s, [?地球椭球体的大小和形状:
% G8 ^$ J8 S+ m8 f6 R& z2 H' v! b9 R常用旋转椭圆的5个基本几何参数表示:长半轴(长半径)a;短半轴(短半径)b;扁率9 O, E6 Y7 c4 X5 M4 ]9 J" Q
α=(a-b)/a;第一偏心率e=(a 2-b 2)0.5/a;第二偏心率e′= (a 2-b 2)0.5
. h* g6 k" B6 u& s. T5 t1 F, E/b' v- L* Z: b: V+ @3 t: E
; O3 S5 V2 `% G?椭球参数根据天文大地测量、重力测量及人造地球卫星观测资料综合确定。4 G/ k* O/ t* A/ \; s
世界各国常用的地球椭球体参数; i( t. K5 }' Q0 U6 O; g
1980国家大地坐标系- v3 b J' a' z+ v
1/298.257 6356755.2886378140.08 U+ @6 b9 k2 Y' w* Y
IUGG (1975) 1954年北京坐标系3 b2 S4 S* R& I
1/298.3* @( I% q @. l. m/ A
63568635 ]: _3 `# D& k9 _
63782458 H# J3 p/ S2 G
克拉索夫斯基 (Krassovsky 1940)" C+ j& c1 h2 L( `: v
1/297 6356912 6378388 海福特(Hayford 1910) 1/293.5 6356515 6378249 克拉克(Clarke 1880) 1/299.15 6356079 6377397 白塞尔(Bessel 1841) 备注
$ ~- u. t. w- v. S扁率 短半轴(米) 长半轴(米) 椭球体名称
! z! @* u1 l" a( }; i/ P 椭球参数是进行坐标转换的基础!
: C4 D \) P4 J# i7 R三、GIS常用坐标系+ p# F% P( r4 R# Q! B% ^. d
1.大地坐标系:又称地理坐标系
4 h8 z; ?9 |3 Z0 a8 `?地理坐标系是以地理极(北极、南极)为极点。. \, D4 F$ M7 i" w
/ `7 m8 D+ Z3 ]* U4 ^0 ^
?通过A点作椭球面的垂线,称之为过A点的法线。2 o, s( c4 Y9 d6 T+ R. _5 p/ l
?法线与赤道面的交角,称为A点的大地纬度,B。
3 ]& C: C3 U: o7 T' Q?过A点的子午面与通过英国格林尼治天文台的子午面所夹
) n: F8 U, j: I; x& |- T的二面角,称为A点的大地经度,L。
9 c ?2 N2 ?+ D6 _& }?A点距大地水准面的高度,称为A点的大地高,H。
% k( t8 Q; {4 i0 a/ h2 [ zA点的地理位置由B、L、H确定. A4 \& {- D5 u* N* S: I
2.空间直角坐标系4 i. ^+ I$ I: }/ j4 l+ E; a3 v
坐标原点:地球椭球(参考椭球)$ U) a. _; E6 z7 w, O
Z轴:与地球自转轴重合,指向北极点。
8 Q* h& F8 Q; i' C" u1 aX轴:指向本初子午线与赤道线的交点! ?5 }* K) \, p
: H3 _5 _. n7 JY轴:垂直于X、Y轴确定的平面,向东为正。
0 Y" @ Y0 f/ O7 Y) @空间任意点的坐标用X、Y、Z表示: a. O5 N! j/ g) y- X9 A# _
3.平面直角坐标系& n) A2 q' c, H. ]0 ~
运用地图投影的方法,建立地球表面和平面上点的函数关
, Z6 [4 }4 N3 Y5 e% G系,使地球表面上任一点由地理坐标(L、B)确定的点,在平面8 F6 H q& @3 h8 O
上必有一个与它相对应的点,平面上任一点的位置可以用直角坐$ U: k$ Y. q4 U, t5 z# g
标表示。+ e& h( H3 }/ W( R: L: u
在平面上选一点O为直角坐标原点,过该点O作相互垂直
0 t3 @9 B- Z* t4 A# ?, _3 _+ H- {/ F的两轴X’OX和Y’OY而建立平面直角坐标系平面直角坐标系。
/ P! s% |% l, }( w( L注意:纵轴为X轴,指向北为正横轴为Y轴,指向东为正
8 q: k) j! u0 l7 Z四、地图投影
" G+ I# D5 S5 `" Y+ R+ b; t为什么要进行投影?8 ]2 L! Z% F9 y' g, t
地图投影实质3 {/ C0 f- u2 `/ [5 L: ?
投影变形/ U! L. t: C( Z/ E( n& r
投影分类" e! y- f6 Y; H5 P" |
投影选择所考虑的因素
9 {1 `6 }1 y0 g我国常用的投影方法
% o, A8 {# Q6 q' s. P7 O) V1.为什么要进行投影?
6 j/ c; k6 e4 ~5 i6 b定义:将地球椭球面上的点映射到平面上的方法,称为地图投影。- j& ]: ~ v" Y6 _& v! `
? 地理坐标为球面坐标,不方便进行距离、方位、面积等参数的量算。2 c! l) Y* n( Z" }4 V& M7 j/ r
? 地球椭球体为不可展曲面。
4 \5 G; M# T- G! X; g3 U% i; z? 地图为平面,符合视觉心理,并易于进行距离、方位、面积等量算和各种空间分析。 2. 地图投影实质1 I7 |2 D) r. n0 q; J1 a
建立地球椭球面上经纬线网和平面上相应经纬线网的数学基础,也就是建立地球椭球面上的点的地理坐标(λ,φ)与平面上对应点的平面直角坐标(x,y)之间的函数关系:当给定不同的具体条件时,将得到不同类型的投影方式。
1 p& d* Z1 V& L5 h地图投影示意! S9 I3 N( J/ x/ D3 w
8 u% S m" ? r$ v( W* B3. 投影变形:
7 F G* _+ `$ M: s1 G# k将不可展的地球椭球面展开成平面,并且不能有断裂,则图形必将在某些地方被拉伸,某些地方被压缩,故投影变形是不可避免的。
k, X1 V" |- Z/ i5 r2 [/ D' T5 z! v长度变形 面积变形 角度变形2 N/ O; W* t* \# a
4. 投影分类
6 [- D8 o7 |7 W% i1 c; A) m! f8 ?$ {- s) l2 N! u
变形分类:5 v; M9 R L$ h- Q
等角投影:投影前后角度不变
- y7 L/ F, m1 W7 X* W5 X% w等面积投影:投影前后面积不变;$ b+ }; N% i! e! ^) L6 [2 o* u1 o
等距投影:投影前后长度不变变形* ?) @) W9 Q( O+ i/ B/ f
投影面:6 ?9 b9 M3 `/ A% i; Y$ a
圆柱投影:投影面为圆柱(包括椭圆柱)
/ V( w$ |/ B! ~$ y+ z圆锥投影:投影面为圆锥/ S) ^3 m; x; Q, K6 V
方位投影:投影面为平面
( X% r/ g3 Z# F; J" j) k# @ K投影面位置:
8 s% Y [8 f2 p/ _1 @" h7 d' @正轴投影:投影面中心轴与地轴相互重合
U9 z- X+ K! _1 I斜轴投影:投影面中心轴与地轴斜向相交
r( n; Q. J# ]# [) [( {; d0 _横轴投影:投影面中心轴与地轴相互垂直
: B Z0 P W& x" @3 d; |相割或相切:1 u: i0 C* T4 K7 X4 o0 {4 n
相切投影:投影面与椭球体相切
% u- g# K9 D# h9 @$ [7 g相割投影:投影面与椭球体相割
* j% _/ q2 T, L7 K$ p5. 影响投影选择的因素7 b* e& U) J' m2 Y
?制图区域的地理位置、形状和范围/ ?9 n/ Z( o3 R. h
?制图比例尺' l* A' o( F- d1 k. Q# K; G
?地图内容
3 S+ w/ b9 G. d?出版方式8 g( d% n) U! g$ d4 q: f' q
地图投影选择,主要指中、小比例尺地图,不包括国家基本比例尺地形图。
5 r( |8 Z% U5 Z1 e国家基本比例尺地形图的投影、分幅等,是由国家测绘主管部门研究制订,不容许任意改变。' ^/ j, l3 L+ i2 ~7 f% ]0 }% j
小区域大比例尺地图,无论采用什么投影,变形都很小。
4 a2 |) r$ ^8 _2 c; B% I& AGIS中,地理数据的显示可根据用户的需要而指定投影方式,但当所显示的地图与国家基本地图系列的比例尺一致时,一般采用国家基本系列地图所用的投影。
% {0 h0 B$ H- j/ h8 m! A: w& i 6. 我国常用的投影方法; w8 `: O- N* C. `) N
高斯-克吕格投影(Gauss-Kruger)——国家基本比例尺地形图(≥1:50万)
+ e- H* K* B' R: C$ r( n1 Q5 K9 n墨卡托投影(Mercator)——海图1 u# S y& A( T3 h
兰勃特投影(Lambert)——国家基本比例尺地形图(1:100万)、大部分省区图
' f0 ~" p7 Z6 o# o% }通用横轴墨卡托投影(UTM)
; [& D) L, J6 o高斯—克吕格投影:2 n- x/ D. b- b& _% T
横轴等角椭圆柱切投影假想一椭圆柱面横套在地球椭球体上,并与某一条子午线(中央子午线)相切,椭圆柱的中心轴通过椭球体中心,将中央子午线两侧各一定经差范围内的地区投影到椭圆柱面上。
/ p' w- c0 V4 R+ D2 y1 U?投影条件:1)中央经线和地球赤道投影成为直线且为投影的对称轴;2)投影后无角度变形;3)中央经线上没有长度变形。
% [; x1 l* S" K* w% ??特点:1)中央经线上无变形;2)同一条纬线上,离中央经线越远,变形越大;3)同一条经线上,纬度越低,变形越大。
2 C' G$ p- T# |3 ]3 X?我国高斯投影的分带:我国规定按6°或3 °分带; f( {9 C! H8 c8 Z( R
6度带:从0度子午线起,自西向东每隔经差6为一投影带,全球分为60带,各带的带号(n)用自然序数1,2,3,…60表示。中央子午线经度L0=6n-36度带3度带;
5 l- O$ y) ^" e n, h3度带:从东经1度30分的经线开始,每隔3度为一带,全球划分为120个投影带,各带的带号(n′)用自然序数1,2,3,…120表示。中央子午线经度L0= 3n′;" K! }4 ?) s, c- R! f
; p% e- I; S- n9 X?我国高斯平面直角坐标系?以中央经线为X轴,赤道为Y轴,两轴的交点为坐标原点。
) H$ F& b4 E8 R7 c) v?为了避免Y坐标出现负值,将各带的坐标纵轴西移500公里,即将所有Y值都加500公里。
5 Q& a3 R G a0 p) ]$ B, e?各带的投影完全相同,某一坐标值(x,y),在每一投影带中均有一个,不能确切表示该点的位置。因此,在Y值前,需冠以带号,称之为通用坐标。$ O0 e& i( x: Y( q4 r
例如大学路操场内某点:X=3995920 Y=21 259470——21带6 w# u4 E C/ l
墨卡托投影:
9 r. Z! P0 y4 t3 D0 v$ m& d等角正轴圆柱(切、割)投影?圆柱与地球椭球的切线或割线称为基准纬线。1 e4 a+ w7 p/ x6 N1 S
特点:1. 没有角度变形-等角航线。2. 经纬线都是平行直线,且相交成直角。3.经线间隔相等,纬线间隔从基准纬线向两极逐渐增大。4.常用作航海图和航空图。
@. H; P& T- ~: _% b8 L?墨卡托投影坐标取零子午线或自定义原点经线与赤道交点的投影为原点,零子午线或自定义原点经线的投影为纵坐标x轴,赤道的投影为横坐标y轴,构成墨卡托平面直角坐标系,此投影标准纬线无变形。
" Q) y& F, L3 T' S: q* ~3 v3 {五、我国常用坐标系统
/ p- B: [) a6 I1 u' {8 t1954年北京坐标系:20世纪50年代初,将我国一等锁与原苏联远东一等锁相连接,然后以连接处呼玛、吉拉宁、东宁基线网扩大边端点的原苏联1942年普尔科沃坐标系的坐标为起算数据,平差我国东北及东部区一等锁,由此传算过来的坐标系定名为1954年北京坐标系: a.属参心大地坐标系; b.采用克拉索夫斯基椭球参数; c. 大地原点在原苏联的普
( a7 ^ x' T& U* { O6 A 尔科沃;
! K% t/ A) G0 }% }: }+ F) U1980年国家大地坐标系
, E. K9 u. U: y. N1 N?又称1980西安坐标系
* }' H' |7 |3 W+ o: c; T?为进行全国天文大地网整体平差而建立(1)属参心坐标系。椭球短轴Z轴平行于地球质心指向地极原点方向,大地起始子午面平行于格林尼治平均天文台子午面;X轴在大地起始子午面内与Z轴垂直指向经度0方向;Y轴与Z、X轴成右手坐标系; (3)采用IUGG 1975椭球参数(3)大地原点在我国中部,具体地点是陕西省径阳县永乐镇;
( O" m/ ?$ g* k) SWGS一84坐标系
- o# U% z0 P/ { K2 V* G?美国国防部研制确定的大地坐标系
$ o( k3 U+ T8 _% T6 k- g* y. w?GPS采用的地心坐标系。
% K. F* e0 O7 O3 t?原点在地球质心,z轴指向BIH 1984.0定义的协议地球极(CTP)方向,X轴指向BIH 1984.0 的零子午面和CTP赤道的交点。Y轴与Z、X轴构成右手坐标系,又称为1984年世界大地坐标系统。?椭球参数采用IUGG第17届大会的推荐值。( U1 u$ D( y) W# Q) g4 i) H
六、坐标系统与投影小结
( t' ? N9 V, U+ V7 z3 H: U坐标系统是由坐标原点位置、坐标轴的指向和椭球参数所定义的。; X8 @% K, S4 j: k8 C+ M. e
任一坐标系统均可分为大地坐标系和空间直角坐标系。而其平面直角坐标系需根据投影方式确定。, O& {+ @, U3 @
由于坐标系统定义所采用的参数不同,地球上的同一点,在不同坐标系统中的大地坐标值不同。$ \6 {% x% K, o; C/ _
在同一个坐标系统里,坐标的转换是严密的,而不同坐标系统之间的坐标转换是不严密的。
* X. p" G, p: A2 {* V8 Y2 d§3.4地理空间数据
# Z v1 | [$ }2 e. n0 w$ {/ k: F8 |一、地理空间数据的定义
, t- l- Q ~8 e' [. t4 E; n9 H6 ^; W观测地球某一位置的地物所得到的、以空间位置为参照的数据称为地理空间数据。
, _) ]* `5 E S7 r地理空间数据描述的主体称为地理实体,任何一种地理实体都可以用以下三种基本实体来抽象描述:1)点实体2)线实体3)面实体。
# S( f/ B4 f: I2 a点实体——?有位置,无宽度和长度;?抽象的点?用坐标对表示(x , y)(L , B)) R- _4 s$ H% H$ I4 u4 O" \: L
线实体——?有位置、长度,但无宽度和高度?用来描述线状实体,通常在网络分析中使用较多?度量实体距离?用坐标对序列表示(x1, y1),(x2, y2)…. (xn, yn)/ i4 h( m8 _9 K2 a" K2 t+ u
面实体——?有位置,具有特定的形状、范围?通常用来表示自然或人工的封闭多边形?用坐标对序列表示的闭合线来表示中国土地利用分布图(不连续面)
! B/ I0 z: _2 d: D二、地理空间数据的内容
3 s9 S# j3 k+ |" _( q0 G/ B地理空间数据一般包括三个方面的内容:1 B" x2 B& Q4 Z) o
1、某个参考坐标系所确定的位置7 _2 R/ O3 x) V
即几何坐标,标识地理实体在某个已知坐标系(如大地坐标系、直角坐标系、极坐标系、自定义坐标系)中的空间位置。
" r2 {0 E9 ]1 }0 P0 L2、实体间的空间相关性——空间拓扑关系
& a6 u2 F2 Y; r2 Y表示点、线、面实体之间的空间联系,如网络结点与网络线之间的枢纽关系,边界线与面实体间的构成关系,面实体与岛或内部点的包含关系等。: O3 t! I8 f6 Y' `0 I8 W
3、与几何位置无关的属性 即非几何属性或简称属性(Attribute),是与地理实体相联系的地理变量或地理意义。定性:名称、类型、特性… 定量:等级、数量…
7 i& }/ y/ d, z0 o9 @" ?地理空间数据三个方面的内容可分为两种数据类型:
' E5 g/ m. J* U! e; l$ a- b7 \" V1、某个参考坐标系所确定的位置——必需的 a' H" T4 k# M, H
Q6 R4 E, W) g, y. _7 v7 ^2 m
. E6 o& {, D. R1 D7 y7 q$ y- g ?, E# D8 V0 o) a' S0 l
三、空间拓扑关系
6 ]* V, O H6 B! G地理实体间的空间相关性,即点、线、面之间的空间关系称为空间拓扑关系。 “拓扑”——源于希腊文,原意:形状的研究 “拓扑学”(Topology)——几何学的分支,研究内容:图形在拓扑变换时能够保持不变
: @8 u- ^5 v5 [! ^) ~- C2 a/ y- |的几何属性——拓扑属性4 y r" O6 O& V+ ^6 b9 w
举例:将橡皮拉伸或压缩
& q/ Y. `1 G( F) g) w% P. `# n) m?空间拓扑关系的分类:7 i. R* l0 A* E' e: O; p3 s5 l8 t
地理实体间的拓扑关系可以从不同的角度来描述:
. C: s2 p6 V% s( l' W5 l( b( K空间位置 实体的构成
3 Y# F, T6 \3 U+ K, m不同种类实体间的相互组合
( _1 k7 Q0 d7 Z) w. M$ R1. 空间位置的角度 a* m: I- O, z4 N+ x) E W
?拓扑邻接:同类实体之间的拓扑关系。
; X" R9 b* f5 A7 X' u; H' \7 N0 z6 I+ q% e+ n
如:N1/N2 、N3 、N4 P1/P3 ; P2/P3 ?拓扑关联:不同类实体之间的拓扑关系。5 D1 R& k' t$ o1 b" ~6 W& z, m
如:N1/е1、е3 、е6 P1/е1、е5 、е6 ?拓扑包含:同类不同级实体之间的拓扑关系。4 c; E$ D/ N- z
如:P3 与P4) c! I0 R6 z+ y' H* g
拓扑包含分类:简单包含:多层包含:等价包含:
- ~& i; O2 k8 ]0 Q2 N2 T2. 实体构成的角度0 i6 k9 z) x0 N0 `
e7 P4 e3 、 e4 、 e6
( \6 ~; J v% F4 kP3 e2 、 e4 、 e5 P2 e1 、 e5 、 e6 P1 弧段 面 …… …… N1 、 N3 e3 N2 、 N3 e2 N1 、 N2 e1 结点 弧段: ]* e6 f( v8 s
7 O8 X6 S1 m/ C) J6 ?3 W面实体的构成 线实体的端点
) x. X4 K- y$ K! G ( P) B1 H- p2 {+ U, k W
e1、e2、e5 N2点
; Y- ^' ^/ r+ a: B! {! ge1、e3、e6 N1 弧段 结点 …… P1. Z" h& N# }: Q [, F
右邻面 右多边形)
) c. l! ~! `; K' e* B0 u…… ……P50 g { m0 F s ~) T) F8 t3 v
e1
+ [1 i g$ G1 B' P' ]8 o左邻面 (左多边形) 弧段
; L- I( i" {; `$ m, `: Z 线实体的结点 面实体的公共边界
2 ^% _- r! ~: Q# y5 F, s3.不同种类实体间的相互组合0 F R9 Y3 T* _* W# y+ q9 U' u
1)点—点关系
5 [" n7 S( s$ V8 d+ F7 ?点和点之间的关系主要有两点(通过某条线)是否相连,两点之间的距离是多少?如城市中某两个点之间可否有通路,距离是多少?这是在实际生活中常见的点和点之间的空间关系问题。
9 u6 n4 U/ p$ J1 B8 I0 `2)点—线关系8 {5 P7 w3 ?0 Y# p- ?
点和线的关系主要表现在点和线的关联关系上。如点是否位于线上,点和线之间的距离等等。+ I9 ?1 D/ \/ G
3)点—面关系" u @+ X! g4 r3 g" T" C
点和面的关系主要表现在空间包含关系上。如某个村子是否位于某个县内?或某个县共有多少个村子?- a) Q3 Q8 V$ ~
4)线—线关系
' q/ s7 B1 W, l" J8 ]线和线是否邻接、相交是线和线关系的主要表现形式。如河流和铁路的相交,两条公路是否通过某个点邻接?9 O; _; Z6 f4 u& S
5)线—面关系7 f; U: X9 F; e) q; p2 W+ U
线和面的关系表现为线是否通过面或和面关联或包含在面之内?
( }- |( Q4 f8 Q$ |$ L6)面—面关系
) l' |1 q& m( o; H, l6 e. w6 u?地理空间数据拓扑关系的意义:. j# S9 N8 R5 v. {. x
(1)确定地理实体间的相对空间位置,它比坐标数据有更大的稳定性。——与坐标投影无关!1 d% _) i3 J0 F( D1 `- U# \5 y1 P1 k6 `
(2)利于空间要素查询。
! H$ J6 k" R/ y5 e' G! N(3)重建地理实体。
2 n8 [* Z/ V7 Q§3.5 系统开发、管理与使用人员/ `7 b+ c: ]0 G) R7 a
系统开发、管理、维护、操作2 g" K+ T3 F' M9 d, r7 N9 T
数据更新、系统扩充完善、应用程序开发人员
+ a+ E" m7 Z# ^/ W9 e& R1 B必备的知识:计算机、测绘、遥感、相关应用学科
) x6 z* X( R+ a- a# ^/ ^/ s. @/ U?GIS的使用年限和费用组成
! H! D* s5 C4 T+ z+ `, N4 B硬件:2-5年,10%- O7 e! j& a; c" t' Z
软件:5年,20%
1 P+ Z& _5 K5 m. l& a地理空间数据:10-20年,65%, o- ]+ p" `7 r; U2 ]- V# @
人员:5%(培训费)% m, q6 P' F) e$ j! v% v; k8 M
第四章地理空间数据的采集与处理: t8 Q# s; v0 Y' f
§4.1 地理空间数据的采集
. k* v' D6 n- g! ?' U2 C一、数据源的种类3 u! `7 v- r6 ^6 x6 w+ F! e$ B
1. 地图各种类型的地图是GIS最主要的数据源。纸质地图、电子地图?
! @" I% ?& x; [地图作为数据源应注意的问题:' L3 r! r1 U7 w+ D4 q9 j
a. 地图存储介质的缺陷——变形(纸质地图)
, `6 Z' E) V P& W: ob. 现势性差——制作周期长、更新慢
. s1 L$ [5 X& K3 x% y9 kc. 不同类型的地图采用的坐标系和投影方式不同——需转换
& W" i) d+ o- u: H: b# q' H2. 遥感数据/ u. m; y- L, Z/ r( T
记录形式:胶片、数字格式——影像合成3 d# I5 z/ `& a+ T- Y; Z
遥感数据的特点:快速地获取大面积、综合性专题信息。" `: P# \$ T+ `5 ?0 l$ |2 q
遥感作为数据源应注意的问题:各有其自身的成像规律、变形规律,应注意影像的分辨率、影像的纠正和解译特征。
! C6 f$ |" Y0 t! U/ Q5 B! w! k3. 测绘实测数据:?经纬仪?全站仪?GPS(单机、自差分、信标差分、广域差分)4 y4 ?) @* X( A( P
?GPS的数据输出格式NMEA-0183格式:文本格式
4 C$ F: P# n9 K& `4. 数字数据:7 r7 ~* u2 l" @
?其他GIS的数据?数字图件?数字资料
" e g3 r( Q! u使用数字数据应注意的问题:a. 数据格式的转换b. 数据的精度、可信度、时限
$ |4 z2 p9 T' ^6 \/ V5 E) T5. 统计数据:如人口数量、人口构成、国民生产总值等% u. u2 q3 c# C- m" R7 T6 s
6. 各种文字报告和立法文件
) h/ i; ]2 ~2 I3 A, r7. 多媒体文件:声音、图片、动画
R- E* h; T) v- \. n0 t二、地理空间数据的采集方法
: ?5 K3 F1 T% Y3 s6 e" \( S1.空间数据(图形数据)的采集
( ?6 d) U! q" o6 G W; U(1)手工输入:通过手工在计算机终端上输入数据,主要是键盘和鼠标输入。键盘:输入坐标数据鼠标:在计算机屏幕上直接点绘图形/ O, I! F+ t. ?: c2 i) s$ j
(2)手扶跟踪数字化输入:硬件设备——数字化仪。使用数字化仪的地图数字化方式,称为手扶跟踪数字化。数字化仪的种类:机械式、超声波式、全电子式。! a# A8 j* G, Y1 C) m# w, e& N- |8 A
生成图形的格式:矢量图形* n T8 i' Y c: p- S
?数字化仪的工作方式:% R0 m9 R9 n' l- X0 n* M
a.点方式:每次定标器的键被按下,感应板发送一对坐标数据到计算机。
4 m. R" M5 t7 E7 r1 db.距离流方式:定标器在感应板上移动一定距离,数字化仪就自动向计算机发送一个坐标数据。! }4 H0 q" L+ F/ C5 s7 w
c.时间流方式:不论定标器的键是否按下,数字化仪每隔一定的时间自动向计算机发送一个坐标数据。
' ]; K. @4 o* J, a9 w" w影响手扶跟踪数字化精度的因素:数字化仪的精度,人为误差,经验、技能工作态度。 (3)扫描数字化输入:. b- o9 ]" X0 V y; i7 m
硬件设备:扫描仪(滚筒式、平板式、CCD直接摄像式)。. c& C" t1 @1 K( u
描扫描模式:二值扫描、灰度扫描、真彩色扫。! d2 K6 o8 x! i
使用扫描仪的地图数字化方式,称为扫描数字化。5 W5 U6 `9 @" Y& v! Z# C
生成图形的格式:栅格图形" t# i) p0 T6 i5 F* u1 ^9 n
?扫描图像的用途:①用作底图显示②以之为地图,进行屏幕矢量化(自动、手动)③用软件自动转成矢量地图R2V——Raster to Vector。
: _# z& O$ ]! \( O (4)遥感影像处理和信息提取, B J& P& P& `/ _! Z" M
从遥感影像上自动提取专题信息。
; l0 H! e# y4 u# i以遥感影像为底图进行人工屏幕矢量化。
% }: v/ ^& C! s+ Z3 h(5)数字数据的转换:数据格式的转换、坐标的转换
4 M4 w; C3 z9 K( E/ b% _(6)定位测量设备的输入:GPS、全站仪等& c# ?7 N' h5 b1 m7 l, k
2.属性数据的采集7 A" d! {4 Y4 e0 O9 H# ]
属性数据即地理实体的特征数据,一般包括名称、等级、数量、代码等多种形式。4 @+ t/ A+ s# W' Q" `) Y& j5 Q5 z
属性数据的存储方式:1)单独输入数据库存储为属性文件,通过关键码(指针、索引)与图形数据相联系。输入方式:键盘输入或其他数据的转换。2)直接记录在栅格或矢量数据文件中。输入方式:首先对属性数据进行编码,然后再输入。
$ W' F# U0 A/ i' Q+ v8 e编码目的:将属性数据变为计算机可以接受的数字或字符形式,便于GIS存储管理。: ?) @# p; q5 X- X b
编码原则:(1)编码的系统性和科学性。(2)编码的一致性。(3)编码的标准化和通用性。(4)编码的简捷性。(5)编码的可扩展性。! F: s7 W3 T1 O ]
§4.2 地理空间数据(采集后)的处理& v% e7 z3 i' A' w$ K
一、数据处理的目的8 Y/ M# L, S: H6 O- v
1. 修正数据输入错误! O, ?. `$ p" f4 Y+ J" a
2. 维护数据的完整性和一致性
3 N; Y' K2 s) @/ {* u0 W) i二、地理空间数据一般性错误
( a8 }; Z+ R8 u9 h+ Z3 P1.数据不完整、重复
9 h! {/ `# |" t0 V2 K" m- H4 W2.空间数据位置不正确+ O3 I Y2 x0 J* }6 d- h6 Z* j
3.空间数据变形9 Q+ q: S4 s7 F! ]2 X! [; A3 ^
4.空间数据和属性数据连接有误8 P( v2 l* X& X+ n& G
5.属性数据不完整% j$ H, n m7 w9 c) R$ |
错误检查的主要方法:7 _4 X+ L" b4 {+ ~9 ?
? 叠合比较法
. x. l6 R8 B! P7 u- X? 目视检查法# G, ?; U: p* h' N
? 逻辑检查法) c8 |+ q1 I2 E/ U7 |( g" _
三、数据处理的主要内容
+ D/ a/ a# E, U" m4 d" H. p, L1、图像纠正
' X7 A! }) i9 q T扫描得到的地形图的纠正
+ O3 }' q! z* _+ d \. P遥感影象的纠正8 s& ^3 l9 [2 z4 |3 w2 q5 I
2、数据格式的转换
& I: {5 y* j7 a0 D不同数据介质之间的转换; g* D! v1 j% c5 ] t
数据结构之间的转换
; E, F/ M; d5 B3、地图投影转换
: Y; ]2 H q+ O8 i% ]! K
2 t4 N* v* k S0 o' ]' l8 ?3 R4 M) F0 i. {# s+ r, Q) x$ ]! ^
投影B(X,Y)- M: g5 U6 w7 {
4、遥感影像的解译* s0 q3 B8 ^& G5 Q
知识要求:
, W8 p6 B+ e$ K% S9 F! E, b% R6 t?研究地理区域的一般知识;
& c& T. e- @+ T- l. Q?掌握影像分析的经验和技能;% O% \0 Y. {& G; ~* ^9 R
?对影像特征的深入理解。
5 \6 \, G* s+ @ 目的:建立典型地物类型的遥感影像解译标志:包括图象的色调或色彩、大小、形状、纹理、阴影、位置及地物之间的相互关系等。8 j* H$ s& o0 H+ r4 s( t `
5、图幅拼接2 x+ ?$ x/ {* W/ U1 [* G& g
在相邻图幅的边缘部分,由于原图本身的数字化误差,使得同一实体的线段或弧段的坐标数据不能相互衔接,或是由于坐标系统、编码方式等不统一,需进行图幅数据边缘匹配处理。) o( [0 e. A5 D3 c
图形接边:实体的连续性; 逻辑接边:属性的一致性, }$ C3 n Z; W( b# r1 }' C$ G
6、拓扑生成
0 }& E+ y, U; H' q(1)图形修整:伪节点、悬挂节点((a)多边形不封闭(b)结点不重合)、碎屑多边形不正规多边形。
* G3 ?0 _& `# v; z+ U' _2 h3 y(2)建立拓扑关系
e3 I: H N9 F$ H5 p% p" x) I, j计算机自动生成:大多数GIS软件提供完善的拓扑功能;
) M& C8 P* u5 e* D( T6 F0 `7 u对计算机创建的拓扑关系进行手工修改。
+ E1 E& K: T( z ?四、GIS的数据质量
7 ~1 B% K/ @4 y& O8 c1.定义:空间位置空间位置、属性特征属性特征以及时间时间是表达现实世界空间变化的三个基本要素。GIS的数据在表达这三个基本要素时,所能够达到的准确性、一致性、完整性,以及它们三者之间统一性的程度,称为GIS的数据质量。
& K4 }! N1 j; q7 ?2 .GIS数据质量的基本内容
6 M( r. ?& Z9 K4 r! U) R3 U(1)位置精度:如数学基础、平面精度、高程精度等,用以描述几何数据的质量。
2 i. q+ ?4 p9 ^/ _6 R- d, B4 Y(2)属性精度:如要素分类的正确性、属性编码的正确性、注记的正确性等,用以反映属性数据的质量。
" J+ x/ H7 ~. V3 |* T, Z(3)逻辑一致性:如多边形的闭合精度、结点匹配精度、拓扑关系的正确性等。& v `1 S5 T$ q7 }0 C
(4)完备性:如数据分类的完备性、实体类型的完备性、属性数据的完备性、注记的完整性等。1 z0 N# Z0 H9 G4 ]! u9 _2 b
(5)现势性:如数据的采集时间、数据的更新时间等+ ~/ B# V5 @, [5 H3 }" N
3 .影响GIS数据质量的误差, _1 m% w: A$ O
(1)误差类型
+ K0 k L( `4 e) j# u5 M/ D! ]源误差:数据源的不准确性以及数据采集和录入中产生的误差
/ R0 b* J7 {5 r% V6 u2 G处理误差:对地理空间数据进行处理时产生的误差- u. A& L8 {" Q' `( B
(2)GIS中的误差传播误8 E1 O1 M, \/ x/ ^/ P- r# h
差传播是指对有误差的数据,经过处理生成的GIS产品也存在着误差。
- C+ }# u: s. U0 a2 x①代数关系下的误差传播:这是指对有误差的数据进行代数运算后,所得结果的误差。
7 A* s) e7 z/ }, a1 b5 {②逻辑关系下的误差传播:即指在GIS中对数据进行逻辑交、并等运算所引起的误差传播,如叠置分析时的误差传播。9 c1 X: k) [* h9 X4 V% ^
③推理关系下的误差传播:这是指不精确推理所造成的误差。
5 r* k2 ]* H2 n) _& S) s五、GIS的空间元数据2 E4 n1 h& O( ?: `/ f( z
元数据(Metadate):描述数据的数据,是关于数据和信息资源的描述性信息。% _' R2 W7 }8 ~: h
空间元数据(GeospatialMetadata):地理空间数据和信息资源的描述性信息。它通过对地理空间数据的内容、质量、条件和其他特征进行描述与说明,以便人们有效地定位、评价、比较、获取和使用与地理相关的数据。
+ ^- ?. w d; z+ e0 |! R9 A?空间元数据的应用
/ I: c" u8 G! P$ l+ o帮助用户获取数据、空间数据质量控制、数据集成
; C) t9 K; b! r1 V8 _: M6 n" J 地理空间数据采集及处理流程
2 o$ E6 v4 {+ m: t
) Q8 E. y6 |9 |# O3 X 第五章GIS的空间数据结构
6 o2 r+ B9 U! Y t- p. J§5.1 空间数据模型及空间数据结构 Q# w1 ?, B$ f# b: p. `4 d
一、空间数据模型7 v0 h( V) E8 D/ t
1.模型(Model)
3 h e$ l' _$ y0 Q h模型是对现实世界的简化表达,是现实世界的某些特征的模拟和抽象。1 s8 Z! f" @; H0 _3 P8 U
模型:5 w5 \' ]5 p# ?6 Z1 B1 W, t
实物模型:飞机模型
6 J5 K% E+ O' Y抽象模型:数学模型(应用数学分析方法建立的数学表达式)、数据模型
4 O; p% i; a3 [' @2.空间数据模型8 r* W9 K7 d! T8 ]- [' [6 t, J
地理空间数据特征的抽象,是对地理数据特征的抽象,是客观事物及其联系的描述。
! F& {# {8 ]6 }: v( ^3 c3.空间数据建模
b( m) ?" t- E?建模过程. B! b# X1 `$ k* ^6 L+ K
(1)选择一种数据模型数据模型来对描述现实世界的数据进行组织。
! t8 c% N$ z0 L(2)选择一种数据结构数据结构来表达该数据模型。
, j1 e8 d9 v- c' P(3)选择适合的文件格式文件格式记录该数据结构。- I, {- ^' [* v' E, W. s
把描述地理空间实体的数据组织成信息世界中计算机能够接受的数据集的过程。4 W- y5 v8 Y$ b) o9 L
数据建模通常可选多种数据模型,一种数据模型可用多种数据结构来表示,一种数据结构又可用多种文件格式来存储。
( f& [/ e$ T0 [$ T( }* ~栅格离散化——栅格模型
) l- {7 [$ ~, y. L矢量离散化——矢量模型
3 k6 w7 x! z! G0 l+ {8 t4 k) v5 p二、空间数据结构
) E4 H& `9 P/ V- Q) v1.数据结构. F: L, x; K" h0 X: |/ B1 z$ P/ S
数据结构是指数据记录的编码格式及数据间的关系描述,即数据的组织形式。
+ X) Q( X( o* L/ k/ O? 抽象数据结构(逻辑结构)
: t! c- L( |9 A" p1 F/ N? 数据存储结构(物理结构)
: \! C) i0 a+ [2.空间数据结构. f" B/ |/ r3 s+ h# |' P
空间数据结构是指地理空间数据在计算机内的组织和编码格式。它是一种适合于计算机存储、管理和分析地理空间数据的逻辑结构。是地理实体的空间排列和相互关系的抽象描述。
* Q$ L( m! c* G7 D6 v# ` A$ Q3 B- f4 A8 K
3.空间数据编码# b. \# k# U0 P3 f
根据选定的数据结构,将GIS的数据源转换为适合于计算机存储和处理的数据的过程,称为空间数据编码。
8 r _2 `5 r% \$ S图形数据编码 属性数据编码——矢量数据结构所采用" `# ]& {: e6 n9 S
图形数据与属性数据混合编码——栅格数据结构所采用/ N8 Y( J3 ]" B: k( }
§5.2 栅格数据结构及其编码+ Z+ N) _3 M6 j% K. c
一、栅格数据结构(Raster Data Structure)
; s* N, }; b t: A 1.栅格数据结构的概念
, i; `) C! j# V; E) q" C0 _6 L. q将地面划分为均匀的网格、每个网格作为一个像元,像元的位置由所在行、列号确定,像元所含有的代码表示其属性或仅是与其属性相联系的指针。用这种离散化方法获得的地理空间数据称为栅格数据。地理空间数据的这种组织方式,称为栅格数据结构。
- l6 V/ q4 S2 f, g) L: [4 F实质:像元阵列
2 o t: |7 [3 i9 ~& O) E4 M特点:定位隐含 属性明显
5 n" X6 `$ C, O# J# S7 ~( N?像元的表示:5 U" k* `/ e* ^. q
位置:像元的行、列号确定
" R! D/ v) q. @2 U7 S7 {: ]属性:像元的代码(像元值)确定% H; z" D! x! K" H" V
?栅格像元的分辨率:, L1 ]/ c3 C J. W9 L
每个栅格像元所代表的实际地块的面积的倒数。0 i! N2 d9 @* C( \
R=1/实际地块的面积; e5 ^( o1 @( @& p, Z: E. v
2 .地理实体的表示方法
. l3 M8 N5 N }" U实体:用一个像元表示。8 I8 L- j3 a! h! }5 @' |2 j. [, E
线实体:用沿线走向的一组相邻像元表示。
; B4 Z; U! m# X. Q) p. m% \4 N面实体:用具有相同属性的相邻像元的集合表示。
! C9 P- H1 b a4 D0 _3 .栅格数据的获取途径
% ]7 \0 \9 P; g1 |% n! G: d0 w) l(1)格网法:在地图上进行均匀网格划分,逐个网格决定其属性代码,形成栅格数字地图文件。' g4 b" |2 M7 n* U( _- r
(2)扫描数字化。6 T! C; L, R$ `9 a- b& i
(3)遥感影像数据。0 N" K% V& K+ o; b/ i# E. _! U& J
(4)矢量数据转换为栅格数据。
8 x6 b5 m" J' } ~, C& y" e8 G% Z(5)DEM数据. Z4 g- j" J; x0 O' R. C4 g4 f
4 .栅格数据的偏差
5 R4 E6 u0 F4 b2 q7 ^几何偏差
2 {9 N. d+ |0 I9 ^属性偏差
% h, l4 A! b% u4 S! m2 `(1)几何偏差(又称形态畸变):(举例:三角形栅格化)& y' E" m1 T+ K$ Q4 ]
降低几何偏差的方法: 提高栅格像元的分辨率
6 Y! x8 @6 |2 j2 j$ v8 i$ C) G(2)属性偏差# M2 \7 L2 d2 y: [2 P' a- e# M4 h: [* f& m
像元的属性是地表相应区域内地物类型的近似值,因而存在属性方面的偏差。: v0 b8 a! P! g; A/ l6 Q6 a
Landsat TM卫星影像的单个像元对应着地表30×30m2的矩形区域,影像记录的光谱数据是每个像元所对应的地表区域内所有地物光谱辐射的总和。: E" e1 m- H) ]2 |
?栅格像元属性值的确定:
0 i3 x$ L# h. H- Y( ?①中心点法
! s* d: R9 q+ z) k7 x2 S; |( h由像元代表的地块的中心点所在的面域的属性来确定。适用于具有连续分布特性的实体要素。8 w5 N" H" E9 @
②面积占优法
( B9 o* ^" W0 d6 ]; i8 h c由像元代表的地块中面积最大的实体的属性来确定。适用于分类较细,地物类别面积较小的实体要素。% l3 N6 u9 J2 S! N v: \
③重要性法
6 |4 Z Y+ _6 t: Z选取像元代表的地块中最重要的地物的类型作为栅格单元的属性值。适用于具有特殊意义而面积较小的实体要素。
1 E+ d/ P! o+ R d* U/ [? 降低属性偏差的方法:
" S- H7 g0 @$ w9 {+ A4 ]# v 使用上述属性值确定方法
1 ]" o0 S- g: e* ^! k; l- ?提高栅格像元的分辨率
: h% Q1 { I# d% Q+ x4 m使用数据分层处理
# D5 q7 b# C1 P二、栅格数据的编码方法- V) `. K( |: w
根据选定的数据结构,将GIS的数据源转换为适合于计算机存储和处理的数据的过程,称为空间数据编码。, m0 V$ W7 @7 Q0 o' [: x
直接栅格编码:栅格矩阵编码) s* r! g* ] i: R% w1 Q
压缩编码:链码、游程长度编码、块码、四叉树编码5 P5 ~9 R; j* }6 \4 y
1.栅格矩阵编码(直接栅格编码)$ e' z8 ~7 ?& Q6 d# o
栅格矩阵编码就是将栅格地图看作一个数据矩阵,逐行(或逐列)逐个记录像元代码,形成栅格矩阵文件。
# g" w8 G9 W8 q- y8 V* c# g栅格矩阵文件是其它压缩编码方法的逻辑原型!
) d7 z+ { S9 S) N R优点:编码简单,信息无损失# Z! }7 }" G" i3 |- i
缺点:数据量大6 c3 D7 M, ?, H8 U: Q! T: }
2.链码/ F2 O* K3 ?# O0 f3 ^$ {
又称为弗里曼链码又称为弗里曼链码(Freeman) Freeman)或边界链码。或边界链码。用一个起点(行、列值)和一系列在基本方向上的单位矢量记录线实体或面实体的边界。
4 E3 ~7 K( D" o7 U1 O0 N/ S链码的特点:
, f2 Y5 H) P" [$ x+ J1 f# w(1)易于估算面积、长度、转折方向的凹凸度
6 z2 {, Y6 w; S3 W! |. w8 D8 M R(2)边界合并和插入等修改比较困难
) u) o6 t4 x, p# k5 L( F/ H(3)相邻区域的边界被重复存储而产生冗余
" S& ^; ~3 W/ C. C3.游程长度编码
' M% d! g; W6 p* u$ q基本思路:
$ g9 {, R0 D7 ]8 c1 b, l对于一幅栅格图像,常常有行(或列)方向上相邻的若干像元具有相同的属性代码,因而可采取某种方法压缩那些重复的记录内容。
. a$ k1 S( q0 B9 c' M$ |. o) N游程长度编码有两种编码方案
. C: Y# ]5 p. d H) U' P4 H, w方案一:只在各行(或列)数据的代码发生变化时依次记录该代码以及相同代码重复的个数。 方案二:" y$ f0 p1 j6 e( k7 G
逐行(或列)记录代码发生变化的位置(列号或行号)和相应代码。
; {9 M/ }. P q% E. v游程长度编码的特点1 J# p* F& c5 X
(1)压缩比的大小与图的复杂程度成反比,在变化多的部分,游程数就多,变化少的部分游程数就少,图件越简单,压缩效率就越高。: E$ j' g% L/ o8 D
(2)游程长度编码在栅格加密时,数据量没有明显增加。
* \: @( t! r0 |, H5 f4.块码
& ~* W/ R, H4 I. C( b: P( W又称为二维游程长度编码。
* Q4 P2 i$ E) d* o0 h/ g0 o( O$ p8 W3 l把栅格矩阵划分成由属性相同的像元组成的正方形,每个正方形称为一个记录单元,正方形行(或列)方向上像元的个数称为记录单元的边长。; |9 p4 W+ t6 \' t) J! T: a- k
记录单元的划分原则:按行(或列)扫描栅格矩阵,视其周围的像元能否组成记录单元。若某个像元已成为记录单元的组成部分,则不再重复划分。栅格文件由记录单元的初, w7 c: f8 W9 _8 c! {9 E* p$ q
始位置(左上角像元的行、列号)和边长,再加上记录单元的代码组成。- S) T- y+ s3 _: u3 R$ |- }1 K
? 其特点与游程长度编码相同。
, f# C+ b* G9 p, }# Y& W5.四叉树编码
1 f0 m6 w- E' @% p! [7 V四叉树编码是将整个栅格图像按照4个象限进行递归分割(2n×2 n,且n>1),直至整个图像区逐步分解为一系列具有单一属性的的方形区域,最小的方形区域为一个栅格象元。
* k9 T& j5 P5 U 分割原则:将图像区域划分为四个大小相同的象限,而每个象限又可根据一定终止规则判断是否继续等分为次一层的四个象限。& v+ c& O c2 P, L- O7 ~
终止规则:不管是哪一层上的象限,只要划分到仅代表一种地物或符合既定要求的少数几种地物时,则不再继续划分,否则一直划分到单个栅格象元为止。
: g6 d% Q7 F% t 为了保证四叉树分解能不断的进行下去,要求图形必须为2n×2n的栅格阵列。
; X! L; v1 V) Q+ g* ??四叉树编码的记录方式:4 d5 Z4 v: J! ?8 X2 B
——以美国马里兰大学四叉树编码为例(线形四叉树)' t4 R7 I! p4 m" Q* U6 V, \$ y
(1)编码中只记录叶结点
# y, E! r+ O9 g$ ~(2)使用32位二进制数记录每个叶结点
' G9 H6 q/ w. Y3 ^9 u. U" E& z1 e(3)每个叶结点只记录其地址和属性代码 S9 y1 I; ^, P- V4 ?1 {
深度:即处于四叉树的第几层上
% J/ x+ b, _/ k9 t& b- U, k路径:每次分割所在的象限代码的组合% W4 p- ?" }+ w: s0 Q8 s
?四叉树编码的特点:5 s; \% S; |! h
①阵列各部分的分辨率是可变的,边界复杂部分四叉树分级多,分辨率高;多边形内部分级少,分辨率低;因而既可精确表示图形结构又可减少存贮量。8 Z' B' D0 F! h6 ^
② 树状表示的结构缺乏统一性。用同一形状和大小的多边形可能得出多种不同的四叉树结构,不利于形状分析和模式识别。
0 x9 G6 U3 S( ?# F f三、栅格数据的压缩性
& i" W. q* D3 o% u1.栅格像元大小的确定
) T! }$ I; U% \设图形(地块)内最小多边形的面积为A,则合理的栅格像元的大小(边长H )为:9 `2 |# b% D* A6 L: M6 O# q
H = A0.5/26 W/ T$ Y! p: o3 M
2 .栅格数据的压缩率; _- [; T7 C- t5 O, a% o6 R/ {/ U! _
设栅格矩阵由i 行、j 列像元组成,像元属性值变化的频率为q,则该栅格矩阵的信息冗余度为Rc:% R# h* D6 }* n, m. b
Rc=1- q/(i×j)
: y$ t1 n! P! g?Rc越大,栅格数据的压缩率越高
7 E( K4 v" X: |5 Q- F?Rc越小,栅格数据的压缩率越低 m5 W7 i4 s" D; h
?栅格数据编码总结:- ], T @$ I0 _% k2 g- ]8 ?
①栅格矩阵法产生的栅格文件是其他压缩编码的逻辑原型。5 [$ Z5 b8 d/ a4 A) e& {$ Q
②数据的压缩以增加运算时间为代价。
" }3 r: n* W1 ^( {! ^( B4 B③链码的压缩效率较高,已经近矢量结构,对边界的运算比较方便,但不具有区域的性质,区域运算困难;
! a/ ~) s0 l$ I$ U/ x- |: {9 m④游程长度编码既可以在很大程度上压缩数据,又最大限度地保留了原始栅格结构,编码解码十分容易;. y% U0 R: o) t' {& u" J' `$ d9 I
⑤块码和四叉树码具有区域性质,又具有可变的分辨率,有较高的压缩效率。0 D+ [2 y5 F) l# Y9 v- G* c9 K
§5.3 矢量数据结构及其编码, w4 v; B' M' B* ?7 Z7 f
一、矢量数据结构(Vector Data Structure)
: i# D( T3 X% y- g7 X% Z! i9 R4 F1.矢量数据结构的概念
6 A. F# v+ Q( ~4 `1 R2 x通过记录坐标的方式尽可能精确地表示点、线、面等地理实体,地理空间数据的这种组织方式,称为矢量数据结构。
J; P2 F! a0 Q( }9 x特点:定位明显、属性隐含
, e+ L0 S) ^8 E# x5 M$ i 2 .地理实体的表示方法* J" q6 o k1 @( W; p+ I
点实体:用参考坐标系下的一个坐标对表示。(x,y)$ o0 E% G- A6 a) q( n
线实体:用一系列足够短的、首尾相连的直线段表示,并依次记录这些直线段的首尾坐标——坐标对序列。(x1,y1)…(xn,yn)/ ^" I0 o/ n' a; A7 L% @
面实体:用闭合的线段圈闭成多边形 s: B, H' l; _0 T& M
?计算长度、面积、形状和图形编辑、几何变换操作中,效率和精度均较高。 c) h" z! G8 c1 w! I1 T( @
?在叠加运算、邻域搜索等操作时则比较困难。
9 g2 |0 [/ @ W3 .矢量数据的获取途径
! T; k. \$ k# n, g. o(1)手工输入2 L: [3 \% \& @3 ~$ n" c, J
(2)手扶跟踪数字化
^/ e% N% f9 Q/ K1 {(3)屏幕矢量化) d" Z8 s0 P9 x+ M* U
(4)栅格数据的转换
8 e8 W- q& A1 R$ m U0 Q(5)其它数据格式的转换% p4 {' s7 y- `$ Q- @" l2 |# A
(6)实测数据
5 M: O& E6 k0 _+ G+ @1 I二、矢量数据的编码方法
0 D' U6 o5 V. Z) ^) I0 x1.点实体的矢量编码0 p& n1 L3 A' ~; D. s
空间位置
V) f) d9 _9 i* t' y) R属性:非几何属性、显示属性0 C- {; H( I A: C) e/ b
点实体矢量编码内容:$ n+ y- J# k5 T, b/ T5 x9 J5 u) ~0 G* {
$ V( x( t% Y- G: ~+ g1 w+ y1 c2.线实体的矢量编码7 f- p% h- A) `1 T- y
, I9 i, @! x) k
3.面实体的矢量编码
* b' O. ^8 @7 D9 c! x4 o7 b q空间位置; }# o8 r3 ^- V& P4 r
拓扑关系(形状、邻域和层次结构)
/ C( { d3 Y, S$ b- n! k$ Q" m属性:非几何属性、显示属性
: @4 d6 t: `" ~! C1 h2 r. L面实体的矢量编码基本要求+ r+ o2 T# Q. R) e" u& Y
①组成地图的每个多边形应有唯一的形状、周长和面积。- y" v6 B2 f# |) E% R0 G8 d1 i
②地理分析要求数据结构应能够记录每个多边形的邻域关系。# Z5 J% j4 q0 @+ z2 S9 R
③专题地图上的多边形并不都是同一等级的多边形,而可能是多边形内的嵌套多边形(次一级)。5 L; C' q; d, B* Z' @
(1)坐标序列法
% F$ [) @% B) M& w+ Y# u7 Z# }) P& G0 f8 G% m2 z; E
2 b) l1 Y; S9 ]& V- v$ W?坐标序列法的特点
3 u& z3 m) R# M5 Ka.坐标序列法文件结构简单,易4 w3 {+ U' M" F
于实现以多边形为单位的运算和显示。2 v$ i; u$ m9 `
b.多边形之间的公共边界被数字化和存储两次,由此产生冗余和碎屑多边形;# h% y1 N: d! K, A. P
c.每个多边形自成体系而缺少邻域信息,难以进行邻域处理,如消除某两个多边形之间的共同边界;
) h5 _- v5 N3 o: |. V/ Y# Ed.岛只作为一个单个的图形建造,没有与外包多边形的联系;9 h/ g+ ?7 z$ D% M" x& t2 Y
e.不易检查拓扑错误。3 R1 d$ Z" K: ?) x2 J6 A
(2)树状索引法
0 a, F1 s e3 S8 b* L4 }$ x# R+ R% N 该法采用树状索引以减少数据冗余并间接表示邻域信息。
* R$ o$ x4 v/ g/ S# p1 s' y; h记录所有边界点的点号和坐标对,由点索引与边界线号相联系,以线索引与各多边形相联系,形成树状索引结构。
+ C# U2 o$ n* Y) t点号和坐标对——点文件
* S' L6 A/ ]2 n; U- l4 b Z3 V2 U点索引——线文件2 I1 z" I9 @9 S
线索引——多边形文件
( e- k$ r; `' X. [ v; K7 W3 J8 ]% e1 l4 z5 l: j
4 c k1 R& v0 f; c7 f?树状索引法的特点% x* \% \* M/ ?7 F
树状索引法消除了相邻多边形边界的数据冗余和不一致的问题。- {7 t& l& s! O( c9 I4 V% r3 N# Y
邻域信息和岛状信息可以通过对多边形文件的线索引处理得到。5 I r4 j$ k* N" s2 \
编码文件需要以人工方式建立,工作量大且容易出错。
, C8 x0 h+ _* [) \$ J6 @9 w(3)拓扑结构法/ z* y1 D# C0 i& Q# O3 y
要求:在进行矢量编码时,明确反映点、线、面的拓扑关系。; Z8 k$ z* h5 ^& u3 ?3 Q( m
典型方法:8 F- R8 k% D$ T+ W; Q9 w0 a
?双重独立地图编码(DIME)
1 Y! D* F4 j* U8 n& u?多边形转换器(POLYVRT)
7 T9 u5 t* a' J" J# K* y?地理编码与参考系统的拓扑集成(TIGER)) r7 O0 Z- X$ ?3 H. ^/ {1 H0 M
' V- x/ P+ z$ i( `0 |) F' i. S
?右图的拓扑结构编码可分两步进行
" Z$ x$ t) \" @6 P/ P① “多边形-边界线-结点”的拓扑编码表(大—7 R% w9 N1 h1 m0 J: X
>小). h" e! E: F" U8 q5 A( `" g( Q% n! X
多边形-边界线的拓扑关系
# n6 f- e- |0 M 多边形 边界线2 p: D$ q; H, [; @- a4 ^
P1c1 ,c3, c7: q5 J+ w8 D' U9 Q! @
… …
{6 N* k7 q+ z% gP4c2 ,c3, c4" K& ^6 B- C: k: [* [) a4 B" @
边界线-结点的拓扑关系! T% q. d' w5 \3 ?* a' i
边界线 起点 终点 边界点坐标
3 Y& \ @# g1 _" fc1N1N2
% |$ P: H: s3 r9 b9 o, b(x11,y11) …(x1n,y1n)
0 e2 a5 U; g8 P- e% i3 R5 ?& C2 l… … … … …
) b5 a7 M5 Q, Q* Q9 M- P; W; gc7N1N4
9 q% t7 u5 v* L9 A2 y9 H(x71,y71) …(x7n,y7n)! R/ {# ~1 J+ P' |4 c/ G
② “结点-边界线-多边形”的拓扑编码表(小—>大)
- J/ p# n! {0 r结点-边界线的拓扑关系
7 k: q# J# }1 R3 F5 l$ R结点 边界线
# O& v% M, T1 W7 JN1* ^1 T) I5 @! F- v) F
c1 ,c5, c75 |. [1 x! t& [/ c5 k
N2, Z7 B1 z q% M
c1 ,c2, c3
7 }" J2 E. G6 [: Z… …
. l4 [0 ]0 E( [4 G# q边界线-多边形的拓扑关系/ d5 M4 S( d4 U u i
边界线 左多边形 右多边形3 J: c* J& B+ E0 S
c1
$ t/ B+ e- U l9 W2 m3 AP1P3. M2 Q' q1 p" c. Y
… … …
7 F# w6 M9 j- Uc7" R5 I' e' q6 v3 ~* Y" @; P
P5P19 T3 h e: x+ ` ?
拓扑结构图:
% g* ^) u; d! V1 q0 V L8 Q2 I6 Z5 T, D6 _+ V( D
?由大到小和由小到大的拓扑关系应该是一致的,否则编码错误!+ U$ \5 [7 e7 _/ j! x' M
?DIME编码
# V& c- w! E7 ]$ L: `——Dual Independent Map Encoding 双重独立地图编码7 @' `1 k6 h: I a& {8 h' r
——技术贡献:拓扑检查* D% k) z: W& \1 @' u7 y
?拓扑结构法的特点
) X3 b7 `8 K+ w8 y6 I4 I6 k' N: b数据结构紧密,拓扑关系明确,有利于空间数据的拓扑查询和拓扑分析。. d3 o$ P R' q i4 e
公共边界只存储一次,降低了数据冗余。' y j3 K. \% W3 L
拓扑关系面向几何要素而不是地理实体,强调几何要素之间的连接关系,而忽视地理实体的完整性和独力性,因此难以表达复杂的地理实体,对实体的操作(增加、删除、修改)效率低,影像对实体的快速查询和复杂的空间分析。7 E( z) y9 E* T! P6 [9 e
数据结构复杂,不便于系统的维护和更新。
) R1 j' Q7 J t4 T8 J* N/ a& P: m
8 t! P# Z7 l! d面实体的一般矢量编码内容
. z6 K3 Y; m( \1 \?矢量数据结构总结/ l, t9 I8 r6 I3 F
① 用离散的点描述空间对象与特征,定位明显,属性隐含; S( w; l& {* P! E; R6 n, H
② 用拓扑关系描述空间对象之间的关系;% _5 P6 Y0 U) q* T, Q
③ 面向目标操作,精度高,数据冗余度小;
2 T& `, O W1 s. i④输出图形质量好,精度高。
# A$ T1 t% m8 ?$ S3 E( r2 P⑤ 与遥感影像难以结合;
- N6 O: Z4 ?2 j§5.4 栅格数据与矢量数据的比较和转换
& q/ ^2 l2 t! a! u一、两种数据结构的比较
& R3 U; _3 I2 _0 P; w3 C: ?" q6 I, v栅格数据 矢量数据( y0 `, _8 Y: R4 B( v4 j! H1 _
优点 ?数据结构及处理算法?数据结构及处理算法均较
: B) n/ Z. y% }! n2 D2 s) O简单;
7 f# J: a3 g- [# j4 y1 q?数学模拟方便;
9 g0 N: ?# \) D6 d( t0 G0 v?叠置分析易实现;
& o! ^0 y" s& g5 C3 u" W2 ~! ~$ b?数据输入费用低;
7 b2 M9 S& ^# k8 J?易于与遥感影像数据结合。6 J1 }7 G" |1 P2 |2 G6 f
均较简单;
. w7 z8 d9 B3 x+ M* ^6 X?数学?数据量小;?能完整描述拓' l9 g% [3 ^2 `4 R8 {& w
扑关系;
# p+ A# u+ u1 s2 {/ |' P?表示地理空间数据的精度高; `( B$ r- K" B5 ^6 g/ c* P
?网络分析易实现。?图形输出精
, s( N! l5 Z3 G" K' o. j( J确美观;: ? Q1 X4 C" {. [1 P# ~
模拟方便;( [% A6 P. N2 \; W: w
缺点 ?叠置分析易实现;?数据量大;
- O3 F7 r8 {- G- s% b# P+ |4 F, A?降低分辨率以减少数据量时,精度和信息量受损;3 O# ~% e& P+ L* ], T
?投影变换困难;& q3 ? d T9 q, H& N" x( J
?难以表达拓扑关系,网络分析难以实施; ?地图输出不美观,线条有锯齿。 ?数据输入?数据结构及处理算法复杂;
( X' r( O. s- L9 {, f" ??数学模拟困难; ?叠置分析难;
9 c" N5 s3 H& _9 c7 _. F' H?数据输入费用高;& ]2 D" h* u8 j4 U7 j' t7 q
?与遥感影像数据结合困难。 ?0 E: ]; v7 |% l6 z) h
费用低;
& i M+ ` x+ ^" g! p ( {) @8 y* H1 D4 B! P8 K7 D$ L
?数据结构选择原则6 F" p! C( {# A3 h1 E6 ^$ J8 s
可获取的数据源种类 所需的位置精度 所需的空间分析类型 输出地图的类型
1 w9 V8 N2 c( ~. H4 P8 SGIS 中需考虑两种数据结构转换的原因:+ C j! o% D# d
?栅格数据与矢量数据各具特点与适用性; ?空间数据库中两种数据类型共存;1 r/ L% Y6 x! Z+ _& ?' v9 Z
?建立同时适应两种数据结构的空间分析模型非常困难; ?GIS的空间分析需基于同一种数据结构。/ A- p/ E, G/ `& _0 b. ~
二、两种数据结构的转换6 g& F' A2 y" |9 o, t
矢量数据结构向栅格数据结构的转换 栅格数据结构向矢量数据结构的转换& G4 g8 Q5 _7 s$ _3 y
1.矢量数据结构向栅格数据结构的转换 确定栅格像元的大小:: T0 l. r# `- J; d: g
设图形(地块)内最小多边形的面积为A,则栅格像元的边长H 为:
9 A. H) @+ x: L& V7 Y3 Z+ HH = A 0.5. z( O) d1 D1 `. x; I: p
/24 s* A9 y' O8 c8 c2 m# S
7 F1 W6 H3 p$ \: ]
?点实体的转换
8 w% y1 K! g9 d) F, l7 c, S
# |: }0 u3 P) N$ ~: K; P! j 3 ^5 f3 z7 y$ D
?线实体的转换4 p1 x& d N9 m. h9 q1 o' ]$ V
线实体由首尾相连的直线段组成,假如某直线段的首尾坐标为(x1,y1)、(x2,y2)两点组& O. {. f. L e( A8 ~
成直线方程:y = ax + b ,以△x为步长(step),扫描直线方程,求出直线上某点坐标(xi,yi),用点实体的转换公式计算该点的行列号I、J,若与前一个点的行列号相同则不计,否则计为新的像元。扫描完所有直线段,则完成整个线实体的转换。: t# l7 Q. e6 Q5 k5 v4 J# `( K
" Z! ]' L5 F& B6 t5 }/ h?面实体的转换
- F1 v2 W. {* e) I实质:多边形填充(将多边形的属性赋予多边形内部的每个栅格像元)
6 B" u4 J/ O: i1 U常用的多边形填充算法:) _ z2 {) K; c! t. Z& i: a' J
1)内部点扩散算法
0 |5 Q$ G! I6 {$ y+ b$ Q" D由每个多边形一个内部点(种子点)开始,向其八个方向的邻点扩散,判断各个新加入点是否在多边形边界上,如果在边界上,则该新加入点不作为种子点,否则把非边界点的邻点作为新的种子点与原有种子点一起进行新的扩散运算,并将该种子点赋以该多边形的编号。重复上述过程直到所有种子点填满该多边形并遇到边界停止为止。& x" Q; Y5 n1 ? `+ m# Y/ O
2)射线算法6 G. e2 m: s7 i+ E/ O/ X2 b- _
由待判点向四周引射线,判断该射线与某多边形所有边界相交的总次数。( A k8 D1 E8 J. E2 `$ o( m G C7 z
如相交偶数次,则待判点在该多边形外部;如为奇数次,则待判点在该
) @ V) } F' y2 v7 D多边形内部,将此多边形的属性赋予该点;
- D3 m# ~# d/ ^6 g$ R, z9 i射线算法特例4 S3 S$ M( c }" J* K Z G
3)扫描算法
1 {% I) |$ k8 v扫描算法是射线算法的改进,将射线改为沿栅格阵列列或行方向扫描线,判断与射线算法相似。扫描算法省去了计算射线与多边形边界交点的大量运算,大大提高了效率。+ \9 b5 ?% M* ^( J i
4)边界代数算法
* `$ Z( R u; t如图所示,多边形编号为a。初始化栅格阵列各栅格值为零,以栅格行列为参考坐标轴,由多边形边界上某点开始顺时针搜索边界线,当边界上行时(图a),位于该边界左侧的具有相同行坐标的所有栅格被减去a;当边界下行时(图b),该边界左边所有栅格点加一个值a,边界搜索完毕则多边形内部所有像元被赋予属性a, 外部仍然为0,完成了多边形的转换。! Z! J" P& G. o' p: P4 ]9 z0 }
" l3 |8 f3 N; o/ k2.栅格数据结构向矢量数据结构的转换0 C7 @1 x' ^9 A c3 ?: C
? 点状实体栅格的矢量化
2 v& y* U( G6 ~* ]4 Z 点状栅格数据的矢量化是将栅格点的中心转换为矢量坐标的过程。5 p! ]3 N; B7 f
对于任何一个栅格点而言,将其行列号I、J转换为其中心点的X、Y的公式:, d9 \" S# z( q9 E1 c5 l# E- O
X = X0 + (J – 0.5)×dx
- W3 X3 t# V5 f- y7 c0 X1 Q4 c K4 SY = Y0 + (I – 0.5) ×dy
5 B2 w R5 V4 }9 ]? 线状实体栅格的矢量化9 \; v4 [! m( N6 @ o1 ~
线状栅格数据向矢量化数据的转换过程是提取弧段栅格序列点中心的矢量坐标的过程。 ? 面状栅格的矢量化9 @3 F( O- R; B7 ?' P: M
面状栅格数据的矢量化过程就是提取具有相同属性编码的栅格集合的矢量边界及边界与边界之间拓扑关系的过程。; f5 Y& [" F5 d" W3 E
? 多边形边界提取) t8 B5 [3 s6 [) T. z
? 边界线追踪$ \- [: ?: k a3 G
? 拓扑关系生成& Q! [6 \3 W) w* v* E/ f' q+ M9 D
? 去除多余点及曲线圆滑
% s. u, @% H. k! Y+ ^# z+ X4 s?面状栅格数据矢量化的一般过程:
6 ^7 ~$ L. B. K' v$ m- p; oa.多边形边界提取
4 u& p/ b$ z% c# @. p' p8 K+ D) M确定边界点和边界结点。) K2 c" m3 m2 a
b.边界线追踪. Y; w5 @; A- }- z2 w ?# y
对每个边界弧段由一个结点向另一个结点搜索,通常对每个已知边界点需沿除了进入方向的其他7个方向搜索下一个边界点,直到连成边界弧段。
7 z/ ^- L2 x" u# z# `c.拓扑关系生成
) @" v3 w% _8 r2 c$ ^( m: ^, C5 ]对于矢量表示的边界弧段数据,判断其与原图上各多边形的空间关系,以形成完整的拓扑结构并建立与属性数据的联系。7 S0 Q; N: e+ ]8 M
d.去除多余点及曲线圆滑8 S- D7 f& a! g
采用一定的插补算法进行光滑处理,常用的算法有:线形迭代法;分段三次多项式插值法等。: h( A1 K U# X7 _! P
第六章GIS的基本功能 o, B1 e5 ^4 ?3 e1 @
§6.1数据的输入功能
1 O: w6 s- N0 l, _GIS的数据输入需考虑的问题:
1 U: F5 [. t4 |! L) W统一的地理数学基础3 N6 m. i) E7 A
图形数据的采集
5 Z% B( a% ~( n* `属性数据的采集4 p4 F3 G1 v6 w" T# b
图形数据与属性数据的连接
" F% ~+ G/ ?+ M; T- l一、统一的地理基础% Y: p5 K% X4 k5 ]* p
——GIS空间数据库内地理空间数据的地理基础的一致性0 Q2 q" W4 I" R/ o( c; ^
统一的地理坐标系统——应与国家基本比例尺地形图的要求一致!
! J1 v2 W, d& O. E统一的地图投影系统——应与国家基本比例尺地形图的要求一致!
. f4 v; r; \' }: k# E5 l统一的地理编码系统——指属性数据编码
8 ?; p8 I, [3 x. X我国基本比例尺地形图:' g* X; p7 ^7 D% F y) x5 q/ B5 U
1:5000、1:1万、1:2.5万、1:5万、1:10万、1:25万、1:50万(高斯投影3°带)、1:100万(高斯投影6°带)# \$ ]8 U1 y6 v0 Y, [/ Y' F7 l
地理坐标系统:1954北京坐标系、1980西安坐标系、兰勃特投影
0 K% Z1 z3 V% c+ O4 l4 j0 t0 T! J影响地图投影选取的因素:
, Y' v5 I3 O4 r制图比例尺
/ k' @! ~/ L* q0 @* O" y区域的地理位置、形状和范围
& E" G# e' f Q/ t, c3 }二、图形数据的采集
& G) W4 f9 |* w: T7 `手工输入 ——键盘、鼠标5 F7 `6 c" d! p2 L% S; e3 o& \% z& \: I
手扶跟踪数字化输入 ——支持数字化仪
7 s+ j* M# ?' g6 N# C" u. @8 b% [9 G扫描数字化输入 ——支持扫描仪
; W% A7 U8 O4 w& Q遥感影像处理和信息提取——与遥感影像格式的兼容 空间数据库
. h& B$ D0 K4 ?" m数字数据的转换 ——与其它数据格式的兼容* z9 R, J. _- F$ G+ a. ]0 X2 v
定位测量设备的输入 ——测量设备接口
- O( V# `$ |0 s三、属性数据的采集+ u2 u U3 Q4 W! `( M
属性数据的存储方式 采集过程6 c+ d1 f# }. ]- f2 f
1、直接记录在栅格或矢量数据文件中 属性数据编码
0 C+ }: v( \/ w+ Q+ d1 e' p$ B空间数据库' J- C! {- O- J4 B( D. G
2、单独输入数据库存储为属性文件 通过关键码与图形数据相联系
z5 k+ a$ o& Y$ q5 E- G四、图形数据与属性数据的连接0 F7 y- x2 J F" \
——地理编码(Geocode): t3 x: X' G, F" ]% H
§6.2数据的编辑功能
' n S0 K! W" Z; k5 i) i5 [图形数据的编辑4 u. K v$ H8 v% ?" {: b
属性数据的编辑! {/ I% r: ~# d+ r& l
数据的更新
/ x% F4 Y$ p4 m% y一、图形数据的编辑
% Z7 N8 S/ c: [ a.误差或错误的检查与编辑
$ k+ ~9 \" ]/ o: e1 sb.拓扑关系的建立
J/ v- h- w) s* N Z* y2 ^7 Bc.图幅拼接% T+ P' z- o" v& e, O
数据采集——>空间数据库——>d.遥感影像的解译 ——>空间数据库
+ m3 L4 X7 ^5 b% qe.地图投影转换数
" P/ h" Q( {6 \* A0 mf.据格式的转换
9 K/ J) v, i/ ]7 @: i: Wg.图像纠正
- f' q* S+ ^9 @" |& D/ T二、属性数据的编辑
# \- I' ?" E7 f& |数据内容的编辑/ I: u, q1 ^! L& F; e7 b
与图形数据连接的编辑, |0 i# A. A4 d+ [& j, P
统计计算与分析报表
- _. |# B; l) p- ?三、数据的更新——保持数据的现势性!
( @" Y: T+ d* ~图形数据的更新、属性数据的更新
9 \+ Y6 g5 O. `. W* R§6.3数据的存储与管理功能' L; w: Q& Z5 e
一、计算机数据管理模式及其发展
6 v$ A3 ^8 h9 L% \1. 程序管理阶段程—20世纪50年代% |; ^1 P0 \) P: T- ?1 {
数据面向应用,程序间不能实现数据共享。
7 N5 D7 u/ ]9 b% Z2. 文件管理阶段—20世纪60年代0 M0 v* b7 X& p( t- h
?由文件系统管理数据,应用程序与数据的逻辑结构有关联;. A2 M$ i. Y6 y- L# q% w: I/ o% Z
?应用程序直接访问所使用的数据文件;
! J0 e% n2 F N& ?8 g?数据共享性差。
& C6 |/ h- `' b! \+ B0 G3. 数据库管理阶段—20世纪70年代+ F% `/ t+ I9 h* E8 V, S# V
?数据由数据库管理系统(DBMS)统一管理和控制;' N# N5 g+ H& j4 B2 }3 q- \/ |- E" o
?DBMS提供方便的用户接口;
5 L K. u: l4 R1 h/ f$ `) a?数据的共享性、独立性高,易扩展,冗余度低% o S! y; M& g2 I/ G0 Z
二、数据库的相关概念
5 I7 t) y. b9 V) j7 c0 K3 J8 [1.数据库
2 N3 i: j6 T: @1 H$ v3 L, g为了一定的目的,在计算机系统中以特定的结构组织、存储和应用的相关联的数据集合。 2.数据库管理系统(DBMS)8 [" i7 v0 i: v" u# o
位于用户与操作系统之间的一层数据管理软件。它提供数据库的访问接口。$ \2 l) R$ ^% A9 b' }
DBMS的功能:数据库的定义、数据操作、数据库的运行管理、数据库的维护、数据库的通讯$ k% }" y* g: [
3.数据库系统(DBS)+ s; H6 P% o- S4 o9 l) F2 Y% G+ D
指使用数据库管理数据的计算机系统。包括:数据库、操作系统、DBMS、应用系统、数据库管理员和用户。
* @1 `% z1 W" y7 v1 o- p# g三、数据库的数据组织的分级# R9 q0 F6 ^$ e! ]' R, t& v
1.数据项:数据项是可以定义数据的最小单位,也叫元素、基本项、字段等。
0 V p' r2 ^- i, d0 d$ A7 q/ O8 y8 F2.记录:记录是由若干相关联的数据项组成,是处理和存储信息的基本单位,是关于一个实体的数据总和,构成该记录的数据项表示实体的若干属性。
- \8 _5 f- d% T- R7 n2 w3.文件:文件是一给定类型的(逻辑)记录的全部具体值的集合,文件用文件名称标识,文件根据记录的组织方式和存取方法可以分为:顺序文件、索引文件、直接文件和倒排文件等。6 r ~. W @1 W M/ K) u
4.数据库:数据库是比文件更大的数据组织,数据库是具有特定联系的数据的集合,也可以看成是具有特定联系的多种类型的记录的集合。数据库的内部构造是文件的集合,这些文件之间存在某种联系,不能孤立存在。, R Q2 |3 `( P' t
四、传统数据库系统的数据模型5 M% ]/ b/ m2 V
数据模型是数据库系统中关于数据和联系的逻辑组织的形式表示。每一个具体的数据库都是由一个相应的数据模型来定义。数据模型的主要任务就是研究记录类型之间的联系。 _3 ]- X* i6 Q4 N
主要有三种类型:层次模型、网状模型、关系模型( V( ^+ c) i3 \
1.层次模型) S" P! w. j ~+ Q( I
层次次模型所表达的基本联系是一对多的关系,它把数据按其自然的层次关系组织起来,以反应数据之间的隶属关系,将数据组织成有向有序的树结构。/ R. U5 Z% k. `7 a7 y9 A( m
优点:模型层次分明、结构清晰,较容易实现,反映了现实世界中实体间的层次关系
# ^. V T- c, c0 ?缺点:?对任何对象的查询必须始于其所在层次结构的根,使得低层次对象的处理效率较低,并难以进行反向查询。?数据独立性较差。删除一个记录其所有下属记录也同时被删除。' d2 e! Q n: U, Y( b' U% J( s
2.网状模型
" r9 X2 Z- w' Q7 E网状模型基本特征是在记录之间没有明确的主从关系,任何一个记录可与任意其他多个记录建立联系。网络模型将数据组织成有向图结构。
; r8 F, m- w4 j% @优点:?与层次模型相比较,大大压缩了数据的存贮量,数据存贮效率高于层次模型;?可以表示实体的多种,关系更为灵活,对确定的数据表示效率高,冗余小;?表示关系复杂的地理数据和具有网络特征的地理实体效果较好。' j p8 U* K1 _! |9 b! \
缺点:?结构复杂,用户查询和定位困难。?不直接支持对于层次结构的表达。?数据更新较为繁琐。3 @" b) D, v* F8 P2 k$ r$ x' q) u
3.关系模型* n& _3 a& ]* c' t
关系模型是将数据的逻辑结构归结为满足一定条件的二维表。实体本身的信息以及实体之间的联系均表现为二维表,把这种二维表叫做“关系”。这些关系表的集合就构成了关系模型。
2 w k, x" e ], M# p优点:?数据结构灵活、清晰,可以通过数学运算进行各种查询、计算和修改;?数据描述具有较强的一致性和独立性。$ W; w' h* {3 M9 h8 F1 ]
缺点:?实现效率不够高,重构数据的操作费时;?不允许嵌套元组和嵌套关系存在;?模拟和操纵复杂对象的能力较弱。: _5 |+ w) B& Q+ X) c) _
五、空间数据库
6 \1 f l% B; Y$ _$ k. `; u: @1.空间数据库特点2 D+ P' z4 d7 t1 J; n6 e
空间数据库与一般数据库相比,具有以下特点:+ J' Z$ C& N( U
(1)数据量特别大。——海量数据
' r; E2 R5 Q4 w" j(2)不仅有地理要素的属性数据(与一般数据库中的数据性质相似),还有大量的空间数据,即描述地理要素空间分布位置的数据,并且这两种数据之间具有不可分割的联系。1 K) Z4 N! q! S0 |' `& Q4 I
(3)数据应用广泛,例如地理研究、环境保护、土地利用与规划、资源开发、生态环境、市政管理、道路建设等。5 \8 W! m# |' h( S v2 T. V
2.传统数据库系统管理地理空间数据的局限性2 m% |2 ?5 Z. ~
(1)传统数据库系统管理的是不连续的、相关性较小的数字和字符;而地理空间数据是连续的,并且具有很强的空间相关性。4 ]$ h' u+ k7 c* y* T# M
(2)传统数据库系统管理的实体类型较少,并且实体类型之间通常只有简单、固定的空间关系;而地理空间数据的实体类型繁多,实体类型之间存在着复杂的空间关系,并且还能产生新的关系(如拓扑关系)。
- \3 S) j, }+ t (3)传统数据库系统存贮的数据通常为等长记录的数据;而地理空间数据通常由于不同空间目标的坐标串长度不定,具有变长记录,并且数据项也可能很大,很复杂。6 S( r N# c- T3 o$ g' k1 o
(4)传统数据库系统只操纵和查询文字和数字信息;而空间数据库中需要有大量的空间数据操作和查询,如相邻、连通、包含、叠加等。
; w2 ~8 @) q. v* h' R3.GIS的空间数据管理方案
" \! _0 l, B# F$ g/ a3 n9 y目前,大多数商品化的GIS软件都不是采取传统的某一种单一的数据模型,也不是抛弃传统的数据模型,而是采用建立在关系数据库管理系统(RDBMS)基础上的综合的数据模型,归纳起来,主要有以下三种管理解决方案:
8 h9 R" s' e2 X% e# K9 P1 F! Y文件关系数据库混合管理方案; [5 i& K9 h. t3 y; J: E# {3 R6 X
全关系式数据库管理方案
. t- G1 l, `# z6 ^7 i" o- P扩展关系数据库管理方案
0 ~; i) o) J5 ?2 g文件关系数据库混合管理方案' p' ]7 N; v1 \* n0 Q
图形数据和属性数据分开存储,数据的完整性有可能遭到破坏。
/ o1 M" d: s$ u, @属性数据建立在RDBMS上,数据存储和检索比较可靠、有效;
. @8 i9 i& |/ w图形数据采用图形文件管理,功能较弱,特别是在数据的安全性、一致性、完整性、并发控制方面,比商用数据库要逊色得多。
1 D/ f1 S: a7 d! O# ?# R' S1 c, RGIS软件:Arc/Info,MGE,SICARD、GENEMAP等。
* |6 z5 d$ b: n9 p) `全关系式数据库管理方案" p, F, z! U& m! V: v
属性数据、图形数据同时采用关系式数据库进行管理' C- R1 F$ |! M. i
空间数据和属性数据不必进行烦琐的连接,数据存取较快
; F5 E A/ U( Q# q8 Z( N属性数据间接存取,效率比DBMS的直接存取慢,特别是涉及空间查询、对象嵌套等复杂的空间操作' ?" ^# y4 @! `
GIS软件:System9,Small World、Geovision等
, V( Y, m* {( m扩展关系数据库管理方案+ i2 r( J) e5 N0 G) G$ ?; h7 L
在标准的关系数据库上增加空间数据管理层,即利用该层将地理结构查询语言(GeoSQL)转化成标准的SQL查询,借助索引数据的辅助关系实施空间索引操作。
+ e5 ~1 ~; D1 @解决了空间数据变长记录的存储问题,由数据库软件商开发,效率较高
5 n! @3 \1 w" u+ L+ [; K: a用户不能根据GIS要求进行空间对象的再定义,因而不能将设计的拓扑结构进行存储 GIS软件:TIGER,Geo++、Geo Tropics等
5 `/ Q" a2 K# D5 eMapInfo的数据库方案——。。。2 _$ F9 J! A1 P: m- T6 b/ k% }& p
对空间数据库的总结:2 s3 e$ ?( m. t7 Q/ r8 C, M
?空间数据库系统必须具备对地理对象(大多为具有复杂结构和内涵的复杂对象)进行模拟和推理的功能。
' c. @2 P( f3 J" L @/ G- e0 h?可将空间数据库技术视为传统数据库技术的扩充;. H4 O7 A L9 x5 z; [; G1 N) _
?空间数据库突破了传统数据库理论(如将规范关系推向非规范关系),其实质性发展必然导致理论上的创新。
( ^, ]( H: h8 W§6.4空间查询与空间分析功能
1 j) D- _0 g* |2 u, {①物体A在何处? ——空间查询
- r7 \2 Q- P, x, o/ P: K ②物体B是做什么用的?——空间查询
" Y' V- d. n i1 j③物体A处于物体B的什么地方? ——空间查询
: g) Q3 h( G. n$ S④具有属性A的地块的面积是多少? ——空间量算
" ]+ R4 V' x- Z⑤B地物的形状怎样(面积周长比)? ——空间量算
8 j! j' I0 h. t& ? {# r⑥各种空间数据的交、并运算结果是什么? ——空间分析
1 m5 o3 d; d4 @) t, P# q⑦各种空间数据叠加运算后结果是什么? ——空间分析3 t# i0 p' X1 L) @4 x" ^, W
⑧具有一定属性组合的物体与什么物体相邻? ——空间分析
4 z/ V1 H0 F1 I! |6 {; }⑨按某属性组合如何将物体重新分类? ——空间分析
' f% m9 Z6 g7 `4 S5 s: U$ j?空间分析和空间查询是GIS系统的重要功能之一,是GIS系统与计算机辅助绘图系统的主要区别。空间查询和分析的对象是一系列跟空间位置有关的数据,这些数据包括空间坐标和非几何属性两部分。& ?" c! t( ~+ z, U
?GIS的空间查询和空间分析能力来自于对地理空间数据的分析、变换功能。. D7 F, @; ]" C& Y
——数据变换与分析模块
2 c6 s% l& p& W; k一、空间量算5 T+ O; t) Y" u& B$ a% T
空间信息的自动化量算是地理信息系统所具有的重要功能,也是进行空间分析的定量化基础。其中的主要量算有:# i! T0 |* y9 Z% [9 B* t Z
?几何量算 ?形状量算 ?质心量算 ?距离量算/ ^% T) @$ u2 c* G1 N& P
1. 几何量算; f. ^) n% m; T
几何量算对点、线、面目标物而言,其内容不同:
# f: d" P* ~5 e6 U" P) t点状目标:坐标;
8 V- ~: `$ ?& N+ s2 Z线状目标:长度、曲率、方向;) X3 R: O3 X2 I7 f" j; U L% @
面状目标:面积、周长等
: [9 h0 i: Z& M5 o( k9 x9 ] ?线的长度计算:……
& K# x" w) V! H面状地物的面积:……4 {2 n* S( b8 R! \
几何交叉处理方法:沿多边形的每个顶点作垂直与X轴的垂线,然后计算每条边、它的两条垂线及这两条垂线所截得X轴部分所包围的面积,向右(左)为正,向左(右)为负,所求出的面积的代数和的绝对值,即为多边形面积。4 n+ B& D, D$ P1 {* D* V3 I
2. 形状量算
$ g( ?( ?0 Y3 M$ @) _地物外形是影像处理中模式识别的一个重要部分。例如海岸线的外形是岛屿的重要特征,森林中不同类型的土地外形对野生生物显得非常重要。目标物的外观是多变的,很难找到一个准确的量对其进行描述。$ b5 z. J! A9 X6 L0 k2 x3 d6 `
基本考虑:空间完整性、多边形形状特征。
& E1 h/ d+ \4 a) \: g第一种量算方法:首先统计线段的总长度得到多边形的周长,然后与同该多边形面积相同的圆的周长进行比较。将多边形周长与圆周长相除得到一个参数,用以比较多边形的不同。如:多边形长、短轴之比,周长面积比,面积长度比等。
; P& `$ _6 e7 E第二种量算方法:如果认为一个标准的圆目标既非紧凑型也非膨胀型,则可定义其形状系数为: r=P/(2*(π ? A)0.5)4 ?) N$ T b. g7 r0 \1 j
其中,P为目标物周长,A为目标物面积。如果r=1,目标物为一标准圆;r>1,目标物为膨胀型。
& c( _$ J8 }0 u( Z3 n. ~+ s2 B9 a" n% r3. 质心量算% O$ }/ m. W7 a. \5 c% }) W
定义:保持目标均匀分布的平衡点,一般为多边形的几何中心或重心。 计算公式:
7 A( i! s1 h4 k ∑∑==i i y )/1( x )/1(N Y N X G G& w$ a% g, p' y9 O$ m5 h X
?某些情况下,质心为加权平均中心& _9 H$ w; b0 f$ {# {: r& s5 q8 n
∑∑∑
- }. y3 t, g7 w# z3 d+ G/ b& A$ S∑=
( B* d7 ^7 ]- _9 ]" |) j=$ e- R4 a9 B- H4 ~' V9 A. s
i
" f5 K6 M( L5 n, \; Ki% U0 \+ G" j" E; E: O3 `
i G, a; v+ V% Q, b) x
i- Q0 h4 r' F( W+ L" p7 @5 D
i i G( q3 P' T' `6 n- Z$ Q9 d
w
5 f9 G( _0 L7 }4 E+ L" y# Gy w Y w x w X# i# L. G: z5 T0 |
应用: a.跟踪某些地理分布的变化,如人口变迁、土地类型变化等。& C. g5 L2 k6 B2 t( P3 g2 m) d {& }
b.简化复杂目标的模型建立等
! a# }0 ^( K: y0 b3 P7 l2 E& {9 O4. 距离量算:. V# y+ N0 X! B: ^* w' I
n 维匀质空间广义距离公式:: h+ ?' e" \1 x2 a2 |0 r" A
()()
9 T2 \% L! s6 H4 m! w' G0 ?; a[
( p6 j4 Z% y) Q]2 Q( M- A; |7 j0 M8 T0 o- Q
k6 _2 K* O4 d7 l: S) I
k j i k2 f2 C" A# U1 L3 f7 {
j i Y Y X X d 1?+?=
0 n; ^ i0 s T* `. [K 为阻力参数
" d5 `) c+ S8 W1 P" K+ U$ k a?当k=1时:
; }* i9 m# y) e1 G8 l6 d/ nj i j i ij y y x x d ?+?=——曼哈顿距离9 U( p; H$ F3 X; x8 r1 L I
?当k=2时:
4 C, A3 `$ q) K$ S; G22)()(j i j i ij y y x x d ?+?=——欧氏距离. j7 y) f$ J" f. w8 @
?当k=0.6时:! O1 G0 }$ r. X, R+ e: U. S$ H
[]
- n9 M) Y9 j c6
+ f( f: @+ r7 e- n7 O.0/168 v- ^9 n7 K3 F0 H; s+ m+ r6 U* q2 E
.06.0)()(j i j i ij y y x x d ?+?=——非标准欧氏距离1 b6 c9 G; |9 q8 n, _
二、空间查询 1. 简单查询& Q9 ~8 J& X" ?! N) ~0 E
?按属性信息的要求来查询定位空间位置,称为“属性查图形”。 ?根据对象的空间位置查询有关属性信息,称为“图形查属性”。 2. 拓扑空间查询面面查询:. j& I/ Q" T) n2 }& ]: e$ ~
面面查询:如与某个多边形相邻的多边形有哪些。 面线查询:如某个多边形的边界有哪些线。 面点查询:如某个多边形内有哪些点状地物。
+ l" Z, E7 E. `' ?4 ]/ o% Z P线面查询:如某条线经过(穿过)的多边形有哪些,某条链的左、右多边形是哪些。 线线查询:如与某条河流相连的支流有哪些,某条道路跨过哪些河流。 线点查询:如某条道路上有哪些桥梁,某条输电线上有哪些变电站。 点面查询:如某个点落在哪个多边形内。 点线查询:如某个结点由哪些线相交而成。 三、空间分析 1. 叠置分析8 `" K+ G/ N" h) W8 ?5 T- G w
在统一的空间参照系下,将两个或多个数据图层进行叠加,以确定地理实体之间的空间% u: l( u9 p6 ~( r' d+ B: ]
对应关系或生成具有多重属性特征的新实体。这种空间分析方法,称为叠置分析。
0 M. f# I' S- @7 \2 w实质:分析具有不同属性实体的相互空间关系和属性关系。
! n7 K. M$ X! E(1)视觉信息叠加
8 H* i$ g( h C7 \将不同侧面的信息内容叠加显示在结果图件或屏幕上,以便研究者判断其相互空间关系,获得更为丰富的空间信息。?点状图,线状图和面状图之间的叠加显示。?面状图区域边界之间或! a o" U& U$ C2 ?5 c1 h: w* R
一个面状图与其他专题区域边界之间的叠加。?遥感影象与专题地图的叠加。?专题地图与数字高程模型(DEM)叠加显示立体专题图。视觉信息叠加不产生新的数据图层!# F/ r/ j3 F7 F4 L5 I
(2)点与多边形叠加计算& F$ S+ a" W# ^1 C" R1 X% i
多边形对点的包含关系。分析点是否落在多边形内或某多边形包含那几个点,然后将点的属性赋予多边形或将多边形的属性赋予点。例如将一个中国政区图(多边形)和一个全国矿产分布图(点)叠加:?将点的属性赋予多边形——查询某省有哪几种矿产;?将多边形的属性赋予点——查询某种矿产在哪些省里有分布。点与多边形叠加不产生新的数据图层!2 R3 N& V8 m2 m) z. R& [
(3)线与多边形叠加( Q4 A% ]- C% @( Y, x; I
判断线与多边形是否相交,只要相交,就产生一个结点,将原线打断成一条条弧段,并将原线和多边形的属性信息一起赋给新弧段。——查询某多边形内指定线种类的总长度线与多边形叠加产生新的数据图层!
9 F9 \# V' u: @7 Z7 v# B% t(4)多边形与多边形叠加" G9 W# e- u5 i
根据不同图层上的多边形的边界的交点,创建有多重属性的新多边形。多边形与多边形叠加产生新的数据图层!
" ~5 k/ w. h5 {7 U2. 缓冲区分析" Q. a% u% a' j% {0 o* R- J
根据地理对象点、线、面的空间特性,自动建立对象周围一定距离的区域范围(缓冲区域),以分析该对象对缓冲区内邻近对象的影响程度。这种空间分析方法,称为缓冲区分析。
# A4 A: j7 s# q( b( @. \2 w- S1 B缓冲区分为:点缓冲区、线缓冲区、面缓冲区! V" b4 v8 g, U( O8 B5 Q
缓冲区分析三要素:* k- I! S4 W; u- E p! }* G1 [
(1)主体
j$ |7 U/ z: Y. n+ y(2)主体的作用条件7 p1 q- V$ k5 b# X7 J. Y
(3)邻近对象
. {7 ~' C+ f: H( |9 i+ V/ a4 U. f应用:分析某地理要素(主体)对邻近对象的影响程度、影响范围! x1 L2 E- g; C: T2 k
3. 网络分析
6 W" N8 }1 H3 `1 | b* I依据网络元素的拓扑关系(线实体之间、线实体与结点之间、结点与结点之间的连接连通关系),考查网络元素的空间及属性数据特征,对网络的性能进行分析评估的分析计算方法,称为网络分析。. a* V" y9 A k4 y( p4 i" a
网络的基本组成部分
0 M/ Y/ ]- }' s" a链(Link):供物体运移的通道
- X1 `3 B: v0 f' t; f节/结点(Node):网络中链的端点或任意两条链的交点。结点中又有下面几种特殊的类型: 中心(Center):接受或分配资源的位置,有源、汇之分9 Z% w0 O3 W- \0 [. g( x
站点(Stop):资源中转站
4 {2 Y' o8 V/ W4 {+ u障碍(Barrier):资源不能通过的节点
, @! k7 k& k% r0 q' n$ Z4 b% X拐点(Turn):资源流向发生改变处,出现在网络链中的分割结点上9 a l: X# z) z- r
网络的主要属性:阻强、资源需求量、资源容量
8 A( G1 s2 K# Q5 s应用:" H8 i) f! \7 {$ p8 ~, B6 \
选择最优路径、资源调配、地址匹配 n2 ^) X% p( u9 s
1)路径分析:求最优路径 ?最优状态:静态、动态;?最优目标:距离最短路径、时间最短路径、耗费最低路径;?最优数量:1条、N条
: ^* Q0 n6 h8 `0 q' f2)资源分配?负荷设计?时间和距离估算" W1 |: y |1 o i) Z7 K) j
3)连通分析1 n( w3 o9 m0 x# q' j; m d, Z
4. 地形分析——DTM & DEM. d# |( L" ` z" _4 e6 U7 G; n
(1)数字地形模型——DTM(Digital Terrain Model)4 `/ J$ k5 b O3 u2 w+ g4 Z8 \
数字地面模型是描述地球表面形态多种属性信息空间分布的有序数值阵列,是通过地表点的空间坐标及属性数据集合表示地形特征的地学模型。
0 {" [9 H& Z, `7 R: Z+ B, CDTM={Zx,y}4 \- G( d1 ]4 a- T2 Q9 a9 Z8 Q
其中,Zx,y为平面坐标为(x,y)的点的地形属性数据。地形属性数据包括:高程、坡度、坡向、坡度变化率。
9 `7 O( B% x ]8 J5 u* P1 s与DTM相近的术语:
. }9 @- L9 `. `/ F德国:DHM (Digital Height Model)
; [7 c! f# i) {, v$ q# y G% L4 P英国:DGM(Digital Ground Model)
& ~* \! T' H9 z4 }1 b0 N美国地质测量局USGS:DTEM (Digital Terrain Elevation Model)
' K4 l; n: i1 Q(2)数字高程模型——DEM (Digital Elevation Model)
1 R2 ~$ K: W5 e) e, x( a# P- r当数字地形模型(DTM)中地形属性数据为高程时,DTM称为数字高程模型。
* E" ]% ^; i1 c9 I. M: Q6 K数字高程模型是描述地球表面高程分布的有序数值阵列,是通过地表点的空间坐标及高程集合表示地形特征的地学模型。
* E$ I* h2 B' q3 ]DEM是建立DTM的基础数据,其它的地形要素可由DEM直接或间接导出!# a; `% J" c7 A3 n
国家地理信息的基础数据:
4 k& k- y8 e' [" v9 p0 T8 o数字线化图(DLG,Digital Linear Graphs)8 {+ ]/ A% d8 e
数字高程模型(DEM,Digital Elevation Models)9 X+ `0 c5 D) N
数字正射影像(DOQ,Digital OrthophotoQuadrangles)* ~; H9 Z4 l) {, s" s. l. l
数字栅格图(DRG,Digital Raster Graphs)
5 Z; ^; c7 h5 i$ {前三种构成了国家空间数据基础设施框架数据!
4 S" |" A$ k9 \! i, _/ N% T(3)DEM的主要表示模型* x7 }) ?$ {, z3 X- N% f* L
①规则格网模型:规则网格通常是正方形,也可以是矩形、三角形等规则网格。规则网格将区域空间切分为规则的格网单元,每个格网单元对应一个数值。数学上可以表示为一个矩阵,在计算机实现中则是一个二维数组。每个格网单元或数组的一个元素,对应一个高程值。$ {2 _2 J+ O b B% y
DEM={Zi,j}i=1,2,…m-l,m;j=1,2,…n-1,n
/ |8 j1 R* R- ?* ~2 u' Z2 y②等高线模型:每一条等高线对应一个已知的高程值,这样一系列等高线集合和它们的高程值一起构成了一种地面高程模型。* Q' U# \4 `" ]) H4 F- Z' y: @
?等高线通常被存成一个有序的坐标点对序列,可以认为是一条带有高程值属性的简单多边形或多边形弧段。0 o( U3 t3 O) ]# C j8 d
③不规则三角网模型(TIN——Triangulated Irregular Network)
- N" j( ^4 m( o) FTIN模型根据区域有限个点集将区域划分为相连的三角面网络,区域中任意点落在三角面的顶点、边上或三角形内。如果点不在顶点上,该点的高程值通常通过线性插值的方法得到.
/ v, B+ x5 T5 S4 V" q; {?TIN可根据地形的复杂程度来确定采样点的密度和位置9 b8 v0 k3 E& V9 \
?能充分表示地形特征点和线
6 Q+ y6 w5 }* r2 w7 ?8 ?0 q?减少了地形较平坦地区的数据冗余。
! K; [8 r# S% M% X, d(4)建立DEM的数据源及采集方法
! m' ]0 g) C/ ?2 J$ D①以航空或航天遥感图像为数据源由航空或航天遥感立体像对,用摄影测量的方法建立空间地形立体模型,量取密集数字高程数据,建立DEM 。
9 {, t3 T" ^) z. n, [?SRTM计划——雷达干涉测量
: C6 I8 Y! g: mSRTM Level-2地形高程数据集:数字商程数据是在任务中通过C波段数据采集而产生 w( h: d# C# p, f$ U! e
的。每一个数据点代表一个高程点,水平基准面是WGS84椭球模型。绝对水平和高程精度分别为20米和16米;高程数据间隔为1弧秒,约30米。
, X5 d7 o& C( N! Q! M②以地形图为数据源主要以比例尺不大于1:1万的国家近期地形图为数据源,从中量取中等密度地面点集的高程数据,建立DTM。8 d9 S- n) g; f; h+ V8 r
手工方法、手扶跟踪数字化仪采集、扫描数字化仪采集* h! L) }( I/ M; Z
③以地面实测记录为数据源利用自动记录的测距经纬仪(常用电子速测经纬仪或全站经纬仪)在野外实测,建立DEM。) M3 l, v6 h, Y) q! N
(5)DTM在地图制图与地学分析中的应用
~* F% H" d$ {4 p0 ~①利用DEM绘制等高线图
( p) B4 W. t2 d' L5 q. h) b0 o利用DEM绘制等高线图,是以格网点高程数据或者将离散的高程数据由栅格追踪法原理转换为矢量等值线所产生的。该方法可以适用于所有的利用格网数据方法绘制等值线图。 i2 V5 h: w% [; K* J
②利用DEM绘制地面晕渲图3 v" t0 |: T% u: P) b
晕渲图是以通过模拟实际地面本影与落影的方法有效反映地形起伏的重要的地图制图学方法。& w, i$ w: v. [
③绘制立体透视图反映地形的立体形态。' e# Q6 ~" K( j" y4 T
④通视分析
% ?0 R4 X ?" n1 M; x* F7 j⑤地形分析:
9 f3 Z0 u, A& \1 [8 F0 I1 c5 x1 Y9 @a.坡度和坡向分析
* i( k% W7 ]! O+ eb.地表粗糙度的计算
' {/ w( E7 l# N- L ~2 s) Q vc.地表曲率(坡度变化率)的计算: L6 \& z. ~$ M: h
d.谷脊特征分析地面坡度图地面坡向图
& X2 I4 |7 }* {⑥计算表面积、土方量、绘制地形剖面图4 x& A: r1 T k0 q
§6.5 地理信息的显示与表达功能
6 V4 a! b" Q6 K! S+ e8 ]! d地理系统信息产品输出
4 `' T0 r1 y; h3 ~+ S0 L6 G屏幕显示
& j' B9 R o* X4 X' {# V* _打印机输出
) n8 G6 K3 L) J5 _绘图机输出
, x1 b3 g/ q: p) g# x电子地图的编制
4 P o. t* L; b8 K8 o% L?信息输出硬件的特点4 h& z) {) W% o: t( [8 o; D! p
设备 图形输出方式 精度特点. W# L8 y3 g! d% W/ h0 E$ v! [
矢量绘图机 矢量线划 高 适合绘制一般的线划地图,还可以进行刻图等
$ q, ~9 h& v; Z0 ]! A! V% J特殊方式的绘图
: {" [, u. ?$ M5 O e8 }8 \喷墨打印机 栅格点阵 高 可制作彩色地图与影像地图等各类精致地图制
& B* F7 t5 U6 t& Y; V品。! h6 ^: U6 O2 L: ^
高分辨彩显 屏幕像元点阵 一般实时显示GIS的各类图形、图像产品。
: Z1 B' g# Q$ r3 z h4 S行式打印机 字符点阵 差 以不同复杂度的打印字符输出各类地图,精度( y4 D6 [6 {+ {* y1 K, ~
差,变形大。0 B# `2 X0 r! S1 T+ U0 y j# Z/ F
胶片拷贝机 光栅 较高可将屏幕图形复制至胶片上,用于制作幻灯片9 O2 ?6 J7 h1 e8 l" T* |: [
或正胶片。1 r$ Z4 }" l, b9 F1 U
?地理信息系统输出产品的类型" G' N* [( J/ w6 ~
地图、图像、统计表、三维多媒体虚拟现实、电子地图、专题地图
. Q7 @3 e) Z# ]专题地图:
4 d7 [9 ~# r' U8 z4 Z 专题地图是突出地表示一种或几种自然现象和社会经济现象的地图。按内容可分为三大类:
) d) g- g3 A+ `- U7 Q自然地图、社会经济地图和其他专题地图。自然地图:表示自然界各种现象的特征、地理分布及其相互关系,如地质图、水文图等;
* K& b8 @& J" t* Y! M社会经济地图:表示各种社会经济现象的特征、地理分布及其相互关系,如人口图、行政区划图等;! H3 j) g4 q: b! D! |
其他专题地图:指不属于上述二类的专题地图,如航海图、航空图等。
9 v0 H/ G) _6 u- `1 l电子地图:+ \, ~# A7 u& r& r& R5 ]+ y
产生于20世纪80年代中期,它随着数字地图及GIS技术的发展和应用,随着计算机视觉化研究的深入应运而生,侧重于空间信息的表现与显示。电子地图与全球定位系统(GPS)集成,在航天、航空、军事领域以及汽车导航中也发挥着十分广泛的作用。
# A( A' g$ b2 z4 p/ _电子地图的特点:( |" g* T, F! @3 [! V) ~
电子地图是以地图数据库为基础,以数字形式存贮在外存贮器上,并能在屏幕上快速显示的可视地图。电子地图的主要特点是:快速、动画显示;生产周期短;检索十分方便编辑修改容易;动态模拟;通讯传送简便;量测自动化。
+ G8 s% ^- K6 q J电子地图的生成:数据采集、数据处理、符号化& g5 p+ O) |' g! n" y2 h
电子地图的显示:目前电子地图的显示,广泛应用光栅扫描图形显示器。屏幕电子地图的显示,实质就是存贮于帧存贮器中的象素矩阵在监视器屏幕上的映射。显示模块包括检索方式、属性查询、静态显示、动态显示、图形缩放、翻页和产品输出等功能。
- M* d) o) `# J, z 第七章GIS工程设计
/ x9 g6 x4 i" V# ^( ?0 u§7.1 GIS设计概述
* e5 @ N8 y1 n! PGIS的分类:
. _7 e* d: W- r& L专题GIS
/ g' C: _9 u! [区域GIS
7 I9 A: u0 x# ?: ?* H5 K6 {+ CGIS工具(GIS TOOLS)——应用型GIS、工具型GIS$ r: S$ U( m0 F9 E9 v
工具型GIS:7 e) Z' b1 U# F& @
通用型GIS,提供通用的GIS操作平台,一般没有地理空间实体,由用户根据目的和需要,以其为基础进行进一步的设计和二次开发,达到解决实际问题的目的。如:ArcGis,MapInfo,MapGis) T8 G: A) q1 c* _ f
应用型GIS:" ?7 o. l. s# i T
以工具型GIS为基础,设计、开发的具有解决实际问题能力的GIS。
1 u% p$ R n8 @9 B0 D! q5 tGIS的特点:5 |# ]7 g; }& D9 G1 R+ O6 c. R+ g: ?
(1)多学科交叉的边缘学科,涉及的专业知识广。
0 \' b, {! K* ?, w(2)以空间数据为主,兼有专业属性数据,数据结构复杂,数据量大。
2 z \* V9 A% ]/ \ U6 h/ W0 e& `(3)以应用为目的,针对不同的领域有不同的GIS。9 G0 y; Z" U9 O
(4)以空间分析、统计处理为手段,辅助决策为目的。
; H! c6 m0 J, U$ P1 a% \8 [' JGIS的设计、开发和应用是一项复杂的系统工程,涉及“系统”的诸多方面:
1 B' _3 Y6 Z5 `* V4 B$ m(1)系统的优化设计4 M2 C1 `" f2 G2 |/ }9 v6 ]7 E+ t
(2)系统的控制运行, z/ U0 M: ^7 `2 g5 L
(3)系统的管理
9 |: F! k0 G$ z X5 ]$ M(4)系统资源(人、财、物)的配置与组织* `; O# K% t% @: m" C( a+ |
——一个完整的GIS的建立过程,不是简单的原理或方法的应用,而是基于系统化思想指导下的工程化的建立过程!8 {/ z- ?7 x+ g L1 J1 o( \. z8 E( v
GIS设计的模式:
+ V0 Q5 @& |2 O& \4 \+ c结构化的系统设计模式——Calkins在1972年地理数据处理学术会议上提出
2 ]8 ?* |' ~1 N5 i5 q4 A: C8 |+ e结构化——有组织,有计划,有规律的一种安排。
. x+ N! s1 q5 ~( k/ J结构化系统分析方法——利用一般系统工程分析法和有关结构概念,把它们应用于地理信息系统的设计,采用自上而下划分模块、逐步求精的系统的分析方法。
1 s4 D2 {4 c& x! _8 T结构化的基本思想:' \) Y2 v4 J; o7 n/ _
系统的观点7 [& Q0 K0 [! \$ A6 f2 X
调研的观点; X( g2 B" c; a! N0 H; r
结构化的观点$ V! ?& t% t+ G& t0 D8 R
面向用户的观点
" Q3 R2 t6 G+ k3 R0 M! v8 N系统的观点:* Y0 X0 Y! M; s! F0 f" K* F: o
首先要从总体出发,考虑全局的问题,然后再自上而下,一层一层的完成系统的研制,这是结构化的思想核心。
7 D* {9 `& m0 t) D6 g' b调研的原则:
# ^$ I2 Q% W# N+ W) m- w只有先进行调查研究,了解用户的需求,掌握必要的数据,才能进行系统分析,做出合理的设计。调查研究是GIS设计的基本原则。
/ F8 ?' ^6 D2 Z+ e5 T$ Q% T结构化的方法:
2 o A$ l! M- b# P: KGIS的逻辑设计,应采用结构化的方法,包括: 分析信息流程,绘制数据流程图; 根
4 S* Y k5 E2 V 据数据的规范,编制数据字典;; 根据概念结构的设计,确定数据文件的逻辑结构; 选择系统执行的结构化语言
; m1 }! L1 t) b* P面向用户的观点:& S' u5 ^) m9 }, R; v
地理信息系统类型的选择,很大程度上取决于使用部门的工作性质、工作领域以及该领域内的应用范围和应用期限。另外,用户的技术选择,用户的财力和人力的费用等,都是设计地理信息系统需要考虑的内容。
3 o, @- [- V4 A& `) M地理信息系统设计要满足的三个基本要求:
/ I; X# F9 v! M7 O- b加强系统的实用性& Q& e, ]: p8 z8 y
降低系统开发和应用的成本
6 j- H# u# g8 b& k8 G3 D8 z提高系统的生命周期2 q4 H) m( Y* v3 ], g" u
§7.2 GIS设计与开发的步骤1 Z( s. V% ]; M
GIS设计与开发的四个阶段:
8 a9 m# u) n; ~: b- N3 w/ U系统分析——“做什么?”" U) o- L) l8 i* S0 Z
系统设计——“怎么做?”) b2 Z# T( b5 R0 x" j8 m
系统实施——“做!”: j/ Q6 o& r/ R
系统维护及评价——“做的怎样?”
3 C. {: i2 r( D7 ~# `& N# D一、系统分析9 d3 |1 H3 s# ?
1、用户需求分析- E8 J+ ]* X" B. C: j
调查用户部门或其它有关部门对相应GIS系统的信息需求情况。
X2 p$ n( H% I1 W7 d从上至下调查用户部门各级机构在目前和将来发展业务上需要些什么信息;/ q( x3 @& m, O
从下自上调查他们完成本部门专业活动所需要的数据和所采用的处理手段,以及为改善本部门工作进行了哪些实践活动等。; e3 R& z0 l6 {0 N& ?" g7 p; X# l
收集他们对本部门的业务活动实现现代化的设想与建议。
$ m7 l) \& A5 c, T1 s可能碰到的用户类型:3 k9 s. c5 \$ U) j, P' P. E- o
(1)能明确提出要求。这类用户希望用GIS来实现现有工作业务的现代化,改善数据采集、分析、表示方法及过程,并用以对工作领域的前景进行评估,以及对现有技术方法更新改造等。4 ` U8 ?" O" u1 u2 d
(2)不十分明确提出要求,提出要求不合理,或者要求与投入矛盾极大。
6 a2 ]5 Z6 j- L; a4 A1 F(3)不能提出具体要求,只是跟风炒作。+ n% ~% r2 ]$ w5 j f
2、系统目的和任务
, ~( A* a- x3 M) p9 U①空间信息管理与制图;* V7 m/ j; [& m$ H
②空间指标量算;
a5 W% _# E" j7 w* S( p }! s. c③空间分析与综合评价;5 _4 I T9 o6 d% L1 i' z
④空间过程模拟。
+ J9 `8 D; B' i; f! n3、数据源调查和评估% b! q( R+ s% z; m
种类、数量、精度、收集方法——工作量估计
1 Z$ s" w- H0 x) I9 W4、评价地理信息系统的年处理工作量、数据库结构和大小、GIS的服务范围、输出形式和质量等。
/ | c9 |# }$ |& C5、系统的支持状况
7 ^$ g0 X# m8 C* a" ?, |部门管理者的态度
+ n6 e; J, g5 i+ O; y) O/ ?( S8 `工作人员对建立GIS的支持情况;包括有多少人力可用于GIS系统,人员的知识结构等;
8 J7 f6 |3 k% e! f! ~8 ~) E" H' X 财力支持情况包括组织部门所能给予的当前的投资额及将来维护GIS的逐年投资额等。
. n, f$ _6 A. b& Y* F根据上述调查结果确定GIS的可行性及GIS的结构形式和规模,估算建立GIS所需投资和人员编制等。——可行性分析9 t5 J7 v7 g- m& p X: J6 w. E
可行性分析就是根据社会、经济和技术条件,确定系统开发的必要性和可能性,主要进行①效益分析;②经费估算;③进度预测;④技术水平的支持能力;⑤有关部门的支持程度等。?系统分析的重要性:
% k: m. H/ q0 U8 M5 O2 n根据美国一些软件公司的统计资料,在后期引入一个变动比在早期引入相同变动所需付出的代价高 2~3个数量级。( T% A7 w6 k: L* w
二、系统设计4 r" K8 n. o# `& z6 M0 W
(1)功能设计:根据系统研制的目标,确定系统必须具备的空间操作功能。
& g9 z& d( Q/ ~3 H2 w/ [(2)数据设计:数据分类和编码,完成空间数据的存储和管理。; W5 T1 W7 f3 T
(3)应用设计:系统的建模和产品的输出。
* ?8 q$ @$ i8 a" H$ `/ N①硬件配置与选购;& X+ z5 p$ C& p
②数据库设计;& k/ P+ r& e5 @' _6 P, v: o
③软件设计/ q: d% p2 }, l9 d9 R) ?, k
1)模块设计2)代码设计3)数据库设计4)数据获取方案设计5)界面设计6)输入输出设计7)程序模块设计8)安全性设计9)实施方案设计
! m4 {+ p6 D( r& m U$ \三、系统实施
7 {+ K0 t) @% q3 @1. 程序编制与调试( n( n1 \1 Y2 L9 I4 Y8 ]0 \- L
2. 数据采集与数据库建立
- o* k$ [% w+ Z3 A m+ A( Z* U3. 人员技术培训,编制用户手册
% T; u' _5 [) h' h: Z3 i4. 系统测试
* s, e, r% l/ z! X/ q$ e四、系统维护及评价2 i7 z/ B- x' `
1 . 维护:〈1〉纠错〈2〉数据更新〈3〉完善与适应性维护〈4〉硬件设备维护
. j, r- @- g k- s8 V2 . 评价:& K: }0 ~, F: N8 o3 t
评价的主要内容:系统效率、系统可靠性、可扩展性、可移植性、系统效益
% z P) w9 v1 T) {& v* O" {系统效率:是通过系统的各种职能指标、技术指标和经济指标反映出来的。例如:系统操作是否方便快捷,系统出错率如何,资源的使用效率高低等等。1 S4 c2 [" [3 ?
系统可靠性:是指系统在运行时的稳定性。可靠性好的系统,一般很少发生事故,即使发生事故也能很快修复。可靠性还包括系统有关的数据文件和程序能否妥善保存,系统是否有后备体系等。& i, o6 k2 o4 T: n) B; d1 i
系统可扩展性:任何系统的开发都是从简单到复杂的不断求精和逐步完善的过程,特别是地理信息系统常常是从清查和汇集空间数据开始,然后逐步演化到从管理到决策的高级阶段,功能逐步增加和提高.因此,要在现行系统上不做大的改动或不影响整个系统的情况下增加功能模块,就必须在系统设计时留有接口,否则,当数据增加或功能增加时,系统就要重新设计,这就是一个没有生命力的系统。
1 J: Q/ x, X. r! y& K; \系统可移植性:是评价地理信息系统的一项重要指标。一个有价值的地理信息系统的软件和数据库,不仅在于它自身结构的合理,还在于他对环境的适应能力。这首先要求系统必须按国家规范标准设计,以保证软件和数据的匹配,交换和共享.此外,在地理信息系统软件的设计中,采用自上而下的结构设计,即首先将与软件,硬件有关的公用子程序列出,在列出与软件,硬件无关的子% L1 ~ N M" j' Z# h
程序,最后组成软件系统,也可提高软件开发的可移植性。
8 O3 u( e L: b* W系统效益:' ^; Y3 D$ `7 v4 b# C+ l' S5 H
经济效益---体现:促进生产和产值的提高,减少盲目投资,降低工时耗费,减轻灾害损失等.4 n) D0 E, Q% j f) S
社会效益---体现:信息共享的效果,地学综合能力的发展,应用系统的模型化,工作时间的缩短等。
- W; F& E- S `( M8 S§7.3 GIS的运行及管理( T2 N Y6 f2 e% Q
1.日常操作
, r) U1 Q# r% i4 f2 .数据更新9 E6 N6 k3 t5 x" y2 m
3 .功能扩充
. r) J3 a* v* ^, {GIS设计与开发周期的各个阶段:
$ Q* K3 T& `5 ~) ?& _; Q7 U
- r, H# }( F1 L/ K6 h' ~3 q5 I§7.4 应用型GIS的开发方法
. w- j/ Q$ V; j0 b5 |' b" i?独立开发
; b4 W, ?" f2 T3 t! H; {?单纯二次开发
, B5 @ y: l4 t8 i. w2 ~7 ]5 f?集成二次开发
; |- k+ \- S$ a6 f% k. K1 x一、独立开发:
1 l8 c" [+ M/ Q指不依赖于任何GIS工具软件,从空间数据的采集、编辑到数据的处理分析及结果输出,所有的算法都由开发者独立设计,然后选用某种程序设计语言,如Visual C++、Delphi等,在一定的操作系统平台上编程实现。这种方式的好处在于无须依赖任何商业GIS工具软件,减少了开发成本,但一方面对于大多数开发者来说,能力、时间、财力方面的限制使其开发出来的产品很难在功能上与商业化GIS工具软件相比,而且在购买GIS工具软件上省下的钱可能还抵不上开发者在开发过程中绞尽脑汁所花的代价。& H& U! {& {; z1 C i( C4 C/ c& S
二、单纯二次开发:
: q7 C! b& r" `! c2 A' {! G 指完全借助于GIS工具软件提供的开发语言进行应用系统开发。GIS工具软件大多提供了可供用户进行二次开发的宏语言,如ESRI的ArcView提供了Avenue语言,MapInfo公司研制的MapInfo Professional提供了MapBasic语言等等。用户可以利用这些宏语言,以原GIS工具软件为开发平台,开发出自己的针对不同应用对象的应用程序。这种方式省时省心,但进行二次开发的宏语言,作为编程语言只能算是二流,功能极弱,用它们来开发应用程序仍然不尽如人意。
% M: n# ?# V3 ^, M& L9 f& M三、集成二次开发 Y. ~9 }( Q' Q1 f, [, p, N/ I
指利用专业的GIS工具软件,如ArcView、MapInfo等,实现GIS的基本功能,以通用软件开发工具尤其是可视化开发工具,如Visual C++、Delphi、Visual Basic、PowerBuilder等为开发平台,进行二者的集成开发。——GIS控件
7 j" g$ C) K5 P9 r集成二次开发的优点: 既可以充分利用GIS工具软件对空间数据库的管理、分析功能,又可以利用其它可视化开发语言具有的高效、方便等编程优点,集二者之所长,不仅能大大提高应用系统的开发效率,而且使用可视化软件开发工具开发出来的应用程序具有更好的外观效果,更强大的数据库功能,而且可靠性好、易于移植、便于维护。这种方法唯一的缺点是前期投入比较大,需要同时购买GIS工具软件和可视化编程软件。
" w0 A$ @) X1 M( I1 Y" K6 h& e% n三种开发方法的对比:& [* W" j) a+ } l9 M! X4 X
应用GIS开发的主流方向由于独立开发难度太大,单纯二次开发受GIS工具提供的编程语言的限制差强人意,因此结合GIS工具软件与当今可视化开发语言的集成二次开发方式就成为GIS应用开发的主流。 |