《一种真正射影像的生成方法.pdf》由会员分享,可在线阅读,更多相关《一种真正射影像的生成方法.pdf(21页珍藏版)》请在专利查询网上搜索。
1、10申请公布号CN104200527A43申请公布日20141210CN104200527A21申请号201410443508922申请日20140902G06T17/05201101G06T7/00200601G01C11/0020060171申请人西安煤航信息产业有限公司地址710054陕西省西安市长胜街78号72发明人刘敏郭永春梁菲姚春雨左涛74专利代理机构西安创知专利事务所61213代理人谭文琰54发明名称一种真正射影像的生成方法57摘要本发明公开了一种真正射影像的生成方法,包括以下步骤一、获取被测量区域的原始测量数据;二、遮蔽检测,过程如下201、遮蔽检测范围确定;202、遮蔽检测对。
2、当前所检测纠正图像上的所有待检测点分别进行遮蔽检测;对任一个待检测点遮蔽检测时,过程如下自适应搜索路径方位确定、自适应搜索路径长度确定、自适应搜索路径上检查点数量及各检查点位置确定和遮蔽判断;三、影像纠正根据步骤二中的遮蔽检测结果对当前所检测纠正图像进行纠正。本发明方法步骤简单、设计合理且生产效率高、使用效果好,能解决现有真正射影像制作方法存在的效率低、实用性差等问题,以解决在制作真正射影像中的精度和速度问题,提高生产效率。51INTCL权利要求书5页说明书14页附图1页19中华人民共和国国家知识产权局12发明专利申请权利要求书5页说明书14页附图1页10申请公布号CN104200527ACN。
3、104200527A1/5页21一种真正射影像的生成方法,其特征在于该方法包括以下步骤步骤一、获取被测量区域的原始测量数据对被测量区域进行航空摄影测量,摄取被测量区域的航空摄影图像,并获得被测量区域的原始测量数据;所述原始测量数据为DSM测量数据或DBM测量数据,所述DSM测量数据为数字表面模型的测量数据,所述DBM测量数据为数字建筑物模型的测量数据;步骤二、遮蔽检测根据步骤一中所述原始测量数据,先采用数据处理器获取步骤一中所述航空摄影图像的纠正图像,所述纠正图像为根据当前所处理航空摄影图像在地面上的投影范围和需生成真正射影像的地面分辨率进行划分所形成的格网图像,所述格网图像的大小与需生成真正。
4、射影像的大小相同且所述格网图像上的各格网点分别与所述真正射影像上的各像素点一一对应,所述真正射影像为对所述航空摄影图像进行数字微分纠正和重采样后获得的影像;之后,再对所述纠正图像进行遮蔽检测,过程如下步骤201、遮蔽检测范围确定先根据当前所处理航空摄影图像的外方位元素和内方位元素,计算得出当前所处理航空摄影图像的四个角点的地面坐标;之后,根据计算得出的四个角点的地面坐标,确定当前所检测纠正图像的遮蔽检测范围;当前所处理航空摄影图像的四个角点分别为点A、点B、点C和点D,点A、点B、点C和点D的地面坐标分别为X1,Y1、X2,Y2、X3,Y3和X4,Y4;并且,根据公式计算得出四个角点的地面坐标。
5、;公式1中,I为正整数且I1、2、3、4;XS,YS,ZS为当前所处理航空摄影图像的所述外方位元素中的摄影中心点坐标,F为所述内方位元素中的一个参数且其为步骤一中对被测量区域进行航空摄影测量时所用航摄仪的焦距;XI,YI为当前所处理航空摄影图像上四个角点的二维平面坐标;Z为被测量区域的地面平均高度;计算得出当前所处理航空摄影图像上四个角点的地面坐标后,根据公式计算出所述纠正图像的大小;公式9中GRIDSIZE为预先设定的所述真正射影像的地面分辨率,所述真正射影像的大小为MN个像素点,所述纠正图像上包括MN个格网点;其中,M和N分别为所述纠正图像中格网点的列数和行数;所述遮蔽检测范围的横坐标为X。
6、MINXMAX且其纵坐标为YMINYMAX,其中XMIN和XMAX分别为四个角点的地面坐标中横坐标的最小值和最大值,YMIN和YMAX分别为四个角点的地面坐标中纵坐标的最小值和最大值;当前所检测纠正图像的旋转矩阵权利要求书CN104200527A2/5页3步骤202、遮蔽检测先将当前所检测纠正图像划分为多个遮蔽检测面元,每个所述检测面元中均包括N1N1个格网点,其中N1为正整数,且N1N11155;将当前所检测纠正图像上的各检测面元均作为一个待检测点分别进行遮蔽检测,并且各待检测点的遮蔽检测方法均相同;其中,对当前所检测纠正图像上的任一个待检测点进行遮蔽检测时,过程如下步骤2021、自适应搜索。
7、路径方位确定根据公式对当前检测点的自适应搜索路径的方位进行确定;其中,当前检测点为当前状态下进行遮蔽检测的待检测点;公式2中,为当前检测点相对于投影中心的方位角,所述投影中心为所述航空摄影图像的摄影中心点在地面上的投影点;当N11时,XP,YP为当前检测点的二维平面坐标;当N125时,XP,YP为当前检测点中所包括N1个格网点的平均二维平面坐标;所述自适应搜索路径为对当前检测点进行遮蔽检测时所采用的搜索路径;步骤2022、自适应搜索路径长度确定先根据公式计算得出当前检测点到投影中心的水平距离D;再根据公式计算得出当前检测点的自适应搜索路径的长度L;公式4中,L为当前检测点到投影中心的搜索长度,。
8、HMAX为当前所检测纠正图像上所有格网点的高程最大值;步骤2023、自适应搜索路径上检查点数量及各检查点位置确定根据公式对自适应搜索路径上的检查点数量Q进行确定;公式5中,GSD为自适应搜索路径上检查点的间距且GSDK1FB1,其中K115,FB1为当前所检测纠正图像的地面分辨率;对自适应搜索路径上Q个检查点的位置进行确定时,根据公式计算得出Q个检查点的二维平面坐标XQ,YQ;公式51中,Q为正整数且Q1、2、Q;步骤2024、遮蔽判断采用步骤2023中所确定的Q个检查点对当前检测点分别进行遮蔽判断当Q个检查点对当前检测点的遮蔽判断结果均为可见时,说明当前检测点为可见点;否则,说明当前检测点为。
9、遮蔽点;采用Q个检查点对当前检测点进行遮蔽判断时,所有检查点对当前检测点的遮蔽判断过程均相同;其中,利用Q个检查点中第Q个检查点对当前检测点进行遮蔽判断时,过程如下权利要求书CN104200527A3/5页4步骤、根据公式计算得出当前所采用检查点到当前检测点的水平距离;步骤、根据公式计算得出当前所采用检查点在当前检测点的摄影光线上的高程增量;公式7中,ZP为当前检测点的高程值;步骤、根据公式HZPDZ8,计算得出当前所采用检查点的遮蔽判断阈值H;步骤、将当前所采用检查点的高程值ZQ与步骤中所述的H进行差值比较当ZQH时,说明当前所采用检查点对当前检测点的遮蔽判断结果为遮蔽;当ZQH时,说明当前。
10、所采用检查点对当前检测点的遮蔽判断结果为可见;步骤2025、多次重复步骤2021至步骤2024,直至完成当前所检测纠正图像上所有待检测点的遮蔽检测过程;步骤三、影像纠正根据步骤二中的遮蔽检测结果,采用所述数据处理器对当前所检测纠正图像进行纠正;对当前所检测纠正图像进行纠正时,根据步骤二中对各待纠正点的遮蔽判断结果,对当前所检测纠正图像上的所有待纠正点分别进行纠正,步骤201中当前所检测纠正图像上的所有待检测点均为待纠正点;其中,当所纠正的待纠正点为可见点时,对待纠正点进行数字微分纠正并对该待纠正点的灰度值进行重采样;当所纠正的待纠正点为遮蔽点时,利用步骤一中航空摄影测量所摄取的与当前所检测纠正。
11、图像存在重叠区域的航空摄影图像,对待纠正点进行修复。2按照权利要求1所述的一种真正射影像的生成方法,其特征在于步骤三中修复完成后,还需对修复区域与其周侧原始影像区域的接边处进行羽化。3按照权利要求1或2所述的一种真正射影像的生成方法,其特征在于步骤201中所述的XI,YI为当前所检测纠正图像上四个角点的像点坐标观测数据;步骤2021中所述的XP,YP为当前检测点的地面坐标,并且按照公式对XP,YP进行计算;公式12中XP,YP为当前检测点的像点坐标观测数据,IP为当前待修复点在当前所检测纠正图像中所处的列数,JP为当前待修复点在当前所检测纠正图像中所处的行数;步骤中所述的ZP和步骤中所述的ZQ。
12、均为步骤一中所述原始测量数据中的测量值。4按照权利要求1或2所述的一种真正射影像的生成方法,其特征在于步骤2024中遮蔽判断完成后,还需根据遮蔽判断结果,对当前检测点进行标记当当前检测点为可见点时,将当前所检测点标记为可见;当当前检测点为遮蔽点时,将当前所检测点标记为遮蔽;步骤202中进行遮蔽检测之前,还需建立用于记录遮蔽检测结果的标记矩阵,所述标记矩阵为MN矩阵,且所述标记矩阵中记录有所述遮蔽检测范围内的MN个像素点的遮蔽判断结果。5按照权利要求1或2所述的一种真正射影像的生成方法,其特征在于步骤一中所摄取航空摄影图像的数量为多张,多张所述航空摄影图像之间存在重叠区域;步骤一中获取被测量区域。
13、的原始测量数据后,需按照步骤二中所述的遮蔽检测方法和步骤三中所述的权利要求书CN104200527A4/5页5影像纠正方法,对多张所述航空摄影图像分别进行遮蔽检测和影像纠正,并获得多张纠正后的航空摄影图像;之后,再对所获得的多张纠正后的航空摄影图像进行影像镶嵌,获得被测量区域的真正射影像。6按照权利要求5所述的一种真正射影像的生成方法,其特征在于步骤一中所摄取多张航空摄影图像的图像大小、摄影比例尺和像素大小均相同;步骤三中进行影像纠正之前,还从步骤一中所摄取的多张所述航空摄影图像中选出当前所检测纠正图像的副影像,所选取副影像的数量为K张,其中K为正整数且K2;当前所检测纠正图像为主影像;对所述。
14、副影像进行选取时,根据公式计算得出当前所选取航空摄影图像与所述主影像的重叠度,公式14为当前所选取航空摄影图像与所述主影像的投影中心之间的距离,IMAGESIZE为所述主影像的图像大小,PHOTOSCALE为所述主影像的摄影比例尺,PIXELSIZE为所述主影像的像素大小;并且,还需按照距离对所选取的多张副影像进行权值排序对任一张副影像的权值进行计算时,根据公式进行计算;公式141中,XJ,YJ为所选取的多张副影像中第J张副影像的摄影中心点坐标,WJ为第J张副影像的权值;之后,按照权值WJ由大到小的顺序对所选取的多张副影像进行排序;步骤三中进行影像纠正时,优先采用权值大的副影像进行修复。7按照。
15、权利要求5所述的一种真正射影像的生成方法,其特征在于步骤三中进行影像纠正之前,从步骤一中所摄取的多张航空摄影图像中,找出与主影像存在重叠区域的航空摄影图像作为副影像,所述副影像的数量为多张;其中,所述主影像为当前所检测纠正图像的原始影像且其为航空摄影图像;步骤三中利用步骤一中航空摄影测量所摄取的与当前所检测纠正图像存在重叠区域的航空摄影图像,对待纠正点进行修复时,过程如下步骤301、纠正用副影像确定先根据当前待修复点在所述主影像上的位置,找出能对当前待修复点进行修复的副影像,所述当前待修复点为当前所修复的待纠正点且其为遮蔽点;之后,对所找出的副影像是否能进行修复进行判断,过程如下步骤3011、。
16、矩形区域划分先获取当前所找出副影像的纠正图像,再在当前所找出副影像的纠正图像上划分出一个矩形区域;所述矩形区域为将当前所找出副影像的投影中心与所述主影像的摄影中心之间的连线作为对角线生成一个矩形区域,所述矩形区域的两个直边分别为当前所找出副影像的纠正图像的X轴和Y轴坐标轴;步骤3012、矩形区域内格网点遮蔽检测按照步骤2021至步骤2024中所述的遮蔽检测方法,对当前所找出副影像的纠正图像上位于所述矩形区域内且与当前待修复点对应的格网点进行遮蔽检测,当对应格网点的遮蔽判断结果为可见时,则选用当前所找出副影像对当前待修复点进行修复;否则,返回步骤3011,选择多张所述副影像中的其它副影像并对所找。
17、出的副影像是否能进行修复进行判断,直至找出对当前待修复点进行修复的副影像;权利要求书CN104200527A5/5页6步骤302、待纠正点修复利用步骤301中所确定的副影像,对当前待修复点进行修复;对当前待修复点进行修复时,按照公式计算得出当前待修复点在所述纠正用副影像上的像点位置,并对当前待修复点的灰度值进行重采样;公式11中XS,YS,ZS为所述纠正用副影像的摄影中心点坐标,ZT为当前待修复点的高程值,XT,YT为当前待修复点的地面坐标且公式13中XT,YT为当前待修复点的像点坐标观测数据,IT为当前待修复点在当前所检测纠正图像中所处的列数,JT为当前待修复点在当前所检测纠正图像中所处的行。
18、数。8按照权利要求7所述的一种真正射影像的生成方法,其特征在于步骤302中待纠正点修复完成后,还需对当前待修复点进行二次修复,过程如下步骤3031、副影像上修复点有效性判断判断步骤301中所确定副影像上对当前待修复点进行修复的当前修复点是否为有效修复点,其中所述当前修复点为步骤301中所确定副影像上对当前待修复点进行修复的修复点;判断当前修复点是否为有效修复点时,先在所述当前修复点投影方向上布设多个遮蔽检测点,并按照步骤2021至步骤2024中所述的遮蔽检测方法对所布设的各遮蔽检测点分别进行遮蔽检测,当且仅当所布设的多个遮蔽检测点均可见时,说明当前修复点是否为有效修复点,当前待修复点的修复过程。
19、完成;否则,需返回步骤3011,重新从多张所述副影像中选择其它副影像并对所找出的副影像是否能进行修复进行判断,直至找出对当前待修复点进行修复的副影像。9按照权利要求7所述的一种真正射影像的生成方法,其特征在于步骤三中对当前所检测纠正图像进行纠正过程中,对当前所检测纠正图像上的遮蔽点进行修复时,先以当前所检测纠正图像的投影中心为中心点,沿圆周方向将当前所检测纠正图像分割成多个图像块,再对多个所述图像块分别进行修复;并且,多个所述图像块的修复方法均相同;对任一个所述图像块进行修复时,均按照步骤301至步骤302中所述的方法对该图像块上的各遮蔽点进行修复。10按照权利要求7所述的一种真正射影像的生成。
20、方法,其特征在于步骤三中进行影像纠正时,还需利用所选取的多张所述副影像的投影中心所在图像区域的像素点对当前所检测纠正图像上的相应区域分别进行修复。权利要求书CN104200527A1/14页7一种真正射影像的生成方法技术领域0001本发明属于航空摄影测量技术领域,具体涉及一种真正射影像的生成方法。背景技术0002真正射影像制作是同时消除地形和地面三维目标投影变形的正射影像处理流程。传统的正射影像利用定向后的影像和数字高程模型DEM通过数字微分纠正技术获取,并未考虑人工建筑物等空间目标,从而导致其偏离真正的位置并遮挡其他地物目标。大比例尺城市地区正射影像中,建筑物倾斜、遮挡问题更为严重,如将GI。
21、S矢量地图叠加到正射影像上,会表现出道路矢量线横穿建筑物以及建筑物矢量轮廓无法与其自身叠合等问题,严重时将会导致其失去地理参考价值,不足以作为底图使用。0003随着摄影测量系统成像技术的不断提高以及我国城市化进程的高速发展,城区建筑物、道路桥梁日益复杂,传统的正射影像已经不能满足应用的需求,因此真正射影像这种高级的正射影像产品得到了越来越多的重视,它很好地解决了正射影像中建筑物倾斜、遮挡等问题,不仅具有地形图的特点,而且包含了丰富的纹理信息。近年来,国外的一些学者对真正射影像进行了较为广泛而深入的研究,出现了一些令人鼓舞的科研成果。我国在2006年的中国测绘学科发展蓝皮书中才提到“真正射影像”。
22、的概念,这种全新的测绘产品才开始引起普遍的关注。因此,研究真正射影像的制作以满足经济社会快速发展的需要对我国具有重要的现实意义。0004从国内外的研究成果上看,真正射影像制作主要利用数字建筑物模型DBM或数字表面模型DSM。真正射影像与正射影像最显著的差异在于纠正的同时分析地物的可见性,也称为遮蔽区域检测,检测遮蔽的能力是真正射影像生成算法的主要技术差别,国内外的研究人员也都将此作为重点议题讨论,遮蔽检测结果的好坏直接影响真正射影像的视觉效果。根据成像时投影射线的几何关系进行的遮蔽检测是理论上最严密的方法,典型的方法有比较重叠像点对应地物点到摄影中心距离的ZBUFFER方法;检测扫描方向上投影。
23、射线与天底方向夹角变化趋势的基于角度的检测方法;检测扫描方向上是否有采样点高于投影射线的基于高程的检测方法。其中,ZBUFFER方法在计算过程中需存储投影中心到目标点距离和与DSM具有同样分辨率和大小的二进制索引图两个主要的矩阵,运算效率低、占用的内存空间大,且对DSM分辨率敏感,存在伪遮蔽与伪可见问题,对于狭窄的垂直结构存在MPORTION问题。基于角度的检测方法可以有效地避免伪遮蔽和伪可见问题,但是扫描全部DSM复杂费时,重复访问现象不可避免,计算量大,效率低。而基于高程的检测方法不受影像分辨率的影响,基于共线方程理论严密精度较高,效率比基于角度的检测方法高,但是从总体上来看计算量还是比较。
24、大。另外,还有学者提出了基于矢量多边形的检测方法,该方法考虑了影像像素的投影变形,避免分辨率差异和地形坡度引起的伪遮蔽和伪可见,但是栅格矢量化工作量大,前向投影迭代计算耗时。因而,上述真正射影像制作方法均不同程度地存在效率低、实用性差等缺陷,不能有效解决真正射影像制作过程中存在的制作精度和速度问题,生产效率较低。说明书CN104200527A2/14页8发明内容0005本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种真正射影像的生成方法,其方法步骤简单、设计合理且生产效率高、使用效果好,能有效解决现有真正射影像制作方法存在的效率低、实用性差等问题。0006为解决上述技术问题,本发。
25、明采用的技术方案是一种真正射影像的生成方法,其特征在于该方法包括以下步骤0007步骤一、获取被测量区域的原始测量数据对被测量区域进行航空摄影测量,摄取被测量区域的航空摄影图像,并获得被测量区域的原始测量数据;所述原始测量数据为DSM测量数据或DBM测量数据,所述DSM测量数据为数字表面模型的测量数据,所述DBM测量数据为数字建筑物模型的测量数据;0008步骤二、遮蔽检测根据步骤一中所述原始测量数据,先采用数据处理器获取步骤一中所述航空摄影图像的纠正图像,所述纠正图像为根据当前所处理航空摄影图像在地面上的投影范围和需生成真正射影像的地面分辨率进行划分所形成的格网图像,所述格网图像的大小与需生成真。
26、正射影像的大小相同且所述格网图像上的各格网点分别与所述真正射影像上的各像素点一一对应,所述真正射影像为对所述航空摄影图像进行数字微分纠正和重采样后获得的影像;之后,再对所述纠正图像进行遮蔽检测,过程如下0009步骤201、遮蔽检测范围确定先根据当前所处理航空摄影图像的外方位元素和内方位元素,计算得出当前所处理航空摄影图像的四个角点的地面坐标;之后,根据计算得出的四个角点的地面坐标,确定当前所检测纠正图像的遮蔽检测范围;0010当前所处理航空摄影图像的四个角点分别为点A、点B、点C和点D,点A、点B、点C和点D的地面坐标分别为X1,Y1、X2,Y2、X3,Y3和X4,Y4;并且,根据公式计算得出。
27、四个角点的地面坐标;公式1中,I为正整数且I1、2、3、4;XS,YS,ZS为当前所处理航空摄影图像的所述外方位元素中的摄影中心点坐标,F为所述内方位元素中的一个参数且其为步骤一中对被测量区域进行航空摄影测量时所用航摄仪的焦距;XI,YI为当前所处理航空摄影图像上四个角点的二维平面坐标;Z为被测量区域的地面平均高度;0011计算得出当前所处理航空摄影图像上四个角点的地面坐标后,根据公式计算出所述纠正图像的大小;公式9中GRIDSIZE为预先设定的所述真正射影像的地面分辨率,所述真正射影像的大小为MN个像素点,所述纠正图像上包括MN个格网点;其中,M和N分别为所述纠正图像中格网点的列数和行数;0。
28、012所述遮蔽检测范围的横坐标为XMINXMAX且其纵坐标为YMINYMAX,其中XMIN和XMAX说明书CN104200527A3/14页9分别为四个角点的地面坐标中横坐标的最小值和最大值,YMIN和YMAX分别为四个角点的地面坐标中纵坐标的最小值和最大值;0013当前所检测纠正图像的旋转矩阵0014步骤202、遮蔽检测先将当前所检测纠正图像划分为多个遮蔽检测面元,每个所述检测面元中均包括N1N1个格网点,其中N1为正整数,且N1N11155;将当前所检测纠正图像上的各检测面元均作为一个待检测点分别进行遮蔽检测,并且各待检测点的遮蔽检测方法均相同;其中,对当前所检测纠正图像上的任一个待检测点。
29、进行遮蔽检测时,过程如下0015步骤2021、自适应搜索路径方位确定根据公式对当前检测点的自适应搜索路径的方位进行确定;其中,当前检测点为当前状态下进行遮蔽检测的待检测点;公式2中,为当前检测点相对于投影中心的方位角,所述投影中心为所述航空摄影图像的摄影中心点在地面上的投影点;当N11时,XP,YP为当前检测点的二维平面坐标;当N125时,XP,YP为当前检测点中所包括N1个格网点的平均二维平面坐标;0016所述自适应搜索路径为对当前检测点进行遮蔽检测时所采用的搜索路径;0017步骤2022、自适应搜索路径长度确定先根据公式计算得出当前检测点到投影中心的水平距离D;再根据公式计算得出当前检测点。
30、的自适应搜索路径的长度L;公式4中,L为当前检测点到投影中心的搜索长度,HMAX为当前所检测纠正图像上所有格网点的高程最大值;0018步骤2023、自适应搜索路径上检查点数量及各检查点位置确定根据公式对自适应搜索路径上的检查点数量Q进行确定;公式5中,GSD为自适应搜索路径上检查点的间距且GSDK1FB1,其中K115,FB1为当前所检测纠正图像的地面分辨率;0019对自适应搜索路径上Q个检查点的位置进行确定时,根据公式计算得出Q个检查点的二维平面坐标XQ,YQ;公式51中,Q为正整数且Q1、2、Q;0020步骤2024、遮蔽判断采用步骤2023中所确定的Q个检查点对当前检测点分别进说明书CN。
31、104200527A4/14页10行遮蔽判断当Q个检查点对当前检测点的遮蔽判断结果均为可见时,说明当前检测点为可见点;否则,说明当前检测点为遮蔽点;0021采用Q个检查点对当前检测点进行遮蔽判断时,所有检查点对当前检测点的遮蔽判断过程均相同;其中,利用Q个检查点中第Q个检查点对当前检测点进行遮蔽判断时,过程如下0022步骤、根据公式计算得出当前所采用检查点到当前检测点的水平距离;0023步骤、根据公式计算得出当前所采用检查点在当前检测点的摄影光线上的高程增量;公式7中,ZP为当前检测点的高程值;0024步骤、根据公式HZPDZ8,计算得出当前所采用检查点的遮蔽判断阈值H;0025步骤、将当前所。
32、采用检查点的高程值ZQ与步骤中所述的H进行差值比较当ZQH时,说明当前所采用检查点对当前检测点的遮蔽判断结果为遮蔽;当ZQH时,说明当前所采用检查点对当前检测点的遮蔽判断结果为可见;0026步骤2025、多次重复步骤2021至步骤2024,直至完成当前所检测纠正图像上所有待检测点的遮蔽检测过程;0027步骤三、影像纠正根据步骤二中的遮蔽检测结果,采用所述数据处理器对当前所检测纠正图像进行纠正;0028对当前所检测纠正图像进行纠正时,根据步骤二中对各待纠正点的遮蔽判断结果,对当前所检测纠正图像上的所有待纠正点分别进行纠正,步骤201中当前所检测纠正图像上的所有待检测点均为待纠正点;其中,当所纠正。
33、的待纠正点为可见点时,对待纠正点进行数字微分纠正并对该待纠正点的灰度值进行重采样;当所纠正的待纠正点为遮蔽点时,利用步骤一中航空摄影测量所摄取的与当前所检测纠正图像存在重叠区域的航空摄影图像,对待纠正点进行修复。0029上述一种真正射影像的生成方法,其特征是步骤三中修复完成后,还需对修复区域与其周侧原始影像区域的接边处进行羽化。0030上述一种真正射影像的生成方法,其特征是步骤201中所述的XI,YI为当前所检测纠正图像上四个角点的像点坐标观测数据;0031步骤2021中所述的XP,YP为当前检测点的地面坐标,并且按照公式对XP,YP进行计算;公式12中XP,YP为当前检测点的像点坐标观测数据。
34、,IP为当前待修复点在当前所检测纠正图像中所处的列数,JP为当前待修复点在当前所检测纠正图像中所处的行数;0032步骤中所述的ZP和步骤中所述的ZQ均为步骤一中所述原始测量数据中的测量值。0033上述一种真正射影像的生成方法,其特征是步骤2024中遮蔽判断完成后,还需根据遮蔽判断结果,对当前检测点进行标记当当前检测点为可见点时,将当前所检测点标说明书CN104200527A105/14页11记为可见;当当前检测点为遮蔽点时,将当前所检测点标记为遮蔽;0034步骤202中进行遮蔽检测之前,还需建立用于记录遮蔽检测结果的标记矩阵,所述标记矩阵为MN矩阵,且所述标记矩阵中记录有所述遮蔽检测范围内的M。
35、N个像素点的遮蔽判断结果。0035上述一种真正射影像的生成方法,其特征是步骤一中所摄取航空摄影图像的数量为多张,多张所述航空摄影图像之间存在重叠区域;步骤一中获取被测量区域的原始测量数据后,需按照步骤二中所述的遮蔽检测方法和步骤三中所述的影像纠正方法,对多张所述航空摄影图像分别进行遮蔽检测和影像纠正,并获得多张纠正后的航空摄影图像;之后,再对所获得的多张纠正后的航空摄影图像进行影像镶嵌,获得被测量区域的真正射影像。0036上述一种真正射影像的生成方法,其特征是步骤一中所摄取多张航空摄影图像的图像大小、摄影比例尺和像素大小均相同;0037步骤三中进行影像纠正之前,还从步骤一中所摄取的多张所述航空。
36、摄影图像中选出当前所检测纠正图像的副影像,所选取副影像的数量为K张,其中K为正整数且K2;0038当前所检测纠正图像为主影像;对所述副影像进行选取时,根据公式计算得出当前所选取航空摄影图像与所述主影像的重叠度,公式14为当前所选取航空摄影图像与所述主影像的投影中心之间的距离,IMAGESIZE为所述主影像的图像大小,PHOTOSCALE为所述主影像的摄影比例尺,PIXELSIZE为所述主影像的像素大小;0039并且,还需按照距离对所选取的多张副影像进行权值排序对任一张副影像的权值进行计算时,根据公式进行计算;公式141中,XJ,YJ为所选取的多张副影像中第J张副影像的摄影中心点坐标,WJ为第J。
37、张副影像的权值;之后,按照权值WJ由大到小的顺序对所选取的多张副影像进行排序;0040步骤三中进行影像纠正时,优先采用权值大的副影像进行修复。0041上述一种真正射影像的生成方法,其特征是步骤三中进行影像纠正之前,从步骤一中所摄取的多张航空摄影图像中,找出与主影像存在重叠区域的航空摄影图像作为副影像,所述副影像的数量为多张;其中,所述主影像为当前所检测纠正图像的原始影像且其为航空摄影图像;0042步骤三中利用步骤一中航空摄影测量所摄取的与当前所检测纠正图像存在重叠区域的航空摄影图像,对待纠正点进行修复时,过程如下0043步骤301、纠正用副影像确定先根据当前待修复点在所述主影像上的位置,找出能。
38、对当前待修复点进行修复的副影像,所述当前待修复点为当前所修复的待纠正点且其为遮蔽点;之后,对所找出的副影像是否能进行修复进行判断,过程如下0044步骤3011、矩形区域划分先获取当前所找出副影像的纠正图像,再在当前所找出副影像的纠正图像上划分出一个矩形区域;所述矩形区域为将当前所找出副影像的投影中心与所述主影像的摄影中心之间的连线作为对角线生成一个矩形区域,所述矩形区域的说明书CN104200527A116/14页12两个直边分别为当前所找出副影像的纠正图像的X轴和Y轴坐标轴;0045步骤3012、矩形区域内格网点遮蔽检测按照步骤2021至步骤2024中所述的遮蔽检测方法,对当前所找出副影像的。
39、纠正图像上位于所述矩形区域内且与当前待修复点对应的格网点进行遮蔽检测,当对应格网点的遮蔽判断结果为可见时,则选用当前所找出副影像对当前待修复点进行修复;否则,返回步骤3011,选择多张所述副影像中的其它副影像并对所找出的副影像是否能进行修复进行判断,直至找出对当前待修复点进行修复的副影像;0046步骤302、待纠正点修复利用步骤301中所确定的副影像,对当前待修复点进行修复;0047对当前待修复点进行修复时,按照公式计算得出当前待修复点在所述纠正用副影像上的像点位置,并对当前待修复点的灰度值进行重采样;0048公式11中XS,YS,ZS为所述纠正用副影像的摄影中心点坐标,ZT为当前待修复点的高。
40、程值,XT,YT为当前待修复点的地面坐标且公式13中XT,YT为当前待修复点的像点坐标观测数据,IT为当前待修复点在当前所检测纠正图像中所处的列数,JT为当前待修复点在当前所检测纠正图像中所处的行数。0049上述一种真正射影像的生成方法,其特征是步骤302中待纠正点修复完成后,还需对当前待修复点进行二次修复,过程如下0050步骤3031、副影像上修复点有效性判断判断步骤301中所确定副影像上对当前待修复点进行修复的当前修复点是否为有效修复点,其中所述当前修复点为步骤301中所确定副影像上对当前待修复点进行修复的修复点;判断当前修复点是否为有效修复点时,先在所述当前修复点投影方向上布设多个遮蔽检。
41、测点,并按照步骤2021至步骤2024中所述的遮蔽检测方法对所布设的各遮蔽检测点分别进行遮蔽检测,当且仅当所布设的多个遮蔽检测点均可见时,说明当前修复点是否为有效修复点,当前待修复点的修复过程完成;否则,需返回步骤3011,重新从多张所述副影像中选择其它副影像并对所找出的副影像是否能进行修复进行判断,直至找出对当前待修复点进行修复的副影像。0051上述一种真正射影像的生成方法,其特征是步骤三中对当前所检测纠正图像进行纠正过程中,对当前所检测纠正图像上的遮蔽点进行修复时,先以当前所检测纠正图像的投影中心为中心点,沿圆周方向将当前所检测纠正图像分割成多个图像块,再对多个所述图像块分别进行修复;并且。
42、,多个所述图像块的修复方法均相同;对任一个所述图像块进行修复时,均按照步骤301至步骤302中所述的方法对该图像块上的各遮蔽点进行修复。0052上述一种真正射影像的生成方法,其特征是步骤三中进行影像纠正时,还需利用所选取的多张所述副影像的投影中心所在图像区域的像素点对当前所检测纠正图像上的相应区域分别进行修复。说明书CN104200527A127/14页130053本发明与现有技术相比具有以下优点00541、方法步骤简单、设计合理且操作简单,实现方便。00552、生成效率较高。00563、所采用的遮蔽检测方法步骤简单、设计合理且检测速度快、检测精度高,采用基于方位关系的自适应路径搜索和面元的遮。
43、蔽判断方法,克服了传统遮蔽检测方法的效率和精度问题。00574、所采用的遮蔽点修复方法步骤简单、设计合理且修复速度快、修复效果好,对遮蔽区域进行修复时,采用两次补偿,第一次采用分块对向补偿方法进行补偿,第二次采用中心补偿方法即基于投影中心的中心补偿方法进行中心补偿。并且,所述中心补偿方法是对分块对向补偿方法的补充,对于经分块对向补偿方法进行修复后,当前所检测纠正图像上存留的遮蔽点进行修复。并且,实际进行分块对向补偿时,采用两次修复相结合的方法,不仅能对各遮蔽点进行有效修复,而且实现方便,省时省力,并能有效克服现有遮蔽点补偿过程中容易出现的“双眼皮”现象,从而大大改善了补偿的效果。另外,对副影像。
44、的选用优先权进行排序,进一步节省了人力物力,花费时间少,并且修复效果好。00585、使用效果好且实用价值高、实用性强,主要包括以下步骤1获取高精度DBM测量数据或DSM测量数据,并计算遮蔽检测范围,确定纠正图像的格网大小和纠正面元大小;2计算待检测点遮蔽判断的自适应搜索路径;3采用基于摄影光线高度的自适应方法对每个小面元的可见性进行遮蔽检测,并作可见或不可见的标记;4对可见的面元做数字微分纠正,对不可见区域进行标记;5对不可见区域,按照方位并结合遮蔽检测的方法利用副影像数据进行首次补偿;6对剩余的不可见区域,利用各个副影像投影中心附近的像素进行二次补偿;7将补偿的遮蔽区域的周围和原始影像之间做。
45、羽化过度,完成真正射影像的制作。本发明具有实际生产效率高、精度高等优点,所生产的真正射影像产品可用于城市影像地图的制作、三维建模、城市规划等领域。0059综上所述,本发明方法步骤简单、设计合理且生产效率高、使用效果好,能有效解决现有真正射影像制作方法存在的效率低、实用性差等问题,以解决在制作真正射影像中的精度和速度问题,提高生产效率。0060下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。附图说明0061图1为本发明的方法流程框图。0062图2为采用本发明进行遮蔽检测时的方法流程框图。0063图3为采用本发明进行遮蔽补偿时的方法流程框图。具体实施方式0064如图1所示的一种真正射影。
46、像的生成方法,包括以下步骤0065步骤一、步骤一、获取被测量区域的原始测量数据对被测量区域进行航空摄影测量,摄取被测量区域的航空摄影图像,并获得被测量区域的原始测量数据;所述原始测量数据为DSM测量数据或DBM测量数据,所述DSM测量数据为数字表面模型的测量数据,所述DBM测量数据为数字建筑物模型的测量数据。说明书CN104200527A138/14页140066其中,数字表面模型DIGITALSURFACEMODEL,缩写DSM是指包含了地表建筑物、桥梁和树木等高度的地面高程模型。0067由于高精度的DSM测量数据是制作真正射影像产品的基础,通过采集或者自动匹配的方法获取DSM测量数据,并对。
47、DSM测量数据中建筑物的边界和顶面进行处理,使其边界闭合顶面特征明显。0068步骤二、遮蔽检测结合图2,根据步骤一中所述原始测量数据,先采用数据处理器获取步骤一中所述航空摄影图像的纠正图像,所述纠正图像为根据当前所处理航空摄影图像在地面上的投影范围和需生成真正射影像的地面分辨率进行划分所形成的格网图像,所述格网图像的大小与需生成真正射影像的大小相同且所述格网图像上的各格网点分别与所述真正射影像上的各像素点一一对应,所述真正射影像为对所述航空摄影图像进行数字微分纠正和重采样后获得的影像,所述航空摄影图像为所述真正射影像的原始影像,所述格网图像上的各格网点分别与所述原始影像上的各像素点一一对应;之。
48、后,再对所述纠正图像进行遮蔽检测,过程如下0069步骤201、遮蔽检测范围确定先根据当前所处理航空摄影图像的外方位元素和内方位元素,计算得出当前所处理航空摄影图像的四个角点的地面坐标;之后,根据计算得出的四个角点的地面坐标,确定当前所检测纠正图像的遮蔽检测范围。0070当前所处理航空摄影图像的四个角点分别为点A、点B、点C和点D,点A、点B、点C和点D的地面坐标分别为X1,Y1、X2,Y2、X3,Y3和X4,Y4;并且,根据公式计算得出四个角点的地面坐标;公式1中,I为正整数且I1、2、3、4;XS,YS,ZS为当前所处理航空摄影图像的所述外方位元素中的摄影中心点坐标,F为所述内方位元素中的一。
49、个参数且其为步骤一中对被测量区域进行航空摄影测量时所用航摄仪的焦距;XI,YI为当前所处理航空摄影图像上四个角点的二维平面坐标;Z为被测量区域的地面平均高度。0071计算得出当前所处理航空摄影图像上四个角点的地面坐标后,根据公式计算出所述纠正图像的大小;公式9中GRIDSIZE为预先设定的所述真正射影像的地面分辨率,所述真正射影像的大小为MN个像素点,所述纠正图像上包括MN个格网点,也就是说,所述纠正图像上的各格网点均对应一个像素点;其中,M和N分别为所述纠正图像中格网点的列数和行数。0072所述遮蔽检测范围的横坐标为XMINXMAX且其纵坐标为YMINYMAX,其中XMIN和XMAX分别为四个角点的地面坐标中横坐标的最小值和最大值,YMIN和YMAX分别为四个角点的地面坐标中纵坐标的最小值和最大值。其中,XMIN和XMAX分别为四个角点的地面坐标中最小的横坐标值和最大的横坐标值,YMIN和YMAX分别为四个角点的地面坐标中最小的纵坐标值和最大说明书CN104200527A149/14页15的纵坐标值。0073当前所检测纠正图像的旋转矩阵0074实际进行航空测量时,所述外方位元素包括外方位角元素和以及摄影中心点坐标XS,YS,ZS;所述内方位元素包括X0、Y0和。