基于暗原色融合和维纳滤波的单幅图像去雾-杨爱萍.pdf

《基于暗原色融合和维纳滤波的单幅图像去雾-杨爱萍.pdf》由会员分享，可在线阅读，更多相关《基于暗原色融合和维纳滤波的单幅图像去雾-杨爱萍.pdf（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第49卷 第6期 2016年6月 天津大学学报(自然科学与工程技术版) Journal of Tianjin University(Science and Technology) Vol.49 No.6Jun. 2016 收稿日期：2014-11-06；修回日期：2015-01-23. 基金项目：国家自然科学基金资助项目(61372145，61472274，61201371). 作者简介：杨爱萍（1977 ），女，博士，副教授. 通讯作者：杨爱萍，. 网络出版时间：2015-03-12. 网络出版地址：http:/ DOI:10.11784/tdxbz201411021 基于暗原色融合和维纳滤

2、波的单幅图像去雾 杨爱萍，刘华平，何宇清，白煌煌，宋曹春洋 (天津大学电子信息工程学院，天津 300072) 摘 要：针对在雾霾等恶劣天气下捕获的户外场景图像对比度降低、颜色失真等问题，对基于暗原色先验的去雾方法进行改进，应用小波变换将块暗原色和点暗原色进行融合后，得到新的透射率估计，并利用自适应维纳滤波细化透射率同时提出了四分加权法重新估计大气光，使得大气光更具鲁棒性实验结果表明，本文方法不仅能有效恢复清晰的无雾图像，而且能够大幅提升运行速度，便于实时应用 关键词：图像去雾；暗原色融合；小波变换；透射率估计；四分加权法 中图分类号：TP391 文献标志码：A 文章编号：0493-2137(2

3、016)06-0574-07 Single Image Dehazing Based on Dark Channel Fusion and Wiener Filtering Yang Aiping，Liu Huaping，He Yuqing，Bai Huanghuang，Songcao Chunyang (School of Electronic Information Engineering，Tianjin University，Tianjin 300072，China) Abstract：Haze is one of the major factors that cause color d

4、istortion and contrast loss of the outdoor image. To re-duce these effects，an improved dehazing method based on dark channel prior is proposed in this paper. A new trans-mission is estimated after the fusion of patch-based dark channel prior and pixel-based dark channel prior using wave-let transfor

5、m，and adaptive wiener filter is introduced to further refine the transmission estimation. Simultaneously，in order to make the result more robust，we propose a method called weighted quadtree subdivision to estimate atmos-pheric light. Comparative experimental results demonstrate that the proposed met

6、hod is effective not only in restoring the clear images，but in speeding up the computation，which is appropriate for real time application. Keywords：image dehazing；dark channel fusion；wavelet transform；transmission estimation；weighted quadtree subdivision 雾天降质图像清晰化一直是计算机视觉领域研究的重要内容，简单有效的图像去雾对提高视觉系统的可

7、靠性和鲁棒性具有非常重要的意义当前的图像去雾方法主要分为两类1-2：一是基于图像增强的方法，该方法不考虑图像退化的真正原因，去雾后的图像容易出现信息丢失和色彩失真；二是基于物理模型的方法，即基于大气散射规律建立图像退化模型，并能充分利用图像的先验知识，具有内在的优越性，是该领域的研究热点 为了适应各种应用场合，基于单幅图像的去雾算法成为当前的研究热点Tan3观察到无雾图像相比有雾图像具有g3101高的对比g2318，g6447g6398g3107大化图像g2124g6537对比g2318g6393到去雾的g4395的g877该方法g2450g4244了g815物理模型g700g3514g132

8、0g6440射g4106，g2409到的无雾图像g7079色g6398于g7119和，g3071g3718g4681g1322g1834g790g4214g3061g6338效应Fattal4g987g6009g6440射g4106g704g5866g7014g2693g2389g2124g6537不相g1057，g6611用g4055立成分分g3169g868计g6440射g4106，g877是该方法对于g7079色信息要g3514g6358高，g1834理g3651雾图像g3024g855g790g4214g6358大失真Tarel g47495利用g4017g3443的g726g984g

9、3801g3592对大气g5105散g1125g2972g6408g5856g868计，并g4934合色g6075g3041射g2051现图像去雾g877该方法g1311g2972g6358g1842，不易g6075g2974，g1358g3024在g3071g3718g4681g1322的g1244域不能g5516g2409g5345g1885的去雾效g3174为了g5964g1092g700g6423g6855g7076，He g47496提出了一种基于g3084原色先验的去雾方法，g1g1g1 2016年6月 杨爱萍等：基于暗原色融合和维纳滤波的单幅图像去雾 575 对一g5321g26

10、49g1840图像g4212g5295g3651雾场g3071g1704具有g5345g1885的去雾效g3174g15 g877是g1060g2124g6924性在于算法g1836g3133g2318高，g1834理g6451g2318g2596g15 g6634对g6407g782g6855g7076，g2021g5095g834相g4952提出了一系g1140g2951g6408算 法7-9g15 g3124g2976在g3084原色先验去雾方法的基g4520g700，g6634对点g3084原色和g1712g3084原色g3514g1320g6440射g4106的g5028g6941，

11、提出g2097二g5095g6408g5856g5815合重g2995g868计g6440射g4106，g3944后g6611用g5292适应g4964g4905g3801g3592对g6440射g4106g6408g5856g4921化，g1358g3024提出g1661分g1194g3134法对大气g1039g6408g5856g868计g2051验g4934g3174g5866g3033，g3124g2976方法不g808g5516g2409了g5345g1885的去雾效g3174，g5098g708能g1844大幅g2318提高算法的g6401g5856g6451g2318 1 暗原色

12、先验去雾算法 1.1 雾天图像退化模型 大气散射模型g2848g6423了雾天图像的退化g6398g46536,10，即 g1g1g1 () ()() 1 ()Ix Jxtx tx A=+g1(1)g2351g726：()Ix为观察到的雾图像强g2318； ()J x 为不受空气介质干扰的g2051际g3071物g1039线强g2318；A为大气g1039；()tx为场g3071g6440射g4106； ()()J xtx为直接衰减项； 1 ( )Atx 为大气g1039成分图像去雾的g4395的是g815已知图像 ()Ixg726恢g1836 ()J x 1.2 基于暗原色先验的去雾方法 g

13、3084原色先验是g6447g6398统计观察大量g2649g1840无雾图像提出的：在大g1842g2972无雾图像的非天空g1244域，某g782像素g5295少有1个g7079色g6447道具有非常小的强g2318g984，g1060g3107小强g2318g984趋于 0，即对于无雾图像 ()J x ，g1060g3084原色先验dark()J x 可g5866示为 g1g1g1dark, , ()( ) min ( min ( ( ) 0ccrgbyxJx Jy=g1(2)g2351g726： ()cJ y 为 ()J x 的某一个g7079色g6447道的亮g2318g984；()

14、 x 为以x为g726心的方形g1244域 首先g987g6009大气g1039A给定，g708g987g6009g1358g1244域g1712内g6440射g4106t恒定，在g1244域 () x 内对g2351(1)两边g6408g5856g3107小g984 g3801g3592，并应用g3084原色先验理论g2409到 g1g1g1, , ()() 1 () min(min( ()ccrgbyxvx A tx I y= = g1(3)由于A总是正g984，g2409到g6440射g4106g868计为 g1g1g1, , ()() ()() 1 min(min( ) 1ccrgby

15、xIy vxtxAA= = g1 (4)g2351g726 ()vx为大气g5105散g1125g2972g1060g726，大气g1039A的g868计方法是：g5516g1320g3084原色g7260.1%,亮g2318g3107大的像素，g6407g782像素所对应的原图 ()Ixg726强g2318g3107大的像素点作为大 气g1039 直接应用g2351(4)g2409到的g6440射g4106恢g1836图像g855g790g4214g3061g6338伪g2389，g2976献6g6611用软抠图方法对g6440射g4106g6408g5856优化，g3107后g2097优化

16、后的g6440射g4106代入g2351(1)，g2409到清晰的无雾图像 2 本文改进算法 2.1 基于暗原色融合和自适应维纳滤波估计透射率 2.1.1 基于暗原色融合的透射率估计 基于g3084原色先验的去雾方法是在g2124g6537图像g1712内g868计g6440射g4106，g6407样g2409到的g6440射g4106在g1712内是恒定的g877在g2051际图像g1834理g726，g1712内的g6440射g4106并不总是恒定不g1322的，尤g1060是在g3718g2318g790g4214大跳跃的边缘，g855导致g6440射g4106图出现严重的g1712效应

17、(见图 1(b)，使恢g1836图像g790g4214g3061g6338伪g2389(见图 1(c)g2976献6g6611用软抠图技术对g6440射g4106g6408g5856优化，虽g3944g2409到了g5345g1885的去雾效g3174，g877计算量巨大，g3024间消g5105严重g2976献5基于g726g984g3801g3592g3514g1320大气g5105散g1125g2972g2051现去雾，虽能减少算法g6401g5856g3024间，g877不能彻底去除g3061g6338伪g2389 g3084原色先验理论以统计g2021为基g4520，当图像g1712

18、g6358大g3024，g1712g726包含g3084像素的概g4106也g855g6358高，g3084原色先验理论成立g877当g1712g6398大g3024，g1712g726g6440射g4106恒定的g987g6009失效，导致色彩畸g1322因此，g1712的大小是g3084原色先验去雾算法g726一个很g1057键的因素11 g1（a）有雾图像 （b）块暗原色估计透射率 （c）图（b）对应清晰图像 （d）点暗原色估计透射率 （e）图（d）对应清晰图像 图1 块效应和晕轮伪影 Fig.1 Block artifacts and halos phenomenons g1g1g1

19、576 天津大学学报(自然科学与工程技术版) 第49卷 第6期 He g47496提出的方法是在图像g2124g6537g1712内g3514g1320g3084原色(即g1712g3084原色)，该方法出现g1712效应的原因是g6407g782g1712的边缘包含一g782错误的高频信息换言之，g6611用g1712g3084原色g2409到的g6440射g4106g726，低频g6537分接近g2051际的g6440射g4106，g877g1358g3024也带来许g1842g1712边缘g4749高频信息(见图 1(b)，g815g5098导致了g1836原图像g726的g3061g6

20、338伪g2389(见图1(c)g1g1g1Tarel g47495提出的方法相当于在单像素点g1834g3514g1320g3084原色(即点g3084原色)，场g3071g726的g4921节g2097g855g2409到保留，g708不g790g4214g1712效应，也就是说使用点g3084原色g2409到的g6440射g4106g726，高频信息接近真g2051的g6440射g4106g877是g6407种情况下，低频g6537分的g868计不g1844精确(见图 1(d)，恢g1836图像g855g790g4214严重的色彩失真(见图1(e) 基于以g700分g3169，g1712

21、g3084原色的低频和点g3084原色的高频分别接近于真g2051g3084原色的低频和高频，g3124g2976分别g3514g1320雾天图像的g1712g3084原色和点g3084原色，g3944后对两种g3084原色图g6408g5856小g3592g1322换，并g2097g1712g3084原色的低频系g2972和点g3084原色的高频系g2972g6408g5856g5815合，作为要恢g1836图像g3084原色的小g3592g1322换，g3107后g6408g5856小g3592反g1322换g2409到空间域g3084原色图，g6408g5098g2409到g6440射g

22、4106g868计具体算法流g4653见图2 g1图2 暗原色融合估计透射率框图 Fig.2 Block diagram of the estimated transmission mapusing the dark channel fusion 2.1.2 基于局部自适应维纳滤波的透射率细化 g3084原色g5815合g2409到的g6440射g4106虽g3944g6358 He g47496提出的方法g3101具鲁棒性，g877在g5815合g6398g4653g726不可避免g855引入一g782错误的g4921节信息因此，g3124g2976g6611用g2124g6537g5292

23、适应g4964g4905g3801g359212对g6440射g4106g6408g5856g4921化 g987g6009大气g1039A已知，由于错误的g4921节信息是在g3084原色g5815合的g6398g4653g726引入的，分g3169g2351(3)可知，g5815合后g2409到的g3084原色可看做是大气g5105散g1125g2972和错误g4921节信息之和 g1g1g1 () () ()dx vx nx=+g1(5)g2351g726： ()dx为g5815合后的g3084原色；()nx为错误的g4921节信息g987g6009 ()vx和 ()nx是相互g4055

24、立的 给定g2351(5)，g6611用g2124g6537g5292适应g4964g4905g3801g3592器来g868计g6611样窗口 ()x 内的 ()vx，有 g1g1g1222()() () () ()()vnvdvxvx x dx xx+=+ g1(6)g2351g726：()vx为大气g5105散g1125g2972g868计； ()vx 、2()vx 分别为 ()vx在g6611样窗口内的g1704g984和方差； ()dx 为 ()dx在g6611样窗口内的g1704g984；2n 为g4921节信息的方差，g987g6009g1060在g2974幅图像g726是恒定的

25、 此g1840，还需要对错误的g4921节信息 ()nx的方差2ng6408g5856g868计g987g6009d和n不相g1057g708ng1704g984为 0，根据g2351(4)可知 g1g1g1222() ()dvnxx=+g1 (7)g2351g7262()dx 为 ()dx在g6611样窗口内的方差g704 Tan 和Fattal 的方法类似，大气g5105散g1125g2972在大的g1320样窗口内是互相g1057的，并g708方差2()vx 很低g987g600922()vnx ，可用g3084原色方差的全g2124平g1704作为g4921节方差的g868计 g1g1

26、g112201 ()Mndx xM=g1(8)g2351g726M为g2974幅图像的像素点g2972 g868计g2409到 ()vx的g1704g984和方差以及g4921节信息的方差后，g6447g6398g2351(6)可g2409到大气g5105散g1125g2972 ()vxg3107优g868计，g3107后计算g6440射g4106为 g1g1g1()() 1vxtxA= g1 (9)g2351g726为常g2972的作用是保留g6537分雾，g6407是因为如g3174彻底去除雾的存在，g2974体还原图像g2097g855不真g2051并g708丢失g3718g2318感g

27、4921化后的g6440射g4106可有效去除g1712效应和g3061g6338伪g2389 图 3 为g6611用g3124g2976方法g2409到的g6440射g4106及去雾效g3174图，g6447g6398对比 He 的算法可以看出，g3124g2976方法不g808可有效去除g1712效应和g3061g6338伪g2389，还能g1844保留图像的真g2051色彩，去雾效g3174g3033显优于He和Tarel 提出的方法 2.2 四分加权法估计大气光 大气g1039的g868计在图像恢g1836g726起到了举足轻重的作用13在常用的去雾算法g726，g6447常g1320

28、图像g726g3107亮点的像素g984作为大气g1039，g877是g2051际g700g3107亮点很g1842g3024候并不是雾g3651g2318g3107大的地方，不能代g5866大气g1039的g984，g6407g782点往往都是场g3071g726的白色物体(如图 4(a)g726矩形g3251g6537分所示)，g1060像素g984高于真g2051的大气g1039，使g2409恢g1836图像g991g3084(如图 4(b)所示)为使大气g1039的g868计g3101具鲁棒性，g3124g2976提出了一种g2995的大气g1039g868计方法g250g250g16

29、61g1g1g1 2016年6月 杨爱萍等：基于暗原色融合和维纳滤波的单幅图像去雾 577 分g1194g3134法 g1（a）He算法估计透射率 （b）图（a）对应的清晰图像 g1（c）本文算法估计透射率 （d）图（c）对应的清晰图像 图3 透射率以及去雾效果图 Fig.3 Transmission and dehazing images g1（a）有雾图像 （b）去雾图像 图4 最亮像素去雾结果 Fig.4 Result of dehazing using the brightest pixel 2.2.1 四分法确定浓雾区域 大气g1039A要在雾g3651g2318g3107高的g12

30、44域确定，g708雾g3651g2318越高的g1244域，像素g984越高，像素间的差g2345越小，g1704g984g704g3221g1101差的差g984就越大，定义g3463个g1244域的差g984为 g1g1g1 () () ()Sn Mn Dn=g1(10)g2351g726： ()Sn为某一矩形g1244域的g1704g984和g3221g1101差之间的差g984； ()Mn和 ()D n 分别为该g1244域的g1704g984和g3221g1101差；n为g1060g726某一g1244域， 1, 2, 3, 4n = 首先，g2097图像g6443g2373g23

31、51地分成4个相g1358大小的矩形g1244域，分别计算 4 个g1244域的差g984，并g6436g1320差g984g3107高的g1244域，g3944后重g1836g700g6423g6398g4653，直到矩形g1244域的大小g3799足g7060先g6009定的g6879g984，g3107后g6436g1320的g1244域就是雾g3651g2318g3107高的g1244域具体g6436g2744g6398g4653如图5所示 g1（a）有雾图像 （b）雾浓度最高的区域 图5 四分法选择图像块的过程 Fig.5 Process of selecting the regio

32、n using quadtree sub-division 2.2.2 加权法估计大气光 g6436定g3651雾g1244域 ()Yx后，在此g1244域内g868计大气g1039现有的方法是g2097图像g726g3107亮点的像素g984作为大气g1039，g877该方法具有一定的g6947机性，g868计大气g1039的点可能是图像g726的白色物体为g5964g1092g700g6423g6855g7076，g3124g2976g6611用g1194g3134法g868计 g2097图像g1712 ()Yxg726所有的像素点分为 2 个g1244域：所有g3877g2318g984

33、大于g1704g984的点g2138于亮g1244，g3877g2318g984小于g1704g984的点g2138于g3084g1244，g6009亮g1244和g3084g1244的像素点个g2972分别为brightN 和darkN 定义 ()Yx在33 的窗口内的g3084原色为Edark()Yx，g5516g1320亮g1244和g3084g1244g3084原色的g3107大g984，分别g5997为brightA 和darkA g987g6009brightA 和darkA 分别在点 ()Ym和 ()Yng1834g1320g2409，g1142大气g1039A计算g2351为

34、g1g1g1bd() ()AWYm WYn=+g1 (11)g2351g726：brightbYYNWMN=；darkdYYNWMN=，并g708g3799足bd1WW+=；YYMN 为图像g1712 ()Yxg2118g2089g6407样，当亮g1244点比g6358g1842g3024，A由 ()Ym主导，当g3084g1244点比g6358g1842g3024，Ag1142由 ()Yn主导，g877是无论g884种情况， ()Ym和 ()Yn都是g1056g1358g1318g2776作用，相互g1155g4893，相互g5865g1001，g2409到一个合理的大气g1039g984

35、，g5516g2409视觉g700g6358g5292g3944的重建图像(如图6所示) g1（a）有雾图像 （b）去雾图像 图6 四分加权法去雾结果 Fig.6 Result of dehazing using weighted quadtree subdi-vision g1g1g1 578 天津大学学报(自然科学与工程技术版) 第49卷 第6期 g1g1g4971g700，g1661分g1194g3134法g868计大气g1039A的具体g3429g7200如下： (1) g2097图像分成4个相g1358大小的矩形g1244域，分别计算g3463个g1244域的差g984 ()Sn；

36、(2) g6436g1320差g984g3107大的g1244域，重g1836g3429g7200(1)，直到g3799足g7060先g6009定的g6879g984，确定雾g3651g2318g3107高的g1244域 ()Yx； (3) g2097g1244域 ()Yx内的像素点分为g3084g1244和亮g1244，g5516g13202个g1244域内g3084原色的g3107大g984； (4) 确定g3107大g984在图像g726的g878g5045，利用g1194g3134法g3514g1320大气g1039A 2.3 复原图像 g1g1图7 是g3124g2976提出方法的流

37、g4653g6372入 1 幅有雾图像，首先g868计大气g1039A，利用g1661分法确定图像g726雾g3651g2318g3107高的矩形g1244域，并在该g1244域内利用g1194g3134法g868计g2409到大气g1039g3944后g5516g1320有雾图像的g1712g3084原色和点g3084原色，g6447g6398g3084原色g5815合g2409到g4835g4244g868计的 ()tx为了g2951g1554 g1图7 本文算法流程 Fig.7 Flow chart of the proposed algorithm in this paper()tx的

38、质量，g1330利用g4964g4905g3801g3592器对g1060g6408g5856g4921化，g3107后g2097g868计的A和 ()tx代入g2351(1)g2409到清晰的无雾图像 g1g1g10()()max( ( ), )Ix AJx Atx t=+g1(12)为避免去雾后的图像包含g1641g1829，g6009定1个g6440射g4106的下g6924g9840t ，g3124g2976g13200.1 3 实验结果的比较与分析 3.1 主观评价 为了验g6012去雾方法的有效性，g3124g2976g704 Tarel g47495提出的g726g984g380

39、1g3592方法以及 He g47496提出的g3084原色先验去雾方法g6408g5856比g6358图 8 给出了g6537分g2056g1840有雾图像的去雾效g3174 对比3g4919g2051验g4934g3174可以看出，g3124g2976方法能g1844有效地去除图像g726的雾气，提高了图像的可见g2318，g704 He提出的方法去雾效g3174相近，某g782情况下优于g1060方法g815g4921节g700分g3169，g704He的方法相比，g47391g4919和g47392g4919图像天空g1244域恢g1836效g3174g3101g1885，g4739

40、3 g4919效g3174相似；g704Tarel 方法相比，3 g4919图像g726，Tarel 方法g7079色g6398于g7119和，g3124g2976方法g3033显优于g1060g4934g3174 图9给出g1330g1840g1110幅图像的去雾效g3174，具有类似的g4934论 （a）雾天图像 （b）He 算法 （c）Tarel 算法 （d）本文算法 图8 雾天图像恢复结果对比(图像从上到下分别为City，Hongkong，Ny1) Fig.8 Comparison of the recovered scenes for hazy images (From top to

41、 bottom，the images are City，Hongkong，Ny1) g1g1g1 2016年6月 杨爱萍等：基于暗原色融合和维纳滤波的单幅图像去雾 579 （a）雾天图像 （b）He算法 （c）Tarel 算法 （d）本文算法 图9 雾天图像恢复结果对比(图像从上到下分别为Trees，House，Tanmen) Fig.9 Comparison of the recovered scenes for hazy images (From top to bottom，the images are Trees，House，Tanmen) 3.2 客观评价 下g7014g815g640

42、1g5856g3024间、可见边g3288g2318g4749方g7014g6408g5856比g6358 He 去雾方法g726g6611用 soft matting 对g6440射g4106g6408g5856优化，g1060g2051质是g6408g5856大规模g4649g4255线性方g4653的g3514g5964，计算g1836g3133g2318g6358高；Tarel 的去雾方法g3024间消g5105在 2 g3410g726g984g3801g3592；g3124g2976方法的g3024间主要消g5105在g5516g1320g1712g3084原色先验g6398g46

43、53g726的g3107小g984g3801g3592以及g6440射g4106优化g6398g4653g726的g4964g4905g3801g3592，g3801g3592g6398g4653近似于线性，g2955g3124g2976方法g704 Tarel方法的g3024间消g5105相近 可见边g3288g2318法14是g4395前去雾领域常用的g4398g6015方法，该方法g6611用有雾图像g704恢g1836图像的可见边g2972g4395比( e )和平g1704g3288g2318比( r )来g2054观g6015g840图像的去雾效g3174即 g1g1g100rnn

44、en= g1(13)g1g1g10rgrg= g1(14)g2351g726：rn 和0n 分别为有雾图像 ()Ix和恢g1836图像 ()J x的可见边g2972g4395；rg 和0g 分别为有雾图像 ()Ix和恢g1836图像 ()J x 的平g1704g3288g2318一g5321情况下，e和r 的g984越大，g5866g3033g1834理后图像边缘强g2318越大，g2972量越g1842 g4934g3165相似性(structural similarity，SSIM)算法15是图像g1834理领域常用的质量g6015g840g3221g1101由于g800类视觉很容易g81

45、5图像g726g2713g1320出g4934g3165信息，2 幅图像的SSIMg984可作为g6015g840图像去雾质量g1885g1709的g3221g1101 g1g1g1()()1222 22SSIMIJ IJIJ IJCC +=+ +g1(15)g2351g726：I 、J 、I 、J 分别为图像 ()Ix、 ()J x 的g1704g984和方差；IJ 为 ()Ix、 ()J x 的g1256方差；1C 、2C 为2 个比g6358小的常g2972SSIM g984越趋近于 0，两幅图g4934g3165相似性越低，SSIM g984越趋近于 1，两幅图的g4934g3165相

46、似性g1越高 g5866 1 给出了图 9 所示 3 幅图像在不g1358去雾方法下各g2760g3221g984的大小，g2051验图像大小分别为 300g102216、350g102262和350g102262 表1 各方法相关指标对比 Tab.1 Comparison results in related measuring index of various methods 图像 算法 时间/s e r SSIM Trees He Tarel本文9.45 0.95 0.79 0.569 0 1.020 3 0.704 6 1.022 8 1.920 6 1.987 3 0.761 0 0

47、.799 6 0.812 6 House He Tarel本文12.410 1.86 1.12 0.254 1 0.129 5 0.271 3 1.211 3 1.198 9 1.464 8 0.967 4 0.795 9 0.885 5 TanmenHe Tarel本文12.220 1.86 1.12 0.241 1 0.260 1 0.267 4 1.215 0 1.313 6 1.458 2 0.880 9 0.884 5 0.889 1 g1g1g1 580 天津大学学报(自然科学与工程技术版) 第49卷 第6期 由g5866 1 可以看出，g3124g2976方法g704 Tarel 方法相比，g3024间g4244g4480，g704 He 方法相比，g6401g5856g6451g2318g2443了g2097近 10g965；g1358g3024g704 He 和 Tarel 方法相比，g3124g2976方法的e和rg1704g2409到很g1885保g2758虽g3944g1834理 Trees 图像g3024