• <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
    欧美黑人欧美精品刺激| 亚洲情色 制服丝袜| 不卡av一区二区三区| 成人av一区二区三区在线看| 亚洲久久久国产精品| 一边摸一边抽搐一进一小说 | 久热爱精品视频在线9| 国产成人系列免费观看| 90打野战视频偷拍视频| 国产成人一区二区三区免费视频网站| 一区二区三区国产精品乱码| 国产深夜福利视频在线观看| 亚洲熟妇熟女久久| 18禁黄网站禁片午夜丰满| 黄色怎么调成土黄色| 久久香蕉精品热| 久久精品aⅴ一区二区三区四区| 性色av乱码一区二区三区2| 日本撒尿小便嘘嘘汇集6| 一级毛片精品| 国产精品影院久久| 黑丝袜美女国产一区| 久久九九热精品免费| 高清在线国产一区| 9热在线视频观看99| 亚洲欧美色中文字幕在线| 老司机影院毛片| av福利片在线| 在线国产一区二区在线| 国产在线一区二区三区精| 国产一区二区激情短视频| 久久九九热精品免费| 两个人看的免费小视频| 少妇粗大呻吟视频| 久久影院123| 精品国产乱子伦一区二区三区| 一级毛片精品| 国产精品98久久久久久宅男小说| 建设人人有责人人尽责人人享有的| 亚洲一区中文字幕在线| 成人特级黄色片久久久久久久| 三上悠亚av全集在线观看| 一边摸一边做爽爽视频免费| 超碰成人久久| 成人18禁高潮啪啪吃奶动态图| 久久影院123| 午夜日韩欧美国产| 高清欧美精品videossex| 国产乱人伦免费视频| 999精品在线视频| 婷婷精品国产亚洲av在线 | 欧美成人免费av一区二区三区 | 一进一出抽搐gif免费好疼 | 宅男免费午夜| 露出奶头的视频| a级毛片在线看网站| 久久久国产成人免费| 亚洲第一欧美日韩一区二区三区| 可以免费在线观看a视频的电影网站| 亚洲av成人不卡在线观看播放网| 精品电影一区二区在线| 一边摸一边做爽爽视频免费| 亚洲国产精品合色在线| 国产精品乱码一区二三区的特点 | 丝袜美腿诱惑在线| 夜夜爽天天搞| 一级作爱视频免费观看| 色综合欧美亚洲国产小说| 91大片在线观看| 黑人操中国人逼视频| 久热爱精品视频在线9| 欧美+亚洲+日韩+国产| 在线永久观看黄色视频| 露出奶头的视频| 一级,二级,三级黄色视频| x7x7x7水蜜桃| 欧美久久黑人一区二区| 亚洲欧美一区二区三区久久| 妹子高潮喷水视频| 精品一区二区三卡| 色94色欧美一区二区| 久久精品91无色码中文字幕| 亚洲aⅴ乱码一区二区在线播放 | 动漫黄色视频在线观看| 午夜影院日韩av| 国产一卡二卡三卡精品| 丝袜美腿诱惑在线| 99香蕉大伊视频| 老司机午夜福利在线观看视频| 在线av久久热| 男女高潮啪啪啪动态图| 久久久国产成人精品二区 | www.精华液| 99久久99久久久精品蜜桃| 久热爱精品视频在线9| 乱人伦中国视频| 人妻 亚洲 视频| 岛国毛片在线播放| 欧美精品人与动牲交sv欧美| 99久久国产精品久久久| 午夜成年电影在线免费观看| 国产99久久九九免费精品| 欧美精品啪啪一区二区三区| 激情视频va一区二区三区| 免费少妇av软件| www.999成人在线观看| 亚洲黑人精品在线| 看免费av毛片| 久久国产乱子伦精品免费另类| 狠狠婷婷综合久久久久久88av| 可以免费在线观看a视频的电影网站| 亚洲一区中文字幕在线| 99在线人妻在线中文字幕 | 日韩欧美一区视频在线观看| 男人操女人黄网站| 一本综合久久免费| 精品国产美女av久久久久小说| 欧美亚洲 丝袜 人妻 在线| 国产成人av激情在线播放| 中文字幕制服av| 国产精品一区二区在线观看99| 国产成人欧美在线观看 | 波多野结衣一区麻豆| 中文字幕av电影在线播放| 十八禁高潮呻吟视频| 在线看a的网站| 免费av中文字幕在线| 黄网站色视频无遮挡免费观看| 91精品三级在线观看| 久久这里只有精品19| 757午夜福利合集在线观看| 日韩欧美一区视频在线观看| 无人区码免费观看不卡| 国产精品免费视频内射| 91成年电影在线观看| 中出人妻视频一区二区| 亚洲精品成人av观看孕妇| www.999成人在线观看| 国产麻豆69| 一本综合久久免费| 精品亚洲成a人片在线观看| 在线天堂中文资源库| 欧美日韩乱码在线| 午夜免费观看网址| 另类亚洲欧美激情| 久久久久久久午夜电影 | av有码第一页| 久久久久久久久久久久大奶| 日本撒尿小便嘘嘘汇集6| 亚洲情色 制服丝袜| 老司机在亚洲福利影院| 精品人妻1区二区| 亚洲视频免费观看视频| 午夜亚洲福利在线播放| 成人黄色视频免费在线看| 国产精品 欧美亚洲| 亚洲中文日韩欧美视频| 亚洲精品中文字幕在线视频| 国产主播在线观看一区二区| 中文字幕人妻丝袜制服| 国产单亲对白刺激| 欧美久久黑人一区二区| 黄色a级毛片大全视频| 好看av亚洲va欧美ⅴa在| 亚洲午夜精品一区,二区,三区| 又黄又爽又免费观看的视频| 久久久久久免费高清国产稀缺| 国产亚洲精品第一综合不卡| 天天添夜夜摸| 丝袜人妻中文字幕| 欧美黄色淫秽网站| 国产99久久九九免费精品| 亚洲欧美日韩另类电影网站| 村上凉子中文字幕在线| 国产免费现黄频在线看| tocl精华| 50天的宝宝边吃奶边哭怎么回事| 少妇 在线观看| 亚洲国产看品久久| 精品免费久久久久久久清纯 | 国产午夜精品久久久久久| 精品欧美一区二区三区在线| 久久久久久人人人人人| 亚洲五月天丁香| 精品少妇久久久久久888优播| 新久久久久国产一级毛片| av超薄肉色丝袜交足视频| 免费观看精品视频网站| 在线观看免费视频网站a站| 久久久国产成人免费| 99国产精品一区二区三区| 国产精品免费一区二区三区在线 | 国产精品久久久久久精品古装| 天天添夜夜摸| 脱女人内裤的视频| 两个人免费观看高清视频| 欧美日韩福利视频一区二区| 国产精品98久久久久久宅男小说| 丝袜美腿诱惑在线| svipshipincom国产片| 中文字幕高清在线视频| 欧美乱妇无乱码| 免费不卡黄色视频| 亚洲人成伊人成综合网2020| 久久九九热精品免费| 一区二区三区精品91| 两性夫妻黄色片| 国产精品偷伦视频观看了| 两性夫妻黄色片| 91国产中文字幕| 黄频高清免费视频| 最近最新中文字幕大全免费视频| 天堂√8在线中文| 99精品欧美一区二区三区四区| 成年动漫av网址| 久久狼人影院| 天天影视国产精品| 国产在线观看jvid| 国产一区有黄有色的免费视频| 免费看a级黄色片| 国产日韩欧美亚洲二区| 亚洲精品在线观看二区| 一级a爱片免费观看的视频| 嫩草影视91久久| 精品视频人人做人人爽| 一边摸一边抽搐一进一小说 | 国产深夜福利视频在线观看| 99国产综合亚洲精品| 精品一区二区三区四区五区乱码| 亚洲 欧美一区二区三区| 免费日韩欧美在线观看| 啦啦啦视频在线资源免费观看| 91九色精品人成在线观看| 色在线成人网| 天天添夜夜摸| av网站在线播放免费| 亚洲av成人一区二区三| 国产aⅴ精品一区二区三区波| 国产视频一区二区在线看| 精品国产一区二区三区四区第35| 亚洲精品粉嫩美女一区| 黄片大片在线免费观看| 大片电影免费在线观看免费| 美女福利国产在线| 一区二区三区激情视频| 亚洲中文字幕日韩| 亚洲精品国产区一区二| 变态另类成人亚洲欧美熟女 | 国产精品自产拍在线观看55亚洲 | 国产人伦9x9x在线观看| 成人18禁在线播放| 日日爽夜夜爽网站| 一级黄色大片毛片| 妹子高潮喷水视频| 国产精品二区激情视频| 91在线观看av| 免费在线观看影片大全网站| 国产午夜精品久久久久久| www.999成人在线观看| 欧美乱妇无乱码| 免费在线观看黄色视频的| 两人在一起打扑克的视频| 免费久久久久久久精品成人欧美视频| 美女国产高潮福利片在线看| 日本欧美视频一区| cao死你这个sao货| www.熟女人妻精品国产| 国产无遮挡羞羞视频在线观看| 亚洲av第一区精品v没综合| 免费观看精品视频网站| 一本一本久久a久久精品综合妖精| 欧美黄色淫秽网站| av天堂在线播放| 中文欧美无线码| 久久久久国产精品人妻aⅴ院 | av不卡在线播放| 人妻一区二区av| 成熟少妇高潮喷水视频| 午夜福利视频在线观看免费| 99热国产这里只有精品6| 在线观看66精品国产| 久久人妻熟女aⅴ| 看黄色毛片网站| 欧美黑人欧美精品刺激| 亚洲情色 制服丝袜| 欧美黑人精品巨大| 亚洲国产中文字幕在线视频| 国产淫语在线视频| 欧美日韩亚洲国产一区二区在线观看 | 国产xxxxx性猛交| 黄色片一级片一级黄色片| 亚洲成人免费电影在线观看| 色精品久久人妻99蜜桃| 国产精品1区2区在线观看. | 欧美 亚洲 国产 日韩一| 国产三级黄色录像| 国产成人免费无遮挡视频| a级片在线免费高清观看视频| 久久久国产一区二区| 午夜影院日韩av| 99热国产这里只有精品6| 动漫黄色视频在线观看| 久久性视频一级片| 久久久久久免费高清国产稀缺| 久久精品91无色码中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 久久人妻熟女aⅴ| 王馨瑶露胸无遮挡在线观看| 一个人免费在线观看的高清视频| 成人精品一区二区免费| 久久精品人人爽人人爽视色| 中文字幕av电影在线播放| 男女床上黄色一级片免费看| 我的亚洲天堂| 亚洲av成人av| av一本久久久久| 日本欧美视频一区| 成熟少妇高潮喷水视频| 巨乳人妻的诱惑在线观看| 99热网站在线观看| 亚洲美女黄片视频| 久久久国产成人精品二区 | 侵犯人妻中文字幕一二三四区| 国产精品免费一区二区三区在线 | 午夜精品在线福利| 成人黄色视频免费在线看| 女人被狂操c到高潮| 怎么达到女性高潮| 激情在线观看视频在线高清 | av超薄肉色丝袜交足视频| 桃红色精品国产亚洲av| 三级毛片av免费| 欧美成人免费av一区二区三区 | 欧美乱码精品一区二区三区| 色老头精品视频在线观看| 免费看a级黄色片| 他把我摸到了高潮在线观看| 亚洲精品美女久久久久99蜜臀| 叶爱在线成人免费视频播放| 欧美成狂野欧美在线观看| 国产单亲对白刺激| 又黄又粗又硬又大视频| 高清视频免费观看一区二区| 日本vs欧美在线观看视频| 国产精品98久久久久久宅男小说| 亚洲精品美女久久av网站| 亚洲精品乱久久久久久| 国产精品98久久久久久宅男小说| 久久精品国产清高在天天线| 桃红色精品国产亚洲av| 欧美久久黑人一区二区| 亚洲专区中文字幕在线| 天堂动漫精品| 久久中文字幕一级| 黄片大片在线免费观看| 久久久久久亚洲精品国产蜜桃av| 精品国产乱子伦一区二区三区| 国产精品99久久99久久久不卡| 老汉色∧v一级毛片| 亚洲第一av免费看| 午夜福利免费观看在线| 亚洲欧美一区二区三区黑人| 在线av久久热| 色综合婷婷激情| 久久狼人影院| 91av网站免费观看| 中文欧美无线码| 久久青草综合色| 欧美av亚洲av综合av国产av| 男人的好看免费观看在线视频 | 老司机影院毛片| 亚洲熟女精品中文字幕| 亚洲在线自拍视频| 中文字幕高清在线视频| 亚洲欧美精品综合一区二区三区| 欧美日韩亚洲高清精品| 日韩欧美一区视频在线观看| 人人妻人人澡人人看| av不卡在线播放| 午夜福利,免费看| 女警被强在线播放| 免费日韩欧美在线观看| 天堂动漫精品| 777米奇影视久久| 999精品在线视频| 91av网站免费观看| 欧美大码av| 99国产极品粉嫩在线观看| 一区二区三区精品91| 无限看片的www在线观看| 男女床上黄色一级片免费看| 99国产精品一区二区蜜桃av | 1024香蕉在线观看| 在线观看日韩欧美| 国产精品亚洲一级av第二区| 日韩视频一区二区在线观看| 伦理电影免费视频| 巨乳人妻的诱惑在线观看| 黄色片一级片一级黄色片| 国产亚洲一区二区精品| 一夜夜www| 久久国产亚洲av麻豆专区| www.999成人在线观看| svipshipincom国产片| 日本五十路高清| 成年人黄色毛片网站| www日本在线高清视频| 久久久久久久午夜电影 | 精品久久久久久电影网| 欧美性长视频在线观看| 宅男免费午夜| 美女扒开内裤让男人捅视频| av欧美777| 亚洲国产中文字幕在线视频| 精品乱码久久久久久99久播| xxxhd国产人妻xxx| 人人妻人人澡人人爽人人夜夜| 女人爽到高潮嗷嗷叫在线视频| 极品教师在线免费播放| av不卡在线播放| 亚洲av日韩在线播放| 成年版毛片免费区| 成人手机av| 国产99白浆流出| 香蕉国产在线看| 欧美日韩亚洲高清精品| 日本wwww免费看| 女人高潮潮喷娇喘18禁视频| 人妻久久中文字幕网| 在线观看免费视频网站a站| 夫妻午夜视频| 捣出白浆h1v1| 成人精品一区二区免费| 黄色视频,在线免费观看| 国产免费男女视频| 亚洲五月天丁香| 亚洲熟女毛片儿| 青草久久国产| 亚洲专区中文字幕在线| 成年女人毛片免费观看观看9 | 三级毛片av免费| 99久久人妻综合| 久久精品国产99精品国产亚洲性色 | 亚洲精品乱久久久久久| av欧美777| 精品第一国产精品| 亚洲精品国产精品久久久不卡| 看片在线看免费视频| 久久中文字幕人妻熟女| 手机成人av网站| 这个男人来自地球电影免费观看| 纯流量卡能插随身wifi吗| 精品无人区乱码1区二区| 亚洲伊人色综图| 精品欧美一区二区三区在线| 午夜福利一区二区在线看| 制服人妻中文乱码| 久久精品人人爽人人爽视色| 又黄又粗又硬又大视频| a级毛片在线看网站| 亚洲精品中文字幕在线视频| 99热国产这里只有精品6| 十八禁高潮呻吟视频| 男女之事视频高清在线观看| 欧美日韩福利视频一区二区| 亚洲中文字幕日韩| 香蕉丝袜av| svipshipincom国产片| 日韩免费av在线播放| 久久精品国产亚洲av高清一级| 成人影院久久| 十八禁高潮呻吟视频| 一级片免费观看大全| 亚洲午夜理论影院| 日本精品一区二区三区蜜桃| 免费av中文字幕在线| 在线观看一区二区三区激情| av网站在线播放免费| 无遮挡黄片免费观看| 热99re8久久精品国产| 成人特级黄色片久久久久久久| 成在线人永久免费视频| 亚洲色图综合在线观看| 亚洲视频免费观看视频| 99久久综合精品五月天人人| 久久人人97超碰香蕉20202| 黄片大片在线免费观看| 激情在线观看视频在线高清 | 亚洲av片天天在线观看| 久久久水蜜桃国产精品网| 欧美大码av| 80岁老熟妇乱子伦牲交| 99久久综合精品五月天人人| 亚洲少妇的诱惑av| 亚洲人成伊人成综合网2020| 欧美+亚洲+日韩+国产| 免费在线观看黄色视频的| 激情视频va一区二区三区| 这个男人来自地球电影免费观看| 97人妻天天添夜夜摸| 亚洲精品中文字幕一二三四区| 精品久久久精品久久久| 日韩三级视频一区二区三区| 黄片播放在线免费| 岛国毛片在线播放| 久久 成人 亚洲| 精品人妻熟女毛片av久久网站| 男女高潮啪啪啪动态图| 成年动漫av网址| 欧美中文综合在线视频| 欧美在线黄色| 可以免费在线观看a视频的电影网站| av天堂在线播放| 很黄的视频免费| ponron亚洲| 无人区码免费观看不卡| 久久久久久人人人人人| 波多野结衣一区麻豆| 国产一区二区三区视频了| 亚洲熟女精品中文字幕| 久久久久久免费高清国产稀缺| 亚洲五月婷婷丁香| 亚洲欧美激情在线| 18禁裸乳无遮挡免费网站照片 | a级毛片黄视频| 中文字幕人妻丝袜制服| 岛国在线观看网站| 少妇裸体淫交视频免费看高清 | 丁香欧美五月| 国产成人av教育| 精品电影一区二区在线| 久99久视频精品免费| 国产一区二区三区视频了| 国产成人精品无人区| 亚洲国产精品合色在线| 交换朋友夫妻互换小说| 变态另类成人亚洲欧美熟女 | 老司机靠b影院| 欧美日韩成人在线一区二区| 手机成人av网站| 国产精品久久视频播放| 国产高清视频在线播放一区| 欧美日韩一级在线毛片| 身体一侧抽搐| 国产极品粉嫩免费观看在线| 飞空精品影院首页| 午夜精品在线福利| 天天躁夜夜躁狠狠躁躁| 久久久国产成人精品二区 | 嫁个100分男人电影在线观看| 欧美日韩av久久| 欧美色视频一区免费| 黄色 视频免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 免费在线观看影片大全网站| 青草久久国产| 亚洲精品在线美女| 成人手机av| 丝瓜视频免费看黄片| av网站在线播放免费| 啦啦啦免费观看视频1| 精品久久蜜臀av无| 国产又爽黄色视频| 日韩欧美在线二视频 | 老司机福利观看| 久久久精品免费免费高清| 国产片内射在线| 久久中文看片网| 国产激情欧美一区二区| 成人免费观看视频高清| 一本大道久久a久久精品| 天堂中文最新版在线下载| 老司机在亚洲福利影院| 18禁裸乳无遮挡免费网站照片 | 不卡av一区二区三区| 亚洲国产精品一区二区三区在线| 欧美精品啪啪一区二区三区| 色精品久久人妻99蜜桃| 国产亚洲av高清不卡| 两性夫妻黄色片| 精品电影一区二区在线| 日本撒尿小便嘘嘘汇集6| 欧美在线黄色| 高清在线国产一区| 国产欧美日韩精品亚洲av| 成人精品一区二区免费| 国内毛片毛片毛片毛片毛片| 99国产精品一区二区三区| 一夜夜www| 欧美日韩瑟瑟在线播放| 亚洲黑人精品在线| 女同久久另类99精品国产91| 国产精品成人在线| 亚洲精品国产区一区二| 香蕉国产在线看| 好男人电影高清在线观看| 十分钟在线观看高清视频www| 欧美精品高潮呻吟av久久| 夜夜夜夜夜久久久久| 欧美在线黄色| 9色porny在线观看| 啦啦啦 在线观看视频| 91在线观看av| 男女下面插进去视频免费观看| 国产精品亚洲一级av第二区| 欧美日韩视频精品一区| av天堂在线播放| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产毛片av蜜桃av| 一级片'在线观看视频| 人妻一区二区av| 一级毛片女人18水好多| 大陆偷拍与自拍| 一本综合久久免费| 超色免费av|