[发明专利]数字正射影像拉花区域的处理方法

权利要求书

1.一种数字正射影像拉花区域的处理方法，其特征在于该方法包括以下步骤：

步骤一、获取被测量区域的原始测量数据：对被测量区域进行航空摄影测量，摄取被测量区域的航空摄影图像，并获得被测量区域的原始测量数据；所述原始测量数据为DEM测量数据，所述DEM测量数据为数字高程数据；

步骤二、正射影像的地面格网范围确定：先根据步骤一中所述原始测量数据，获取步骤一中当前所处理航空摄影图像的纠正图像，所述纠正图像为根据当前所处理航空摄影图像在地面上的投影范围和需生成正射影像的地面分辨率进行划分所形成的格网图像，所述格网图像的大小与当前所处理航空摄影图像的正射影像的大小相同且所述格网图像上的各格网点分别与所述正射影像上的各像素点一一对应，所述正射影像为对所述航空摄影图像进行数字微分纠正和重采样后获得的影像，所述航空摄影图像为所述正射影像的原始影像；之后，根据当前所处理航空摄影图像的外方位元素和内方位元素，计算得出当前所处理航空摄影图像的四个角点的地面坐标；然后，根据计算得出的四个角点的地面坐标，确定所述纠正图像的格网范围；

当前所处理航空摄影图像的四个角点分别为点A、点B、点C和点D，点A、点B、点C和点D的地面坐标分别为(X₁，Y₁)、(X₂，Y₂)、(X₃，Y₃)和(X₄，Y₄)；根据公式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；(X_S，Y_S，Z_S)为当前所处理航空摄影图像的所述外方位元素中的摄影中心点坐标，f为所述内方位元素中的一个参数且其为步骤一中对被测量区域进行航空摄影测量时所用航摄仪的焦距；(x_i,y_i)为当前所处理航空摄影图像上四个角点的二维平面坐标；Z为被测量区域的地面平均高度；

当前所处理航空摄影图像的旋转矩阵R=a1a2a3b1b2b3c1c2c3;]]>

对所述纠正图像的格网范围进行确定时，根据公式m=Xmax-Xmingridsizen=Ymax-Ymingridsize]]>(2)，计算得出所述纠正图像的格网范围；公式(2)中X_min和X_max分别为四个角点的地面坐标中横坐标的最小值和最大值，Y_min和Y_max分别为四个角点的地面坐标中纵坐标的最小值和最大值；gridsize为预先设定的所述正射影像的分辨率，所述正射影像的大小为m×n个像素点，所述纠正图像上包括m×n个格网点；其中，m和n分别为所述纠正图像中格网点的列数和行数；

步骤三、拉花检测：结合步骤一中所述的原始测量数据，对当前所处理纠正图像上的各格网点分别进行拉花检测，并且各格网点的拉花检测方法均相同；对任一个格网点进行拉花检测时，过程如下：

步骤301、摄影光线方位确定：根据公式α=arctg(Ys-YpXs-Xp)Xs-Xp>0α=π+arctg(Ys-YpXs-Xp)Xs-Xp<0---(3),]]>对当前检测点的投影光线方位进行确定；其中，当前检测点为当前状态下进行拉花检测的格网点；公式(3)中，α为当前检测点的投影光线相对于投影中心的方位角，所述投影中心为当前所处理航空摄影图像的摄影中心点在地面上的投影点；(X_p,Y_p)为当前检测点的地面坐标；

步骤302、拉花判断：判断当前所处理纠正图像中沿当前检测点的投影光线上，是否存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点：当判断得出存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点时，说明当前检测点为“拉花”点；反之，当判断得出不存在其它与当前检测点在所述原始影像上对应的像素点相同的格网点时，说明当前检测点为未“拉花”点；

当前检测点的投影光线为步骤301中所确定当前检测点的摄影光线在所述纠正图像上的投影线；

步骤303、多次重复步骤301至步骤302，直至完成当前所处理纠正图像上所有格网点的拉花检测过程；

步骤四、影像纠正：根据步骤三中的拉花检测结果，对当前所处理航空摄影图像的纠正图像进行纠正，获得纠正后的正射影像；

对当前所处理航空摄影图像的纠正图像进行纠正时，根据步骤三中对各格网点的拉花检测结果，对所述纠正图像上的各像素点分别进行数字微分纠正，并对各像素点的灰度值进行重采样。