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

    基于高階正則與非光滑數(shù)據(jù)擬合項(xiàng)的圖像邊緣檢測(cè)模型①

    2020-01-15 06:44:50陳靜思王鵬彥
    關(guān)鍵詞:正則高階算子

    李 春,陳靜思,王鵬彥,李 健,羅 澤

    1(中國(guó)科學(xué)院 計(jì)算機(jī)網(wǎng)絡(luò)信息中心,北京 100190)

    2(中國(guó)科學(xué)院大學(xué),北京 100049)

    3(云南財(cái)經(jīng)大學(xué) 云南省經(jīng)濟(jì)社會(huì)大數(shù)據(jù)研究院,昆明 650221)

    4(四川臥龍國(guó)家級(jí)自然保護(hù)區(qū)管理局,臥龍 623006)

    1 圖像分割簡(jiǎn)介

    圖像分割和邊緣檢測(cè)是計(jì)算機(jī)視覺(jué),特別是圖像分析中一個(gè)基礎(chǔ)而重要的問(wèn)題(如:目標(biāo)識(shí)別和圖像解析),其目的是把一幅給定圖像根據(jù)其圖像特征(如:邊緣,顏色,紋理,運(yùn)動(dòng),角點(diǎn)特征等)劃分為不同區(qū)域[1].近年來(lái),隨著計(jì)算機(jī)技術(shù)發(fā)展,使得圖像分割成文現(xiàn)代科學(xué)領(lǐng)域的研究熱點(diǎn),圖像分割在醫(yī)學(xué)圖像分析,遙感,監(jiān)測(cè)等領(lǐng)域得到廣泛應(yīng)用[2–11].在諸如物體識(shí)別,分類和分割的應(yīng)用中尋找物體邊緣是非常重要的一步.因此,使用和設(shè)計(jì)何種邊緣檢測(cè)算法直接影響這些應(yīng)用性能.如:2019年Yuan等[12]提出了一個(gè)量子圖像邊緣檢測(cè)算法,將邊緣檢測(cè)算法靈活應(yīng)用于量子計(jì)算.圖像邊緣檢測(cè)在雷達(dá)圖像檢測(cè)中也有重要應(yīng)用,Ma等[13]提出通過(guò)稀疏表示進(jìn)行SAR圖像邊緣檢測(cè),作者提出了一種基于去噪算法的新型合成孔徑雷達(dá)(SAR)圖像檢測(cè)算法,該算法通過(guò)稀疏表達(dá)和一種新的形態(tài)學(xué)邊緣檢測(cè)器.首先,作者將Shearlet變換應(yīng)用于SAR圖像以獲得圖像的稀疏表示.然后,將帶方向的形態(tài)學(xué)邊緣檢測(cè)器應(yīng)用于Shearlet的方向子帶系數(shù),其通過(guò)迭代去噪處理來(lái)恢復(fù).2019年Sert等[14]提出了一個(gè)新的基于最大范數(shù)熵的中智學(xué)邊緣檢測(cè)方法.同時(shí),邊緣檢測(cè)算法在醫(yī)學(xué)圖像領(lǐng)域也有重要應(yīng)用,隨著醫(yī)學(xué)領(lǐng)域成像技術(shù)的快速發(fā)展,不同模態(tài)的醫(yī)學(xué)圖像具有不同的成像原理,反映了人體生理信息的不同重點(diǎn)和缺陷.邊緣檢測(cè)是醫(yī)學(xué)超聲圖像處理的關(guān)鍵步驟,檢測(cè)結(jié)果將直接影響醫(yī)生對(duì)疾病的診斷.圖像邊緣檢測(cè)可以看作邊緣點(diǎn)和非邊緣點(diǎn)的分類問(wèn)題.基于此,Song等于2019年[15]提出了基于改進(jìn)差分進(jìn)化算法和Prewitt算子的醫(yī)學(xué)圖像邊緣檢測(cè).Churchill等于2019年[16]提供了理論證據(jù)來(lái)解釋邊緣位置在CT影像重建中的潛在重要性.然后開(kāi)發(fā)一種實(shí)用的方案來(lái)利用這一理論.最后,作者通過(guò)實(shí)驗(yàn)證實(shí)了提出的使用邊緣遮蔽正則化可以提高CT重建的準(zhǔn)確性和重建速度.

    近年來(lái),基于變分法和偏微分方程(PDES)的圖像分割方法得到了廣泛研究[17–24].早期圖像分割模型代表為活動(dòng)輪廓模型,活動(dòng)輪廓模型根據(jù)圖像邊緣信息和區(qū)域信息對(duì)給定圖像進(jìn)行分割.基于邊緣信息和區(qū)域信息考慮,Kass等于1988年提出基于邊緣信息的能量最小化方法,即,Snake模型[25].其后,受到Snake活動(dòng)輪廓模型啟發(fā),Caselles等提出了綜合活動(dòng)輪廓長(zhǎng)度信息,提出了幾何活動(dòng)輪廓模型[26],與此同時(shí),Mumford和Shah提出了基于區(qū)域信息圖像分割變分框架,即,著名的Munmford-Shah (MS)模型[27].但是,由于MS模型能量泛函的非凸性,導(dǎo)致在當(dāng)時(shí)對(duì)該模型直接求解十分困難.所以,Chan和Vese為了克服MS模型所存在的弊端,并整合了變分法和水平集方法[17],提出了Chan-Vese (CV)模型[28].

    由于CV模型應(yīng)用廣泛,所以在過(guò)去幾十年里,該模型一直受到研究學(xué)者青睞.同時(shí),各研究學(xué)者開(kāi)發(fā)各種對(duì)該模型求解的有效方法,例如:為克服初始化輪廓敏感度問(wèn)題,Li等[29]修改CV模型,提出了利用懲罰函數(shù)作為約束且無(wú)需初始化過(guò)程的圖像分割模型.2014年Duan等[30]綜合了變量分離法,對(duì)偶方法,Bregman迭代和增廣的拉格朗日方法提出了更快更有效更的算法.這些算法在高計(jì)算精度和無(wú)需初始化方面取得重要成果.

    隨著計(jì)算機(jī)技術(shù)的發(fā)展,深度學(xué)習(xí)在圖像分割中取得重要進(jìn)展,最顯著的進(jìn)展是深度卷積神經(jīng)網(wǎng)絡(luò)(CNNS)不斷被應(yīng)用到場(chǎng)景解析中.典型的網(wǎng)絡(luò)包括:DeepLab (V1-V3)[31,32],Refinenet[33],PSP-Net[34]和DANet[35],然而,這些模型的缺點(diǎn)是計(jì)算代價(jià)相當(dāng)大,需要學(xué)習(xí)的參數(shù)特別多,為了降低這些模型在場(chǎng)景解析中的計(jì)算復(fù)雜度,過(guò)去幾年里,很多學(xué)者在降低計(jì)算復(fù)雜度方面做了大量工作.例如:ENet[36],ESPNet[37]大量減少了需要學(xué)習(xí)的參數(shù).

    但是,上述所提的變分優(yōu)化模型中,基本上都是假設(shè)圖像所含噪聲為高斯噪聲,幾乎是用L2范數(shù)對(duì)高斯噪聲加以擬合.然而在現(xiàn)實(shí)生活中,情況相對(duì)較為復(fù)雜,往往圖像在成像時(shí)還有來(lái)自外加噪聲的污染,幸運(yùn)的是,外加噪聲可以用L1范數(shù)得以很好擬合[38,39].利用L1范數(shù)作為數(shù)據(jù)擬合項(xiàng)作為分段常數(shù)圖像分割模型,并利用二階正則函數(shù)(TV2)對(duì)目標(biāo)函數(shù)加以懲罰,使得所提出模型能更好地逼近目標(biāo)函數(shù)的同時(shí)還能能更好地分割低對(duì)比度和含有外加噪聲的圖像.再者,本文通過(guò)變量分離方法把一復(fù)雜非凸問(wèn)題轉(zhuǎn)換為若干個(gè)簡(jiǎn)單凸子問(wèn)題進(jìn)行求解,從而巧妙地處理了不可導(dǎo)項(xiàng)和正則項(xiàng).

    2 預(yù)備知識(shí)

    在本節(jié)中,因?yàn)锳lternating Direction Method of Multipliers (ADMM)或Split Brgman迭代[40]方法是求解線性約束凸二次規(guī)劃問(wèn)題的有效方法,所以在此對(duì)該方法做如下簡(jiǎn)單介紹.考慮如下線性約束凸優(yōu)化問(wèn)題:

    其中,通過(guò)上述目標(biāo)函數(shù),可得如下增廣拉格朗日泛函:

    當(dāng)上述系統(tǒng)方程滿足某些收斂條件時(shí),ADMM方法在(2)中固定v更新u,固定u更 新v.因此,ADMM算法可以總結(jié)如算法1.

    算法1.ADMM算法γ>0u0v0 b0=0初始化:設(shè)k=0,選擇 ,,and ;

    for do(1)解問(wèn)題argmin k=1,2,3,···,K u E1(u)+γ 2A1u+A2v?f?bk2 2;(2)求解問(wèn)題argmin u E2(v)+γ 2A1u+A2v?f?bk2 2;(3)拉格朗日乘子更新bk+1=bk?(A1uk+1+A2vk+1?f);end

    3 模型提出與數(shù)值求解

    在本節(jié)中將介紹所提出的高階正則圖像輪廓檢測(cè)算法,隨后用變量分離法和ADMM算法對(duì)該模型詳細(xì)優(yōu)化過(guò)程和部分收斂性分析加以描述.

    設(shè)f:[0,1]2→R?L1(?)為觀測(cè)到信號(hào)(圖像),且滿足如下數(shù)學(xué)表達(dá)式,f=u+n,其中u為未受噪聲污染信號(hào)(圖像),n為外加噪聲(IN).為了重建u,可以考慮如下優(yōu)化問(wèn)題:

    其中,?u為u的 分布導(dǎo)數(shù),∥ ?u∥1為著名的total variation(TV)[41]有界變差空間中的半范數(shù),λ>0為權(quán)參數(shù).若用L2-數(shù)據(jù)擬合項(xiàng) ∥f?u∥2替代∥f?u∥1數(shù)據(jù)擬合項(xiàng),從而可得經(jīng)典的Rudin-Osher-Fatemi (ROF)[41]模型.從而圖像分割模型數(shù)學(xué)表達(dá)式如下:

    其中,Γ 為u的 不連續(xù)集合,H1為一維的Hausdorff測(cè)度,即,Γ的長(zhǎng)度.

    假設(shè)區(qū)域 [0,1]2被劃分為如下的N個(gè)區(qū)域,?1,?2,···,?N且同時(shí)滿足∪i=1,···,N?i=[0,1]2,?i∩?j=?,ij,且對(duì)于所有的i具 有ci∈R.特別地,當(dāng)N=2可以得到除了數(shù)據(jù)擬合項(xiàng)之外,其他都類似于CV模型的模型.

    由于 ∥f?u∥1不 可微,所以可以通過(guò)迭代得出c1c2近似解.在圖像恢復(fù)中,盡管基于TV[41]正則已經(jīng)取,得了很好的效果,但是在圖像恢復(fù)過(guò)程中TV正則會(huì)引起階梯效應(yīng).而高階正則不僅有可以去除階梯效應(yīng)的功能,而且高階正則還有光滑區(qū)域,很好地保持物體邊緣的效果,在本模型的實(shí)驗(yàn)過(guò)程中,假設(shè)圖像都是帶有噪聲的,又因?yàn)楦唠A正則比低階正則更具有光滑性和非線性性,所以為了更高的提取邊緣信息,綜合Jung等[38]提出的模型,我們提出如下的高階正則圖像邊緣檢測(cè)模型:

    通過(guò)觀察上述目標(biāo)函數(shù)可以發(fā)現(xiàn),為了對(duì)該目標(biāo)泛函進(jìn)行求解,(c1,c2,u),即,兩步驟算法:第1步:先固定u求解c1,c2,第2步:固定c1,c2求解u.同時(shí)并對(duì)離散的散度算子做了如下定義,且散度算子div:(RM×N)2→RM×N具有如下的共軛性質(zhì).?div.u=p.?u,?u∈RN×M,p∈(RM×N)2.

    故此,離散散度算子可以作如下定義:對(duì)

    其中,Dx和Dy為如下的向前向后分算子.

    以此類推,離散二階散度算子定義如下:div2:(RN×M)4→RN×M.且具有如下共軛性質(zhì):div2q.u=q.?u,?u∈RM×N,q∈(RN×M)4.對(duì)于,q=(q11,q12,q21,q22)∈(RN×M)4,我們定義,div2q(i,j)=Dxxq11(i,j)+Dxyq12(i,j)+Dyyq21(i,j)+Dyxq22(i,j).

    在實(shí)際應(yīng)用中,學(xué)者經(jīng)常用算子分離法(split Bregman iteration和ADMM)對(duì)(6)進(jìn)行求解.為達(dá)到求解目的,本文先把(6)無(wú)約束優(yōu)化問(wèn)題轉(zhuǎn)化成一個(gè)有約束優(yōu)化問(wèn)題,從而利用Bregman iteration進(jìn)行求解,引入中間變量的好處在于:通過(guò)引入如下輔助變量,不僅把無(wú)約束優(yōu)化問(wèn)題變成一個(gè)有約束優(yōu)化問(wèn)題,引入輔助變量使用變量分離法把一個(gè)復(fù)雜問(wèn)題分解成若干簡(jiǎn)單子問(wèn)題,使得分別對(duì)各子問(wèn)題求解相對(duì)容易.

    其中,v=(vx,vy),w=(wxx,wyy,wxy,wyx).這里需要進(jìn)行說(shuō)明的是該符號(hào)不是對(duì)v,w分別求導(dǎo),而是表示坐標(biāo).顯而易見(jiàn),無(wú)約束優(yōu)化問(wèn)題(6)和有約束優(yōu)化問(wèn)題(7)等價(jià),從而可以利用ADMM或Split Bregman迭代進(jìn)行求解.其中,問(wèn)題(7)等價(jià)于:

    Bregmen迭代算法總結(jié)如算法2.

    算法2.Bregman迭代圖像邊緣檢測(cè)算法u1=u0,v1=?u0,w1=?u0 bk1=0,bk2=0初始化:設(shè)k=0,選擇 and ;for do 計(jì)算(8)的最優(yōu)解;bk+1 k=1,2,3,···,K 1 =bk1+vk+1??uk+1 更新:;2 =bk2+wk+1??uk+1 更新:;bk+1 end

    直接對(duì)(8)求解時(shí)十分困難,從而可把上式分離成多個(gè)單變量,即,把一個(gè)復(fù)雜問(wèn)題分解成幾個(gè)子問(wèn)題求解.Split Bregman iteration第一次被Goldstain和OSher[40]提出,它在很多寬松條件下等價(jià)于ADMM算法.該方法被廣泛地應(yīng)用到圖像噪聲去除,圖像去模糊,圖像修補(bǔ)等相關(guān)領(lǐng)域.

    3.1 u-子問(wèn)題

    對(duì)u-子問(wèn)題進(jìn)行求解時(shí)首先固定r=|f?c1|+|f?c2|.得到如下目標(biāo)泛函:

    此優(yōu)化問(wèn)題可以通過(guò)其優(yōu)化條件加以求解,在連續(xù)情況下,其解為四階線性PDE,離散情況,其優(yōu)化條件如下:

    通過(guò)FFT (快速傅里葉變換)可以得到(10)如下近似解.其中,F 代表傅里葉變換,F?1代表傅里葉的逆變換.

    3.2 v-子問(wèn)題

    v-子問(wèn)題等價(jià)于求解如下優(yōu)化問(wèn)題:

    上述優(yōu)化問(wèn)題的解可以使用收縮算子[42]表達(dá),其解析表達(dá)如下:

    其中,收縮算子shrink(x)[42]定義如下:

    對(duì)于任意點(diǎn)α ∈[0,1]2,為了方便起見(jiàn),令

    3.3 w-子問(wèn)題

    w-子問(wèn)題滿足如下優(yōu)化條件:

    類似于v-子問(wèn)題,上述優(yōu)化問(wèn)題解為具有收縮性質(zhì)的軟閾算子.

    為了說(shuō)明收縮算子的表達(dá)效果和設(shè)置依據(jù)我們考慮如下的一維最小值問(wèn)題加以說(shuō)明:

    其中,x,h都為標(biāo)量.令可得:

    當(dāng)h≥λα?xí)r ,Q(x)的 最小值為h?λα≥0;當(dāng)h≤?λα?xí)r,Q(x)的最小值為h+λα≤0 ;當(dāng)? λα

    收縮算子對(duì)稀疏解是有利的,最終大量的h值會(huì)趨近于區(qū)間[ ?λα,λα]的 某個(gè)值從而使得x變 得稀疏.當(dāng)λ 值越大,區(qū)間[ ?λα,λα]就 越大,從而使得x的解就越稀疏.

    對(duì)于本模型來(lái)說(shuō)若圖像的大小為N×N,以式(13)為例,其線性時(shí)間復(fù)雜度為O (N2),所以,v,w-子問(wèn)題可以被快速求解,綜上所述,對(duì)于v,w-子問(wèn)題我們?cè)O(shè)置了如式(14)的收縮算子是合理的.

    3.4 c 1,c2-子問(wèn)題

    現(xiàn)考慮如下最小值問(wèn)題,從而得到c1,c2的近似解.先固定u,得到如下最小值問(wèn)題:

    其中,當(dāng)i=1時(shí) ,h=u,當(dāng)i=2 時(shí) ,h=1?u.

    首先我們定義區(qū)域 [0,1]2上的特征函數(shù),χ=從而常數(shù)值ci可以用如下的表達(dá)式來(lái)代替χci,從而可以重寫最小值(20):

    其中,

    為使用變量分離法,故此引入如下輔助變量,

    式(23)可以通過(guò)ADMM或Split Bregman迭代進(jìn)行求解.

    從而,cki+1和eki+1的近似解解析表達(dá)式為:

    故此ADMM圖像邊緣檢測(cè)算法總結(jié)如算法3.

    算法3.ADMM圖像邊緣檢測(cè)算法λη1η2α β>0γ1,γ2>0,bk1=v0=0,bk2=w0=0,e0i=d0i=0 u0=1 ??[0,1]2初始化:設(shè)k=0,選擇,,,,,在某些區(qū)域 ,0其他;k=1,2,3,···,K for do hk=uk i=1hk=(1?uk),i=2 if ,,ck+1 i ck+1 i i=1,2 計(jì)算 :用式(27)計(jì)算 ,;ek+1 i ek+1 i i=1,2 計(jì)算 :用式(28)計(jì)算 ,;dk+1 i dk+1 i i=1,2f?ck+1 計(jì)算 :用 ;更新 :用式(26)計(jì)算 ,;rk+1 rk+1=f?ck+1 1 2 vk+1 vk+1 計(jì)算 :用式(13)計(jì)算 ;wk+1 wk+1 計(jì)算 :用式(16)計(jì)算 ;uk+1 uk+1 計(jì)算 :用式(11)計(jì)算 ;1 bk+11 =bk1+vk+1??uk+1 更新 :;bk+1 2 bk+12 =bk2+wk+1??uk+1 更新 :;bk+1 end

    4 算法部分收斂性分析

    本節(jié)中,受到文獻(xiàn)[38]啟發(fā),本文將給出本算法的部分收斂性分析,最小值問(wèn)題(6)可以通過(guò)引入變量分離辦法重寫成有約束條件:

    其對(duì)應(yīng)的拉格朗日泛函如下:

    其中,μ1,μ2,μ3為對(duì)偶變量.

    設(shè)X?=(c?1,c?2,u?,e?1,e?2,v?,w?,μ?1,μ?2,μ?3,μ?4)為 問(wèn) 題(29)的KKT點(diǎn),則點(diǎn)X?滿足如下的KKT條件:

    其中,x⊙y為coordi nate-wise乘法.

    定理.設(shè)Xk=(ck1,ck2,uk,ek1,ek2,vk,wk,d1k,d2k,bk1,bk2)為算法3中的迭代,并且,設(shè):

    假設(shè),

    證明:首先,對(duì)于u-子問(wèn)題等價(jià)于求解最小化問(wèn)題(9),從而得出如下優(yōu)化條件:

    其次,從優(yōu)化條件和ADMM算法對(duì)各變量的表達(dá)式可得:

    5 實(shí)驗(yàn)結(jié)果

    為了說(shuō)明本文提出算法優(yōu)越性,本文將給出該模型在灰度圖像,真實(shí)圖像,含噪聲圖像,CT圖像的分割實(shí)驗(yàn)結(jié)果,在討論結(jié)果之前,本文先作以下說(shuō)明,該算法的迭代終止條件應(yīng)滿足如下的相對(duì)容忍度,即滿足下列條件:

    此后,用該模型和Chan-Vese (CV)做比較.該算法與CV模型不同之處在于計(jì)算常數(shù)值ci,i=1,2,對(duì)于CV模型而言,ci,i=1,2在每次迭代過(guò)程中都有具體迭代表達(dá)式.該模型另外一個(gè)優(yōu)點(diǎn)是很少依賴于初始化位置和微調(diào)參數(shù),在本算法中,η1,η2,α ,β >0 ,γ1,γ2>0,所有參數(shù)都相對(duì)很固定.因?yàn)樵撃P褪且粋€(gè)局部極小值問(wèn)題,不同初始輪廓可能會(huì)收斂到不到能量局部極小值.我們之所以選擇和CV模型作比較,主要原因在于想體現(xiàn)高階正則模型與低階正則模型對(duì)圖像邊緣檢測(cè)不同效果影響.

    圖1為CV模型和新模型對(duì)飛機(jī)圖像的分割效果圖,其中,第一,二行分別代表CV模型和高階正則模型對(duì)飛機(jī)圖像的分割效果圖;第一列代表原始圖像,第二列分別代表CV模型和高階正則模型對(duì)飛機(jī)圖像的初始化輪廓;第三列分別代表CV模型對(duì)飛機(jī)圖像迭代800次時(shí)的邊緣檢測(cè)情況和高階正則模型對(duì)飛機(jī)圖像迭代80次時(shí)的邊緣檢測(cè)情況;第四列分別代表CV模型對(duì)飛機(jī)圖像迭代800次時(shí)的實(shí)際分割效果圖和高階正則模型對(duì)飛機(jī)圖像迭代80次時(shí)的實(shí)際分割效果圖.

    圖1 CV模型和高階正則模型對(duì)飛機(jī)圖像的檢測(cè)效果圖

    圖2表示CV模型對(duì)自然風(fēng)景圖像的分割效果圖,圖3表示高階正則模型對(duì)自然風(fēng)景圖像的分割效果圖.其中圖2的4幅圖分別為:第一幅代表原始圖像,第二幅代表CV模型對(duì)自然風(fēng)景圖像的初始化輪廓,第三幅代表CV模型對(duì)自然風(fēng)景這幅圖像迭代800次時(shí)的邊緣檢測(cè)圖,第四幅代表CV模型對(duì)自然風(fēng)景這幅圖迭代800次時(shí)的實(shí)際分割效果圖.圖3的4幅圖分別代表的意義為:第一幅為原始圖像,第二幅為高階正則模型對(duì)自然風(fēng)景圖像的初始化輪廓圖,第三幅圖像表示高階正則模型對(duì)自然風(fēng)景圖迭代200次時(shí)的邊緣檢測(cè)圖,第三幅圖像表示高階正則模型對(duì)自然風(fēng)景圖迭代200次時(shí)的實(shí)際分割效果圖.

    圖2 CV模型對(duì)自然風(fēng)景圖像的分割效果圖

    從圖4可以看出,CV模型對(duì)自然風(fēng)景這一幅圖像迭代30次后基本收斂,而高階正則模型需要迭代到40次左右才收斂,但是值得注意的是,從兩個(gè)模型對(duì)同一幅圖像Loss函數(shù)曲線圖可知,CV模型基本收斂后Loss函數(shù)曲線圖還有小幅度跳躍情況,然而,高階正則模型收斂后Loss函數(shù)曲線圖相對(duì)更平滑,說(shuō)明其分割效果也很好,更能精確地逼近物體邊界.

    下面我們給出一個(gè)例子說(shuō)明該高階正則模型在醫(yī)學(xué)圖像領(lǐng)域的應(yīng)用,圖像分割在醫(yī)學(xué)圖像領(lǐng)域的應(yīng)用就是通過(guò)發(fā)展大量的自動(dòng)或者半自動(dòng)圖像分割方法,精準(zhǔn)分割醫(yī)學(xué)圖像,從而替代醫(yī)生的大量手工標(biāo)注,從而使得醫(yī)生從大量繁重的體力勞動(dòng)中解放出來(lái).

    我們用本模型在病人肺部CT圖上做分割實(shí)驗(yàn),如圖5所示,從其分割效果圖得知,高模型不僅能對(duì)病人肺部做精準(zhǔn)分割,CV模型和高階正則模型Loss函數(shù)曲線收斂圖如圖6,可知,新模型Loss函數(shù)曲線圖相對(duì)CV模型而言相對(duì)較光滑,這更能說(shuō)明該模型比CV模型來(lái)更能精確地逼近物體邊界.

    需要說(shuō)明的是,在分割過(guò)程中,CV模型的模型參數(shù)設(shè)置均為原始論文的參數(shù),考慮到CV模型的數(shù)據(jù)項(xiàng)數(shù)擬合高斯噪聲,而本文所設(shè)計(jì)的模型數(shù)據(jù)保真項(xiàng)模擬的是外加噪聲,所以,在實(shí)驗(yàn)過(guò)程中,對(duì)于CV模型而言,并沒(méi)有人工地加入任何噪聲.但是,對(duì)于高階正則模型而言,由于本文假設(shè)其成像過(guò)程中可能含有外加噪聲(如:椒鹽噪聲),所以在實(shí)驗(yàn)過(guò)程中我們用Matlab中的加噪聲函數(shù)imnoise(Im0,'salt & pepper' 0.6);每幅原始圖像都人為的加了參數(shù)為0.6的椒鹽噪聲.新模型對(duì)所有圖像分割的迭代總次數(shù)均設(shè)置為200,CV迭代總次數(shù)均設(shè)為800,迭代停止條件為stoppling_tol=–10,λ =20,α =β=1,η1=η2=5 0,γ1=γ1=50.

    圖3 高階正則模型對(duì)自然風(fēng)景圖像的分割效果圖

    圖4 CV模型和高階正則模型對(duì)自然風(fēng)景圖像分割的Loss函數(shù)曲線收斂圖

    圖5 高階正則模型對(duì)病人肺部CT圖像的分割效果圖

    圖6 CV模型和高階正則模型對(duì)病人肺部CT圖像分割的Loss函數(shù)曲線收斂圖

    6 總結(jié)與展望

    本文提出一個(gè)修改的Chan-Vese模型,并引入高階正則函數(shù)對(duì)目標(biāo)函數(shù)進(jìn)行懲罰,然后對(duì)新模型用ADMM或者Split Bregmen迭代算法進(jìn)行數(shù)值求解.再者,本文還分析了該模型的數(shù)學(xué)性質(zhì),并給出該模型的部分收斂性分析,實(shí)驗(yàn)結(jié)果表明,通過(guò)引入高階正則函數(shù)后,該模型不僅能分割對(duì)比對(duì)低的物體,而且還可以探測(cè)帶噪聲的物體.大量實(shí)驗(yàn)表明,該算法在各領(lǐng)域具有廣泛應(yīng)用價(jià)值,在未來(lái)的工作中,我們會(huì)綜合一些深度學(xué)習(xí)的模塊,去探測(cè)物體的landmarks點(diǎn),從而使得算法更自動(dòng),更魯棒.

    猜你喜歡
    正則高階算子
    擬微分算子在Hp(ω)上的有界性
    有限圖上高階Yamabe型方程的非平凡解
    高階各向異性Cahn-Hilliard-Navier-Stokes系統(tǒng)的弱解
    滾動(dòng)軸承壽命高階計(jì)算與應(yīng)用
    哈爾濱軸承(2020年1期)2020-11-03 09:16:02
    各向異性次Laplace算子和擬p-次Laplace算子的Picone恒等式及其應(yīng)用
    剩余有限Minimax可解群的4階正則自同構(gòu)
    一類Markov模算子半群與相應(yīng)的算子值Dirichlet型刻畫
    類似于VNL環(huán)的環(huán)
    Roper-Suffridge延拓算子與Loewner鏈
    基于Bernstein多項(xiàng)式的配點(diǎn)法解高階常微分方程
    午夜免费成人在线视频| 亚洲国产精品合色在线| 两个人视频免费观看高清| 国产三级在线视频| 成年版毛片免费区| 免费观看精品视频网站| 亚洲人成伊人成综合网2020| 大型av网站在线播放| 国产成人精品久久二区二区91| 亚洲自拍偷在线| 久久婷婷成人综合色麻豆| 99精品在免费线老司机午夜| 多毛熟女@视频| 少妇裸体淫交视频免费看高清 | 一边摸一边做爽爽视频免费| 日韩一卡2卡3卡4卡2021年| 亚洲av熟女| 免费高清视频大片| 精品一区二区三区视频在线观看免费| 啦啦啦韩国在线观看视频| 国产成人欧美在线观看| 麻豆国产av国片精品| av有码第一页| 丰满的人妻完整版| 亚洲精品在线美女| 亚洲va日本ⅴa欧美va伊人久久| 人人澡人人妻人| 真人一进一出gif抽搐免费| 欧美老熟妇乱子伦牲交| 亚洲人成网站在线播放欧美日韩| 看黄色毛片网站| 首页视频小说图片口味搜索| 午夜福利视频1000在线观看 | 一本久久中文字幕| 久久午夜亚洲精品久久| 一个人免费在线观看的高清视频| 午夜影院日韩av| 9热在线视频观看99| 中文字幕av电影在线播放| 多毛熟女@视频| 国产精品久久视频播放| 亚洲第一av免费看| 免费无遮挡裸体视频| 在线国产一区二区在线| 最近最新中文字幕大全电影3 | 好男人电影高清在线观看| 曰老女人黄片| 一本综合久久免费| 高潮久久久久久久久久久不卡| 又大又爽又粗| 日本免费a在线| 一进一出抽搐gif免费好疼| 天天一区二区日本电影三级 | 国产精品久久久av美女十八| 国产一区二区三区视频了| 无限看片的www在线观看| 50天的宝宝边吃奶边哭怎么回事| 欧美成人免费av一区二区三区| 黄片大片在线免费观看| 中文字幕另类日韩欧美亚洲嫩草| 国产成年人精品一区二区| 纯流量卡能插随身wifi吗| 国产成人精品久久二区二区免费| 在线av久久热| 香蕉国产在线看| 久久人妻福利社区极品人妻图片| 日韩欧美一区二区三区在线观看| 国产高清视频在线播放一区| 天天一区二区日本电影三级 | 亚洲精品一区av在线观看| 久久香蕉精品热| 老司机深夜福利视频在线观看| 国产在线精品亚洲第一网站| 久久这里只有精品19| 99国产精品一区二区蜜桃av| 欧美日韩中文字幕国产精品一区二区三区 | 黄色丝袜av网址大全| 十分钟在线观看高清视频www| 在线观看免费日韩欧美大片| 欧美日韩乱码在线| 久久狼人影院| 亚洲精品国产区一区二| 不卡av一区二区三区| 99精品在免费线老司机午夜| 国产精品 欧美亚洲| 久9热在线精品视频| 国产精品久久视频播放| 精品国产乱码久久久久久男人| 精品国产美女av久久久久小说| 美女国产高潮福利片在线看| 久久精品亚洲熟妇少妇任你| 久久人人爽av亚洲精品天堂| 成人三级做爰电影| 人人妻人人澡人人看| 嫁个100分男人电影在线观看| а√天堂www在线а√下载| 国产精品一区二区三区四区久久 | 中文字幕色久视频| 久久久国产成人免费| 中亚洲国语对白在线视频| 叶爱在线成人免费视频播放| 欧美黑人欧美精品刺激| 久久精品国产亚洲av香蕉五月| 久久人人精品亚洲av| 黑人巨大精品欧美一区二区mp4| 欧美性长视频在线观看| 波多野结衣巨乳人妻| 欧美久久黑人一区二区| 国内久久婷婷六月综合欲色啪| 国内精品久久久久久久电影| 一级作爱视频免费观看| 欧美乱妇无乱码| 久久精品成人免费网站| 欧美成人性av电影在线观看| 露出奶头的视频| 女生性感内裤真人,穿戴方法视频| 丰满的人妻完整版| 亚洲自拍偷在线| 午夜福利高清视频| 亚洲三区欧美一区| 午夜激情av网站| 国产1区2区3区精品| 黄片播放在线免费| 国产aⅴ精品一区二区三区波| 欧美日韩亚洲国产一区二区在线观看| 热99re8久久精品国产| 欧美日韩福利视频一区二区| 午夜精品在线福利| 一本久久中文字幕| 男人的好看免费观看在线视频 | 级片在线观看| 色尼玛亚洲综合影院| 欧美一级毛片孕妇| 国产97色在线日韩免费| 国产精品美女特级片免费视频播放器 | 一区二区三区激情视频| 亚洲人成网站在线播放欧美日韩| 激情视频va一区二区三区| 成人国产综合亚洲| 又黄又爽又免费观看的视频| 中文字幕精品免费在线观看视频| 女人爽到高潮嗷嗷叫在线视频| 久久 成人 亚洲| 欧美日本视频| 黑人欧美特级aaaaaa片| 午夜福利成人在线免费观看| 人妻久久中文字幕网| 成人国语在线视频| 可以免费在线观看a视频的电影网站| 大码成人一级视频| 黄色丝袜av网址大全| 国产精品电影一区二区三区| 欧美日韩亚洲国产一区二区在线观看| 亚洲天堂国产精品一区在线| 50天的宝宝边吃奶边哭怎么回事| 精品国内亚洲2022精品成人| 免费在线观看日本一区| 日日爽夜夜爽网站| 国产亚洲精品综合一区在线观看 | 97碰自拍视频| 老熟妇仑乱视频hdxx| 在线观看免费视频网站a站| 国产免费av片在线观看野外av| 中文字幕人成人乱码亚洲影| 日本欧美视频一区| 国产免费av片在线观看野外av| 午夜福利,免费看| 日本精品一区二区三区蜜桃| 国产亚洲欧美在线一区二区| 日日干狠狠操夜夜爽| 久久亚洲真实| 中国美女看黄片| 99精品久久久久人妻精品| 午夜老司机福利片| bbb黄色大片| 熟女少妇亚洲综合色aaa.| 日韩大尺度精品在线看网址 | av中文乱码字幕在线| 巨乳人妻的诱惑在线观看| 男女下面进入的视频免费午夜 | 9191精品国产免费久久| 99国产精品99久久久久| 欧美激情久久久久久爽电影 | 19禁男女啪啪无遮挡网站| 成人亚洲精品一区在线观看| 99国产综合亚洲精品| 69av精品久久久久久| 一边摸一边抽搐一进一小说| bbb黄色大片| 亚洲情色 制服丝袜| 在线国产一区二区在线| 国产成+人综合+亚洲专区| 国语自产精品视频在线第100页| 亚洲中文字幕日韩| 每晚都被弄得嗷嗷叫到高潮| 国产亚洲精品第一综合不卡| 黄色女人牲交| 中文字幕人妻丝袜一区二区| 亚洲中文av在线| av在线天堂中文字幕| 久久精品91无色码中文字幕| 啪啪无遮挡十八禁网站| 桃红色精品国产亚洲av| 人人妻人人澡人人看| 极品人妻少妇av视频| 欧美日韩亚洲国产一区二区在线观看| 免费观看精品视频网站| 免费在线观看影片大全网站| 亚洲片人在线观看| 国产黄a三级三级三级人| 最近最新免费中文字幕在线| 如日韩欧美国产精品一区二区三区| 97超级碰碰碰精品色视频在线观看| 大码成人一级视频| 黄片播放在线免费| 国产精品亚洲美女久久久| 亚洲色图av天堂| 亚洲 国产 在线| 国产又爽黄色视频| 欧美黑人欧美精品刺激| 国产黄a三级三级三级人| 精品免费久久久久久久清纯| 叶爱在线成人免费视频播放| 一个人观看的视频www高清免费观看 | 91大片在线观看| 91精品国产国语对白视频| 久久精品亚洲精品国产色婷小说| 亚洲中文字幕日韩| 免费看a级黄色片| 国产免费av片在线观看野外av| 久久精品人人爽人人爽视色| 日本免费a在线| a在线观看视频网站| e午夜精品久久久久久久| 久久久久精品国产欧美久久久| 国内精品久久久久精免费| 色尼玛亚洲综合影院| 国产亚洲欧美精品永久| 国内精品久久久久精免费| 无限看片的www在线观看| 麻豆成人av在线观看| 18禁观看日本| 午夜影院日韩av| 女生性感内裤真人,穿戴方法视频| 日本免费a在线| 一级作爱视频免费观看| 多毛熟女@视频| 国产精品乱码一区二三区的特点 | 亚洲性夜色夜夜综合| 两个人视频免费观看高清| 多毛熟女@视频| 日本三级黄在线观看| 脱女人内裤的视频| 国语自产精品视频在线第100页| 伊人久久大香线蕉亚洲五| 在线永久观看黄色视频| 香蕉丝袜av| 久久香蕉激情| 中文字幕久久专区| 亚洲中文字幕日韩| 亚洲,欧美精品.| 无人区码免费观看不卡| 亚洲国产日韩欧美精品在线观看 | 脱女人内裤的视频| 成人精品一区二区免费| 婷婷精品国产亚洲av在线| 91av网站免费观看| 久久久国产成人精品二区| 搡老熟女国产l中国老女人| 久久人人97超碰香蕉20202| 久久精品国产清高在天天线| 国产成人精品久久二区二区91| 麻豆国产av国片精品| √禁漫天堂资源中文www| 日韩欧美国产在线观看| 国产亚洲精品一区二区www| 亚洲精品中文字幕在线视频| 午夜福利欧美成人| 久久久精品欧美日韩精品| 欧美乱妇无乱码| 中文字幕av电影在线播放| 精品不卡国产一区二区三区| 天天躁夜夜躁狠狠躁躁| 宅男免费午夜| 叶爱在线成人免费视频播放| 三级毛片av免费| 啦啦啦观看免费观看视频高清 | 久热这里只有精品99| 嫩草影院精品99| 久久久国产成人免费| 成人18禁高潮啪啪吃奶动态图| 女性被躁到高潮视频| 巨乳人妻的诱惑在线观看| 久久天躁狠狠躁夜夜2o2o| xxx96com| 神马国产精品三级电影在线观看 | 真人做人爱边吃奶动态| 亚洲电影在线观看av| 搡老岳熟女国产| 日韩三级视频一区二区三区| 久久精品国产亚洲av高清一级| 亚洲人成伊人成综合网2020| 最近最新中文字幕大全电影3 | 欧美成狂野欧美在线观看| 精品久久久久久成人av| 99精品在免费线老司机午夜| 亚洲国产高清在线一区二区三 | 女人爽到高潮嗷嗷叫在线视频| 精品国产超薄肉色丝袜足j| 国产av一区二区精品久久| av在线天堂中文字幕| 国产av又大| 久久久久久亚洲精品国产蜜桃av| 黄色丝袜av网址大全| 久久精品亚洲精品国产色婷小说| 国产精品综合久久久久久久免费 | 亚洲一码二码三码区别大吗| 午夜老司机福利片| 欧美最黄视频在线播放免费| 一级黄色大片毛片| 国产在线观看jvid| 亚洲精品美女久久久久99蜜臀| 夜夜爽天天搞| 亚洲成国产人片在线观看| √禁漫天堂资源中文www| 一本大道久久a久久精品| 久久久久久人人人人人| 黄片小视频在线播放| 啪啪无遮挡十八禁网站| 十分钟在线观看高清视频www| 99精品在免费线老司机午夜| 亚洲五月色婷婷综合| 中文字幕人成人乱码亚洲影| 一个人免费在线观看的高清视频| a在线观看视频网站| 老司机福利观看| 欧美乱色亚洲激情| 97碰自拍视频| 99在线人妻在线中文字幕| 两个人视频免费观看高清| 欧美色欧美亚洲另类二区 | 久久狼人影院| 香蕉丝袜av| e午夜精品久久久久久久| 久久国产亚洲av麻豆专区| 亚洲中文字幕日韩| 久久精品aⅴ一区二区三区四区| 国产亚洲精品久久久久5区| 精品久久蜜臀av无| 国语自产精品视频在线第100页| 变态另类丝袜制服| 亚洲人成伊人成综合网2020| 国内精品久久久久精免费| 亚洲av五月六月丁香网| 99国产极品粉嫩在线观看| 日韩三级视频一区二区三区| 麻豆国产av国片精品| 一级毛片女人18水好多| 国产午夜福利久久久久久| 亚洲成人国产一区在线观看| 在线观看午夜福利视频| 亚洲精品国产色婷婷电影| 97人妻天天添夜夜摸| 国产又爽黄色视频| 变态另类丝袜制服| 亚洲成av人片免费观看| 黑人欧美特级aaaaaa片| 中文亚洲av片在线观看爽| 国产91精品成人一区二区三区| 日本免费a在线| 精品熟女少妇八av免费久了| 国产99久久九九免费精品| 精品国产超薄肉色丝袜足j| 人人妻人人爽人人添夜夜欢视频| 国产精品久久视频播放| 多毛熟女@视频| 成年版毛片免费区| 99re在线观看精品视频| 狠狠狠狠99中文字幕| 亚洲第一欧美日韩一区二区三区| 一本大道久久a久久精品| av欧美777| 99香蕉大伊视频| 国产精品亚洲美女久久久| 男女床上黄色一级片免费看| 欧洲精品卡2卡3卡4卡5卡区| 国产av在哪里看| 一区二区三区精品91| 日韩国内少妇激情av| 亚洲无线在线观看| 正在播放国产对白刺激| 变态另类成人亚洲欧美熟女 | 黄色片一级片一级黄色片| 色播在线永久视频| 国产三级在线视频| 黄频高清免费视频| www国产在线视频色| 欧美日韩瑟瑟在线播放| 亚洲av美国av| 精品欧美一区二区三区在线| 日本免费一区二区三区高清不卡 | 国产精品精品国产色婷婷| 欧美另类亚洲清纯唯美| 日韩精品中文字幕看吧| 亚洲熟妇中文字幕五十中出| 欧美一级毛片孕妇| 亚洲精品在线美女| 女性生殖器流出的白浆| 色综合欧美亚洲国产小说| 国产亚洲精品久久久久久毛片| 国产黄a三级三级三级人| 国产精品 欧美亚洲| 不卡一级毛片| 亚洲人成77777在线视频| 亚洲第一欧美日韩一区二区三区| 老司机在亚洲福利影院| 亚洲,欧美精品.| 欧美日韩一级在线毛片| 热re99久久国产66热| 亚洲国产精品999在线| 夜夜爽天天搞| 亚洲 欧美一区二区三区| 亚洲全国av大片| 久久中文字幕人妻熟女| 亚洲第一电影网av| 男人舔女人的私密视频| av超薄肉色丝袜交足视频| 黑人欧美特级aaaaaa片| 午夜a级毛片| 又大又爽又粗| 国产精品野战在线观看| 久久人人精品亚洲av| 一个人观看的视频www高清免费观看 | 色哟哟哟哟哟哟| 变态另类丝袜制服| 日本vs欧美在线观看视频| 日韩国内少妇激情av| 精品国产美女av久久久久小说| 亚洲色图 男人天堂 中文字幕| 可以在线观看毛片的网站| 久久国产乱子伦精品免费另类| 69av精品久久久久久| av天堂久久9| 久久久久久人人人人人| 欧美激情高清一区二区三区| 一区二区三区激情视频| 国产精品香港三级国产av潘金莲| 一本久久中文字幕| 亚洲午夜理论影院| 黄色女人牲交| 亚洲一卡2卡3卡4卡5卡精品中文| 黑人操中国人逼视频| av在线播放免费不卡| 久久精品成人免费网站| 亚洲精品av麻豆狂野| 日韩 欧美 亚洲 中文字幕| www.www免费av| 欧美乱色亚洲激情| 国产精品,欧美在线| 午夜老司机福利片| 亚洲最大成人中文| 亚洲久久久国产精品| 国产精品亚洲av一区麻豆| 国产视频一区二区在线看| 黄色视频,在线免费观看| 日本在线视频免费播放| 久久国产乱子伦精品免费另类| 男人的好看免费观看在线视频 | 亚洲人成电影观看| 午夜日韩欧美国产| 国产一卡二卡三卡精品| 99久久精品国产亚洲精品| 国产乱人伦免费视频| 国产午夜精品久久久久久| 色av中文字幕| 悠悠久久av| 青草久久国产| 一区二区三区激情视频| 国产亚洲精品久久久久5区| 91老司机精品| 色在线成人网| 好男人在线观看高清免费视频 | 欧美 亚洲 国产 日韩一| 亚洲熟妇中文字幕五十中出| 久久青草综合色| 精品久久久久久,| 精品欧美国产一区二区三| 亚洲自偷自拍图片 自拍| 男女床上黄色一级片免费看| 精品久久蜜臀av无| 午夜福利高清视频| 夜夜爽天天搞| 国产亚洲精品综合一区在线观看 | 国产人伦9x9x在线观看| 久久亚洲精品不卡| 色av中文字幕| 国产真人三级小视频在线观看| 三级毛片av免费| 动漫黄色视频在线观看| 香蕉丝袜av| 黄色毛片三级朝国网站| 久久精品国产清高在天天线| 女人精品久久久久毛片| 老司机靠b影院| 中国美女看黄片| 亚洲一区中文字幕在线| 黄色视频,在线免费观看| 午夜免费鲁丝| 丰满的人妻完整版| 成年女人毛片免费观看观看9| 国产高清视频在线播放一区| 日韩大码丰满熟妇| 日韩精品中文字幕看吧| 国产精品自产拍在线观看55亚洲| 精品第一国产精品| 中文字幕人成人乱码亚洲影| 性少妇av在线| 999久久久国产精品视频| 黄片大片在线免费观看| 日韩精品中文字幕看吧| 老鸭窝网址在线观看| 黄色成人免费大全| 麻豆av在线久日| 亚洲三区欧美一区| 午夜福利视频1000在线观看 | 两个人视频免费观看高清| 欧美乱色亚洲激情| 国产精品自产拍在线观看55亚洲| 精品国产国语对白av| 久久久国产精品麻豆| 精品国内亚洲2022精品成人| 天堂动漫精品| 窝窝影院91人妻| 女性生殖器流出的白浆| 可以在线观看毛片的网站| 在线天堂中文资源库| 亚洲色图综合在线观看| 一二三四社区在线视频社区8| 久久热在线av| 精品乱码久久久久久99久播| 午夜久久久久精精品| 国产精品免费一区二区三区在线| 午夜久久久久精精品| 国产91精品成人一区二区三区| 国产精品一区二区在线不卡| 亚洲第一青青草原| 国产精品野战在线观看| 欧美日韩亚洲综合一区二区三区_| 香蕉久久夜色| 丝袜美腿诱惑在线| 搞女人的毛片| 男女下面插进去视频免费观看| 久久精品亚洲熟妇少妇任你| 精品免费久久久久久久清纯| 岛国在线观看网站| 久久久久久久久久久久大奶| 国产精品免费视频内射| 1024视频免费在线观看| 97人妻天天添夜夜摸| 51午夜福利影视在线观看| 亚洲一区中文字幕在线| 亚洲精华国产精华精| 欧美成人午夜精品| 大香蕉久久成人网| 亚洲精品在线美女| 一区二区三区精品91| 亚洲国产高清在线一区二区三 | 精品久久久久久,| 首页视频小说图片口味搜索| 国产成人系列免费观看| 国产精品亚洲av一区麻豆| av在线播放免费不卡| 国产av一区在线观看免费| 色综合亚洲欧美另类图片| 精品第一国产精品| 国产一区二区三区综合在线观看| 久久 成人 亚洲| 90打野战视频偷拍视频| 久久精品国产综合久久久| 极品教师在线免费播放| 超碰成人久久| 电影成人av| 黑人巨大精品欧美一区二区mp4| 十八禁人妻一区二区| 亚洲熟女毛片儿| 高清毛片免费观看视频网站| 精品少妇一区二区三区视频日本电影| 久久久久九九精品影院| 久久久久久免费高清国产稀缺| 香蕉国产在线看| 国产欧美日韩一区二区三区在线| 99久久国产精品久久久| 午夜免费成人在线视频| 欧美黄色淫秽网站| 国产成人系列免费观看| 精品久久久久久成人av| 久久精品国产清高在天天线| 制服人妻中文乱码| 在线免费观看的www视频| 人妻丰满熟妇av一区二区三区| 变态另类丝袜制服| 99国产精品99久久久久| 91成年电影在线观看| 亚洲专区字幕在线| 日韩av在线大香蕉| 侵犯人妻中文字幕一二三四区| 亚洲五月婷婷丁香| 亚洲中文字幕日韩| 男人操女人黄网站| 成人国产一区最新在线观看| 亚洲欧美激情综合另类| 男女午夜视频在线观看| 99久久综合精品五月天人人| 久久精品aⅴ一区二区三区四区| 国产片内射在线|