第13章-数字摄影测量基础.pptVIP

第十三章数字摄影测量根底;§13-1概述;§13-2数字影像及数字影像的获取;二、采样及量化

像片上的像点是连续分布的，但在数字化过程中不可能将每一个连续的像点全部数字化，而只能是每隔一个间隔Δ，读取一个点的灰度值，这个过程称为采样，所谓采样就是确定需要观察的数据的瞬时点。其实质是将连续的函数模型离散化。

通过上述得到每个点的灰度值一般不是整数，应将各点的灰度值取为整数，这一过程称为影像灰度的量化。所谓量化就是把像片上点上的灰度值转换成整数数字量。其方法是将航摄像片上有可能出现的最大灰度范围进行等分，分为假设干个“灰度等级”，每个点的灰度值在其相应的灰度级内，取整的原那么是四舍五入。由于计算机中数字均用二进制，因此，灰度等级一般为：

i=2M (M=1，2，……，8) (13-2)

当M=1时，i=2，只有黑、白两个灰度值；当M=8时，i=256，即有256个灰度值，其级数介于0~255之间的一个整数，0为黑，255为白，每个像元素的灰度值占8bit，即一个字节。;三、扫描坐标与像片坐标

在摄影测量中常取以像主点为坐标原点，框标连线为坐标轴的像片坐标系，由于在像片扫描的数字化过程中，其扫描坐标系一般不与像片坐标系平行且原点不同，所以同一像点的像片坐标x，y与其扫描坐标不相等，需要加以换算，这种换算称为数字影像的内定向。这项工作可借助于像片上框标的像片坐标和相应的扫描坐标来完成。

假定框标的像片坐标为(x，y)，它的扫描阵列坐标为()，用仿射变换公式进行，即

(13-3)

式中，h0、h1、h2、k0、k1、k2称为内定向参数，其数值由像片上四个框标点的扫描坐标及其相应的像平面坐标(视为理论值)组成误差方程式，平差运算求得。;四、数字影像重采样

数字影像是个规那么排列的灰度格网序列，但当对数字影像进行几何处理时，如对核线的排列、数字纠正等，由于所求得的像点不一定恰好落在原始像片上像元素的中心，要获得该像点的灰度值，就要在原采样的根底上再一次采样，即重采样。???时就必须采用适当的方法，把该点周围整数点位上灰度值对该点的灰度奉献积累起来，构成该点位新的灰度值。重采样是数字摄影测量的重要工具。常用的重采样方法有双线性插值法、双三次卷积法和最邻近像元法。;1、双线性插值法

双线性插值法的卷积核(权函数)是一个三角函数，即：

W(x)=1-|x| 0≤|x|≤1 (13-4)

此时需要待重采样点P附近的四个原始影像灰度值参加计算。如图13-1所示，11、12、21、22为相邻像元中心，像元间隔为一个单位，它们的灰度值分别为g11、g12、g21、g22，P为待重采样点的位置。计算出四个原始点为点P所作奉献的权值，以构成一个2×2的二维卷积核W(权矩阵)，把它与四个原始像元素灰度值构成的2×2灰度矩阵g作哈达玛(Hadmard)积运算，即可得到重采样点P的灰度值gP：

(13-5); W11=W(x1)W(y1) W12=W(x1)W(y2)

W21=W(x2)W(y1) W22=W(x2)W(y2)

g(i，j)*W(i，j)为两个矩阵的哈达玛积，它是这两个矩阵各对应元素的乘积之和。

W(x1)=1-Δx W(x2)=Δx

W(y1)=1-Δy W(y2)=Δy

Δx=x–int(x) Δy=y–int(y)

∴g=(1-Δx)(1-Δy)g11+(1-Δx)Δyg12+Δx(1-Δy)g21+ΔxΔyg22 (13-6);2、双三次卷积法

双三次卷只法是以三次样条函数为卷积核，其函数表达式为：;;3、最邻近像元法

最邻近像元法是取离重采样点位置最近的像元(N)的灰度值作为重采样点的灰度值，即：

gP=gN (13-9)

式中：N为最邻近点，其影像坐标值为：

xN=int(x+0.5) yN=int(y+0.5) (13-10)

以上三种方法中，最邻近像元法计算最简单、速度最快，且不破坏原始影像的灰度信息，但其几何精度较差，最大误差可达0.5个像素；双三次卷积法精度高，但计算量最大；双线性插值法既能获得较好的精度，也能到达较快的速度，是一种普遍被采用的方法。;§13-3数字影像相关;二、基于灰度的数字影像相关

所谓基于灰