地理配准与矢量配准
地理数据矫正
- 栅格数据通常是通过扫描纸质地图或采集航空及卫星照片获得 。
- 通过扫描获取的影像不包含定义其地理空间位置所需的信息。
- 航空及卫星照片所使用的坐标系统相对于通用GIS平台软件所使用的坐标系统是独立的。
- 为了能够将这些影像数据与其它的数据集成,以便进行分析, 就必须对其进行处理:用户需要事先将这些数据校准(配准)到一个指定的地图坐标系
地图数据来源主要来自于扫描仪,在扫描过程中会产生各种各样的误差,但所有这些扫描的误差引起的几何变形**,可看成平移、旋转、缩放、仿射、弯曲以及各种更高变形的综合作用结果**。在实际操作过程中,很难对这些误差一一进行变形改正,只能综合考虑它们的影响,综合校正。
矫正变换方式
输入到计算机中的图形,实际上都是通过其位置坐标(x,y)来表示,因此校正过程实质上是找一种数学关系(或函数关系),描述变换前图形坐标(x,y)与变换后图形坐标(x′,y′)之间的换算,其数学关系一般描述为
x
’
=
f
1
(
x
,
y
)
;
y
’
=
f
2
(
x
,
y
)
x’=f1(x,y); y’=f2(x,y)
x’=f1(x,y);y’=f2(x,y)
多项式变换
这个数学关系常表示为二元多项式一次、二次或三次及更高次表达式.
{
x
′
=
a
1
x
+
a
2
y
+
a
11
x
2
+
a
12
x
y
+
a
22
y
2
+
A
,
y
′
=
b
1
x
+
b
2
y
+
b
11
x
2
+
b
12
x
y
+
b
22
y
2
+
B
\begin{cases} x'=a_1x+a_2y+a_{11}x^2+a_{12}xy+a_{22}y^2+A,\\ y'=b_1x+b_2y+b_{11}x^2+b_{12}xy+b_{22}y^2+B \end{cases}
{x′=a1x+a2y+a11x2+a12xy+a22y2+A,y′=b1x+b2y+b11x2+b12xy+b22y2+B
上式是高次变换方程,在进行高次变换时,需要有6对以上控制点的坐标和理论值,才能求出待定系数。
最少控制点计算公式:
(
t
+
1
)
(
t
+
2
)
2
(t+1)(t+2)\over2
2(t+1)(t+2)
- 1次方方程最少需要3个控制点
- 2次方方程最少需要6个控制点
- 3次方方程最少需要10个控制点
当不考虑高次变换方程中的A和B时,则变成二次变换方程,称为二次变换。二次变换适用于原图有非线性变形的情况,至少需要5对控制点的坐标及其理论值,才能求出待定系数。
仿射变换
仿射变换是使用最多的一种几何纠正方式,只考虑到x和y方向上的变形,仿射变换的特性是:
- 直线变换后仍为直线;
- 平行线变换后仍为平行线;
- 不同方向上的长度比发生变化。
{ x ′ = a 1 x + a 2 y + a 3 y ′ = b 1 x + b 2 y + b 3 \begin{cases} x'=a_1x+a_2y+a_3\\ y'=b_1x+b_2y+b_3 \end{cases} {x′=a1x+a2y+a3y′=b1x+b2y+b3
对于仿射变换,只需知道不在同一直线上的三对控制点的坐标及其理论值,就可求得待定系数。但在实际使用时,往往利用4个以上的点进行纠正,利用最小二乘法处理,以提高变换的精度。
地理配准操作流程
校准栅格数据 (选择控制点)
地理配准的基本过程是在栅格图像中选取一定数据的控制点,将它们的坐标指定为矢量数据中对应点的坐标(在空间数据中,这些点的坐标是已知的,坐标系统为地图坐标系)
控制点
-
在配准中我们需要知道一些特殊点的坐标,即控制点。
-
控制点的选取:控制点可以是经纬 线网格的交点、公里网格的交点或者一些典型地物的坐标。
-
控制点的坐标:
- 如果我们知道这些点在我们矢量坐标系内坐标, 则直接输入控制点的坐标值;
- 如果不知道它们的坐标,则可以采用间接方法获取-从矢量数据中选取。
-
选取控制点:
- 控制点的数目取决于 使用哪一种数学方法来实现坐标转换.
- 但是,过多的控制点并不一定能够保证高精度的配准 。要尽可能使控制点均匀分布于整个格格图像,而不是只在图像的某个较小区域 选择控制点。
- 通常,先在图像的四个角选择4个控制点,然后在中间的位置有规律地选择一些控制点能得到较好的效果
坐标变换
- 一旦选取了足够的控制点, 就可以将栅格数据变换(或转换)到地图坐标系统下. 转换(Transformation) 运用一种数学变换方法来重新确定栅格数据中每个像元的灰度值。
- 一次多项式:仿射(affine)—变换可以将栅格数据平移, 缩放, 及旋转.
- 栅格图像中的一条直线变换后仍然为直线。矩形和正方形变换后为平行四边形
检查均方差
- 坐标转换的准确程度可以通过比较某一点在地图中的实际坐标与根据变换公式得到的坐标来判断。
- 这两个点之间的距离之差称为残差(residual error). 通过计算均方差(RMS)获取控制点总误差。
- 均方差(RMS) 的大小描述了变换公式在不同控制点间的一致性。
- 可以将残差特别大控制点删除,然后添加新的控制点。
- RMS 比较小时,说明控制点的选取是比较准确的,但也要注意有可能存在残差非常大的控制点,但由于其它控制点是很精确的,所以总的RMS误差比较小的情况。
矫正栅格数据-重采样
在地理配准过程中,将基于地图坐标生成一个“空的” 矩阵, 矩阵中每个元素的值(表示颜色)将通过重采样重新计算。
有三种通用的图像重采样技术
- 最小紧邻
- 双线性内插
- 立方卷积