数字正射影像拉花区域的处理方法技术

本发明专利技术公开了一种数字正射影像拉花区域的处理方法，包括步骤：一、获取被测量区域的原始测量数据；二、正射影像的地面格网范围确定；三、拉花检测：对当前所处理航空摄影图像的纠正图像上的各格网点分别进行拉花检测，过程如下：301、投影光线方位确定；302、拉花判断：判断沿所确定投影光线上是否存在其它与当前检测点在原始影像上对应的像素点相同的格网点：当存在时，说明当前检测点为“拉花”点；反之，说明当前检测点为未“拉花”点；四、影像纠正：根据拉花检测结果，对当前所处理航空摄影图像进行纠正。本发明专利技术方法步骤简单、设计合理且效率高、使用效果好，能解决现有正射影像的拉花处理方法存在的效率低、使用效果较差等问题。

【技术实现步骤摘要】
【专利摘要】本专利技术公开了一种，包括步骤：一、获取被测量区域的原始测量数据；二、正射影像的地面格网范围确定；三、拉花检测：对当前所处理航空摄影图像的纠正图像上的各格网点分别进行拉花检测，过程如下：301、投影光线方位确定；302、拉花判断：判断沿所确定投影光线上是否存在其它与当前检测点在原始影像上对应的像素点相同的格网点：当存在时，说明当前检测点为“拉花”点；反之，说明当前检测点为未“拉花”点；四、影像纠正：根据拉花检测结果，对当前所处理航空摄影图像进行纠正。本专利技术方法步骤简单、设计合理且效率高、使用效果好，能解决现有正射影像的拉花处理方法存在的效率低、使用效果较差等问题。【专利说明】
本专利技术属于航空摄影测量
，具体涉及一种。
技术介绍
正射影像制作是指同时消除地形起伏和相片倾斜引起的影像变形的过程。正射影像是根据有关参数和数字高程模型(DEM)，利用相应的构像方程式，计算地面点对应的像点坐标，并对原始影像进行灰度重采样，使得成为一幅既有正确平面位置又有丰富纹理信息的影像的技术。 正射影像是对原始影像进行微分纠正和灰度重采样的结果图，而原始影像由于中心投影和地形起伏的影响并不能保证地面上的每个位置都能在影像中成像。因此，在数字微分纠正中，重采样时对于摄影信息充足的区域采样较为稀疏，对于摄影信息匮乏区域采样过密或者重复采样，从而导致图像出现颗粒感或者拉伸的现象，沿一个方向拉伸过度时会出现区域纹理失真的现象，我们称之为“拉花”现象。根据中华人民共和国测绘行业标准《CH/T9008.3-2010》中“基础地理信息数字成果1:500、1:1000、1:2000数字正射影像图”的规定，数字正射影像的质量元素主要有空间参考系、位置精度、逻辑一致性、时间精度、影像质量、元数据质量、表征质量和附件质量。而“拉花”问题的存在严重影响了质量元素的位置精度和影像质量。 目前，在国内外文献和专利中，还没有针对“拉花”提出的自动化解决方案，在通常的生产过程中，当出现拉花现象的正射影像时，一般都要通过利用人工的方法对“拉花”区域进行寻找，并到相应的纠正影像图或者原图当中去定位，寻找相应的纹理对“拉花”区域进行替补，以此来解决“拉花”问题，但这种人工的处理方法非常费时费力，特别是在“拉花”区域比较多的情况下，生产效率非常低。
技术实现思路
本专利技术所要解决的技术问题在于针对上述现有技术中的不足，提供一种，其方法步骤简单、设计合理且效率高、使用效果好，能有效解决现有正射影像的拉花处理方法存在的效率低、使用效果较差等问题。 为解决上述技术问题，本专利技术采用的技术方案是:一种，其特征在于该方法包括以下步骤: 步骤一、获取被测量区域的原始测量数据:对被测量区域进行航空摄影测量，摄取被测量区域的航空摄影图像，并获得被测量区域的原始测量数据；所述原始测量数据为DEM测量数据，所述DEM测量数据为数字高程数据； 步骤二、正射影像的地面格网范围确定:先根据步骤一中所述原始测量数据，获取步骤一中当前所处理航空摄影图像的纠正图像，所述纠正图像为根据当前所处理航空摄影图像在地面上的投影范围和需生成正射影像的地面分辨率进行划分所形成的格网图像，所述格网图像的大小与当前所处理航空摄影图像的正射影像的大小相同且所述格网图像上的各格网点分别与所述正射影像上的各像素点一一对应，所述正射影像为对所述航空摄影图像进行数字微分纠正和重采样后获得的影像，所述航空摄影图像为所述正射影像的原始影像；之后，根据当前所处理航空摄影图像的外方位元素和内方位元素，计算得出当前所处理航空摄影图像的四个角点的地面坐标；然后，根据计算得出的四个角点的地面坐标，确定所述纠正图像的格网范围； 当前所处理航空摄影图像的四个角点分别为点A、点B、点C和点D，点A、点B、点C和点D的地面坐标分别为(X1, Y1)、(X2, Y2)、(X3, Y3)和(X4，Y4)；根据公式 【权利要求】1.一种，其特征在于该方法包括以下步骤: 步骤一、获取被测量区域的原始测量数据:对被测量区域进行航空摄影测量，摄取被测量区域的航空摄影图像，并获得被测量区域的原始测量数据；所述原始测量数据为DEM测量数据，所述DEM测量数据为数字高程数据； 步骤二、正射影像的地面格网范围确定:先根据步骤一中所述原始测量数据，获取步骤一中当前所处理航空摄影图像的纠正图像，所述纠正图像为根据当前所处理航空摄影图像在地面上的投影范围和需生成正射影像的地面分辨率进行划分所形成的格网图像，所述格网图像的大小与当前所处理航空摄影图像的正射影像的大小相同且所述格网图像上的各格网点分别与所述正射影像上的各像素点一一对应，所述正射影像为对所述航空摄影图像进行数字微分纠正和重采样后获得的影像，所述航空摄影图像为所述正射影像的原始影像；之后，根据当前所处理航空摄影图像的外方位元素和内方位元素，计算得出当前所处理航空摄影图像的四个角点的地面坐标；然后，根据计算得出的四个角点的地面坐标，确定所述纠正图像的格网范围； 当前所处理航空摄影图像的四个角点分别为点A、点B、点C和点D,点A、点B、点C和点D的地面坐标分别为(X1, Y1)、(X2, Y2)、(X3, Y3)和(X4，Y4);根据公式(1)，计算得出四个角点的地面坐标；公式(I)中，I为正整数且i = 1、2、3、4;(XS，Ys, Zs)为当前所处理航空摄影图像的所述外方位元素中的摄影中心点坐标，f为所述内方位元素中的一个参数且其为步骤一中对被测量区域进行航空摄影测量时所用航摄仪的焦距；(Xi，yi)为当前所处理航空摄影图像上四个角点的二维平面坐标；Z为被测量区域的地面平均高度；当前所处理航空摄影图像的旋转矩阵对所述纠正图像的格网范围进行确定时，根据公式(2)，计算得出所述纠正图像的格网范围；公式(2)中Xmin和Xmax分别为四个角点的地面坐标中横坐标的最小值和最大值，Yniin和Yniax分别为四个角点的地面坐标中纵坐标的最小值和最大值；gridsize为预先设定的所述正射影像的分辨率，所述正射影像的大小为mXn个像素点，所述纠正图像上包括mXn个格网点；其中，m和η分别为所述纠正图像中格网点的列数和行数； 步骤三、拉花检测:结合步骤一中所述的原始测量数据，对当前所处理纠正图像上的各格网点分别进行拉花检测，并且各格网点的拉花检测方法均相同；对任一个格网点进行拉花检测时，过程如下: 步骤3 O 1、 摄影光线方位确定:根据公式,对当前检测点的投影光线方位进行确定；其中，当前检测点为当前状态下进行拉花检测的格网点；公式(3)中，α为当前检测点的投影光线相对于投影中心的方位角，所述投影中心为当前所处理航空摄影图像的摄影中心点在地面上的投影点；(Xp，Yp)为当前检测点的地面坐标； 步骤302、拉花判断:判断当前所处理纠正图像中沿当前检测点的投影光线上，是否存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点:当判断得出存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点时，说明当前检测点为“拉花”点；反之，当判断得出不存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点时，说明当前检测点为未“拉花”点；本文档来自技高网...

【技术保护点】
一种数字正射影像拉花区域的处理方法，其特征在于该方法包括以下步骤：步骤一、获取被测量区域的原始测量数据：对被测量区域进行航空摄影测量，摄取被测量区域的航空摄影图像，并获得被测量区域的原始测量数据；所述原始测量数据为DEM测量数据，所述DEM测量数据为数字高程数据；步骤二、正射影像的地面格网范围确定：先根据步骤一中所述原始测量数据，获取步骤一中当前所处理航空摄影图像的纠正图像，所述纠正图像为根据当前所处理航空摄影图像在地面上的投影范围和需生成正射影像的地面分辨率进行划分所形成的格网图像，所述格网图像的大小与当前所处理航空摄影图像的正射影像的大小相同且所述格网图像上的各格网点分别与所述正射影像上的各像素点一一对应，所述正射影像为对所述航空摄影图像进行数字微分纠正和重采样后获得的影像，所述航空摄影图像为所述正射影像的原始影像；之后，根据当前所处理航空摄影图像的外方位元素和内方位元素，计算得出当前所处理航空摄影图像的四个角点的地面坐标；然后，根据计算得出的四个角点的地面坐标，确定所述纠正图像的格网范围；当前所处理航空摄影图像的四个角点分别为点A、点B、点C和点D，点A、点B、点C和点D的地面坐标分别为(X1，Y1)、(X2，Y2)、(X3，Y3)和(X4，Y4)；根据公式Xi=Xs+(Z-Zs)a1×xi+a2×yi-a3×fc1×xi+c2×yi-c3×fYi=Ys+(Z-Zs)b1×xi+b2×yi-b3×fc1×xi+c2×yi-c3×f---(1),]]>计算得出四个角点的地面坐标；公式(1)中，i为正整数且i＝1、2、3、4；(XS，YS，ZS)为当前所处理航空摄影图像的所述外方位元素中的摄影中心点坐标，f为所述内方位元素中的一个参数且其为步骤一中对被测量区域进行航空摄影测量时所用航摄仪的焦距；(xi,yi)为当前所处理航空摄影图像上四个角点的二维平面坐标；Z为被测量区域的地面平均高度；当前所处理航空摄影图像的旋转矩阵R=a1a2a3b1b2b3c1c2c3;]]>对所述纠正图像的格网范围进行确定时，根据公式m=Xmax-Xmingridsizen=Ymax-Ymingridsize]]>(2)，计算得出所述纠正图像的格网范围；公式(2)中Xmin和Xmax分别为四个角点的地面坐标中横坐标的最小值和最大值，Ymin和Ymax分别为四个角点的地面坐标中纵坐标的最小值和最大值；gridsize为预先设定的所述正射影像的分辨率，所述正射影像的大小为m×n个像素点，所述纠正图像上包括m×n个格网点；其中，m和n分别为所述纠正图像中格网点的列数和行数；步骤三、拉花检测：结合步骤一中所述的原始测量数据，对当前所处理纠正图像上的各格网点分别进行拉花检测，并且各格网点的拉花检测方法均相同；对任一个格网点进行拉花检测时，过程如下：步骤301、摄影光线方位确定：根据公式α=arctg(Ys-YpXs-Xp)Xs-Xp>0α=π+arctg(Ys-YpXs-Xp)Xs-Xp<0---(3),]]>对当前检测点的投影光线方位进行确定；其中，当前检测点为当前状态下进行拉花检测的格网点；公式(3)中，α为当前检测点的投影光线相对于投影中心的方位角，所述投影中心为当前所处理航空摄影图像的摄影中心点在地面上的投影点；(Xp,Yp)为当前检测点的地面坐标；步骤302、拉花判断：判断当前所处理纠正图像中沿当前检测点的投影光线上，是否存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点：当判断得出存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点时，说明当前检测点为“拉花”点；反之，当判断得出不存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点时，说明当前检测点为未“拉花”点；当前检测点的投影光线为步骤301中所确定当前检测点的摄影光线在所述纠正图像上的投影线；步骤303、多次重复步骤301至步骤302，直至完成当前所处理纠正图像上所有格网点的拉花检测过程；步骤四、影像纠正：根据步骤三中的拉花检测结果，对当前所处理航空摄影图像的纠正图像进行纠正，获得纠正后的正射影像；对当前所处理航空摄影图像的纠正图像进行纠正时，根据步骤三中对各格网点的拉花检测结果，对所述纠正图像上的各像素点分别进行数字微分纠正，并对各像素点的灰度值进行重采样。...

【技术特征摘要】

【专利技术属性】
技术研发人员：刘敏，郭永春，梁菲，姚春雨，左涛，
申请(专利权)人：西安煤航信息产业有限公司，
类型：发明
国别省市：陕西;61