• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    結(jié)合局部頻率估計(jì)的小波域InSAR相位濾波新方法

    2015-06-05 15:33:40李芳芳胡東輝丁赤飚
    關(guān)鍵詞:子帶條紋小波

    李芳芳,林 雪,胡東輝,丁赤飚

    (1.中國(guó)科學(xué)院電子學(xué)研究所,北京100190;2.中國(guó)科學(xué)院空間信息處理與應(yīng)用系統(tǒng)技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京100190)

    結(jié)合局部頻率估計(jì)的小波域InSAR相位濾波新方法

    李芳芳1,2,林 雪1,2,胡東輝1,2,丁赤飚1,2

    (1.中國(guó)科學(xué)院電子學(xué)研究所,北京100190;2.中國(guó)科學(xué)院空間信息處理與應(yīng)用系統(tǒng)技術(shù)重點(diǎn)實(shí)驗(yàn)室,北京100190)

    干涉相位噪聲的存在直接影響相位解纏的效果及最終干涉測(cè)量的精度。針對(duì)這一問(wèn)題,本文提出了一種結(jié)合局部頻率估計(jì)的小波域干涉合成孔徑雷達(dá)(interferometric synthetic aperture radar,InSAR)相位濾波方法。根據(jù)局部頻率估計(jì)可以判斷出干涉相位小波系數(shù)中包含有用信息的子帶。結(jié)合VisuShrink和NeighShrink兩種小波閾值收縮方法分別具有去噪效果好和細(xì)節(jié)保持能力強(qiáng)的特點(diǎn),對(duì)有用信息所在子帶的小波系數(shù)利用NeighShrink方法進(jìn)行閾值收縮,而對(duì)其他子帶的小波系數(shù)則利用VisuShrink方法進(jìn)行閾值收縮,從而盡可能地濾除噪聲,同時(shí)保持干涉條紋的細(xì)節(jié)信息不被破壞。仿真和實(shí)測(cè)數(shù)據(jù)實(shí)驗(yàn)驗(yàn)證了本文提出的濾波方法的有效性。

    干涉合成孔徑雷達(dá);相位濾波;小波變換;局部頻率估計(jì);閾值收縮

    0 引 言

    合成孔徑雷達(dá)(synthetic aperture radar,SAR)干涉測(cè)量技術(shù)通過(guò)從同一場(chǎng)景多次觀測(cè)的SAR復(fù)圖像中獲取相位信息,從而反演地表的高程信息和地形變化信息[1]。干涉測(cè)量的精度和可靠性在很大程度上取決于干涉相位圖的質(zhì)量。然而,在實(shí)際系統(tǒng)中,受熱噪聲去相干、時(shí)間去相干、基線去相干、信號(hào)處理去相干等多種去相干因素的影響,干涉相位噪聲難以避免[2]。相位噪聲不僅給相位解纏帶來(lái)困難,而且影響最終的干涉測(cè)量精度。因此,在相位解纏前必須對(duì)干涉相位進(jìn)行濾波,從而獲得較為準(zhǔn)確的干涉相位估計(jì)值。

    目前,干涉相位的濾波方法大體可以分為空間域?yàn)V波和頻率域?yàn)V波兩類。圓周期均值或中值濾波[3]是最常用的空間域?yàn)V波方法,它實(shí)現(xiàn)簡(jiǎn)單,運(yùn)算速度快,但是濾波窗口難以確定,在條紋密集時(shí)容易破壞相位細(xì)節(jié),降低分辨率。在此基礎(chǔ)上,發(fā)展出了基于局部頻率補(bǔ)償?shù)淖赃m應(yīng)濾波方法[45]和自適應(yīng)選擇窗口大小和方向的濾波方法[67]。近年來(lái),發(fā)展出了基于非局部的濾波方法[89],并有學(xué)者將其應(yīng)用到了干涉相位處理中[1011],與傳統(tǒng)基于局部窗口的濾波方法不同,該方法基于整幅圖像,利用圖像的紋理相似性來(lái)濾除噪聲,可以更好地保留圖像細(xì)節(jié),但該方法的計(jì)算復(fù)雜度過(guò)高,運(yùn)算效率較低。頻率域?yàn)V波方法中,Goldstein濾波[12]應(yīng)用最為廣泛,但該方法受分塊大小和濾波參數(shù)的影響較大,在信噪比很低時(shí),濾波效果較差。另外,文獻(xiàn)[13]提出了基于信號(hào)子空間投影的干涉相位估計(jì)方法,可以在配準(zhǔn)存在誤差的情況下獲得較好的估計(jì)結(jié)果,然而在相干性較低時(shí),噪聲子空間的維數(shù)難以判斷,影響估計(jì)的準(zhǔn)確性。

    小波變換由于其良好的時(shí)頻分析特性和多分辨率特性,在信號(hào)分析和圖像處理領(lǐng)域具有廣泛的應(yīng)用。Carlos等人提出利用離散小波變換(discrete wavelet transform,DWT)進(jìn)行干涉相位濾波[14],通過(guò)增大小波系數(shù)中的信號(hào)成分來(lái)實(shí)現(xiàn),該方法能夠很好地保持干涉條紋的細(xì)節(jié)信息,并在一定程度上提高了圖像的信噪比。但由于對(duì)噪聲信息沒(méi)有進(jìn)行有效的抑制,使得去噪效果不夠好。因此,本文提出一種結(jié)合局部頻率估計(jì)和小波閾值收縮的干涉相位濾波方法,通過(guò)局部頻率估計(jì)結(jié)果區(qū)分干涉相位小波變換后條紋信息及噪聲所在的子帶,分別利用NeighShrink和VisuShrink方法進(jìn)行閾值收縮,從而在濾除噪聲的同時(shí)保持干涉條紋的細(xì)節(jié)信息不被破壞。仿真和實(shí)測(cè)數(shù)據(jù)實(shí)驗(yàn)驗(yàn)證了本文算法的有效性和適應(yīng)性。

    本文的內(nèi)容安排如下:第1節(jié)對(duì)小波域中的干涉相位模型進(jìn)行簡(jiǎn)單介紹;第2節(jié)闡述了兩種小波閾值收縮方法的閾值選取原則及各自的優(yōu)缺點(diǎn);第3節(jié)在介紹局部頻率估計(jì)方法的基礎(chǔ)上,給出了本文提出的新的濾波方法流程;第4節(jié)利用仿真和實(shí)際數(shù)據(jù)對(duì)本文提出的干涉相位濾波方法進(jìn)行了驗(yàn)證;最后對(duì)文章進(jìn)行了總結(jié)。

    1 小波域干涉相位模型

    在實(shí)數(shù)域中,干涉相位噪聲符合加性噪聲模型的特征[7],可以表示為

    式中,φz為實(shí)際的干涉相位;φx為理想的干涉相位;v為相位噪聲,均值為零。受三角測(cè)量的限制,干涉相位纏繞在(-π,π]之間,為了避免相位纏繞可能導(dǎo)致的誤差,大多數(shù)干涉相位濾波算法在復(fù)數(shù)域進(jìn)行,即有

    根據(jù)式(1)可知,式(2)的實(shí)部和虛部[14]可以分別表示為

    由此可見(jiàn),Nc值與干涉相干性有關(guān)。當(dāng)相干性高時(shí),小波系數(shù)中的有用相位信息應(yīng)占主導(dǎo),此時(shí)由于Nc值較大,故小波系數(shù)的幅度值也較大;當(dāng)相干性差時(shí),小波系數(shù)中的噪聲占主要成分,由于Nc值較小,故小波系數(shù)的幅度值也較小。因此,通過(guò)設(shè)置閾值可以區(qū)分相位信息和噪聲對(duì)應(yīng)的小波系數(shù),分別進(jìn)行處理,從而實(shí)現(xiàn)干涉相位濾波[15]。

    小波變換的多分辨率特性使得不同的子帶對(duì)應(yīng)不同的頻率范圍。圖1給出了尺度為2時(shí)的二維DWT示意圖,在進(jìn)行第i尺度的變換后,各個(gè)子帶對(duì)應(yīng)的頻率范圍如下:

    圖1 尺度為2的二維DWT示意圖

    2 小波閾值收縮方法

    小波變換具有將信號(hào)或圖像進(jìn)行稀疏表達(dá)的能力。也就是說(shuō),含噪聲的信號(hào)或圖像經(jīng)過(guò)小波變換后,由有用信息產(chǎn)生的小波系數(shù)僅集中在較少的位置和尺度上,且通常具有較大的幅值;而由噪聲產(chǎn)生的小波系數(shù)在時(shí)頻空間分布較廣,但幅值很小。小波閾值收縮算法就是基于這一特征,通過(guò)抑制低于一定閾值的小波系數(shù),而保留高于該閾值的小波系數(shù),從而達(dá)到去除噪聲的目的。選擇合適的閾值對(duì)于小波閾值收縮算法至關(guān)重要,本文將結(jié)合VisuShrink和NeighShrink兩種閾值收縮算法進(jìn)行干涉相位濾波,下面分別介紹其閾值選取原則。

    2.1 VisuShrink方法

    由Donoho和Johnstone提出的VisuShrink方法是目前最為經(jīng)典的閾值收縮方法[16],也稱為通用閾值收縮方法。對(duì)于二維圖像而言,該方法的實(shí)現(xiàn)過(guò)程為:首先通過(guò)小波變換獲得噪聲圖像的小波系數(shù),其中i為小波尺度,m,n為小波系數(shù)的位置。然后根據(jù)式(8)計(jì)算通用閾值T:

    式中,N為信號(hào)長(zhǎng)度或圖像尺寸;σ為噪聲標(biāo)準(zhǔn)差,按式(9)進(jìn)行估計(jì)

    由于高頻子帶上的小波系數(shù)主要由噪聲引起,因此這里僅利用第一層小波分解后的對(duì)角高頻子帶上的小波系數(shù)估計(jì)噪聲標(biāo)準(zhǔn)差。得到通用閾值后,利用軟閾值函數(shù)對(duì)小波系數(shù)進(jìn)行閾值收縮,即有

    式中,sgn(·)為符號(hào)函數(shù),下標(biāo)“+”表示保持正值不變,將負(fù)值置零。最后,將收縮后的小波系數(shù)進(jìn)行小波重構(gòu)得到去噪后的圖像。

    VisuShrink方法實(shí)現(xiàn)簡(jiǎn)單,能夠有效地去除噪聲。然而,由于通用閾值較大,該方法容易導(dǎo)致圖像的過(guò)度平滑,使細(xì)節(jié)保持效果不好。而且VisuShrink方法僅對(duì)高頻子帶的小波系數(shù)進(jìn)行閾值收縮,而保留低頻子帶的小波系數(shù),因而難以選擇合適的小波變換尺度。當(dāng)小波變換尺度較小時(shí),低頻子帶中的噪聲無(wú)法濾除,小波變換尺度過(guò)大時(shí),又容易破壞細(xì)節(jié)。

    2.2 NeighShrink方法

    小波變換具有能力集中和系數(shù)聚簇的特點(diǎn)。也就是說(shuō),一個(gè)小波系數(shù)若包含有用信息,具有較大幅值時(shí),則其鄰域的小波系數(shù)包含有用信息的可能性也很大。根據(jù)小波變換的這一性質(zhì),文獻(xiàn)[17]提出了利用鄰域信息的小波系數(shù)收縮方法。文獻(xiàn)[18]將其擴(kuò)展到二維,用于圖像去噪,提出了NeighShrink方法。該方法以當(dāng)前要處理的小波系數(shù)為中心,選取大小合適的窗口,令

    小波系數(shù)的收縮方法為

    式中,T為通用閾值,計(jì)算方法見(jiàn)式(8)。

    由此可見(jiàn),與VisuShrink方法在閾值處理時(shí)僅考慮當(dāng)前待處理的小波系數(shù)不同,NeighShrink方法充分利用了其周圍小波系數(shù)的分布特點(diǎn)并將其引入收縮策略,具有比VisuShrink方法較低的閾值,能夠更好地保持圖像的細(xì)節(jié)。選擇的鄰域窗口越大,收縮閾值越低,細(xì)節(jié)保持效果更好,但去噪能力相應(yīng)會(huì)下降。為保證噪聲濾除的效果,通常選取的窗口大小為3×3或者5×5。

    3 結(jié)合局部頻率估計(jì)的小波域干涉相位濾波方法

    根據(jù)上節(jié)對(duì)兩種閾值收縮方法的分析可以看出,VisuShrink方法的閾值較高,去噪能力較好,而NeighShrink方法考慮了當(dāng)前小波系數(shù)的鄰域信息,閾值低于VisuShrink方法,細(xì)節(jié)保持能力更強(qiáng)。根據(jù)以上特點(diǎn),本文希望結(jié)合VisuShrink和NeighShrink方法進(jìn)行干涉相位濾波,即對(duì)于干涉相位有用信息所在頻率子帶的小波系數(shù)利用NeighShrink方法進(jìn)行閾值收縮,而對(duì)于其他主要包含噪聲的頻率子帶的小波系數(shù)則利用VisuShrink方法進(jìn)行閾值收縮,這樣可以在盡可能地濾除噪聲的同時(shí)保持干涉條紋的細(xì)節(jié)信息不被破壞。因此,在進(jìn)行小波系數(shù)閾值收縮之前,需要首先確定有用信息所在的子帶。局部頻率估計(jì)可以獲得干涉條紋的頻率范圍,根據(jù)式(7)給出的小波變換后不同子帶的頻率范圍即可確定干涉條紋信息所在的子帶。因此本節(jié)首先介紹局部頻率估計(jì)的方法,之后給出本文提出的濾波處理算法流程。

    3.1 局部頻率估計(jì)

    干涉相位局部頻率估計(jì)的方法有多種。其中,最大似然頻率估計(jì)方法是應(yīng)用最為廣泛的方法之一,它可以通過(guò)二維快速傅里葉變換(fast Fourier transform,F(xiàn)FT)實(shí)現(xiàn)。在復(fù)干涉相位圖的一個(gè)小窗口內(nèi),地形的坡度可以認(rèn)為是一個(gè)常數(shù),即復(fù)干涉相位圖中僅包含一個(gè)主要頻率,可以用一個(gè)二維復(fù)正弦信號(hào)表示。在一個(gè)矩形窗內(nèi),干涉相位可寫為

    φ(m+k,n+l)=2πk·fa+2πl(wèi)·fr+φ(m,n)(13)式中,φ(m,n)為窗口中心像素的相位;k和l為窗口中的像素相對(duì)于中心像素沿方位向和距離向上的偏移;fa和fr分別為方位向和距離向的干涉相位頻率。最大似然頻率估計(jì)方法按下式進(jìn)行

    為保證窗口內(nèi)信號(hào)的平穩(wěn)性,窗口的選擇不能過(guò)大,這樣在具體實(shí)現(xiàn)過(guò)程中,由于二維FFT的量化效應(yīng),會(huì)使估計(jì)結(jié)果的方差較大。因此,在不增加窗口大小的前提下,為了提高估計(jì)精度,可以通過(guò)對(duì)窗口內(nèi)的信號(hào)補(bǔ)零后,再進(jìn)行二維FFT,從而達(dá)到減小頻率采樣間隔的目的。根據(jù)二維局部頻率估計(jì)方差的克拉美羅界[19],可以利用式(15)近似計(jì)算出補(bǔ)零后的窗口大小

    式中,La,Lr為補(bǔ)零后方位向和距離向的窗口長(zhǎng)度,rs/n為干涉相位圖的信噪比,不等式左側(cè)為離散傅里葉變換后輸出頻率分辨率的平方,右側(cè)為局部頻率估計(jì)方差的克拉美羅界。當(dāng)信噪比為0 dB,截取干涉相位圖的窗口大小為9×9時(shí),根據(jù)式(15)可計(jì)算出需要進(jìn)行256×256點(diǎn)的補(bǔ)零FFT。然而,對(duì)干涉相位圖中的每一像素都進(jìn)行如此多點(diǎn)的補(bǔ)零FFT是非常耗時(shí)的。實(shí)際上,只有在頻譜的主瓣內(nèi)需要做內(nèi)插,這樣可以大大提高計(jì)算效率。而Chirp-Z變換(Chirp-Z transform,CZT)能夠在任意指定的頻率范圍內(nèi)調(diào)整頻率采樣間隔。已知x(n),0≤n≤N-1,則它的CZT為

    式中,M表示欲分析的復(fù)頻譜點(diǎn)數(shù),不一定等于N;A和W均為任意復(fù)數(shù),A為復(fù)頻譜的起始取樣點(diǎn),W為取樣點(diǎn)之間的間隔。因此,可以首先利用16×16點(diǎn)的補(bǔ)零FFT找出

    這樣就達(dá)到了256×256點(diǎn)補(bǔ)零FFT的估計(jì)精度,但相比直接補(bǔ)零FFT而言,計(jì)算效率提高了約30倍。

    3.2 濾波算法流程

    根據(jù)上述分析,本文提出了結(jié)合小波閾值收縮和局部頻率估計(jì)的干涉相位濾波方法,算法流程如圖2所示。

    圖2 本文提出相位濾波方法流程圖

    具體實(shí)現(xiàn)步驟如下:

    步驟1利用最大似然方法對(duì)干涉相位圖進(jìn)行局部頻率估計(jì),得到干涉條紋所在的頻率范圍。

    步驟2對(duì)式(2)中的實(shí)部和虛部分別利用DWT進(jìn)行分解,得到各自的小波系數(shù)。根據(jù)式(7)和步驟1的局部頻率估計(jì)結(jié)果,可以判斷出干涉相位中有用信息所在的子帶。

    步驟3對(duì)于有用信息所在子帶的小波系數(shù),利用NeighShrink方法進(jìn)行閾值收縮;而對(duì)于其他主要為噪聲的子帶中的小波系數(shù),則利用VisuShrink方法進(jìn)行閾值收縮。

    步驟4將閾值收縮后的小波系數(shù),利用IDWT進(jìn)行重構(gòu),分別得到濾波后的復(fù)干涉相位的實(shí)部和虛部。

    步驟5計(jì)算出濾波后的干涉相位為=arctan(I︵m/)。

    4 實(shí)驗(yàn)結(jié)果分析

    4.1 仿真實(shí)驗(yàn)

    本小節(jié)利用Matlab中的peaks函數(shù)生成理想的干涉相位圖如圖3(a)所示,數(shù)據(jù)大小為800×800。根據(jù)單視時(shí)的干涉相位概率密度函數(shù)[2],仿真相干系數(shù)為0.85時(shí)的干涉相位如圖3(b)所示。干涉相位方位向和距離向的局部頻率估計(jì)結(jié)果分別如圖3(c)和圖3(d)所示,從中得到干涉相位圖的頻率范圍為[0,0.133π]×[0,0.094π]。實(shí)驗(yàn)中,小波變換的尺度為5,由于第i個(gè)尺度分解后的低頻子帶范圍為[0,2-iπ],因此本文方法對(duì)前兩個(gè)尺度分解后的高頻子帶中的小波系數(shù)均利用VisuShrink方法進(jìn)行閾值收縮,而其他子帶的小波系數(shù)利用NeighShrink方法進(jìn)行閾值收縮,濾波結(jié)果如圖3(i)所示。圖3(e)~圖3(h)分別給出了WInPF濾波[14]、Lee濾波[7]、單獨(dú)用VisuShrink方法及NeighShrink方法閾值收縮的濾波結(jié)果。為了更清楚地顯示濾波后干涉條紋細(xì)節(jié)保持效果,將圖3(e)~圖3(i)中白色方框的區(qū)域進(jìn)行放大顯示,如圖3(j)~圖3(n)所示??梢钥闯觯琖InPF和Lee濾波后,干涉相位圖中還殘余一定的噪聲。VisuShrink方法在條紋密集處一定程度上會(huì)破壞條紋的結(jié)構(gòu),而NeighShrink方法的細(xì)節(jié)保持能力較好但殘余了部分噪聲。相比而言,本文方法獲得了最好的濾波效果。為了進(jìn)一步量化評(píng)價(jià)不同方法的濾波效果,表1給出了不同方法濾波后的殘差點(diǎn)數(shù)目以及濾波后干涉相位與理想干涉相位之間的均方根誤差(root-meam-square error,RMSE),可以看出本文方法較其他方法均有較大的改善。

    另外,表1還給出了不同濾波方法處理該仿真數(shù)據(jù)的運(yùn)行時(shí)間。可以看出,本文方法由于需進(jìn)行局部頻率估計(jì)比WInPF濾波、VisuShrink濾波及NeighShrink濾波運(yùn)行時(shí)間長(zhǎng),但是相比Lee濾波而言運(yùn)算效率有較大提高,能夠處理較大維數(shù)的數(shù)據(jù)。

    表1 不同濾波方法濾波效果及運(yùn)行時(shí)間對(duì)比

    圖3 仿真干涉相位濾波結(jié)果

    4.2 實(shí)測(cè)數(shù)據(jù)實(shí)驗(yàn)

    本小節(jié)利用SIR-C/X-SAR獲取的意大利Etna火山口干涉相位進(jìn)行相位濾波實(shí)驗(yàn)。實(shí)際的干涉相位如圖4(a)所示,數(shù)據(jù)大小為512×512??梢钥闯?,該干涉相位圖條紋非常密集,而且信噪比較低。分別利用WInPF濾波、Lee濾波、VisuShrink濾波、NeighShrink濾波以及本文方法進(jìn)行處理,濾波結(jié)果如圖4(b)~圖4(f)所示??梢钥闯觯琖InPF濾波和Lee濾波的去噪效果較差,Neigh-Shrink方法去噪效果稍好,而VisuShrink方法對(duì)條紋的細(xì)節(jié)信息破壞嚴(yán)重。本文濾波方法的小波變換尺度仍為5,根據(jù)圖4(g)和圖4(h)的頻率估計(jì)結(jié)果得到干涉相位圖的頻率范圍為[0,0.273π]×[0,0.438π],因而對(duì)第一個(gè)尺度分解后的高頻子帶中的小波系數(shù)用VisuShrink方法處理,而其他子帶的小波系數(shù)用NeighShrink方法處理。從圖4(f)可以看出,條紋密集處的干涉相位仍得到了很好的保持。表1顯示了不同方法濾波后干涉相位的殘差點(diǎn)數(shù),可以看出本文方法濾波后殘差點(diǎn)數(shù)目最少。

    圖4 實(shí)際干涉相位濾波結(jié)果

    5 結(jié) 論

    本文提出了一種局部頻率估計(jì)與小波閾值收縮相結(jié)合的干涉相位濾波新方法。該方法針對(duì)VisuShrink和NeighShrink兩種小波閾值收縮算法各自的優(yōu)點(diǎn)和缺陷,將二者結(jié)合使用進(jìn)行干涉相位濾波。該方法的關(guān)鍵在于通過(guò)最大似然估計(jì)得到干涉相位的局部條紋頻率,根據(jù)估計(jì)結(jié)果,對(duì)估計(jì)出的頻率范圍所在子帶的小波系數(shù)利用NeighShrink方法進(jìn)行閾值收縮,而對(duì)其他子帶的小波系數(shù)則利用VisuShrink方法進(jìn)行閾值收縮,從而可以有效地濾除干涉相位噪聲,并保證干涉條紋結(jié)構(gòu)盡量不被破壞。本文方法繼承了基于小波變換的濾波方法效率高的優(yōu)點(diǎn),可以用于較大維數(shù)數(shù)據(jù)的處理。仿真和實(shí)測(cè)數(shù)據(jù)實(shí)驗(yàn)驗(yàn)證了本文方法的有效性。

    [1]Rosen P A,Hensley S,Joughin I R,et al.Synthetic aperture radar interferometry[J].Proceedings of the IEEE,2000,88(3):333 -382.

    [2]Hanssen R F.Radar interferometry,data interpretation and error analysis[M].The Netherlands:Kluwer Academic Publishers,2001.

    [3]Lanari R.Generation of digital elevation models by using SIR-C/X-SAR multifrequency two-pass interferometry:the Etna case study[J].IEEE Trans.on Geoscience and Remote Sensing,1996,34(5):1097- 1114.

    [4]Suo Z,Li Z,Bao Z.A new strategy to estimate local fringe frequencies for InSAR phase noise reduction[J].IEEE Geoscience and Remote Sensing Letters,2010,7(4):771- 775.

    [5]Cai B,Liang D,Dong Z.A new adaptive multiresolution noise-filtering approach for SAR interferometric phase images[J].IEEE Geoscience and Remote Sensing Letters,2008,5(2):266- 270.

    [6]Wu N,F(xiàn)eng D,Li J.A locally adaptive filter of interferometric phase images[J].IEEE Geoscience and Remote Sensing Letters,2006,3(1):73- 77.

    [7]Lee J S,Papathanassiou K P,Ainsworth T L,et al.A new technique for noise filtering of SAR interferometric phase images[J].IEEE Trans.on Geoscience and Remote Sensing,1998,36(5):1456- 1465.

    [8]Buades A,Coll B,Morel J M.Nonlocal image and movie denoising[J].International Journal of Computer Vision,2008,76(2):123- 140.

    [9]Parrilli S,Poderico M,Angelino C V,et al.A non-local SAR image denosing algorithm based on LLMMSE wavelet shrinkage[J].IEEE Trans.on Geoscience and Remote Sensing,2012,50(2):606- 616.

    [10]Deledalle C A,Denis L,Tupin F.NL-InSAR:non-local interferogram estimation[J].IEEE Trans.on Geoscience and Remote Sensing,2011,49(4):1441- 1452.

    [11]Lin X,Li F F,Meng D D,et al.Nonlocal SAR interferometric phase filtering through higher order singular value decomposition[J].IEEE Geoscience and Remote Sensing Letters,2015,12(4):806- 810.

    [12]Goldstein R M,Werner C L.Radar interferogram filtering for geophysical application[J].Geophysical Research Letters,1998,25(21):4035- 4038.

    [13]Li Z,Bao Z,Li H,et al.Image autocoregistration and InSAR interferogram estimation using joint subspace projection[J].IEEE Trans.on Geoscience and Remote Sensing,2006,44(2):288- 297.

    [14]Lopez-Martinez C,F(xiàn)abregas X.Modeling and reduction of SAR interferometric phase noise in the wavelet domain[J].IEEE Trans. on Geoscience and Remote Sensing,2002,40(12):2553- 2566.

    [15]Wang P,Wang Y F,Zhang B C,et al.InSAR interferogram filtering based on Bayesian threshold in stationary wavelet domain[J].Journal of Electronics&Information Technology,2007,29(11):2706- 2710.(汪沛,王巖飛,張冰塵,等.基于貝葉斯門限的靜態(tài)小波域干涉相位圖濾波[J].電子與信息學(xué)報(bào),2007,29(11):2706- 2710.)

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

    [17]Cai T T,Silverman B W.Incorporating information on neighboring coefficients into wavelet estimation[J].The Indian Journal of Statistics,Series B,2001,63(2):127- 148.

    [18]Chen G Y,Bui T D,Krzyak A.Image denoising with neighbor dependency and customized wavelet and threshold[J].Pattern Recognition,2005,38(1):115- 124.

    [19]Apagnolini U.2-D phase unwrapping and instantaneous frequency estimation[J].IEEE Trans.on Geoscience and Remote Sensing,1995,33(3):579- 589.

    Novel InSAR phase filtering method in wavelet domain integrated with local frequency estimation

    LI Fang-fang1,2,LIN Xue1,2,HU Dong-hui1,2,DING Chi-biao1,2
    (1.Institute of Electronics,Chinese Academy of Sciences,Beijing 100190,China;2.Key Laboratory of Technology in Geo-spatial Information Processing and Application System,Chinese Academy of Sciences,Beijing 100190,China)

    Phase noise in the interferogram directly affects the result of phase unwrapping and the precision of interferometirc measurement.A novel approach combining the local frequency estimation with wavelet transform is presented to reduce interferometric phase noise for interferometric synthetic aperture radar(InSAR).According to the local frequency estimation result,the subband of wavelet coefficients containing useful information can be found out.Based on the fact that the VisuShrink method has effective noise reduction ability and the NeighShrink method performs well in detail preservation,the wavelet coefficients containing useful information are shrunk by the NeighShrink method while the others are shrunk by the VisuShrink method.As a result,the noisy interferogram can be efficiently filtered without destroying the fringes.The experiments using the simulated data and the real data demonstrate the validity of the proposed method.

    interferometric synthetic aperture radar(InSAR);phase filtering;wavelet transform;local frequency estimation;threshold shrinkage

    TP 391

    A

    10.3969/j.issn.1001-506X.2015.12.09

    李芳芳(1986- ),女,助理研究員,博士,主要研究方向?yàn)楦缮鍿AR信號(hào)處理。

    E-mail:iecas_lab8@126.com

    林 雪(1986-),女,博士研究生,主要研究方向?yàn)楦缮鍿AR信號(hào)處理。

    E-mail:linxue862002@163.com

    胡東輝(1970-),男,副研究員,碩士,主要研究方向?yàn)镾AR/ISAR成像。

    E-mail:dhhu@m(xù)ails.ie.ac.cn

    丁赤飚(196-9- ),男,研究員,博士,主要研究方向?yàn)檫b感信息處理和應(yīng)用系統(tǒng)。

    E-mail:cbding@m(xù)ails.ie.ac.cn

    1001-506X(2015)12-2719-06

    2014- 12- 02;

    2015- 05- 29;網(wǎng)絡(luò)優(yōu)先出版日期:2015- 07- 07。

    網(wǎng)絡(luò)優(yōu)先出版地址:http://www.cnki.net/kcms/detail/11.2422.TN.20150707.1350.001.html

    國(guó)家自然科學(xué)基金(61401428)資助課題

    猜你喜歡
    子帶條紋小波
    一種基于奇偶判斷WPT的多音干擾抑制方法*
    構(gòu)造Daubechies小波的一些注記
    子帶編碼在圖像壓縮編碼中的應(yīng)用
    電子制作(2019年22期)2020-01-14 03:16:24
    基于MATLAB的小波降噪研究
    電子制作(2019年13期)2020-01-14 03:15:32
    誰(shuí)是窮橫條紋衣服的人
    小福爾摩斯(2019年2期)2019-09-10 07:22:44
    別急!丟了條紋的斑馬(上)
    別急!丟了條紋的斑馬(下)
    基于改進(jìn)的G-SVS LMS 與冗余提升小波的滾動(dòng)軸承故障診斷
    基于虛擬孔徑擴(kuò)展的子帶信息融合寬帶DOA估計(jì)
    條紋,條紋,發(fā)現(xiàn)啦
    丰满人妻熟妇乱又伦精品不卡| 黄色片一级片一级黄色片| 久久精品国产99精品国产亚洲性色 | 亚洲午夜理论影院| 成人影院久久| 日韩大码丰满熟妇| 天堂中文最新版在线下载| 午夜福利欧美成人| 国产熟女xx| 欧美国产精品va在线观看不卡| 9热在线视频观看99| 九色亚洲精品在线播放| 12—13女人毛片做爰片一| 亚洲第一av免费看| 亚洲精品粉嫩美女一区| 久久人妻熟女aⅴ| 久久天躁狠狠躁夜夜2o2o| 免费在线观看亚洲国产| 性欧美人与动物交配| 久久草成人影院| 国内久久婷婷六月综合欲色啪| 亚洲国产欧美网| 在线观看www视频免费| 99久久人妻综合| 精品一区二区三卡| 亚洲精品国产精品久久久不卡| 亚洲七黄色美女视频| 午夜精品在线福利| 国产麻豆69| 欧美丝袜亚洲另类 | av在线播放免费不卡| 亚洲久久久国产精品| 精品久久久久久久毛片微露脸| 亚洲一码二码三码区别大吗| tocl精华| 法律面前人人平等表现在哪些方面| 亚洲av片天天在线观看| 黄网站色视频无遮挡免费观看| 成年人免费黄色播放视频| 国产欧美日韩精品亚洲av| 中文字幕人妻丝袜制服| 久久精品成人免费网站| 亚洲人成电影免费在线| 国产精品 欧美亚洲| 国产成+人综合+亚洲专区| 亚洲一区二区三区色噜噜 | 亚洲欧美一区二区三区黑人| 欧美日韩视频精品一区| 日韩欧美在线二视频| 成人精品一区二区免费| 变态另类成人亚洲欧美熟女 | 国产一卡二卡三卡精品| 欧美日韩中文字幕国产精品一区二区三区 | 老汉色∧v一级毛片| 另类亚洲欧美激情| 亚洲性夜色夜夜综合| 视频在线观看一区二区三区| 欧美一级毛片孕妇| 一级毛片女人18水好多| 国产亚洲精品久久久久久毛片| 村上凉子中文字幕在线| 深夜精品福利| 国产成人精品久久二区二区91| 在线观看免费日韩欧美大片| 精品久久蜜臀av无| 两个人看的免费小视频| 欧美精品啪啪一区二区三区| 国产一区二区三区视频了| 国产xxxxx性猛交| 亚洲精品美女久久久久99蜜臀| 男女下面插进去视频免费观看| 制服诱惑二区| 免费久久久久久久精品成人欧美视频| 99久久99久久久精品蜜桃| 一区二区三区激情视频| 亚洲精品一区av在线观看| 他把我摸到了高潮在线观看| 一区二区日韩欧美中文字幕| 国产成人精品在线电影| 99在线视频只有这里精品首页| 级片在线观看| 伦理电影免费视频| 女人爽到高潮嗷嗷叫在线视频| 久久久久国产精品人妻aⅴ院| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久久av美女十八| 99在线人妻在线中文字幕| 电影成人av| 亚洲人成伊人成综合网2020| 成人18禁高潮啪啪吃奶动态图| 男男h啪啪无遮挡| 亚洲av五月六月丁香网| 国产无遮挡羞羞视频在线观看| 欧美日本亚洲视频在线播放| 黄片小视频在线播放| 国产亚洲精品久久久久5区| 日本撒尿小便嘘嘘汇集6| 怎么达到女性高潮| 老熟妇仑乱视频hdxx| 欧美日韩中文字幕国产精品一区二区三区 | 丝袜美腿诱惑在线| 免费av中文字幕在线| 一边摸一边做爽爽视频免费| 国产三级黄色录像| 日本撒尿小便嘘嘘汇集6| 超碰成人久久| 亚洲,欧美精品.| 国产av又大| av国产精品久久久久影院| 999久久久国产精品视频| 91成年电影在线观看| 欧美中文日本在线观看视频| 在线观看免费视频日本深夜| 琪琪午夜伦伦电影理论片6080| 久久这里只有精品19| 好看av亚洲va欧美ⅴa在| 国产真人三级小视频在线观看| 久久亚洲真实| 真人做人爱边吃奶动态| 日韩欧美一区视频在线观看| www.www免费av| 50天的宝宝边吃奶边哭怎么回事| 精品久久久久久,| cao死你这个sao货| 淫妇啪啪啪对白视频| 97碰自拍视频| 亚洲五月天丁香| 国产精品一区二区在线不卡| 午夜视频精品福利| 日日干狠狠操夜夜爽| 亚洲欧美精品综合久久99| 亚洲精品一卡2卡三卡4卡5卡| 岛国视频午夜一区免费看| 久99久视频精品免费| 韩国av一区二区三区四区| 国产在线观看jvid| 亚洲色图 男人天堂 中文字幕| 日日摸夜夜添夜夜添小说| 91成人精品电影| 久久99一区二区三区| 久久久久精品国产欧美久久久| 国产99白浆流出| 亚洲熟妇熟女久久| 91国产中文字幕| 在线观看一区二区三区| 亚洲国产毛片av蜜桃av| 天堂√8在线中文| 国产黄色免费在线视频| 欧美中文日本在线观看视频| 日韩大尺度精品在线看网址 | 欧美大码av| 亚洲va日本ⅴa欧美va伊人久久| 每晚都被弄得嗷嗷叫到高潮| 国产精品国产av在线观看| 极品教师在线免费播放| 国产熟女午夜一区二区三区| 国产黄色免费在线视频| 性欧美人与动物交配| 91精品国产国语对白视频| 水蜜桃什么品种好| 国产三级在线视频| 成人亚洲精品一区在线观看| 午夜精品久久久久久毛片777| 黑人操中国人逼视频| 男女下面进入的视频免费午夜 | 91精品三级在线观看| av片东京热男人的天堂| 97人妻天天添夜夜摸| 99国产精品免费福利视频| 最新在线观看一区二区三区| 老汉色∧v一级毛片| av在线天堂中文字幕 | 亚洲国产精品sss在线观看 | 欧美精品啪啪一区二区三区| 老汉色av国产亚洲站长工具| 国产激情久久老熟女| 国产99白浆流出| 国内毛片毛片毛片毛片毛片| 亚洲在线自拍视频| av在线播放免费不卡| 成人黄色视频免费在线看| 香蕉久久夜色| 天堂√8在线中文| 可以免费在线观看a视频的电影网站| 国产亚洲av高清不卡| 成人免费观看视频高清| 欧美大码av| 嫩草影视91久久| 日日干狠狠操夜夜爽| 9色porny在线观看| 99re在线观看精品视频| 欧美精品一区二区免费开放| 久久久久国产精品人妻aⅴ院| 亚洲精品美女久久久久99蜜臀| 亚洲成国产人片在线观看| 国产亚洲av高清不卡| 日韩视频一区二区在线观看| 夫妻午夜视频| 乱人伦中国视频| 这个男人来自地球电影免费观看| 免费在线观看视频国产中文字幕亚洲| 免费高清视频大片| 成熟少妇高潮喷水视频| 神马国产精品三级电影在线观看 | 国产精品一区二区免费欧美| 在线观看一区二区三区| 如日韩欧美国产精品一区二区三区| 老司机福利观看| 亚洲三区欧美一区| 日韩欧美一区二区三区在线观看| 18禁裸乳无遮挡免费网站照片 | 性少妇av在线| 亚洲精品中文字幕在线视频| 亚洲国产精品一区二区三区在线| 中文字幕高清在线视频| 免费在线观看完整版高清| 首页视频小说图片口味搜索| 亚洲精品在线美女| 一级片免费观看大全| 黄色毛片三级朝国网站| 99热只有精品国产| 欧美一级毛片孕妇| www.自偷自拍.com| 国产av又大| 久久久久国产精品人妻aⅴ院| 99在线人妻在线中文字幕| 熟女少妇亚洲综合色aaa.| 麻豆久久精品国产亚洲av | 久久精品91蜜桃| 成人亚洲精品一区在线观看| 在线播放国产精品三级| 日本欧美视频一区| 日韩一卡2卡3卡4卡2021年| 变态另类成人亚洲欧美熟女 | 成人国语在线视频| 18禁观看日本| 18禁裸乳无遮挡免费网站照片 | 久久草成人影院| 亚洲成人久久性| 免费搜索国产男女视频| netflix在线观看网站| 亚洲激情在线av| 夜夜夜夜夜久久久久| 亚洲 欧美 日韩 在线 免费| 一级片免费观看大全| 9191精品国产免费久久| 欧美激情极品国产一区二区三区| 久久草成人影院| 日韩免费av在线播放| 18美女黄网站色大片免费观看| 国产熟女xx| 深夜精品福利| 国产高清国产精品国产三级| 男女下面进入的视频免费午夜 | 99精国产麻豆久久婷婷| 日日夜夜操网爽| 精品人妻在线不人妻| 国产精品秋霞免费鲁丝片| 日韩中文字幕欧美一区二区| 午夜福利在线观看吧| 国产成+人综合+亚洲专区| 久久精品国产清高在天天线| 久9热在线精品视频| 好看av亚洲va欧美ⅴa在| 国产av一区二区精品久久| 中国美女看黄片| 国产高清国产精品国产三级| 亚洲中文字幕日韩| 亚洲精品久久成人aⅴ小说| 久久久久久久久久久久大奶| 精品久久久久久成人av| 手机成人av网站| 他把我摸到了高潮在线观看| 1024视频免费在线观看| 91大片在线观看| 久热爱精品视频在线9| 国产av一区二区精品久久| 久久久久九九精品影院| 国产片内射在线| 国产亚洲av高清不卡| 国产精品一区二区三区四区久久 | 午夜成年电影在线免费观看| 国产免费男女视频| 两个人免费观看高清视频| 久久久国产成人免费| 亚洲欧美一区二区三区久久| 国产高清激情床上av| 久久香蕉激情| 久久精品亚洲熟妇少妇任你| 亚洲片人在线观看| tocl精华| 黄色成人免费大全| 国产男靠女视频免费网站| 欧美在线一区亚洲| 99国产精品99久久久久| avwww免费| 国产熟女午夜一区二区三区| 欧美不卡视频在线免费观看 | 亚洲五月天丁香| 免费av中文字幕在线| 久久久久久免费高清国产稀缺| 久久香蕉国产精品| 国产成人欧美在线观看| 成年版毛片免费区| 亚洲精品一二三| 国产精品美女特级片免费视频播放器 | 91精品三级在线观看| 国产精品偷伦视频观看了| 人人妻人人添人人爽欧美一区卜| 12—13女人毛片做爰片一| 女性生殖器流出的白浆| 成人国语在线视频| 色婷婷久久久亚洲欧美| 国产精品国产av在线观看| 午夜精品在线福利| 热re99久久国产66热| 新久久久久国产一级毛片| 满18在线观看网站| 国产亚洲欧美精品永久| 精品少妇一区二区三区视频日本电影| 国产蜜桃级精品一区二区三区| 欧美乱妇无乱码| 美国免费a级毛片| 国产在线观看jvid| 一边摸一边抽搐一进一小说| 国产精品1区2区在线观看.| 亚洲av第一区精品v没综合| 亚洲国产欧美网| 国产一区二区激情短视频| 黑人欧美特级aaaaaa片| 日本精品一区二区三区蜜桃| 制服人妻中文乱码| 18禁国产床啪视频网站| 91大片在线观看| 黄色女人牲交| 搡老乐熟女国产| 两个人免费观看高清视频| 久久久国产一区二区| 搡老岳熟女国产| 久久久久久亚洲精品国产蜜桃av| 村上凉子中文字幕在线| 可以免费在线观看a视频的电影网站| 女人高潮潮喷娇喘18禁视频| 黄色片一级片一级黄色片| 国产精品综合久久久久久久免费 | 桃色一区二区三区在线观看| 无遮挡黄片免费观看| 免费在线观看亚洲国产| 久久人人精品亚洲av| 无人区码免费观看不卡| 精品人妻1区二区| 国产精品秋霞免费鲁丝片| 好男人电影高清在线观看| 99精品久久久久人妻精品| 国产一区二区三区视频了| 91国产中文字幕| 亚洲人成网站在线播放欧美日韩| 日韩大码丰满熟妇| 久久精品91蜜桃| 久久精品亚洲av国产电影网| 999精品在线视频| 91国产中文字幕| 日本黄色视频三级网站网址| 麻豆国产av国片精品| 丝袜在线中文字幕| bbb黄色大片| 成人国产一区最新在线观看| 操出白浆在线播放| 99久久人妻综合| x7x7x7水蜜桃| 久久人人精品亚洲av| 亚洲免费av在线视频| 亚洲少妇的诱惑av| 欧美中文日本在线观看视频| av在线播放免费不卡| 国产亚洲欧美精品永久| 另类亚洲欧美激情| 淫妇啪啪啪对白视频| 亚洲精品美女久久久久99蜜臀| 亚洲欧美日韩高清在线视频| 久久久水蜜桃国产精品网| 99精品在免费线老司机午夜| 淫秽高清视频在线观看| 欧美乱码精品一区二区三区| 亚洲七黄色美女视频| 欧美国产精品va在线观看不卡| 最近最新中文字幕大全免费视频| 午夜激情av网站| 久久精品91无色码中文字幕| 久久久久久免费高清国产稀缺| 香蕉国产在线看| 每晚都被弄得嗷嗷叫到高潮| 久久久精品国产亚洲av高清涩受| 亚洲自拍偷在线| 精品久久蜜臀av无| 最好的美女福利视频网| 波多野结衣一区麻豆| 国产乱人伦免费视频| 女生性感内裤真人,穿戴方法视频| 午夜精品久久久久久毛片777| 免费日韩欧美在线观看| 久久香蕉国产精品| 国产成人啪精品午夜网站| 亚洲全国av大片| 中国美女看黄片| 看片在线看免费视频| 天天躁狠狠躁夜夜躁狠狠躁| 一区在线观看完整版| 9热在线视频观看99| 国产精品爽爽va在线观看网站 | 国产精品久久久av美女十八| 国产欧美日韩综合在线一区二区| 变态另类成人亚洲欧美熟女 | 中国美女看黄片| 一本大道久久a久久精品| 在线天堂中文资源库| 正在播放国产对白刺激| 一级a爱片免费观看的视频| 国产国语露脸激情在线看| 亚洲 国产 在线| 性少妇av在线| 男人舔女人下体高潮全视频| 欧美大码av| 国产一区二区三区综合在线观看| netflix在线观看网站| 久久人妻熟女aⅴ| 9热在线视频观看99| 18禁裸乳无遮挡免费网站照片 | 一级毛片精品| 视频区欧美日本亚洲| 精品一区二区三区av网在线观看| 亚洲欧美一区二区三区久久| 欧美日韩av久久| 99国产综合亚洲精品| 啦啦啦在线免费观看视频4| 久久人人精品亚洲av| bbb黄色大片| 久久久久久免费高清国产稀缺| 人成视频在线观看免费观看| 免费一级毛片在线播放高清视频 | 精品一区二区三卡| 两人在一起打扑克的视频| 黄色 视频免费看| 黄片大片在线免费观看| 亚洲av五月六月丁香网| 国产精品av久久久久免费| 国产99久久九九免费精品| 91麻豆精品激情在线观看国产 | 午夜免费成人在线视频| 亚洲国产精品sss在线观看 | 成人精品一区二区免费| 免费高清视频大片| 91大片在线观看| 美女大奶头视频| 日韩有码中文字幕| 精品人妻1区二区| 欧美日韩亚洲国产一区二区在线观看| 国产精品国产av在线观看| av天堂久久9| 亚洲av成人av| 在线观看舔阴道视频| 免费女性裸体啪啪无遮挡网站| 日韩欧美在线二视频| 国产又色又爽无遮挡免费看| 搡老乐熟女国产| 啦啦啦 在线观看视频| 精品国产国语对白av| 亚洲七黄色美女视频| 在线观看免费视频日本深夜| 精品熟女少妇八av免费久了| 国产成人系列免费观看| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一二三| 国产1区2区3区精品| 欧美日韩亚洲综合一区二区三区_| 9热在线视频观看99| 男人的好看免费观看在线视频 | 久久草成人影院| 80岁老熟妇乱子伦牲交| 99精品在免费线老司机午夜| 香蕉丝袜av| 亚洲午夜精品一区,二区,三区| 高清毛片免费观看视频网站 | 在线观看午夜福利视频| 老熟妇仑乱视频hdxx| 最新在线观看一区二区三区| 欧美午夜高清在线| 日韩大尺度精品在线看网址 | 女人被狂操c到高潮| 男人舔女人下体高潮全视频| 日韩大尺度精品在线看网址 | 97超级碰碰碰精品色视频在线观看| 亚洲第一av免费看| 亚洲在线自拍视频| 一个人免费在线观看的高清视频| 97人妻天天添夜夜摸| 久久精品亚洲av国产电影网| 男女午夜视频在线观看| 精品久久久久久成人av| 黄色视频不卡| 免费在线观看视频国产中文字幕亚洲| 亚洲五月婷婷丁香| 黑人巨大精品欧美一区二区蜜桃| 久久久国产一区二区| 亚洲精品美女久久久久99蜜臀| 久久精品国产综合久久久| 色精品久久人妻99蜜桃| 黄色 视频免费看| 性欧美人与动物交配| 精品国产美女av久久久久小说| 在线观看免费高清a一片| 亚洲欧美一区二区三区黑人| 丰满饥渴人妻一区二区三| 很黄的视频免费| 午夜a级毛片| 日韩大尺度精品在线看网址 | x7x7x7水蜜桃| 男女做爰动态图高潮gif福利片 | 亚洲自偷自拍图片 自拍| 精品福利永久在线观看| 老熟妇仑乱视频hdxx| 首页视频小说图片口味搜索| 亚洲男人天堂网一区| 亚洲av片天天在线观看| 国产精品久久视频播放| 午夜视频精品福利| 亚洲男人的天堂狠狠| 99国产精品免费福利视频| 老司机深夜福利视频在线观看| 国产精品久久视频播放| 久久久久亚洲av毛片大全| 亚洲精品国产一区二区精华液| 成年人免费黄色播放视频| 欧美日韩瑟瑟在线播放| 波多野结衣一区麻豆| 免费少妇av软件| 人人澡人人妻人| 性色av乱码一区二区三区2| 在线观看www视频免费| 亚洲成a人片在线一区二区| 91麻豆精品激情在线观看国产 | 一级片免费观看大全| 自线自在国产av| 午夜久久久在线观看| 国产三级黄色录像| 一边摸一边抽搐一进一出视频| 在线观看日韩欧美| 久热这里只有精品99| av国产精品久久久久影院| 欧美在线一区亚洲| 99久久久亚洲精品蜜臀av| 高清av免费在线| 99在线人妻在线中文字幕| 热99国产精品久久久久久7| 久久久国产成人免费| 国产真人三级小视频在线观看| 国产激情欧美一区二区| 亚洲欧美一区二区三区黑人| 久久婷婷成人综合色麻豆| 高清av免费在线| 欧美日韩中文字幕国产精品一区二区三区 | 久久热在线av| 俄罗斯特黄特色一大片| 999久久久精品免费观看国产| 国产成人欧美在线观看| 亚洲精品美女久久av网站| 老司机在亚洲福利影院| 黄色片一级片一级黄色片| 免费av毛片视频| 午夜日韩欧美国产| 十八禁人妻一区二区| 国产成人系列免费观看| 精品国产乱子伦一区二区三区| cao死你这个sao货| 91字幕亚洲| 正在播放国产对白刺激| 欧美日韩乱码在线| 国产99白浆流出| 亚洲精品国产一区二区精华液| 久久久久久大精品| 一边摸一边做爽爽视频免费| 狂野欧美激情性xxxx| 91精品三级在线观看| 欧美日韩av久久| 18禁黄网站禁片午夜丰满| 1024香蕉在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美人与性动交α欧美软件| tocl精华| 日韩欧美一区视频在线观看| 精品卡一卡二卡四卡免费| 久久性视频一级片| 日本欧美视频一区| 看片在线看免费视频| 亚洲人成电影观看| 亚洲自拍偷在线| 窝窝影院91人妻| 在线观看午夜福利视频| 久久久国产欧美日韩av| 午夜免费鲁丝| 日韩av在线大香蕉| 亚洲全国av大片| 国产人伦9x9x在线观看| 国产伦人伦偷精品视频| 纯流量卡能插随身wifi吗| 国内久久婷婷六月综合欲色啪| 国产精品av久久久久免费| 久久久久久免费高清国产稀缺| 国产精品久久电影中文字幕| 俄罗斯特黄特色一大片| 色哟哟哟哟哟哟| 久久久久精品国产欧美久久久| 精品免费久久久久久久清纯| 日韩免费av在线播放| 丰满迷人的少妇在线观看|