您的当前位置:首页正文

基于局部Hough变换的机场跑道检测方法

2024-09-16 来源:欧得旅游网
2017年 6月 装备学院学报 June 2O17 Vol_28 No.3 第28卷第3期 Journal of Equipment Academy 基于局部Hough变换的机场跑道检测方法 焦 姣 , 巩向武 , 郭皇皇 (1.装备学院研究生管理大队,北京101416; 2.装备学院复杂电子系统仿真实验室,北京101416) 摘 要 针对可见光遥感图像中的机场跑道检测问题,提出了一种基于方 向预判的局部Hough变换方法。通过频谱能量楔向采样,确定频谱图中对能量贡献 起主导作用的系数分布角度范围与机场跑道的对应关系,得到跑道的方向预判结果; 利用方向信息对Hough变换角度的取值范围进行压缩,通过局部Hough变换进行 直线检测得到机场跑道。实验结果表明:该方法能够较好地消除复杂背景中地物的 干扰,快速准确地检测出机场跑道。 关 键词 频谱分析;边缘检测;局部Hough变换;跑道检测 TP751 A 中图分类号文献标志码文章编号 2095—3828(2017)03—0032—06 DOI 10.3783/j.issn.2095—3828.2017.03.006 Airport Runway Detection Based on Local Hough Transform JIAO Jiao , GONG Xiangwu。, GUO Huanghuang (1.Department of Graduate Management,Equipment Academy,Beijing 10141 6,China; 2.Science and Technology on Complex Electronic System Simulation Laboratory,Equipment Academy,Beijing 101416,China) Abstract In this paper,a local Hough transform method based on direction prediction is pro— posed for the airport runway detection of visible remote sensing image.The runway direction predic— tion is achieved by determining the correspondence between the coefficient distribution angle range with a dominant function in energy contribution and the airport runway.With the use of the direction information of the runway,Hough transformation range is reduced,and the airport runway is ob— tained via the line detection through local Hough transformation.Experiments show that,with this method,the disturbance from the ground obj ects in complex background can be better removed,and the airport runway can be detected with great speed and accuracy. Keywords spectrum analysis;edge detection;local Hough transformation;airport runway de tection 遥感图像中的机场跑道具有较为明显的特 取出跑道。文献E2]首先利用形态学方法进行图 征:跑道通常表现为平直且边沿较为齐整的带状, 长宽具有一定的范围,边缘线一般为平行线对;跑 像增强,然后采用最大类间方差法进行阈值分割, 提取得到机场连通域,经过细化之后通过2次 Hough变换提取跑道。文献[3]根据机场跑道在 道区域的灰度级通常高于附近背景区域,具有与 周围环境较为明显的差异。这些特征可应用于跑 道检测与提取。基于此,许多学者进行了研究。 各波段所具有的不同灰度特征及其直线特征进行 跑道定位,再对跑道关键点进行区域生长。文献 [4IN用最大连通域提取得到机场区域,然后通过 边缘提取和Hough变换确定跑道。这些特征在 复杂背景下存在的损失容易造成目标漏警和虚 文献Eli首先对图像进行分割,然后利用各向异性 扩散和Frangi滤波去除噪声干扰,检测得到可能 的跑道信息,最后利用目标的形状和颜色特征提 收稿日期 2O17一O4一l7 基金项目 部委级资助项目 作者简介焦姣(1988),女,博士研究生,主要研究方向为空间信息处理和遥感图像处理。jiaojiao nk@163.corn 笫3期 熊 姣。等:基于肘部H(mgh变换的机场跑道检测方法 :{:{ 警,而传统Hough变换的计算量也非常大。 本文首先通过频域能量分布分析确定机场跑 道的大致方向,然后将Hough变换中的角度取值 限定在较小的范围内,减少了Hough变换的计算 量,同时准确枪测出_rl机场跑道。l本文提出的方 法包括以下几个步骤:(1)采用改进的中值滤波 策略对图像进行降噪处理;(2)通过阈值分割得 到二值图像;(3)采用Canny算子进行边缘榆测; (4)将图像变换至频域通过能量谱楔状采样峰值 分析得到跑道的大致方向;(5)刈’边缘图像进行 局部Hough变换得到机场跑道的检测结果。整 个流程如图1所示。 l 机场跑道检测流秤 l 图像预处理 1.1图像降噪 遥感图像在拍摄的过程 }J难免会存在一些噪 声,一般是突发噪声点。为了更好地提取日标,首 先需要对图像进行降噪处理。目前.图像降噪有 许多成熟的滤波器,其中中值滤波对随机噪声具 有非常理想的降噪能力,但其效果依赖于滤波窗 口大小——窗U太大会导致边缘模糊,太小则去 噪效果不理想。本文采用一种改进的中值滤波策 略,即对图像进行扫描,判断像素点是否是滤波窗 口覆盖下邻域像素的极小值或者极大值。如果 是,依然采用普通中值滤波进行处理;如果不足, 则不予处理。由于噪声点基本郁属于邻域极值. 利川改进的方法能够在几乎不对边缘造成影响的 情况下有效去除突发噪声点,尤其是消除椒盐 噪声。 1.2阈值分割 一幅图像可以看作由前景与背景构成,可以 通过统计学方法获取分割阈值,使前景与背景尽 可能分离。一般而言,机场具有与背景较为明显 的灰度差异,可以采用一种自适应闽值分割方法 对遥感图像进行二值化操作,即完成机场 背景 的分割。本文采用()tsu方法。。。()tsu阈值分割 方法又称最大类问方差法,该方法计算简单, 一 定条件下不受图像亮度与对比度变化的影响。 假设待分割图像仃L个灰度级 (0,l,…,L一1),通过阈值,将图像分为(、 和 (、: 2类,即具有肤度级{0。l,…,了、}的像索分为 (、 ,{7、+1,了 _厂2,…, 一1}的像素分为(、 设 ( )表示 值为f时的类问方差,通过求最 大值叮得到最 阈值1’。 T、一arg Illax{ l{( ) (1) l§(,) _- l(,)(, l(f)一川 ) +/J z( )(7, 2(f)一,, g) (2) 式【{1:,,j ( ) , ( )分别是(、 和(、 的像素平均 值; ( ),,) (,)分别是柑应的概率密度; 是 全局像素均值。 如图2为耿度图像的直方图,使用()tsu法求 得丁一0.5/19,转换在[0,255]之问为1 40.0.闽值 丁位于2个峰值之问的波谷位置。 ● 瓣 5 4 3 2 l 啪 0 _l cI  l魏g霹 壤 h f j ih lJ 1: . 0 50 l00 150 200 250 灰度级 J割2 版度l习1黎rl:力I冬I 1.3边缘检测 Canny足常川的边缘检测算子,它能够在噪 卢抑制和边缘精确定位之 取得最佳折 f 。本文 采用Canny算子进行边缘检测,其原理是:在 、 y方向求一阶导数,将其组合为q个方向的导数, 这些方向的导数达到局部最大值的点就足构成边 缘的候选像素点。具体算法如下“:(1)用高斯 滤波器对图像进行平滑滤波;(2)用一阶偏导有 限差分来计算梯度的幅值l不1l方 ;(3)对梯度幅 值进行非极大值抑制;(4)刷双阈值算法对边缘 进行检测和连接。 遥感图像背景通常较为复杂,存在多种地物。 完成 像分割操作之后依然存在背景残留,导致 边缘提取结果【_}l存在部分背景连通体边缘。【六I 此,本文对边缘提取结果进行简单的连通体筛选. 滤除较小的连通体边缘。 2基于频谱能量分析的方向预判 图像的傅● 叶变换是将I兜l像从审问域变换到 频率域,也就是将图像的光谱能量分布./’( ,j,) 变换为频率能 分布F(“, ),频谱图能够很好 装 备学院学报 地反映 不同的频率成分,【_i土反映了 像能垣随 步页率的变化情况 。在图像的傅里叶频谱中,亮 点的叫峭表 了陔点与邻域点差异强弱,亮度强 的点梯度值较人。『』I】果频潜【冬1中亮度强的点较 多. 像通常较为尖锐,具有较为剧烈的灰度变化 离中心的区域属高频部分.与原 像r{ 灰度级变 化较快的分量奉}{对应,常衷征图像的细 ,如边缘 秆1噪声等。 2.1频谱能量分析 将图像进行二维傅咀ll1‘变换, I 3所示 。 趋势,边 分明并I{.边界处像素差异较大。通常 图像t 仃在较多的水平与垂直办向线条, 此频 谱 ll 点多集・} 丁水平 乖 方向。原 像中 傅里口1‘变换得到的频谱数摧并非以零频为【fI心进 行排列,而是按照原始的讣算川贝序. 此零频位于 频谱图角上。利用频潜周期性特点.通过平移可 的水平边缘 频谱 中的 直亮线相刘‘应。『『ii竖 以将零频移至图像中心位置。低频分}I{=从频谱 直边缘!』{Ij与水平亮线相对应。频潜图中靠近巾心 r『1. 向四周扩散;同时.高频分垃!J!lJ分散在四周. 的 域为低频部分,对应原图像的慢变化分艟;远 彳丁利于财图像特征进行分析。 亩流中 I 直流中心 : 维图像 高频成分 低频成分 3 l矧像频潸坐标移位示惑 2.2跑道方向预判 分划分为6个小区域.分别做6个 的能{}}统 遥感图像lII.各类地物特征线条表现出多 计,划分角度为:[6()。。9()。) ㈤。.6O。).[o。,30。). 方 .且随机分 于图像『f】符个位置。根据傅 [一30。,0。).[ 60。.一3O。).1 9()。.一60。).&IJ 里叶能f1谱的白眦准性,具有ll{=卡}1同方向的特征线 4所示。 不受位 约束.其能量贡献会替置在一起,并在频 0。 jU 30。 IIJ集中表征为通过中心且垂直于原特征线方 IP]的条带。若原图像中存在大 的某一角度的边 ∞ 缘线,则其频谱能量的楔状采样在 +丌 2 有较 大的峰ffIf。遥感 像中机场线条较为 一,分布 在 像『fI的某一个位置,通过分析楔向能量分布 的峰值化置即叮大敛确定图像巾机场的方向。 采川极坐标系进行频谱能 E( , )分析, 似没(r.0)为( .-、,)平面上的极坐标 l 能 i}}楔状采样 r f。 F,(0)一f F( . )dl 一I 『F( 。 )l。dr 3基于方向信息的Hough变换 J¨ √l l(3) 3.1 Hough变换基本原理 迎过积分町以得到各扇形 域中的频潜能量 Hough变换是一种常用的提取 线的 法. 分布。离散二维 像可用求和米替代积分形式。 实质上是一种投票机制,埘参数 M『lfl的点进仃 在相J、 n勺扇 埘, 求和 投票,投票值超过阈值则认为有足够多的点化于 陔参数点所决定的直线l 。该方法具仃受亩线是 ,( ,)一>:F,(0), 一1,2....川 (4) 11 存在断点和噪声影响小的优点。经过Hough 考虑经过平移的频谱图[fj频谱能量系数关于 变换,可以得到直线长度、端点、夹角.以及直线 原点『 I刈称分布,只收上半部分的频谱系数进行 垂直距离等参数。Hough变换的 本思想是点 楔 分 能量采样.根据对称性将频谱 卜半部 线的对偶性.它把图像 中直线的榆测问题变 第3期 焦 姣,等:基于局部Hough变换的机场跑道检测方法 35 换到参数空间中点的检测问题。变换关系为 D COS +ysin (5) 由于背景粘连、成像条件等原因容易出现边缘断 裂。本文对Hough变换检测出的直线进行扫描, 式中:p为坐标原点到直线的距离; 为直线的法 线与 轴的夹角。图像空间r.y中的每个点与参 数空间p-O中的一条正弦/余弦曲线相对应,且同 一根据平行线对距离和局部灰度分布信息对部分短 线段进行区域生长,得到较为完整的跑道。 直线上的点所对应参数空间中的曲线均交于一 4实验结果及分析 如图6所示,选取多幅机场遥感图像进行实 验验证,统一图像大小为800×800。实验硬件环 境为Intel Core i5—4200U,CPU主频1.60GHz,软 点。将参数空间离散化,统计过参数空间每个点 [fD, ]的曲线条数, 值相同则表征2条平行的 直线。如图5所示,空间中有3个点,分别对各点 求过该点的N个方向直线的 , ]坐标。通常 N—l8O,对应的检测角度精度为l。。图中3个点 件环境为MATI AB R2014a。这3处机场具有 较好的代表性和差异性,机场背景主要以农田和 村庄为主:图6a)中有一条跑道,跑道一侧有一条 滑行道;图6b)中机场较为复杂,有2条平行跑 道;图6c)中有2条“X”形跑道。 的[10, ]曲线相交于点[1D ,],该交点即表征 经过这3个点的直线。Hough变换直线检测对 经过同一点的曲线数进行累加,超过一定阈值的 [J0, ]坐标,即可判定此处存在一条直线。 a)北京两苑机场 b)参数空间 图5 Hough变换原理图 b)台湾桃园机场 3.2 跑道检测 遥感图像中的跑道一般为具有一定宽度的平 行线对,且具有显著的方向性,与其他地物目标有 明显区别 l。根据跑道边缘线特征,结合方向预 判信息,对边缘检测结果进行局部Hough变换。 传统Hough变换参数臼的取值范围是[一9O。, 90。),经过边缘检测提取到的边缘像素点都要参 与运算,具有计算量大的缺点。本文方法利用频 谱能量分布分析得到机场跑道方向信息,将 的 取值范围压缩到跑道方向邻域内的局部小区间 [ ,臼。)内,将对 的全局投票简化为局部投票。 能有效减小计算量并准确检测到机场跑道。 机场跑道边缘一般为基本等长的平行直线, 4.1 跑道方向预判结果及分析 c)芷江机场 图6机场遥感原图像 首先针对原始图像,进行图像降噪以及阈值 装 衔学院学报 分割;然后采川('army算子进行边缘检测,滤掉 为楔向能皱累计峰值,与机场跑道方向的角度埘 过小的连通体;llIf通过傅 叶变换将降噪处理后 血关系满足 { 2 的灰度 像, 换仝频域。提取 嘴 中对能 贯 根据能 分布检测剑I殳I 6; )对应的能 统汁 』{jJ:起t 作川的系数分 角 范m.得到跑道的 图中条柱l 峰值。卡H应角发范围为L㈤ -90 )・ ,J‘ 预划结 ..刈‘ 6中3 I 像分别进行处 则机场跑道力‘ 大敛在[一3O .0。);图 )对J、¨内 ,得剑 f1ll_分荆结果、边缘榆测结果、对心的频 能量统训 叶1榆测到条柱5为峰值,卡}1应们度 频 6个 向的楔 能 统汁结果以及 罔为[一6()。,一30。),则跑道力‘向大致仡l: )”, 、机场跑道力‘ 预州结果. 表1J听,Ji。表1叫|第 60。); 6c)跑道为“X”形.埘心的能 统汁 Il 1)q列为条仆 .条佳l~6分圳对嘘[6() ・90 )、 检测到条f 2为峰值.条 3能量累计值 j条 t 6(J。)、『()。. 。)、l一 。.0。)、 一6()。.一30。)、 2相筹0.: .卡Ij应的角发范 分圳为i .㈤”)、 1一㈨。,一60。)6个角度范 。除 频谱图巾的水 [so。。6O。)。!J!IJ跑道方向分圳 [一㈨ 。一6() )、 与垂卣凫线之外,一般还仔 通过频谱 心的 r 6()。, .{(】。).据此l1』以通过 值设置检测 划 斜IhJ亮色条带,j 仃突 的能 集-II.该条带表现 双峰值。提取“X”形机场跑道。 表I 图6a)~图6c)机场跑道方向预判结果  l,I。I 阀 ‘}削 粜 边缘俭洲 顷 颂 【 l恻I lj能:i}境il 机场跑i盟,J…J ■l,r ● -■I-■E _ 暖瞄臻鼹翳鞣勰 ) 2 I 翟 暖 霍I I  _ _-●■■。。‘。一 2未3. l4 50   I !拦蚵1滕  .5  _删 】.1I 1 1 2 3 的局部Hough变换和传统Hough变换方法进fj 僻ylj厅…颅削信息后 将 6中3幅 像的 跑道榆测,结 如 7~ )所示。跑道愉洲2f 边缘 像作为输入.分别利用 __r方向预削信息 果以红色线条作为标记。 )Hough变换食 投票 c{)lt()ugh 笾换局部擞 7 )机场跑道卡=;=测结 第3期 焦 姣.等:丛于局部Hougl 变换的机场跑道榆测力‘法 全局}tough变换检测结粜 b)水文方法检测结果 c、)Hough变换全局投票d)Hm g1、变换^ 邾投票 图8 )机场跑道检测结果 臼 )术史方法榆测结果 c、)Hough变换伞硒投票 d)Hough变换硒 投燥 9 6c)机场跑道卡令测结果 图7~图9中的图像a)分别为图6中机场 机场跑道方向的对应关系。将Hough变换f}1全局 投票简化为局部投票,大大减少了Hough变换的 计算量,快速准确地检测出机场跑道。下一步工 像的全局Hough变换检测结果,图像b)分别为 局部Hough变换的检测结果,图像c)和图像d) 分别为Hough变换全局和局部投票结果。从上 述结果u丁以看出:图6a)中山于全局变换的缺点, 容易榆测出虚假直线;罔6b)中机场跑道结构较 为复杂,缩小检测范同可以得到更为精细的跑道 检测结果;图6c)中机场存在“X”形跑道,全局检 作需要研究自适应闽值分割技术,消除复杂背景 干扰,进一步提高跑道的检测率和算法的鲁棒性。 参考文献 (References) L 1 FRI1’A rHI A K.SWARUP S.Sba1) and color{eattlYeS 测仪能检测到一个方向的跑道。而局部检测可以 检测剑“X”形跑道。下面对Hough变换在时问 效率f 的改进进行验证。分别对局部与全局 Hough变换跑道检测所用时『开】、预处理时问以及 based airport rtlnway detec'tion L(、j//IEEE Im{、rnational Advance('om1)uling(’onference(1AL、(、).( baziabad.India: 20l 3:836 8 1 1. [2]董书驹.吴巍.冲p遥感罔像巾机场跑道的提取力法 -lj.武汉 理工夫学学报(信息与管理【程版). )06.28(7):l 2—1 55. [3] 博.哭秀清.刨拥军.统汁的}tougb变换在机场遥感图像分 割巾的廊用[J].汁算机12程,20{)2,28(8):261—265. 总耗时进行统计。结果如表2所示。 表2跑道检测耗时测试结果 (图像尺寸800×800 pixc1.分辨率要求0.1。) [4]陈旭光.林卉.遥感 像巾机场日标的识别方法 I_.、}’算机l一 程与应用.2()1 2.18(25):l 94 l97. 5]阳波.摹于R大炎问力‘差遗传算法的『 I像分割力法[J]。湖南 师范大学自然科学学报.2003.26(1):32—36. :6]IxJ萨需斯.数字l 像处理[M].3版.阮秋琦.译.北寐:电子工 分析表2测试结果可知,局部Hough变换引 入的预处理开销略大于全局变换,但是通过参数 业f“版社.2()J 2:2粕 1. [7]吴桂平,肖鹏峰.冯学钾.等.利用频谱能 进行高分辨率遥 感图像地物 {驯_J].武汉大学学报(信息利学版).2()l1.36 (11):1 29 I 1 2 )7. 范隔压缩,大大降低_r Hougt 变换本身的计算 量.使得跑道检测时间明显优于全局Hough变 换;同时,综合性能也较全局变换更优。综上所 述.本文方法可快速准确地检测出机场跑道。 [8]吴桂平.高分辨牢遥感 像频渚能鲢分析 典7 地物特征识 别研究[I)].}_}j : 求火学.2Oll:lO j6. :9]艾淑芳.闩钧华,李夫甫.等.遥感图像巾的机场跑道检测算 法[J].电光 j控制.20l 7(2):13 l6. 5 结束语 本文通过频域能量楔向采样确定能量峰值与 (编辑:李江涛) 

因篇幅问题不能全部显示,请点此查看更多更全内容