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

    基于交互式分割技術和決策級融合的SAR圖像變化檢測

    2012-01-04 08:01:08萬紅林焦李成辛芳芳
    測繪學報 2012年1期
    關鍵詞:變化檢測鄰域像素

    萬紅林,焦李成,辛芳芳

    西安電子科技大學智能感知與圖像理解教育部重點實驗室,陜西西安710071

    基于交互式分割技術和決策級融合的SAR圖像變化檢測

    萬紅林,焦李成,辛芳芳

    西安電子科技大學智能感知與圖像理解教育部重點實驗室,陜西西安710071

    為免去降斑預處理及克服選擇分布模型的限制,結合差異圖的特點和一種不涉及分布模型的交互式分割方法,產(chǎn)生不同“種子點”下的變化檢測結果后,再利用投票策略進行決策級融合給出最終的變化檢測結果。分割中,將每個像素的特征設置為由差異圖及靜態(tài)小波變換分解差異圖再丟棄高頻系數(shù)后重構得到的各層表示內(nèi),對應位置上的灰度值構成的矢量。此特征及決策級融合的策略使本文變化檢測技術對SAR圖像中的斑點噪聲具有一定的抗差性。在無需對SAR圖像做預處理的情況下,對真實SAR圖像數(shù)據(jù)集的變化檢測結果證實了方法的有效性。

    交互式分割技術;差異圖的特點;決策級融合;靜態(tài)小波變換;SAR圖像;變化檢測

    1 引 言

    合成孔徑雷達(synthetic aperture radar,SAR)圖像的變化檢測,由于數(shù)據(jù)的內(nèi)在復雜性,既要求一個可有效去除斑點噪聲的預處理,又要求一個能有效處理乘性斑點噪聲的分析技術,這使得相對光學圖像以SAR圖像為研究對象的變化檢測技術較少。然而,與光學遙感相比,SAR具有的不受氣候和太陽光照條件影響的特點,使之具有潛在的實際應用價值,比如,SAR能有規(guī)律地監(jiān)測常出現(xiàn)云雨天氣的地區(qū),而光學遙感則易被中斷。目前,SAR圖像的變化檢測技術已成功用于地震災區(qū)的標識[1]、水災監(jiān)控[2]和水稻栽種面積的評估[3]等方面。

    一般地,若檢測由水災、火災或地震造成的突變性變化,只需事件發(fā)生前后兩時相的圖像,而檢測如綠地覆蓋率、濕地變化等逐漸發(fā)生的變化,則需一系列不同時刻獲取的圖像[4]。本文檢測由水災引起的突變性變化,兩時相圖像即可。對已配準、已校正的兩時相SAR圖像,無監(jiān)督的變化檢測技術一般基于以下3步:預處理、圖像比較和閾值分割[5]。預處理的目的是去除斑點噪聲,現(xiàn)有的自適應濾波器,如Lee或Gamma MAP濾波器等,均可用于去除斑點噪聲,但都需事先設定濾波程度。圖像比較是以對應位置有變化發(fā)生的可能性大小作為像素的灰度值生成差異圖。由于SAR圖像所含斑點噪聲的乘性特點,最廣泛使用的方法是比值法[6],但該法只利用了單個像素的信息?,F(xiàn)有基于空間信息的差異圖構造方法,如:交叉互相關測度[7]、區(qū)域均值比等[8],這種比較策略能較好地抑制因輻射等原因引起的灰度值突變的點。本文即采用區(qū)域均值比構造差異圖。

    差異圖生成后,比較經(jīng)典的方法是通過閾值分割差異圖確定變化類和非變化類。一般的,基于閾值的差異圖分析技術在基于像素間相互獨立的假設下,利用單個像素的信息人工[9]或自動[10-13]地給出決策閾值。然而,真實圖像中鄰域像素間具有很大相關性。為利用空間關系,文獻[12]通過馬爾可夫隨機場(Markov random field,MRF)模型在決策階段將其引入,而后,文獻[13]利用鄰域像素間的空間關系糾正廣義高斯模型下EM算法的分類結果。結果表明變化檢測性能有所提高,但與文獻[10—11]中的方法一樣,變化檢測結果仍受所選統(tǒng)計模型與真實分布函數(shù)間擬合度的影響。為避免分布模型的假設,文獻[14]提出基于DCT特征抽取,文獻[15]提出基于紋理分類器的差異圖分析技術。然而,上述方法的變化檢測性能均依賴于預處理,這意味著濾波程度的設定問題仍存在。為此,文獻[5,16—18]利用多尺度分析工具代替濾波器,免去了預處理操作,但文獻[16—18]中的技術仍涉及分布模型的假設,而文獻[5]涉及多個參數(shù)的組合優(yōu)化。

    為了免去預處理及克服分布模型假設的限制,本文借助靜態(tài)小波變換,結合差異圖的特點及交互式的分割技術,給出了一種新的SAR圖像變化檢測方法。利用靜態(tài)小波變換將每個像素的特征設置為對噪聲具有一定抗差性的由灰度值構成的矢量后,再由一種不涉及分布模型的交互式分割技術[19]結合由差異圖的特點給出的“種子點”,并通過像素間的競爭(特征間的比較)將差異圖分成兩類,得到不同“種子點”下的變化檢測結果;而后,再對不同“種子點”的變化檢測結果進行決策級融合,完成變化檢測。該方法無需對已配準、已校正的兩時相SAR圖像作降斑預處理,能產(chǎn)生優(yōu)于文獻[5]中最好方法的變化檢測結果,并克服文獻[10—18]需為類分布選擇合適的統(tǒng)計模型的限制。對真實SAR圖像試驗數(shù)據(jù)的變化檢測結果及比較,說明本文方法優(yōu)于其他相關技術。

    2 差異圖的分析

    對已配準、已校正的兩時刻SAR圖像It1和It2,構造差異圖DA

    式中,N(p)表示位于(i,j)的像素p的鄰域像素構成的集合,p′為其中的一個元素??紤]到窗口越大,噪聲去除得越干凈,但變化區(qū)域擴大得也越厲害,本文依據(jù)實際情況選取大小為3×3的窗口。給定差異圖后,若把變化檢測看成是一個對差異圖二分類的問題,則利用交互式的分割技術生成變化檢測結果是可行的,但存在以下兩個問題:如何給出交互式分割技術所需的“種子點”及克服差異圖中噪聲對其分割性能的影響。對交互式分割技術所需的“種子點”,本文根據(jù)差異圖的灰度分布特點,給出分別屬于變化類和非變化類的子集,并將其作為分割技術所需的種子點;對差異圖中殘留的部分噪聲,通過設置對噪聲具有一定抗差性的特征矢量來減少殘留的部分噪聲對分割性能的影響。最后,由于不同“種子點”對應不同的變化檢測結果,對不同“種子點”下得到的變化檢測結果,采用投票競爭策略對其進行決策級的融合生成最終的變化檢測結果。決策級融合的策略使最終的變化檢測結果對“種子點”的選取不敏感的同時,也加強了對噪聲的抗差性。流程圖如圖1所示。

    圖1 本文方法的框架圖Fig.1 The general scheme of our proposed method

    2.1 交互式分割技術對差異圖的分析

    將交互式分割技術用于分析差異圖DA時,其實質(zhì)為對差異圖中每個像素賦予一個屬性向量后,由指定的“種子點”開始,根據(jù)競爭規(guī)則δ,通過鄰域N逐漸向外擴張,直到所有像素的狀態(tài)子集s不再變化。其中每個像素p的屬性向量由標簽lp、強度θp、特征Vp、鄰域系統(tǒng)N及競爭規(guī)則δ構成,狀態(tài)子集由標簽lp、強度θp、特征矢量Vp構成,即sp={lp,θp,Vp}。對給定的“種子點”,迭代結束后,所有像素的標簽即構成該“種子點”下的變化檢測結果。

    如前所述,迭代前需解決種子點的選取和像素特征的設置問題。對分割方法所需的種子點,根據(jù)差異圖的分布特點[13]給出。利用差異圖的直方圖hD計算出兩個閾值Tn和Tc后,由該閾值確定分別屬于變化類和非變化類的“種子”。若兩閾值為Tn=MD(1-α),Tc=MD(1+α),則變化類的種子有Sc={DA(p)|DA(p)>Tc},非變化類的種子有Sn={DA(p)|DA(p)<Tn},其中MD為直方圖hD的中間值,即MD=[max{DA}-min{DA}]/2;α∈(0,1)是人為設定的一個參數(shù),表征灰度值介于以MD為中心,MD×α為鄰域半徑內(nèi)的像素,其類別是有待確定的。

    關于設置像素p的特征Vp問題,其關鍵在于部分殘留在差異圖中的斑點噪聲對分割技術的性能有影響,因此需設置對噪聲具有一定的抗差性的特征。受文獻[5]啟發(fā),根據(jù)靜態(tài)小波變換對差異圖分解后再重構所得的對差異圖的不同層表示,其內(nèi)噪聲被去除程度不同的特點設置特征Vp。利用靜態(tài)小波變換獲得對差異圖的不同層表示后,將像素p在差異圖及其各層表示內(nèi)對應位置上的灰度值構成的矢量作為特征Vp,以消除斑點噪聲對分割技術性能的影響。靜態(tài)小波變換分解差異圖時,每層都產(chǎn)生4個大小與原圖一致的圖像,其中1個表示低頻信息,3個表示高頻信息。鑒于斑點噪聲都處在高頻系數(shù)內(nèi),將分解后表示高頻信息的變換系數(shù)置0,只利用低頻系數(shù)進行逆靜態(tài)小波變換,得到對差異圖的某個分解層數(shù)下的表示。隨著分解層數(shù)的增加,該層表示內(nèi)噪聲去除的越多,但同時丟失的細節(jié)信息也越多??紤]到噪聲去除和細節(jié)信息保持間的兩難,分別求得對差異圖的1層和2層表示,將特征設置為由三維灰度值構成的矢量。若Dk表示對差異圖的k層表示,則對位置(i,j)上的像素p,其特征矢量Vp為Vp=[DA(i,j) D1(i,j) D2(i,j)]。

    初始時,給定一個α,確定出變化類wc的種子,非變化類wn的種子及類別不確定的像素集Sα后,將差異圖中每個像素的狀態(tài)子集設置為

    并開始迭代。迭代中,標簽不為零的種子點像素,通過鄰域系統(tǒng)N(選取二階鄰域),根據(jù)競爭規(guī)則δ更新差異圖中每個像素的狀態(tài)子集,直到Sα中所有像素的狀態(tài)子集不再變化。其中競爭規(guī)則δ為若當前像素自身的強度大于鄰域N 內(nèi)所有像素的強度,則該像素狀態(tài)子集保持不變;反之,該像素的狀態(tài)子集由其鄰域像素中強度最強的像素狀態(tài)決定,即標簽更新為與其一致,強度更新為最強強度與一個函數(shù)值的乘積,其中函數(shù)值由兩像素特征的相似度確定;若假設像素p的鄰域內(nèi)強度最強像素為q,則lp=lq,θp=g(x)·θq,其中

    g(x)是值域在[0,1]間的單調(diào)遞減函數(shù),能確保算法的收斂性;x表示像素p和q的特征Vp和Vq間的相似度,可由歐氏距離度量得到;maxc為歸一化常數(shù),選取三維的灰度值作為特征矢量,maxc=441.673。相關的偽代碼可參見文獻[19]。

    2.2 決策級融合生成變化檢測圖

    不同的α代表選取不同的種子點,而不同的種子點產(chǎn)生不同的變化檢測結果,這表明變化檢測結果的效果與α的取值有關。為了使最終的變化檢測結果較穩(wěn)定,同時增強其對斑點噪聲的抗差性,對不同種子點下的變化檢測結果進行決策級融合。將每個α對應的變化檢測結果都作為一個專家的判斷,根據(jù)少數(shù)服從多數(shù)的原則,通過投票策略[20]決定最終的變化檢測結果。對位置(i,j)上的像素p,其類別為獲取投票數(shù)最多的一類。若以Vwk(i,j),k=c,n表示專家認為該像素屬于變化類(或非變化類)的票數(shù),則差異圖DA中每個像素的類別可按式(4)確定

    若給屬于變化類wc的像素賦予標簽1,而屬于非變化類wn的像素賦0,則對位置(i,j)上的像素p,其標簽Lp為最終,標簽圖L中,標簽值為1所在的位置即為檢測到的變化區(qū)域。

    3 試驗設計及試驗結果對比

    3.1 試驗數(shù)據(jù)描述

    第一組真實的試驗數(shù)據(jù)是關于Bern城市水災的SAR圖像,由歐洲遙感2號衛(wèi)星SAR傳感器獲得,大小為301×301像素,如圖2所示。圖2(a)、(b)分別為水災前1999年4月的圖像和水災后1999年5月的圖像,圖2(c)為變化參考圖,其中真實變化目標數(shù)為1155個像素單元。要檢測的是1999年5月水災引起的變化。

    圖2 第1組真實的有關Bern城市的試驗數(shù)據(jù)Fig.2 The first real dataset to be experimented on Bern city

    第2組真實的試驗數(shù)據(jù)是Ottawa地區(qū)水災的Radarsat SAR圖像,大小為290×350像素,如圖3所示。圖3(a)為正處雨季,1997年5月的圖像,圖3(b)為雨季過后,1997年8月的圖像,圖3(c)為變化參考圖,其中變化目標數(shù)為16 049個像素單元。兩幅圖像已經(jīng)過配準,要檢測的是因為雨季來臨,洪水泛濫造成的變化。

    圖3 第2組真實的有關Ottawa地區(qū)的試驗數(shù)據(jù)Fig.3 The second real dataset on Ottawa district

    3.2 試驗設計

    雖然已知α的值域為(0,1),且已知決策級融合可以消除不同的α對變化檢測性能的影響,但融合哪些α下的變化檢測結果得到最終的變化檢測結果,即存在α該如何取值的問題。為簡單起見,本文采取均分[0.05,0.95]的取值方式。這樣,最終變化檢測結果轉化為受取值間隔Δα的影響。為驗證本文方法的有效性,先分析參數(shù)Δα對變化檢測結果性能的影響,而后,從本文方法具有的3個特點組織了以下的比較試驗證明本文方法的有效性:① 與無監(jiān)督方式的變化檢測技術GGKIT[10]和LN-GKIT(或WR-GKIT)[11]比較;②與上下文敏感分析技術(EM+MRF)[13]比較;③與無分布假設及利用上下文信息的變化檢測技術(FFL-ARS)[5]比較。其中,本文方法的變化檢測結果為投票策略融合α={0.05:Δα:0.95},Δα=0.05下的變化檢測結果所得。

    為了更客觀地說明變化檢測結果的效果,給出4個定量評價檢測結果的分析指標:正確檢測數(shù)(right detections(RD),變化的像素被判斷為變化的像素個數(shù))、漏檢數(shù)(missed alarms(MA),變化的像素被判為未變化的像素個數(shù))、虛警(false alarms(FA),未變化的像素被判為變化的像素個數(shù),也叫錯誤檢測數(shù))和總錯誤檢測數(shù)(overall error(OE),即漏檢數(shù)+虛警)。一個好的變化檢測技術在于使正確檢測數(shù)盡可能地多,而總錯誤檢測數(shù)盡可能地少。

    3.3 參數(shù)分析

    為分析參數(shù)α的取值間隔Δα對變化檢測結果的影響,以Bern試驗數(shù)據(jù)為例,給出以不同取值間隔Δα={0.01,0.02,0.03,0.04,0.05,0.06,0.07,0.08,0.09,0.1}均分區(qū)間[0.05,0.95]下各α對應的變化檢測結果,并在圖4給出其定量分析指標隨參數(shù)Δα的變化曲線。

    圖4 本文方法的變化檢測結果受參數(shù)Δα的影響曲線Fig.4 Behaviors of quantitative analysis terms changing with the parameterΔα

    給出不同Δα下變化檢測結果的定量分析指標后,得出正確檢測數(shù)在90.4±1.562 2個像素單元內(nèi)波動,而虛警數(shù)在289.1±5.557 1個像素單元范圍內(nèi)波動,這說明本文方法具有較好的變化檢測精度及較穩(wěn)定的性能。圖4展示的是取值間隔Δα對變化檢測結果性能的影響曲線,也表明本方法的變化檢測性能受取值間隔Δα的影響較小。

    3.4 試驗結果及分析

    對Bern試驗數(shù)據(jù),各差異圖分析技術的變化檢測結果在圖5展示,其中圖5(a)、(b)分別為文獻[10—11]給出的對Bern試驗數(shù)據(jù)最好的變化檢測結果圖,圖5(e)為本文方法的變化檢測結果圖。表1給出了圖5中各變化檢測結果的定量分析指標。圖5(a)為對數(shù)比值法差異圖,假設其類分布模型符合廣義高斯分布(GGKIT)產(chǎn)生的變化檢測結果圖,圖5(b)為對比值法差異圖,假設其分布函數(shù)符合對數(shù)正態(tài)分布模型(LN-GKIT)得到的變化檢測結果;比較GGKIT和LN-GKIT法,GGKIT法的漏檢數(shù)比本文的多33個像素單元,LN-GKIT法的漏檢數(shù)比本文的多175個像素單元,這一方面是因為GGKIT和LN-GKIT法均需降斑預處理,而濾波程度設置得過小漏檢數(shù)較多,反之,則錯誤檢測數(shù)較多;另一方面是因為GGKIT和LN-GKIT法均涉及統(tǒng)計模型的選擇,對相同的試驗數(shù)據(jù),相同的差異圖構造,不同的分布模型對應不同的變化檢測結果,即變化檢測結果的好壞取決于所選與真實的分布模型間的擬合度,若選取的分布模型不合適的話,變化檢測效果會不如本文方法的。

    圖5 差異圖分析技術對Bern試驗數(shù)據(jù)均值比差異圖的變化檢測結果Fig.5 The results on Bern dataset of different techniques

    表1 圖5中各變化檢測結果的定量分析指標Tab.1 The quantitative analysis terms on results shown in Fig.5

    圖5(c)和圖5(d)分別為上下文敏感技術EM+MRF和FFL-ARS法的變化檢測結果。比較EM+MRF法與本文方法,雖然EM+MRF法利用了鄰域信息,但EM+MRF法多3085個像素單元。這主要是因為:① 差異圖中的斑點噪聲影響了EM+MRF法的變化檢測性能;②暗含地認為差異圖的類分布符合高斯模型,而這與Bern試驗數(shù)據(jù)的差異圖中兩類的真實分布不太符合。比較FFL-ARS法與本文方法,F(xiàn)FL-ARS法的總錯誤檢測數(shù)比本文方法多107個像素單元。雖然FFL-ARS法不涉及分布模型的選擇且利用了鄰域信息,但該方法丟棄了差異圖自身所含的信息,只利用了由靜態(tài)小波變換對差異圖進行分解后再丟棄高頻系數(shù)重構的圖像數(shù)據(jù)(即對差異圖的k(k≥1)層表示)信息,而無論k取值多少,細節(jié)信息難免會有所丟失。本文方法利用了差異圖自身的信息,在一定程度上抑制了細節(jié)信息的丟失。比較兩方法的漏檢數(shù),本文方法的漏檢數(shù)比FFL-ARS方法的少303個像素單元,這也定量地說明了FFL-ARS法細節(jié)信息丟失的較多。

    通過上述比較,結果顯示本文方法優(yōu)于上述變化檢測技術。本文方法的變化檢測結果擁有較多的正確檢測數(shù)的同時,含有的總錯誤檢測數(shù)也較少,這說明在總錯誤檢測數(shù)和正確檢測數(shù)這個兩難的問題上達到了一個較好的平衡。究其原因,主要是因為:① 本文方法不涉及降斑預處理,避免了因濾波程度設置不當引起效果降低的可能性,且設置特征時,利用了差異圖本身,在一定程度上抑制了信息的丟失;② 不涉及分布模型的選擇,變化檢測結果只與初始選取的“種子點”有關,而通過決策級融合后,變化檢測結果已相對穩(wěn)定,且進一步克服了斑點噪聲對變化檢測結果的影響;③ 在變化檢測過程利用鄰域內(nèi)像素的競爭完成變化檢測,將鄰域信息引入提高了變化檢測精度。圖5(f)為給定差異圖下能獲取的最優(yōu)的變化檢測結果,即MTEP技術的變化檢測結果,與之相比進一步說明本文方法效果的有效性。

    對Ottawa試驗數(shù)據(jù),文獻[10]中GGKIT法的變化檢測結果較好,即對數(shù)比值法構造的差異圖,假設其類條件分布符合廣義高斯分布更合適,如圖6(a)所示;對比值法構造的差異圖,假設其類分布符合“Weibull-Ratio”模型更合適,即WRGKIT法的變化檢測結果較好,如圖6(b)所示。圖6(c)為上下文敏感技術EM+MRF法的變化檢測結果圖,而圖6(d)為上下文敏感和無分布假設的FFL-ARS法的變化檢測圖;本文方法及MTEP技術的變化檢測結果分別在圖6(e)和圖6(f)展示。表2給出了各變化檢測結果的定量分析結果。

    圖6 差異圖分析技術對Ottawa試驗數(shù)據(jù)均值比差異圖的變化檢測結果Fig.6 The results on Ottawa dataset of different techniques

    表2 圖6中各變化檢測結果的定量分析指標Tab.2 The quantitative analysis terms on results shown in Fig.6

    分析表2的數(shù)據(jù)可知:與Bern試驗數(shù)據(jù)相同的是,對Ottawa試驗數(shù)據(jù),本文方法相比于其他技術,變化檢測效果有所改善。對Ottawa試驗數(shù)據(jù),本文方法變化檢測結果的總錯誤檢測數(shù)為1199個像素單元,比GGKIT法少1508個像素單元,比WR-GKIT法的少1394個像素單元;比上下文敏感的EM+MRF法少547個像素單元;比同時具有上下文敏感和無分布假設的FFL-ARS法少2002個像素單元。與MTEP技術的變化檢測結果相比,可進一步說明本文方法的有效性。而與Bern試驗數(shù)據(jù)不同的是,對Ottawa試驗數(shù)據(jù),文獻[11]中WR-GKIT法的變化檢測效果較好,而對Bern試驗數(shù)據(jù),相同方法構造的差異圖,LN-GKIT法的變化檢測效果較好,這進一步說明本文提出的變化檢測技術,不涉及分布模型假設的變化檢測技術的有效性。

    4 結 論

    本文利用差異圖的特點給出種子點后,結合種子點和一種基于鄰域像素競爭的交互式分割技術及決策級融合策略,給出了一種無監(jiān)督、上下文敏感及無分布假設SAR圖像變化檢測技術。該方法無需要求對SAR圖像作降斑預處理,且性能優(yōu)于其他相關的變化檢測技術,但當處理的圖像所含噪聲嚴重時,可以利用靜態(tài)小波變換對差異圖作多次分解后再重構,將特征設置為由灰度值構成的高維矢量,以增強對斑點噪聲的抗差性;也可以先作預處理去除部分斑點噪聲再由本文技術產(chǎn)生變化檢測結果,其中,濾波程度不夠強、仍有殘留的斑點噪聲對變化檢測結果有影響的情況可以不考慮,因為本文方法可以較好地處理含噪聲不嚴重的圖像。

    分割技術中鄰域像素間競爭實為特征間的相似性比較,本文以灰度為特征、歐式范數(shù)為測度。實際上,使用者可以考慮任何具有辨別力的特征,諸如紋理、方差、局部直方圖等以及相應的相似性度量測度,如自適應的距離測度[22]等。

    [1] RIDD M,LIU J.A Comparison of Four Algorithms for Change Detection in an Urban Environment[J].Remote Sensing of Environment,1998,63(2):95-100.

    [2] MAS J F.Monitoring Land-cover Changes:a Comparison of Change Detection Techniques[J].International Journal of Remote Sensing,1999,20(1):139-152.

    [3] WEYDAHL D J.Detecting Seasonal Changes at Nordic Latitudes Using ERS-1SAR Images[C]∥Proceedings of IEEE International Geoscience and Remote Sensing Symposium,IGARSS’93.Tokyo:IEEE Computer Society,1993:1462–1464.

    [4] COPPIN P,JONCKHEERE I,NACKAERTS K,et al.Digital Change Detection Methods in Ecosystem Monitoring:a Review[J].International Journal of Remote Sensing,2004,25(9):1565-1596.

    [5] BOVOLO F,BRUZZONE L.A Detail-preserving Scaledriven Approach to Change Detection in Multitemporal SAR Images[J].IEEE Transactions on Geoscience and Remote Sensing,2005,43(12):2963-2972.

    [6] RIGNOT E J M,VAN ZYL J J.Change Detection Techniques form ERS-1SAR Data[J].IEEE Transactions on Geoscience and Remote Sensing,1993,31(4):896-906.

    [7] ONANA V P,TROUVE E,MAURIS G,et al.Change Detection in Urban Context with Multitemporal ERS-SAR Images by Using Data Fusion Approach[C]∥Proceedings of IEEE International Geoscience and Remote Sensing Symposium,IGARSS’03.Toulouse:IEEE Computer Society,2003:3650-3652.

    [8] WU F,WANG C,Zhang H,et al.Change Detection and Analysis with Radarsat-1SAR Image[C]∥Proceedings of IEEE International Geoscience and Remote Sensing Symposium,IGARSS’07.Barcelona:IEEE Computer Society,2007:2601-2604.

    [9] DEKKER R J.Speckle Filtering in Satellite SAR Change Detection Imagery[J].International Journal of Remote Sensing,1998,19(6):1133-1146.

    [10] BAZI Y,BRUZZONE L,MELGANI F.An Unsupervised Approach Based on the Generalized Gaussian Model to Automatic Change Detection in Multitemporal SAR Images[J].IEEE Transactions on Geoscience and Remote Sensing,2005,43(4):874-887.

    [11] MOSER G,SERPICO S B.Generalized Minimum-error Thresholding for Unsupervised Change Detection from SAR Amplitude Imagery[J].IEEE Transactions on Geoscience and Remote Sensing,2006,44(10):2972-2983.

    [12] BRUZZONE L,PRIETO D F.Automatic Analysis of the Difference Image for Unsupervised Change Detection[J].IEEE Transactions on Geoscience and Remote Sensing,2000,38(3):1171-1182.

    [13] BAZI Y,MELGANI F,BRUZZONE L,et al.A Genetic Expectation-maximization Method for Unsupervised Change Detection in Multitemporal SAR Imagery[J].International Journal of Remote Sensing,2009,30(24):6591-6610.

    [14] YU H T,LIN A P,CHAI G L,et al.A DCT-based Change Detection Method for Multi-temporal SAR Images[C]∥Proceedings of IEEE International Geoscience and Remote Sensing Symposimu,IGARSS’08,Boston:IEEE Computer Society,2008:166-169.

    [15] ZHANG S Q,LU H Q.Learning Texture Classifier for Flooded Region Detection in SAR Images[C]∥Proceedings of the International Conference on Computer Graphics,Imaging and Visualization,CGIV’05.Washington DC:IEEE Computer Society,2005:93-98.

    [16] CELIK T,MA K K.Unsupervised Change Detection for Satellite Images Using Dual-tree Complex Wavelet Transform[J].IEEE Transactions on Geoscience and Remote Sensing,2010,48(3):1199-1210.

    [17] CELIK T.A Bayesian Approach to Unsupervised Multiscale Change Detection in Synthetic Aperture Radar Images[J].Signal Processing,2010,90(5):1147-1471.

    [18] HUANG Shiqi,LIU Daizhi,HU Mingxing,et al.Multi-temporal SAR Images Change Detection Technique Based on Wavelet Transform[J].Acta Geodaetica et Cartographica Sinica,2010,39(2):180-186.(黃世奇,劉代志,胡明星,等,基于小波變換的多時相SAR圖像變化檢測技術[J].測繪學報,2010,39(2):180-186.)

    [19] VEZHNEVETS V,KONOUCHINE V.“GrowCut”—Interactive Multi-label N-D Image Segmentation by Cellular Automata[C]∥Proceedings of Graphicon.Novosibirsk:[s.n.],2005.

    [20] KITTLER J,HATEF M,DUIN R P W,et al.On Combining Classifier[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,1998,20(3):226-239.

    [21] KITTLER J,ILLINGWORTH J.Minimum Error Thresholding[J].Journal of Pattern Recognition,1986,19(1):41-47.

    [22] DAVIS J V,KULIS B,JAIN P,et al.Information—Theoretic Metric Learning[C]∥Proceedings of the 24th International Conference on Machine learning,ICML’07.New York:ACM,2007:209-216.

    Interactive Segmentation Technique and Decision-level Fusion Based Change Detection for SAR Images

    WAN Honglin,JIAO Licheng,XIN Fangfang
    The Key Laboratory of Intelligent Perception and Image Understanding of Ministry of Education,Xidian University,Xi’an 710071,China

    To avoid reducing the speckle noise,and meanwhile overcome the limitation of selecting the distribution model,first,characteristics of the difference image(DI)are integrated with an interactive segmentation method not referring to any distribution assumption,to generate change detection maps corresponding to different“seeds”,then a voting competition strategy is used to fuse those results to give the final change detection map.During segmenting,the feature of each pixel is set as a vector consisted of the corresponding intensities in the DI and each scale representation of the DI given by stationary wavelet transform(SWT).This kind of features and the decision level fusion make our proposed method robust to the speckle noise.Results on real SAR datasets obtained under the situation that there is no despeckling preprocessing of SAR images confirm the effectiveness of our method.

    interactive segmentation method;characteristics of the difference image;decision-level fusion;stationary wavelet transform;SAR images;change detection

    WAN Honglin(1983—),female,PhD candidate,majors in SAR image segmentation,object detection and change detection.

    WAN Honglin,JIAO Licheng,XIN Fangfang.Interactive Segmentation Technique and Decision-level Fusion Based Change Detection for SAR Images[J].Acta Geodaetica et Cartographica Sinica,2012,41(1):74-80.(萬紅林,焦李成,辛芳芳.基于交互式分割技術和決策級融合的SAR圖像變化檢測[J].測繪學報,2012,41(1):74-80.)

    TP751

    A

    1001-1595(2012)01-0074-07

    國家863計劃(2008AA01Z125;2009AA12Z210);國家自然科學基金(60702062;60803097;60972148;60971128;60970066);國家部委科技項目(XADZ2008159;51307040103);教育部博士點基金(200807010003);教育部重點項目(108115);高等學校學科創(chuàng)新引智計劃(111計劃)(B07048);中央高?;究蒲袠I(yè)務費專項(JY10000902001;JY10000902038;JY10000902032)

    雷秀麗)

    2010-07-15

    2011-03-15

    萬紅林(1983—),女,博士生,主要研究方向為SAR圖像分割、目標檢測及變化檢測等。

    E-mail:hlwan@m(xù)ail.xidian.edu.cn

    猜你喜歡
    變化檢測鄰域像素
    趙運哲作品
    藝術家(2023年8期)2023-11-02 02:05:28
    用于遙感圖像變化檢測的全尺度特征聚合網(wǎng)絡
    像素前線之“幻影”2000
    基于多尺度紋理特征的SAR影像變化檢測
    稀疏圖平方圖的染色數(shù)上界
    “像素”仙人掌
    基于稀疏表示的視網(wǎng)膜圖像對變化檢測
    基于鄰域競賽的多目標優(yōu)化算法
    自動化學報(2018年7期)2018-08-20 02:59:04
    基于Landsat影像的黃豐橋林場森林變化檢測研究
    關于-型鄰域空間
    国产伦人伦偷精品视频| 亚洲av成人不卡在线观看播放网| 18禁裸乳无遮挡免费网站照片 | 亚洲五月婷婷丁香| 亚洲av熟女| 在线天堂中文资源库| 亚洲国产毛片av蜜桃av| 成人特级黄色片久久久久久久| 丝袜美足系列| 午夜成年电影在线免费观看| 给我免费播放毛片高清在线观看| 91麻豆精品激情在线观看国产| 亚洲成人久久性| 成人国产综合亚洲| 搡老熟女国产l中国老女人| 少妇的丰满在线观看| 一区二区三区精品91| 久久久久久大精品| 黄色 视频免费看| 丰满人妻熟妇乱又伦精品不卡| 日本a在线网址| 久久精品国产清高在天天线| 国产精品免费视频内射| 成人特级黄色片久久久久久久| 欧美+亚洲+日韩+国产| 久久欧美精品欧美久久欧美| 侵犯人妻中文字幕一二三四区| 欧美日本视频| 久久精品国产综合久久久| 国产精品野战在线观看| 人人澡人人妻人| 欧美一级a爱片免费观看看 | 自线自在国产av| 黑人操中国人逼视频| 午夜福利成人在线免费观看| 老鸭窝网址在线观看| 亚洲精品美女久久久久99蜜臀| 91字幕亚洲| 妹子高潮喷水视频| 日韩精品中文字幕看吧| 淫妇啪啪啪对白视频| 欧美一级毛片孕妇| 一区在线观看完整版| 国产成人影院久久av| 大香蕉久久成人网| 国产精品久久久久久人妻精品电影| 久久人人97超碰香蕉20202| 国产亚洲精品第一综合不卡| 亚洲 国产 在线| ponron亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 久久亚洲精品不卡| 日本精品一区二区三区蜜桃| 国产精品久久视频播放| 国产精品亚洲美女久久久| 老司机午夜福利在线观看视频| 国产成人影院久久av| 99在线视频只有这里精品首页| 午夜精品在线福利| 禁无遮挡网站| 后天国语完整版免费观看| 国产欧美日韩精品亚洲av| 日韩三级视频一区二区三区| 午夜福利免费观看在线| 俄罗斯特黄特色一大片| av有码第一页| 久久久久久久午夜电影| 黄色视频不卡| 精品久久久久久,| 曰老女人黄片| 中文字幕久久专区| 亚洲国产中文字幕在线视频| 欧美日韩福利视频一区二区| 久久精品成人免费网站| 大陆偷拍与自拍| 久久天堂一区二区三区四区| 欧美日本亚洲视频在线播放| 亚洲成人免费电影在线观看| 国产免费男女视频| 黄片播放在线免费| 国产亚洲欧美98| 十八禁人妻一区二区| 麻豆av在线久日| 视频区欧美日本亚洲| 久久天堂一区二区三区四区| 国产成人系列免费观看| 国产av在哪里看| 亚洲精品一区av在线观看| 亚洲 国产 在线| 色在线成人网| 国产亚洲精品第一综合不卡| 欧美亚洲日本最大视频资源| 国产精品99久久99久久久不卡| 亚洲一区二区三区不卡视频| 久久国产精品人妻蜜桃| 操美女的视频在线观看| av网站免费在线观看视频| 亚洲精品国产精品久久久不卡| 欧美激情 高清一区二区三区| 88av欧美| 极品教师在线免费播放| av超薄肉色丝袜交足视频| 69精品国产乱码久久久| 三级毛片av免费| 中文字幕色久视频| 欧美久久黑人一区二区| 亚洲av片天天在线观看| 美女国产高潮福利片在线看| 亚洲一区中文字幕在线| 又黄又爽又免费观看的视频| 一级黄色大片毛片| 亚洲中文字幕一区二区三区有码在线看 | 岛国在线观看网站| 免费一级毛片在线播放高清视频 | 成人国产一区最新在线观看| 黄网站色视频无遮挡免费观看| 又黄又爽又免费观看的视频| 国产激情欧美一区二区| 正在播放国产对白刺激| 少妇的丰满在线观看| 一进一出抽搐动态| 亚洲五月婷婷丁香| 国产精品 国内视频| 欧美在线一区亚洲| 日韩欧美免费精品| 亚洲九九香蕉| 日韩三级视频一区二区三区| 亚洲精品久久成人aⅴ小说| 不卡一级毛片| 不卡一级毛片| 99国产精品一区二区蜜桃av| 天堂动漫精品| 久久国产精品人妻蜜桃| 18美女黄网站色大片免费观看| 国内毛片毛片毛片毛片毛片| 在线天堂中文资源库| 国产精品二区激情视频| 午夜福利欧美成人| 97碰自拍视频| 婷婷精品国产亚洲av在线| 日韩中文字幕欧美一区二区| 国产乱人伦免费视频| 欧美乱色亚洲激情| 久久久久久久久中文| 日韩 欧美 亚洲 中文字幕| 丁香欧美五月| 久久亚洲精品不卡| 国产欧美日韩精品亚洲av| or卡值多少钱| 午夜亚洲福利在线播放| 久久久久国内视频| 日韩有码中文字幕| 黄色片一级片一级黄色片| 精品无人区乱码1区二区| 如日韩欧美国产精品一区二区三区| 国产成人一区二区三区免费视频网站| 亚洲人成电影免费在线| 亚洲一卡2卡3卡4卡5卡精品中文| 久久香蕉激情| 国产成人一区二区三区免费视频网站| 欧美黄色淫秽网站| 最新在线观看一区二区三区| 中文字幕色久视频| 十分钟在线观看高清视频www| 精品久久蜜臀av无| 国产精品 欧美亚洲| 色尼玛亚洲综合影院| 日韩大尺度精品在线看网址 | 每晚都被弄得嗷嗷叫到高潮| 亚洲伊人色综图| av天堂在线播放| 亚洲国产精品999在线| 日韩欧美一区二区三区在线观看| 天堂动漫精品| 久久久久久国产a免费观看| 欧美日本中文国产一区发布| 免费高清视频大片| 国产一区二区三区在线臀色熟女| 国产av精品麻豆| 99riav亚洲国产免费| 成人国语在线视频| 怎么达到女性高潮| 97人妻天天添夜夜摸| www.自偷自拍.com| 在线永久观看黄色视频| 视频区欧美日本亚洲| www.www免费av| www国产在线视频色| 欧美久久黑人一区二区| 成人永久免费在线观看视频| 12—13女人毛片做爰片一| 午夜福利一区二区在线看| 一a级毛片在线观看| 一级a爱片免费观看的视频| 国产xxxxx性猛交| 国产不卡一卡二| 日韩 欧美 亚洲 中文字幕| 国产高清激情床上av| 欧美黑人精品巨大| 精品人妻在线不人妻| 久久婷婷人人爽人人干人人爱 | 欧美成人一区二区免费高清观看 | 女性生殖器流出的白浆| 国产高清有码在线观看视频 | 亚洲人成电影免费在线| 欧美一级毛片孕妇| 99精品在免费线老司机午夜| 成人欧美大片| 在线观看舔阴道视频| 日韩精品中文字幕看吧| 亚洲成人免费电影在线观看| 国产精品久久久久久亚洲av鲁大| 日韩欧美国产一区二区入口| 欧美av亚洲av综合av国产av| 黄色女人牲交| 国产精品,欧美在线| 成人免费观看视频高清| 电影成人av| 亚洲成a人片在线一区二区| 免费在线观看日本一区| 热re99久久国产66热| 日本撒尿小便嘘嘘汇集6| 精品久久久久久久人妻蜜臀av | 嫩草影视91久久| 国产不卡一卡二| 日日夜夜操网爽| 国产亚洲精品av在线| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲午夜精品一区,二区,三区| 久久欧美精品欧美久久欧美| 国产精品一区二区三区四区久久 | av欧美777| 久久人妻福利社区极品人妻图片| 香蕉国产在线看| 成人亚洲精品av一区二区| 啦啦啦观看免费观看视频高清 | 91成年电影在线观看| 日本a在线网址| 丝袜人妻中文字幕| 精品欧美国产一区二区三| x7x7x7水蜜桃| 精品一区二区三区视频在线观看免费| 动漫黄色视频在线观看| 国产午夜福利久久久久久| www.自偷自拍.com| 99re在线观看精品视频| 丰满的人妻完整版| 欧美日本中文国产一区发布| 极品教师在线免费播放| 母亲3免费完整高清在线观看| 99在线视频只有这里精品首页| 夜夜夜夜夜久久久久| 日韩欧美在线二视频| 国语自产精品视频在线第100页| 国产精品久久久av美女十八| 老司机在亚洲福利影院| 国产成人欧美| 999久久久精品免费观看国产| 国产精品永久免费网站| 人人澡人人妻人| 久久精品亚洲熟妇少妇任你| 亚洲欧美日韩另类电影网站| 欧美一区二区精品小视频在线| 日韩成人在线观看一区二区三区| 国产主播在线观看一区二区| 午夜福利,免费看| 19禁男女啪啪无遮挡网站| 日本免费a在线| 男人的好看免费观看在线视频 | 国产1区2区3区精品| 欧美 亚洲 国产 日韩一| 日韩中文字幕欧美一区二区| 久久人人精品亚洲av| 日韩欧美三级三区| 亚洲黑人精品在线| 国产熟女午夜一区二区三区| 国产高清videossex| 看黄色毛片网站| 999精品在线视频| 日韩精品免费视频一区二区三区| or卡值多少钱| 久久午夜亚洲精品久久| 法律面前人人平等表现在哪些方面| 91精品国产国语对白视频| 亚洲色图av天堂| 丁香欧美五月| videosex国产| 一级片免费观看大全| 国产区一区二久久| 人妻丰满熟妇av一区二区三区| 国产精品免费视频内射| 人人妻人人澡欧美一区二区 | 精品日产1卡2卡| 亚洲精品av麻豆狂野| 亚洲精品一卡2卡三卡4卡5卡| 搡老岳熟女国产| 日本精品一区二区三区蜜桃| 少妇裸体淫交视频免费看高清 | 免费看十八禁软件| 欧美成人免费av一区二区三区| 最近最新中文字幕大全电影3 | 国产高清有码在线观看视频 | 欧美人与性动交α欧美精品济南到| 日本vs欧美在线观看视频| 久久青草综合色| 亚洲欧美激情在线| 欧美激情极品国产一区二区三区| 欧美日韩一级在线毛片| 精品第一国产精品| 大码成人一级视频| 亚洲精品国产一区二区精华液| 757午夜福利合集在线观看| 日本 av在线| 日韩欧美一区二区三区在线观看| 久久国产乱子伦精品免费另类| 一个人免费在线观看的高清视频| 性少妇av在线| 色播在线永久视频| 精品人妻1区二区| 中文字幕另类日韩欧美亚洲嫩草| 亚洲免费av在线视频| 国产精品亚洲美女久久久| 国产精品久久久av美女十八| 最近最新中文字幕大全免费视频| 校园春色视频在线观看| 午夜福利影视在线免费观看| 无遮挡黄片免费观看| а√天堂www在线а√下载| 午夜成年电影在线免费观看| 日韩中文字幕欧美一区二区| 长腿黑丝高跟| 亚洲人成伊人成综合网2020| 亚洲天堂国产精品一区在线| 欧美国产精品va在线观看不卡| ponron亚洲| 99热只有精品国产| 亚洲av熟女| 最新在线观看一区二区三区| 高潮久久久久久久久久久不卡| 国产精品免费视频内射| 村上凉子中文字幕在线| 国产精品自产拍在线观看55亚洲| 十八禁人妻一区二区| 色尼玛亚洲综合影院| 亚洲最大成人中文| 久久久国产欧美日韩av| 亚洲av成人av| 国产精华一区二区三区| www.精华液| 亚洲精品中文字幕在线视频| ponron亚洲| 午夜福利,免费看| 久久精品成人免费网站| 免费久久久久久久精品成人欧美视频| netflix在线观看网站| 无遮挡黄片免费观看| 一个人观看的视频www高清免费观看 | 亚洲国产精品久久男人天堂| 精品高清国产在线一区| 老汉色∧v一级毛片| 88av欧美| 久久午夜亚洲精品久久| 窝窝影院91人妻| 99国产极品粉嫩在线观看| 亚洲av电影在线进入| 女人精品久久久久毛片| 国产成+人综合+亚洲专区| 色在线成人网| 18禁美女被吸乳视频| 亚洲免费av在线视频| 9191精品国产免费久久| 大型黄色视频在线免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品美女久久久久99蜜臀| 成人手机av| 日韩大尺度精品在线看网址 | 免费在线观看黄色视频的| 国产成人免费无遮挡视频| 色在线成人网| 十八禁人妻一区二区| 国产1区2区3区精品| 搡老岳熟女国产| 日韩大码丰满熟妇| 日日摸夜夜添夜夜添小说| 村上凉子中文字幕在线| 国产精品久久视频播放| 看片在线看免费视频| 无遮挡黄片免费观看| 成人手机av| 免费在线观看视频国产中文字幕亚洲| 精品一品国产午夜福利视频| 淫秽高清视频在线观看| 黄色片一级片一级黄色片| 日日爽夜夜爽网站| 老司机靠b影院| av有码第一页| 日本一区二区免费在线视频| 久久国产亚洲av麻豆专区| www.自偷自拍.com| 亚洲性夜色夜夜综合| 亚洲一区高清亚洲精品| 又黄又粗又硬又大视频| 亚洲自拍偷在线| 极品教师在线免费播放| 在线av久久热| 天天一区二区日本电影三级 | 精品熟女少妇八av免费久了| 久久 成人 亚洲| 涩涩av久久男人的天堂| 国产高清有码在线观看视频 | 在线十欧美十亚洲十日本专区| 国产av一区二区精品久久| 久久久久国内视频| 人妻丰满熟妇av一区二区三区| 在线观看66精品国产| 丰满的人妻完整版| 一区二区三区精品91| netflix在线观看网站| av在线播放免费不卡| 日韩大码丰满熟妇| 久久午夜综合久久蜜桃| 精品少妇一区二区三区视频日本电影| 亚洲精品粉嫩美女一区| 可以免费在线观看a视频的电影网站| 国产伦一二天堂av在线观看| 男人舔女人下体高潮全视频| 麻豆国产av国片精品| 一本综合久久免费| 亚洲,欧美精品.| 最新在线观看一区二区三区| 久久香蕉精品热| 丝袜美腿诱惑在线| 我的亚洲天堂| 国产又爽黄色视频| 欧美日韩福利视频一区二区| 欧美激情久久久久久爽电影 | 国内久久婷婷六月综合欲色啪| 免费在线观看影片大全网站| av视频在线观看入口| 免费观看精品视频网站| 97人妻精品一区二区三区麻豆 | 久久久久亚洲av毛片大全| 精品一品国产午夜福利视频| 高清黄色对白视频在线免费看| 两人在一起打扑克的视频| 窝窝影院91人妻| 亚洲成人精品中文字幕电影| 色播亚洲综合网| 国产在线精品亚洲第一网站| 大型黄色视频在线免费观看| 少妇裸体淫交视频免费看高清 | 91九色精品人成在线观看| 国产蜜桃级精品一区二区三区| 国产一区二区激情短视频| 国产91精品成人一区二区三区| 欧美色欧美亚洲另类二区 | 一个人免费在线观看的高清视频| 国产成人精品久久二区二区91| 最新在线观看一区二区三区| 国产精品99久久99久久久不卡| 国产一区二区三区综合在线观看| x7x7x7水蜜桃| 欧美激情久久久久久爽电影 | 欧美人与性动交α欧美精品济南到| 男男h啪啪无遮挡| 巨乳人妻的诱惑在线观看| 亚洲五月婷婷丁香| 一级毛片精品| 日韩精品免费视频一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 真人做人爱边吃奶动态| 波多野结衣高清无吗| 这个男人来自地球电影免费观看| 精品第一国产精品| 亚洲av成人不卡在线观看播放网| 国产一区二区三区在线臀色熟女| 色尼玛亚洲综合影院| 日韩成人在线观看一区二区三区| 精品欧美国产一区二区三| 成年版毛片免费区| 亚洲九九香蕉| 久久青草综合色| 少妇被粗大的猛进出69影院| 女人被狂操c到高潮| 一区二区三区高清视频在线| 免费看十八禁软件| 少妇熟女aⅴ在线视频| 91成年电影在线观看| 黄片播放在线免费| 午夜福利一区二区在线看| 欧美黄色淫秽网站| 成人国产一区最新在线观看| 久久九九热精品免费| 曰老女人黄片| 免费看美女性在线毛片视频| 热re99久久国产66热| 亚洲精品美女久久av网站| 精品一区二区三区av网在线观看| 一区福利在线观看| 极品人妻少妇av视频| 男人舔女人的私密视频| 一区在线观看完整版| 亚洲电影在线观看av| 制服丝袜大香蕉在线| 精品一区二区三区av网在线观看| 亚洲国产毛片av蜜桃av| 如日韩欧美国产精品一区二区三区| 欧美性长视频在线观看| 国产欧美日韩精品亚洲av| 国产av一区在线观看免费| 人人妻人人澡欧美一区二区 | 国产成人精品久久二区二区免费| 18禁国产床啪视频网站| 亚洲一区高清亚洲精品| 老汉色∧v一级毛片| 国产成人精品久久二区二区91| 免费一级毛片在线播放高清视频 | 亚洲欧美激情综合另类| 亚洲电影在线观看av| 免费女性裸体啪啪无遮挡网站| 国产99久久九九免费精品| 久久久国产欧美日韩av| 亚洲人成电影免费在线| 亚洲专区国产一区二区| 久久人妻福利社区极品人妻图片| 欧美一级a爱片免费观看看 | 中出人妻视频一区二区| 久久久国产欧美日韩av| 成人国语在线视频| 精品久久久久久,| 国产97色在线日韩免费| 成人国产一区最新在线观看| 丝袜人妻中文字幕| 欧美绝顶高潮抽搐喷水| 精品乱码久久久久久99久播| 99久久综合精品五月天人人| 欧美+亚洲+日韩+国产| 制服诱惑二区| av免费在线观看网站| 国产伦人伦偷精品视频| 国产精品野战在线观看| 久久天堂一区二区三区四区| 大码成人一级视频| 国产精品影院久久| 亚洲色图 男人天堂 中文字幕| 老司机在亚洲福利影院| 香蕉丝袜av| 成人精品一区二区免费| 欧美日本中文国产一区发布| av天堂久久9| 久久狼人影院| 夜夜夜夜夜久久久久| 好男人电影高清在线观看| 欧美日韩福利视频一区二区| 一区二区三区激情视频| 一级片免费观看大全| 超碰成人久久| 国内毛片毛片毛片毛片毛片| 精品高清国产在线一区| 亚洲精品美女久久av网站| 精品第一国产精品| 国产av一区二区精品久久| 老司机午夜福利在线观看视频| 好男人电影高清在线观看| 久久九九热精品免费| 国产精品免费视频内射| а√天堂www在线а√下载| 国产亚洲精品综合一区在线观看 | 成年女人毛片免费观看观看9| 伊人久久大香线蕉亚洲五| 国产熟女午夜一区二区三区| 中文亚洲av片在线观看爽| 一个人观看的视频www高清免费观看 | 精品国产一区二区三区四区第35| 色av中文字幕| 巨乳人妻的诱惑在线观看| 成人国语在线视频| 好男人电影高清在线观看| 女警被强在线播放| 纯流量卡能插随身wifi吗| 99精品久久久久人妻精品| 国产一级毛片七仙女欲春2 | a在线观看视频网站| 少妇 在线观看| 麻豆av在线久日| 在线观看免费视频日本深夜| 久久天躁狠狠躁夜夜2o2o| 97人妻天天添夜夜摸| 中亚洲国语对白在线视频| 午夜两性在线视频| 十八禁人妻一区二区| 成人特级黄色片久久久久久久| 亚洲五月婷婷丁香| 亚洲精品国产区一区二| 亚洲人成电影观看| 亚洲人成网站在线播放欧美日韩| 久久精品91无色码中文字幕| 国产精品1区2区在线观看.| 免费看十八禁软件| 国产精品亚洲一级av第二区| 一级片免费观看大全| 亚洲午夜精品一区,二区,三区| 久久久久久久精品吃奶| 色老头精品视频在线观看| 搡老妇女老女人老熟妇| 亚洲专区字幕在线| 亚洲av熟女| 国产精品1区2区在线观看.| 中国美女看黄片| 一级黄色大片毛片| 麻豆av在线久日| 黄网站色视频无遮挡免费观看| 亚洲色图综合在线观看| 国产欧美日韩一区二区三| 国产aⅴ精品一区二区三区波|