基于萤火虫算法的农业遥感图像增强研究.pdf

《基于萤火虫算法的农业遥感图像增强研究.pdf》由会员分享，可在线阅读，更多相关《基于萤火虫算法的农业遥感图像增强研究.pdf（6页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、浙江农业学报Acta Agriculturae Zhejiangensis, 2016,28(6): 1076 -1081 http:/ / www. zjnyxb. cn郭红山,张慧宁.基于萤火虫算法的农业遥感图像增强研究J.浙江农业学报,2016,28(6): 1076 -1081.DOI: 10. 3969/ j. issn.1004-1524. 2016. 06. 27收稿日期:2015-10-12基金项目:河南省科技厅鉴定项目(预科鉴委字2014第664号)作者简介:郭红山(1978 ),男,河南宁陵人,讲师,硕士,从事计算机技术研究。 E-mail:ghs3171412 foxma

2、il. com基于萤火虫算法的农业遥感图像增强研究郭红山,张慧宁(黄河水利职业技术学院,河南开封475004)摘 要:农业遥感图像增强有利于图像信息的提取与分析,萤火虫算法是近年来较为新颖的智能仿生算法,目前国内外关于其能否用于农业遥感图像增强的研究未见报道。文章首先利用非完全Beta函数建立农业遥感图像增强模型,结合人眼最小灰度分辨力函数进行图像细节增强,将每个输入区间的像素灰度值变换到适当的输出灰度级区间,最终生成对比度均衡的图像;然后通过萤火虫优化算法在其动态决策域半径进行伪差分操作更新;最后确定最佳参数的收敛条件,给出了算法流程。试验仿真结果表明,萤火虫算法的农业遥感图像检测在图像细节

3、增强评价指标、相位一致性指标、通用质量评价指标等方面与直方图算法、Retinex算法、小波变换算法、模糊聚类算法等相比数据较优,能够用于农业遥感图像增强。关键词:萤火虫算法;Beta函数;农业遥感图像;增强;灰度中图分类号:TP393文献标志码:A文章编号:1004-1524(2016)06-1076-06Agricultural remote sensing image enhancement based on firefly algorithmGUO Hong-shan,ZHANG Hui-ning(Yellow River Conservancy Technical Institute,

4、 Kaifeng 475004, China)Abstract: Agricultural remote sensing image enhancement is advantageous to extraction and analysis of image infor-mation. Firefly algorithm is new intelligent bionic algorithm in recent years. Domestic and foreign research on wheth-er it can be used for agricultural remote sen

5、sing image enhancement has not been reported at present. Firstly, incom-plete Beta function was used to establish agricultural remote sensing image enhancement model, which was combinedwith human minimum gray level resolution function for image detail enhancement, each input interval of the pixelgra

6、y value was transformed to the appropriate output gray scale interval, and generated image contrast equalization.Secondly, difference operation was updated in the dynamic decision domain radius based on firefly algorithm. Final-ly, the convergence condition of the optimal parameter was determined, a

7、nd the algorithm flow was given. The exper-imental simulation results showed that the agricultural remote sensing image enhancement based on firefly algorithmdata was better than histogram algorithm, Retinex algorithm, wavelet transform, fuzzy clustering algorithm about en-hancement measure evaluati

8、on, phase congruency and universal image quality index, so that it can be used for agri-cultural remote sensing image enhancement.Key words: firefly algorithm; Beta function; agricultural remote sensing image; enhancement; gray通过遥感图像可对农田、森林、草场等进行监测,也可以对土壤侵蚀、沙尘暴等进行评估。但是遥感图像在获取的过程中由于受到大气的散射、反射、折射或者天气等

9、的影响,获得的图像难免会带有噪声,如对比度不够、图像模糊;有时总体效果较好,但所需要的信息不够突出,如线状地物或地物的边缘部分。通过图像增强技术,改善遥感图像质量,提高图像目视效果,使遥感应用者易于从经过增强处理的遥感影像上获得所感兴趣的有用信息,快速实现从遥感数据向有用信息的转化1 -2。当前,遥感图像增强有直方图、小波变换、Retinex理论、模糊聚类等算法。国内的研究有:直方图算法能够从全局显著提高图像对比度,在图像增强中取得了不错的应用效果,但仅能描述图像中的奇异点的特性和位置,无法准确描述图像的“沿”边缘特征,会使图像在低亮度区域和高亮度区域对比度过饱和3;小波变换能有效抑制图像的噪

10、声4,但是在实际应用中仍然存在不能保持图像的细节及边缘信息,图像清晰度低等缺陷,只能捕获图像有限的方向信息。国外的研究有:美国物理学家Land等提出的基于人类视觉的光亮度和色彩感知模型的Retinex理论5,使图像具有锐化、动态范围压缩大或颜色恒常性等优点,但容易放大噪声、图像明暗对比度强烈易产生光晕现象以及图像过增强等现象;模糊聚类算法设计简单6,易于实现,但是不能对无中心的像素数据进行聚类,如果针对这种情况改进,时间复杂度和空间复杂度都增加很多。为了解决这个问题,国内外学者研究了仿生学算法,提高了寻优效率,但处理速度仍较慢。农业遥感图像增强算法不仅需要提高图像的对比度,有效地抑制图像噪声,

11、而且应当使低亮度区域和高亮度区域的像素失真减到最小。萤火虫算法是近年来提出的群体智能算法,其在高光谱影像分类、遥感波段选择等领域已经有部分研究,目前国内外关于其能否用于农业遥感图像增强的研究未见报道。本研究拟采用萤火虫优化算法利用非完全Beta函数通过实验仿真研究萤火虫算法能否用于遥感图像增强,为萤火虫算法在农业遥感领域的应用提供理论与算法支持。1 农业遥感图像增强模型1.1 非完全Beta函数采用非完全Beta函数来实现遥感图像的增强7,归一化的非完全Beta函数F(u)定义为:F(u) = B -1 ,( ) u0t -1(1 - t) -1dt。 (1)B ,( ) = 10t -1(1

12、 - t) -1dt。 (2)其中:B ,( )为Beta函数,0 10,0 10,和的取值控制着F(u)的形状,如图1所示。设yij表示像素(i,j)的原始灰度值,yij为处理后的灰度值,非完全Beta函数增强算法的步骤如下。1)对原始图像每个像素进行归一化变换:yij = yij - Lmin( )/ Lmax - Lmin( )。 (3)其中:yij为像素(i,j)的归一化灰度值,Lmax和Lmin分别为原始图像的最大和最小灰度值。2)对归一化图像每个像素进行增强处理:yij = F yij( )。 (4)其中:函数F为非完全Beta函数。3)根据图像灰度值范围,对每个像素进行反变换得到

13、结果图像:yij = Lmax - Lmin( )yij + Lmin。 (5)其中:Lmax和Lmin分别为图像灰度值范围中的最大和最小灰度值,与结果图像的位数有关,对于8位结果图像,Lmax =255,Lmin =0。只有和获得最优取值时,遥感图像的增强效果才能最佳,因此可将参数,选取问题转换为优化问题进行处理。图1 Beta函数图形Fig.1 Beta function graphics7701郭红山,等.基于萤火虫算法的农业遥感图像增强研究1.2 细节增强为了适应人眼的视觉效果8,在低灰度和高灰度级区域,将灰度间隔拉伸;在中等灰度区域,对灰度级间隔过大的部分缩短其间隔,使人眼更好分辨。

14、人眼最小灰度分辨力函数I(k)随灰度k的函数关系,用分段函数表示为:I(k) =- k8 +6,0 k 32- k32 +3,32 k 64k128 +12 ,64 k 129k64 -1,129 k 255。 (6)设输入图像X的灰度值分布为p(x),其中x X,将之视为一个用高斯混合模型模拟的由N个函数线性组合而成的统计密度函数:p(x) = Nn =1p(wn)p(x wn)。 (7)其中:p(x wn)为第n个成分的密度,p(wn)为高斯成分wn的先验概率。则高斯模型中的成分密度函数:p(x wn) = 122wnexp - x - ( )222wn。 (8)其中:,2wn分别为第n个

15、成分的均值和方差。2 萤火虫算法的农业遥感图像增强模型寻优2.1 萤火虫算法描述萤火虫算法包括相对荧光亮度和吸引度2个要素9 -10,相对荧光亮度计算为:l = l0 e - rij。 (9)其中:rij为萤火虫i与j之间的空间距离,l0为萤火虫在rij = 0处的最大萤光亮度,l0与适应度函数有关,适应度函数越优自身亮度越高,为光强吸收系数,亮度体现了萤火虫所处位置的优劣并决定其移动方向。吸引度计算为:q = q0 e - r2ij。 (10)其中:q0为最大吸引度,吸引度决定了萤火虫移动的距离。萤火虫i被吸引向萤火虫j移动位置为:x i = xi + q (xj - xi) + h (ra

16、nd -0.5)。(11)其中,xi、xj为萤火虫i和j所处的空间位置,为x i更新后位置,h 0,1 为步长因子,rand为0,1 上服从均匀分布的随机因子。2.2 基于伪差分算法的迭代策略通过伪差分操作进行迭代,这样萤火虫通过多次移动后,所有个体都将聚集在亮度最高的萤火虫的空间位置上,从而实现寻优。荧光亮度值更新:li(t) = li(t -1) + J(xi(t)。 (12)其中:J(xi(t)为每只萤火虫i在t次迭代位置xi(t)对应的适应度函数值;li(t)代表t次迭代萤火虫i的荧光亮度值;为适应度提取比例,荧光挥发因子 0,1 。每个萤火虫在其动态决策域半径内选择荧光亮度值比自己高

17、的个体组成其邻域集合11 -14,设t次迭代选择移向萤火虫i的领域集合Ni(t)内萤火虫个体j的概率pij(t):pij(t) = lj(t) - li(t)k Ni(t)(lk(t) - li(t)。 (13)其中:领域集Ni(t) = j:dij(t) rid(t);li(t) lj(t),0 rid(t) rs,rs为萤火虫个体的感知最大半径,dij(t)为t次迭代萤火虫个体i和j之间的距离,rid (t)为萤火虫i在t次迭代的决策范围。位置更新:xi(t +1) = 1s xi(t) + s xj(t) - xi(t)xj(t) - xi(t)。(14)其中:s为移动步长。动态决策域半

18、径更新:rid (t + 1) = minrs,max0,rid(t) + ni(t) - Ni(t)( )。 (15)其中:动态决策域更新率 0,1 ,ni(t)为t次迭代萤火虫个体i邻域集内包含的萤火虫数目阈值。伪差分操作更新为:x i(t) =xi(t) +(1 -)(xbest(t) - xi(t) +(xm(t) - xn(t)。 (16)8701浙江农业学报 第28卷 第6期其中:为0,1区间的随机数,xbest(t)为当前种群中最优的个体,xm(t)和xn(t)分别为种群中任意的2个个体,为变异因子,这样减少算法寻优过程的盲目性,加快算法寻优速度。2.3 非完全Beta函数最佳参

19、数选择利用萤火虫寻优算法确定非完全Beta函数的最佳参数,首先找出萤火虫中局部最好位置,( )b,按动态伪差分扰动对每个萤火虫最差解,( )w进行更新,直至算法满足对农业遥感图像增强的收敛要求,最终确定,( )的最优值,( )g,算法收敛是在给定的精度 =10 -5下,如果满足,( )g - ,( ) ,则认为算法收敛,否则不收敛, ,( )g为求得的最优解, ,( )为理论最优解。算法流程:(1)输入遥感图像,初始化萤火虫算法参数;(2)计算每个萤火虫所对应的遥感图像像元的灰度值,将适应度函数值作为其灰度值,并排序,得到亮度最大的萤火虫位置;(3)按照伪差分算法进行迭代,在决策域半径范围内更

20、新;(4)获得,参数,满足算法收敛要求,进行步骤(5),否则转(3);(5)输出遥感图像。3 结果与分析为了检验农业遥感图像增强的效果,选择2幅大小为512 512,灰度级为256的遥感图像作为测试样本。在实验中,依次进行直方图、Ret-inex算法、小波变换、模糊聚类算法、萤火虫算法对比实验,实验PC配置为CPU 3. 6 GHz、内存4GB、IntelH61主板,集成显卡,Matlab 7. 0实现仿真。3.1 视觉分析本研究选取2幅不同的遥感图像,图像由于光照条件不好或有薄雾,图像模糊不清,对比度差,需要进行图像增强处理,实验结果如图2、图3所示。在图2和图3的结果可以看出:萤火虫算法能

21、够增强遥感图像中的边缘信息片段,且边缘定位准确,可以看到细节和灰度拉伸至人眼敏感的中等区域上,适合人眼观察。直方图算法、Ret-inex算法、小波变换算法、模糊聚类算法在增强图像的同时,将一些杂乱背景也增强了,且导致了细节的下降,细节部分如图2中的居民点和农田被归为同一灰阶,无法细分,同时有过增强的现象。3.2 评价指标a,图像1;b,直方图增强结果;c, Retinex算法增强结果;d,小波变换增强结果;e,模糊聚类算法增强结果;f,萤火虫算法增强结果图2 图像1不同算法增强结果Fig.2 Different algorithms enhancement results of image N

22、O.19701郭红山,等.基于萤火虫算法的农业遥感图像增强研究a,图像2;b,直方图增强结果;c, Retinex算法增强结果;d,小波变换增强结果;e,模糊聚类算法增强结果;f.萤火虫算法增强结果图3 图像2不同算法增强结果Fig.3 Different algorithms enhancement results of image NO.23.2. 1 EME指标图像细节增强评价指标(enhancement meas-ure evaluation, EME),其原理是将图像分成k1 k2个区域,通过计算区域中灰度最大值和最小值之比的对数均值作为评价结果。其物理意义是表现出图像局域灰度的变化

23、程度,当局域灰度变化越强烈,细节表现越丰富,值越大。其表达式:EMEk1,k2 = 1k1k2k2l =1k1k =120log Imax(k,l)Imin(k,l) 。 (17)其中:Imax(k,l),Imin(k,l)分别为图像区域(k,l)中的灰度最大值和最小值。3.2. 2 PC指标相位一致性(phase congruency,PC)源于局部能量模型,其实质是图像各个频率成分的相位相似度的一种度量方式,即使在对比度非常弱的图像局部区域,相位一致性也能够很好地检测图像轮廓细节特征。相位一致性是一种频率域的测度值,是一个无量纲参数,取值范围0,1,越大说明某处的特征越明显,表达式为:PC

24、(x) = E x( )nAn x( )。 (18)其中:E x( )为局部能量,An x( )为傅里叶分量在位置x处的振幅。3.2. 3 DI指标偏差指数(difference index,DI)用来比较增强图像和标准图像偏离的程度,定义重建图像的偏差指数的表达式为:DI = 1m n mi =1 nj =1|yij - xij |xij 。 (19)式中:m,n为图像宽和高;xij,yij分别为原始图像和增强图像上对应(i,j)点的灰度值。偏差指数的大小反映了增强图像对标准图像的保持程度,偏差指数越大,增强效果越差。增强评价指标如表1所示,其中每个测试指表1 增强评价指标Table 1 E

25、nhanced evaluation index算法图像1EME PC DI图像2EME PC DI直方图0. 7929 0. 8102 0. 7919 0. 8152 0. 8059 0. 7082Retinex算法0. 7987 0. 8301 0. 9218 0. 8519 0. 8305 0. 8901小波变换0. 8138 0. 8408 0. 8913 0. 8418 0. 8269 0. 9001模糊聚类0. 8621 0. 8860 0. 6068 0. 8362 0. 8363 0. 7621萤火虫算法0. 9991 0. 9933 0. 3565 0. 9793 0. 99

26、56 0. 34470801浙江农业学报 第28卷 第6期标进行20次的蒙特卡洛测试,再取均值。从增强评价指标可以看出,萤火虫算法的评价指标结果较优,因此,其他的农业遥感图像增强方法无法兼顾细节信息,而萤火虫算法却能较好地满足细节增强的需求。4 小结本研究建立农业遥感图像增强模型,通过萤火虫优化算法在其动态决策域半径进行数据更新,对模型求解,实验仿真显示萤火虫算法与直方图算法、Retinex算法、小波变换算法、模糊聚类算法相比,图像细节增强评价指标值较大,相位一致性指标值接近于1,通用质量评价指标值接近于1,数据指标较优,萤火虫算法具有能够增强遥感图像的细节信息的性能。因此,为遥感图像目标的检

27、测提供了一种新思路。参考文献:1 鲍程辉,贺新光,蒋卫国.非下采样Contourlet变换和脉冲耦合神经网络相结合的遥感图像融合方法J.遥感信息,2015,30(2):52 -56.2 吴乐芹,姜春,陈进栋,等.面向城市道路目视识别的遥感图像滤波增强方法J.地理空间信息,2015,13(2):104 -106,109.3 刘陶胜,李沛鸿,李辰风.结合像素频率分布特征的遥感图像自适应线性增强J.江西理工大学学报,2014,35(5):40 -44.4 伊力哈木亚尔买买提,谢丽蓉,孔军.基于PCA变换与小波变换的遥感图像融合方法J.红外与激光工程,2014,43(7):2335 -2340.5 李

28、建存,周丽梅,刘军.基于多尺度Retinex理论的遥感图像增强算法J.西安工业大学学报,2014,34(1):27 -33.6 杨波,贾振红,覃锡忠,等.基于改进模糊算法和相对熵的遥感图像增强J.激光杂志,2014,35(2):37 -38.7 姜建国,周佳薇,周润生,等.一种采用改进细菌觅食优化算法的图像增强方法J.控制与决策,2015,30(3):461-466.8 王才子,孙士保,段建辉,等.基于人眼视觉的遥感图像广义直方图增强算法J.合肥工业大学学报,2012,35(9):1208 -1211,1292.9 刘洲洲,王福豹.基于离散萤火虫压缩感知重构的无线传感器网络多目标定位J.光学精

29、密工程,2014,22(7):1904 -1911.10 BAYKASOLU A, OZSOYDAN F B. Adaptive firefly algo-rithm with chaos for mechanical design optimization problemsJ. Applied Soft Computing,2015,36(11):152 -164.11 吴鹏.萤火虫算法优化最大熵的图像分割方法J.计算机工程与应用,2014,50(12):115 -119.12 AYDIN I. A new approach based on firefly algorithm for vi

30、-sion-based railway overhead inspection systemJ. Measure-ment,2015,74(10):43 -55.13 KOUGIANOS E, MOHANTYA S P. A nature-inspired fireflyalgorithm based approach for nanoscale leakage optimal RTLstructureJ. Integration,the VLSI Journal,2015,51(9):46-60.14 YU S H, ZHU S L, MA Y, et al. A variable step size fireflyalgorithm for numerical optimization J. Applied Mathemat-ics and Computation,2015,263(15):214 -220.(责任编辑 张 韵)1801郭红山,等.基于萤火虫算法的农业遥感图像增强研究