结合模糊分割和遗传算法的航空影像车辆提取方法

赵泉华，常 波，王 玉，李 玉

(辽宁工程技术大学测绘与地理科学学院遥感科学与应用研究所，辽宁 阜新 123000 )

结合模糊分割和遗传算法的航空影像车辆提取方法

赵泉华，常 波，王 玉，李 玉

(辽宁工程技术大学测绘与地理科学学院遥感科学与应用研究所，辽宁 阜新 123000 )

利用高分辨率遥感影像实现道路车辆监测在交通管理系统中占有重要地位，然而受成像条件及车辆本身颜色多样性的影响，很难从遥感影像中直接提取道路车辆。为此，本文提出了一种结合遗传算法和数学形态学的遥感影像道路车辆提取方法。首先利用模糊聚类方法建立影像分割模型，并利用遗传算法对其进行求解，从而实现影像中道路、车辆及环境的精确分割；然后对分割结果中不同层进行赋值,实现分割结果的二值化;最后利用数学形态学操作从分割的道路—车辆的二值结果中剔除道路，从而实现道路车辆的提取。利用本文方法对航空影像进行了道路车辆提取试验。试验结果表明，本文方法不仅能够准确识别地面车辆位置，还能精确提取其轮廓线。

目标提取；车辆边缘；遗传算法；航摄影像；模糊分割

作为图像处理的重要内容，目标提取已有广泛的应用，如房屋[1]、道路[2-3]、植被[4]等的提取。道路车辆提取作为目标提取中的重要任务之一，在监测道路车流量、车流密度等交通管理方面具有独特的潜在应用价值。近年来，随着无人机航测技术的发展，获取高分辨率遥感影像的成本越来越低，相比于利用固定放置的传统道路信息采集系统获取交通信息，利用遥感手段的交通信息获取方式具有更好的应用前景。目前，已有很多遥感道路车辆提取方法，主要分为基于显式[5]和隐式模型[6-8]两类。

针对小型车辆在遥感影像中细节信息不充分的问题，Jeon等[5]提出了一种基于显式模型的目标匹配方法，即选取影像的部分车辆作为模板，并设置匹配度，筛选出符合模板信息的区域标记为车辆，从而实现车辆提取。该方法虽能提取影像中车辆数量及其位置信息，但由于背景环境复杂，影响匹配结果，导致提取精度较低。

显式模型目标提取方法不能大范围提取车辆目标。为此，研究者提出通过检测车辆与背景环境间纹理及光谱特性的差别实现了目标提取的隐式模型方法。Stefan等[9]利用层次模型描述车辆在不同层次中的细节特征，同时结合车辆与其周围地物的邻域关系，对大范围城市遥感影像进行车辆提取。该方法虽然能提取车辆信息，但需要预先获取道路的光谱反射信息。

目前，还没有可精确提取高分辨率航摄影像道路车辆的有效方法。为此，本文提出一种结合遗传算法和数学形态学的目标提取方法。首先，采用遗传算法对航摄影像进行优化分割，利用其较强的鲁棒性弱化影像噪声对算法稳定性的影响，从而提高分割精度；然后，运用数学形态学方法处理其分割影像，平滑无关信息；最后利用一阶微分算子，精确提取道路车辆的边界轮廓。

1 算法描述

1.1 初始分割

首先对影像进行分割,设Z={zi(xi,yi),i=1,2，…,n}为彩色航空影像，其中，i为像素索引，zi=(zir,zig,zib)为像素i的彩色矢量，zir、zig、zib分别为zi的红、绿、蓝分量，(xi,yi)∈D为像素i的格点位置坐标，D为影像域，n为像素数。

为了刻画像素对目标类的类属性，其模糊隶属度定义为

(1)

式中，β为模糊因子，β越大表示算法模糊程度越大；μ(uj,zi)为聚类中心uj与zi的相似度

(2)

式中，t1和t2为给定常数；d(uj,zi)为向量uj和zi间的距离；θ(uj,zi)为两向量间夹角，表示如下

(3)

(4)

式中，l({r,g,b}为彩色分量索引。

为了实现最优模糊分割，需建立约束E和U的模糊目标函数

(5)

最小化式(5)，可求得最优模糊类属矩阵和聚类中心集，以实现图像的初始分割，即

(6)

1.2 基于遗传算法的目标函数解

1.2.1 构建初始种群

用遗传算法迭代求解目标函数的重点在于给定初始种群Po(0)。本文首先采用颜色直方图获取图像的候选聚类中心集，再由候选聚类中心集构建初始种群Po(0)[13]，如图1所示。

图1 颜色空间划分

首先获取图像的颜色直方图，将RGB颜色空间沿R、G、B轴离散化为N(N=NR×NG×NB)个颜色块bin，图1是离散化的颜色空间，每个bin在离散化空间中的位置为(R,G,B)。按照R→G→B的顺序对每个bin进行编号，则任意bin的索引k表示为

(7)

式中，R=0,1,…,NR-1；G=0,1,…,NG-1；B=0,1,…,NB-1，表示bin在离散化颜色空间的位置。则图像中值为zi(zir,zig,zib)的像素所在bin为

(8)

(9)

由此得到表征彩色图像颜色分布的离散直方图H。假定图像类别数为c，认为聚类中心存在于H中峰值最高的c个bin中，在这c个bin中分别随机选取一个彩色矢量作uj为聚类中心，构成聚类中心集U={uj:j=1,…,c}。

图2 c个聚类中心矢量组成的染色体Um

1.2.2 遗传操作

为了获得最优种群，对任意代种群Po(τ)需经遗传操作，演变出下一代种群Po(τ+1)。每次迭代遗传演变由3种遗传操作完成，即选择、交叉和变异操作。

(10)

1.2.3 图像分割

上述遗传操作经过T次迭代后，生成的种群为Po(T)。从Po(T)中选取适应度最大的染色体作为最佳聚类中心集

(11)

(12)

式中，Li为像素i的类属标号，而L={Li:i=1,2，…,n}表征图像Z的最优分割。

1.3 基于形态学的车辆提取

1.3.1 道路提取

将分割影像L中道路及其误分像素赋值为0，其他像素赋值为1，得到二值道路分割层影像L1。

为了恢复车辆在影像分割过程中所丢失的细节信息，对影像进行简单的二值形态学膨胀，使车辆恢复原始形态,再采用二值数学形态学中的开运算，消除车辆导致的道路分割层中几何噪声。由于车辆多为矩形目标，故设置方形结构元的数学形态学操作更符合车辆的形态特性。对L1进行膨胀运算得到车辆信息完整的影像L2,对其进行二值开运算，得到去噪二值影像L3，道路上的车辆被视为几何噪声被滤掉，保留完整的道路体，即所需的道路层影像。

1.3.2 车辆提取

道路车辆在分割结果的道路层影像L2中表现为几何噪声，因此L2包含了道路车辆的完备信息，为了提取道路车辆不必关心对其他地物类型的分割情况。L2消噪后道路层影像为L3，两者的差异仅仅在于车辆部分，因此将L2和L3相减即可获取车辆影像。

1.3.3 车辆轮廓线提取

由于获取的车辆层二值影像已经完成了对噪声的滤除，并且车辆与背景间具有明显颜色差异，故直接利用Roberts算子提取边缘[15-16]。

2 结果及分析

为了验证该算法的可行性及有效性，对航摄影像进行道路车辆提取试验。图3(a)为130×130像素、分辨率为0.15 m的城区道路的原始图像，人为给定聚类数c为2。设置参数t1、t2分别为经验值0.000 1和0.2，概率pc为0.5，变异概率pm为经验值0.2，最大迭代次数为30，种群大小K为40，利用颜色直方图优化的遗传算法对图3(a)进行图像分割，其结果如图3(b)所示。通过图3(b)可以看出，全部车辆可以与周围地物区分开，且分割精度较高。将分割结果进行二值化，得到二值影像，如图3(c)所示。在图3(c)中，除车辆和道路信息外，还存在部分斑点噪声。图3(d)为对二值影像进行膨胀处理的结果，消除了斑点噪声，并恢复了车辆的原始形态。图3(e)为将车辆看作噪声进行开运算处理的结果，可以看出，车辆在道路体上造成的空洞被成功填补且道路及周围背景都保留了的原有特性。将图3(d)与图3(e)相减，得到如图3(f)所示的二值影像，其中黑色区域为提取车辆，白色区域为无关背景。利用Roberts算子提取出车辆的轮廓线，结果如图3(g)所示，可以看出，车辆边界都被成功捕获且无遗漏。为对提出算法进行定性评价，将车辆轮廓(电子版中为红色线)叠加到原始影像(图3(a))上，结果如图3(h)所示，可以看出，绝大多数车辆轮廓与原图匹配成功，且能完美拟合。

为进一步验证提出算法的优越性，利用Canny边缘检测算法和ENVI软件中基于规则的目标提取方法对试验图像进行车辆目标提取对比试验，试验结果如图3(i)—(j)所示。两种对比算法在提取车辆轮廓的同时，并不能有效去除道路的影响，并且提取出的车辆轮廓线也不能很好地表述车辆形态。对3种算法提取的相应车辆信息进行定量精度分析，包括识别目标总数、识别车辆数、识别车辆比例、实际车辆数、车辆提取精度，并且以人工获取的车辆轮廓为模板图像，进行提取车辆轮廓的重合精度分析，其中重合精度是试验提取的车辆轮廓与模板中车辆轮廓重叠像素数与提取的像素总数的比值。该值不仅可以反映提取车辆的准确性，又考虑了无关地物的影响，能对提取精度进行充分的评价。结果见表1。通过表1可以看出，传统的Canny算法和基于规则的目标提取方法的轮廓重合率较低，且都识别了过多无关地物目标，导致识别车辆比例略低于本文算法；而利用本文算法不仅能识别所有车辆，并且较好地去除无关地物的影响，与模板重合率在85%以上，更能准确描述车辆的位置信息。通过对比，可以说明本文算法的优越性。

图3 试验影像

3 结 语

遥感影像作为交通管理的信息来源，具有观察范围广、观测周期短等优点。利用车辆信息的自动提取技术可以快速测量路段车流密度，获取高速路段出现拥堵或事故情况时的交通信息，以及根据车辆长宽比自动识别不同类型车辆。针对遥感影像的道路车辆提取问题，本文提出一种结合遗传算法和数学形态学的目标提取方法。相比于传统遗传算法随机生成初始种群的方法，本文利用影像的颜色直方图构建初始种群，以实现聚类中心优化，从而大大缩减了目标函数收敛的时间，提高了算法的效率，再采用数学形态学方法对车辆进行提取。该方法不仅能消除影像中的细小噪声，还弥补了影像分割所缺失的信息，有效且高精度地提取出车辆信息，为智能交通管理提供高效准确的数据来源。

表1 精度分析

[1] KHUMALO P P, TAPAMO J R. Rotation Invariant Texture Feature Algorithms for Urban Settlement Classification [J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 24(8): 511-514.

[2] 唐伟, 赵书河. 基于GVF和Snake模型的高分辨率遥感图像四元数空间道路提取[J]. 遥感学报, 2011,15(5): 1040-1052.

[3] 王科, 黄智, 钟志华. 基于多特征融合的道路理解方法[J]. 中国公路学报,2013, 26(4): 176-183.

[4] 杨胜天, 李茜, 刘昌明,等. 应用“北京一号”遥感数据计算官厅水库库滨带植被覆盖度[J]. 地理研究, 2006, 25(4): 570-578.

[5] JEON S H, LEE K, KWON B D. Application of Template Matching Method to Traffic Feature Detection Using KOMPSAT EOC Imagery [J]. IEEE International Geoscience and Remote Sensing Symposium, 2005, 180(2): 845-867.

[6] EIKVIL L, AURDAL L, KOREN H. Classification-based Vehicle Detection in High-resolution Satellite Images [J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2009, 64(1): 65-72.

[7] LARSEN S O, KOREN H, SOLBERG R. Traffic Monitoring Using Very High Resolution Satellite Images [J]. Photogrammetric Engineering and Remote Sensing, 2009, 7(75): 859-869.

[8] LEITLOFF J, HINZ S, SITLLA U. Vehicle Detection in Very High Resolution Satellite Images of City Areas [J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 7(48): 2795-2806.

[9] HINZ S, BAUMGARTHER A. Vehicle Detection in Aerial Images Using Generic Features, Grouping, and Context [M]. Berlin: Springer Berlin Heidelberg, 2001:45-52.

[10] BEZDEK J C. Pattern Recognition with Fuzzy Objective Function Algorithms [M]. New York: Plenum Press, 1981.

[11] POTTS J C， TERRI D. The Development and Evaluation of An Improved Genetic Algorithm Based on Migration and Artificial Selection [J]. IEEE Transactions on Systems, Man, and Cybernetics, 1994, 24(1): 73-86.

[12] 刘正军, 王长耀, 张继贤. 基于小波与遗传算法的特征提取与特征选择[J]. 遥感学报, 2005, 9(2): 176-185.

[13] 王向阳, 杨红颖, 郑宏亮,等. 基于视觉权值的分块颜色直方图图像检索算法[J]. 自动化学报, 2010, 36(10): 1489-1492.

[14] 邓素娟, 蒲石. 赌轮盘算法的实现及分析[J]. 乐山师范学院学报, 2014, 29(12): 63-82.

[15] 李晖, 肖鹏峰, 冯学智,等. 基于向量场模型的多光谱遥感图像多尺度边缘检测[J]. 测绘学报, 2012, 41(1): 100-107.

[16] ACCAME M, FRANCESCO G B. Edge Detection by Point Classification of Canny Filtered Images [J]. Signal Processing, 1997, 60(1): 11-22.

Fuzzy Segmentation and Genetic Algorithm Based Road Vehicle Extraction Method from High Resolution Aerial Image

ZHAO Quanhua,CHANG Bo,WANG Yu,LI Yu

(The Institute for Remote Sensing School of Geomatics，Liaoning Technical University, Fuxin 123000,China)

Nowadays, vehicle object extraction of remote sensing image and its potential applications become a hot topic. However, the optical information of vehicles is seriously weakened by a series of factors, such as environment and the velocity of the vehicles, so that the extraction accuracy of vehicles is accordingly reduced. To this end，a vehicle extraction method of remote sensing image using genetic algorithm and mathematical morphological operation has been proposed. First, a pre-clustering operation of remote sensing image is done by using histogram-based technique, and optimizing the initial cluster centers applied the genetic algorithm to improve the accuracy of segmentation results. Then, the mathematical morphological operation is used to extract the edges of vehicles. The aerial images are tested to extract vehicles using the proposed approach. The experiment results show that this method can identify the situation of vehicles on the ground as well as extract the contour line accurately.

object extraction; vehicle extraction; genetic algorithm; aerial images; fuzzy segmentation

赵泉华，常波，王玉,等.结合模糊分割和遗传算法的航空影像车辆提取方法[J].测绘通报，2017(8)：62-66.

10.13474/j.cnki.11-2246.2017.0255.

2016-12-20;

2017-03-03

国家自然科学基金(41301479)；大学生创新创业训练项目(201510147005)；辽宁省自然科学基金(2015020090)

赵泉华(1978—)，女，博士，副教授，主要研究方向为随机几何、空间统计学、模糊集理论等在遥感图像建模、解译及海洋环境遥感中的应用。E-mail: zqhlby@163.com

常 波。E-mail:953416733@qq.com

P237

A

0494-0911(2017)08-0062-05