一种高分辨率影像的镶嵌处理方法技术

本发明专利技术公开了一种高分辨率影像的镶嵌处理方法，利用航摄区范围内的道路、河流地理实体中心线数据作为区域整体镶嵌线的骨架，结合地表覆盖数据叠加相幅面数据进行重叠边缘计算，获取相幅重叠边缘区相交的完整地表覆盖地理实体范围数据，在保证地理实体完整性的情况下对相幅面进行拓扑重构计算，并将计算结果与区域道路骨架线再次进行拓扑运算，得到最小分割区域作为图斑成果数据，最后对重复图斑数据进行像主点距离判断，取舍最小距离图斑，从而消除重复图斑数据，最终计算得到覆盖全航摄区的镶嵌分割小图斑数据。试验结果表明本发明专利技术的方法具有较好的数据处理效果，方法简单且易于实现，重点应用可在高分辨率航空影像镶嵌处理中推广。

【技术实现步骤摘要】
一种高分辨率影像的镶嵌处理方法
本专利技术属于遥感影像处理
，具体涉及一种高分辨率影像的镶嵌处理方法。
技术介绍
在遥感影像处理中为了获取更大范围、更为完整的影像数据，需将多幅(景)遥感图像进行拼接处理，在这个过程中，影像镶嵌是最为关键的一步。遥感影像数据镶嵌处理是将两幅或多幅遥感影像数据按照镶嵌规则拼接形成整体遥感影像数据的技术过程。由于遥感影像数据源采用航空、航天遥感等不同方式获取，成像方法和数据特征均有很大不同，因此对数据源进行裁剪、去重、融合过程中需要制定适合的镶嵌规则，尽量避免出现因灰度或颜色差异而产生的拼接缝，减少因不同成像条件产生的影像投影差与变形。目前高分辨率数字正射影像数据(DOM)生产主要采用数字摄影测量系统进行。采用数字摄影测量系统生产航空正射影像数据时，影响影像数据成果精度的主要因素包括航摄相机性能、航摄比例尺、航摄质量、航片扫描质量、控制点精度、DEM的平面精度等。常规航空数字摄影测量法生产流程如图1所示，由于在航空影像获取过程中的各种环境因素不同，使得每条航带里面的影像和航带间相互重叠、连接的影像存在色差、亮度差、构筑物投影差等多方面不同程度的差异。因此在生产正射影像数据中需要用专业化的遥感影像处理软件来对影像数据进行处理，特别是在影像镶嵌处理环节，由于镶嵌线的计算不合理，常常会导致大量的手工镶嵌处理作业，极大影响了镶嵌处理效率。
技术实现思路
专利技术目的：针对现有技术中存在的不足，本专利技术的目的是一种高分辨率影像的镶嵌处理方法，利用道路、地理国情覆盖数据，编写算法实现遥感影像镶嵌线自动计算、影像自动裁切与拼接、标准分幅影像制作等自动化处理流程。本文研究的影像数据源主要是影像分辨率优于20公分的高分辨率影像数据。技术方案：为了实现上述专利技术目的，本专利技术采用的技术方案为：一种高分辨率影像的镶嵌处理方法，利用航摄区范围内的交通道路地理实体数据作为区域整体镶嵌线的骨架，结合地表覆盖数据叠加相幅面数据进行重叠边缘计算，获取相幅重叠边缘区相交的完整地表覆盖地理实体数据，在保证地理实体完整性的情况下对相幅面进行拓扑重构计算，并将计算结果与区域道路骨架线再次进行拓扑运算，得到最小分割区域作为图斑成果数据，最后对重复图斑数据进行像主点距离判断，取舍最小距离图斑，从而消除重复图斑数据，最终计算得到覆盖全航摄区的镶嵌分割小图斑数据。所述的高分辨率影像的镶嵌处理方法，包括以下步骤：1)获取航摄相机参数和航线航带信息；2)生成单片相幅面图形信息；3)利用航摄区范围内的交通道路地理实体数据作为区域整体镶嵌线的骨架，结合地表覆盖数据叠加相幅面数据进行重叠边缘计算，获取相幅重叠边缘区相交的完整地表覆盖地理实体数据；4)拓扑计算后进行构面处理：当拓扑面没有出现一个完整的地理实体对象时，关联相幅边线检索、航摄区地表覆盖地理实体，进入步骤3)运算；当拓扑面中出现完整地理实体对象，且不同相幅面与同一地理实体构出了唯一完整面，输出镶嵌线构面成果数据；不同相幅面与同一地理实体构出了多个完整面，进行像主点与地理实体距离计算，选距离最短、投影差最小的构面数据，输出镶嵌线构面成果数据。步骤3)中，重叠边缘计算获取相幅重叠边缘区相交的完整地表覆盖范围数据的具体过程为：(1)按照ID号由小到大的顺序从道路、河流地理实体数据构成的骨架的多边形集合X中每次选择1个多边形b，再按照单片相幅面的多边形集合Y的ID号顺序由大到小每次选择1个多边形a1，先判断多边形b是否完全包含在多边形a1内；如果包含，将a1记录到完全包含b的多边形集合Ab1中，如果不包含，则判断多边形b与多边形a1是否相交，如果相交，将a1记录到与b相交的多边形集合Bb1中，如果既不包含也不相交，判断结束；按照ID顺序选择多边形集合Y中下一个多边形a2，继续判断包含或者相交，直到多边形集合Y中的所有多边形遍历完成，得到与多边形b有包含关系或者相交关系的多边形的集合Ab1与集合Bb1；按照ID由小到大的顺序从多边形集合X中选择下一个多边形b2，继续遍历判断b2与多边形集合Y中的所有多边形的关系，得到相应的集合Ab2与集合Bb2；按照ID号由小到大的顺序逐个判断完集合X中所有的多边形，直至得到多边形集合Abn与Bbn；(2)按照步骤(1)得到的集合AB{Ab1、Ab2、......Abn}，从由小到大的ID顺序依次取出Ab1、Ab1......Abn集合，判断Ab1中的多边形记录值与1的大小关系；如果等于1，则得到唯一包含多边形b的相幅an，即相幅包含的唯一完整地表覆盖范围，记入多边形集合P1，其对应的相幅记入多边形集合P2；如果大于1，则按照记录ID由大到小的顺序每次从集合Abn中取1个多边形，分别计算其对应的像主点与多边形b1的平面几何中心点距离集合S，按照排序取min{S}对应的相幅an，同时记录对应多边形b，此时完全包含条件下的裁切线就是多边形b的范围，其对应的相幅就是min{S}，即an。步骤4)中，将计算结果与区域道路骨架线再次进行拓扑运算的具体过程为：(1)按照ID由小到大的顺序遍历集合{Bb1、Bb2、......Bbn}中的多边形，逐一计算与b1的交集，从而得到新的多边形集合{Bb11、Bb21、......Bbn1}，取集合{Bb11、Bb21、......Bbn1}中面积最大且唯一多边形记为MaxAreaPolygon；(2)取地理国情图斑数据，逐一与多边形MaxAreaPolygon计算包含关系，得到所有包含在多边形MaxAreaPolygon内部的图斑多边形数据，重新合并组成新的多边形TuB；(3)将多边形TuB按照第3条中步骤(1)、步骤(2)的方法计算，得到TuB相应的唯一包含相幅an，即相幅重叠包含的唯一完整地表覆盖范围，此时将TuB记入多边形集合P1，其对应相幅记入多边形集合P2；(4)运算到Ab1中的多边形记录值为1，即所有多边形TuBn与相幅都是包含关系结束；此时得到最终多边形集合P1与P2。步骤5)中，像主点与地理实体距离计算的具体过程为：(1)计算地理实体几何中心点point_A；(2)采用空间距离计算公式计算像主点point_B与point_A的空间距离；空间距离计算公式为：有益效果：与现有技术相比，本专利技术的高分辨率影像的镶嵌处理方法，针对航空正摄影像数据镶嵌处理环节进行了研究，在充分考虑数据源的成像方式和特征的情况下，利用航摄区道路与地表覆盖地理实体数据，使用镶嵌线生成算法，将航带航向、像幅、像主点、地理实体等参数作为影像重叠区数据取舍的依据，在确保投影变形最小的情况下实现影像镶嵌的自动化处理；本专利技术通过大范围航空影像数据进行实际试验，结果表明本专利技术的方法具有较好的数据处理效果，方法简单且易于实现，重点应用可在高分辨率(优于20公分)航空影像镶嵌处理中推广。附图说明图1是现有的空数字摄影测量法生产流程图；图2是利用航摄参数和地理实体计算镶嵌线的算法流程图；图3是影像单片面重叠情况图；图4是影像单片面与道路、地表覆盖地理实体数据拓扑计算完成后的图斑成果图；图5是利用图斑成果自动裁切、镶嵌后的正摄影像成果图。具体实施方式下面结合具体实施例对本专利技术做进一步的说明。实施例11、影像镶嵌、裁切处理环节改进本实施例使用的影像数据源本文档来自技高网...

【技术保护点】
1.一种高分辨率影像的镶嵌处理方法，其特征在于，利用航摄区范围内的道路、河流地理实体数据作为区域整体镶嵌线的骨架，结合地表覆盖数据叠加相幅面数据进行重叠边缘计算，获取相幅重叠边缘区相交的完整地表覆盖地理实体数据，在保证地理实体完整性的情况下对相幅面进行拓扑重构计算，并将计算结果与区域道路骨架线再次进行拓扑运算，得到最小分割区域作为图斑成果数据，最后对重复图斑数据进行像主点距离判断，取舍最小距离图斑，从而消除重复图斑数据，最终计算得到覆盖全航摄区的镶嵌分割小图斑数据。

【技术特征摘要】
1.一种高分辨率影像的镶嵌处理方法，其特征在于，利用航摄区范围内的道路、河流地理实体数据作为区域整体镶嵌线的骨架，结合地表覆盖数据叠加相幅面数据进行重叠边缘计算，获取相幅重叠边缘区相交的完整地表覆盖地理实体数据，在保证地理实体完整性的情况下对相幅面进行拓扑重构计算，并将计算结果与区域道路骨架线再次进行拓扑运算，得到最小分割区域作为图斑成果数据，最后对重复图斑数据进行像主点距离判断，取舍最小距离图斑，从而消除重复图斑数据，最终计算得到覆盖全航摄区的镶嵌分割小图斑数据。2.根据权利要求1所述的高分辨率影像的镶嵌处理方法，其特征在于，包括以下步骤：1)获取航摄相机参数和航线航带信息；2)生成单片相幅面图形信息；3)利用航摄区范围内的道路、河流地理实体数据作为区域整体镶嵌线的骨架，结合地表覆盖数据叠加相幅面数据进行重叠边缘计算，获取相幅重叠边缘区相交的完整地表覆盖地理实体范围数据；4)在保证地理实体完整性的情况下对相幅面进行拓扑重构计算，并将计算结果与区域道路骨架线再次进行拓扑运算，得到最小分割区域作为图斑成果数据；5)拓扑计算后进行构面处理：当拓扑面没有出现一个完整的地理实体对象时，关联相幅边线检索、航摄区地表覆盖地理实体，进入步骤3)运算；当拓扑面中出现完整地理实体对象，且不同相幅面与同一地理实体构出了唯一完整面，输出镶嵌线构面成果数据；不同相幅面与同一地理实体构出了多个完整面，进行像主点与地理实体距离计算，选距离最短、投影差最小的构面数据，输出镶嵌线构面成果数据。3.根据权利要求1或2所述的高分辨率影像的镶嵌处理方法，其特征在于，步骤3)中，重叠边缘计算获取相幅重叠边缘区相交的完整地表覆盖范围数据的具体过程为：(1)按照ID号由小到大的顺序从道路、河流地理实体数据构成的骨架的多边形集合X中每次选择1个多边形b，再按照单片相幅面的多边形集合Y的ID号顺序由大到小每次选择1个多边形a1，先判断多边形b是否完全包含在多边形a1内；如果包含，将a1记录到完全包含b的多边形集合Ab1中，如果不包含，则判断多边形b与多边形a1是否相交，如果相交，将a1记录到与b相交的多边形集合Bb1中，如果既不包含也不相交，判断结束；按照ID顺序选择多边形集合Y中下一个多边形a2，继续判断包含或者相交，直到多边形集合Y中的所有多边形遍历完成，得到与多边形b有包含关系或者相交关系的多边形的集合Ab1与集合...

【专利技术属性】
技术研发人员：宋玉兵，张光伟，王朝辉，羌鑫林，
申请(专利权)人：江苏省测绘工程院，
类型：发明
国别省市：江苏,32