数字摄影测量.pptx

《数字摄影测量.pptx》由会员分享，可在线阅读，更多相关《数字摄影测量.pptx（56页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、2023/2/121数字摄影测量数字摄影测量Digital Photogrammetry第1页/共56页2023/2/122（一）数字摄影测量概述 （二）影像数字化及影像重采样 （三）基于灰度的影像相关 （四）同名核线的确定与核线相关（五）基于特征的影像匹配（六）数字摄影测量系统（DPSDPS）数字摄影测量数字摄影测量第2页/共56页2023/2/123为能利用摄影测量方法获取地面目标的几何信息，首先必须对影像进行量测。为能利用摄影测量方法获取地面目标的几何信息，首先必须对影像进行量测。无论在模拟摄影测量阶段，还是解析摄影测量阶段，量测工作均需要人工进无论在模拟摄影测量阶段，还是解析摄影测量阶

2、段，量测工作均需要人工进行。行。（一）数字摄影测量概述 v例如：在立体坐标量测仪上进行像点坐标立体量测，或在模拟测图仪上进行定向、测绘地物与地貌时，都要人眼在立体观察情况下使左右测标对准左右同名像点。随着现代科学技术的发展，摄影测量工作者一直在研究如何用计算机代替人随着现代科学技术的发展，摄影测量工作者一直在研究如何用计算机代替人完成一些摄影测量任务，如同名像点的量测及建立立体模型等。这就是自动完成一些摄影测量任务，如同名像点的量测及建立立体模型等。这就是自动化测图的重要内容。化测图的重要内容。早期的自动化测图（早期的自动化测图（1950年美国研制了第一台）是将像片上的灰度信号转化年美国研制了

3、第一台）是将像片上的灰度信号转化为电信号，利用电子相关技术实现自动化量测。为电信号，利用电子相关技术实现自动化量测。第3页/共56页2023/2/124随着计算机技术的发展，将电信号转变成数字信号，由计算机来实现相关运随着计算机技术的发展，将电信号转变成数字信号，由计算机来实现相关运算，这种方法称为算，这种方法称为数字相关数字相关（Digital Correlation）。v20世纪60年代初美国研制的自动解析测图仪SA-11B-X及RASTER均利用数字相关技术。到了80年代，对数字相关的研究占了统治地位。利用数字灰度信号，采用数字相关技术量测同名像点，在此基础上通过解析计算，进行相对定向和

4、绝对定向，建立数字立体模型，从而建立DEM、绘制等高线、制作DOM以及为GIS提供基础信息等。这就是数字摄影测量。由于整个过程都是以数字形式在计算机中完成，因而又称为由于整个过程都是以数字形式在计算机中完成，因而又称为全数字化摄影测全数字化摄影测量量。第4页/共56页2023/2/125实现数字影像自动测图的系统称为实现数字影像自动测图的系统称为数字摄影测量系统数字摄影测量系统DPS（Digital Photogrammetric System）或或数字摄影测量工作站数字摄影测量工作站DPW（Digital Photogrammetric Workstation）。v这种系统不同于其它自动化测

5、图仪，它即没有模拟测图仪，也没有解析测图仪，它实质上是一个普通的计算机影像数据处理系统，其硬设备就是影像数字化装置、影像或图形输出装置及电子计算机，而由事先编制好的置于计算机中的软件系统来完成各种摄影测量处理工作。因此，这种系统的功能主要取决于计算机及软件系统的功能。19881988年京都第1616届ISPRSISPRS，标志着数字摄影测量在迅速发展。19921992年华盛顿第1717届ISPRSISPRS，标志着数字摄影测量正在走向实用，并步入摄影测量生产。19961996年维也纳召开的第1818届ISPRSISPRS，实现了摄影测量几何处理的自动化，并把GPSGPS技术引入摄影测量，以确定

6、摄影时的方位元素。第5页/共56页2023/2/126数字摄影测量系统数字摄影测量系统DPS之所以受到广泛重视，是因为它的功能优越于解析测之所以受到广泛重视，是因为它的功能优越于解析测图仪，更远远超过模拟测图仪。图仪，更远远超过模拟测图仪。v模拟时代的摄影测量作业大部分依赖于各种各样的摄影测量仪器。如坐标量测仪、纠正仪、正射投影仪等等。每种仪器均各具不同的用途，其设计、结构、通用性、灵活性、精度、输入输出等亦各异。v解析测图仪则前进了一大步：在一台解析测图仪上便可完成多种不同类型的任务。v而数字摄影测量系统则又在解析测图仪的基础上前进了一大步：集各种功能于一体，且其应用范围也将大大扩展。第6页

7、/共56页2023/2/127数字摄影测量工作站数字摄影测量工作站DPW除了能胜任解析测图仪的一切任务外，尚具有许多除了能胜任解析测图仪的一切任务外，尚具有许多新的功能，如：新的功能，如：影像位移的去除、任意方式的纠正、反差的扩展、多幅影像的比较分析、图像识别、影像数字相关、以及数据库的管理等等；通过显示器还可观察数字图像以及框标、控制点、连接点、DEM及其它所需特征；易于实现用于整体检查和质量控制的图形显示或叠合，甚至进行立体显示；可对图像自动进行所需要的特征提取；具有实时数据获取和处理的能力。等等。第7页/共56页2023/2/128目前数字摄影测量系统仍处于发展的时期，其自动化功能仅限于

8、几何处理，目前数字摄影测量系统仍处于发展的时期，其自动化功能仅限于几何处理，即可进行自动内定向、相对定向，自动建立即可进行自动内定向、相对定向，自动建立DEM、制作、制作DOM。其它工作要。其它工作要作采用自动或人工的方式进行。作采用自动或人工的方式进行。特别是地物的测绘，目前全部是人工交互方式。特别是地物的测绘，目前全部是人工交互方式。数字摄影测量系统数字摄影测量系统第8页/共56页2023/2/129数字影像的产生：直接从空间飞行器中的直接从空间飞行器中的扫描式传感器扫描式传感器产生；产生；利用利用影像数字化器影像数字化器对摄取的像片通过影像数字化过程获得。对摄取的像片通过影像数字化过程获

9、得。（二）影像数字化与影像重采样 Vexcel Ultra5000专业航片扫描仪第9页/共56页2023/2/1210一、影像数字化仪一、影像数字化仪一、影像数字化仪一、影像数字化仪有多种，主要有：有多种，主要有：电子光学扫描器电子光学扫描器和和固体阵列式数字化器固体阵列式数字化器等。等。电子光学扫描器电子光学扫描器主要分：滚筒式和平台式两类主要分：滚筒式和平台式两类。xy光源第10页/共56页2023/2/1211一、影像数字化仪一、影像数字化仪一、影像数字化仪一、影像数字化仪在一条线上可以排列到在一条线上可以排列到2048个传感器。亦可由多组个传感器。亦可由多组2048个传感个传感器串联成

10、更长的线阵列。在一个面积上可以排列成面阵列式传器串联成更长的线阵列。在一个面积上可以排列成面阵列式传感器。感器。固体阵列式数字化器固体阵列式数字化器是使用在一条线上或者是在一个面积上排列的半导体传感器（电荷耦合装置CCD），对像片进行数字化，而无需使用扫描头的移动。第11页/共56页2023/2/1212二、影像数字化过程二、影像数字化过程二、影像数字化过程二、影像数字化过程影像数字化影像数字化：将透明正片（或负片）放在影像数字化器上，把像片上将透明正片（或负片）放在影像数字化器上，把像片上像点的灰度值用数字形式记录下来的过程。像点的灰度值用数字形式记录下来的过程。Vexcel Ultra50

11、00专业航片扫描仪影像的灰度影像的灰度又称为光学密度。透明像片（正片或负片）上影像的灰度值，反映了又称为光学密度。透明像片（正片或负片）上影像的灰度值，反映了它透明的程度，即透光的能力。它透明的程度，即透光的能力。第12页/共56页2023/2/1213二、影像数字化过程二、影像数字化过程二、影像数字化过程二、影像数字化过程影像的灰度影像的灰度设投影在透明像片上的光通量为设投影在透明像片上的光通量为F0，而透过透明像片后的光通，而透过透明像片后的光通量为量为F，则，则F与与F0之比称为透过率之比称为透过率T，F0与与F之比称为不透过率之比称为不透过率O：因此，像点愈黑，则透过的光通量愈小，不透

12、过率愈大，所以，因此，像点愈黑，则透过的光通量愈小，不透过率愈大，所以，透过率和不透过率都可以说明影像黑白的程度。透过率和不透过率都可以说明影像黑白的程度。第13页/共56页2023/2/1214二、影像数字化过程二、影像数字化过程二、影像数字化过程二、影像数字化过程影像的灰度影像的灰度但是，人眼对明暗程度的感觉是按对数关系变化的。但是，人眼对明暗程度的感觉是按对数关系变化的。为了适应为了适应人眼的视觉，在分析影像的性能时，不直接用透过率或不透过人眼的视觉，在分析影像的性能时，不直接用透过率或不透过率表示其黑白程度，而用不透过率的对数值表示：率表示其黑白程度，而用不透过率的对数值表示：D称为影

13、像的灰度称为影像的灰度，当光线全部透过时，即透过率为，当光线全部透过时，即透过率为1，则影像，则影像的灰度为的灰度为0；当光通量仅透过百分之一，即不透过率是；当光通量仅透过百分之一，即不透过率是100时，时，则影像的灰度为则影像的灰度为2，实际的航空底片的灰度一般在至范围之内。，实际的航空底片的灰度一般在至范围之内。第14页/共56页2023/2/1215二、影像数字化过程二、影像数字化过程二、影像数字化过程二、影像数字化过程采样：采样：在像片上像点是连续分布的，但在影像数字化过程中不在像片上像点是连续分布的，但在影像数字化过程中不可能将每一个连续的像点全部数字化，而只能每隔一个间隔可能将每一

14、个连续的像点全部数字化，而只能每隔一个间隔（）读一个点的灰度值，这个过程称为采样，）读一个点的灰度值，这个过程称为采样，称为采样间称为采样间隔。隔。采样与量化采样与量化影像数字化过程包括影像数字化过程包括采样采样与与量化量化两项内容两项内容。采样后是不连续的等间隔序列，采样过程会给影像的灰度采样后是不连续的等间隔序列，采样过程会给影像的灰度带来误差。所以采样间隔越小越好。但是采样间隔越小，带来误差。所以采样间隔越小越好。但是采样间隔越小，数据量越大。数据量越大。第15页/共56页2023/2/1216二、影像数字化过程二、影像数字化过程二、影像数字化过程二、影像数字化过程量化：量化：通过采样过

15、程得到的每个点的灰度值不是整数，不方便通过采样过程得到的每个点的灰度值不是整数，不方便计算。为此，应将各点的灰度值取为整数，这个一过程称为计算。为此，应将各点的灰度值取为整数，这个一过程称为影影像灰度的量化像灰度的量化。采样与量化采样与量化影像数字化过程包括影像数字化过程包括采样采样与与量化量化两项内容两项内容。其方法是将透明像片有可能出现的最大灰度变化范围进行等分，其方法是将透明像片有可能出现的最大灰度变化范围进行等分，等分的数目称为等分的数目称为“灰度等级灰度等级”，然后将每个点的灰度值在其相应的灰度等级内取整，取整的原然后将每个点的灰度值在其相应的灰度等级内取整，取整的原则是四舍五入。则

16、是四舍五入。第16页/共56页2023/2/1217二、影像数字化过程二、影像数字化过程二、影像数字化过程二、影像数字化过程量化：量化：采样与量化采样与量化影像数字化过程包括影像数字化过程包括采样采样与与量化量化两项内容两项内容。由于数字计算机中数字均用二进制表示，因此灰度等级一般都由于数字计算机中数字均用二进制表示，因此灰度等级一般都取为取为2m（m是正整数）。是正整数）。当当m=1时，灰度只有黑白两级，当时，灰度只有黑白两级，当m=8时，则得时，则得256个灰度级。个灰度级。其中，级数是介于其中，级数是介于0255之间，之间，0为黑，为黑，255为白，每个像元素为白，每个像元素的长度值占的

17、长度值占8bit，即一个字节。，即一个字节。第17页/共56页2023/2/1218三、数字影像的构成三、数字影像的构成三、数字影像的构成三、数字影像的构成经过影像数字化以后得到的数字影像是一个二维的数字矩阵：经过影像数字化以后得到的数字影像是一个二维的数字矩阵：矩阵的每个元素称一个元素矩阵的每个元素称一个元素(Pixel=Picture element)。对各像元素所。对各像元素所赋予的灰度值赋予的灰度值g(m,n)一般在一般在0255之间某个整数，代表像元素的黑之间某个整数，代表像元素的黑白程度。矩阵的每一行对应于一个扫描行。白程度。矩阵的每一行对应于一个扫描行。第18页/共56页2023

18、/2/1219三、数字影像的构成三、数字影像的构成三、数字影像的构成三、数字影像的构成第19页/共56页2023/2/1220四、数字影像内定向四、数字影像内定向四、数字影像内定向四、数字影像内定向在摄影测量中常取用在摄影测量中常取用以像主点为原点的以像主点为原点的像平面坐标像平面坐标来建立像点与地面来建立像点与地面点的坐标关系。点的坐标关系。数字影像内定向：由于在像片扫描的数字化过程中，像片的扫描坐标系与像平面由于在像片扫描的数字化过程中，像片的扫描坐标系与像平面坐标系坐标系不平行不平行，且，且坐标原点不同坐标原点不同，所以同一像点的像平面坐标，所以同一像点的像平面坐标x，y与其扫描与其扫描

19、坐标坐标x，y不相等不相等，需要加以换算，这种换算称为，需要加以换算，这种换算称为数字影像内定向数字影像内定向。第20页/共56页2023/2/1221四、数字影像内定向四、数字影像内定向四、数字影像内定向四、数字影像内定向 其数值由像片上四个框标的扫描坐标及其相应的像平面坐标（视为理论值）组成误差方程式，用平差运算求得。第21页/共56页2023/2/1222立体测图的关键是要已知同名像点在左、右像片上的位置，无论是在模拟测立体测图的关键是要已知同名像点在左、右像片上的位置，无论是在模拟测图仪上还是在解析仪上作业，都需要作业员通过人眼的立体观测，不断地从图仪上还是在解析仪上作业，都需要作业员

20、通过人眼的立体观测，不断地从左、右像片上搜索同名像点，也是就探求影像的相关。左、右像片上搜索同名像点，也是就探求影像的相关。对于全数字化摄影测量，在没有人眼的立体观测的情况下，如何从左、右数对于全数字化摄影测量，在没有人眼的立体观测的情况下，如何从左、右数字影像中寻找同名像点，亦即数字影像相关，是其核心问题。字影像中寻找同名像点，亦即数字影像相关，是其核心问题。对这个问题的研究，从分析影像的灰度特性入手，提出了许多各具特色的数对这个问题的研究，从分析影像的灰度特性入手，提出了许多各具特色的数字影像相关方法，如：协方差法、相关系数法、高精度最小二乘相关等等。字影像相关方法，如：协方差法、相关系数

21、法、高精度最小二乘相关等等。这些方法有一个共同的特点，即：它们都是基于待相关点所在的一个小区域这些方法有一个共同的特点，即：它们都是基于待相关点所在的一个小区域内的影像灰度。内的影像灰度。（三）基于灰度的影像相关 第22页/共56页2023/2/1223一、相关系数法一、相关系数法一、相关系数法一、相关系数法假设在假设在左片上左片上有一个有一个目标点目标点，为了搜索它在，为了搜索它在右片上的同名点右片上的同名点，须以，须以这这目标点目标点为中心取其周围为中心取其周围nn个像元素的灰度序列组成一个个像元素的灰度序列组成一个目标区目标区，如，如图：图：在在目标区目标区中任意一个像元素的灰度值设为中

22、任意一个像元素的灰度值设为gij（i,j1,2,n），一般），一般取取n为奇数，其为奇数，其中心点即为目标点中心点即为目标点。根据根据左片上目标点的坐左片上目标点的坐标标概略地估计出它在概略地估计出它在右右片上的近似点位片上的近似点位，并以，并以此为中心取其周围此为中心取其周围lm个个影像灰度序列（影像灰度序列（l,mn），），组成一个组成一个搜索区搜索区。第23页/共56页2023/2/1224一、相关系数法一、相关系数法一、相关系数法一、相关系数法在搜索区内有（在搜索区内有（l-n+1）（m-n1）个与目标区等大的区域，称为）个与目标区等大的区域，称为相关窗口相关窗口，窗口内任意一点的灰度

23、值设为，窗口内任意一点的灰度值设为gi+k,j+h（k=0,1,l-n，h=0,1,m-n）。）。为了在右像片上的搜索区内寻找同名点，须按下式计算相关系数值：为了在右像片上的搜索区内寻找同名点，须按下式计算相关系数值：第24页/共56页2023/2/1225上式可计算出（l-n+1）（m-n1）个相关系数。当取最大值时，对应的相关窗口的中点即被认为是待定点的同名像点。第25页/共56页2023/2/1226二、协方差法二、协方差法二、协方差法二、协方差法用协方差法搜索同名像点的过程与相关系数法基本相同，不同之处仅用协方差法搜索同名像点的过程与相关系数法基本相同，不同之处仅在于采用的相似性判据不

24、同。在于采用的相似性判据不同。这里采用协方差值作为相似性判据。为了搜索同名像点须按下式计算这里采用协方差值作为相似性判据。为了搜索同名像点须按下式计算协方差值：协方差值：这样，共有（l-n+1）（m-n1）个协方差值，取其最小者对应的相关窗口的中心，即为目标点的同名像点。第26页/共56页2023/2/1227略略三、高精度最小二乘相关三、高精度最小二乘相关三、高精度最小二乘相关三、高精度最小二乘相关第27页/共56页2023/2/1228在影像相关方法中，无论是目标区，还是搜索区，都是一个二维的影像窗口，在影像相关方法中，无论是目标区，还是搜索区，都是一个二维的影像窗口，在进行相关计算时，其

25、计算量是相当大的。在进行相关计算时，其计算量是相当大的。由摄影测量基本知识可知，通过摄影基线所作的任意一个与像片面相交的平由摄影测量基本知识可知，通过摄影基线所作的任意一个与像片面相交的平面，与像片对相交，就会在左、右像片上获得一对面，与像片对相交，就会在左、右像片上获得一对同名核线同名核线，由核线的几何，由核线的几何关系确定了同名点必然位于同名核线上。关系确定了同名点必然位于同名核线上。这样，利用核线这样，利用核线(epipolar line)的概念就能将的概念就能将沿着沿着x，y方向搜索同名点的二方向搜索同名点的二维相关问题维相关问题，改成为，改成为沿同名核线的一维相关问题沿同名核线的一维

26、相关问题。（把二维变到一维上。）（把二维变到一维上。）（四）同名核线的确定与核线相关 第28页/共56页2023/2/1229但是在影像数字化过程中，像元素按矩阵形式规则排列，扫描行不是核线方但是在影像数字化过程中，像元素按矩阵形式规则排列，扫描行不是核线方向。向。因此，欲进行核线相关，必须先找到核线，而且核线确定的精度直接影响到因此，欲进行核线相关，必须先找到核线，而且核线确定的精度直接影响到影像相关的精度。影像相关的精度。当两个比例尺不同的影像相关时，可能同名点的灰度相关很大，灰度相关就当两个比例尺不同的影像相关时，可能同名点的灰度相关很大，灰度相关就很难实现，所以只能基于特征相关。很难实

27、现，所以只能基于特征相关。第29页/共56页2023/2/1230一、核线几何关系解析一、核线几何关系解析一、核线几何关系解析一、核线几何关系解析确定同名核线的方法很多，原理上最简单的一种方法是：确定同名核线的方法很多，原理上最简单的一种方法是：基于数字影基于数字影像几何纠正的方法像几何纠正的方法。同名核线在航摄像片上是互相不平行的，它们交于一个同名核线在航摄像片上是互相不平行的，它们交于一个点（点（核点核点）。但是，如果将像片上的核线但是，如果将像片上的核线投影到一对投影到一对“相对水平相对水平”的像片的像片（平行于平行于投影基线的像片对投影基线的像片对），则），则核线相互平行核线相互平行。

28、如下图：。如下图：（b）水平像片（b）实际像片核线的几何关系第30页/共56页2023/2/1231正是由于正是由于相对水平像片对相对水平像片对具有这一特性，我们就有可能在相对水平像片上建具有这一特性，我们就有可能在相对水平像片上建立规则的格网，它的行就是立规则的格网，它的行就是核线核线，核线上像元素（坐标为核线上像元素（坐标为xt，yt）的灰度可由它对应的实际像片上的像元素（坐）的灰度可由它对应的实际像片上的像元素（坐标为标为x，y）的灰度求得，即）的灰度求得，即g(xt，yt)g(x，y)。下右图代表通过摄影基线下右图代表通过摄影基线SS=B和某一个摄影射线和某一个摄影射线SA所构成的平面

29、，亦即通所构成的平面，亦即通过像点过像点a的核平面，的核平面，图中图中P代表位在左方的航摄像片代表位在左方的航摄像片。第31页/共56页2023/2/1232t代表相应的平行于摄影基线的水平像片，代表相应的平行于摄影基线的水平像片，at为为A点在点在左水平像片上的相应像左水平像片上的相应像点点。设设a，at在各自的像片坐标系中的坐标分别为（在各自的像片坐标系中的坐标分别为（x，y）和（）和（xt，yt），则：），则：第32页/共56页2023/2/1233式中式中a1，a2，c3为左片的方向余弦，为这张像片相对于摄影基线的角为左片的方向余弦，为这张像片相对于摄影基线的角方位元素的函数，方位元素

30、的函数，f 为像片主距。为像片主距。显然在相对水平像显然在相对水平像片上，当片上，当yt常数，常数，则为则为核线核线。将。将yt=c代入上式，经整理代入上式，经整理得：得：第33页/共56页2023/2/1234若以等间隔取一系列的若以等间隔取一系列的xt值值k，(k+1)，(k+2)，即求得一系列的像，即求得一系列的像片坐标片坐标(x1，y1),(x2，y2),(x3，y3)，这些像点就位于倾斜像片，这些像点就位于倾斜像片P的核线的核线上，将这些像点的灰度上，将这些像点的灰度 g(x1，y1)，g(x2，y2)，直接赋给相对水平像片直接赋给相对水平像片上相应的像点，即：上相应的像点，即：g0

31、(k,c)=g(x1，y1)，g0(k+1),c)=g(x2，y2)就能获得相对水平像片上的核线影像。就能获得相对水平像片上的核线影像。第34页/共56页2023/2/1235由于在相对水平像片上，同名核线的由于在相对水平像片上，同名核线的yt坐标值相等，因此，将同样的坐标值相等，因此，将同样的yt=c代代入右片共线方程：入右片共线方程：式中式中a1，a2，c3为右片的方向余弦，分别是右片相对于摄影基线的为右片的方向余弦，分别是右片相对于摄影基线的角方位元素的函数，角方位元素的函数，f 为像片主距。为像片主距。由此即得右片上的同名核线。由此即得右片上的同名核线。第35页/共56页2023/2/

32、1236数字影像是个规则排列的灰度格网序列，为了获取核线的灰度序列，数字影像是个规则排列的灰度格网序列，为了获取核线的灰度序列，必须对原始影像灰度重采样。必须对原始影像灰度重采样。二、沿核线重采样二、沿核线重采样二、沿核线重采样二、沿核线重采样右图为相对水平像片上沿核线排右图为相对水平像片上沿核线排列的规则格网，为求得每个格网列的规则格网，为求得每个格网的灰度，须按下式依次将每个格的灰度，须按下式依次将每个格网点的坐标（网点的坐标（xt，yt）反算到原始）反算到原始像片上，得到相应的坐标（像片上，得到相应的坐标（x，y）。）。第36页/共56页2023/2/1237但是，由于所求得的像点不一定

33、恰好落在原始像片上像元素的中心，因此，但是，由于所求得的像点不一定恰好落在原始像片上像元素的中心，因此，就必须进行灰度重采样。就必须进行灰度重采样。如下右图，如下右图，为原始像片数字化时的采样间隔，待求的像元素的灰度为原始像片数字化时的采样间隔，待求的像元素的灰度g可由其可由其周围四个像元的灰度值周围四个像元的灰度值g1，g2，g3，g4经双线性内插求得：经双线性内插求得：第37页/共56页2023/2/1238与二维相关不同的是，其目标区和搜索区分别位于左、右同名核线上，与二维相关不同的是，其目标区和搜索区分别位于左、右同名核线上，均为一维的影像窗口，现以相关系数法为例，介绍核线相关的过程。

34、均为一维的影像窗口，现以相关系数法为例，介绍核线相关的过程。为了沿同名核线搜索同名点，在左核线上建立一个目标区，该目标区为了沿同名核线搜索同名点，在左核线上建立一个目标区，该目标区中心就是目标点，目标区的长度为中心就是目标点，目标区的长度为n个像元素（个像元素（n为奇数）；为奇数）；在右片上沿同名核线建立搜索区，其长度为在右片上沿同名核线建立搜索区，其长度为m个像元素，如下图：个像元素，如下图：三、核线相关三、核线相关三、核线相关三、核线相关第38页/共56页2023/2/1239为找同名点，须按下式计算相关系数值：为找同名点，须按下式计算相关系数值：共计算（共计算（m-n+1）个相关系数，判

35、别其中最大的一个，设）个相关系数，判别其中最大的一个，设kk0时相时相关系数取得最大值，则同名像点在搜索区内的序号为关系数取得最大值，则同名像点在搜索区内的序号为k0+(n+1)/2。第39页/共56页2023/2/1240前面介绍了几种基于灰度的影像相关方法，它们均是直接以待定点为中心的前面介绍了几种基于灰度的影像相关方法，它们均是直接以待定点为中心的窗口内影像灰度值为依据进行同名点的搜索。窗口内影像灰度值为依据进行同名点的搜索。一般而言，这些方法大多不能顾及图像的总体结构，而是机械地按照某种或一般而言，这些方法大多不能顾及图像的总体结构，而是机械地按照某种或几种相似性判据逐像元以一定大小的

36、窗口顺序进行相关搜索。几种相似性判据逐像元以一定大小的窗口顺序进行相关搜索。虽然在某些情况下也能获得较高的相关精度，但是对于信息贫乏区域，或相虽然在某些情况下也能获得较高的相关精度，但是对于信息贫乏区域，或相关影像之间存在着较大比例尺差异或扭曲的区域，相关则难免失败。关影像之间存在着较大比例尺差异或扭曲的区域，相关则难免失败。（五）基于特征的影像匹配 第40页/共56页2023/2/1241而而基于特征的影像匹配（基于特征的影像匹配（feature based image matching）是另一类搜索同是另一类搜索同名像点的方法。名像点的方法。基本思想基本思想是：是：首先用某种特征提取算子提

37、取影像中的特征（点、线、面）；然后对提取的特征进行参数描述；最后以特征的参数值为依据进行同名特征的搜索，继而获得同名像点。第41页/共56页2023/2/1242一、特征提取一、特征提取一、特征提取一、特征提取用于影像匹配的特征分为用于影像匹配的特征分为点特征点特征和和线特征线特征两种。在数字图像处理中通两种。在数字图像处理中通常将常将线特征线特征称为称为边缘边缘。用于用于提取点特征的算子提取点特征的算子称为称为有利算子有利算子或或兴趣算子兴趣算子，提取的特征点提取的特征点称为称为有利点有利点或或兴趣点兴趣点；用于用于提取边缘的算子提取边缘的算子称为称为边缘检测算子边缘检测算子或或边缘检测器边

38、缘检测器。一般而言，点特征（如明显地物点）具有较高的匹配精度。一般而言，点特征（如明显地物点）具有较高的匹配精度。特别是当图像的方位元素未知时，往往需要首先匹配少量点求解特别是当图像的方位元素未知时，往往需要首先匹配少量点求解图像的相对方位元素，这时点特征就显示其重要性。图像的相对方位元素，这时点特征就显示其重要性。但孤立的特征点匹配亦可能出现多义性。但孤立的特征点匹配亦可能出现多义性。第42页/共56页2023/2/1243一、特征提取一、特征提取一、特征提取一、特征提取如果已知图像的方位元素，则可以方便地引用核线约束条件。如果已知图像的方位元素，则可以方便地引用核线约束条件。这时与核线正交

39、或近似正交的边缘特征不仅可以达到高的匹配精这时与核线正交或近似正交的边缘特征不仅可以达到高的匹配精度，同时还可用作核线与核线之间匹配信息的传递媒介。度，同时还可用作核线与核线之间匹配信息的传递媒介。即：将已匹配成功的核线上的边缘信息沿着边缘传递到待匹配的即：将已匹配成功的核线上的边缘信息沿着边缘传递到待匹配的核线，从而提高核线匹配的精度、速度和可靠性。核线，从而提高核线匹配的精度、速度和可靠性。第43页/共56页2023/2/12441 1 1 1、用于提取点特征的有利算子、用于提取点特征的有利算子、用于提取点特征的有利算子、用于提取点特征的有利算子Moravec算子算子是通过逐像元量测与其邻

40、元的灰度差，搜索相邻像元之间具有高反差是通过逐像元量测与其邻元的灰度差，搜索相邻像元之间具有高反差的点，的点，具体方法有：具体方法有：计算各像元的有利值计算各像元的有利值IV（Interest Value）。）。给一定阈值，确定特定的有利点。给一定阈值，确定特定的有利点。抑制局部非最大。抑制局部非最大。第44页/共56页2023/2/1245Forstner算子算子是通过逐像元计算其是通过逐像元计算其Roberts梯度和协方差矩阵寻找具有尽可能小而梯度和协方差矩阵寻找具有尽可能小而圆的误差随圆的点作为特征点。圆的误差随圆的点作为特征点。具体方法有：具体方法有：以以55或更大的窗口逐像元计算或更

41、大的窗口逐像元计算Roberts梯度和协方差矩阵，以梯度和协方差矩阵，以及有利值及有利值q和和w。确定待选的有利窗口。确定待选的有利窗口。抑制局部非最大。抑制局部非最大。在最佳窗口中确定加权中心作为最后所需的有利点，即特征点。在最佳窗口中确定加权中心作为最后所需的有利点，即特征点。第45页/共56页2023/2/1246LY算子算子。具体方法略。具体方法略。第46页/共56页2023/2/12472 2 2 2、边缘检测算子、边缘检测算子、边缘检测算子、边缘检测算子直方图法直方图法，具体方法略。具体方法略。第47页/共56页2023/2/1248二、基于特征的影像匹配二、基于特征的影像匹配二、

42、基于特征的影像匹配二、基于特征的影像匹配由于提取的特征有由于提取的特征有点特征点特征和和边缘特征边缘特征之分，而且对于边缘特征，又由之分，而且对于边缘特征，又由于描述特征的参数不同，继而在匹配过程中采用的相似性基元也不同，于描述特征的参数不同，继而在匹配过程中采用的相似性基元也不同，从而形成了许多具有不同特点的特征匹配算法。从而形成了许多具有不同特点的特征匹配算法。Barnard-Thompson算法算法Greenfeld-Schenk算法算法跨接法影像匹配方法跨接法影像匹配方法第48页/共56页2023/2/1249数字摄影测量系统的任务是基于数字影像或数字化影像完成摄影测量作业。数字摄影测

43、量系统的任务是基于数字影像或数字化影像完成摄影测量作业。原则上，数字摄影测量系统是对影像进行自动量测与识别的系统。原则上，数字摄影测量系统是对影像进行自动量测与识别的系统。（六）数字摄影测量系统 JX4-4AVirtuoZo NT第49页/共56页2023/2/1250一、主要功能一、主要功能一、主要功能一、主要功能AT：空中三角测量；：空中三角测量；DAT：半自动数字化空中三角测量；：半自动数字化空中三角测量；DC：数字制图；：数字制图；DM数字测图：特征提取与更新、属性化及绘图；数字测图：特征提取与更新、属性化及绘图；DTM：人工或自动化获取数字高程模型，包括影像相关功能；：人工或自动化获

44、取数字高程模型，包括影像相关功能；GIS：GIS功能，通常用于栅格数据的分析、叠置和表达；功能，通常用于栅格数据的分析、叠置和表达；GO：几何操作，坐标量测，影像配准；：几何操作，坐标量测，影像配准；第50页/共56页2023/2/1251GPS：GPS辅助空中三角测量；辅助空中三角测量；IP：影像处理；：影像处理；OM：正射影像图的生成；：正射影像图的生成；OR：内定向和外定向；：内定向和外定向；RS：遥感功能，主要是多光谱影像分类；：遥感功能，主要是多光谱影像分类；SCA：像片扫描数字化；：像片扫描数字化；TIN：不规则三角网；：不规则三角网；VI：可视化，特别是地形三维显示；：可视化，特

45、别是地形三维显示；一、主要功能一、主要功能一、主要功能一、主要功能第51页/共56页2023/2/1252二、硬件二、硬件二、硬件二、硬件主要由两部分组成：主要由两部分组成：一是数字影像获取装置及成果输出设备；一是数字影像获取装置及成果输出设备；另（也是核心部分）是计算机及其它外部设备；另（也是核心部分）是计算机及其它外部设备；第52页/共56页2023/2/1253三、软件三、软件三、软件三、软件包括操作系统和应用软件两大部分：包括操作系统和应用软件两大部分：操作系统：操作系统：DOS，UNIX，WINDOWS，NT等；等；应用软件：实际上是应用软件：实际上是解析摄影测量软件解析摄影测量软件

46、和和数字图形与图像处理软数字图形与图像处理软件件的集合。的集合。如果：如果：定向参数计算：定向参数计算：内定向；外定向，包括相对定向和绝对定向；内定向；外定向，包括相对定向和绝对定向；空中三角测量空中三角测量形成按核线方向排列的立体影像形成按核线方向排列的立体影像影像相关或特征匹配影像相关或特征匹配第53页/共56页2023/2/1254建立数字高程模型建立数字高程模型自动绘制等高线自动绘制等高线制作正射影像制作正射影像等高线与正射影像叠加，制作带等高线的正射影像图等高线与正射影像叠加，制作带等高线的正射影像图制作透视图和景观图制作透视图和景观图数字影像的机助量测，如地物、地貌元素的量测数字影像的机助量测，如地物、地貌元素的量测地图编辑与注记地图编辑与注记第54页/共56页2023/2/1255四、几种数字摄影测量系统简介四、几种数字摄影测量系统简介四、几种数字摄影测量系统简介四、几种数字摄影测量系统简介VirtuoZo NT：原武汉测绘科技大学研制：原武汉测绘科技大学研制JX-4A(DPW)：中国测绘科学研究院研制：中国测绘科学研究院研制JX4-4AVirtuoZo NT第55页/共56页2023/2/1256感谢您的观看！第56页/共56页