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

    基于Kirsch算子和多級(jí)閾值分割的河道輪廓識(shí)別和提取算法

    2023-08-02 14:53:32周少君佘海龍
    人民珠江 2023年7期
    關(guān)鍵詞:輪廓算子灰度

    周少君,佘海龍

    (肅南裕固族自治縣水務(wù)局,甘肅 張掖 734400)

    隨著遙感成像技術(shù)的不斷發(fā)展,其在水文、水利等研究領(lǐng)域中的運(yùn)用越來(lái)越廣泛。依據(jù)遙感圖像可以對(duì)水旱災(zāi)害進(jìn)行檢測(cè)和評(píng)估,為水利工程的規(guī)劃建設(shè)與運(yùn)行、水文資源調(diào)查及水質(zhì)監(jiān)測(cè)提供服務(wù)[1-2]。對(duì)水利設(shè)施、水體、河道等遙感圖像進(jìn)行深入分析,有助于排除水利工程中存在的風(fēng)險(xiǎn),保證設(shè)施的完好運(yùn)行[3-4],河道輪廓的提取是其中的重要一環(huán)。然而面對(duì)大數(shù)據(jù)量的遙感圖像,單純依賴(lài)人工分析處理難以滿(mǎn)足要求?,F(xiàn)階段需要借助一些智能化分析及處理手段實(shí)現(xiàn)對(duì)河道的提取。

    對(duì)于遙感圖像的目標(biāo)檢測(cè)智能化分析手段主要分為基于手工設(shè)計(jì)特征和基于深度學(xué)習(xí)兩類(lèi)。前者通過(guò)提取候選區(qū)域,針對(duì)目標(biāo)人工設(shè)計(jì)特征,在結(jié)合分類(lèi)器進(jìn)行分類(lèi)確定目標(biāo)類(lèi)別[5]。這種針對(duì)特定目標(biāo)的特征進(jìn)行目標(biāo)檢測(cè)能夠快速準(zhǔn)確地識(shí)別特定目標(biāo),因此得到了廣泛的應(yīng)用。后者隨著人工神經(jīng)網(wǎng)絡(luò)技術(shù)的快速發(fā)展,能夠提取圖像的深層特征,具有更強(qiáng)的語(yǔ)義表達(dá)和判別性[6]。但是河道的遙感圖像具有尺寸大、形狀變化大、背景復(fù)雜等特點(diǎn),這極大地增加了后者的開(kāi)發(fā)和使用難度,因此本文主要采用前者實(shí)現(xiàn)對(duì)河道輪廓的提取。

    目前針對(duì)遙感圖像中河道輪廓的檢測(cè)和提取算法主要分為兩類(lèi)。一類(lèi)是采用灰度或色彩特征的水體區(qū)域識(shí)別。該方法主要是基于先驗(yàn)知識(shí)的指導(dǎo),認(rèn)為水體的光譜或電磁波反射相對(duì)于背景或其他的目標(biāo)具有明顯差異。因此通過(guò)對(duì)灰度信息或光譜特征建模并尋找出和背景間的差異便能夠識(shí)別出相應(yīng)的河道。但是在復(fù)雜背景的條件下,與水體反射率相近的物體將會(huì)被錯(cuò)誤識(shí)別,缺乏普遍性。另一類(lèi)是依賴(lài)河道本身的形態(tài)學(xué)特征,該方法根據(jù)同背景目標(biāo)間輪廓及邊緣形狀的差異,識(shí)別出相應(yīng)的河道。該方法可以擺脫水體反射率多變對(duì)目標(biāo)提取的影響,適用范圍更加廣泛。但是河道輪廓形態(tài)的復(fù)雜性嚴(yán)重影響到該類(lèi)特征提取的準(zhǔn)確度,在復(fù)雜背景下河道提取效果較差[7]。將兩者相結(jié)合,利用各自?xún)?yōu)點(diǎn)成為解決河道輪廓準(zhǔn)確提取的有效方式。

    灰度閾值分割是最為常用的、實(shí)現(xiàn)簡(jiǎn)單的閾值分割方法,能夠有效地將部分水體與其他地物分割。其中1979年由日本學(xué)者提出OSTU閾值分割法就是一種典型的閾值分割方式。由于其能夠自適應(yīng)的尋找出合適的分割閾值對(duì)圖像進(jìn)行分割,被廣泛使用[8]并不斷改進(jìn)[9]。形態(tài)學(xué)方法是利用圖像的形狀或特征,以便進(jìn)一步進(jìn)行圖像分析和目標(biāo)識(shí)別的方法。常規(guī)的形態(tài)學(xué)處理方法利用結(jié)構(gòu)元對(duì)圖像進(jìn)行膨脹、腐蝕、聯(lián)通等基本操作。形態(tài)學(xué)方法在遙感圖像識(shí)別[10]等方面具有廣泛應(yīng)用。王鵬程等[11]基于河流的形態(tài)學(xué)特征,依據(jù)遙感圖像,對(duì)河流輪廓進(jìn)行識(shí)別提取。朱賀等[7]結(jié)合閾值分割和河流條帶狀輪廓形態(tài),較好的提取出了簡(jiǎn)單河道的輪廓。但是在實(shí)際情況中,河流在彎折處、分叉處等并不是完美的條帶狀,形態(tài)更為復(fù)雜,條帶形假設(shè)并不一定成立。對(duì)于形狀更為復(fù)雜的河流輪廓提取,需要更改形態(tài)學(xué)處理的方式。另外,其處理的河道周?chē)牡匚镄畔⒉⒉皇謴?fù)雜,通過(guò)多級(jí)的閾值分割便可較好解決。但是在實(shí)際情況中,河流周邊可能為城市、農(nóng)田、山區(qū)等,地物信息更為復(fù)雜,背景噪聲需要通過(guò)合理的手段進(jìn)行抑制。

    本文針對(duì)高分辨率圖像(圖像分辨率小于等于10 m)中的河道輪廓識(shí)別和提取過(guò)程中的建模和背景噪聲抑制難等問(wèn)題,提出了基于Kirsch算子和多級(jí)閾值分割的算法。其總體框架見(jiàn)圖1。首先輸入河道遙感圖像,將其轉(zhuǎn)化為灰度圖像。利用Kirsch算子獲取圖像主要地物的邊緣輪廓,再使用OSTU閾值分割算法對(duì)圖像進(jìn)行簡(jiǎn)單的一次分割,從而對(duì)于復(fù)雜背景進(jìn)行抑制。之后依據(jù)河道的主要形態(tài)學(xué)特征提出新的形態(tài)學(xué)判據(jù)對(duì)分割后的圖像進(jìn)行處理,獲取圖像中的河流區(qū)域。最后對(duì)河流區(qū)域再次采用OSTU的閾值分割算法,獲取準(zhǔn)確的河流輪廓,提取出河流區(qū)域。

    圖1 河道輪廓識(shí)別和提取算法總框架

    1 背景模糊處理

    高分辨率遙感圖像中所包含的數(shù)據(jù)是相同面積的中、低分辨率圖像中所包含數(shù)據(jù)的百倍以上,原本較低分辨率圖像中的孤立點(diǎn)在高分辨率圖像中會(huì)變成目標(biāo)。因而高分辨率圖像中將會(huì)包含更多更加豐富的地物信息、紋理及細(xì)節(jié)信息。但是同時(shí)也放大了在中、低分辨率遙感影響中比較弱、甚至可以忽略的干擾噪聲。在河道識(shí)別時(shí),豐富的背景信息將會(huì)對(duì)河道的準(zhǔn)確識(shí)別產(chǎn)生較大干擾。本文將結(jié)合Kirsch算子和閾值分割算法,去除復(fù)雜背景的影響。

    1.1 Kirsch算子

    Kirsch算子是R.Kirsch提出的一種邊緣檢測(cè)算法。對(duì)于包含復(fù)雜背景噪聲的河道遙感圖像,相較于其他傳統(tǒng)算子,Kirsch算子具有更好的適用性。Kirsch算子的原理為:采用8個(gè)模板對(duì)圖像上的每個(gè)像素點(diǎn)進(jìn)行卷積求導(dǎo)。這8個(gè)模板代表8個(gè)方向,計(jì)算出8個(gè)方向的梯度幅值和方向,并以最大的卷積值作為該點(diǎn)的灰度值。圖2為8個(gè)梯度方向,式(1)為采用的Kirsch算子。

    圖2 Kirsch算子的8個(gè)梯度方向

    (1)

    算例均以中國(guó)主要河流及其支流的遙感衛(wèi)星圖像為例,原始圖像均為RGB三通道圖像。利用相關(guān)地圖軟件,全波段的遙感衛(wèi)星地形圖可較為容易取得。長(zhǎng)江某平原地區(qū)河流的原始圖像見(jiàn)圖3a。圖形大小為1 024×1 024,圖像分辨率為5 m。將其灰度化后得到灰度(圖3b)。采用Kirsch算子對(duì)圖像中所有像素點(diǎn)進(jìn)行卷積求導(dǎo)后,得到的主要地物的邊緣圖像見(jiàn)圖3c。由圖可知采用Kirsch算子后,主要地物尤其式河流的邊緣形貌已經(jīng)被提取出來(lái),背景的復(fù)雜程度也得到了降低。

    a)原始RGB三通道衛(wèi)星遙感

    1.2 閾值分割

    采用Kirsch算子對(duì)原始灰度圖進(jìn)行處理后,雖然背景一定程度得到了簡(jiǎn)化,但是仍比較復(fù)雜。此時(shí)圖像為連續(xù)的灰度圖,可以采用閾值分割算法,將背景進(jìn)一步簡(jiǎn)化,并將其轉(zhuǎn)化為二值圖像。本文采用OSTU閾值分割算法,其基本原理是將灰度分布的圖像分為前景和背景兩部分,前景為分割出的部分。分割前景和背景的分隔值就是所求解分割閾值。算法基本原理如下。

    設(shè)圖像灰度為L(zhǎng),灰度級(jí)為i的像素點(diǎn)為ni,則直方圖分布為式(2):

    (2)

    式中,N為總的像素點(diǎn),pi為灰度級(jí)為i的像素點(diǎn)的直方分布。圖像總的灰度值的均值記為ut,見(jiàn)式(3):

    (3)

    按照灰度級(jí)用閾值t將總灰度L劃分為兩類(lèi)C0(0,1,…,t)和C1(t+1,t+2,…,L-1)。則C0和C1部分像素占總像素的比例分別為w0和w1,見(jiàn)式(4)、(5):

    (4)

    w1=1-w0

    (5)

    圖像中C0和C1部分灰度值的均值分別記為u0和u1,見(jiàn)式(6)、(7):

    (6)

    (7)

    依據(jù)OSTU法給出類(lèi)間方差的定義,見(jiàn)式(8):

    σB=w0(u0-ut)2+w1(u1-ut)2

    (8)

    通過(guò)運(yùn)算選取出合適的閾值t,使得類(lèi)間方差σB的數(shù)值最大。在OSTU閾值分割法中分割閾值t可以通過(guò)計(jì)算機(jī)本身的窮舉計(jì)算得到,對(duì)于不同的情況可以動(dòng)態(tài)的自適應(yīng)性選取閾值。此時(shí)按照閾值t,對(duì)圖3b進(jìn)行閾值分割。圖4為一次閾值分割后的圖像,從圖中可以看出,圖中河流和其他地物的分界線(xiàn)更加清晰明確,背景也得到了進(jìn)一步的簡(jiǎn)化。

    圖4 采用OSTU分割后的二值圖像

    2 河流輪廓識(shí)別和提取

    經(jīng)過(guò)Kirsch算子提取地物邊緣和一次閾值分割后,原始的灰度圖背景噪聲得到了簡(jiǎn)化。但是從圖4可以看出,一些零星的背景噪聲仍無(wú)法完好地去除。此時(shí)需要結(jié)合河流的基本形態(tài)特征對(duì)圖像進(jìn)行一些基本的降噪處理。

    2.1 形態(tài)學(xué)處理

    觀(guān)察圖中的零星噪聲,可以發(fā)現(xiàn)不屬于河流的噪聲為一個(gè)個(gè)非連通的孤立區(qū)域??梢酝ㄟ^(guò)一定原則從圖中將這些小對(duì)象去除。去除原則為:當(dāng)孤立區(qū)域的面積S小于設(shè)定的閾值Sc時(shí),將其去除。參考的Sc取值見(jiàn)式(9):

    Sc=k·l1·l2

    (9)

    式中,l1、l2為圖像的長(zhǎng)度、寬度;k為面積系數(shù),取值盡量小,本文推薦k=0.000 5。當(dāng)k設(shè)置為0.000 5時(shí),由于此值足夠小,河流在圖像中的面積必然大于這一數(shù)值,因此不會(huì)將河流區(qū)域誤刪除。由于Sc隨著圖像面積進(jìn)行變化,可以保證在不同分辨率圖像中較小的噪聲被去除。故對(duì)于不同分辨率的圖像,這一閾值具有普適性。

    圖5為刪除圖中小型非聯(lián)通區(qū)域后的圖像。從圖中可以看出此時(shí)河流的輪廓已經(jīng)基本上清晰。但是遺憾的是,仍存在著一些較大的非聯(lián)通區(qū)域無(wú)法去除,如圖5中的紅色虛線(xiàn)所示。此時(shí)若增大k值雖然可以去除這些區(qū)域,但是可能去除一些細(xì)小的河流分支,使得識(shí)別出的河流輪廓不完整。

    圖5 刪除小型非聯(lián)通區(qū)域

    此時(shí)可以利用河流的形態(tài)特征:河流長(zhǎng)度較長(zhǎng),且順?biāo)鞣较虻拈L(zhǎng)度遠(yuǎn)大于河流寬度。此時(shí)可以檢查余下的區(qū)域,對(duì)于每一個(gè)非聯(lián)通區(qū)域,當(dāng)面積較小,且2個(gè)方向的長(zhǎng)度d1和d2的比值較小,則可以將該區(qū)域去除。設(shè)定如下2個(gè)去除條件,2個(gè)條件同時(shí)滿(mǎn)足時(shí),則將該區(qū)域去除。

    (10)

    式中,k2為面積系數(shù),其值相對(duì)較大,本文推薦k2=0.01;n為河流長(zhǎng)寬比的限定值,本文推薦n=4。當(dāng)閾值k2設(shè)置為0.01時(shí),這一數(shù)值也相對(duì)較小,一般情況下河流在圖像中的面積會(huì)大于這一數(shù)值。當(dāng)然若是河流被截?cái)嗪竺娣e較小,長(zhǎng)寬比的限制條件將會(huì)發(fā)揮作用,保證部分河流不被誤刪除。因此該閾值對(duì)于不同情況也具有一定的普適性。在后文中2.3節(jié)的算例,也證明了這點(diǎn)。

    圖6為刪除無(wú)河流特征的連通域后的圖像,從圖中可以看出一些較大的噪聲也已經(jīng)基本上被刪除,僅留下河流附近區(qū)域的圖像。在此過(guò)程中河流的一些細(xì)小的部分也可能被誤刪除。圖6所示,圖中圓形虛線(xiàn)區(qū)域所指的區(qū)域原本是聯(lián)通的,卻在形態(tài)學(xué)處理時(shí)被誤刪除。

    圖6 刪除無(wú)河流特征的較大型連通域

    這種誤刪除的錯(cuò)誤將會(huì)導(dǎo)致識(shí)別出的河流區(qū)域出現(xiàn)殘缺。此時(shí)可以利用基本的形態(tài)學(xué)處理方式將一些斷開(kāi)的河流區(qū)域進(jìn)行連接。具體方法為:首先對(duì)圖像進(jìn)行膨脹操作,之后對(duì)距離相近的非聯(lián)通域進(jìn)行連接,具體的實(shí)施過(guò)程可參考王康等[12]關(guān)于裂紋識(shí)別的處理。圖7為對(duì)圖像進(jìn)行膨脹和連接操作之后的圖像,從圖中可以看出紅色虛線(xiàn)指示區(qū)域河流已連接。

    圖7 膨脹和連接的形態(tài)學(xué)運(yùn)算

    2.2 二次閾值分割

    對(duì)河流圖像進(jìn)行基本形態(tài)學(xué)處理后,雖然河流周?chē)脑肼曇鸦颈磺宄?但是河流輪廓和形貌的提取依然存在問(wèn)題。例如圖7中紅色矩形框是河流與田地的混雜的區(qū)域,但是卻完全被劃分為河流。因此仍需要進(jìn)行二次閾值分割,以獲取更為準(zhǔn)確的河流輪廓形貌。

    圖7中的白色區(qū)域已基本上將河流包含在內(nèi),因此可以?xún)H對(duì)圖7中的白色區(qū)域進(jìn)行第二次閾值分割,不考慮圖7中的黑色區(qū)域。為此將原始的灰度圖中圖像映射至白色區(qū)域中,得到圖8,以作為二次閾值分割的原始圖形。

    圖8 第二次閾值分割的區(qū)域

    采用1.2節(jié)的OSTU中的閾值分割算法,對(duì)剩余圖像進(jìn)行二次閾值分割,得到圖9a閾值分割結(jié)果。此時(shí)圖形會(huì)存在一些毛刺和細(xì)小的噪聲,采用腐蝕膨脹等形態(tài)學(xué)處理手段,得到相對(duì)光滑的河流輪廓和形貌,見(jiàn)圖9b。

    a)原始灰度圖

    2.3 適用性驗(yàn)證

    若提取河道輪廓時(shí),未去除背景噪聲的影響,直接采用閾值分割。受到復(fù)雜背景的影響,提取內(nèi)容不僅包含河道,還包含一些背景中的地物,效果極差,見(jiàn)圖10。

    圖10 直接閾值剖分的提取效果

    為驗(yàn)證本文方法對(duì)于不同情況下的適用性,更換不同地形區(qū)域的河道遙感圖像,查看提取效果。分別選取平原地區(qū)、山區(qū)和丘陵地區(qū)河流圖像各20張,總計(jì)60張,以驗(yàn)證本文算法的適用性。圖11為不同地區(qū)典型的河道遙感圖像及提取效果,可以看出對(duì)于不同地形下的河流遙感圖像,河道輪廓均被準(zhǔn)確的提取出來(lái)。

    a) 平原復(fù)雜河道灰度

    為更為客觀(guān)地評(píng)價(jià)河道提取的準(zhǔn)確性,采用像素精度指標(biāo)pixel accuracy (PA)和均像素精度指標(biāo)mean pixel accuracy (MPA)評(píng)價(jià)該方法識(shí)別的準(zhǔn)確性。其中像素精度含義為:對(duì)于每一類(lèi)像素,正確分類(lèi)的像素個(gè)數(shù)之和與該類(lèi)所有像素的個(gè)數(shù)之和的比值。均像素精度含義為:每一類(lèi)像素的精度平均值,即先求出每一類(lèi)像素的PA,再取平均值。表1分別為演示算例和平原地區(qū)、山區(qū)和丘陵地區(qū)每類(lèi)的PA值。從表1中可以看出在不同地形下,背景像素的識(shí)別精度達(dá)到了98%以上,河流像素的識(shí)別精度也均超過(guò)了90%。平均像素精度均超過(guò)了95%,誤差在5%以?xún)?nèi)。

    表1 河道提取的像素精度值

    當(dāng)然對(duì)于遙感信息中河流信息提取精度的評(píng)價(jià),還可以采用生產(chǎn)者精度(Producer’s accuracy)和使用者精度(User’s Accuracy)來(lái)表征。對(duì)于相同的算例,表2中列出了生產(chǎn)者精度和使用者精度。其中河流信息中的使用者精度與像素精度指標(biāo)的定義相同,因此數(shù)值相同。從表2中可以看出河流像素識(shí)別的生產(chǎn)者精度和使用者精度相近,均超過(guò)90%。

    表2 河流信息提取的生產(chǎn)者精度和使用者精度

    2種評(píng)價(jià)指標(biāo)的結(jié)果均表明:本文方法對(duì)于河道提取結(jié)果較好,識(shí)別精度較高,能夠滿(mǎn)足高精度河道提取的需要。

    對(duì)于較大范圍的河流,本文給出的提取方案仍具有適用性。圖12a為2018年拍攝的穿過(guò)重慶主城區(qū)的長(zhǎng)江和嘉陵江圖像,圖12b提取出的河道輪廓基本上與河流形貌吻合。當(dāng)然由于河流上橋梁較多,而Kirsch算子提取的邊緣輪廓較粗,這將會(huì)導(dǎo)致提取出的河道輪廓存在斷點(diǎn),見(jiàn)圖12b中的虛線(xiàn)。這需要在之后的研究中加以改善。

    a)重慶主城區(qū)嘉陵江河道灰度

    圖13a為1967年采用鎖眼衛(wèi)星拍攝的相同地點(diǎn)處的遙感圖像,此時(shí)衛(wèi)星圖像的分辨率低(分辨率大于10 m)。采用本文的算法得到的河道輪廓見(jiàn)圖13b,可以看出此時(shí)提取出的河道輪廓與實(shí)際情況相吻合。因此本文提出的算法對(duì)中等分辨率的圖像也具有適用性。

    a)重慶主城區(qū)嘉陵江河道灰度

    綜合2.3節(jié)中的算例,可以發(fā)現(xiàn)本文提出的輪廓提取算法對(duì)于不同地區(qū),不同范圍以及不同分辨率的衛(wèi)星遙感圖像,均具有較好的適用性,能夠較好地區(qū)分出河流與山丘、城市、農(nóng)田等其余地物,從而準(zhǔn)確地提取出河流的輪廓。由于多次閾值分割均采用OSTU法,而OSTU法閾值本身具有的動(dòng)態(tài)自適應(yīng)特性,一定程度上也增加了方法的普適性。

    3 結(jié)論

    主要解決現(xiàn)階段遙感圖像分辨率提高,導(dǎo)致在復(fù)雜背景噪聲下河道輪廓難以準(zhǔn)確識(shí)別和提取的問(wèn)題。為此提出了一種基于Kirsch算子、多級(jí)閾值分割和形態(tài)學(xué)處理的河道提取方法。采用Kirsch算子對(duì)高清灰度圖進(jìn)行預(yù)加工,提取出主要地物的輪廓,之后采用OSTU的閾值分割方法對(duì)預(yù)處理后的圖像進(jìn)行預(yù)分割,去除了部分背景噪聲。再利用河流面積和長(zhǎng)寬比等形態(tài)學(xué)特征,去除大部分的復(fù)雜背景噪聲。最后對(duì)圖像局部區(qū)域采用二次閾值分割,準(zhǔn)確地提取出河道的信息。

    該方法能夠去除遙感圖像中的復(fù)雜背景噪聲,在不同地形條件下實(shí)現(xiàn)對(duì)河道輪廓的準(zhǔn)確識(shí)別和提取。通過(guò)生產(chǎn)者精度等客觀(guān)的評(píng)價(jià)指標(biāo),對(duì)方法的適用性和準(zhǔn)確性進(jìn)行了驗(yàn)證。結(jié)果表明在不同地形條件下,河流區(qū)域識(shí)別的準(zhǔn)確度均超過(guò)90%,圖像整體的識(shí)別精度超過(guò)95%。因此該方法能夠滿(mǎn)足高精度河道提取的需要,為智能化的提取河道信息提供了參考。

    當(dāng)然由于Kirsch算子提取的邊緣輪廓較粗,當(dāng)河流上存在橋梁時(shí),會(huì)導(dǎo)致提取的河道存在斷點(diǎn),需要在之后的研究中對(duì)此加以改善。由于本文主要采用全波段的衛(wèi)星遙感圖像,對(duì)于單波段或者多波段的遙感圖像未作深入的分析和研究。在下一階段的研究中,可以考慮不同波段的選擇對(duì)于提取效果的影響,以尋找更為合適的遙感波段,進(jìn)一步提高河道提取的準(zhǔn)確性。

    猜你喜歡
    輪廓算子灰度
    采用改進(jìn)導(dǎo)重法的拓?fù)浣Y(jié)構(gòu)灰度單元過(guò)濾技術(shù)
    擬微分算子在Hp(ω)上的有界性
    基于灰度拉伸的圖像水位識(shí)別方法研究
    OPENCV輪廓識(shí)別研究與實(shí)踐
    各向異性次Laplace算子和擬p-次Laplace算子的Picone恒等式及其應(yīng)用
    基于實(shí)時(shí)輪廓誤差估算的數(shù)控系統(tǒng)輪廓控制
    一類(lèi)Markov模算子半群與相應(yīng)的算子值Dirichlet型刻畫(huà)
    基于最大加權(quán)投影求解的彩色圖像灰度化對(duì)比度保留算法
    Roper-Suffridge延拓算子與Loewner鏈
    基于灰度線(xiàn)性建模的亞像素圖像抖動(dòng)量計(jì)算
    精品一区在线观看国产| a级片在线免费高清观看视频| 如何舔出高潮| 麻豆乱淫一区二区| 王馨瑶露胸无遮挡在线观看| 欧美日韩视频高清一区二区三区二| 午夜免费鲁丝| 久久99热6这里只有精品| 乱码一卡2卡4卡精品| 亚洲,欧美,日韩| 精品少妇黑人巨大在线播放| 视频在线观看一区二区三区| 欧美精品国产亚洲| 大香蕉97超碰在线| 高清在线视频一区二区三区| 久久久国产精品麻豆| 中文天堂在线官网| xxx大片免费视频| 91精品伊人久久大香线蕉| 51国产日韩欧美| 新久久久久国产一级毛片| 久久这里有精品视频免费| 涩涩av久久男人的天堂| 日产精品乱码卡一卡2卡三| 久久女婷五月综合色啪小说| 亚洲国产看品久久| 国产深夜福利视频在线观看| 女人被躁到高潮嗷嗷叫费观| 国产精品久久久久久久久免| 免费黄色在线免费观看| 国产精品久久久久久精品电影小说| 美女国产高潮福利片在线看| 国产黄色视频一区二区在线观看| 下体分泌物呈黄色| 日日摸夜夜添夜夜爱| 午夜福利影视在线免费观看| 丝袜喷水一区| 亚洲欧美精品自产自拍| 亚洲精品色激情综合| 嫩草影院入口| 国产精品一二三区在线看| 久久狼人影院| 亚洲一码二码三码区别大吗| 国产成人精品婷婷| 欧美人与性动交α欧美软件 | 午夜福利在线观看免费完整高清在| 制服诱惑二区| 天天躁夜夜躁狠狠久久av| 久久人人97超碰香蕉20202| 咕卡用的链子| 性色avwww在线观看| 久久午夜福利片| 亚洲欧美日韩卡通动漫| 国产精品人妻久久久影院| 国产麻豆69| 成人国产麻豆网| 日韩伦理黄色片| 少妇人妻 视频| 国产国语露脸激情在线看| 日本午夜av视频| 菩萨蛮人人尽说江南好唐韦庄| 69精品国产乱码久久久| 人人妻人人澡人人看| 日韩伦理黄色片| 国产欧美日韩综合在线一区二区| 桃花免费在线播放| 国产精品三级大全| 99九九在线精品视频| av线在线观看网站| 老司机影院毛片| 亚洲伊人久久精品综合| www日本在线高清视频| 黄色 视频免费看| 两个人免费观看高清视频| 天堂中文最新版在线下载| xxx大片免费视频| 美女福利国产在线| 精品99又大又爽又粗少妇毛片| 亚洲精品美女久久久久99蜜臀 | 色婷婷久久久亚洲欧美| 久久婷婷青草| 亚洲国产精品一区二区三区在线| av又黄又爽大尺度在线免费看| 国产免费一区二区三区四区乱码| av线在线观看网站| 精品国产露脸久久av麻豆| 搡女人真爽免费视频火全软件| 天堂8中文在线网| 国产黄色视频一区二区在线观看| 国产又色又爽无遮挡免| 如日韩欧美国产精品一区二区三区| 中国国产av一级| 国产 精品1| 男人添女人高潮全过程视频| freevideosex欧美| 精品第一国产精品| 亚洲国产成人一精品久久久| 最近中文字幕2019免费版| 极品人妻少妇av视频| 久久 成人 亚洲| 欧美成人精品欧美一级黄| 久久久久久久久久成人| 国产熟女午夜一区二区三区| 美女福利国产在线| 午夜影院在线不卡| 男女边摸边吃奶| 一级毛片电影观看| 少妇 在线观看| 人妻 亚洲 视频| 成人毛片a级毛片在线播放| 大片电影免费在线观看免费| 狂野欧美激情性bbbbbb| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 视频在线观看一区二区三区| 色视频在线一区二区三区| 国产精品久久久久久精品古装| 亚洲精品美女久久av网站| 69精品国产乱码久久久| 夜夜骑夜夜射夜夜干| 久久久国产欧美日韩av| 这个男人来自地球电影免费观看 | 国产男人的电影天堂91| 免费黄色在线免费观看| 大香蕉久久网| 99久久中文字幕三级久久日本| 99久久精品国产国产毛片| 欧美精品高潮呻吟av久久| 久久人妻熟女aⅴ| 又粗又硬又长又爽又黄的视频| 国产精品久久久久久精品电影小说| 人妻 亚洲 视频| 丝袜人妻中文字幕| 最近中文字幕高清免费大全6| 亚洲国产欧美日韩在线播放| 国产精品一区二区在线观看99| 国产成人91sexporn| 黑人欧美特级aaaaaa片| 青春草亚洲视频在线观看| 亚洲成色77777| 久久精品国产鲁丝片午夜精品| a级片在线免费高清观看视频| 免费看不卡的av| av有码第一页| 国产免费一级a男人的天堂| 中文欧美无线码| 美女主播在线视频| 精品久久蜜臀av无| 欧美国产精品一级二级三级| 久久精品久久精品一区二区三区| 亚洲欧美清纯卡通| 欧美精品亚洲一区二区| 久久人人97超碰香蕉20202| 美女中出高潮动态图| 色5月婷婷丁香| 亚洲中文av在线| 天天操日日干夜夜撸| 亚洲伊人久久精品综合| 久久青草综合色| 亚洲精品国产av蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕最新亚洲高清| 欧美成人精品欧美一级黄| 亚洲人成77777在线视频| 日韩中文字幕视频在线看片| 久久99热6这里只有精品| 亚洲中文av在线| 爱豆传媒免费全集在线观看| 日韩一区二区三区影片| 久久国产精品大桥未久av| 少妇人妻 视频| 最近最新中文字幕免费大全7| 在线天堂最新版资源| 亚洲熟女精品中文字幕| 丰满乱子伦码专区| 五月天丁香电影| 九色亚洲精品在线播放| 久久久国产欧美日韩av| 久久精品国产亚洲av涩爱| 女人久久www免费人成看片| 午夜福利网站1000一区二区三区| 男的添女的下面高潮视频| 国产国语露脸激情在线看| 国产高清国产精品国产三级| 国语对白做爰xxxⅹ性视频网站| 最近的中文字幕免费完整| 国产高清三级在线| 日韩熟女老妇一区二区性免费视频| 一级片免费观看大全| 最新中文字幕久久久久| 中文字幕最新亚洲高清| 春色校园在线视频观看| 久久精品国产a三级三级三级| av网站免费在线观看视频| 成人18禁高潮啪啪吃奶动态图| 久久久久久人妻| 中文天堂在线官网| 啦啦啦视频在线资源免费观看| 亚洲精品国产av蜜桃| 午夜日本视频在线| 亚洲国产最新在线播放| 亚洲在久久综合| av又黄又爽大尺度在线免费看| 麻豆精品久久久久久蜜桃| 女人久久www免费人成看片| 亚洲国产欧美在线一区| 国产精品久久久久久久久免| 1024视频免费在线观看| 国产免费一区二区三区四区乱码| 精品亚洲成国产av| 免费大片黄手机在线观看| 中文欧美无线码| 亚洲欧美日韩卡通动漫| 国产xxxxx性猛交| a级毛片在线看网站| 美女xxoo啪啪120秒动态图| 婷婷成人精品国产| 中文字幕另类日韩欧美亚洲嫩草| 免费大片18禁| 少妇精品久久久久久久| 亚洲精品美女久久久久99蜜臀 | 日韩不卡一区二区三区视频在线| 久久国内精品自在自线图片| 深夜精品福利| 亚洲国产av新网站| 97人妻天天添夜夜摸| 日本欧美视频一区| 街头女战士在线观看网站| 欧美日韩视频高清一区二区三区二| 亚洲国产毛片av蜜桃av| 久久99蜜桃精品久久| 亚洲高清免费不卡视频| 日韩在线高清观看一区二区三区| 亚洲综合色网址| 18禁动态无遮挡网站| 美女大奶头黄色视频| 91精品伊人久久大香线蕉| 国产又色又爽无遮挡免| 99久久中文字幕三级久久日本| 97人妻天天添夜夜摸| 国产成人精品在线电影| 国产视频首页在线观看| videos熟女内射| 日本av免费视频播放| 黄色毛片三级朝国网站| 乱码一卡2卡4卡精品| av在线观看视频网站免费| 一区二区三区乱码不卡18| 中国美白少妇内射xxxbb| 久久99蜜桃精品久久| 咕卡用的链子| 久久精品aⅴ一区二区三区四区 | 高清黄色对白视频在线免费看| 91aial.com中文字幕在线观看| 欧美少妇被猛烈插入视频| 国产精品久久久久成人av| 久久人人97超碰香蕉20202| 飞空精品影院首页| 久久久a久久爽久久v久久| 精品亚洲成国产av| 亚洲成色77777| 免费在线观看黄色视频的| 一区二区三区四区激情视频| 久久久久国产精品人妻一区二区| 亚洲精品一区蜜桃| 九草在线视频观看| 亚洲欧美成人综合另类久久久| √禁漫天堂资源中文www| 久久国内精品自在自线图片| 亚洲三级黄色毛片| 国产精品一区二区在线观看99| 十分钟在线观看高清视频www| 好男人视频免费观看在线| tube8黄色片| 国产午夜精品一二区理论片| 夫妻性生交免费视频一级片| 看非洲黑人一级黄片| 国产日韩一区二区三区精品不卡| 丰满少妇做爰视频| 精品亚洲乱码少妇综合久久| 欧美日本中文国产一区发布| 免费观看a级毛片全部| 极品少妇高潮喷水抽搐| 巨乳人妻的诱惑在线观看| 国产淫语在线视频| 免费看av在线观看网站| 欧美日本中文国产一区发布| 男的添女的下面高潮视频| 有码 亚洲区| 丰满少妇做爰视频| 中文字幕制服av| 插逼视频在线观看| 国产一区二区在线观看av| 国产精品一区二区在线不卡| a级毛色黄片| 午夜福利视频在线观看免费| 亚洲成色77777| 国产精品秋霞免费鲁丝片| 夜夜骑夜夜射夜夜干| 久久久国产欧美日韩av| 全区人妻精品视频| 亚洲成国产人片在线观看| 国语对白做爰xxxⅹ性视频网站| 一级黄片播放器| 99国产精品免费福利视频| 国产国语露脸激情在线看| 国产精品嫩草影院av在线观看| 久久久欧美国产精品| 日韩制服骚丝袜av| 女人被躁到高潮嗷嗷叫费观| 欧美日韩成人在线一区二区| a级毛片黄视频| 日韩精品免费视频一区二区三区 | 亚洲成人一二三区av| 国产成人精品一,二区| 亚洲人成77777在线视频| 成人毛片60女人毛片免费| 搡老乐熟女国产| 久久人人爽人人片av| 九九爱精品视频在线观看| 国产一区二区三区综合在线观看 | 99九九在线精品视频| 婷婷色av中文字幕| 欧美激情 高清一区二区三区| 精品福利永久在线观看| 免费黄频网站在线观看国产| 亚洲欧美一区二区三区黑人 | 亚洲av国产av综合av卡| 国产欧美日韩综合在线一区二区| 91在线精品国自产拍蜜月| 涩涩av久久男人的天堂| 美女主播在线视频| 欧美成人精品欧美一级黄| 91午夜精品亚洲一区二区三区| 一级毛片我不卡| 国产一区二区激情短视频 | 亚洲精品自拍成人| 毛片一级片免费看久久久久| 亚洲性久久影院| 亚洲色图 男人天堂 中文字幕 | 男人添女人高潮全过程视频| 精品一区二区三区四区五区乱码 | 97在线人人人人妻| 99香蕉大伊视频| 最近的中文字幕免费完整| 99视频精品全部免费 在线| 亚洲av电影在线观看一区二区三区| videos熟女内射| 亚洲国产精品专区欧美| 久久久精品区二区三区| 国产女主播在线喷水免费视频网站| 国内精品宾馆在线| 波野结衣二区三区在线| 亚洲精品久久午夜乱码| 少妇高潮的动态图| 搡老乐熟女国产| 中文字幕制服av| 亚洲成人av在线免费| 亚洲欧洲国产日韩| 在线观看免费高清a一片| 欧美成人午夜免费资源| 国产精品一区二区在线观看99| 女人精品久久久久毛片| 又黄又粗又硬又大视频| 免费黄频网站在线观看国产| 制服丝袜香蕉在线| 国产成人精品婷婷| 韩国av在线不卡| 在线精品无人区一区二区三| 最新中文字幕久久久久| 丰满少妇做爰视频| 亚洲国产毛片av蜜桃av| 精品人妻偷拍中文字幕| 精品一区在线观看国产| 婷婷色av中文字幕| 久久久久精品人妻al黑| 极品人妻少妇av视频| a 毛片基地| 国产探花极品一区二区| 精品一区在线观看国产| 在线观看免费视频网站a站| 亚洲综合色惰| 内地一区二区视频在线| 26uuu在线亚洲综合色| 黄网站色视频无遮挡免费观看| 欧美xxⅹ黑人| 精品一区二区三区视频在线| 大陆偷拍与自拍| 久久97久久精品| 免费看不卡的av| 国产精品一区二区在线不卡| 久久精品aⅴ一区二区三区四区 | 国产一区二区在线观看日韩| 激情视频va一区二区三区| 制服诱惑二区| 激情视频va一区二区三区| 久久久亚洲精品成人影院| av国产精品久久久久影院| 男人添女人高潮全过程视频| 最近最新中文字幕大全免费视频 | 国产色爽女视频免费观看| 国产精品嫩草影院av在线观看| 少妇的逼水好多| www.av在线官网国产| 欧美bdsm另类| 久久人人爽人人片av| av.在线天堂| 日本爱情动作片www.在线观看| 久久毛片免费看一区二区三区| 亚洲色图 男人天堂 中文字幕 | 亚洲欧美成人精品一区二区| av一本久久久久| 在线观看一区二区三区激情| 美女国产视频在线观看| 中文字幕亚洲精品专区| 午夜日本视频在线| www日本在线高清视频| 制服丝袜香蕉在线| 国产亚洲精品久久久com| 捣出白浆h1v1| 国产精品欧美亚洲77777| 国产又爽黄色视频| 大码成人一级视频| 国产日韩欧美在线精品| 日韩三级伦理在线观看| 久久久a久久爽久久v久久| 日韩精品有码人妻一区| 美女视频免费永久观看网站| 大香蕉久久成人网| 一本大道久久a久久精品| 香蕉丝袜av| 丰满少妇做爰视频| 国产高清不卡午夜福利| 国产成人a∨麻豆精品| 久久亚洲国产成人精品v| 国产精品欧美亚洲77777| 亚洲精品第二区| 亚洲精品自拍成人| 在线免费观看不下载黄p国产| 国产高清不卡午夜福利| 青青草视频在线视频观看| 日韩成人伦理影院| 精品午夜福利在线看| 少妇精品久久久久久久| 国产 一区精品| 麻豆乱淫一区二区| 欧美xxⅹ黑人| 王馨瑶露胸无遮挡在线观看| 亚洲av电影在线进入| 激情五月婷婷亚洲| 内地一区二区视频在线| 久久国内精品自在自线图片| 国产av一区二区精品久久| 国产一区二区在线观看日韩| 秋霞在线观看毛片| 国产黄色免费在线视频| 一区二区三区乱码不卡18| 午夜福利乱码中文字幕| 亚洲伊人久久精品综合| 亚洲精品中文字幕在线视频| 香蕉国产在线看| 男人爽女人下面视频在线观看| 国产在线免费精品| 男女高潮啪啪啪动态图| 国产色婷婷99| 欧美日韩精品成人综合77777| 国产黄色视频一区二区在线观看| 久久久国产精品麻豆| 男女边吃奶边做爰视频| 久久精品国产a三级三级三级| 精品人妻偷拍中文字幕| www.熟女人妻精品国产 | 狂野欧美激情性bbbbbb| 国产成人精品久久久久久| 成人18禁高潮啪啪吃奶动态图| 丁香六月天网| 一区二区三区精品91| 一本久久精品| 久久精品国产a三级三级三级| 亚洲一级一片aⅴ在线观看| 91成人精品电影| 国产国语露脸激情在线看| 成年av动漫网址| xxx大片免费视频| 日韩精品免费视频一区二区三区 | 日日爽夜夜爽网站| 99视频精品全部免费 在线| 亚洲精品国产av成人精品| 成年美女黄网站色视频大全免费| 99热全是精品| av播播在线观看一区| 国产 精品1| 日韩中字成人| 九九爱精品视频在线观看| 少妇的逼水好多| 制服人妻中文乱码| 亚洲av电影在线进入| 性色avwww在线观看| 你懂的网址亚洲精品在线观看| 在线免费观看不下载黄p国产| 欧美亚洲 丝袜 人妻 在线| 国产精品久久久久久精品古装| av片东京热男人的天堂| 丝袜喷水一区| 一边亲一边摸免费视频| 久久精品国产亚洲av涩爱| 亚洲精品av麻豆狂野| 免费av不卡在线播放| a级毛片黄视频| 亚洲精品自拍成人| 国产亚洲一区二区精品| 日日撸夜夜添| 午夜影院在线不卡| 爱豆传媒免费全集在线观看| 国产成人一区二区在线| 91精品国产国语对白视频| 人成视频在线观看免费观看| 日韩 亚洲 欧美在线| 蜜臀久久99精品久久宅男| 国产深夜福利视频在线观看| 亚洲精品乱久久久久久| 日本91视频免费播放| 亚洲欧洲精品一区二区精品久久久 | 久久久久久人妻| av片东京热男人的天堂| 亚洲三级黄色毛片| 中文字幕av电影在线播放| 色94色欧美一区二区| 好男人视频免费观看在线| 男女边摸边吃奶| 欧美3d第一页| 大码成人一级视频| 久久国产精品男人的天堂亚洲 | av女优亚洲男人天堂| 免费观看a级毛片全部| 三级国产精品片| 91午夜精品亚洲一区二区三区| 国产精品蜜桃在线观看| av免费在线看不卡| 在线观看人妻少妇| 青青草视频在线视频观看| 五月开心婷婷网| 亚洲精品日本国产第一区| 一区在线观看完整版| 极品少妇高潮喷水抽搐| 两个人看的免费小视频| 一级毛片黄色毛片免费观看视频| 丝袜脚勾引网站| 看非洲黑人一级黄片| 亚洲丝袜综合中文字幕| 成人漫画全彩无遮挡| 少妇高潮的动态图| 亚洲精品av麻豆狂野| 最近中文字幕2019免费版| 国内精品宾馆在线| 男的添女的下面高潮视频| 亚洲精品日本国产第一区| 日韩伦理黄色片| 亚洲伊人色综图| 午夜激情久久久久久久| 日韩视频在线欧美| 在线观看免费高清a一片| 成人免费观看视频高清| 欧美日韩视频高清一区二区三区二| 成年人免费黄色播放视频| 成人18禁高潮啪啪吃奶动态图| 在线观看免费视频网站a站| 日日撸夜夜添| 黄色 视频免费看| 日韩伦理黄色片| 在线观看人妻少妇| 一区二区日韩欧美中文字幕 | 人妻 亚洲 视频| 国产免费一区二区三区四区乱码| 国产精品久久久久久久电影| 97人妻天天添夜夜摸| 日本色播在线视频| 亚洲精品久久成人aⅴ小说| 色视频在线一区二区三区| 在线看a的网站| 在线观看www视频免费| 亚洲av国产av综合av卡| 一个人免费看片子| 欧美另类一区| 你懂的网址亚洲精品在线观看| 一二三四在线观看免费中文在 | 亚洲av.av天堂| 国产毛片在线视频| 日韩欧美精品免费久久| 午夜av观看不卡| 国产亚洲av片在线观看秒播厂| 成人漫画全彩无遮挡| 欧美日韩综合久久久久久| 国产精品久久久久久精品古装| 桃花免费在线播放| 亚洲精品美女久久av网站| av不卡在线播放| 少妇被粗大猛烈的视频| 午夜免费男女啪啪视频观看| 男男h啪啪无遮挡| 国产精品国产三级专区第一集| 免费黄网站久久成人精品| 精品第一国产精品| 色网站视频免费| 久久久精品区二区三区| 十分钟在线观看高清视频www| 自线自在国产av| 久久人人爽av亚洲精品天堂| 精品国产一区二区三区久久久樱花| 美女视频免费永久观看网站| 日韩中文字幕视频在线看片| 成人免费观看视频高清| 丝袜脚勾引网站| 日韩中字成人| a级片在线免费高清观看视频| 久久精品国产综合久久久 |