亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于噪聲類型及強度估計的狹葉錦雞兒葉切片圖像盲去噪

        2017-07-07 00:43:38王海超王春光宗哲英殷曉飛張文霞王曉蓉張海軍
        農(nóng)業(yè)工程學(xué)報 2017年10期
        關(guān)鍵詞:錦雞兒標準差高斯

        王海超,王春光,宗哲英,殷曉飛,張文霞,王曉蓉,張海軍,李 靖,劉 濤

        (內(nèi)蒙古農(nóng)業(yè)大學(xué)機電工程學(xué)院,呼和浩特 010018)

        基于噪聲類型及強度估計的狹葉錦雞兒葉切片圖像盲去噪

        王海超,王春光※,宗哲英,殷曉飛,張文霞,王曉蓉,張海軍,李 靖,劉 濤

        (內(nèi)蒙古農(nóng)業(yè)大學(xué)機電工程學(xué)院,呼和浩特 010018)

        狹葉錦雞兒葉切片顯微圖像在獲取過程中不可避免的受到噪聲污染,會對后續(xù)處理造成不良影響。針對現(xiàn)有噪聲類型未知,去噪算法存在速度慢、效果不理想等問題,該文提出圖像噪聲類型估計-強度估計-去噪這一處理過程,實現(xiàn)對狹葉錦雞兒葉切片顯微圖像降噪目的。首先采用平滑區(qū)直方圖重構(gòu)和擬合法確定噪聲類型;然后在此基礎(chǔ)上,應(yīng)用基于圖像塊的SVD(singular valuable decomposition,SVD)域圖像噪聲強度估計法對噪聲標準差進行估計;最后在確定噪聲類型和強度基礎(chǔ)上,采用幾何均值濾波(geometric mean filtering,GMF)和三維塊匹配濾波(block-matching and 3-D filtering,BM3D)對圖像進行聯(lián)合去噪。試驗結(jié)果表明:該文噪聲類型估計法估計出切片圖像噪聲類型為加性高斯噪聲,高斯函數(shù)對隨機選取的 15幅狹葉錦雞兒葉切片圖像平滑區(qū)域直方圖數(shù)據(jù)點擬合優(yōu)度2R均值為 0.996,平均均方根誤差RMSE(root mean squared error,RMSE)為0.144 6;采用該文噪聲強度估計法估計出的切片圖像噪聲標準差處理標準差較小噪聲,該文算法處理精度、運行速度和穩(wěn)定性等方面存在明顯優(yōu)勢;GMF-BM3D算法在較好去除圖像噪聲同時,極大的保留了圖像紋理、邊緣和細節(jié)等信息,同時極大的提高了算法運行速度,處理后的圖像 BRISQUE(blind/referenceless image spatial quality evaluator,BRISQUE)值為10左右,相當于原圖BRISQUE值的1/2左右。與傳統(tǒng)BM3D算法相比,去噪效果相當,但耗時約相當于傳統(tǒng)BM3D算法的1/9。與小波去噪算法(wavelet threshold,WT)算法相比,雖速度相對較慢,但去噪后圖像BRISQUE值比使用WT法低4左右。因此,該算法較好實現(xiàn)了對狹葉錦雞兒葉切片圖像準確降噪,為其后續(xù)處理提供了可靠技術(shù)支持。

        圖像處理;高斯噪聲;評估;切片圖像;狹葉錦雞兒;噪聲強度;圖像去噪

        王海超,王春光,宗哲英,殷曉飛,張文霞,王曉蓉,張海軍,李 靖,劉 濤. 基于噪聲類型及強度估計的狹葉錦雞兒葉切片圖像盲去噪[J]. 農(nóng)業(yè)工程學(xué)報,2017,33(10):229-238. doi:10.11975/j.issn.1002-6819.2017.10.000 http://www.tcsae.org

        Wang Haichao, Wang Chunguang, Zong Zheying, Yin Xiaofei, Zhang Wenxia, Wang Xiaorong, Zhang Haijun, Li Jing, Liu Tao.Blind image denoising of microscopic slices image ofCaragana stenophylla Pojarkbased on noise type and intensity estimation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(10): 229-238. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2017.10.030 http://www.tcsae.org

        0 引 言

        狹葉錦雞兒(Caragana stenophylla Pojark)是荒漠草原植被群落的優(yōu)勢物種,具有良好的防風(fēng)固沙和保持水土能力,且飼用價值較高,在內(nèi)蒙古荒漠草原帶分布廣泛,其對當?shù)夭菰哪卫砭哂兄匾饔谩*M葉錦雞兒葉切片試驗的主要目的是認識葉片內(nèi)部結(jié)構(gòu),揭示其內(nèi)部結(jié)構(gòu)特征指標對干旱脅迫下荒漠草原退化梯度的適應(yīng)機制,為狹葉錦雞兒保護利用和荒漠草原退化防治提供理論依據(jù)[1]。但狹葉錦雞兒葉切片圖像在獲取、傳輸和存儲等過程中會受到各種因素干擾。其中噪聲是影響圖像最主要和最常見因素,其不但會降低圖像視覺效果,而且會給圖像特征提取、拼接修復(fù)等后續(xù)處理造成不良影響。

        圖像去噪一直是圖像處理領(lǐng)域研究的重點,近年來,國內(nèi)外學(xué)者提出了各種高效的圖像去噪算法,并取得了良好去噪效果。具有代表性的基于空間域去噪算法有:He等[2]提出的引導(dǎo)濾波,此法在已知質(zhì)量較好引導(dǎo)圖像情況下去噪效果較好;Buades等[3-4]提出的非局部均值濾波(non local means,NLM),此法對圖像進行全局搜索,利用圖像中所有像素點進行加權(quán)平均,獲得較好去噪效果。具有代表性的基于變換域去噪算法有:Silva等[5]提出的基于自適應(yīng)小波變換圖像邊緣保持去噪算法,去噪結(jié)果較好兼顧了噪聲消除和邊界保持;Dabov等[6-7]提出的三維塊匹配濾波(block-matching and 3-D filtering,BM3D)結(jié)合了非局部和變換域思想,同時還用到了線性變換閾值、維納濾波和稀疏表示等算法,從而達到良好去噪效果。

        上述各種去噪算法中,執(zhí)行效率較高去噪法對圖像邊緣、紋理等細節(jié)特征不能全面兼顧,效果較差。而效果較好的NLM和BM3D算法計算量相當龐大,執(zhí)行效率較低,限制其實際應(yīng)用。BM3D去噪法需預(yù)先知道噪聲水平[8],以噪聲方差作為先驗輸入?yún)?shù)進行去噪。因此,通常需要利用噪聲水平估計算法指導(dǎo)BM3D算法參數(shù)輸入,以實現(xiàn)自然圖像去噪[9-11]。

        然而,對于自然圖像,事先準確得知噪聲方差十分困難,目前大致噪聲強度估計方法可以分為:基于平滑濾波估計法、基于分塊方差估計法和基于變換域估計法3類?;谄交瑸V波估計法是將濾波前后圖像差值作為圖像噪聲數(shù)據(jù)源,進而達到估計噪聲的目的。Tai等[12]提出的基于 Sobel-Laplacian確定和計算圖像結(jié)構(gòu)邊緣響應(yīng)并估計噪聲標準差,此法優(yōu)于傳統(tǒng)濾波算法。然而,此法濾波參數(shù)的選擇對結(jié)果影響較大,濾波太強會使差值圖像中包含圖像結(jié)構(gòu)信息,導(dǎo)致噪聲估計偏大。濾波太弱又無法完全去除噪聲,導(dǎo)致噪聲估計不足;基于分塊方差估計法是將圖像分割成大量圖像塊,用相對平滑圖像塊的方差估計噪聲強度。Salmeri等[13]采用分塊統(tǒng)計圖像特征,模糊分類法劃分區(qū)域并進行方差加權(quán),最終確定噪聲方差。Tian等

        [14]采用蟻群優(yōu)化算法確定圖像中平滑子塊,以各子塊方差均值作為圖像噪聲方差估計值。然而,此法中塊區(qū)選擇會直接影響計算精度,且計算量較大,執(zhí)行效率較低;基于變換域估計法是假設(shè)圖像經(jīng)變換后,較大系數(shù)由圖像信息決定,而較小系數(shù)對應(yīng)圖像噪聲。Donoho等[15]提出以小波高頻系數(shù)絕對值的中值除以 0.674 5作為噪聲標準差為圖像小波分解后高頻對角子帶,D為HH系數(shù))。Uss等[16]對含噪圖像進行正交余弦變換,針對高頻部分估計噪聲。Liu等[17-18]通過對含噪圖像奇異值特性研究,提出了較精確的SVD(singular valuable decomposition,SVD)域高斯噪聲估計算法,然而圖像噪聲和信息不可能完全分離,故圖像類型、噪聲強度和顏色變化均對算法的有效性產(chǎn)生影響。

        狹葉錦雞兒葉切片圖像具有較多紋理、平滑區(qū)域、邊緣等特征,對其去噪需考慮邊緣紋理等細節(jié)的保護,但目前國內(nèi)外對其去噪的研究卻很少。因此,本文在分析總結(jié)前人算法基礎(chǔ)上,首先針對采集的狹葉錦雞兒葉切片圖像采用平滑區(qū)直方圖重構(gòu)和擬合法估計圖像噪聲類型,然后采用本文提出的基于圖像塊的SVD域圖像噪聲估計法實現(xiàn)圖像噪聲強度估計,最后采用本文提出的基于均值濾波和三維塊匹配濾波聯(lián)合去噪算法實現(xiàn)圖像去噪。

        1 圖像獲取及噪聲類型識別

        1.1 圖像獲取

        將經(jīng)過標準固定液 FAA(formalin–acetic acid–alcohol,F(xiàn)AA)固定處理的狹葉錦雞兒葉片制成石蠟切片(軟化—脫水和透明—浸蠟、包埋和修塊—切片、粘片和烘片—染色—封片),具體石蠟切片制作過程參照李正理(1987)和李和平(2009)的方法[19-20]。脫水和透明及染色步驟十分關(guān)鍵,根據(jù)狹葉錦雞兒葉片特點,脫水和透明步驟為:

        染色步驟為:

        其中E為酒精,D為二甲苯。采用YD-1508R輪轉(zhuǎn)切片機進行切片,經(jīng)顏色處理,葉片韌皮部和木質(zhì)部被染成淡紅色,其他結(jié)構(gòu)被染成綠色。葉切片圖像的獲取由YYS-80E型光學(xué)生物顯微鏡上CM1400攝像機拍照獲得,CM1400攝像機采用彩色CMOS逐行掃描圖像傳感器。由于人工制片染色溫度、染液濃度、攝像機曝光增益比和白平衡等因素影響,使得圖像存在色差和噪聲,部分圖像如圖1所示。

        圖1 不同狹葉錦雞兒葉切片色差和含噪圖像(×4)Fig.1 Colour aberration and noise leaf slices images of differentCaragana stenophyllaPojark(×4)

        1.2 圖像噪聲類型估計

        目前多數(shù)情況下,圖像噪聲被近似為加性高斯噪聲處理。然而自然圖像中噪聲是否服從其分布,并未具體驗證,噪聲模型選擇恰當與否直接影響噪聲強度估計和后續(xù)去噪算法的選擇。因此,本文采用平滑區(qū)直方圖重構(gòu)和擬合法對狹葉錦雞兒葉切片圖像噪聲類型進行估計。首先將被測圖像轉(zhuǎn)化為灰度圖像,人工交互選取平滑區(qū)域,選定區(qū)域如圖1(①、②和③)所示。選定區(qū)域亮度值變化較小,灰度幾乎不變,造成直方圖的波動主要由噪聲變化引起。以區(qū)域①為例,此區(qū)域直方圖如圖2a所示。由圖知,直方圖近似服從高斯分布。為驗證直方圖是否服從高斯分布,依選定區(qū)域灰度均值=214,方差s2=9,采用高斯函數(shù)重構(gòu)區(qū)域①直方圖,結(jié)果如圖2b所示。圖2a和2b十分相似,故認為圖像噪聲類型為加性高斯噪聲。為進一步驗證圖像噪聲類型為高斯噪聲,采用高斯函數(shù)對隨機選取的15幅狹葉錦雞兒葉切片圖像平滑區(qū)域直方圖數(shù)據(jù)點進行擬合,擬合優(yōu)度R2均值為0.996,平均均方根誤差RMSE(root mean squared error,RMSE)為 0.144 6,①區(qū)域直方圖高斯曲線擬合結(jié)果如圖2c所示。因此可確定狹葉錦雞兒葉切片圖像噪聲類型為加性高斯噪聲,噪聲幅度概率密度函數(shù)可表示為

        式中z表示灰度值;表示灰度均值; 表示z的標準差,對于均值= 0的加性高斯噪聲,唯一待定的噪聲參數(shù)是 。

        圖2 選定平滑區(qū)域(見圖1a中①)噪聲類型確定Fig.2 Noise-type determining of selected smooth areas(① in Fig. 1a)

        2 圖像噪聲強度估計

        傳統(tǒng)SVD域圖像高斯噪聲強度估計算法對噪聲強度較低,細節(jié)過多或過少的圖像估計精度均不高。這是由于當噪聲強度較低細節(jié)過多時,用來估計噪聲強度的奇異值尾部數(shù)據(jù)可能被噪聲和圖像信號共同控制;當噪聲強度較低細節(jié)過少時,如背景區(qū)域,忽略亮度變化情況下,用來估計噪聲強度的奇異值全部數(shù)據(jù)幾乎都被噪聲控制。可行的解決辦法是根據(jù)噪聲強度合理確定奇異值尾部數(shù)據(jù)長度,保證用來估計噪聲強度的奇異值盡量由噪聲控制。另外該算法對全圖進行 3次奇異分解,對于a′b(a為圖像高度(pix),b為圖像寬度(pix))圖像,此法計算量級為相對較大。因此,本文在SVD域噪聲水平估計法基礎(chǔ)上,為提高算法對小噪聲圖像噪聲強度估計準確性和運算速率,采用基于圖像塊的SVD域圖像噪聲強度估計算法對狹葉錦雞兒葉切片圖像噪聲強度進行準確估計。

        2.1 基于圖像塊的SVD域圖像噪聲強度估計

        SVD域噪聲強度估計法主要針對高斯噪聲,其原理為:含噪圖像奇異分解后,圖像絕大部分信息集中在前面較大少數(shù)奇異值上,噪聲信息分布在所有奇異值上。圖像能量越大,噪聲對前端奇異值影響相對越小,相較圖像信息,噪聲主要影響奇異值尾部數(shù)據(jù),故采用尾部數(shù)據(jù)對噪聲強度進行估計,其表達式為

        式中圖像尺寸為a′b;PM為尾部M個奇異值的平均值;r為奇異值數(shù)量;1s和2s為加入噪聲標準差,二者差值不能過??;P1M和P2M為加入已知強度噪聲圖像的奇異值矩陣尾部M個奇異值的平均值;固定尺寸圖像的線性函數(shù)的斜率 已知,算法詳見文獻[17]和[18]。

        基于圖像塊的SVD域圖像噪聲強度估計算法主要分3步,具體實現(xiàn)過程如圖3所示。

        圖3 基于圖像塊的SVD域狹葉錦雞兒葉切片圖像噪聲強度估計流程圖Fig.3 Flow chart of noise level estimation in SVD domain forCaragana stenophyllaPojarkimages based on image block

        第一步:圖像分塊

        將狹葉錦雞兒葉切片彩色圖像轉(zhuǎn)化為灰度圖像I,然后將圖像等分為4×4大小圖像塊,采用公式(4)計算各塊熵值。由于原SVD域圖像噪聲強度估計算法對細節(jié)較少圖像估計更精確,故取最小熵值圖像塊進行后續(xù)處理,此圖塊細節(jié)最少,隨機性最小。

        式中zmi為灰度隨機變量;L為灰度級,對于8 bit圖像為256;p(zmi)為灰度zmi出現(xiàn)的概率,取值范圍為[0,1]。

        第二步:SVD域圖像噪聲估計

        對熵值最小圖像塊按原SVD域圖像噪聲強度估計算法進行估計,初始M=3r/4。

        第三步:修正M值輸出結(jié)果

        式中0.02這個系數(shù)是經(jīng)驗值,不排除經(jīng)過精細調(diào)節(jié)得到更好的結(jié)果,但對于本文圖像,修正后所得計算結(jié)果已十分準確。

        圖 4為本文算法對狹葉錦雞兒葉切片圖像分塊求熵值示例,圖塊I5熵值最小,僅為3.847 9。采用傳統(tǒng)SVD域噪聲強度估計法對熵值最小圖塊I5進行3次奇異分解,本文算法計算量級為計算速度提升顯著。

        圖4 狹葉錦雞兒葉切片圖像塊熵值示例Fig.4 Entropy to the leaf slices image blocks ofCaragana stenophyllaPojark

        2.2 噪聲強度估計試驗與結(jié)果分析

        2.2.1 試驗方法

        為評價本文所提噪聲強度估計算法精度、穩(wěn)定性和速度,從美國德克薩斯大學(xué)奧斯汀分校圖像測試庫中選取標準無噪圖像。由前述知,狹葉錦雞兒葉切片圖像噪聲類型為加性高斯噪聲,且強度相對較小,故對測試圖像分別加入強度較小加性高斯噪聲,為測試算法適用范圍,同時對加入強度較大加性高斯噪聲圖像進行處理。分別采用本文算法、基于圖像分割的噪聲方差估計法[21]和傳統(tǒng)SVD域噪聲強度估計法對噪聲強度進行估計。試驗采用Windows7旗艦版64位系統(tǒng)、主頻2.40G Hz、8 G內(nèi)存Asus筆記本電腦,軟件采用Matlab R2014a,具體試驗分為2部分:

        1)選取復(fù)雜度、尺寸與狹葉錦雞兒葉切片圖像相似的monarch圖像,分別向monarch圖像中加入標準差 為1、3、5、7、9、20、30、40、50和 60的高斯白噪聲各100次,采用上述3種算法估計噪聲強度均值,以測試3種算法精度、穩(wěn)定性和速度;

        2)為測試算法適應(yīng)能力,選擇復(fù)雜度和尺寸各異圖像100幅,分別向100幅各異圖象中加入標準差 為1、3、5、7、9、20、30、40、50和60的高斯白噪聲1次,采用上述3種算法估計噪聲強度均值,以測試3種算法精度、穩(wěn)定性和速度。

        2.2.2 結(jié)果與分析

        1)monarch圖像加入不同強度高斯白噪聲估計結(jié)果如表1所示。由表1知,當噪聲標準差小于10時,3種算法有效性都會相對降低,本文算法降低程度相對較小。本文算法估計精度總體上優(yōu)于圖像分割和SVD法(僅在噪聲標準差等于1時差于圖像分割噪聲強度估計法);算法穩(wěn)定性與圖像分割噪聲強度估計法相當,優(yōu)于SVD法;由于本文算法采用分塊處理方式,故算法運算速度得到顯著提升。

        表1 不同算法對同一幅圖像疊加100次高斯噪聲估計結(jié)果對比Table 1 Different algorithms estimation results on Gaussian noise averaged over 100 times for same image

        本文算法以最小熵值圖塊(細節(jié)較少)SVD域估計噪聲強度,當圖像噪聲標準差小于10時,噪聲對圖像影響較小,圖塊經(jīng)奇異分解后,奇異值幾乎全被噪聲控制,此時修正的尾部奇異值數(shù)量M大于奇異值總量的4/5;當圖像噪聲標準差大于10時,噪聲對圖像影響較大,此時M取值范圍較寬泛,為奇異值數(shù)量)。由結(jié)果對比知,從精度上看,當噪聲強度標準差低于40時,本文算法精度較高。當噪聲強度標準差高于40時,3種算法精度相當。本文算法估計噪聲的均方根誤差RMSE值為0.314 1,低于其它2種算法,故總體精度較高;從穩(wěn)定性上看,本文算法與圖像分割噪聲強度估計相當,優(yōu)于傳統(tǒng)SVD域噪聲強度估計法;從處理速度上看,本文算法平均耗時約相當于圖像分割噪聲強度估計法耗時的1/3,約相當于SVD域噪聲強度估計法耗時的1/2,速度提升明顯。

        2)對 100幅復(fù)雜度和尺寸各異測試圖像加入不同強度高斯白噪聲估計結(jié)果如表2所示。由表2知,當噪聲標準差小于 10時,本文算法精度較高;當噪聲標準差大于 10時,本文算法精度與圖像分割噪聲強度估計法相當,優(yōu)于 SVD法。從穩(wěn)定性上看,本文算法與圖像分割噪聲強度估計法相當,傳統(tǒng)SVD法穩(wěn)定性較差,當噪聲標準差為60時,傳統(tǒng)SVD法估計值標準差相當于其他 2種算法的10倍左右。從速度上看,本文算法占有明顯優(yōu)勢。

        表2 不同算法對100幅不同圖像疊加高斯噪聲估計結(jié)果對比Table 2 Different algorithms estimation results on Gaussian noise averaged over 100 different images

        2.3 本文算法對狹葉錦雞兒葉切片圖像噪聲估計

        從拍攝的×4狹葉錦雞兒葉切片圖像中隨機選取 100幅,采用本文算法對其進行噪聲強度估計,每幅測10次求均值,測試結(jié)果如圖 5所示。由圖知,圖像噪聲標準差準確的先驗強度標準差估計值,為后續(xù)噪聲去除奠定了基礎(chǔ)。

        圖5 狹葉錦雞兒葉切片圖像噪聲強度估計Fig.5 Noise level estimation onCaraganastenophyllaPojark leaf slices images

        3 圖像去噪

        BM3D算法已經(jīng)十分接近基于圖像塊去噪算法的最優(yōu)界[22-23],經(jīng)擴展后,被廣泛應(yīng)用到其他圖像處理問題中[24-25]。但正如前述,BM3D算法屬于非局部多點型算法中的多模型方法,其計算量和耗時較大,當圖像分辨率較高時耗時更大,制約了該算法的實際應(yīng)用。因此,本文在分析和總結(jié)BM3D算法基礎(chǔ)上,在保證較優(yōu)去噪效果前提下最大限度提高算法處理速度,嘗試把幾何均值濾波(geometric mean filtering,GMF)和BM3D算法相融合,簡稱GMF-BM3D去噪法,實現(xiàn)對狹葉錦雞兒葉切片圖像聯(lián)合去噪。

        3.1 基于GMF-BM3D去噪算法

        GMF濾波器是比較簡單的非線性均值濾波器,其與BM3D去噪法相比,速度較快,但去噪效果差距很大。BM3D算法去噪效果較好,但處理效率較低。BM3D算法基本原理為:首先將圖像I窗口化,將參考塊與其他相似塊形成三維矩陣;然后用閾值收縮變換域系數(shù)對三維矩陣進行線性變換,從而減少噪聲,通過反變換得到所有組中圖像塊估計;最后通過式(6)對重疊估計部分進行加權(quán)平均,實現(xiàn)多估計聚合,得到最終去噪圖像。該算法詳細介紹請參考文獻[6]。

        基于GMF-BM3D去噪算法主要分2步,具體實現(xiàn)過程見圖6。

        圖6 基于GMF- BM3D狹葉錦雞兒葉切片圖像去噪流程圖Fig.6 Denoising flow chart ofCaragana stenophyllaPojarkleaf slices images based on GMF-BM3D

        1)區(qū)域分割

        采用區(qū)域歸并(statistical region merging,SRM)算法[26]將圖像分為 2部分:背景區(qū)域和目標區(qū)域,不同區(qū)域采取不同去噪方法,任意兩區(qū)域歸并條件為

        式中P(R,R¢)為鄰域合并概率,值為0或1;R、R¢為相鄰兩區(qū)域,灰度均值;b(R)表示合并閾值(b(R¢)類似); 為灰度級;分別為區(qū)域R和圖像I像素數(shù);Q為復(fù)雜度參數(shù)(Q越大區(qū)域合并越嚴格,分割結(jié)果越細,反之越寬容,分割結(jié)果越粗);d為修正系數(shù)。

        試驗中采用復(fù)雜度Q=8的SRM算法對狹葉錦雞兒葉切片圖像進行分割,計算圖像各像素點與其 4鄰域間灰度差值,依式(8)對升序排列差值點對進行判斷,將滿足限定條件的區(qū)域歸并標號。根據(jù)狹葉錦雞兒圖像特征,像素數(shù)最大標號區(qū)為背景區(qū)域(如圖6中①),其余部分為目標區(qū)域。

        2)分區(qū)域去噪

        ① 背景區(qū)域去噪

        采用GMF算法對背景區(qū)域進行去噪,為消除邊緣效應(yīng),對整幅圖像做 3×3的幾何均值濾波,然后將特征區(qū)域置為0。由于幾何均值濾波算法屬經(jīng)典算法,運算速度較快,耗時相對較小,可忽略不計;

        ② 目標區(qū)域去噪

        特征區(qū)域采用BM3D法進行去噪(此處使用的BM3D算法已經(jīng)過Lebrun簡化處理),同樣為消除邊緣效應(yīng),舍棄目標區(qū)域像素數(shù)小于 100的小區(qū)域(避免切片雜質(zhì)干擾),對目標區(qū)域最小外接矩形區(qū)進行BM3D濾波處理,然后將矩形框中背景區(qū)域置為0。

        最后將 2個去噪法區(qū)域合并,從而實現(xiàn)對狹葉錦雞兒葉切片圖像GMF-BM3D算法聯(lián)合去噪。

        3.2 圖像質(zhì)量評價標準

        對去噪算法性能客觀評價按有無參考圖像一般分為全參考、半?yún)⒖己蜔o參考法[27]。當存在標準圖像時,通常采用全參考峰值信噪比(peak signal to noise ratio,PSNR)和結(jié)構(gòu)相似性(structural similarity image measurement,SSIM)對算法去噪性能進行評價;當無標準圖像時,通常采用方法噪聲(method noise,MS)和盲/無參考圖像質(zhì)量評價法(blind/referenceless image spatial quality evaluator,BRISQUE)對算法去噪性能進行評價。

        1)PSNR評價法:衡量原始圖像與降噪后圖像對應(yīng)像素灰度值相似程度,取值范圍為[0,100],值越大,圖像降噪效果越好。具體表達式為

        式中fij和fij¢分別為原始圖像和降噪圖像在像素點(i,j)處灰度值;MN表示圖像總個數(shù);L表示圖像最大灰度級,對于8 bits圖像,L=255。

        2)SSIM評價法:衡量原始圖像與降噪后圖像亮度、對比度和結(jié)構(gòu)相似性的度量,取值范圍為[0,1],值越大,圖像降噪效果越好,具體表達式為

        式中X和Y分別表示原始圖像和降噪后圖像相同位置處的圖像信息;Xm和Ym分別表示X和Y圖像的灰度均值,是圖像結(jié)構(gòu)信息亮度的度量;sX和Ys分別表示X和Y的標準差,是圖像結(jié)構(gòu)信息對比度的度量;sXY為X和Y的協(xié)方差,是圖像結(jié)構(gòu)信息結(jié)構(gòu)屬性的度量;N為像素數(shù);C1和C2為常數(shù),通常對于8 bits圖像,L=255。

        3)MS評價法:衡量自然圖像和降噪后圖像差值圖像所含結(jié)構(gòu)信息成分的評價方法,其差值圖像所含圖像特征結(jié)構(gòu)成分越少(成分多少一般需要主觀判斷),圖像降噪效果越好,具體表達式為

        式中u為含噪圖像;Dh濾波算子(h為濾波參數(shù))。

        4)BRISQUE評價法:依據(jù)自然圖像亮度值歸一化趨于高斯分布特性,圖像失真會造成歸一化系數(shù)統(tǒng)計特征的改變,將2個尺度下提取的4個廣義高斯分布特征和32個相鄰系數(shù)相關(guān)性特征利用訓(xùn)練好的SVR模型計算出圖像質(zhì)量。質(zhì)量得分范圍為[0,100],得分越低,圖像質(zhì)量越好[28-29]。算法中支持向量機回歸得到的SVR模型訓(xùn)練采用的是LIVE(Laboratory for image & video engineering)圖像質(zhì)量評價數(shù)據(jù)庫。

        3.3 去噪試驗與結(jié)果分析

        3.3.1 試驗方法

        為檢驗本文去噪算法性能,試驗共分為 2部分,試驗軟件和硬件與噪聲估計使用相同。

        1)從美國德克薩斯大學(xué)奧斯汀分校圖像測試庫中選取與狹葉錦雞兒葉切片復(fù)雜度和尺寸相近標準無噪圖像15幅,向圖象加入標準差為3的高斯白噪聲(狹葉錦雞兒葉切片噪聲強度分別采用傳統(tǒng) BM3D算法(窗口半徑4像素、搜索窗口半徑32像素、滑步3像素)、小波去噪算法(wavelet threshold,WT,使用sym5小波2層分解、軟閾值法、處罰閾值,迭代25次)[30-32]與本文算法對圖像進行降噪處理。使用 PSNR、SSIM、MS和BRISQUE評價法對去噪結(jié)果進行評價。

        2)從基于圖像塊的SVD域圖像噪聲強度估計后的×4狹葉錦雞兒葉切片圖像中隨機選取15幅,使用上述算法對圖像進行降噪處理,采用MS和BRISQUE評價法對去噪結(jié)果進行評價。

        3.2.2 結(jié)果與分析

        1)對15幅標準測試圖像部分處理結(jié)果如圖7所示,幾種去噪算法去噪結(jié)果客觀評價如表 3所示。其中,圖7a為無噪monarch圖像,為更好觀察去噪效果,圖7b為圖7a中A、B兩矩形框區(qū)域的放大顯示,圖7c為A、B兩區(qū)域加入標準差為3的高斯白噪聲圖像,圖7d、7e和7f為各算法去噪結(jié)果和對應(yīng)MS圖像。由圖7d上面2幅圖可以看出,傳統(tǒng)BM3D法能夠在去除噪聲同時很好地保護圖像細節(jié)、紋理和邊緣,對應(yīng)7d下面2幅MS圖像與白噪聲信號十分接近;圖7e上面2幅圖表明,WT算法對含噪圖像具有一定去噪效果,但仍會對圖像紋理細節(jié)等較多部分造成模糊,對應(yīng)7e下面2幅MS圖像中所含邊緣和紋理等信息較多;圖7f表明,本文算法與BM3D算法處理效果相當,很好地兼顧了去噪和細節(jié)、紋理和邊的評價結(jié)論基本一致,去噪后圖像PSNR值越大, SSIM緣保持的優(yōu)點。由表3可以看出,各評價指標對圖像質(zhì)量值越大,對應(yīng)BRISQUE值越小。本文算法平均PSNR值為35.99 dB,與BM3D算法相當,比WT算法高5 dB左右;平均SSIM值為0.82,與BM3D算法相當,比WT算法高0.05左右;平均BRISQUE值為13.12,與BM3D算法相當,比WT算法低4左右。相較傳統(tǒng)BM3D算法,本文算法PSNR值和SSIM值略小,BRISQUE值略大。主要是由于目標區(qū)域采用BM3D算法去噪,背景部分采用GMF法去噪所致,若單獨對目標區(qū)域進行評價,評價結(jié)果應(yīng)與傳統(tǒng)BM3D算法一致。本文算法耗時相對于傳統(tǒng)BM3D算法的1/9。綜合衡量,本文算法較好。

        圖7 本文算法與其他去噪算法對測試圖像monarch圖像去噪結(jié)果示例Fig.7 Example of denoising results of tested images monarch image by proposed algorithm in contrast toother algorithms

        表3 本文算法與其他去噪算法對測試圖像去噪結(jié)果的客觀評價Table 3 Objective evaluation of denoising results for tested images by proposed algorithm and other algorithms

        2)對15幅×4狹葉錦雞兒葉切片圖像部分處理結(jié)果如圖8所示,幾種去噪算法去噪結(jié)果客觀評價如表4所示。其中,圖8a為葉切片原始自然圖像,為更好觀察去噪效果,圖8b為圖8a中A、B兩矩形框區(qū)域的放大顯示,圖8c、8d和8e為各算法去噪結(jié)果和對應(yīng)MS圖像。由圖8知,本文算法去噪效果與 BM3D算法幾乎一致,優(yōu)于WT算法。由表4知,15幅×4狹葉錦雞兒葉切片圖像噪聲標準差?s均值為3.21,受噪聲影響,原圖BRISQUE均值為30.95。經(jīng)本文算法處理后,BRISQUE均值為17.24,約相當于原圖BRISQUE均值的1/2,平均耗時8.04s,約相當于傳統(tǒng)BM3D算法耗時的1/9。本文算法能夠較好實現(xiàn)狹葉錦雞兒葉切片圖像去噪。

        圖8 本文算法與其他去噪算法對切片圖像去噪結(jié)果示例Fig.8 Denoising results of section images by proposed algorithm and other algorithms

        表4 本文算法與其他去噪算法對切片圖像去噪結(jié)果的客觀評價Table 4 Objective evaluation of denoising results for section images by proposed algorithm and other algorithms

        4 結(jié) 論

        本文針對狹葉錦雞兒葉切片圖像進行去噪研究,得出以下結(jié)論:1)平滑區(qū)直方圖重構(gòu)和擬合法能夠?qū)崿F(xiàn)對狹葉錦雞兒葉切片圖像噪聲類型估計,高斯函數(shù)對平滑區(qū)域直方圖數(shù)據(jù)擬合優(yōu)度2R均值高達0.996,平均均方根誤差RMSE(root mean squared error,RMSE)為 0.144 6,故切片圖像噪聲類型為加性高斯噪聲;2)本文提出的基于圖像塊的SVD域噪聲強度估計法,估計過程主要分2步進行,先對熵值最小圖塊 SVD(singularly valuable decomposition,SVD)域估計得到一個初始估計,然后自適應(yīng)的通過初始估計修正算法參數(shù),在第二次估計中得到更準確的估計。處理精度方面,當噪聲標準差小于 10時,本文算法精度最高,當噪聲標準差大于10時,本文算法精度與圖像分割噪聲強度估計法相當,優(yōu)于傳統(tǒng)SVD法。本文算法穩(wěn)定性與圖像分割噪聲強度估計法相當,優(yōu)于傳統(tǒng)SVD法。采用本文算法估計出的狹葉錦雞兒葉切片圖像噪聲標準差本文提出的基于GMF-BM3D聯(lián)合去噪法不僅能較好去除噪聲,而且能完好保留圖像邊緣、紋理等細節(jié)特征。相較傳統(tǒng) BM3D(block-matching and 3-D filtering,BM3D)算法,去噪效果相當,但耗時約相當于BM3D算法的1/9。相較WT(wavelet threshold,WT)算法,本文算法去噪后圖像BRISQUE(block-matching and 3-D filtering,BM3D)值比使用WT法去噪低4左右,圖像質(zhì)量提升明顯。本文算法處理后的狹葉錦雞兒葉切片圖像BRISQUE值僅為10左右,相當于含噪圖像1/2,因此本文采用噪聲類型估計—強度估計—去噪過程及算法對狹葉錦雞兒葉切片圖像去噪取得了較好效果。

        [1]王海超,王春光,賀曉,等.狹葉錦雞兒葉片解剖結(jié)構(gòu)對干旱脅迫下荒漠草原退化梯度的響應(yīng)[J].生態(tài)環(huán)境學(xué)報,2016,25(5):744-751.Wang Haichao, Wang Chunguang, He Xiao, et al. The response ofCaragana stenophyllaPojark leaf anatomical structure to desert grassland deterioration gradient under the drought stress[J]. Ecology and Environmental Sciences, 2016,25(5): 744-751. (in Chinese with English abstract)

        [2]He K, Sun J, Tang X. Guided image filtering[M]. Springer Berlin Heidelberg: Computer Vision-Eccv, 2010: 1-14.

        [3]Buades A, Coll B, Morel J M. A review of image denoising algorithms, with a new one[J]. Multiscale Modeling &Simulation, 2005, 4(2): 490-530.

        [4]Buades A, Coll B, Morel J M. A non-local algorithm for image denoising[C]//IEEE Computer Society Conference on Computer Vision and Pattern Recognition. IEEE Computer Society, 2005: 60-65.

        [5]Silva R D, Minetto R, Schwartz W R, et al. Adaptive edge-preserving image denoising using wavelet transforms[J].Pattern Analysis and Applications, 2013, 16(4): 567-580.

        [6]Dabov K, Foi A, Katkovnik V, et al. Image denoising by sparse 3-D transform-domain collaborative filtering[J]. IEEETransactions on Image Processing, 2007, 16(8): 2080-2095.

        [7]Dabov K, Foi A, Katkovnik V, et al. BM3D image denoising with shape-adaptive principal component Analysis[C]//Proc Workshop on Signal Processing with Adaptive Sparse Structured Representation Saint-malo, 2009.

        [8]Lebrun M. An analysis and implementation of the BM3D image denoising method[J]. Image Processing Online, 2012,2(25): 175-213.

        [9]Wen Y, Ng M K, Huang Y. Efficient total variation minimization methods for color image restoration[J]. IEEE Transactions on Image Processing, 2008, 17(11): 2018-2088.

        [10]Tang Q X, Jiao L C. Image denoising with geometrical thresholds[J]. Electronics Letters, 2009, 45(8): 405-406.

        [11]Portilla J, Strela V, Wainwright M J. Image denoising using scale mixtures of Gaussians in the wavelet domain[J]. IEEE Transactions on Image Processing, 2003, 12(11): 1338-1351.

        [12]Tai S C, Yang S M. A design framework for hybrid approaches of image noise estimation and its application to noise reduction[J]. Journal of Visual Communication &Image Representation, 2012, 23(5): 812-826.

        [13]Salmeri M, Mencattini A, Ricci E, et al. Noise estimation in digital images using fuzzy processing[C]//International Conference on Image Processing, 2001. Proceedings. IEEE,2001: 517-520.

        [14]Tian J, Chen L. Image noise estimation using a variationadaptive evolutionary approach[J]. IEEE Signal Processing Letters, 2012, 19(7): 395-398.

        [15]Donoho D L, Johnstone J M. Ideal spatial adaption by wavelet shrinkage[J]. Biometrika, 1994, 81(3): 425-455.

        [16]Uss M, Vozel B, Lukin V. Image informative maps for estimating noise standard deviation and texture parameters[J].Eurasip Journal on Advances in Signal Processing, 2011(1):1-12.

        [17]Liu W, Lin W. Additive white Gaussian noise level estimation in SVD domain for images[J]. IEEE Transactions on Image Processing A Publication of the IEEE Signal Processing Society, 2012, 22(3): 872-883.

        [18]柳薇.SVD域的圖像高斯噪聲強度估計[J].中國圖象圖形學(xué)報,2012,17(8):923-933.Liu Wei. Gaussian noise level estimation in SVD domain for images[J]. Journal of Image and Graphics, 2012, 17(8): 923-933. (in Chinese with English abstract)

        [19]李正理.植物制片技術(shù)(第二版)[M].北京:科學(xué)出版社,1987:3-60.

        [20]李和平.植物顯微技術(shù)(第二版)[M].北京:科學(xué)出版社,2009:2-47.

        [21]王志明.基于圖像分割的噪聲方差估計[J].工程科學(xué)學(xué)報,2015,9(37):1218-1224.Wang Zhiming. Noise variance estimation based on image segmentation[J]. Chinese Journal of Engineering, 2015, 9(37):1218-1224. (in Chinese with English abstract)

        [22]Chatterjee P, Milanfar P. Patch-based near-optimal image denoising[J]. IEEE Trans on Image Processing, 2012, 21(4):1635-1649.

        [23]Levin A, Nadler B. Natural image denoising: Optimality and inherent bounds[C]//IEEE Conference on Computer Vision and Pattern Recognition. IEEE Computer Society, 2011:2833-2840.

        [24]Danielyan A, Katkovnik V, Egiazarian K. BM3D frames and variational image deblurring[J]. IEEE Transactions on Image Processing: A Publication of the IEEE Signal Processing Society, 2012, 21(4): 1715-1728.

        [25]Maggioni M, Katkovnik V, Egiazarian K, et al. Nonlocal transform-domain filter for volumetric data denoising and reconstruction[J]. IEEE Transactions on Image Processing: A Publication of the IEEE Signal Processing Society, 2013,22(1): 119-133.

        [26]Nock R, Nielsen F. Statistical region merging[J]. IEEE Transactions on Pattern Analysis & Machine Intelligence,2004, 26(11): 1452-1458.

        [27]王志明.無參考圖像質(zhì)量評價綜述[J].自動化學(xué)報,2015,6(41):1062-1079.Wang Zhiming. Review of no-reference image quality assessment[J]. Acta Automatica Sinica, 2015, 6(41): 1062-1079. (in Chinese with English abstract)

        [28]Mittal A, Moorthy A K, Bovik A C. Blind/referenceless image spatial quality evaluator[C]// Signals, Systems and Computers, 2011: 723-727.

        [29]Mittal A, Moorthy A K, Bovik A C. No-reference image quality assessment in the spatial domain[J]. IEEE Transactions on Image Processing, 2012, 21(12): 4695-4708.

        [30]李麗,張楠楠,梅樹立,等.基于自適應(yīng)小波偏微分方程的蝗蟲切片圖像去噪[J].農(nóng)業(yè)工程學(xué)報,2015,31(20):172-177.Li Li, Zhang Nannan, Mei Shuli, et al. Image de-noising of locust sections based on adaptive wavelet and partial differential equation method[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(20): 172-177. (in Chinese with English abstract)

        [31]曹晏飛,滕光輝,余禮根,等.含風(fēng)機噪聲的蛋雞聲音信號去噪方法比較[J].農(nóng)業(yè)工程學(xué)報,2014,30(2):212-218.Cao Yanfei, Teng Guanghui, Yu Ligen, et al. Comparison of different de-noising methods in vocalization environment of laying hens including fan noise[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(2): 212-218. (in Chinese with English abstract)

        [32]楊福增,田艷娜,楊亮亮,等.基于雜交小波變換的農(nóng)產(chǎn)品圖像去噪算法[J].農(nóng)業(yè)工程學(xué)報,2011,27(3):172-178.Yang Fuzeng, Tian Yanna, Yang Liangliang, et al. Agricultural image de-noising algorithm based on hybrid wavelet transform[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2011,27(3): 172-178. (in Chinese with English abstract)

        Blind image denoising of microscopic slices image ofCaragana stenophylla Pojarkbased on noise type and intensity estimation

        Wang Haichao, Wang Chunguang※, Zong Zheying, Yin Xiaofei, Zhang Wenxia,Wang Xiaorong, Zhang Haijun, Li Jing, Liu Tao
        (College of Mechanical and Electrical Engineering,Inner Mongolia Agricultural University, Hohhot010018,China)

        The microscopic images from leaf slices ofCaragana stenophylla Pojarkare inevitably corrupted by noise in obtaining, which will have a negative effect on its subsequent processes. The unpitched sound is usually treated as an additive Gaussian sound in order to denoise the natural images. It has not yet been known whether the unpitched sound and the additive Gaussian sound display the same distribution. The choice of the unpitched sound model will determine the estimation on its noise level and its subsequent choice of the denoising algorithm. BM3D (block-matching and three-dimensional filtering) is universally regarded as the best denoising algorithm, since it keeps the maximum information about texture, border, and other details of the image while denoising effectively. But the speed of BM3D is limited owing to too much calculation and its practical use is also limited because it needs the noise variation as a pretested input. SVD (singular valuable decomposition) is one of the best denoising algorithms but the noise of the image and its information can not be completely separated. The type of image, the level of the noise, and the change of color all affect the validity of the algorithm. This paper proposed a process,involving the estimation on the image type, the noise level and the algorithm of denoising, to solve the problems mentioned above and to denoise the microscopic images from leaf slices ofCaragana stenophylla Pojark. The noise type was recognized as Gaussian noise by the rebuilding and fitting of the smooth histogram domain. The values, abstracted from the smooth histogram domain of 15 randomly selected microscopic images from leaf slices ofCaragana stenophylla Pojark, were fitted with the Gaussian function, with the R2 value of 0.996 and the RMSE (root mean squared error) value of 0.144 6. And then the standard deviation of the noise was estimated via the estimation on noise level with the help of SVD domain image block. In the algorithm estimation, an initial estimative value was obtained from the SVD domain image block with the smallest entropy,and a more accurate value was got from its subsequent self-adapting correcting of the algorithm values. A two-step simulation experiment proved that this algorithm was much better in speed and accuracy in denoising than the other two estimation algorithms. When the standard deviation of noise was lower than 10, the proposed algorithm had the best accuracy; when the standard deviation of noise was higher than 10, this algorithm had the same accuracy with the estimation method of sliced image noise level, but had a better accuracy than SVD. This paper got a standard deviation of the image noise of between 2.5 and 4.0 for its estimation on the microscopic images from leaf slices ofCaragana stenophylla Pojark. The image was denoised via GMF-BM3D, a joint algorithm by GMF (geometric mean filtering) and BM3D. GMF was used to denoise the background of images and then BM3D was applied to denoise the target area. Experiments proved the algorithm introduced in this paper was good at denoising and remaining detailed information of the image, such as border and texture, with a much greater speed.The average value of the BRISQUE (blind/referenceless image spatial quality evaluator) of the processed image was 10,equivalent to half of that of the original image. The algorithm also had an equivalent effect in denoising with the traditional BM3D, but eight ninths of time was saved. It was better than WT (wavelet threshold) denoising in BRISQUE, which was lowered by about 4. Therefore, the algorithm in this paper accurately denoises the microscopic images from leaf slices ofCaragana stenophylla Pojarkand provides a reliable technical support for its subsequent image processing.

        image processing; Gaussian noise; estimation; slices images;Caragana stenophyllaPojark; noise level; image denoising

        10.11975/j.issn.1002-6819.2017.10.030

        TP391.41

        A

        1002-6819(2017)-10-0229-10

        2016-10-21

        2017-04-23

        內(nèi)蒙古自治區(qū)高等學(xué)??茖W(xué)研究項目(NJZY070);內(nèi)蒙古自治區(qū)博士研究生科研創(chuàng)新項目(B20151012902Z);內(nèi)蒙古自治區(qū)自然科學(xué)基金項目(2017MS0514);內(nèi)蒙古自治區(qū)自然科學(xué)基金項目(2017MS0361)

        王海超,男,河北承德人,博士生,研究方向為荒漠草原典型植物切片圖像特征與草原早期退化相關(guān)性研究。呼和浩特 內(nèi)蒙古農(nóng)業(yè)大學(xué)機電工程學(xué)院,010018。Email:wanghaichao1129@163.com

        ※通信作者:王春光,男,內(nèi)蒙古鄂爾多斯人,教授,博士生導(dǎo)師,研究方向為圖像與數(shù)字化研究。呼和浩特 內(nèi)蒙古農(nóng)業(yè)大學(xué)機電工程學(xué)院,010018。

        Email:jdwcg@imau.edu.cn

        猜你喜歡
        錦雞兒標準差高斯
        小高斯的大發(fā)現(xiàn)
        用Pro-Kin Line平衡反饋訓(xùn)練儀對早期帕金森病患者進行治療對其動態(tài)平衡功能的影響
        8種野生錦雞兒在烏魯木齊的引種試驗
        防護林科技(2020年5期)2020-07-21 01:53:32
        17種錦雞兒屬植物葉片解剖結(jié)構(gòu)及抗旱性分析
        天才數(shù)學(xué)家——高斯
        不同溫度對變色錦雞兒種子發(fā)芽特征的影響
        錦雞兒屬植物種子性狀和幼苗生長特征比較研究
        對于平均差與標準差的數(shù)學(xué)關(guān)系和應(yīng)用價值比較研究
        有限域上高斯正規(guī)基的一個注記
        醫(yī)學(xué)科技論文中有效數(shù)字的確定
        国产人成视频免费在线观看| 天天躁狠狠躁狠狠躁夜夜躁| 欧美日韩久久久精品a片| 亚洲狼人社区av在线观看| 人妻秘书被社长浓厚接吻| 亚洲一区二区女搞男| 日韩人妻无码一区二区三区| 一本无码av一区二区三区| 麻豆av在线免费观看精品| 久久精品国产亚洲av麻豆瑜伽| 久久99精品久久久久久琪琪| 变态 另类 欧美 大码 日韩| 成人动漫久久| 色av色婷婷18人妻久久久| 亚洲视频在线中文字幕乱码| 爆操丝袜美女在线观看| 国产在线精品一区二区三区| 亚洲国产精品久久亚洲精品| 九九99久久精品午夜剧场免费 | 精品国产麻豆一区二区三区| 日韩精品一区二区三区乱码| 末成年女a∨片一区二区| 亚洲免费视频播放| 国产精品很黄很色很爽的网站| 久久精品亚洲成在人线av乱码| 亚洲中文字幕无码爆乳app| 国产成人无码精品久久二区三区| 97在线观看播放| 亚洲h视频| 久久精品国产91久久性色tv | 精品不卡视频在线网址| 隔壁老王国产在线精品| 日日噜噜夜夜狠狠久久无码区 | 欧美成人免费观看国产| 国产精品久久婷婷六月| 性欧美丰满熟妇xxxx性久久久| 最近中文字幕大全在线电影视频| 亚洲AV成人无码久久精品老人| 无码免费午夜福利片在线| av在线免费观看麻豆| 国产日产欧洲系列|