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

    基于多尺度邊緣保持分解與PCNN 的醫(yī)學(xué)圖像融合

    2021-03-18 08:04:10郭淑娟秦品樂王麗芳
    計算機(jī)工程 2021年3期
    關(guān)鍵詞:濾波器尺度邊緣

    郭淑娟,高 媛,秦品樂,王麗芳

    (中北大學(xué)大數(shù)據(jù)學(xué)院山西省生物醫(yī)學(xué)成像與影像大數(shù)據(jù)重點實驗室,太原 030051)

    0 概述

    圖像融合是指將不同模態(tài)的相關(guān)圖像合并為一幅新圖像。在此過程中,要求生成的融合圖像具有良好的清晰度,并使不同模態(tài)的圖像特征盡可能融合到一幅圖像中。圖像融合在多模態(tài)醫(yī)學(xué)圖像中的應(yīng)用是一個熱門的研究方向。隨著現(xiàn)代醫(yī)學(xué)成像技術(shù)的發(fā)展,結(jié)合不同的醫(yī)學(xué)成像方式能夠獲得器官和組織的全面信息,例如計算機(jī)斷層掃描(Computed Tomography,CT)圖像僅顯示高密度組織,磁共振(Magnetic Resonance,MR)圖像僅顯示高水分組織,正電子發(fā)射斷層掃描(Positron Emission Tomography,PET)圖像僅反映不同人體之間組織強(qiáng)度的差異[1],而將這些醫(yī)學(xué)圖像的互補(bǔ)特征融合到一幅圖像中,能夠提高醫(yī)生的診斷效率和準(zhǔn)確性。

    圖像融合技術(shù)可分為空間域方法和變換域方法兩類,其中,變換域方法可以更好地表征圖像特征,彌補(bǔ)空間域方法在細(xì)節(jié)提取方面的不足。在變換域方法中,多尺度分析方法通過將源圖像映射至不同的尺度空間提取潛在的重要信息,符合人眼視覺系統(tǒng)的生理機(jī)制,被認(rèn)為是一種主流的圖像融合方法[2]?,F(xiàn)有的多尺度分析方法主要有基于離散小波變換[3]、基于剪切波變換[4]、基于輪廓波變換[5]、基于非下采樣輪廓波變換(Non-Subsampled Contourlet Transform,NSCT)[6]和基于非下采樣剪切波變換(Non-Subsampled Shearlet Transform,NSST)[7]的算法,其中,基于NSCT 和NSST 的算法具有平移不變性,能夠很好地抑制偽吉布斯現(xiàn)象。

    脈沖耦合神經(jīng)網(wǎng)絡(luò)(Pulse Coupled Neural Network,PCNN)是一種通過研究哺乳動物視覺皮層而獲得的具有仿生機(jī)制的單層神經(jīng)網(wǎng)絡(luò)模型[8],其不需要學(xué)習(xí)或者訓(xùn)練,且具有全局耦合性和脈沖發(fā)射特性,能夠提取圖像的全局特征并且增強(qiáng)圖像細(xì)節(jié),有利于圖像的實時處理。為提高融合圖像的視覺感知效果并充分利用圖像的空間信息,一些多尺度分析與PCNN 相結(jié)合的圖像融合方法相繼被提出[9],如NSCT-SF-PCNN[10]和NSST-PAPCNN[11]。然而,傳統(tǒng)的多尺度分析方法大多采用線性濾波器,由于無法保留邊緣從而導(dǎo)致分解階段的強(qiáng)邊緣處出現(xiàn)模糊,使得邊緣處產(chǎn)生光暈[12]。

    為解決該問題,研究人員將保邊濾波器引入圖像融合領(lǐng)域。保邊濾波器具有平移不變性、空間一致性和邊緣保護(hù)性能,而且計算效率較高。在此方面:文獻(xiàn)[13]提出一種多尺度雙邊濾波器對源圖像進(jìn)行分解;文獻(xiàn)[14]提出一種基于引導(dǎo)濾波的圖像融合方案,其對源圖像進(jìn)行兩尺度分解,通過引導(dǎo)濾波優(yōu)化權(quán)重實現(xiàn)了空間一致性;文獻(xiàn)[15]提出以多尺度極值濾波對源圖像進(jìn)行多尺度分解,采用對比度融合規(guī)則使融合圖像更符合人眼視覺感知特性;文獻(xiàn)[16]提出一種加權(quán)最小二乘濾波與引導(dǎo)濾波相結(jié)合的新方案,與只基于引導(dǎo)濾波的方法相比,該方案可以獲得更好的效果。盡管上述方法在一定程度上保護(hù)了邊緣信息,但是基礎(chǔ)層主要包含粗尺度結(jié)構(gòu)信息,直接對其融合可能在之后的融合過程中丟失或模糊保邊濾波保留的顯著的邊緣特征[17]。

    為保留圖像完整清晰的邊緣信息,提高融合圖像的視覺感知效果,本文提出一種新的醫(yī)學(xué)圖像融合方法。將多尺度邊緣保持分解方法引入圖像分解過程,對已配準(zhǔn)的源圖像進(jìn)行分解得到低頻層和高頻層。在此基礎(chǔ)上,利用改進(jìn)的空間頻率以及區(qū)域能量激勵PCNN 對兩者進(jìn)行融合,并通過逆變換獲得最終的融合圖像。

    1 相關(guān)理論

    1.1 加權(quán)最小二乘濾波與高斯濾波

    文獻(xiàn)[18]提出一種基于加權(quán)最小二乘(Weighted Least Square,WLS)優(yōu)化框架的非線性濾波器。給定一個輸入圖像g,WLS 濾波器的目的是使濾波后的平滑圖像u盡可能與輸入圖像g近似,則濾波后的圖像u可以表示為以下目標(biāo)函數(shù):

    其中,下標(biāo)p表示像素的空間位置,(up-gp)2用以保證濾波后的平滑圖像u更接近輸入圖像g,是一個正則化項,通過求u的偏導(dǎo)數(shù)來實現(xiàn)平滑,λ是正則化參數(shù),用以調(diào)整兩項之間的平衡,ax,p和ay,p代表平滑因子。ax,p和ay,p的計算公式如下:

    其中,l是輸入圖像g的對數(shù)亮度通道,指數(shù)α確定對g的梯度的敏感度,而ε是一個非常小的常數(shù),通常為0.000 1,可防止在g恒定的區(qū)域被零除。

    高斯濾波是一種線性濾波,由于沒有考慮相鄰像素的影響,其對圖像濾波時去除高頻信息的同時也模糊了邊緣,因此在圖像的平滑和去噪中被廣泛應(yīng)用。高斯濾波公式可表示為:

    1.2 非下采樣方向濾波器組

    文獻(xiàn)[19]提出一種方向濾波器組(Directional Filter Bank,DFB)。DFB 將圖像劃分為具有方向性的子帶,捕捉圖像幾何特征,但是其中的下采樣操作使其缺乏平移不變性。非下采樣方向濾波器組(Non-subsampled Directional Filter Bank,NDFB)[20]與DBF 類似,可以對圖像進(jìn)行多方向分解,提取圖像方向細(xì)節(jié)特征,但是其進(jìn)行上采樣操作,因而能夠獲得位移不變的方向擴(kuò)展。NDFB 是由扇形濾波器和棋盤濾波器構(gòu)成的四方向濾波器組,基本結(jié)構(gòu)如圖1 所示。

    圖1 NDFB 基本結(jié)構(gòu)Fig.1 Basic structure of NDFB

    1.3 脈沖耦合神經(jīng)網(wǎng)絡(luò)

    圖2 PCNN 的簡化模型Fig.2 Simplified model of PCNN

    脈沖耦合神經(jīng)網(wǎng)絡(luò)(PCNN)是一個M×N的二維網(wǎng)絡(luò),其中每一個神經(jīng)元對應(yīng)于圖像的一個特定的像素。每個神經(jīng)元由三部分組成,分別是接收域、調(diào)制域和脈沖產(chǎn)生器。一種PCNN 的簡化模型和數(shù)學(xué)表達(dá)式如圖2 和式(5)所示。其中:下標(biāo)L表示圖像解等級;Fij(n)和Lij(n)分別是PCNN 的饋送輸入和鏈接輸入,處于PCNN 模型的接收域,F(xiàn)ij(n)通常接收的是圖像像素的歸一化灰度值,Lij(n)則通過突觸權(quán)重與8 個相鄰神經(jīng)元的先前放電狀態(tài)相關(guān);aL和VL是Lij(n)的時間衰減常數(shù)和振幅增益;Uij(n)代表內(nèi)部活動,處于PCNN 模型的調(diào)制域,通過Fij(n)和Lij(n)進(jìn)行非線性調(diào)制得到,其中參數(shù)β是連接強(qiáng)度。

    脈沖發(fā)生域控制脈沖的發(fā)生:將Uij與動態(tài)閾值θij進(jìn)行比較,若Uij>θij,Yij=1 即神經(jīng)元點火,產(chǎn)生脈沖輸出;反之Yij=0,aθ和Vθ是θij的時間衰減常數(shù)和振幅增益。將上述過程多次迭代直到滿足設(shè)定的迭代條件,把此時由神經(jīng)元點火形成的點火映射圖作為輸出結(jié)果。

    2 本文醫(yī)學(xué)圖像融合方法

    對已配準(zhǔn)的源圖像進(jìn)行圖像融合的主要步驟包括圖像分解、圖像融合和圖像重構(gòu)。本文在圖像分解過程中引入多尺度邊緣保持分解方法,對已配準(zhǔn)的源圖像進(jìn)行分解得到低頻層和高頻層,而在圖像融合過程中針對得到的低頻層和高頻層不同的特點,分別采用不同的融合規(guī)則進(jìn)行融合,在圖像重構(gòu)中獲得最終的融合的圖像。融合過程如圖3所示。

    圖3 醫(yī)學(xué)圖像融合過程Fig.3 Process of medical image fusion

    2.1 多尺度邊緣保持分解

    WLS 濾波器可以有效地在模糊和銳化之間進(jìn)行最佳折衷,更好地保留邊緣,解決雙邊濾波在邊緣處容易產(chǎn)生梯度反轉(zhuǎn)的問題,彌補(bǔ)雙邊濾波無法很好地在多尺度上提取到細(xì)節(jié)信息的不足,從而更準(zhǔn)確地提取多尺度信息。由于高斯濾波器只考慮像素的空間分布而沒有考慮鄰近像素的影響,因此具有對圖像進(jìn)行濾波時能夠剔除細(xì)節(jié)信息、模糊邊緣的特點。結(jié)合兩種濾波的特性,本文對源圖像進(jìn)行多尺度分解。分解過程由以下步驟組成:

    1)對輸入圖像進(jìn)行WLS 濾波分解得到細(xì)節(jié)層和基礎(chǔ)層。使用WLS 濾波對輸入圖像進(jìn)行保邊平滑處理得到基礎(chǔ)層,細(xì)節(jié)層則由輸入圖像與基礎(chǔ)圖像做差得到?;A(chǔ)層按照式(1)計算。設(shè)Cinput是輸入圖像,則細(xì)節(jié)層表示為:

    其中,CD_layer和CB_layer是經(jīng)過WLS 濾波分解得到的細(xì)節(jié)層和基礎(chǔ)層。由于WLS 濾波的保邊特性,因此基礎(chǔ)層含有強(qiáng)邊緣信息。

    2)為避免光譜損失,使用高斯濾波器再次分解基礎(chǔ)層提取邊緣信息,得到低頻層和邊緣層。低頻層圖像依據(jù)式(4)計算,則邊緣層計算公式為:

    其中,CL_layer和CE_layer是基礎(chǔ)層經(jīng)過高斯濾波分解得到的低頻層和邊緣層?;赪LS 濾波和高斯濾波的分解單元結(jié)構(gòu)如圖4 所示。

    圖4 基于WLS 濾波和高斯濾波的分解單元結(jié)構(gòu)Fig.4 Decomposition unit structure based on WLS filtering and Gaussian filtering

    3)將邊緣層CE_layer和細(xì)節(jié)層CD_layer進(jìn)行疊加,構(gòu)建圖像的高頻層CH_layer。同時為提取高頻層系數(shù)中存在的大量幾何特征,對高頻層使用NSDFB 進(jìn)行方向分析,捕獲高頻層系數(shù)不同方向的特征。邊緣保持分解框架如圖5 所示。

    圖5 邊緣保持分解框架Fig.5 Frame of edge-preserving decomposition

    將得到的低頻層CL_layer作為下一級分解的輸入圖像,經(jīng)過多級分解實現(xiàn)多尺度分解。通過上述分解模式最終將源圖像分解為一個低頻層和一系列不同方向的高頻層方向子帶。將源圖像A 和源圖像B通過多尺度分解得到的低頻層表示為和,不同方向的高頻層方向子帶表示為和,分別表示源圖像A 和源圖像B 在第i級別第d方向上的高頻層。

    2.2 圖像融合規(guī)則

    2.2.1 低頻層融合規(guī)則

    低頻層圖像是源圖像去除高頻信息之后的近似圖像,包含圖像大量的能量信息,其中子帶相鄰系數(shù)之間存在區(qū)域相關(guān)性[21]。區(qū)域能量將鄰域內(nèi)每個中心像素的能量表征為其自身和鄰域內(nèi)像素值的平方和,考慮了圖像局部特征,減少了融合圖像中的不連續(xù)性,邊緣信息在該融合模式下能夠很好地得以保留。因此,本文將區(qū)域能量作為PCNN的外部激勵融合低頻子帶。首先計算低頻層的區(qū)域能量作為PCNN 的外部激勵,然后使用表征圖像清晰度的平均梯度作為PCNN 的連接強(qiáng)度,最后通過PCNN 輸出的點火圖進(jìn)行融合。低頻層融合過程如下:

    其中,W是一個3×3 的滑動窗口。為體現(xiàn)低頻層子帶的平滑特性,選取W為:

    其中,η是一個大于0 的常數(shù),用于調(diào)節(jié)鏈接強(qiáng)度的值,本文設(shè)置為0.2(i,j)是計算所得的平均梯度。的計算公式如下:

    3)將上述計算的區(qū)域能量作為PCNN 模型的外部激勵,平均梯度構(gòu)建的鏈接系數(shù)作為PCNN 模型的鏈接強(qiáng)度,通過PCNN 模型得到低頻子帶相應(yīng)的點火圖

    2.2.2 高頻層融合規(guī)則

    高頻層反映了圖像的紋理細(xì)節(jié)和邊緣信息,對最終的融合結(jié)果有著至關(guān)重要的影響。改進(jìn)的空間頻率能夠反映圖像的邊緣和細(xì)節(jié)信息,同時表示圖像灰度值的變化,因此,本文采用改進(jìn)的空間頻率激PCNN 對高頻層進(jìn)行融合。首先計算每一級不同方向的高頻層方向子帶改進(jìn)的空間頻率作為PCNN 的外部刺激,然后使用表征圖像清晰度的平均梯度作為PCNN 的連接強(qiáng)度,最后通過PCNN 輸出的點火圖進(jìn)行融合。高頻層融合過程如下:

    3)將上述計算的改進(jìn)的空間頻率作為PCNN 模型的外部激勵,將平均梯度構(gòu)建的鏈接系數(shù)作為PCNN 模型的鏈接強(qiáng)度,通過PCNN 模型得到高頻方向子帶相應(yīng)的點火圖

    通過NDFB 逆變換,可以得到第i層的高頻層融合結(jié)果

    2.3 圖像重構(gòu)

    將得到的高頻層融合結(jié)果和低頻層融合結(jié)果重構(gòu)為最終的融合結(jié)果:

    3 實驗結(jié)果與分析

    3.1 實驗參數(shù)設(shè)置

    本文實驗數(shù)據(jù)選取自哈佛大學(xué)醫(yī)學(xué)院的公開數(shù)據(jù),配準(zhǔn)的醫(yī)學(xué)圖像是大小為256 像素×256 像素、顏色深度為8 bit 的灰度圖。

    在圖像分解階段,基于WLS 分解的參數(shù)主要包括平滑參數(shù)λ和邊緣保留參數(shù)α。λ是一個正則化參數(shù),用于保持?jǐn)?shù)據(jù)項和平滑度項之間的平衡,增加λ可產(chǎn)生更平滑的圖像。α通過非線性縮放梯度來確定對梯度的敏感程度,增加α可保留更清晰的邊緣。本文參考文獻(xiàn)[16],設(shè)置λi=[0.075,0.6,4.8],α=0.6,其中,i=1,2,…,L?;诟咚篂V波分解的主要參數(shù)是用于控制圖像平滑效果的高斯分布的標(biāo)準(zhǔn)差σ,第1 級設(shè)置為σ1=1.1,其他級設(shè)置為σi+1=kσi,其中,k是相鄰級之間的分解比例因子,設(shè)置為k=2。圖像分解等級設(shè)置為L=1,2,3,方向分析的分解水平分別設(shè)置為第1 層方向數(shù)為16 方向,第2 層和第3 層分解分別為8 方向和4 方向,采用的方向濾波器為“vk”,將第L級的得到的低頻圖像作為下一級分解的輸入圖像以實現(xiàn)多尺度分解。

    在圖像融合階段,本文參考文獻(xiàn)[9],將用于融合圖像的PCNN 的參數(shù)設(shè)置為:aL=10,aθ=0.815,VL=1.0,Vθ=10,n=200

    3.2 實驗結(jié)果分析

    3.2.1 不同融合方法的比較

    本文采用多發(fā)性腦梗塞MR-T1 和MR-T2 圖像、腦弓形蟲CT 和MR 圖像、正常腦部腦CT 和MR 圖像進(jìn)行3 組仿真實驗,在64 位Windows10 操作系統(tǒng)、Matlab2014b 平臺上進(jìn)行仿真。選取以下5 種融合方法作為對比方法:1)基于NSCT 的圖像融合方法(NSCT);2)基于多尺度變換和稀疏表示的圖像融合方法(MST-SR)[22];3)基于NSCT域下空間頻率激勵PCNN的圖像融合方法(NSCT-SF-PCNN)[10];4)基于NSST 與自適應(yīng)PCNN 的圖像融合方法(NSST-PAPCNN)[11];5)基于引導(dǎo)濾波的圖像融合方法(GFF)[14]。NSCT采用高頻子帶區(qū)域能量取大、低頻子帶加權(quán)平均的融合規(guī)則,MST-SR 的多尺度分解方法選取為LP,分解等級為4。6 種方法的實驗結(jié)果如圖6~圖8 所示。

    圖6 不同方法的多發(fā)性腦梗塞MR-T1/MR-T2圖像融合結(jié)果Fig.6 Fusion results for MR-T1/MR-T2 image of multiple cerebral infarction by different methods

    圖7 不同方法的腦弓形蟲CT/MR 圖像融合結(jié)果Fig.7 Fusion results for CT/MR image of toxoplasma gondii by different methods

    圖8 不同方法的正常腦部CT/MR 圖像融合結(jié)果Fig.8 Fusion results for CT/MR image of normal brain by different methods

    由實驗結(jié)果可以看出:NSCT 和NSCT-SF-PCNN的融合結(jié)果丟失了骨質(zhì)信息,同時NSCT 的融合圖像整體亮度偏低;NSCT-SF-PCNN 在第3 組實驗中的融合結(jié)果MR 信息嚴(yán)重缺失并且出現(xiàn)邊緣模糊的現(xiàn)象;NSST-PAPCNN 的融合結(jié)果雖然較好地保存了骨質(zhì)信息,但是在邊緣處產(chǎn)生了邊緣模糊和光暈現(xiàn)象;MST-SR 能夠較好地保存骨質(zhì)信息和邊緣信息,但是在第2 組實驗中MST-SR 融合圖像MR信息有所缺失,存在對比度下降和細(xì)節(jié)丟失的問題;GFF 的方法能夠有效保護(hù)邊緣信息,但是第1 組和第2 組實驗結(jié)果存在明顯的細(xì)節(jié)丟失現(xiàn)象。

    為更客觀地評價圖像的融合結(jié)果,選取邊緣評價因子QAB/F、平均梯度、互信息、標(biāo)準(zhǔn)差和空間頻率5 種指標(biāo)進(jìn)行客觀評價。3 組實驗的客觀評價指標(biāo)如表1~表3 所示,其中加粗?jǐn)?shù)據(jù)表示最優(yōu)值??梢钥闯觯涸? 種評價指標(biāo)中,各組實驗中本文方法的QAB/F均達(dá)到最大值,驗證了其在邊緣保持方面的可行性和有效性;此外,除第1 組實驗的空間頻率和第2 組、第3 組實驗的標(biāo)準(zhǔn)差之外,本文方法其余各實驗指標(biāo)均達(dá)到最大值,表明其能較好地保留原圖像信息,體現(xiàn)圖像的紋理細(xì)節(jié)。

    表1 不同方法的多發(fā)性腦梗塞MR-T1/MR-T2 圖像融合評價指標(biāo)Table 1 Evaluation indexes of MR-T1/MR-T2 image fusion of multiple cerebral infarction by different methods

    表2 不同方法的腦弓形蟲CT/MRI 圖像融合評價指標(biāo)Table 2 Evaluation indexes of CT/MRI image fusion of toxoplasma gondii by different methods

    表3 不同方法的正常腦部CT/MRI 圖像融合評價指標(biāo)Table 3 Evaluation indexes of CT/MRI image fusion of normal brain by different methods

    通過對以上主觀視覺和客觀指標(biāo)的分析可知,本文方法能夠很好地保留邊緣信息,所得圖像具有較高的對比度和豐富的細(xì)節(jié)信息,較對比方法能夠得到更好的視覺效果。

    3.2.2 不同濾波器的比較

    在本文方法中使用不同的保邊濾波器進(jìn)行實驗,選取引導(dǎo)濾波器(Guided Filter,GF)和雙邊濾波器(Bilateral Filter,BF)與上文使用的WLS 濾波器進(jìn)行對比。雙邊濾波器參數(shù)設(shè)置為δr=2,4,8,δs=0.1,0.2,0.4,引導(dǎo)濾波器參數(shù)設(shè)置為r=2,4,8,ε=0.12,0.22,0.42。表4~表6 列出了使用不同濾波器的客觀定量指標(biāo),其中,加粗?jǐn)?shù)據(jù)表示最優(yōu)值??梢钥闯?,使用WLS 濾波器效果較使用其他濾波器更好。

    表4 本文方法不同濾波器下的多發(fā)性腦梗塞MR-T1/MR-T2圖像融合評價指標(biāo)Table 4 Evaluation indexes of MR-T1/MR-T2 image fusion of multiple cerebral infarction by the proposed method using different filters

    表5 本文方法不同濾波器下的腦弓形蟲CT/MR 圖像融合評價指標(biāo)Table 5 Evaluation indexes of CT/MR image fusion oftoxoplasma gondii by the proposed method using different filters

    表6 本文方法不同濾波器下的正常腦部CT/MR 圖像融合評價指標(biāo)Table 6 Evaluation indexes of CT/MR image fusion of normal brain by the proposed method using different filters

    4 結(jié)束語

    傳統(tǒng)基于多尺度分析的圖像融合方法因無法保留邊緣特征而導(dǎo)致邊緣處產(chǎn)生光暈。針對該問題,本文提出一種多尺度邊緣保持分解與PCNN 結(jié)合的醫(yī)學(xué)圖像融合方法。在分解過程保護(hù)邊緣的同時進(jìn)行多尺度多方向分解,并采用PCNN 作為融合規(guī)則,從而改善融合圖像的視覺感知效果。實驗結(jié)果表明,該方法能夠突出醫(yī)學(xué)圖像的邊緣輪廓并增強(qiáng)圖像細(xì)節(jié),有利于將更多的顯著特征從源圖像分離并轉(zhuǎn)移到融合圖像中。本文方法使用了大量參數(shù),而參數(shù)設(shè)置將影響融合結(jié)果的質(zhì)量,因此,下一步將從參數(shù)自適應(yīng)角度著手對其進(jìn)行優(yōu)化。

    猜你喜歡
    濾波器尺度邊緣
    基于無擾濾波器和AED-ADT的無擾切換控制
    財產(chǎn)的五大尺度和五重應(yīng)對
    從濾波器理解卷積
    電子制作(2019年11期)2019-07-04 00:34:38
    開關(guān)電源EMI濾波器的應(yīng)用方法探討
    電子制作(2018年16期)2018-09-26 03:26:50
    一張圖看懂邊緣計算
    宇宙的尺度
    太空探索(2016年5期)2016-07-12 15:17:55
    基于TMS320C6678的SAR方位向預(yù)濾波器的并行實現(xiàn)
    9
    室外雕塑的尺度
    雕塑(2000年2期)2000-06-22 16:13:30
    在邊緣尋找自我
    雕塑(1999年2期)1999-06-28 05:01:42
    夫妻午夜视频| 99久久精品国产国产毛片| 亚洲va在线va天堂va国产| 天天躁日日操中文字幕| 夫妻午夜视频| 少妇被粗大猛烈的视频| 亚洲av日韩在线播放| 国产黄片美女视频| 亚洲国产精品国产精品| 亚洲精品一区蜜桃| 国产探花在线观看一区二区| 国产色婷婷99| 日本一本二区三区精品| 亚洲aⅴ乱码一区二区在线播放| 男人舔女人下体高潮全视频| 欧美最新免费一区二区三区| 日韩视频在线欧美| 天堂av国产一区二区熟女人妻| 伦精品一区二区三区| 国产男女超爽视频在线观看| 91精品一卡2卡3卡4卡| 成人鲁丝片一二三区免费| 精品久久久精品久久久| 久久午夜福利片| 超碰av人人做人人爽久久| 婷婷色综合www| 99热这里只有是精品在线观看| 一本一本综合久久| 亚洲av成人av| 久久久亚洲精品成人影院| 高清午夜精品一区二区三区| 亚洲经典国产精华液单| 久久久成人免费电影| 高清av免费在线| 亚洲av男天堂| 成人毛片60女人毛片免费| 日韩成人伦理影院| 国产欧美日韩精品一区二区| 亚洲丝袜综合中文字幕| 欧美xxⅹ黑人| 蜜臀久久99精品久久宅男| 亚洲欧美精品自产自拍| 亚洲在线自拍视频| 成人一区二区视频在线观看| 久久久久久久国产电影| 超碰97精品在线观看| 久久99热6这里只有精品| 亚洲欧美中文字幕日韩二区| 国产精品女同一区二区软件| 欧美日韩精品成人综合77777| 婷婷色综合大香蕉| 欧美日韩在线观看h| 亚洲成人中文字幕在线播放| 日日干狠狠操夜夜爽| 国产在视频线精品| 亚洲精品久久午夜乱码| 国产精品三级大全| 亚洲欧洲国产日韩| 国内揄拍国产精品人妻在线| 2021天堂中文幕一二区在线观| 日韩强制内射视频| 一个人观看的视频www高清免费观看| 亚洲国产欧美人成| 少妇丰满av| 亚洲真实伦在线观看| 汤姆久久久久久久影院中文字幕 | 少妇的逼好多水| 日本免费a在线| 超碰av人人做人人爽久久| 美女黄网站色视频| 黄片无遮挡物在线观看| 1000部很黄的大片| 天堂中文最新版在线下载 | 黄片无遮挡物在线观看| 日韩精品青青久久久久久| 日韩成人伦理影院| 国产在视频线在精品| 国产在线一区二区三区精| 久热久热在线精品观看| 免费看不卡的av| 国产免费一级a男人的天堂| 亚洲欧美精品自产自拍| 国产精品.久久久| 色5月婷婷丁香| 国产精品人妻久久久久久| 国产精品不卡视频一区二区| 麻豆成人av视频| 日韩av不卡免费在线播放| 国产成年人精品一区二区| 国产视频首页在线观看| 亚洲精品日韩av片在线观看| 成人毛片a级毛片在线播放| 高清午夜精品一区二区三区| 麻豆久久精品国产亚洲av| 2021天堂中文幕一二区在线观| 超碰av人人做人人爽久久| 99热6这里只有精品| 欧美3d第一页| 久99久视频精品免费| 免费观看的影片在线观看| 成人无遮挡网站| 日本欧美国产在线视频| 亚洲精品,欧美精品| 久久精品人妻少妇| 只有这里有精品99| 深夜a级毛片| 亚洲欧美日韩无卡精品| 国内少妇人妻偷人精品xxx网站| 亚洲电影在线观看av| 两个人的视频大全免费| 黄片无遮挡物在线观看| 97热精品久久久久久| 免费黄频网站在线观看国产| 2022亚洲国产成人精品| 亚洲国产精品成人综合色| 午夜激情欧美在线| 亚洲av不卡在线观看| 午夜老司机福利剧场| 久久99热这里只频精品6学生| 一夜夜www| 成人无遮挡网站| 九九爱精品视频在线观看| 精品酒店卫生间| 午夜日本视频在线| 亚洲国产精品sss在线观看| 国产免费一级a男人的天堂| 欧美最新免费一区二区三区| 在线天堂最新版资源| 亚洲电影在线观看av| 久久精品国产自在天天线| 一区二区三区免费毛片| 看黄色毛片网站| 深夜a级毛片| 国产av国产精品国产| 久久久久性生活片| 成年免费大片在线观看| 在现免费观看毛片| 18禁裸乳无遮挡免费网站照片| 春色校园在线视频观看| 在线免费十八禁| 日韩欧美国产在线观看| 男女边摸边吃奶| 七月丁香在线播放| 成人综合一区亚洲| 亚洲av中文字字幕乱码综合| 看黄色毛片网站| 大又大粗又爽又黄少妇毛片口| 97超碰精品成人国产| 亚洲欧美日韩东京热| 国产老妇女一区| 日本色播在线视频| 精品欧美国产一区二区三| 18+在线观看网站| 欧美bdsm另类| 男插女下体视频免费在线播放| 亚洲人成网站高清观看| 免费av不卡在线播放| freevideosex欧美| 国产精品福利在线免费观看| 日韩欧美精品v在线| 国产不卡一卡二| 久久精品综合一区二区三区| 亚洲精品国产av成人精品| 久久久久久久亚洲中文字幕| 免费观看av网站的网址| 久久久久精品久久久久真实原创| 美女国产视频在线观看| 亚洲内射少妇av| 淫秽高清视频在线观看| 午夜老司机福利剧场| 美女高潮的动态| 精品久久久久久久人妻蜜臀av| 国产爱豆传媒在线观看| 欧美成人a在线观看| 欧美97在线视频| av在线观看视频网站免费| 国产成人a∨麻豆精品| 成人综合一区亚洲| 91av网一区二区| 亚洲婷婷狠狠爱综合网| 国产在线一区二区三区精| 国产成人精品福利久久| 国产有黄有色有爽视频| 亚洲av国产av综合av卡| 亚洲av成人精品一二三区| 亚洲最大成人av| 欧美成人a在线观看| 国产精品伦人一区二区| 成人一区二区视频在线观看| 岛国毛片在线播放| 亚洲欧美一区二区三区黑人 | 91精品一卡2卡3卡4卡| 国产熟女欧美一区二区| 不卡视频在线观看欧美| 十八禁网站网址无遮挡 | 成人亚洲精品一区在线观看 | 老司机影院毛片| 精品久久久久久久久av| 日韩亚洲欧美综合| 中文字幕制服av| 午夜精品国产一区二区电影 | 免费看美女性在线毛片视频| 国产精品久久久久久av不卡| 久久久精品免费免费高清| 91av网一区二区| 国内精品宾馆在线| 啦啦啦韩国在线观看视频| 亚洲成人av在线免费| 国内少妇人妻偷人精品xxx网站| 亚洲av国产av综合av卡| av国产免费在线观看| 一边亲一边摸免费视频| 插阴视频在线观看视频| a级毛色黄片| 久久这里有精品视频免费| 成人毛片60女人毛片免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 乱系列少妇在线播放| 一级黄片播放器| 亚洲美女搞黄在线观看| .国产精品久久| av又黄又爽大尺度在线免费看| 男人和女人高潮做爰伦理| 免费在线观看成人毛片| 亚洲不卡免费看| 国产成人精品一,二区| 能在线免费观看的黄片| 99久久精品国产国产毛片| 欧美成人精品欧美一级黄| 国产日韩欧美在线精品| 国产探花极品一区二区| 欧美区成人在线视频| 久久久久精品久久久久真实原创| 国产精品精品国产色婷婷| 非洲黑人性xxxx精品又粗又长| 女人久久www免费人成看片| 国产欧美另类精品又又久久亚洲欧美| 亚洲综合精品二区| 午夜爱爱视频在线播放| 高清毛片免费看| 亚洲丝袜综合中文字幕| 在线天堂最新版资源| 欧美日韩精品成人综合77777| 深爱激情五月婷婷| 国产色爽女视频免费观看| 91aial.com中文字幕在线观看| 最近手机中文字幕大全| 亚洲经典国产精华液单| 免费观看性生交大片5| 久久精品国产鲁丝片午夜精品| 亚洲经典国产精华液单| 久久久久久久久大av| 一级毛片电影观看| 亚洲欧美清纯卡通| 建设人人有责人人尽责人人享有的 | 高清欧美精品videossex| 久久久久久久久久成人| 狂野欧美激情性xxxx在线观看| 国产 一区 欧美 日韩| 国产精品福利在线免费观看| 国产中年淑女户外野战色| 内射极品少妇av片p| 又大又黄又爽视频免费| 一区二区三区四区激情视频| 午夜福利视频1000在线观看| 成人二区视频| 成人综合一区亚洲| 精品国产露脸久久av麻豆 | 日日干狠狠操夜夜爽| 80岁老熟妇乱子伦牲交| 日韩大片免费观看网站| 天堂俺去俺来也www色官网 | 国产老妇女一区| 国产极品天堂在线| 免费不卡的大黄色大毛片视频在线观看 | 丝袜美腿在线中文| 日本黄色片子视频| 99久久九九国产精品国产免费| 男女国产视频网站| 国产激情偷乱视频一区二区| 国产成年人精品一区二区| 蜜桃久久精品国产亚洲av| 国产黄频视频在线观看| 日韩av在线大香蕉| 尤物成人国产欧美一区二区三区| 欧美一级a爱片免费观看看| 一个人看视频在线观看www免费| 亚洲精品aⅴ在线观看| 国产精品人妻久久久影院| 最后的刺客免费高清国语| 亚洲精品乱码久久久v下载方式| av在线观看视频网站免费| 国语对白做爰xxxⅹ性视频网站| 亚洲在线自拍视频| 欧美xxⅹ黑人| 深夜a级毛片| 久久久久久久久久黄片| 最近中文字幕高清免费大全6| 建设人人有责人人尽责人人享有的 | 日韩人妻高清精品专区| 国模一区二区三区四区视频| 波野结衣二区三区在线| 最近中文字幕高清免费大全6| 丝瓜视频免费看黄片| 国产男女超爽视频在线观看| 自拍偷自拍亚洲精品老妇| 亚洲内射少妇av| av国产久精品久网站免费入址| 久久精品久久久久久久性| 免费看不卡的av| 国内揄拍国产精品人妻在线| 欧美日韩亚洲高清精品| 国产三级在线视频| 久久午夜福利片| av线在线观看网站| 国产亚洲精品av在线| 亚洲成人av在线免费| 91av网一区二区| videossex国产| 成人二区视频| 国产成人精品福利久久| av在线播放精品| 国产亚洲精品av在线| 一级a做视频免费观看| 丝袜喷水一区| 爱豆传媒免费全集在线观看| 十八禁国产超污无遮挡网站| av免费观看日本| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站在线播| 十八禁网站网址无遮挡 | 午夜福利高清视频| 舔av片在线| 97精品久久久久久久久久精品| 欧美人与善性xxx| 丰满乱子伦码专区| 天堂av国产一区二区熟女人妻| 国产成人精品婷婷| av黄色大香蕉| 男人舔奶头视频| av在线播放精品| 91精品国产九色| av在线亚洲专区| av网站免费在线观看视频 | 国产免费视频播放在线视频 | 人人妻人人澡欧美一区二区| 精品人妻熟女av久视频| 日韩一本色道免费dvd| 男人爽女人下面视频在线观看| 成年女人看的毛片在线观看| 在现免费观看毛片| 永久网站在线| 一级毛片久久久久久久久女| 国产在视频线在精品| 国产高清不卡午夜福利| 亚洲av男天堂| 男女视频在线观看网站免费| 亚洲人成网站在线播| 日本免费在线观看一区| 国产成人福利小说| 免费黄网站久久成人精品| 五月伊人婷婷丁香| 欧美成人精品欧美一级黄| 国产黄色免费在线视频| 午夜福利在线观看吧| 欧美潮喷喷水| 日本-黄色视频高清免费观看| 亚洲国产av新网站| 国产黄片视频在线免费观看| 亚洲精品国产av蜜桃| 寂寞人妻少妇视频99o| 淫秽高清视频在线观看| 高清毛片免费看| 久久久午夜欧美精品| 一级黄片播放器| 国内少妇人妻偷人精品xxx网站| 亚洲欧美中文字幕日韩二区| 婷婷色麻豆天堂久久| 亚洲欧美中文字幕日韩二区| 一个人看的www免费观看视频| 97在线视频观看| 超碰97精品在线观看| 亚洲成色77777| 美女黄网站色视频| 18禁在线播放成人免费| 国产精品国产三级国产专区5o| 大香蕉久久网| 亚州av有码| 久久99热6这里只有精品| 免费观看av网站的网址| 天天躁日日操中文字幕| a级一级毛片免费在线观看| 免费观看a级毛片全部| 两个人视频免费观看高清| 午夜福利成人在线免费观看| 国产男女超爽视频在线观看| xxx大片免费视频| 成人av在线播放网站| 欧美+日韩+精品| 一边亲一边摸免费视频| 免费观看精品视频网站| 亚洲欧洲日产国产| 精品国内亚洲2022精品成人| 菩萨蛮人人尽说江南好唐韦庄| 在线免费观看不下载黄p国产| 熟妇人妻久久中文字幕3abv| 热99在线观看视频| 大片免费播放器 马上看| 日韩视频在线欧美| 亚洲精品久久久久久婷婷小说| 大片免费播放器 马上看| 可以在线观看毛片的网站| 婷婷色综合大香蕉| 亚洲成人久久爱视频| 国产乱人视频| 一级av片app| 男女下面进入的视频免费午夜| 高清毛片免费看| 国产成人免费观看mmmm| 亚洲在线自拍视频| 国产亚洲午夜精品一区二区久久 | 国产在线一区二区三区精| 三级国产精品欧美在线观看| 综合色丁香网| 极品教师在线视频| 在线a可以看的网站| 麻豆国产97在线/欧美| 免费看av在线观看网站| 午夜福利网站1000一区二区三区| 欧美成人一区二区免费高清观看| 欧美xxⅹ黑人| 国产精品嫩草影院av在线观看| 成人亚洲精品一区在线观看 | 青青草视频在线视频观看| 女的被弄到高潮叫床怎么办| 国产精品综合久久久久久久免费| 一二三四中文在线观看免费高清| 国产淫语在线视频| 久久久成人免费电影| 丝袜喷水一区| 精品一区在线观看国产| 欧美人与善性xxx| 久久韩国三级中文字幕| 插逼视频在线观看| 亚洲国产精品成人久久小说| 永久网站在线| 少妇猛男粗大的猛烈进出视频 | 99热全是精品| 亚洲经典国产精华液单| 国产精品爽爽va在线观看网站| 精品亚洲乱码少妇综合久久| 国产三级在线视频| 午夜爱爱视频在线播放| 日韩一本色道免费dvd| 超碰97精品在线观看| 热99在线观看视频| 观看免费一级毛片| 男女下面进入的视频免费午夜| 真实男女啪啪啪动态图| 国产av码专区亚洲av| 大香蕉97超碰在线| 午夜日本视频在线| 99热这里只有精品一区| 免费av毛片视频| 亚洲精品国产av蜜桃| 国产三级在线视频| 男人爽女人下面视频在线观看| 熟妇人妻久久中文字幕3abv| 一级片'在线观看视频| 中文字幕亚洲精品专区| 国产乱人视频| 国语对白做爰xxxⅹ性视频网站| 国内精品一区二区在线观看| 男女国产视频网站| 欧美高清成人免费视频www| 一级a做视频免费观看| 午夜精品一区二区三区免费看| 99热网站在线观看| 如何舔出高潮| 国产女主播在线喷水免费视频网站 | 国产亚洲精品久久久com| 你懂的网址亚洲精品在线观看| 能在线免费观看的黄片| 婷婷色综合大香蕉| 天堂网av新在线| 99九九线精品视频在线观看视频| 夫妻性生交免费视频一级片| 乱系列少妇在线播放| 精华霜和精华液先用哪个| 别揉我奶头 嗯啊视频| 久久精品久久久久久噜噜老黄| 欧美一区二区亚洲| 亚洲成人一二三区av| 边亲边吃奶的免费视频| 国产老妇伦熟女老妇高清| 麻豆久久精品国产亚洲av| 禁无遮挡网站| 大陆偷拍与自拍| 亚洲精品影视一区二区三区av| 熟女电影av网| 国产乱来视频区| 久久99热这里只有精品18| 日韩一区二区视频免费看| 最近中文字幕高清免费大全6| 国产av码专区亚洲av| 非洲黑人性xxxx精品又粗又长| 又大又黄又爽视频免费| 国产欧美日韩精品一区二区| 69av精品久久久久久| 国产女主播在线喷水免费视频网站 | 搞女人的毛片| 亚洲av免费高清在线观看| 亚洲乱码一区二区免费版| 日韩欧美国产在线观看| 国产精品精品国产色婷婷| 色网站视频免费| 国产精品久久久久久精品电影小说 | 国产精品99久久久久久久久| 禁无遮挡网站| 亚洲熟女精品中文字幕| 日韩视频在线欧美| 69人妻影院| 麻豆久久精品国产亚洲av| 亚洲精品456在线播放app| 亚洲18禁久久av| 日韩国内少妇激情av| 男人和女人高潮做爰伦理| 亚洲熟女精品中文字幕| 黄色欧美视频在线观看| 国产高清有码在线观看视频| 三级国产精品欧美在线观看| 国产精品久久久久久久久免| 高清av免费在线| 国产成人精品福利久久| 水蜜桃什么品种好| 九九久久精品国产亚洲av麻豆| 赤兔流量卡办理| 国产精品久久视频播放| 日产精品乱码卡一卡2卡三| 极品少妇高潮喷水抽搐| 国产精品精品国产色婷婷| 全区人妻精品视频| 男人狂女人下面高潮的视频| 又粗又硬又长又爽又黄的视频| 人妻制服诱惑在线中文字幕| 午夜爱爱视频在线播放| 亚洲精品日韩在线中文字幕| 少妇熟女欧美另类| 亚洲乱码一区二区免费版| 亚洲精品中文字幕在线视频 | 国产黄色免费在线视频| av免费观看日本| 啦啦啦韩国在线观看视频| 亚洲国产精品成人综合色| 日韩av不卡免费在线播放| av线在线观看网站| 国产欧美另类精品又又久久亚洲欧美| 亚洲av成人精品一区久久| 国产精品一区二区三区四区久久| 国产高清三级在线| 久久99热这里只有精品18| 高清在线视频一区二区三区| 亚洲精品乱久久久久久| 日韩国内少妇激情av| 亚洲国产色片| 婷婷六月久久综合丁香| 久久亚洲国产成人精品v| 免费观看的影片在线观看| 日本午夜av视频| 久久这里有精品视频免费| 好男人视频免费观看在线| freevideosex欧美| 丝袜喷水一区| 两个人视频免费观看高清| 精品人妻偷拍中文字幕| av黄色大香蕉| 最近视频中文字幕2019在线8| 熟妇人妻久久中文字幕3abv| 69av精品久久久久久| h日本视频在线播放| 久久午夜福利片| 美女脱内裤让男人舔精品视频| 丰满少妇做爰视频| 日韩av不卡免费在线播放| 日韩欧美国产在线观看| 一个人免费在线观看电影| 国产av码专区亚洲av| 能在线免费观看的黄片| 亚洲欧美日韩无卡精品| 欧美97在线视频| 大片免费播放器 马上看| 22中文网久久字幕| 蜜桃亚洲精品一区二区三区| 一级毛片 在线播放| 国产成人一区二区在线| 在线免费十八禁| 九九久久精品国产亚洲av麻豆| 久久久久久久亚洲中文字幕| 麻豆国产97在线/欧美| 中文字幕av在线有码专区| 成人漫画全彩无遮挡| 欧美人与善性xxx| 日韩av在线大香蕉| 纵有疾风起免费观看全集完整版 | 街头女战士在线观看网站| 联通29元200g的流量卡| 97热精品久久久久久| 亚洲成人中文字幕在线播放| 看黄色毛片网站| 国产精品熟女久久久久浪| 日本wwww免费看| 内射极品少妇av片p| 狂野欧美白嫩少妇大欣赏| 亚州av有码| 人人妻人人澡人人爽人人夜夜 |