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

    基于改進MC算法的腦圖譜三維重建

    2016-05-09 07:19:04樊令仲李海芳
    計算機應用與軟件 2016年4期
    關鍵詞:體素等值立方體

    朱 愷 樊令仲 李海芳*

    基于改進MC算法的腦圖譜三維重建

    朱 愷1樊令仲2,3李海芳1*

    1(太原理工大學計算機科學與技術學院 山西 太原 030024)

    2(中國科學院自動化研究所腦網(wǎng)絡組研究中心 北京 100190)

    3(中國科學院自動化研究所模式識別國家重點實驗室 北京 100190)

    隨著腦圖譜在神經(jīng)外科手術的教學與應用中起到越來越重要的作用,腦圖譜的三維重建成為了腦圖譜研究的一個主要方向。針對MarchingCubes算法在提取腦圖譜數(shù)據(jù)等值面效率低下的問題,通過對MarchingCubes算法的遍歷方法和拓撲結(jié)構的研究,提出一種適用于腦圖譜三維重建的改進算法,采用該改進算法實現(xiàn)腦圖譜三維重建系統(tǒng)。并以Brainnetome Atlas腦圖譜數(shù)據(jù)作為實驗數(shù)據(jù)對系統(tǒng)效率進行測試,實驗結(jié)果表明,相對于原始算法,改進算法確實提高了三維重建系統(tǒng)處理腦圖譜數(shù)據(jù)的效率。

    等值面 MarchingCubes 三維重建 腦分區(qū)圖譜

    0 引 言

    腦圖譜不僅是臨床神經(jīng)外科手術中的重要工具,而且在神經(jīng)解剖的教育教學中也扮演著重要的角色。早期的腦圖譜均為紙質(zhì)印刷的二維圖像。近年來,隨著大腦圖譜研究的不斷發(fā)展,人們開始尋求將腦圖譜進行數(shù)字化處理,并且從二維圖譜逐漸的向三維圖譜發(fā)展。通過將二維的腦部切片在計算機中重新整合并最終完成腦部結(jié)構的三維重建,從而使人們不僅僅能夠觀察大腦的不同斷面也能夠以三維的形式對腦結(jié)構進行觀察。

    在腦圖譜的三維重建系統(tǒng)中最關鍵的技術就是體數(shù)據(jù)的三維重建技術,體數(shù)據(jù)的三維重建算法直接決定了整個系統(tǒng)的重建效率以及重建準確率。醫(yī)學體數(shù)據(jù)三維可視化方法分為兩大類[1]:面繪制方法和體繪制方法。面繪制計算量小,處理速度快;體繪制能良好地保持醫(yī)學圖像的細節(jié),但是計算量比較大。因此針對于應用背景選擇一種合適的算法就成為了三維腦圖譜繪制的一個核心問題。

    基于對傳統(tǒng)體數(shù)據(jù)三維重建算法的分析,本文提出一種基于移動立方體算法的快速體數(shù)據(jù)三維重建算法,提高了三維重建的速度,并應用于三維腦圖譜繪制系統(tǒng)當中。

    1 三維重建算法

    1.1 重建算法的選擇

    三維重建算法主要分為體繪制和面繪制兩種,其中體繪制的常用算法中有很多基于光線投射算法[2]。例如等值面光線投射繪制、最大密度投影算法和合成體法光線投射算法。面繪制算法常用的有Marching Cubes算法。為了提高系統(tǒng)的運行效率就必須選擇一種合適的算法對腦圖譜數(shù)據(jù)進行三維重建。為此本文選擇了在腦圖譜三維重建中經(jīng)常用到的幾種文件規(guī)模對以上這些算法的效率進行了測試,測試所用文件列表如表1所示。

    表1 測試文件列表

    以上這三個文件中AIPS_IP1為單一腦區(qū)數(shù)據(jù)文件,cpbet為本文中使用的腦模板文件,cpbetter則為一種更為精細的腦模板文件,該文件并沒有在本文中使用,這里作為一種極端情況來測試算法效率。

    算法效率測試不對原始數(shù)據(jù)文件進行二次抽樣,也不對生成的圖形進行優(yōu)化。測試包含三個指標:(1) 數(shù)據(jù)處理的時間;(2) 處理完成后圖形渲染顯示的時間;(3) 所用時間的綜合。測試結(jié)果如表2所示。

    表2 算法效率測試結(jié)果

    通過表2我們可以看出在這些三維重建算法中,體繪制算法中的合成體法光線投射算法效果是最好的,但是其渲染效率卻沒有Marching Cubes算法高。在本文中用戶對生成的三維圖像每進行一次交互操作,都要進行一次渲染,因此渲染時間的長短是算法是否適合本系統(tǒng)的一個重要指標,過長的渲染時間將極大地影響用戶的體驗。此外體繪制在用戶進行視角變換操作時也要進行重新繪制,這大大影響了使用效果。對于移動立方體算法只要提高數(shù)據(jù)處理的效率就能使其整體效率達到最高。因此本文最后選用Marching Cubes作為系統(tǒng)的三維重建算法。

    1.2 MC算法

    Lorensond等人于1987年提出的Marching Cubes算法是一種流行的從體數(shù)據(jù)中提取等值面的算法[3]。其原理是通過對體數(shù)據(jù)中各個體素進行處理,找到含有等值面的體素并將這些體素中的等值面對應到一系列等值面穿過體素所在立方體的模式中。然后將這些對應好的模式以某種拓撲結(jié)構連接成三角面片得到等值面。體素中等值面的模式經(jīng)過反轉(zhuǎn)、旋轉(zhuǎn)等操作后有很多模式是重復的,去除重復的模式最終體素含有的等值面有15種模式,如圖1所示。在此圖中深色的頂點代表體素中該頂點處于等值面的閥值之下,也就是等值面在該點之上。根據(jù)立方體8個頂點不同的狀態(tài)決定了該體素中等值面所處的狀態(tài)。

    圖1 等值面在立方體中的15種模式

    得到等值面在體素內(nèi)的狀態(tài)后,就要根據(jù)等值面狀態(tài)來得到等值面與體素立方體各條棱相交的頂點,假設有一條棱與等值面相交于點c(x,y,z),等值面閾值為T。該棱的兩個頂點分別為a(x,y,z)和b(x,y,z),兩點的值分別為va與vb,那么交點c的位置可以由下式求出:

    (1)

    由此類推能將其他交點頂點的位置求出,計算出等值面各個頂點之后,還要計算等值面每個三角面片的法向量,為了計算該法向量就要先計算體素立方體頂點的梯度值,假設頂點梯度值為G(i,j,k),則梯度值能夠由下式求出:

    (2)

    其中D(i,j,k)為k層(i,j)點的灰度值,在本文應用背景下其實D就可以取為腦圖譜數(shù)據(jù)k層(i,j)點的值,而△x、△y、△z則分別取為xyz方向上的掃描間隔。

    MC算法作為一種經(jīng)典算法本身也存在一些缺點,例如等值面存在二義性、繪制精度不高、處理效率不高等。針對這些缺點Nielson等提出了解決二義性的方法[4];Lopes提出了提高繪制精度的方法[5];張迎平等提出的區(qū)間樹硬件加速索引的方法[6],通過使用GPU提高了MC算法的處理效率;王旭等初提出的最近鄰移動立方體方法[7]則通過建立等分點避免了大量的插值計算;熊邦書等提出的基于連通性的快速曲面重建算法[8]則通過減少遍歷體素的個數(shù)來提高算法的處理效率。

    1.3 MC算法的改進

    傳統(tǒng)的MC算法是通過對體數(shù)據(jù)所形成的體素進行逐個遍歷來生成小立方體內(nèi)的等值面,并最終合稱為體數(shù)據(jù)的表面。但是體數(shù)據(jù)中大多數(shù)體素并不與體數(shù)據(jù)的等值面相交,因此對體素進行逐個遍歷將浪費很多時間在空的體素處理上。腦圖譜的腦區(qū)都是一組體數(shù)據(jù)中的小部分數(shù)據(jù),傳統(tǒng)的MC算法在處理腦圖譜數(shù)據(jù)時重建一個較小的腦區(qū)也要遍歷龐大的體數(shù)據(jù)的每個體素。但是腦區(qū)表面只占整體數(shù)據(jù)的很小一部分,這使得傳統(tǒng)MC算法在重建腦圖譜數(shù)據(jù)時效率比一般的數(shù)據(jù)下降更為嚴重。為了解決這種問題我們提出了一種基于體素立方體頂點監(jiān)測的改進算法。本算法不采用遍歷的方式生成等值面,而是判斷出種子體素相鄰的體素中哪些體素中包含了等值面,再通過區(qū)域增長法將包含有等值面的體素全部提取出來進行等值面的繪制。

    改進算法的核心在于對種子體素相鄰的體素中包含等值面體素的監(jiān)測。我們通過對體素構成的立方體8個頂點的監(jiān)測來判斷與該體素接鄰的體素是否包含等值面。體素構成的立方體包含有8個頂點6個面,每個面包含4個頂點,在屬于同一個面的4個頂點中只要這四個頂點不同時大于或小于等值面閾值,則與種子體素共同包含該面的體素也包含有等值面。例如圖2所示左右兩個體素均包含同一個面,線段ab為該面與左邊的體素立方體中等值面的交線,因為該面同屬于右面的體素,所以在右面體素中的等值面也必然與此面相交于該交線。由此證明在右面的體素中也是含有等值面的。

    圖2 包含有共同面的兩個體素立方體

    而體素立方體各個面上有沒有與體素內(nèi)等值面的交線通過檢測該面的四個頂點就可以確定。四個頂點中我們假設它們都有兩種狀態(tài)即高于閾值和低于閾值狀態(tài),則一個面有1個點處于不同狀態(tài),2個點處于不同狀態(tài),和全部處于同一狀態(tài)三種情況。其中1點狀態(tài)不同有一種交線形式,2點有兩種,全部處于同一狀態(tài)則沒有交線出現(xiàn),如圖3所示。綜上所述通過對種子體素立方體一個方向的面的四個頂點進行監(jiān)測就可以判斷這個方向的接鄰體素是否含有等值面。

    圖3 體素內(nèi)等值面與體素表面相交的幾種情況

    通過上面的規(guī)則,只要給出相應的種子體素我們就能夠?qū)⒄麄€體數(shù)據(jù)中包含有等值面的體素提取出來進行處理,避免將體素逐個遍歷從而節(jié)省大量時間。

    1.4 改進算法的實現(xiàn)

    實現(xiàn)該算法首先要解決的問題是種子體素的選取問題,我們采用了二次取樣的方法進行種子點的提取,即遍歷時對每個維度增長速度進行加倍。這樣能夠提高遍歷的速度,一般二次取樣都會降低體數(shù)據(jù)的精度,但是本算法二次取樣后得到的體素在進行區(qū)域增長獲取邊界體素后并不會影響精度。在獲得了種子體素后就可以進行鄰接體素的檢測了,本算法使用了一個棧和一個哈希表來輔助算法計算。

    算法流程為:

    (1) 將種子體素入棧,監(jiān)測接鄰體素8個頂點,判斷哪些接鄰體素具有等值面,將具有等值面的接鄰體素入棧,并加入哈希表。

    (2) 監(jiān)測棧是否為空,為空則算法結(jié)束,不為空則將棧頂部的體素彈出。

    (3) 將彈出體素的八個頂點進行監(jiān)測,判斷哪些接鄰體素據(jù)有等值面,監(jiān)測具有等值面的接鄰體素是否在哈希表中,不在表中的接鄰體素入棧。

    (4) 為彈出的體素進行等值面繪制,轉(zhuǎn)到第(2)步。

    算法核心編程:

    //算法通過棧和哈希表提取含有等值面的體素立方體

    while (!readyCubes.empty())

    {

    idx=readyCubes.top();

    readyCubes.pop();

    //從圖譜數(shù)據(jù)文件取出組成一個立方體8個頂點的值

    s[0] = scalars[idx];

    s[1] = scalars[idx+1];

    s[2] = scalars[idx+1 + dims[0]];

    s[3] = scalars[idx + dims[0]];

    s[4] = scalars[idx + sliceSize];

    s[5] = scalars[idx+1 + sliceSize];

    s[6] = scalars[idx+1 + dims[0] + sliceSize];

    s[7] = scalars[idx + dims[0] + sliceSize];

    //向棧加入鄰接立方體,并為當前立方體添加哈希標記

    //監(jiān)測上方鄰接立方體在哈希表中是否有記錄以及是否有等值

    //面,如果有等值面且哈希表無記錄則將其入棧

    if (!(s[0] < min && s[1] < min && s[2] < min && s[3] < min||s[0] > max && s[1] > max && s[2] > max && s[3] > max))

    {

    if((idx-sliceSize)<0)

    continue;

    if (!CheckedCubes.find(idx-sliceSize))

    {

    readyCubes.push(idx-sliceSize);

    CheckedCubes.insert(idx-sliceSize);

    }

    }

    //監(jiān)測下方鄰接立方體在哈希表中是否有記錄以及是否有等值

    //面,如果有等值面且哈希表無記錄則將其入棧

    if (!(s[4] < min && s[5] < min && s[6] < min && s[7] < min||s[4] > max && s[5] > max && s[6] > max && s[7] > max))

    {

    if((idx+sliceSize)<0)

    continue;

    if (!CheckedCubes.find(idx+sliceSize))

    {

    readyCubes.push(idx+sliceSize);

    CheckedCubes.insert(idx+sliceSize);

    }

    }

    //監(jiān)測左方鄰接立方體在哈希表中是否有記錄以及是否有等值

    //面,如果有等值面且哈希表無記錄則將其入棧

    if (!(s[0] < min && s[3] < min && s[4] < min && s[7] < min||s[0] > max && s[3] > max && s[4] > max && s[7] > max))

    {

    ?

    }

    //監(jiān)測右方鄰接立方體在哈希表中是否有記錄以及是否有等值

    //面,如果有等值面且哈希表無記錄則將其入棧

    ?

    //監(jiān)測前方……

    ?

    //監(jiān)測后方……

    ?

    }

    2 VTK應用

    2.1 VTK簡介

    VTK(visualization toolkit)[9,10]是一個開源的免費軟件系統(tǒng),我們可以直接從官網(wǎng)上獲得最新版本。它廣泛應用于三維計算機圖形學、圖像處理和可視化。VTK具有強大的三維圖形功能,在支持體繪制的同時也保留了傳統(tǒng)的面繪制。其獨特的數(shù)據(jù)管道化處理使得VTK在處理大量數(shù)據(jù)時不用考慮內(nèi)存資源的限制。它將我們在可視化開發(fā)中常用的一些算法封裝起來,提供給我們相關的類,通過將開發(fā)中的細節(jié)屏蔽起來大大簡化了可視化開發(fā)中的復雜程度。此外VTK還能夠跨平臺并且支持數(shù)據(jù)并行處理,它能夠在任何Unix和Windows平臺上使用,而且具有良好的可移植性。

    2.2 VTK可視化管線

    VTK對數(shù)據(jù)進行可視化采用的是獨特的管線處理機制[11],通過管線機制VTK能夠讀取數(shù)據(jù)、對數(shù)據(jù)進行處理、將處理后的數(shù)據(jù)傳入渲染引擎進行顯示。這種機制減小了編程的復雜度,使數(shù)據(jù)的輸入與數(shù)據(jù)的處理過程相分離。將數(shù)據(jù)與處理算法引入后只需要保證輸入與處理過程的接口一致就可以了,而不用去考慮除了管線以外其他的因素。并且在對同一數(shù)據(jù)用不同的算法進行處理時,只需要更換相關處理過程的類就可以了。

    處理管線一般包括數(shù)據(jù)源、過濾器、映射器、表演者、渲染器和窗口。另外用戶接口和控制器雖然不是處理管線中必要的部分,但是作為應用程序用戶交互一部分也是十分重要的。處理管線根據(jù)數(shù)據(jù)類型的不同具有不同的數(shù)據(jù)處理過程,但是其本質(zhì)是相同的,流程如圖4所示。

    圖4 VTK可視化處理管線

    改進算法融入到VTK可視化管線的方法如下:

    (1) 讀入腦圖譜數(shù)據(jù)

    vtkImageReader *imreader = vtkImageReader::New();

    GetImageReader(imreader);

    (2) 建立數(shù)據(jù)處理過濾器,并將改進算法引入其中

    //其中MarchingCubesReduce為改進MC算法的封裝類

    MarchingCubesReduce *extractor = MarchingCubesReduce::New();

    extractor->SetInputConnection(imreader->GetOutputPort());

    extractor->SetValue(0,3);

    extractor->Update();

    (3) 建立映射器

    vtkPolyDataMapper *mapper = vtkPolyDataMapper::New();

    mapper->SetInputConnection(extractor->GetOutputPort());

    mapper->SetColorModeToMapScalars();

    mapper->SetScalarRange(1.0,20.0);

    (5) 添加表演者、渲染窗口以及控制器

    vtkActor *actor = vtkActor::New();

    actor->GetProperty()->SetColor(0.5,0.0,0.0);

    actor->SetMapper(mapper);

    vtkRenderer *ren = vtkRenderer::New();

    vtkRenderWindow *renWin = vtkRenderWindow::New();

    renWin->AddRenderer(ren);

    vtkRenderWindowInteractor *iren = vtkRenderWindowInteractor::New();

    iren->SetRenderWindow(renWin);

    在引入改進MC算法之后,該算法能夠提高過濾器處理腦圖譜數(shù)據(jù)的效率。

    2.3 生成圖像性能及質(zhì)量的優(yōu)化

    考慮到系統(tǒng)最后的運行效果,我們要對三維重建后得到的圖形進行優(yōu)化,減小它的復雜度、提高顯示效果,加快系統(tǒng)的運行速度。對生成圖像進行優(yōu)化主要從三方面入手:對原始數(shù)據(jù)進行精簡;選擇產(chǎn)生更簡單圖形的算法;對最終生成圖像采用優(yōu)化算法。

    VTK類庫中有多種方法能對原始數(shù)據(jù)進行簡化,縮小原始數(shù)據(jù)的規(guī)模。其中vtkImageResample類通過線性插值方式對輸入的數(shù)據(jù)進行重新采樣,重新采樣的結(jié)果可以更致密或者相對稀疏一些。通過該類我們能夠?qū)⒃紨?shù)據(jù)的質(zhì)量降低,從而降低數(shù)據(jù)的復雜程度,提高處理速度。另外VTK中的vtkImageShrink3D類能夠?qū)崿F(xiàn)二次抽樣,可以加速數(shù)據(jù)的輸入。通過該類的相關函數(shù)能夠控制抽樣的精度從而得到速度和質(zhì)量平衡的效果。該類搭配vtkImageResample類能夠大幅度地提高數(shù)據(jù)的處理速度。

    選擇性能更優(yōu)良的處理算法能夠使最后獲得的圖像有更好的效果,并提高處理速度。本文使用了MC算法計算腦圖譜數(shù)據(jù)的等值面,該算法得到的最后結(jié)果經(jīng)過優(yōu)化與未經(jīng)過優(yōu)化的結(jié)果相比其最后圖像所包含的三角面片大大減少。如圖5所示。優(yōu)化前腦區(qū)三維重建包含25 296個單元,優(yōu)化后包含15 412個單元。

    圖5 三維重建優(yōu)化前后比較

    如果對最后生成的圖像面片數(shù)量還是不滿意還可以通過使用VTK提供的vtkDecimatePro 類,它能夠有效地減少生成圖像的三角面片。不過削減比例過大有可能破壞圖形原本的拓撲結(jié)構。

    對圖形進行削減后雖然提高了性能但是有可能使圖像變得參差不齊,這時我們可以通過圖像平滑類vtkSmoothPoly-DataFilter對圖像進行處理,提高其顯示質(zhì)量。該類通過對生成圖像各個點的坐標進行修改,從而使構成圖形的網(wǎng)格更疏松,各個定點的分布更加均勻,最終使圖形的各個單元獲得更好的形狀。

    優(yōu)化圖像步驟如下:

    (1) 重新采樣縮小數(shù)據(jù)規(guī)模

    vtkImageResample *irreduce =vtkImageResample::New();

    Irreduce->SetInput(imreader.GetOutput());

    Irreduce->SetAxisMagnificationFactor(0, 0.5);

    Irreduce->SetAxisMagnificationFactor(1, 0.5);

    Irreduce->SetAxisMagnificationFactor(2, 0.5)

    (3) 使用vtkDecimatePro 類削減面片

    vtkDecimatePro *DecimatePro = vtkDecimatePro::New();

    DecimatePro->SetInputConnection(extractor->GetOutputPort() );

    DecimatePro->SetTargetReduction(0.3);

    DecimatePro->PreserveTopologyOn();

    (3) 對生成的圖像進行平滑處理

    vtkSmoothPolyDataFilter *smoother = vtkSmoothPolyDataFilter::New();

    smoother->SetlnputConnection(DecimatePro->GetOutput());

    smoother->SetNumberOflterations(200);

    (4) 計算面片的法向量

    vtkPolyDataNormals *Normals = vtkPolyDataNormals::New();

    Normals->SetInputConnection(smoother->GetOutputPort());

    Normals->FlipNormalsOn();

    3 系統(tǒng)結(jié)構及人機界面

    3.1 系統(tǒng)結(jié)構設計

    本文制作的腦圖譜三維重建系統(tǒng)基于VTK以及在醫(yī)學圖像處理平臺3D Slicer 開發(fā)中曾用到的KWWidgets[12]類庫。通過人機交互界面控制腦圖譜數(shù)據(jù)的輸入,以及腦圖譜數(shù)據(jù)的處理。根據(jù)不同的顯示要求對數(shù)據(jù)處理調(diào)用相應的算法和優(yōu)化處理過程,最終為腦圖譜的三維顯示、二維顯示、腦圖譜各個腦區(qū)的偵測和標記提供處理好的數(shù)據(jù)。系統(tǒng)結(jié)構如圖6所示。

    圖6 系統(tǒng)結(jié)構設計圖

    3.2 人機交互界面設計

    系統(tǒng)圖形界面使用KWWidgets界面庫制作,該類庫也是由VTK的開發(fā)公司開發(fā)的。因此與VTK庫能良好的配合使用,同時這套類庫提供了很多可視化編程要用到的控件,能極大地節(jié)省軟件開發(fā)的時間,此類庫還支持TCL/TK腳本語言,方便進行開發(fā)。

    整個系統(tǒng)的用戶界面分為三部分:(1) 主控制面板;(2) 主顯示窗口;(3) 二維圖像顯示窗口。如圖7所示。

    圖7 系統(tǒng)界面

    (1) 主控面板包含有相關文件的輸入輸出功能,為各個腦區(qū)建立標簽的功能,調(diào)整腦區(qū)三維圖形光照的功能,以及調(diào)整腦區(qū)三維圖形的顏色與透明度的功能。

    (2) 主顯示窗口用來顯示腦圖譜的三維重建圖像,用戶可以通過鼠標點擊拖拽和滾輪滾動控制三維圖像的角度變換,以及視角的拉近和推遠。

    (3) 二維圖像顯示窗口顯示二維的腦圖譜數(shù)據(jù),用戶可以在任意一個顯示窗口點擊,從而交互的改變其他窗口顯示的腦切面的位置。

    4 實驗及結(jié)果分析

    4.1 實驗數(shù)據(jù)介紹

    本系統(tǒng)要求的輸入文件為包含腦圖譜數(shù)據(jù)的NIfTI[13]格式文件。在本文中的實驗數(shù)據(jù)為來自于中國科學院自動化研究所腦網(wǎng)絡組研究中心的腦網(wǎng)絡組圖譜(Brainnetome Atlas)中的腦區(qū)分區(qū)結(jié)果[14](http://atlas.brainnetome.org)。文件維度為181×271×181。本圖譜文件包含有兩個大腦半球的39個腦區(qū),圖譜文件各個分區(qū)所包含的體素value值都為其所屬分區(qū)的編號值。如表3所示。

    表3 Brainnetome Atlas腦圖譜數(shù)據(jù)

    續(xù)表3

    此外由于實驗數(shù)據(jù)的腦圖譜文件不包含腦模板。因此本文在實驗中輸入文件還包含了與圖譜文件匹配的的cpbet腦模板文件,用來形成腦模板的三維圖像,使重構后的圖譜文件三維圖像能夠在模板上顯示。

    4.2 系統(tǒng)效率分析

    本系統(tǒng)在CPU為酷睿雙核2 GHz,3 GB內(nèi)存,顯卡為GT130m的平臺上對實驗數(shù)據(jù)進行處理。實驗得到的MC原始算法與改進后的效率對比如表4所示。

    表4 原始算法與改進算法效率對比

    從對比數(shù)據(jù)中可以看出原始MC算法重建整個腦圖譜表面比改進算法多進行了400多倍的遍歷次數(shù)。但是由于棧和哈希表所占用的時間耗費,實際處理時間雖然沒有遍歷次數(shù)下降的幅度大但是也能夠從12 814 ms加快到1852 ms,減少了85%的時間耗費。重建時間分為遍歷時間和繪制時間。繪制時間是對體素的等值面進行繪制的時間這個時間兩種算法基本持平。在遍歷時間上改進算法根據(jù)重建表面所占的體素數(shù)來決定耗費時間,原始算法則遍歷整個體數(shù)據(jù)空間,因此原始算法的遍歷時間基本上是不變的,改進算法提高效率就是源于遍歷時間相對于原始算法大大減少。圖8給出了隨著處理體素數(shù)量增長兩種算法的效率對比。

    圖8 改進算法與原始算法基于體素數(shù)的效率對比

    4.3 三維重建效果

    圖9給出了兩種算法重建的單個腦區(qū)圖形,改進算法與原始算法產(chǎn)生的圖形相差無幾。證明使用改進算法進行加速不會影響重建效果。圖10則給出了二維腦圖譜圖像與三維重建的對比,其中二維圖像為單層圖像看不見下層的腦區(qū),三維重建后所有的腦區(qū)都能夠方便地進行觀察。

    圖9 單個腦區(qū)的重建對比

    圖10 三維重建結(jié)果與二維圖譜對比

    5 結(jié) 語

    腦圖譜一直以來都是研究腦部結(jié)構和功能的重要手段。目前,腦圖譜相關的研究工作正在國內(nèi)外腦科學研究單位如火如荼的開展。醫(yī)學三維重建技術作為一門已經(jīng)成熟的技術應用在腦圖譜三維重構上極大地提升了腦圖譜的可視化效果。本文對腦圖譜三維重建顯示過程中處理算法的選擇、圖形的優(yōu)化以及用戶交互進行設計,并完成了一個處理速度快、顯示質(zhì)量達到應用要求的腦圖譜三維重建和可視化系統(tǒng)。并且使用Brainnetome Atlas的臨床數(shù)據(jù)對系統(tǒng)進行測試,基于改進MC算法的三維重建系統(tǒng)確實在腦圖譜數(shù)據(jù)的處理上展現(xiàn)了高于原始算法的效率。

    在本文課題進行過程中,VTK類庫得到了大量的應用。它作為一款開源的可視化軟件,其功能十分強大,而且可以根據(jù)個人需要對其進行改進,這方面是其他軟件做不到的。

    以上工作只是腦圖譜三維重建及可視化的一個初步探索,很多功能還尚待進一步實現(xiàn)。隨著對體數(shù)據(jù)三維重建算法在并行計算和GPU處理等方面改進的實現(xiàn),我們將進一步豐富該系統(tǒng)的功能,提高其重建效率,改善其可視化效果。

    [1] 常旖旎,魯雯,聶生東.醫(yī)學圖像三維可視化技術及其應用[J].中國醫(yī)學物理學雜志,2012,29(2):3254-3258.

    [2] 解立志,周明全,田沄,等.基于光線投射算法的腦血管體繪制技術[J].系統(tǒng)仿真學報,2012,24(9):1-5.

    [3] 陳卿.醫(yī)學影像三維可視化Marching Cubes算法研究[D].成都:電子科技大學,2011.

    [4] Nielson G M.On marching cubes[J].Visualization and Computer Graphics,IEEE Transactions on,2003,9(3):283-297.

    [5] Lopes A,Brodlie K.Improving the robustness and accuracy of the marching cubes algorithm for isosurfacing[J].Visualization and Computer Graphics,IEEE Transactions on,2003,9(1):16-29.

    [6] 張迎平,高國賢,陸一峰,等.基于區(qū)間樹硬件加速索引的Marching Cubes算法[J].計算機輔助設計與圖形學學報,2012,24(7):871-878.

    [7] 王旭初,王贊.基于最近鄰Marching Cubes的醫(yī)學圖像三維重建[J].計算機工程與應用,2012,48(18):154-158.

    [8] 熊邦書,何明一,俞華璟.基于空間連通性的快速曲面重建算法[J].系統(tǒng)仿真學報,2005,17(1):75-78.

    [9] 洪濤,潘志方,林立本,等.VTK醫(yī)學圖像三維重建應用及實現(xiàn)[J].計算機系統(tǒng)應用,2011,20(4):127-130.

    [10] Sun S,He J,Ma L.3D visualization system for medical image based on VTK and OK Series Image Board[C]//Computer Science and Service System (CSSS),2011 International Conference on.IEEE,2011:951-954.

    [11] Caban J J,Joshi A,Nagy P.Rapid development of medical imaging tools with open-source libraries[J].Journal of digital imaging,2007,20(1):83-93.

    [12] Fang L,Wan Z,Zeng W,et al.Research on Internationalization of 3D Slicer[C]//Computer Sciences and Applications (CSA),2013 International Conference on.IEEE,2013:469-473.

    [13] Larobina M,Murino L.Medical Image File Formats[J].Journal of digital imaging,2014,27(2):200-206.

    [14] Jiang T,Zhou Y,Liu B,et al.Brainnetome-wide association studies in schizophrenia:The advances and future[J].Neuroscience & Biobehavioral Reviews,2013,37(10):2818-2835.

    3D RECONSTRUCTION OF BRAIN ATLAS BASED ON MODIFIED MARCHING CUBES ALGORITHM

    Zhu Kai1Fan Lingzhong2,3Li Haifang1*

    1(SchoolofComputerScienceandTechnology,TaiyuanUniversityofTechnology,Taiyuan030024,Shanxi,China)2(BrainnetomeCenter,InstituteofAutomation,ChineseAcademyofSciences,Beijing100190,China)3(NationalLaboratoryofPatternRecognition,InstituteofAutomation,ChineseAcademyofSciences,Beijing100190,China)

    Along with that brain atlas plays more and more important role in teaching and application of neurosurgery operation, 3D reconstruction of brain atlas has become a major research direction of brain atlas. In order to improve the efficiency of MarchingCubes algorithm in extracting isosurface of brain atlas data, we propose an improved algorithm suitable for 3D reconstruction of brain atlas through the study on the traversal method and topological structure of MarchingCubes algorithm, and realise the brain atlas 3D reconstruction system based on the improved algorithm. We use brain atlas data of Brainnetome Atlas as the experimental data to test the efficiency of the system, experimental results show that relative to original algorithm the improved algorithm really improves the efficiency of 3D reconstruction system in processing brain atlas data.

    Isosurface MarchingCubes 3D reconstruction Brain atlas

    2014-09-06。國家自然科學基金項目(61373101);山西省自然科學(青年科技研究)基金項目(2014021022-5);山西省科技攻關項目(20130311037-4)。朱愷,碩士生,主研領域:智能信息處理。樊令仲,助理研究員。李海芳,教授。

    TP391.7

    A

    10.3969/j.issn.1000-386x.2016.04.042

    猜你喜歡
    體素等值立方體
    基于超體素聚合的流式細胞術自動門控方法
    疊出一個立方體
    基于多級細分的彩色模型表面體素化算法
    異步電動機等值負載研究
    防爆電機(2020年5期)2020-12-14 07:03:50
    運用邊界狀態(tài)約束的表面體素加密細分算法
    基于體素格尺度不變特征變換的快速點云配準方法
    圖形前線
    立方體星交會對接和空間飛行演示
    太空探索(2016年9期)2016-07-12 09:59:53
    折紙
    電網(wǎng)單點等值下等效諧波參數(shù)計算
    午夜免费男女啪啪视频观看 | 久久久久国内视频| 99热这里只有是精品在线观看| 老女人水多毛片| 不卡视频在线观看欧美| 亚洲自拍偷在线| 久久久久久国产a免费观看| 国产成人a∨麻豆精品| 免费看av在线观看网站| 日韩一区二区视频免费看| 人人妻人人澡人人爽人人夜夜 | 国产精品无大码| 好男人在线观看高清免费视频| 一个人观看的视频www高清免费观看| 亚洲最大成人av| 欧美一区二区国产精品久久精品| 国产一区二区三区在线臀色熟女| 国产一区二区三区在线臀色熟女| 偷拍熟女少妇极品色| 欧美一区二区国产精品久久精品| 久久久精品94久久精品| 女生性感内裤真人,穿戴方法视频| 日本精品一区二区三区蜜桃| 免费搜索国产男女视频| 一级av片app| 看黄色毛片网站| 亚洲一区二区三区色噜噜| 男女做爰动态图高潮gif福利片| 国产精品女同一区二区软件| 老司机午夜福利在线观看视频| 亚洲国产日韩欧美精品在线观看| 午夜影院日韩av| 淫秽高清视频在线观看| 国产单亲对白刺激| 国产精品国产高清国产av| 少妇人妻一区二区三区视频| 成人午夜高清在线视频| 中文字幕免费在线视频6| 成人永久免费在线观看视频| 国产一区二区三区av在线 | 综合色丁香网| 麻豆精品久久久久久蜜桃| 男女那种视频在线观看| 亚洲人成网站在线播放欧美日韩| 一进一出好大好爽视频| videossex国产| 免费看日本二区| 波多野结衣巨乳人妻| 不卡一级毛片| 黑人高潮一二区| 亚洲婷婷狠狠爱综合网| 男人狂女人下面高潮的视频| 熟妇人妻久久中文字幕3abv| 日韩欧美精品v在线| 久久精品综合一区二区三区| 乱系列少妇在线播放| 最近最新中文字幕大全电影3| 国产毛片a区久久久久| 蜜臀久久99精品久久宅男| 免费高清视频大片| 别揉我奶头~嗯~啊~动态视频| 久久久久久久亚洲中文字幕| 国产成人精品久久久久久| 美女免费视频网站| 亚洲七黄色美女视频| 最近2019中文字幕mv第一页| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久亚洲av鲁大| 精品久久久久久久久av| 久久精品人妻少妇| 成人无遮挡网站| 我要搜黄色片| 卡戴珊不雅视频在线播放| 我的女老师完整版在线观看| 自拍偷自拍亚洲精品老妇| 久久欧美精品欧美久久欧美| 村上凉子中文字幕在线| 久久亚洲精品不卡| 久久精品国产自在天天线| 日韩中字成人| 精华霜和精华液先用哪个| 成人国产麻豆网| 少妇被粗大猛烈的视频| 日韩精品有码人妻一区| 亚洲国产精品sss在线观看| 国产精品三级大全| 国产人妻一区二区三区在| 国产又黄又爽又无遮挡在线| 亚洲色图av天堂| 91久久精品电影网| 性色avwww在线观看| 变态另类丝袜制服| 日本免费a在线| 天美传媒精品一区二区| 欧美日韩乱码在线| 国产高清有码在线观看视频| 日本免费a在线| 桃色一区二区三区在线观看| 天天一区二区日本电影三级| 丰满的人妻完整版| 亚洲欧美清纯卡通| 午夜精品一区二区三区免费看| 日本-黄色视频高清免费观看| 久久久久久久久大av| 六月丁香七月| 精品久久久久久久久久久久久| 搞女人的毛片| 亚洲激情五月婷婷啪啪| 少妇丰满av| 久久久久久久午夜电影| 成人无遮挡网站| 日本一本二区三区精品| 久久人人精品亚洲av| 搡老妇女老女人老熟妇| 亚洲天堂国产精品一区在线| 男人舔奶头视频| 激情 狠狠 欧美| 麻豆精品久久久久久蜜桃| 欧美绝顶高潮抽搐喷水| 国产精品,欧美在线| 久久久色成人| 欧美人与善性xxx| 3wmmmm亚洲av在线观看| 99久国产av精品国产电影| 亚洲人成网站高清观看| 插逼视频在线观看| 精品少妇黑人巨大在线播放 | 尾随美女入室| 色综合色国产| 不卡一级毛片| 91av网一区二区| 床上黄色一级片| 搡老妇女老女人老熟妇| 亚洲四区av| 特大巨黑吊av在线直播| 久久久精品大字幕| 精品少妇黑人巨大在线播放 | 晚上一个人看的免费电影| 人妻制服诱惑在线中文字幕| 日本精品一区二区三区蜜桃| 国产v大片淫在线免费观看| 国产黄片美女视频| 亚洲一级一片aⅴ在线观看| 好男人在线观看高清免费视频| 国产精品人妻久久久久久| 搡老岳熟女国产| 18禁在线无遮挡免费观看视频 | 免费观看人在逋| 午夜精品国产一区二区电影 | 自拍偷自拍亚洲精品老妇| 卡戴珊不雅视频在线播放| 亚洲人成网站高清观看| 亚洲自拍偷在线| 亚洲精品影视一区二区三区av| 真实男女啪啪啪动态图| 中文字幕av成人在线电影| 久久综合国产亚洲精品| 精品久久久噜噜| 最好的美女福利视频网| 一本久久中文字幕| 国产精品人妻久久久久久| 一区福利在线观看| 69av精品久久久久久| 中文字幕免费在线视频6| 亚洲中文字幕一区二区三区有码在线看| 色哟哟·www| 如何舔出高潮| 一区二区三区高清视频在线| 卡戴珊不雅视频在线播放| 国内精品久久久久精免费| 亚洲精品在线观看二区| 禁无遮挡网站| 国产精品99久久久久久久久| 亚洲熟妇中文字幕五十中出| 成人亚洲精品av一区二区| 九九在线视频观看精品| 丰满乱子伦码专区| 老司机福利观看| 国产一级毛片七仙女欲春2| 狂野欧美激情性xxxx在线观看| 男人狂女人下面高潮的视频| 久久热精品热| 嫩草影院入口| 亚洲国产日韩欧美精品在线观看| 国产精品,欧美在线| 韩国av在线不卡| 少妇高潮的动态图| 噜噜噜噜噜久久久久久91| 成人亚洲精品av一区二区| 亚洲色图av天堂| 国产黄a三级三级三级人| 如何舔出高潮| 国产免费一级a男人的天堂| 久久精品综合一区二区三区| 亚洲不卡免费看| 精品一区二区三区人妻视频| 日日摸夜夜添夜夜爱| 最近最新中文字幕大全电影3| 久久精品人妻少妇| 插阴视频在线观看视频| 亚洲精华国产精华液的使用体验 | 亚洲欧美日韩卡通动漫| 麻豆一二三区av精品| 日韩中字成人| 乱系列少妇在线播放| 欧美不卡视频在线免费观看| 99九九线精品视频在线观看视频| 亚洲国产精品久久男人天堂| 中文字幕av成人在线电影| 日韩av在线大香蕉| 日日摸夜夜添夜夜添av毛片| 日本 av在线| 99在线视频只有这里精品首页| a级一级毛片免费在线观看| 色哟哟哟哟哟哟| 国产 一区 欧美 日韩| 国产成人freesex在线 | 午夜老司机福利剧场| 麻豆乱淫一区二区| 日韩在线高清观看一区二区三区| 久久欧美精品欧美久久欧美| 九色成人免费人妻av| 校园人妻丝袜中文字幕| 国产黄色视频一区二区在线观看 | 国产精品电影一区二区三区| 国产成人91sexporn| 午夜激情福利司机影院| 国产免费一级a男人的天堂| 91在线精品国自产拍蜜月| 久久人妻av系列| 最好的美女福利视频网| 精品欧美国产一区二区三| 久久精品国产亚洲av涩爱 | 99在线视频只有这里精品首页| 人妻制服诱惑在线中文字幕| 欧美性感艳星| 免费人成在线观看视频色| 一区福利在线观看| 欧美中文日本在线观看视频| 日本免费一区二区三区高清不卡| 99riav亚洲国产免费| 尾随美女入室| 少妇的逼好多水| 97超级碰碰碰精品色视频在线观看| 又黄又爽又刺激的免费视频.| 国产高清有码在线观看视频| 在线免费十八禁| 国产在线精品亚洲第一网站| 搡女人真爽免费视频火全软件 | 久久久久久久久久成人| 国产国拍精品亚洲av在线观看| 欧美成人免费av一区二区三区| 天天躁夜夜躁狠狠久久av| 一边摸一边抽搐一进一小说| 成年版毛片免费区| 精品免费久久久久久久清纯| 18+在线观看网站| 国内揄拍国产精品人妻在线| 97超视频在线观看视频| 国产亚洲精品久久久久久毛片| 欧美bdsm另类| 99热6这里只有精品| 女同久久另类99精品国产91| 亚洲欧美日韩卡通动漫| 特级一级黄色大片| 亚洲精品国产成人久久av| 国产乱人偷精品视频| 国产av在哪里看| 美女cb高潮喷水在线观看| 一个人看视频在线观看www免费| 国产亚洲91精品色在线| 乱码一卡2卡4卡精品| 亚洲欧美日韩无卡精品| 久久人人爽人人片av| 夜夜看夜夜爽夜夜摸| 嫩草影院新地址| 亚洲丝袜综合中文字幕| 长腿黑丝高跟| 亚洲天堂国产精品一区在线| www.色视频.com| 麻豆av噜噜一区二区三区| 国产精品一区二区性色av| 免费搜索国产男女视频| 国产精品亚洲美女久久久| 国产欧美日韩精品亚洲av| 亚洲第一电影网av| 国产伦精品一区二区三区视频9| 国产视频内射| 亚洲真实伦在线观看| 老司机影院成人| 三级经典国产精品| 99久久无色码亚洲精品果冻| 国产精品乱码一区二三区的特点| a级毛片a级免费在线| 亚洲在线自拍视频| 国产精品一及| 最近最新中文字幕大全电影3| 少妇丰满av| 国产69精品久久久久777片| 国产高潮美女av| 春色校园在线视频观看| 黄色日韩在线| 一区二区三区免费毛片| 国产高清不卡午夜福利| 免费搜索国产男女视频| 三级经典国产精品| 嫩草影院入口| 免费av不卡在线播放| 日本免费a在线| 国产精品人妻久久久久久| 天堂√8在线中文| av黄色大香蕉| 日韩欧美一区二区三区在线观看| 亚洲国产精品国产精品| 久久这里只有精品中国| 一卡2卡三卡四卡精品乱码亚洲| av中文乱码字幕在线| 国产黄色视频一区二区在线观看 | 国产亚洲精品综合一区在线观看| 国产国拍精品亚洲av在线观看| 国产综合懂色| 免费不卡的大黄色大毛片视频在线观看 | 三级男女做爰猛烈吃奶摸视频| 日韩成人av中文字幕在线观看 | 丝袜美腿在线中文| 人妻制服诱惑在线中文字幕| 一夜夜www| 国产欧美日韩精品一区二区| 九色成人免费人妻av| 久久久久久久亚洲中文字幕| 久久精品国产清高在天天线| 国内久久婷婷六月综合欲色啪| 午夜影院日韩av| 精品无人区乱码1区二区| 国产精品亚洲一级av第二区| 中文字幕熟女人妻在线| 黄色欧美视频在线观看| 国产精品国产三级国产av玫瑰| 一进一出抽搐动态| 女同久久另类99精品国产91| 插阴视频在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 日韩中字成人| 精品无人区乱码1区二区| 日本免费一区二区三区高清不卡| 精品人妻熟女av久视频| 国产精品久久久久久亚洲av鲁大| 桃色一区二区三区在线观看| 人人妻人人澡欧美一区二区| 久久精品国产99精品国产亚洲性色| 久久欧美精品欧美久久欧美| 久久精品国产亚洲网站| 日产精品乱码卡一卡2卡三| 欧美xxxx性猛交bbbb| 国产精品一区二区性色av| 国内精品宾馆在线| 久久精品91蜜桃| 久久久久国内视频| 亚洲中文字幕日韩| 国产欧美日韩一区二区精品| 亚洲第一区二区三区不卡| 性色avwww在线观看| 小说图片视频综合网站| 高清毛片免费观看视频网站| 久久精品91蜜桃| 久久国内精品自在自线图片| 国产成人精品久久久久久| 日本欧美国产在线视频| 两个人的视频大全免费| 能在线免费观看的黄片| 亚洲av中文av极速乱| 一级黄片播放器| 日韩欧美 国产精品| 国产精品嫩草影院av在线观看| 寂寞人妻少妇视频99o| 久久九九热精品免费| 寂寞人妻少妇视频99o| 久久国内精品自在自线图片| 91久久精品电影网| 国产麻豆成人av免费视频| 成人毛片a级毛片在线播放| 精品一区二区三区视频在线| 国产乱人视频| 麻豆av噜噜一区二区三区| 蜜臀久久99精品久久宅男| 国产亚洲欧美98| 日日啪夜夜撸| 干丝袜人妻中文字幕| 国产一区二区三区av在线 | 国产大屁股一区二区在线视频| 精品少妇黑人巨大在线播放 | 狂野欧美激情性xxxx在线观看| 亚洲第一电影网av| 国产男靠女视频免费网站| 91久久精品国产一区二区成人| 老司机午夜福利在线观看视频| 亚洲中文字幕日韩| 欧美性感艳星| 少妇的逼水好多| 丰满乱子伦码专区| 搡老熟女国产l中国老女人| 成人漫画全彩无遮挡| 国产精品99久久久久久久久| 女人十人毛片免费观看3o分钟| 成人高潮视频无遮挡免费网站| 日日干狠狠操夜夜爽| 黄色配什么色好看| 日韩高清综合在线| 久久精品91蜜桃| 国产片特级美女逼逼视频| 国产一级毛片七仙女欲春2| 别揉我奶头~嗯~啊~动态视频| 久久久久久九九精品二区国产| 午夜激情欧美在线| 亚洲自拍偷在线| 亚洲欧美精品综合久久99| 搡老岳熟女国产| 欧美高清性xxxxhd video| 日韩制服骚丝袜av| 久久久久久久亚洲中文字幕| 一级毛片电影观看 | 麻豆乱淫一区二区| 日本色播在线视频| 国产久久久一区二区三区| 亚洲美女视频黄频| 天堂网av新在线| 少妇的逼好多水| 欧美高清成人免费视频www| 成人漫画全彩无遮挡| 特级一级黄色大片| 欧美成人一区二区免费高清观看| 色哟哟哟哟哟哟| 在线免费观看不下载黄p国产| 国产三级中文精品| 简卡轻食公司| 在线天堂最新版资源| 精品一区二区三区人妻视频| 亚洲av免费在线观看| 99热这里只有是精品在线观看| 国产精品福利在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 久久久久久伊人网av| 久久久久久久亚洲中文字幕| 日韩强制内射视频| 欧美绝顶高潮抽搐喷水| 亚洲婷婷狠狠爱综合网| АⅤ资源中文在线天堂| 99精品在免费线老司机午夜| 偷拍熟女少妇极品色| 又黄又爽又免费观看的视频| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品久久久久久毛片| 亚洲av免费高清在线观看| 12—13女人毛片做爰片一| 狂野欧美白嫩少妇大欣赏| 亚洲自拍偷在线| 内地一区二区视频在线| 日韩欧美三级三区| 国产激情偷乱视频一区二区| 免费看美女性在线毛片视频| 亚洲精品一卡2卡三卡4卡5卡| 高清毛片免费观看视频网站| 国产乱人视频| 晚上一个人看的免费电影| 亚洲美女视频黄频| 免费观看精品视频网站| 美女xxoo啪啪120秒动态图| 熟妇人妻久久中文字幕3abv| videossex国产| 精品一区二区三区人妻视频| 久久人人爽人人片av| 精华霜和精华液先用哪个| 国产av一区在线观看免费| 亚洲自偷自拍三级| 人人妻人人看人人澡| 亚洲av成人av| 又黄又爽又刺激的免费视频.| 免费无遮挡裸体视频| 国产三级在线视频| 亚洲欧美清纯卡通| 少妇高潮的动态图| 男女之事视频高清在线观看| 亚洲欧美日韩高清在线视频| 国产蜜桃级精品一区二区三区| 国产中年淑女户外野战色| 寂寞人妻少妇视频99o| 少妇的逼好多水| 黄色欧美视频在线观看| 99在线人妻在线中文字幕| 99热这里只有是精品50| 亚洲成人中文字幕在线播放| 一卡2卡三卡四卡精品乱码亚洲| 成人国产麻豆网| 亚洲天堂国产精品一区在线| www.色视频.com| 校园春色视频在线观看| 男女下面进入的视频免费午夜| 99热全是精品| 最好的美女福利视频网| 麻豆av噜噜一区二区三区| 97人妻精品一区二区三区麻豆| 国产精品久久久久久久久免| 色在线成人网| 日本撒尿小便嘘嘘汇集6| 日产精品乱码卡一卡2卡三| 国产人妻一区二区三区在| 中文亚洲av片在线观看爽| 偷拍熟女少妇极品色| 搡老熟女国产l中国老女人| 我的女老师完整版在线观看| 蜜桃亚洲精品一区二区三区| 亚洲真实伦在线观看| 久99久视频精品免费| 成人特级黄色片久久久久久久| 嫩草影院精品99| 人人妻人人看人人澡| 精品久久久久久久久av| 十八禁国产超污无遮挡网站| 久久亚洲精品不卡| 99热这里只有精品一区| 日韩精品有码人妻一区| 色av中文字幕| 看非洲黑人一级黄片| 欧美精品国产亚洲| 亚洲国产精品国产精品| 亚洲av.av天堂| 看片在线看免费视频| 亚洲七黄色美女视频| 精品熟女少妇av免费看| 干丝袜人妻中文字幕| 久久精品国产清高在天天线| 中国美白少妇内射xxxbb| 日韩高清综合在线| 久久久午夜欧美精品| 99久久成人亚洲精品观看| 欧美一区二区国产精品久久精品| 久久久久国产精品人妻aⅴ院| 好男人在线观看高清免费视频| 亚洲成人精品中文字幕电影| 精品一区二区三区av网在线观看| 久久综合国产亚洲精品| 国产熟女欧美一区二区| 干丝袜人妻中文字幕| 亚洲精品色激情综合| 色5月婷婷丁香| 99久久精品国产国产毛片| 亚洲自拍偷在线| 观看免费一级毛片| 免费观看人在逋| 寂寞人妻少妇视频99o| 久久久国产成人精品二区| 国产精品人妻久久久影院| 国产精品不卡视频一区二区| 成人无遮挡网站| 亚洲精品一区av在线观看| av黄色大香蕉| 一级毛片久久久久久久久女| 寂寞人妻少妇视频99o| 99热全是精品| 亚洲一级一片aⅴ在线观看| 特级一级黄色大片| 亚洲av免费高清在线观看| 精品久久久久久久久亚洲| 联通29元200g的流量卡| 欧美激情在线99| 99久久久亚洲精品蜜臀av| 成人亚洲精品av一区二区| 男女边吃奶边做爰视频| 免费无遮挡裸体视频| 最近2019中文字幕mv第一页| 精品一区二区免费观看| 搡老熟女国产l中国老女人| 俺也久久电影网| 成人三级黄色视频| 一级毛片aaaaaa免费看小| 少妇猛男粗大的猛烈进出视频 | 日本 av在线| 国产精品野战在线观看| 欧美激情在线99| 亚洲精品亚洲一区二区| 老女人水多毛片| 欧美日韩乱码在线| 成人毛片a级毛片在线播放| 国产精品久久电影中文字幕| 亚洲最大成人av| 日本-黄色视频高清免费观看| 男人狂女人下面高潮的视频| 国产国拍精品亚洲av在线观看| 国语自产精品视频在线第100页| 麻豆乱淫一区二区| 日韩欧美三级三区| 久久久久国产精品人妻aⅴ院| 美女 人体艺术 gogo| 日韩制服骚丝袜av| 欧美+日韩+精品| 欧美色视频一区免费| 日韩精品有码人妻一区| 日韩成人伦理影院| 给我免费播放毛片高清在线观看| 亚洲自拍偷在线| 成人欧美大片| 国产成人a∨麻豆精品| 神马国产精品三级电影在线观看| 又黄又爽又刺激的免费视频.| 最近中文字幕高清免费大全6| 又粗又爽又猛毛片免费看| 久久人妻av系列| 日日干狠狠操夜夜爽| 日韩一本色道免费dvd| 国产高清三级在线| 国产高清有码在线观看视频| 此物有八面人人有两片| 国产在视频线在精品| 欧美日韩乱码在线| 中文字幕av成人在线电影|