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

    SAR圖像海岸線檢測的區(qū)域距離正則化幾何主動輪廓模型

    2016-11-07 02:37:28姜大偉范劍超黃鳳榮
    測繪學(xué)報(bào) 2016年9期
    關(guān)鍵詞:海岸線輪廓邊緣

    姜大偉,范劍超,黃鳳榮

    1. 遼寧師范大學(xué)城市與環(huán)境學(xué)院,遼寧 大連116029; 2. 國家海洋環(huán)境監(jiān)測中心海域資源環(huán)境動態(tài)監(jiān)管中心,遼寧 大連 116023

    ?

    SAR圖像海岸線檢測的區(qū)域距離正則化幾何主動輪廓模型

    姜大偉1,2,范劍超2,黃鳳榮1

    1. 遼寧師范大學(xué)城市與環(huán)境學(xué)院,遼寧 大連116029; 2. 國家海洋環(huán)境監(jiān)測中心海域資源環(huán)境動態(tài)監(jiān)管中心,遼寧 大連 116023

    合成孔徑雷達(dá)(SAR)衛(wèi)星遙感圖像可以極大地提高全國海岸線覆蓋頻率,然而受到海洋波浪所引起的隨機(jī)海水表面粗糙度的影響,海岸目標(biāo)與海水背景邊界易混淆不清,因此本文提出了基于區(qū)域距離正則化幾何主動輪廓模型(RDRGAC),引入距離正則項(xiàng),解決重復(fù)初始化水平集函數(shù)為符號距離函數(shù)的問題,提高了算法收斂速度。此外,將區(qū)域面積項(xiàng)系數(shù)與SAR圖像等效視數(shù)(ENL)建立非線性擬合關(guān)系,實(shí)現(xiàn)RDRGAC模型根據(jù)不同SAR遙感圖像的自適應(yīng)調(diào)整,改善海岸線自動提取精度。通過河北省北戴河和大連市金州灣SAR數(shù)據(jù)海岸線提取對比試驗(yàn),驗(yàn)證了所提方法的有效性。

    相干斑噪聲;合成孔徑雷達(dá);非線性擬合;幾何主動輪廓模型;等效視數(shù)

    海岸線是劃分海洋與陸地管理區(qū)域的基準(zhǔn)線,其所處地帶憑借豐富的自然資源和優(yōu)越的地理位置成為人類開發(fā)利用的重要區(qū)域[1]。隨著海洋經(jīng)濟(jì)的迅速發(fā)展和用地需求的快速增長,圍海造地成為緩解沿海地區(qū)土地供需矛盾、擴(kuò)大發(fā)展空間的有效途徑,致使海岸線不斷收縮或擴(kuò)張。國家海域本體庫建設(shè)需要對全國海岸線及海島開發(fā)利用情況進(jìn)行精確測繪調(diào)查,因此高效準(zhǔn)確的海岸線測繪,可以實(shí)現(xiàn)對海岸帶資源的持續(xù)利用和科學(xué)管理,具有十分重要的意義。

    傳統(tǒng)海岸線檢測的方法主要是利用歷史海圖、地形圖、剖面監(jiān)測等[2],這些方法費(fèi)時費(fèi)力,并且動態(tài)連續(xù)性差。遙感技術(shù)憑借自身宏觀、快速、動態(tài)連續(xù)以及低成本等優(yōu)點(diǎn)在海岸線檢測的工作中受到廣泛應(yīng)用。但是可見光遙感受光照條件限制,例如在沿海多陰雨多云天氣,大多數(shù)情況下,可見光遙感圖像無法每年完全覆蓋岸線一次,因此使得應(yīng)用微波波段的合成孔徑雷達(dá)(synthetic aperture radar, SAR)遙感圖像受到廣泛關(guān)注,并逐漸成為海域資源測繪研究的熱點(diǎn)。SAR是主動式成像傳感器,其具有全天時、全天候和高分辨率的數(shù)據(jù)獲取能力[3],在資源探測、地理測繪和海域監(jiān)測等領(lǐng)域有廣泛的應(yīng)用。然而SAR圖像中的相干斑噪聲使圖像的信噪比下降,降低了圖像的有效分辨率,從而給SAR圖像的分割和解譯帶來困難。

    傳統(tǒng)用于遙感圖像海岸線提取的方法有閾值分割[4]、區(qū)域生長[5]等圖像分割法和Sobel算子、Laplace算子、Canny算子、Susan算子等經(jīng)典邊緣檢測算子法[6-8]。這些方法易受噪聲影響,而且檢測出的邊緣不連續(xù)、難以提取到海岸線真實(shí)輪廓。文獻(xiàn)[9]提出邊界追蹤算法,利用圖像的梯度特征,通過閾值分割,實(shí)現(xiàn)了海岸線輪廓提取。這是一種較早提出的SAR圖像海岸線輪廓提取的方法。隨后,文獻(xiàn)[10]提出Markovian分割法。文獻(xiàn)[11]提出Snake(主動輪廓)模型方法,在近目標(biāo)處設(shè)定一個初始輪廓,并將輪廓曲線的運(yùn)動看作力合成的結(jié)果,從而設(shè)計(jì)曲線的能量函數(shù),最后通過能量函數(shù)的最小化使得初始輪廓逼近真實(shí)目標(biāo)邊界。文獻(xiàn)[12]將水平集方法成功引入主動輪廓模型,將二維平面內(nèi)的輪廓曲線隱含地表達(dá)為更高維曲面的一個水平集,提出了幾何活動輪廓(geometric active contour, GAC)模型。平均密度法是文獻(xiàn)[13]在2001年根據(jù)水平集方法和Mumford-Shah函數(shù)法的思想提出的。這個方法的原理是給出初始輪廓線,通過曲線的運(yùn)動不斷逼近目標(biāo)的輪廓邊界,但是與主動輪廓法不同的是,它不是通過檢測圖像梯度而是通過計(jì)算曲線內(nèi)外的平均灰度來驅(qū)使曲線運(yùn)動。文獻(xiàn)[14]通過引入?yún)^(qū)域信息對Snake模型的驅(qū)動力修正,利用粒子群優(yōu)化算法的全局優(yōu)化特性和數(shù)值穩(wěn)定性對模型分割結(jié)果優(yōu)化。在傳統(tǒng)Chan-Vese(C-V)模型的基礎(chǔ)上,文獻(xiàn)[15]充分考慮SAR圖像的斑點(diǎn)噪聲,引入距離正則項(xiàng),避免了水平集函數(shù)的周期性更新,解決了傳統(tǒng)水平集方法計(jì)算量大、分割速度慢的問題。文獻(xiàn)[16]提出基于雙重輪廓演化曲線的圖形分割水平集模型,通過內(nèi)部能量泛函,避免對符號距離函數(shù)的重新初始化,以及對初始輪廓線的過度依賴。文獻(xiàn)[17]提出一種基于G0分布的局部統(tǒng)計(jì)活動輪廓模型對SAR圖像海岸線檢測,并加入水平集懲罰項(xiàng),消除了重新初始化過程,實(shí)現(xiàn)了精確海岸線檢測,但缺點(diǎn)在于模型分割產(chǎn)生的極小值點(diǎn)問題會受到初始輪廓線位置的影響。

    針對SAR圖像包含大量不規(guī)則相干斑噪聲,及受到海洋波浪引起的隨機(jī)海水表面粗糙度干擾等原因,本文提出基于區(qū)域的距離正則化幾何主動輪廓模型(regional distance regularized geometric active contour models, RDRGAC),克服周期性重復(fù)初始化水平集函數(shù)為符號距離函數(shù)的缺點(diǎn),并將區(qū)域面積項(xiàng)系數(shù)與等效視數(shù)(equivalent number of looks, ENL)建立非線性擬合關(guān)系,根據(jù)每幅SAR遙感圖像中所包含相干斑噪聲的水平實(shí)現(xiàn)參數(shù)自動設(shè)定,提高海岸線自動提取精度和計(jì)算速度。

    1 研究區(qū)域SAR數(shù)據(jù)特征和試驗(yàn)環(huán)境分析

    1.1研究區(qū)域概況

    本文主要采用北戴河和金州灣地區(qū)進(jìn)行仿真試驗(yàn)。北戴河位于河北省東北部,西距首都北京279 km,東距秦皇島港18 km、山海關(guān)機(jī)場25 km。東北與秦皇島市海港區(qū)毗鄰,西部、西北部與撫寧縣接壤,東、南面臨渤海灣。東西長11.2 km,南北寬10.15 km,總面積70.14 km2。

    金州灣位于大連市金州鎮(zhèn)西2 km渤海海域,海域?qū)掗煶蕶E圓形,南北長28 km,寬約15 km。岸線北自金州區(qū)大魏家鎮(zhèn)蕎麥山、葫蘆套一帶,南至甘井子區(qū)黃龍尾咀,長約74 km,面積約7453 km2,是大連市地區(qū)渤海沿岸最大的海灣。目前在建的大連市金州灣國際機(jī)場,采取離岸填海建造人工島方式建設(shè),因此實(shí)現(xiàn)該區(qū)域海岸線檢測具有十分重要的意義。

    1.2SAR數(shù)據(jù)分析

    河北省北戴河SAR遙感圖像采用2015年5月20日RADARSAT-2衛(wèi)星數(shù)據(jù),大連市金州灣試驗(yàn)數(shù)據(jù)采用2015年11月28日RADARSAT-2衛(wèi)星遙感圖像,C波段,空間分辨率為8 m,幅寬為50 km×50 km。通過全極化成像模式對海岸線進(jìn)行分析,選用HV交叉極化模式,更有利于區(qū)分獲得海岸線。數(shù)據(jù)級別為SSG地理編碼系統(tǒng)校正產(chǎn)品,該產(chǎn)品在SGF的基礎(chǔ)上進(jìn)行了地圖投影(UTM)校正。

    1.3試驗(yàn)分析環(huán)境

    為了驗(yàn)證算法的分割性能,采用真實(shí)SAR圖像數(shù)據(jù)進(jìn)行試驗(yàn)分析。試驗(yàn)環(huán)境:CPU為Inter(R) Core i3-4150,主頻為3.50 GHz,內(nèi)存為4 GB,軟件環(huán)境為Matlab(R2013b)。SAR圖像數(shù)據(jù)預(yù)處理過程包括裁剪、鑲嵌和配準(zhǔn)均采用ERDAS IMAGINE 9.2版本來實(shí)現(xiàn)。

    2 算法提出

    2.1經(jīng)典幾何主動輪廓模型C-V模型

    Chan-Vese根據(jù)Mumford-Shah(M-S)模型提出了簡化M-S模型C-V模型[18]。假設(shè)定義域?yàn)棣傅膱D像I(x,y)被輪廓線C劃分為內(nèi)部平均灰度c1和外部平均灰度c2,曲線C將圖像分成了外部Ωext和內(nèi)部Ωint兩部分,則目標(biāo)最小化能量泛函可表示為

    F(C,c1,c2)=μL(C)+vA(C)+λ1∫inside(c)

    (1)

    式中,L(C)為閉合輪廓線C的長度項(xiàng);A(C)為閉合輪廓線C面積項(xiàng);μ、v、λ1、λ2是能量項(xiàng)系數(shù);顯然,當(dāng)閉合輪廓線運(yùn)動到圖像邊緣時,能量函數(shù)達(dá)到最小值。

    引入水平集函數(shù)得到以水平集函數(shù)φ表達(dá)的偏微分方程如下

    (2)

    式中

    (3)

    (4)

    2.2RDRGAC模型

    經(jīng)典水平集方法中,水平集在演化過程中會出現(xiàn)非正則性,這會導(dǎo)致數(shù)值錯誤并影響演化的穩(wěn)定性,因此需要周期性的重新初始化水平集函數(shù)為符號距離函數(shù),這會大大影響演化效率,RDRGAC模型與傳統(tǒng)水平集方法相比引入距離正則項(xiàng),不必重復(fù)此步驟,精確搜索輪廓邊緣并大幅度縮短迭代時間。定義水平集函數(shù)為φ為

    初始水平集函數(shù)定義如下

    (5)

    式中,c0>0,是一個常量;R0是圖像域Ω中的區(qū)域。一般情況下,c0≥1,通常設(shè)置c0=2。所提距離正則化水平集能量函數(shù)定義為

    ε(φ)=μRp(φ)+λLg(φ)+αAg(φ)

    (6)

    式中,λ、α都是能量項(xiàng)系數(shù),λ>0,α∈R,當(dāng)初始輪廓設(shè)置在目標(biāo)區(qū)外部,區(qū)域面積項(xiàng)系數(shù)α>0,輪廓線在迭代過程中收縮至岸線邊緣,當(dāng)初始輪廓設(shè)置在目標(biāo)區(qū)內(nèi)部,系數(shù)α<0,輪廓線在迭代過程中擴(kuò)張至岸線邊緣,μ>0是常量。Rp(φ)是水平集正則項(xiàng),定義如下

    (7)

    (8)

    Ag(φ)=∫ΩgH(-φ)dΧ

    (9)

    式中,Lg(φ)計(jì)算沿著φ的零水平集的函數(shù)g的線性積分,當(dāng)φ的零水平集輪廓到達(dá)目標(biāo)邊緣時,能量函數(shù)Lg(φ)會呈現(xiàn)最小化,Ag(φ)計(jì)算權(quán)重區(qū)域面積Ωint:φ(x,y)<0。g為邊緣指示器函數(shù),定義為

    (10)

    (11)

    (12)

    (13)

    (14)

    (15)

    α∫ΩgHε(-φ)dX

    (16)

    通過變分法求解能量泛函極值,得到關(guān)于φ的梯度矢量流如下

    (17)

    水平集函數(shù)更新方程

    (18)

    3 RDRGAC模型改進(jìn)

    在RDRGAC模型中,區(qū)域面積項(xiàng)系數(shù)α與岸線邊緣強(qiáng)弱有關(guān),通過建立參數(shù)α與ENL的關(guān)系,自動設(shè)定參數(shù)α,從而在盡量減少相干斑噪聲影響的同時,最大程度逼近海岸線。由于參數(shù)α設(shè)置與圖像噪聲強(qiáng)度有關(guān),本文通過使不同噪聲的SAR圖像的ENL同相應(yīng)α建立相關(guān)關(guān)系,實(shí)現(xiàn)α值自動設(shè)定。

    等效視數(shù)ENL度量了圖像區(qū)分不同后向散射特性區(qū)域的能力,是衡量一幅SAR圖像斑點(diǎn)噪聲相對強(qiáng)度的一種指標(biāo),等效視數(shù)越大,表明圖像上斑點(diǎn)越弱。定義為[19-21]

    (19)

    式中,μ1為像素均值;σ1為像素標(biāo)準(zhǔn)差。

    以下是不同區(qū)域的SAR圖像ENL,以及在海岸線自動提取過程中設(shè)定的α值,如圖1所示。

    圖1 針對SAR不同噪聲強(qiáng)度的α取值Fig.1 Different values of α in different SAR image

    根據(jù)多景SAR圖像的統(tǒng)計(jì)試驗(yàn),建立α值與相應(yīng)圖像ENL的非線性擬合關(guān)系,如圖2所示。

    圖2 參數(shù)α取值與ENL擬合關(guān)系Fig.2 Regression relation between α and ENL

    本文選用8景不同區(qū)域SAR圖像ENL來研究ENL同噪聲參數(shù)α關(guān)系以測定參數(shù)α的變化趨勢。并選用其中4景典型特例來說明上述關(guān)系。由圖2可以看出,ENL與參數(shù)α建立的相關(guān)關(guān)系接近對數(shù)關(guān)系,即

    (20)

    對大多數(shù)SAR圖像來講,參數(shù)α取值范圍為0<α≤5。針對弱邊緣圖像即ENL≥7.6,取0<α≤1.5,針對強(qiáng)邊緣圖像即ENL<7.6,取1.5<α≤5之間。由于本文獲得SAR圖像樣本數(shù)量的關(guān)系,根據(jù)ENL求得的α值會有少許誤差,隨著樣本數(shù)量的增加,α與ENL會建立更加精確的對應(yīng)關(guān)系。

    4 試驗(yàn)結(jié)果與分析

    為了驗(yàn)證所提改進(jìn)的區(qū)域距離正則化幾何主動輪廓模型對海岸線自動提取的有效性,將本文所提模型與主動輪廓模型-梯度矢量流模型(gradientvectorflow,GVF-Snake)模型[5]、C-V模型[18]和區(qū)域面積項(xiàng)系數(shù)α人工設(shè)定模型進(jìn)行比較,采用真實(shí)SAR遙感圖像進(jìn)行試驗(yàn)。

    GVF-Snake參數(shù)設(shè)置為:α=0.5,β=0,γ=1,κ=0.6。C-V模型參數(shù)設(shè)置為:λ1=λ2=1,Δt=1.0,ε=1,μ=0.02×2552。本文方法參數(shù)設(shè)置λ=5,Δt=5,ε=1.5,μ=0.2。參數(shù)α原始值為α=1.5。

    4.1河北省北戴河海岸線自動提取試驗(yàn)

    本試驗(yàn)采用河北省北戴河SAR遙感圖像,大小為289×297像素,HV極化模式,如圖3所示。為了保證岸線提取結(jié)果的穩(wěn)健性,初始輪廓均設(shè)在岸線邊緣,本試驗(yàn)初始輪廓為(1:260,1:280),即像素為260×280的長方形區(qū)域。

    SAR圖像同可見光遙感圖像相比,最大特點(diǎn)是存在相干斑噪聲,即SAR圖像發(fā)射相干電磁波,因此各理想點(diǎn)目標(biāo)回波是相互干涉的,相干電磁波照射實(shí)際目標(biāo)時,其散射回來的總回波并不完全由地物目標(biāo)的散射系數(shù)決定,而是圍繞這些散射系數(shù)有很大的隨機(jī)起伏。由于海陸信號的混雜作用以及大陸架暗礁區(qū)域的影響,在海岸線提取過程中極易受到相干斑噪聲影響。

    圖3(a)為北戴河岸線初始輪廓,圖3(b)為GVF-Snake模型岸線分割結(jié)果,其從圖像中獲得一個梯度向量場,即梯度矢量流,依賴于梯度信息作為曲線演化的停止條件,該方法對于噪聲比較敏感,出現(xiàn)了岸線誤分割現(xiàn)象,在SAR圖像中,相干斑噪聲使岸線邊緣比較模糊,在同質(zhì)區(qū)域中也存在一定的灰度起伏和斑點(diǎn)噪聲,造成岸線邊緣處梯度不一定最大,能量泛函在邊緣處不一定達(dá)到最小值,是造成岸線許多誤分割的主因。圖3(c)為C-V模型提取結(jié)果,輪廓線不能精確收斂于真實(shí)岸線邊緣,主要原因在于C-V模型的區(qū)域能量項(xiàng)都是基于歐氏距離,在乘性相干斑噪聲下,歐氏距離不適用于度量SAR圖像像素的相似度,以致在能量最小化過程中易陷入局部極小值。因此造成岸線少許誤分割。圖3(d)是固定參數(shù)α=1.5的RDRGAC方法分割結(jié)果,岸線基本收斂于岸線邊緣,但有少許“冒頂”現(xiàn)象,即輪廓線穿越真實(shí)岸線,由于引入了距離正則項(xiàng),使其具有保證水平集函數(shù)穩(wěn)定演變的特性,不必周期性地重新初始化水平集函數(shù)為符號距離函數(shù),大幅度減少迭代次數(shù)和計(jì)算時間,但由于參數(shù)α取值與岸線強(qiáng)弱有關(guān),針對較弱岸線參數(shù)α不能自行設(shè)定,導(dǎo)致“冒頂”現(xiàn)象。圖3(e)中本文方法分割結(jié)果,在原始方法基礎(chǔ)上,引入了ENL,用于測定SAR圖像噪聲強(qiáng)度,并與參數(shù)α建立非線性擬合關(guān)系,對參數(shù)α自動設(shè)定,實(shí)現(xiàn)了海岸線的精確檢測。圖3(f)為圖3(e)最終水平集。

    為了驗(yàn)證本文方法在SAR圖像海岸線提取中的有效性,各模型迭代次數(shù)均為100次,性能對比如表1所示。從表1中可以發(fā)現(xiàn)GVF-Snake分割速度快,迭代時間短,但對岸線噪聲比較敏感,出現(xiàn)了岸線誤分割,傳統(tǒng)幾何主動輪廓C-V模型由于引入了水平集函數(shù),需要在曲線演化中進(jìn)行耗時的重新初始化步驟,計(jì)算效率低,本文方法引入了距離正則項(xiàng),不必重新初始化水平集函數(shù)為符號距離函數(shù),提升了演化速率,在經(jīng)過ENL自動設(shè)定參數(shù)α后,岸線提取精度方面也有了明顯提高。

    圖3 北戴河岸線提取結(jié)果Fig.3 Coastline extraction result of Beidaihe SAR image

    Tab.1Efficiencies comparison from four segmentation approaches of Beidaihe SAR image

    分割算法迭代時間/s收斂效果GVF-Snake[5]6.07不好,含噪聲較多C-V[18]28.10好,含較少噪聲本文方法(α為原始值)11.26較好,但輪廓線在濾過噪聲的同時越過了岸線邊緣本文方法(α自動設(shè)定)10.96好,輪廓線濾過噪聲剛好收斂于岸線邊緣

    4.2大連市金州灣海岸線自動提取試驗(yàn)

    該試驗(yàn)采用大連市金州灣SAR圖像,大小為707×782像素,HV極化模式,如圖4所示。為了保證岸線提取結(jié)果的抗差性,初始輪廓均設(shè)在岸線邊緣,本試驗(yàn)初始輪廓為(5:705,10:780),即像素為700×770的長方形區(qū)域。

    SAR圖像模糊邊界分為兩類:一是真實(shí)岸線的模糊邊緣,二是由SAR圖像固有的乘性相干斑造成的虛假岸線,在SAR圖像中呈現(xiàn)為弱化邊緣。在海陸信號的共同作用及大陸架暗礁區(qū)域的影響下,第2類情況造成的虛假岸線極易被噪聲敏感的模型誤分割為真實(shí)岸線,從而造成真實(shí)岸線的邊界遺漏,相干斑噪聲導(dǎo)致邊界模糊主要是第2類原因造成的。

    圖4(a)為金州灣岸線初始輪廓,圖4(b)為GVF-Snake方法岸線提取結(jié)果,可以看出,相比于噪聲含量低的SAR圖像,岸線提取效果不佳的特性更加明顯,因?yàn)閰?shù)主動輪廓模型用參數(shù)方程表示輪廓曲線,在二維平面內(nèi)處理參數(shù)曲線的演化,很難處理曲線的拓?fù)渥兓?分裂、合并)等,并且對初始輪廓位置敏感。圖4(c)為傳統(tǒng)幾何主動輪廓C-V模型岸線提取結(jié)果,引入了水平集函數(shù)解決輪廓線拓?fù)浞至押喜⒌膯栴},但由于在SAR圖像中乘性相干斑噪聲不可避免,對噪聲不均勻圖像仍會陷入局部極小值。圖4(d)是固定參數(shù)α=1.5的RDRGAC方法分割結(jié)果,中心區(qū)域?yàn)榇筮B市金州灣飛機(jī)場,對于強(qiáng)噪聲區(qū)的岸線提取效果并不理想,這是因?yàn)樵趶?qiáng)噪聲區(qū),相干斑噪聲使原本具有常數(shù)后向散射系數(shù)均勻同質(zhì)區(qū)域變得不均勻,出現(xiàn)許多斑點(diǎn),固定參數(shù)α將噪聲區(qū)默認(rèn)為岸線邊緣,經(jīng)試驗(yàn)表明,經(jīng)過更多的迭代次數(shù),仍不能精確檢測真實(shí)岸線。圖4(e)中本文方法參數(shù)α經(jīng)自動設(shè)定分割結(jié)果,精確收斂于岸線邊緣。試驗(yàn)證明,參數(shù)α經(jīng)自動設(shè)定能根據(jù)SAR圖像噪聲強(qiáng)度,有效濾過相干斑噪聲,收斂于岸線邊緣,實(shí)現(xiàn)模型自適應(yīng)調(diào)整。圖4(f)為圖4(e)最終水平集。

    圖4 金州灣岸線提取結(jié)果Fig.4 Coastline extraction result of Jinzhou Bay SAR image

    大連市金州灣海岸線提取4種分割方法效率對比如表2所示,各模型迭代次數(shù)均為400次。相比于北戴河SAR圖像,金州灣SAR圖像受相干斑噪聲影響更加嚴(yán)重,對比結(jié)果表明,本文方法對乘性相干斑噪聲的抗差性更強(qiáng),對海岸線輪廓的收斂速度更快,且對初始輪廓不敏感,經(jīng)ENL自動設(shè)定參數(shù)的本文方法對不同噪聲強(qiáng)度SAR圖像的海岸線提取針對性更強(qiáng),且有一定的實(shí)用性。

    表2金州灣SAR圖像4種分割方法性能比較

    Tab.2Efficiencies comparison from four segmentation approaches of Jinzhou Bay SAR image

    分割算法迭代時間/s收斂效果GVF-Snake[5]181.26不好,含噪聲較多C-V[18]571.88不好,含較少噪聲,漏掉部分岸線本文方法(α為原始值)231.01較好,但輪廓線默認(rèn)部分強(qiáng)噪聲區(qū)為岸線邊緣本文方法(α自動設(shè)定)225.27好,輪廓線濾過噪聲剛好收斂于岸線邊緣

    5 結(jié) 語

    本文針對SAR圖像中相干斑噪聲和海表隨機(jī)粗糙度造成海岸目標(biāo)與海水背景模糊的特點(diǎn),提出了基于區(qū)域距離正則化幾何主動輪廓模型,實(shí)現(xiàn)了SAR圖像海岸線自動提取。為避免重新初始化水平集函數(shù)為符號距離函數(shù)的耗時步驟,引入距離正則項(xiàng),提高算法穩(wěn)定性及收斂速度。并將區(qū)域面積項(xiàng)系數(shù)α與SAR圖像ENL建立了非線性擬合關(guān)系,實(shí)現(xiàn)了對包含不同噪聲水平SAR圖像的參數(shù)α自動設(shè)定。通過多種主動輪廓模型方法對河北省北戴河和大連市金州灣SAR數(shù)據(jù)做試驗(yàn)對比,表明本文方法在受相干斑噪聲影響的SAR圖像岸線提取方面有顯著的優(yōu)越性。

    [1]馬小峰, 趙冬至, 張豐收, 等. 海岸線衛(wèi)星遙感提取方法研究進(jìn)展[J]. 遙感技術(shù)與應(yīng)用, 2007, 22(4): 575-580.

    MA Xiaofeng, ZHAO Dongzhi, ZHANG Fengshou, et al. An Overview of Means of Withdrawing Coastline by Remote Sensing[J]. Remote Sensing Technology and Application, 2007, 22(4): 575-580.

    [2]嚴(yán)海兵, 李秉柏, 陳敏東. 遙感技術(shù)提取海岸線的研究進(jìn)展[J]. 地域研究與開發(fā), 2009, 28(1): 101-105.

    YAN Haibing, LI Bingbai, CHEN Mindong. Progress of Researches in Coastline Extraction Based on RS Technique[J]. Areal Research and Development, 2009, 28(1): 101-105.

    [3]李雨謙, 皮亦鳴, 王金峰. 基于水平集的SAR圖像與光學(xué)圖像的配準(zhǔn)[J]. 測繪學(xué)報(bào), 2010, 39(3): 276-282.

    LI Yuqian, PI Yiming, WANG Jinfeng. The Registration between SAR and Optical Image Based on Level Set[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(3): 276-282.

    [4]周奇年, 王廷波, 李文書. 區(qū)域信息和水平集方法的圖像分割[J]. 中國圖象圖形學(xué)報(bào), 2011, 16(11): 2002-2008.

    ZHOU Qinian, WANG Tingbo, LI Wenshu. Image Segmentation Based on Region and Level Set Method[J]. Journal of Image and Graphics, 2011, 16(11): 2002-2008.

    [5]胡倩, 張豐, 劉仁義, 等. 基于梯度矢量流和主動輪廓模型的海陸邊界提取[J]. 浙江大學(xué)學(xué)報(bào)(理學(xué)版), 2013, 40(6): 698-704.

    HU Qian, ZHANG Feng, LIU Renyi, et al. Land Sea Boundary Extraction Based on Gradient Vector Flow and Active Contour Model[J]. Journal of Zhejiang University (Science Edition), 2013, 40(6): 698-704.

    [6]孫陽光, 蔡超, 周成平, 等. R-Snake: 一種基于邊緣與區(qū)域信息的圖像主動輪廓模型提取[J]. 電子學(xué)報(bào), 2009, 37(8): 1810-1815.

    SUN Yangguang, CAI Chao, ZHOU Chengping, et al. R-Snake: A Snake Model Using Both Boundary and Region Information[J]. Acta Electronica Sinica, 2009, 37(8): 1810-1815.

    [7]SUN Guangcai, XING Mengdao, XIA Xianggen, et al. Multichannel Full-aperture Azimuth Processing for Beam Steering SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2013, 51(9): 4761-4778.

    [8]劉鵬程. 基于水平集理論的海岸線輪廓特征提取[J]. 國土資源遙感, 2015, 27(2): 75-79.

    LIU Pengcheng. Feature Extraction of Coastline Contour Based on Level Set Theory[J]. Remote Sensing for Land & Resources, 2015, 27(2): 75-79.

    [9]YANG Xi, GAO Xinbo, TAO Dacehng, et al. An Efficient MRF Embedded Level Set Method for Image Segmentation[J]. IEEE Transactions on Image Processing, 2015, 24(1): 9-21.

    [10]歐陽越, 種勁松. SAR圖像海岸線檢測算法綜述[J]. 國土資源遙感, 2006, 18(2): 1-3.

    OUYANG Yue, ZHONG Jinsong. A Review of Coastline Detection on SAR Imagery[J]. Remote Sensing for Land & Resources, 2006, 18(2): 1-3.

    [11]涂松, 李禹, 粟毅. 基于主動輪廓模型的SAR圖像分割方法綜述[J]. 系統(tǒng)工程與電子技術(shù), 2015, 37(8): 1754-1766.

    TU Song, LI Yu, SU Yi. Overview of SAR Image Segmentation Based on Active Contour Model[J]. Systems Engineering and Electronics, 2015, 37(8): 1754-1766.

    [12]LI Chunming, XU Chenyang, GUI Changfeng, et al. Distance Regularized Level Set Evolution and Its Application to Image Segmentation[J]. IEEE Transactions on Image Processing, 2010, 19(12): 3243-3254.

    [13]COLGAN T J, HAGNESS S C, VAN VEEN B D. A 3-D Level Set Method for Microwave Breast Imaging[J]. IEEE Transactions on Biomedical Engineering, 2015, 62(10): 2526-2534.

    [14]任繼軍, 何明一. 基于粒子群優(yōu)化算法的改進(jìn)Snake模型的圖像分割方法[J]. 中國圖像圖形學(xué)報(bào), 2008, 13(9): 1727-1732.REN Jijun, HE Mingyi. Image Segmentation Using Improved Snake Model Based on Particle Swarm Optimization[J]. Journal of Image and Graphics, 2008, 13(9): 1727-1732.

    [15]東野長磊, 鄭永果, 蘇杰. 基于Chan-Vese模型的SAR圖像分割[J]. 數(shù)據(jù)采集與處理, 2012, 27(2): 151-155.DONGYE Changlei, ZHENG Yongguo, SU Jie. SAR Image Segmentation Based on Chan-Vese Model[J]. Journal of Data Acquisition & Processing, 2012, 27(2): 151-155.

    [16]王相海, 李明. 雙重輪廓演化曲線的圖像分割水平集模型[J]. 中國圖象圖形學(xué)報(bào), 2014, 19(3): 373-380.

    WANG Xianghai, LI Ming. Level Set Model for Image Segmentation Based on Dual Contour Evolutional Curve[J]. Journal of Image and Graphics, 2014, 19(3): 373-380.

    [17]黃魁華, 張軍. 局部統(tǒng)計(jì)活動輪廓模型的SAR圖像海岸線檢測[J]. 遙感學(xué)報(bào), 2011, 15(4): 737-749.

    HUANG Kuihua, ZHANG Jun. A Coastline Detection Method Using SAR Images Based on the Local Statistical Active Contour Model[J]. Journal of Remote Sensing, 2011, 15(4): 737-749.

    [18]李傳龍, 李穎, 蘭國新. 一種新的窄帶快速區(qū)域水平集C-V模型圖像分割方法[J]. 計(jì)算機(jī)科學(xué), 2011, 38(12): 17-19, 35.

    LI Chuanlong, LI Ying, LAN Guoxin. New Method of the Fast Narrow Brand C-V Level Set Model for Image Segmentation[J]. Computer Science, 2011, 38(12): 17-19, 35.

    [19]吳良斌. SAR圖像處理與目標(biāo)識別[M]. 北京: 航空工業(yè)出版社, 2013: 21-22.

    WU Liangbin. SAR Image Processing and Target Recognition[M]. Beijing: Aviation Industry Press, 2013: 21-22.

    [20]崔瑞, 薛磊, 汪波. 基于等效視數(shù)的ISAR干擾效果評估方法[J]. 系統(tǒng)工程與電子技術(shù), 2008, 30(5): 887-888.

    CUI Rui, XUE Lei, WANG Bo. Evaluation Method of Jamming Effect on ISAR Based on Equivalent Number of Looks[J]. Systems Engineering and Electronics, 2008, 30(5): 887-888.

    [21]王曉軍, 孫洪. SAR圖像相干斑濾波性能評估模型[J]. 雷達(dá)科學(xué)與技術(shù), 2008, 6(1): 29-34, 38.

    WANG Xiaojun, SUN Hong. Measurement Model and Performance Evaluation for Speckle Filter of SAR Images[J]. Radar Science and Technology, 2008, 6(1): 29-34, 38.

    (責(zé)任編輯:張艷玲)

    SAR Image Coastline Detection Based on Regional Distance Regularized Geometric Active Contour Models

    JIANG Dawei1, 2,FAN Jianchao2,HUANG Fengrong1

    1. School of Urban Environment, Liaoning Normal University, Dalian 116029, China; 2. Marine Resources and Environment Monitoring Center, National Marine Environmental Monitoring Center, Dalian 116023, China

    Synthetic aperture radar (SAR) satellite remote sensing images can greatly increase the frequency of the coastline coverage all over the country. However, due to the influence of the random sea surface roughness caused by waves, it’s a challenge to distinguish the coastline and sea boundary. To solve this problem, this paper proposes regional distance regularized geometric active contour models (RDRGAC), in which the distance regularized term is introduced to avoid periodically initializing the degraded function with a signed distance function, accelerating the speed of convergence. Besides, it establishes the nonlinear regression relationship between the regional term parameters and equivalent number of looks (ENL) in SAR images, leading to the adaptive setting of RDRGAC model with different SAR images, which could improve the accuracy of coastline automatic detection in return. In the experiments, SAR images in Beidaihe and Dalian Jinzhou Bay respectively are adopted to detect the coastline, verifying the effective of the proposed method.

    speckle noise; synthetic aperture radar; nonlinear regression relationship; geometric active contour model; equivalent number of looks

    The National Natural Science Foundation of China (No.61273307); The China Postdoctoral Science Foundation (No.2014M551082); High Resolution Special Research (No.41-Y30B12-9001-14/16).

    JIANG Dawei(1991—),male,postgraduate,majors in marine RS SAR image processing, and marine area use dynamic monitoring.

    FAN Jianchao

    姜大偉,范劍超,黃鳳榮.SAR圖像海岸線檢測的區(qū)域距離正則化幾何主動輪廓模型[J].測繪學(xué)報(bào),2016,45(9):1096-1103.

    10.11947/j.AGCS.2016.20160039.

    JIANG Dawei, FAN Jianchao, HUANG Fengrong.SAR Image Coastline Detection Based on Regional Distance Regularized Geometric Active Contour Models[J]. Acta Geodaetica et Cartographica Sinica,2016,45(9):1096-1103. DOI:10.11947/j.AGCS.2016.20160039.

    P237

    A

    1001-1595(2016)09-1096-08

    國家自然科學(xué)基金(61273307);中國博士后面上基金(2014M551082);高分辨率對地觀測系統(tǒng)重大專項(xiàng)(41-Y30B12-9001-14/16)

    2016-02-24

    姜大偉(1991—),男,碩士生,研究方向?yàn)楹Q筮b感SAR圖像處理、海域使用動態(tài)監(jiān)測。

    E-mail: 18840817436@163.com

    范劍超

    E-mail: jcfan@nmemc.org.cn

    修回日期: 2016-06-16

    猜你喜歡
    海岸線輪廓邊緣
    OPENCV輪廓識別研究與實(shí)踐
    基于實(shí)時輪廓誤差估算的數(shù)控系統(tǒng)輪廓控制
    徒步拍攝英國海岸線
    徒步拍攝英國海岸線
    徒步拍攝英國海岸線
    一張圖看懂邊緣計(jì)算
    徒步拍攝英國海岸線
    在線學(xué)習(xí)機(jī)制下的Snake輪廓跟蹤
    創(chuàng)造早秋新輪廓
    舒適廣告(2008年9期)2008-09-22 10:02:48
    在邊緣尋找自我
    雕塑(1999年2期)1999-06-28 05:01:42
    亚洲国产av新网站| 亚洲精品国产色婷婷电影| 色视频在线一区二区三区| 中文字幕av电影在线播放| 欧美日韩国产mv在线观看视频| 一二三四在线观看免费中文在| 精品国产超薄肉色丝袜足j| av在线app专区| 人人妻人人爽人人添夜夜欢视频| 亚洲国产欧美一区二区综合| 国产一区亚洲一区在线观看| 中文字幕精品免费在线观看视频| 新久久久久国产一级毛片| 99国产精品免费福利视频| 51午夜福利影视在线观看| 日韩伦理黄色片| 男女边吃奶边做爰视频| 色婷婷久久久亚洲欧美| 丝袜喷水一区| 成在线人永久免费视频| 视频区图区小说| 成人影院久久| 又紧又爽又黄一区二区| 国产免费又黄又爽又色| 一级毛片黄色毛片免费观看视频| 国产麻豆69| 久久久精品94久久精品| 欧美日韩av久久| 国产野战对白在线观看| 久久人妻熟女aⅴ| 99热国产这里只有精品6| 久久午夜综合久久蜜桃| 日韩 亚洲 欧美在线| 丝袜喷水一区| 国产不卡av网站在线观看| 十八禁网站网址无遮挡| 久久精品国产亚洲av高清一级| 一级,二级,三级黄色视频| 亚洲精品乱久久久久久| 久久精品亚洲熟妇少妇任你| 日韩视频在线欧美| 欧美精品一区二区大全| 欧美老熟妇乱子伦牲交| 在线观看免费高清a一片| 一级a爱视频在线免费观看| 亚洲国产毛片av蜜桃av| 99国产综合亚洲精品| 美女主播在线视频| 精品一区在线观看国产| 一边摸一边做爽爽视频免费| 午夜福利,免费看| 日韩一本色道免费dvd| 成年女人毛片免费观看观看9 | 欧美日韩精品网址| 国产精品久久久人人做人人爽| 亚洲欧美一区二区三区黑人| 人人妻人人澡人人看| 亚洲精品av麻豆狂野| 热re99久久精品国产66热6| 考比视频在线观看| 十八禁人妻一区二区| 国产精品三级大全| 狂野欧美激情性bbbbbb| 宅男免费午夜| 久久99精品国语久久久| 久久久久精品国产欧美久久久 | 精品少妇一区二区三区视频日本电影| 高清视频免费观看一区二区| 麻豆国产av国片精品| 亚洲av电影在线观看一区二区三区| 午夜激情av网站| 久久亚洲精品不卡| 免费高清在线观看视频在线观看| 多毛熟女@视频| 男人爽女人下面视频在线观看| tube8黄色片| 色网站视频免费| 亚洲欧洲国产日韩| 2021少妇久久久久久久久久久| 国产日韩欧美亚洲二区| 满18在线观看网站| 激情五月婷婷亚洲| 日韩 亚洲 欧美在线| 考比视频在线观看| 观看av在线不卡| 色视频在线一区二区三区| 亚洲精品国产区一区二| 国产一区二区三区综合在线观看| 一区在线观看完整版| 天堂8中文在线网| 亚洲三区欧美一区| 中文字幕人妻丝袜一区二区| 亚洲欧洲国产日韩| 老司机影院成人| 午夜福利视频在线观看免费| 国产成人精品无人区| 精品欧美一区二区三区在线| 国产精品99久久99久久久不卡| 精品国产一区二区三区四区第35| 天天躁夜夜躁狠狠久久av| 国产精品免费大片| 日韩一区二区三区影片| 日本黄色日本黄色录像| 女警被强在线播放| 国产精品熟女久久久久浪| 国产xxxxx性猛交| cao死你这个sao货| 久久ye,这里只有精品| 久久久久久亚洲精品国产蜜桃av| 一本大道久久a久久精品| 日本91视频免费播放| 久久精品人人爽人人爽视色| 中文精品一卡2卡3卡4更新| 精品久久久精品久久久| 在线观看人妻少妇| 欧美变态另类bdsm刘玥| videosex国产| 老鸭窝网址在线观看| 桃花免费在线播放| 亚洲欧美中文字幕日韩二区| 欧美黄色片欧美黄色片| 精品熟女少妇八av免费久了| 一区二区三区乱码不卡18| 极品少妇高潮喷水抽搐| 香蕉丝袜av| 一本色道久久久久久精品综合| 亚洲欧美清纯卡通| 80岁老熟妇乱子伦牲交| 成人黄色视频免费在线看| 一级毛片女人18水好多 | 久久久久视频综合| 国产黄色视频一区二区在线观看| 中文字幕人妻丝袜一区二区| 桃花免费在线播放| 秋霞在线观看毛片| 国产亚洲欧美在线一区二区| 可以免费在线观看a视频的电影网站| 女性被躁到高潮视频| 亚洲精品久久成人aⅴ小说| 丝袜人妻中文字幕| 国产成人免费无遮挡视频| www.自偷自拍.com| 亚洲一码二码三码区别大吗| 一级毛片女人18水好多 | 日韩av免费高清视频| 水蜜桃什么品种好| 中文字幕人妻丝袜一区二区| 菩萨蛮人人尽说江南好唐韦庄| 久久精品国产综合久久久| 亚洲成人免费av在线播放| 男女之事视频高清在线观看 | 亚洲av成人不卡在线观看播放网 | 少妇人妻久久综合中文| 少妇 在线观看| 一本一本久久a久久精品综合妖精| 午夜av观看不卡| 日日爽夜夜爽网站| 成年动漫av网址| 亚洲精品久久久久久婷婷小说| xxxhd国产人妻xxx| 亚洲国产毛片av蜜桃av| 欧美成狂野欧美在线观看| 国产99久久九九免费精品| www.自偷自拍.com| 亚洲一码二码三码区别大吗| 久久综合国产亚洲精品| 91字幕亚洲| 中文字幕人妻丝袜制服| 亚洲欧洲精品一区二区精品久久久| 国产精品国产三级国产专区5o| 99国产精品免费福利视频| 婷婷色综合大香蕉| 欧美性长视频在线观看| 国产欧美日韩一区二区三 | 大型av网站在线播放| 成人影院久久| 久久久国产欧美日韩av| 免费不卡黄色视频| av不卡在线播放| 国产成人系列免费观看| 日韩av不卡免费在线播放| 精品一区在线观看国产| www.自偷自拍.com| 亚洲精品久久久久久婷婷小说| 久久青草综合色| 亚洲伊人色综图| 男女之事视频高清在线观看 | 精品人妻一区二区三区麻豆| 欧美精品一区二区大全| 男男h啪啪无遮挡| 乱人伦中国视频| 在线看a的网站| 激情五月婷婷亚洲| 精品亚洲乱码少妇综合久久| 国产黄色免费在线视频| 丝袜美足系列| 午夜免费观看性视频| 18在线观看网站| 看免费av毛片| 男女无遮挡免费网站观看| 亚洲一区二区三区欧美精品| 亚洲成人手机| kizo精华| 99热国产这里只有精品6| 色婷婷久久久亚洲欧美| 黄色视频不卡| 大片电影免费在线观看免费| 999久久久国产精品视频| 国产伦人伦偷精品视频| 人人澡人人妻人| 午夜视频精品福利| 免费观看人在逋| 午夜久久久在线观看| 婷婷成人精品国产| 亚洲国产精品一区二区三区在线| 久久国产精品人妻蜜桃| www.自偷自拍.com| 欧美日韩成人在线一区二区| 王馨瑶露胸无遮挡在线观看| 欧美黑人精品巨大| 亚洲欧洲国产日韩| 欧美成人午夜精品| av有码第一页| 欧美日韩亚洲国产一区二区在线观看 | 久久影院123| 男人舔女人的私密视频| 中国美女看黄片| av国产精品久久久久影院| 999久久久国产精品视频| 中文字幕另类日韩欧美亚洲嫩草| 在线观看国产h片| 每晚都被弄得嗷嗷叫到高潮| 男女下面插进去视频免费观看| 校园人妻丝袜中文字幕| av在线老鸭窝| 汤姆久久久久久久影院中文字幕| 精品一区二区三区av网在线观看 | 熟女少妇亚洲综合色aaa.| 国产精品偷伦视频观看了| av天堂久久9| 黄网站色视频无遮挡免费观看| 一二三四社区在线视频社区8| 亚洲国产av新网站| 在线观看www视频免费| 叶爱在线成人免费视频播放| 黄色a级毛片大全视频| av在线app专区| 麻豆av在线久日| av在线播放精品| 国产野战对白在线观看| 人人妻人人澡人人看| 777久久人妻少妇嫩草av网站| 亚洲男人天堂网一区| 色婷婷av一区二区三区视频| 亚洲熟女毛片儿| 欧美日韩综合久久久久久| 亚洲情色 制服丝袜| 国产亚洲精品第一综合不卡| 两性夫妻黄色片| 国语对白做爰xxxⅹ性视频网站| 一区二区日韩欧美中文字幕| 日本午夜av视频| 亚洲国产精品成人久久小说| 精品国产乱码久久久久久男人| 久久影院123| 精品人妻1区二区| 男的添女的下面高潮视频| 高潮久久久久久久久久久不卡| 新久久久久国产一级毛片| 免费人妻精品一区二区三区视频| 亚洲欧美中文字幕日韩二区| 老司机影院毛片| 欧美久久黑人一区二区| h视频一区二区三区| 一区二区av电影网| 国产一级毛片在线| 亚洲精品久久久久久婷婷小说| 人妻人人澡人人爽人人| 国产成人精品在线电影| 我的亚洲天堂| 久久av网站| 色综合欧美亚洲国产小说| 中文字幕亚洲精品专区| 日韩大码丰满熟妇| 一区二区三区乱码不卡18| 在线天堂中文资源库| 午夜久久久在线观看| 国产在线一区二区三区精| 中文欧美无线码| 亚洲精品日本国产第一区| 欧美激情极品国产一区二区三区| 高潮久久久久久久久久久不卡| 男女无遮挡免费网站观看| 国产成人一区二区三区免费视频网站 | 亚洲第一av免费看| 一区在线观看完整版| 97人妻天天添夜夜摸| 叶爱在线成人免费视频播放| 国产成人a∨麻豆精品| 少妇被粗大的猛进出69影院| 国产男女内射视频| 国产成人系列免费观看| 亚洲av男天堂| 亚洲图色成人| 成年人黄色毛片网站| 亚洲精品一卡2卡三卡4卡5卡 | 伊人亚洲综合成人网| 男男h啪啪无遮挡| 欧美黄色片欧美黄色片| 国产视频首页在线观看| 日韩 欧美 亚洲 中文字幕| 免费不卡黄色视频| 亚洲国产欧美日韩在线播放| 久久久久久久久久久久大奶| 亚洲自偷自拍图片 自拍| 老汉色av国产亚洲站长工具| 久久综合国产亚洲精品| 赤兔流量卡办理| 一级片'在线观看视频| 少妇粗大呻吟视频| 亚洲av国产av综合av卡| 男人爽女人下面视频在线观看| 美女脱内裤让男人舔精品视频| 亚洲情色 制服丝袜| 久久国产精品男人的天堂亚洲| 午夜av观看不卡| 欧美日韩av久久| 免费观看av网站的网址| av电影中文网址| 国产亚洲欧美在线一区二区| 日本av手机在线免费观看| avwww免费| 亚洲av男天堂| 丝瓜视频免费看黄片| 制服诱惑二区| 亚洲五月色婷婷综合| 欧美日韩综合久久久久久| 一区二区三区四区激情视频| 一二三四社区在线视频社区8| 亚洲精品国产色婷婷电影| 国产精品麻豆人妻色哟哟久久| 麻豆国产av国片精品| 久久久久国产精品人妻一区二区| 久久久久久久久久久久大奶| 亚洲欧美精品综合一区二区三区| 大型av网站在线播放| av国产久精品久网站免费入址| 欧美激情 高清一区二区三区| 国产成人欧美在线观看 | 欧美日韩亚洲国产一区二区在线观看 | 捣出白浆h1v1| 青青草视频在线视频观看| 国产主播在线观看一区二区 | 久热这里只有精品99| 最新在线观看一区二区三区 | 亚洲五月婷婷丁香| 一级a爱视频在线免费观看| 国产精品久久久久久人妻精品电影 | 日本一区二区免费在线视频| 妹子高潮喷水视频| 婷婷色麻豆天堂久久| 老司机深夜福利视频在线观看 | 中文字幕人妻熟女乱码| 亚洲av综合色区一区| 丝袜人妻中文字幕| 午夜两性在线视频| 欧美亚洲 丝袜 人妻 在线| 一级a爱视频在线免费观看| 欧美精品人与动牲交sv欧美| 成年人午夜在线观看视频| 狠狠婷婷综合久久久久久88av| 美女福利国产在线| 五月开心婷婷网| 校园人妻丝袜中文字幕| 国产精品二区激情视频| 色视频在线一区二区三区| 亚洲国产精品999| 成人影院久久| 99re6热这里在线精品视频| 伦理电影免费视频| 宅男免费午夜| 精品少妇久久久久久888优播| 国产欧美日韩一区二区三 | bbb黄色大片| 国产无遮挡羞羞视频在线观看| 男女边摸边吃奶| a级毛片黄视频| 久久精品国产亚洲av高清一级| 亚洲精品国产av成人精品| 人体艺术视频欧美日本| 无限看片的www在线观看| 99国产综合亚洲精品| 天天躁日日躁夜夜躁夜夜| 欧美激情极品国产一区二区三区| 精品卡一卡二卡四卡免费| 久久久国产精品麻豆| 国产免费视频播放在线视频| 色94色欧美一区二区| 男男h啪啪无遮挡| 免费高清在线观看日韩| 99九九在线精品视频| 老司机影院毛片| 一级毛片我不卡| 99精国产麻豆久久婷婷| 欧美精品啪啪一区二区三区 | 亚洲人成电影观看| 丝袜在线中文字幕| 少妇人妻 视频| 亚洲少妇的诱惑av| 777米奇影视久久| 日韩精品免费视频一区二区三区| 欧美精品av麻豆av| 精品人妻熟女毛片av久久网站| 一区二区三区四区激情视频| 老司机影院毛片| 婷婷色综合www| 最近最新中文字幕大全免费视频 | 日韩大码丰满熟妇| 国产老妇伦熟女老妇高清| 美女高潮到喷水免费观看| 国产精品二区激情视频| 精品一区二区三区四区五区乱码 | 亚洲国产欧美日韩在线播放| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久成人av| 日韩,欧美,国产一区二区三区| 午夜久久久在线观看| 亚洲av综合色区一区| 人体艺术视频欧美日本| 欧美日韩亚洲国产一区二区在线观看 | 亚洲av电影在线进入| 日本wwww免费看| 国产成人a∨麻豆精品| 亚洲av成人不卡在线观看播放网 | 国产淫语在线视频| 欧美在线一区亚洲| 国产精品免费大片| 九色亚洲精品在线播放| 18禁国产床啪视频网站| 青草久久国产| 久久国产亚洲av麻豆专区| tube8黄色片| 在线观看免费午夜福利视频| 性色av一级| 国产精品 国内视频| 一本大道久久a久久精品| 九草在线视频观看| 91老司机精品| 久久久精品免费免费高清| 各种免费的搞黄视频| 久久久久久人人人人人| 一区二区三区精品91| 国产老妇伦熟女老妇高清| 国产有黄有色有爽视频| a级毛片黄视频| 婷婷色麻豆天堂久久| 少妇粗大呻吟视频| e午夜精品久久久久久久| 巨乳人妻的诱惑在线观看| 日韩av免费高清视频| 国产欧美日韩一区二区三 | 2021少妇久久久久久久久久久| 亚洲黑人精品在线| 一边摸一边做爽爽视频免费| 亚洲欧美成人综合另类久久久| 久久久国产欧美日韩av| 青春草亚洲视频在线观看| 91成人精品电影| 啦啦啦在线观看免费高清www| 午夜福利乱码中文字幕| 91老司机精品| 国产亚洲欧美在线一区二区| www日本在线高清视频| 久久久久久久久久久久大奶| 欧美成人精品欧美一级黄| 亚洲av成人精品一二三区| 精品人妻1区二区| 2018国产大陆天天弄谢| 蜜桃在线观看..| 50天的宝宝边吃奶边哭怎么回事| 90打野战视频偷拍视频| 欧美性长视频在线观看| 国产精品成人在线| 老熟女久久久| 免费久久久久久久精品成人欧美视频| 国产欧美日韩精品亚洲av| 啦啦啦在线免费观看视频4| 国产亚洲精品久久久久5区| 777久久人妻少妇嫩草av网站| 国产精品成人在线| 丁香六月天网| 亚洲精品国产色婷婷电影| 黄色一级大片看看| 午夜福利影视在线免费观看| 亚洲精品自拍成人| 久久久久久人人人人人| av又黄又爽大尺度在线免费看| 黄片小视频在线播放| 91精品伊人久久大香线蕉| 成年美女黄网站色视频大全免费| 啦啦啦啦在线视频资源| 黄色片一级片一级黄色片| 国产淫语在线视频| 成人影院久久| 大陆偷拍与自拍| 午夜激情av网站| 麻豆国产av国片精品| 欧美av亚洲av综合av国产av| 国产精品一区二区在线不卡| 天天添夜夜摸| 9热在线视频观看99| 久久性视频一级片| 极品少妇高潮喷水抽搐| 久久久久精品国产欧美久久久 | cao死你这个sao货| 亚洲精品在线美女| 女人精品久久久久毛片| 久久亚洲国产成人精品v| 国产99久久九九免费精品| 大话2 男鬼变身卡| 水蜜桃什么品种好| 欧美日韩一级在线毛片| avwww免费| 一级毛片女人18水好多 | 妹子高潮喷水视频| 亚洲免费av在线视频| 欧美黑人精品巨大| 中文字幕高清在线视频| 伊人亚洲综合成人网| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩亚洲国产一区二区在线观看 | 免费高清在线观看视频在线观看| 999久久久国产精品视频| 国产男人的电影天堂91| 看免费成人av毛片| 超碰成人久久| 婷婷色综合大香蕉| 欧美97在线视频| 啦啦啦在线免费观看视频4| 大话2 男鬼变身卡| 国产成人欧美| 熟女少妇亚洲综合色aaa.| 深夜精品福利| 亚洲自偷自拍图片 自拍| 在线观看免费日韩欧美大片| videos熟女内射| 男女高潮啪啪啪动态图| 色网站视频免费| 一边摸一边做爽爽视频免费| 一本久久精品| 性高湖久久久久久久久免费观看| 热re99久久国产66热| 大陆偷拍与自拍| 又黄又粗又硬又大视频| 男女免费视频国产| 亚洲国产精品成人久久小说| 这个男人来自地球电影免费观看| 午夜两性在线视频| 如日韩欧美国产精品一区二区三区| 三上悠亚av全集在线观看| 久久久精品国产亚洲av高清涩受| 亚洲少妇的诱惑av| 各种免费的搞黄视频| 久久精品aⅴ一区二区三区四区| 国产亚洲欧美精品永久| 国产真人三级小视频在线观看| 国产精品免费大片| a级片在线免费高清观看视频| av一本久久久久| 亚洲国产精品999| 国产成人a∨麻豆精品| 香蕉丝袜av| 国产精品免费视频内射| 亚洲天堂av无毛| 国产日韩欧美亚洲二区| 又粗又硬又长又爽又黄的视频| 国产精品香港三级国产av潘金莲 | 一本一本久久a久久精品综合妖精| 久久青草综合色| 波野结衣二区三区在线| av不卡在线播放| 亚洲精品日韩在线中文字幕| 晚上一个人看的免费电影| 极品人妻少妇av视频| 成年动漫av网址| 欧美久久黑人一区二区| 人人妻人人澡人人爽人人夜夜| 高清欧美精品videossex| 满18在线观看网站| 精品国产超薄肉色丝袜足j| 国产国语露脸激情在线看| 日本91视频免费播放| 亚洲欧洲精品一区二区精品久久久| 尾随美女入室| 久久精品国产综合久久久| 99国产精品一区二区蜜桃av | 久久这里只有精品19| 亚洲视频免费观看视频| 老司机深夜福利视频在线观看 | 在线观看免费高清a一片| 亚洲精品自拍成人| 欧美av亚洲av综合av国产av| 日日摸夜夜添夜夜爱| 高清黄色对白视频在线免费看| 黄色 视频免费看| 欧美日韩亚洲国产一区二区在线观看 | 黄色片一级片一级黄色片| 午夜免费男女啪啪视频观看| av不卡在线播放| 亚洲天堂av无毛| videos熟女内射| 国产高清国产精品国产三级| 国产精品九九99| 建设人人有责人人尽责人人享有的| 久久久久久亚洲精品国产蜜桃av|