一种基于线段分布的机场线性特征提取方法与流程

1.本发明属于遥感图像处理技术领域，具体涉及一种基于线段分布的机场线性特征提取方法。


背景技术：

2.随着遥感技术的快速发展，遥感影像的空间和时间分辨率越来越高，使得遥感的实效性和实用性大大提高，被广泛地应用于各个领域。机场作为一种重要的基础设施，利用遥感图像识别其位置和轮廓信息具有非常重要的研究意义。
3.机场在高分辨率遥感影像上具有较为特殊光谱和纹理特征，通常情况下机场由大面积规则不透水层组成，导致其灰度值比周围地物大，且由多条跑道组成，直线特性明显。许多研究基于机场亮度特征和直线特征，重点利用显著性和线段联合构建模型提取机场，取得一定可观的效果。虽然利用机场跑道直线和平行的特点可以实现快速的目标提取，但其对直线特征利用往往不是很充分，导致许多具有线段的非机场目标也会被检测出来。因此，如何在高分辨率遥感影像上提取有效的机场线段特征成为一个亟待解决的研究方向。


技术实现要素：

4.有鉴于此，本发明提供一种基于线段分布的机场线性特征提取方法，其可利用线段的空间位置及方向特性筛选出机场所属线段，为机场识别提供可靠的特征层支撑。
5.为了实现上述目的，本发明采用的技术方案为：
6.一种基于线段分布的机场线性特征提取方法，包括以下步骤：
7.步骤1，对多光谱遥感影像进行处理，提取可见光的三个波段，并进行线性拉伸；
8.步骤2，采用线段检测器lsd在遥感影像上提取线段，得到线段集合l
set
；
9.步骤3，依据各线段的端点坐标计算线段集合l
set
中每一条线段的长度li和倾角θi；
10.步骤4，以线段倾角θi作为依据，对线段集合l
set
进行分组，形成线段子集
11.步骤5，将线段子集中的线段抽象为具有质量的质点pi，并计算所有质点构成质点系的质心cmk；
12.步骤6，计算线段子集中质点pi到质心cmk的距离，并剔除显著离群线段，得到更新后的线段子集
13.步骤7，对更新后的线段子集计算线段长度的变异系数
14.步骤8，先计算更新后的线段子集的质心scmk，再计算线段质心平均距离
15.步骤9，计算更新后的线段子集的机场隶属度pk：
[0016][0017]
其中，为加权调节因子，nor()表示对参数值进行最大值最小值归一化处理；
[0018]
步骤10，选择机场隶属度pk最大的线段子集作为最终的机场区域线段l
ap
。
[0019]
进一步地，步骤1中，对波段bi进行线性拉伸的计算方式如下：
[0020][0021]
式中，表示多光谱遥感图像第i个拉伸后的波段，和分别表示波段bi的最小灰度值和最大灰度值。
[0022]
进一步地，步骤3的具体方式为：
[0023]
获取线段linei的两个端点坐标和
[0024]
每条线段的长度li计算方式如下：
[0025][0026]
每条线段的倾角θi计算方式如下：
[0027][0028][0029]
其中，ki表示线段linei的斜率，arctan( )表示反正切函数，θi以弧度的形式表示。
[0030]
进一步地，步骤4中，分组的方式如下：
[0031]
设定倾角间隔θ
gap
，θ
gap
能够被180整除；将0-180度范围划分为180/θ
gap
个区间，将倾角在同一区间内的线段划为一组；即：
[0032][0033]
θ
gap
×
(k-1)＜θj≤θ
gap
×
k and linej∈l
set
[0034]
k＝1,2,3,...,180/θ
gap
[0035]
其中，θj为线段linej的倾角。
[0036]
进一步地，步骤5中，以线段linei的中点坐标为质点pi的位置，以线段linei的长度为质点pi的质量；
[0037]
中所有质点构成质点系的质心cmk坐标为：
[0038][0039][0040]
其中，和分别为质心cmk的横纵坐标。
[0041]
进一步地，步骤6中，线段子集中每个质点pi到质心cmk的距离为：
[0042][0043]
其中，d(pi,cmk)表示求质点pi到质心cmk的距离；
[0044]
设线段子集中有n条线段，则依据长度li对所有线段进行排序，并获取长度上
下四分位数和则更新后的线段子集为：
[0045][0046][0047]
其中，m为更新后的线段子集包含的线段总数。
[0048]
进一步地，步骤7中，更新后的线段子集中线段长度的变异系数为：
[0049][0050]
其中，和分别表示中线段长度的均值和标准差。
[0051]
进一步地，步骤8中，更新后的线段子集的线段质心平均距离计算方式如下：
[0052][0053]
其中，d(pi,scmk)表示线段质点pi到更新后的线段子集的质心scmk的距离。
[0054]
本发明具有如下有益效果：
[0055]
(1)本发明方法基于线段分布提取机场的线性特征，可以在不需要训练样本的情况下，有效提取遥感图像中的机场线段特性。
[0056]
(2)应用本发明方法可以有效提取线段的空间分布、长度和方向特征，为机场检测提供关键的线段特性描述。
附图说明
[0057]
图1为机场线段提取方法的原理示意图。
具体实施方式
[0058]
下面结合附图及具体实施方式对本发明做进一步详细的说明。
[0059]
如图1所示，一种基于线段分布的机场线性特征提取方法，包括以下步骤：
[0060]
步骤1，对多光谱遥感影像处理，提取可见光三个波段，并对多光谱各波段bi拉伸，计算方式如下：
[0061][0062]
式中，表示多光谱遥感图像第i个拉伸后的波段，和分别表示波段bi的最小灰度值和最大灰度值。
[0063]
步骤2，采用lsd算法在遥感影像上提取线段得到线段集合l
set
；
[0064]
步骤3，依据各线段的端点坐标计算线段集合l
set
中每一条线段的长度li和倾角θi。首先需要获取线段linei的两个端点坐标分别为和则每条线段的长度li计算方式如下：
[0065][0066]
每条线段的倾角θi计算方式如下：
[0067][0068][0069]
其中，ki表示线段linei的斜率，arctan()表示反正切函数，θi以弧度的形式表示。
[0070]
步骤4，以线段倾角θi作为依据，对线段集合l
set
进行分组，形成线段子集其分组规则如下：
[0071][0072]
θ
gap
×
(k-1)＜θj≤θ
gap
×
k and linej∈l
set
[0073]
k＝1,2,3,...,180/θ
gap
[0074]
其中，θj为linej的倾角，θ
gap
为线段倾角间隔。
[0075]
步骤5，将线段子集中线段linei抽象为具有质量的质点pi，以中点为质点位置，长度为质点的质量，其位置计算方式如下：
[0076][0077][0078]
其中，和分别为质点pi的横纵坐标。
[0079]
中所有线段质点构成质点系的质心cmk，其坐标计算方式如下；
[0080][0081][0082]
其中，和分别为质心cmk的横纵坐标。
[0083]
步骤6，计算线段子集中每个质点pi到质心cmk的距离为：
[0084][0085]
其中，d(pi,cmk)表示求质点pi到质心cmk的距离；
[0086]
设线段子集中有n条线段，则依据长度li对所有线段进行排序，并获取长度上下四分位数和利用长度筛选规则剔除显著离群线段，则更新后的新线段子集为：
[0087][0088][0089]
其中，m为新线段子集包含的线段总数。
[0090]
步骤7，对新线段子集计算线段长度的变异系数
[0091][0092]
其中，和分别表示中线段长度的均值和标准差。
[0093]
步骤8，先计算新线段子集的质心scmk，再计算线段质心平均距离
[0094][0095]
其中，m为新线段子集包含的线段总数，d(pi,scmk)表示线段质点pi到质心scmk的距离。
[0096]
步骤9，计算线段子集的机场隶属度pk：
[0097][0098]
其中，为加权调节因子，nor()表示对参数值进行最大值最小值归一化处理；
[0099]
归一化函数nor()为：
[0100]
nor(w)＝(w-minw)/(max
w-minw)
[0101]
其中，w为输入的待归一化的数组，minw和maxw分别表示输入数组w的最小值和最大值。
[0102]
步骤10，选择隶属度pk最大的线段子集作为最终的机场区域线段l
ap
：
[0103][0104]
其中，o表示机场最大隶属度的线段组编号。
[0105]
总之，本发明提出的机场线段特征提取算法，可以为高分辨率遥感影像的机场提取提供较好的线段特征，该方法无需选择训练样本，采用自动阈值选择方法，无需人工干预识别，为便捷、快速和精确的目标线段特征提取提供技术支撑。