点击上方“蓝字”关注我们
第三章 空间数据处理
(一)坐标变换
实质是建立两个平面点之间的一一对应关系,包括几何纠正和投影转换,他们是空间数据处理的
基本内容之一。
几何纠正是对数据坐标转换和图纸变形误差的纠正。
投影变换是指投影方式的变换
(二)仿射变换
仿射变换在几何上定义为两个向量空间之间的一个仿射变换或者仿射映射,由一个线性变换接
上一个平移组成。
是GIS数据处理中使用最多的一种几何纠正方法。
它的主要特性为:同时考虑到因地突变形而引起的实际比例尺在x方向和y方向上的变形,因此
纠正后的坐标数据在不同方向上的长度比将发生变化。
注:一般的GIS软件都有仿射变换、相似变换和二次变换等几何纠正功能
(三)大地基准面
大地基准面(Geodeticdatum),设计用为最密合部份或全部大地水准面的数学模式。它由椭球体
本身及椭球体和地表上一点视为原点间之关系来定义。
此关系能以6个量来定义,通常(但非必然)是大地纬度、大地经度、原点高度、原点垂线偏差之两
分量及原点至某点的大地方位角。
每个国家或地区均有各自的基准面,我们通常称谓的北京54坐标系、西安80坐标系,指的就是
两个大地基准面。
我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐
标系。
1978年采用国际大地测量协会推荐的1975地球椭球体(IAG75)建立了我国新的大地坐标系-
-西安80坐标系。
目前大地测量基本上仍以北京54坐标系作为参照,北京54与西安80坐标之间的转换可查阅国
家测绘局公布的对照表。
WGS1984基准面采用WGS84椭球体,它是一地心坐标系,即以地心作为椭球体中心,目前GPS
测量数据多以WGS1984为基准。
椭球体与基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体
不能代表基准面,同样的椭球体能定义不同的基准面。
(四)地图投影
地图投影就是指建立地球表面(或其他星球表面或天球面)上的点与投影平面(即地图平面)上
点之间的一一对应关系的方法。
即建立地表点与投影点之间的数学转换公式。
将这个公式作为一个不可展平的曲面即地球表面投影到一个平面的基本方法,保证了空间信息
在区域上的联系与完整。
(五)高斯克吕格投影
(1)原理
这个投影是由德国数学家、物理学家、天文学家高斯于19世纪20年代拟定,后经德国大地测量
学家克吕格于1912年对投影公式加以补充,故称为高斯-克吕格投影。
高斯克吕格投影即等角横切椭圆柱投影。
假想用一个椭圆柱横切于地球椭球体的某一经线上,这条与圆柱面相切的经线,称中央经线。
以中央经线为投影的对称轴,将东西各3°或1°30′的两条子午线所夹经差6°或3°的带状地区按
数学法则、投影法则投影到圆柱面上,再展开成平面,即高斯-克吕格投影,简称高斯投影。
这个狭长的带状的经纬线网叫做高斯-克吕格投影带。
(2)基本概念:
假想有一个椭圆柱面横套在地球椭球体外面,并与某一条子午线(此子午线称为中央子午线或轴
子午线)相切,椭圆柱的中心轴通过椭球体中心,然后用一定投影方法,将中央子午线两侧各一定经差
范围内的地区投影到椭圆柱面上,再将此柱面展开即成为投影面,如图2所示,此投影为高斯投影。
高斯投影是正形投影的一种。
(3)变形规律
这种投影,将中央经线投影为直线,其长度没有变形,与球面实际长度相等。
其余经线为向极点收敛的弧线,距中央经线愈远,变形愈大。
赤道线投影后是直线,但有长度变形。
除赤道外的其余纬线,投影后为凸向赤道的曲线,并以赤道为对称轴。
经线和纬线投影后仍然保持正交。
所有长度变形的线段,其长度变形比均大于1.随远离中央经线,面积变形也愈大。
若采用分带投影的方法,可使投影边缘的变形不致过大。
我国各种大、中比例尺地形图采用了不同的高斯-克吕格投影带。
其中大于1:1万的地形图采用3°带;1:2.5万至1:5万的地形图采用6°带。
(4)分带投影n
高斯投影6度带:自0子午线起每隔经差6自西向东分带,依次编号1,2,3,…。我国6带中央子
午线的经度,由75起每隔6而至135,共计11带(13~23带)。
高斯投影3度带:它的中央子午线一部分同6带中央子午线重合,一部分同6带的分界子午线重
合。我国带共计22带(24~45带)。
(5)高斯平面直角坐标系
在投影面上,中央子午线和赤道的投影都是直线,并且以中央子午线和赤道的交点0作为坐标原
点,以中央子午线的投影为纵坐标x轴,以赤道的投影为横坐标y轴。
在我国x坐标都是正的,y坐标的最大值(在赤道上)约为330km。为了避免出现负的横坐标,可
在横坐标上加上500OOOm。此外还应在坐标前面再冠以带号。这种坐标称为国家统一坐标。例如,
有一点y=39623456.789m,该点位在39带内,位于中央子午线以东,其相对于中央子午线而言的横
坐标则是:首先去掉带号,再减去500000m,最后得=-123456.789m。
(6)高斯平面投影的特点
①中央子午线无变形;
②无角度变形,图形保持相似;
③离中央子午线越远,变形越大。
(六)UTM(横轴墨卡托)
UTM投影全称为“通用横轴墨卡托投影”UNIVERSALTRANSVERSEMERCATOLPROJECTION,
是一种“等角横轴割圆柱投影”。
椭圆柱割地球于南纬80度、北纬84度两条等高圈,投影后两条相割的经线上没有变形,而中央
经线上长度比0.9996。
该投影角度没有变形,中央经线为直线,且为投影的对称轴,中央经线的比例因子取0.9996是为
了保证离中央经线左右约330km处有两条不失真的标准经线。
该投影可以改善高斯克吕格投影,我国的卫星影像资料常用这种投影。
(七)地图投影的变形
(1)等角投影———投影前后的角度相等,但长度和面积有变形;如:高斯投影(等角横切椭圆柱投
影),墨卡托投影(等角正切椭圆柱投影),UTM投影(UniversalTransverseMercator,通用横轴等角割圆
柱墨卡托投影),兰勃特投影(等角正轴割圆锥投影)。
(2)等距投影———投影前后的长度相等,但角度和面积有变形;
(3)等积投影———投影前后的面积相等,但角度和长度有变形。如:阿尔伯特投影(正轴等面积割
圆锥投影)。
(八)矢量、栅格数据的相互转换
一、矢量向栅格的转换
重点:确定栅格单元的大小。X、Y为坐标,I、J为单元行列值。已知:Xmax、Xmin和Ymax、Ymin;
I,J。
求:栅格单元的边长ΔX和ΔY。
它们之间的关系为:
ΔY=(Xmax-Xmin)/J
ΔX=(Ymax-Ymin)/I
1.点的栅格化
设矢量坐标点(x,y),转换后的栅格单元行列值(I,J),则有:
I=1+INT[(Ymax-Y)/ΔY]
J=1+INT[(X-Xmin)/ΔX]
式中INT表示取整函数。栅格点的值用点的属性表示。
2.线的栅格化(扫描线法)
线的栅格化可以分解成对组成折线的每个线段的栅格化。
设线段的两端点坐标分别为(x1,y1)和(x2,y2)。其中有y2>y1。
根据点栅格化方法,求得两端点的行列号为(I1,J1)和(I2,J2)。△I= |I2-I1|,△J= |J2-
J1|,分两种情况处理:△I>△J和△J>△I。
当△I>△J,平行于x轴的做每一行的中心扫描线,再求每一条扫描线与线段的交点,按点的栅
格化犯法将交点转化为栅格坐标。
当△J>△I,平行于y轴的做每一列的中心扫描线,再求每一条扫描线与线段的交点,按点的栅
格化犯法将交点转化为栅格坐标。
如果△I>△J,扫描线与线段的交点的坐标为
Yi=Ymax-ΔY·(I-1/2)
Xi=[(X2-X1)/(Y2-Y1)](Yi-Y1)+X1
如果△J>△I,扫描线与线段的交点的坐标为
Xi=Xmax-ΔX·(J-1/2)
Xi=[(Y2-Y1)/(X2-X1)](Xi-X1)+Y1
例:设两个端点的行、列号已经求出,其行号为3和7,则中间网格的行号必为4、5、6。
3.面的栅格化
(1)基于弧段数据的栅格化方法
首先计算所有弧段结点或中间点所在的格网位置,并赋予该结点正确的属性,然后根据下面的算
法完成弧段的栅格化。
如图,利用弧段的数据列与格网的行列线相交,以得到正确的栅格化结果。
实际计算时,需逐段处理弧段中的局部直线段,待处理完某一局部线段后,再进行下一局部线段
的处理(如处理完AB后,再处理BC段),直至完成整条弧段的处理。局部线段与行列线求交后,存储
交点坐标,并对x或y从小到大排序。根据排序结果,相邻交点所构成线段通过的格网需赋予属
性值。
(九)矢量、栅格数据的相互转换
(2)基于多边形数据的栅格化方法
方法一:内点填充法
1.线栅格化边界(赋属性值)
2.边界内部任取一个栅格为起点
3.填充起点邻域8个栅格(赋属性值)
4.以新的邻域栅格为起点,继续重复第3步
5.直到所有的邻域栅格都填充为止
(2)基于多边形数据的栅格化方法
方法二:边界代数法
1.全部栅格都赋值0
2.线栅格化边界
3.沿着边界走,栅格向上,则左边到边界所有栅格数值减去属性值
4.沿着边界走,栅格向下,则左边到边界所有栅格数值加上属性值
方法三:对每个栅格单元,逐个判定其是否包含在某个多边形之内,并给其赋以包含该栅格单元
的多边形属性值。
“点在多边形内”的判定方法
① 检验夹角之和的方法(0或2π)2π则在其内。
② 检验交点数的方法(铅垂线法)(交点为0或偶数则在其外)
边界代数法定义:
边界代数多边形填充算法是一种基于积分思想的矢量格式向栅格格式转换算法,它适合于记录
拓扑关系的多边形矢量数据转换为栅格结构。它不是逐点判断与边界的关系完成转换,而是根据边
界的拓扑信息,通过简单的加减代数运算将边界位置信息动态地赋给各栅格点,实现了矢量格式到栅
格格式的高速转换,而不需要考虑边界与搜索轨迹之间的关系,因此算法简单、可靠性好,各边界弧段
只被搜索一次,避免了重复计算。
二、栅格向矢量的转换
把栅格单元中的空间信息转换为几何图形的过程叫矢量化。
(一)基于图像的矢量化过程如下:
(1)二值化:一般情况下,栅格数据是按0~255的不同灰度值表达的。为了简化追踪算法,需把
256个灰阶压缩为2个灰阶,即0和1两级。为此,假设任一格网的灰度值为,阀值为,那么,根据下
式就可以得到二值图。
将原来栅格单元的灰度值G(i,j)(255级)压缩为两个灰阶,即0和1两级。
(2)细化:细化是消除线划横断面栅格数的差异,使得每一条线只保留代表其轴线或周围轮廓线
(对多边形而言)位置的单个栅格的宽度。
对于栅格线划的细化方法,一般采用"剥皮法"。
剥皮法的实质是剥掉等于一个栅格宽的一层,直到最后留下彼此连通的由单个栅格点组成的
图形。
方法:① 定义阈值T, Rmin<T<Rmax;
②二值化,ifR>=T then R=1;else R=0;
(3)边界跟踪:跟踪的目的是把细化后的栅格数据整理为从结点出发的线段或闭合的线条,并以
矢量形式加以存储。
跟踪时,根据人为规定的搜索方向(如沿图幅边界的顺时针或逆时针方向),从起始点开始,在保
证趋势的情况下对八个邻域进行搜索,依次得到相邻点,最终得到完整的弧段或多边形。
(4)去除多余点及曲线光滑:由于搜索是逐个栅格进行的,所以,弧段或多边形的数据列十分
密集。
为了减少存储量,在保证线段精度的情况下可以删除部分数据点。
(5)拓扑关系的生成:判断弧段与多边形间的空间关系,以形成完整的拓扑结构并建立与属性数
据的关系。
(二)栅格数据的矢量化方法:
按行扫描,产生只记录类型边界栅格值的文件;
建立对类型边界栅格单元的追踪算法,计算边界坐标,整理成有序的坐标数组;
处理相邻类型的公共边界转换为按线段链建立的数据结构。
(十)多元空间数据的融合
一、遥感与GIS数据的融合
遥感影像与数字线化图(DLG)的融合
遥感影像与数字地形模型(DEM)的融合
遥感影像与数字栅格图(DRG)的融合
4D产品:DLG、DOM(数字正射影像图)、DEMandDRG
二、不同格式数据的融合
基于转换器的数据融合:通过交换格式进行;
基于数据标准的的数据融合:统一的数据标准;
基于公共接口的数据融合:接口是大家都遵守并达成统一的标准;
基于直接访问的数据融合:一个软件中提供了多种格式以便直接调用。
(十一)空间数据压缩
空间数据压缩:即从所取得的数据集合S中抽出一个子集A,这个子集作为一个新的信息源,在
规定的精度范围内最好地逼近原集合,而又取得尽可能大的压缩比。
比如:游程编码
(十二)空间数据内插的方法
空间数据内插:即通过已知点或分区的数据,推求任意点或分区数据的方法。
一、点的内插(3个过程)
数据取样:数据点的选取和坐标的确定;
数据内插:以数据点为控制基础,用某一数学模型来模拟地表面,进行内插计算,确定三角网或格
网结点处的特征值;
数据记录:将内插的结果,按一定的数据结构形式,存储于磁性介质中,供系统应用。
二、区域的内插:是研究根据一组分区的已知数据来推求同一地区另一组分区未知数据的内插
方法。
叠置法:将目标区与源区叠加
比重法:根据平滑密度原理将源区的数据从同性质变为非同性质
(十三)拓扑关系与拓扑结构
拓扑关系:是指网结构元素结点、弧段、面域之间的空间关系。
主要表现为拓扑邻接、拓扑关联、拓扑包含。
根据拓扑关系,不需要利用坐标或距离,可以确定一种地理实体相对于另一种地理实体的位置关
系,拓扑数据也有利于空间要素的查询。
拓扑结构:为在点、线和多边形之间建立关联,以及彻底解决邻域和岛状信息处理问题而必须建
立的数据结构。
这种结构应包括以下内容:唯一标识,多边形标识,外包多边形指针,邻接多边形指针,边界链接,
范围(最大和最小x、y坐标值)。
喜欢就点个在看再走吧