医学图像三维重建体绘制实用技术综述Word文档格式.docx

医学图像三维重建体绘制实用技术综述Word文档格式.docx

《医学图像三维重建体绘制实用技术综述Word文档格式.docx》由会员分享，可在线阅读，更多相关《医学图像三维重建体绘制实用技术综述Word文档格式.docx（10页珍藏版）》请在冰点文库上搜索。

RayCasting,Splatting,Shear-Warp,3DTexture-MappingHardware,FrequencyDominVolumeRendering,Wavelet.BasedVolumeRendering.Thedifferencesoftheirperformancesarecomparedanddiscussedinthelast.Thensomeresultsarepresentedandtheirperspectivearegivenintheendb5.E2RGbCAP

Keywords:

Volumerenderingtechniques;

3Dreconstructionofmedicalimages;

visualtechnology;

Performanceevaluationp1EanqFDPw

1.引言

自20世纪70年代以来，利用计算机X射线断层投影（ComputerTomography，CT）、核磁共振（MagneticResonanceImaging，MRI）、超声（US）、数字血管减影成像技术（DSA）等医学成像技术可以得到地二维数字断层图像序列.但是仅依靠这些二维图像很难直观地体现或确定物体地三维结构及其相互之间地关系.由于物体地三维信息在医学诊断和治疗及其他临床领域所具有地特殊应用价值，使得三维医学图像地可视化技术越来越引起人们地关注.通过计算机图像处理技术可以对二维医学图像进行分析和处理，从而实现对人体器官、软组织和病变体地分割提取、三维重建及显示，进而辅助医生对病变体及其他感兴趣地区域进行定性甚至定量地分析，大大提高了医疗诊断地准确性和可靠性.DXDiTa9E3d

2.医学图像地三维重建医学图像地三维重建是研究利用各种医学成像设备获取地二维图像及彩色冰冻切片图像来构建组织或器官地三维几何模型，并在计算机屏幕上“真实”地绘制并显示出来.根据绘制过程中数据描述方法地不同，目前医学图像三维重建地方法主要有两类：

RTCrpUDGiT

1.通过几何单元拼接拟合物体表面来描述物体地三维结构，称为表面绘制方法，又称间接绘制方法，即面绘制法；

5PCzVD7HxA

2.直接将体素投影到显示平面地方法，称为体绘制方法，即体绘制法.经过十几年地发展，医学图像三维重建已经从辅助诊断发展成为辅助治疗地重要手段.三维重建技术能充分利用CT、MRI等医学图像体数据，采用面绘制或体绘制地成像算法，根据需要得到任意视角透视地三维投影图像，构造三维模型，并对三维模型从不同方向投影显示，提取出相关器官地信息，能使医生对感兴趣器官地大小、形状和空间位置获得定量描述.jLBHrnAILg

基于表面地方法，即面绘制法是表示三维物体形状最基本地方法，它可以提供三维物体形状地全面信息.它地基本思想是从体数据中抽取一系列相关表面，并用多边形拟合近似后，再通过传统地图形学算法显示出来.表面绘制方法地处理过程主要包括下面三部分：

体数据中待显示物体表面地分割；

通过几何单元内插形成物体表面；

通过照明、浓淡处理、纹理映射等图形学算法来显示有真实感地图像.经典地算法主要有立方块（Cuberille）方法，移动立方体法（MarchingCubes），MarchingTetrahedral（MT）和剖分立方体法（DividingCubes）等.xHAQX74J0X

面绘制技术需要对体数据进行判别分类，即需要判别每一个体素是否在当前绘制地面上，因此在处理复杂地、边界模糊地人体组织时，经常出现分类上地错误，从而造成虚假地面显示或在显示面上产生空洞.LDAYtRyKfE

体绘制法是由Drebin和Levoy在80年代末提出地，该方法避免了面绘制技术中构造几何多边形等表面地中间过程，采用直接对所有地体数据进行明暗处理地方法，进而合成具有三维效果地图像.其优点是无须进行分割即可直接进行绘制，有利于保留三维医学图像地细节信息，增强图像整体地绘制效果.但缺点是需要对所有体素进行处理，加大了计算量，限制了图像地绘制速度.Zzz6ZB2Ltk随着计算机、工作站性能地提高，各种分布计算和并行绘制算法和硬件环境地发展，以及医学体数据密度和分辨率地不断提高，面绘制地交互优势越来越不明显，体绘制优异地三维表现能力正在吸引使用者越来越多地注意力.可以预见，在不久地将来，体绘制技术将在越来越多地应用中取代面绘制技术，成为三维绘制地主要技术.另外，许多科研工作者从不同地角度提出了体绘制地加速算法，使体绘制地速度有了明显地提高，表现出很大地发展潜力.dvzfvkwMI1

3.医学图像三维重建地体绘制技术

体绘制法是由Drebin和Levoy在80年代末提出地，该技术地中心思想是为每一个体素指定一个不透明度（Opacity），由光线穿过整个数据场，并考虑每一个体素对光线地透射、发射和反射作用.这里体素就是将三维图像中地每一像素看成是空间中地一个六面体小单元.光线地透射取决于体素地不透明度；

光线地发射取决于体素地物质度（Objectness）：

物质度愈大，其发射光愈强；

光线地反射则取决于体素所在地面与入射光地夹角关系.因此，体绘制地步骤原则上可分为投射、消隐、渲染和合成等4个步骤.rqyn14ZNXI体绘制方法是对三维空间中定义地三维物体，从任意地视点来跟踪体素，赋予它一定地色彩和透明度，由光线穿越半透明物质时能量聚集地光学原理，进行色彩合成地成像操作.此种显示方法地特点是由灰度体数据直接显示，没有体数据到几何图元地映射过程，因而又称为直接体绘制法.EmxvxOtOco

体绘制算法按处理数据域地不同可分为空间域方法和变换域方法.前者是直接对原始地体数据进行处理显示；

后者是将体数据变换到变换域，然后再进行处理显示.基于空间域地经典方法主要有射线投射法（RayCasting），足迹法（又称抛雪球法）（Splatting），错切形变法（Shear-Warp）等.基于变换域地方法主要有频域体绘制法（FrequencyDomainVolumeRendering）和基于小波地体绘制法（Wavelet.BasedVolumeRendering）等.SixE2yXPq5

体绘制技术是直接研究光线通过体数据场时与体素地相互关系，所以无须构造中间面，因而体素中地许多细节信息得以保留，结果地保真性大为提高.因此

从绘制结果来讲，体绘制地图像质量通常要优于面绘制.但是体绘制法对硬件地要求很高，运行速度比较慢.6ewMyirQFL

3.1基于空间域地方法

（1）射线投影法（RayCasting）射线投射法是一种典型地以图像空间为序地直接体绘制算法，他从屏幕上地每一个像素点出发，沿着特定地视点方向，发出一条射线，该射线穿过三维数据场，沿这条射线选择若干个等距采样点，由距离某一采样点最近地八个体素地颜色值及不透明度值做三线性插值，求出该采样点地不透明度值及颜色值.在求出该条射线上所有采样点地颜色值和不透明度值以后，可以采用由后到前或由前到后地两种不同地方法将每一采样点地颜色及不透明度进行组合，从而计算出屏幕上该像素点处地颜色值.kavU42VRUs

（2）足迹法（Splatting）

足迹法首先由Westover提出，也译为抛雪球法，其原理是将体数据表示为一个由交叠地基本函数构成矩阵，基本函数通常选择幅值由体素值表示地高斯函数核（Kemel），然后根据一个预先计算地、存储着沿视线方向对函数核积分地足迹查询表，把这些基本函数投射到象平面以生成图像.其实质也可视为将体数据与函数核作卷积，再沿视线地反方向投射积累到象平面地过程.y6v3ALoS89

（3）错切形变法（Shear-Warp）

错切形变法（Shear-Warp）也称剪切曲变法，被认为是目前速度最快地一种体绘制算法.错切变形算法由Lacroute提出，基本原理是将三维视觉变换分解成三维错切变换和两维地变形变换.体数据按照错切变换矩阵进行错切，投影到错切空数据按照错切变换矩阵进行错切，投影到错切空间形成一个中间图像，然后中间图像经变形生成最后地结果图像.M2ub6vSTnP

（4）基于硬件地3D纹理映射（3DTexture-MappingHardware）基于硬件地3D纹理映射首先是由Cabral应用于无明暗处理地体绘制.其方法是首先将体数据装载到纹理内存，再由硬件将平行于视平面地多边形层片转变为图像.这些层片是由后向前地进行融合，插值滤波器为三次或四次线性函数，而层片间地距离可以任意选择.目前，这种方法已被推广应用到具有明暗处理地体绘制中.0YujCfmUCw

3.2基于变换域地方法

（1）频域体绘制法（FrequencyDomainVolumeRendering）变换域体绘制算法地理论基础是1993年Malzbender提出地傅里叶切片投影

定理.频域体绘制法地基本原理是首先用三维傅立叶变换将空间域地体数据

f（X）变换到频域得到离散频谱F（S），然后沿着经过原点并与视正交地抽取平面对离散频谱F（S）进行插值，插值后地频谱再经过重新采样，得到一个二维地频谱，对其作二维傅立叶反变换即可得到该视方向上地空间域投影图.频域体绘制方法又可细分为傅里叶体绘制和哈特里体绘制.eUts8ZQVRd

（2）基于小波地体绘制法（Wavelet.BasedVolumeRendering）基于小波地体绘制是一种通过对体数据进行三维离散小波变换，以构成体数据地多分辨率表示，然后代入到体绘制方程中来生成三维图像地方法.它包括小波足迹法（WaveletBaseSplatting）和小波域射线投射法（RayTracinginWaveletSpace）.前者地基本思想是将体数据地三维离散小波变换地近似结果直接代入到体绘制方程中求解，从本质上来看它是射线投射法在小波域地实现；

后者地基本思想是利用傅立叶频域绘制先得到每个小波和尺度函数地足迹，再通过小波系数加权得到投影图像.sQsAEJkW5T

4.各种算法地比较分析射线投射法是采用为每个体素分配不透明度和光强地方法来合成图像，因此有利于保留图像地细节，绘制高品质地图像，特别适用于绘制区域特征模糊、体素特征相关性高地三维图像.但是因为需要对每一个体素都进行操作，所以极大地限制了绘制速度.为此人们对射线投射法提出了各种地加速算法，如可以略过三维图像空区域地空间跳跃算法（Space-Leaping）、累积透明度接近于数值1时终止射线投射地射线提前终止法（Early-RayTermination）等.GMsIasNXkA足迹法最大地优点是只有与图像有关地体素才会被映射到象平面，从而可以大大地减少需要处理和存储地数据量.足迹法用经过一个足迹样条地采样均值来代替射线投射法地点采样，这相当于引入了一个低通滤波器对信号进行平滑滤波这样一来，一方面有利于克服图像地失真或混叠，但另一方面有时也会使图像地高频分量受到衰减.足迹法也有应用类似于射线提前终止法概念地加速算法:

基于动态计算屏幕阻塞图（ScreenOcclusionMap）地足迹提前终止法.TIrRGchYzg在错切形变法算法中，体素行程是按不透明度地初分类进行RLE编码地.这就需要在三个主视方向上构造出独立地编码体素；

因为采样插值仅仅出现在剪切出地体素层，因此直接体绘制积分地间距是与视相关地，且不能任意改变来满足沿射线进行密集采样（或过采样），这样在非主视方向上可能不能遵守采样定理，特别是在视域地分辨率高于体素分辨率时，会导致所绘制地图像品质明显下降.7EqZcWLZNX

3D纹理映射地实现需要昂贵地专用图形硬件，对于较大规模地医学体素据则需要在有限地纹理内存中（在较小地机器中通常为4M）反复搬动数据.限制3D纹理影射绘制图像质量地主要因素归结于其画面帧存储器（Framebuffer）有限地位分辨率（8-12bits），这远低于软件算法中使用地浮点数地精度.特别是有限地位分辨率严重制约了不透明度加权颜色（或亮度）方法地使用，因为低不透明度地体素经不透明度加权后，其体素地颜色（或亮度）值会降低到帧存储器分辨率以下，这就限制了具有低不透明度、低密度体素区域地合成绘制.一种解决方法是按比例地放大体素颜色（亮度）值，但这样有时会造成其它区域值地饱和.lzq7IGf02E

频域体绘制法利用快速傅立叶变换（FFT）可以达到较快地绘制速度.由于可以在频域内灵活地根据不同地视对离散频谱抽取平面，因此可以便捷地得到不同视角地图像.但是由于不能保留体素地光吸收特性，因此不能得到半透明地图像.小波域射线投射法是射线投射法在小波域地实现，因此它具有射线投射法地许多优点，如高品质地图像、半透明地视觉效果以及可以加入各种明暗处理等，但缺点是计算量大，绘制速度慢.小波足迹法地最大优点是绘制速度快，可以达

到网络交互级地速率，但缺点是绘制地图像质量呈现一种指数型地自阻塞特

性.zvpgeqJ1hk

表1各种算法地比较分析

体绘制算法

图像品质

绘制速度

算法特点

空

间

域

射线投射法

最高

慢

无需分割，

可以利用不透明度得到整体地层次结构

占用内存大

足迹法

高

中等

占用内存小，可渐进

显示

错切形变法

最快

占用内存小

3D纹理映射

较低

快

图像品质依赖帧存储

器位分辨率

变换域

频域体绘制法

较高

X光片效果，利用FFT算法简洁

小

波

可利用不透明度得到整体地层次

结构，占用内存大

较快

X光片效果,可渐进显示,局部细节添加

5.研究展望

随着现代医学影像技术日新月异地发展和各种医学数字影象设备技术地发展及广泛应用，医学图像三维重建地体绘制技术已经成为研究地热点之一.尽管目前体绘制技术地绘制速度仍然受到计算机运算速度地限制，但随着计算机技术地发展和体绘制算法研究地深入，可以预见体绘制地绘制速度很快能达到实时交互地速率.NrpoJac3v1当前体绘制技术地研究可以根据现有地硬件条件与具体地应用要求，寻求图像品质与绘制速度之间地最佳方案.在各类算法中，基于小波地体绘制技术表现出较好地应用前景，可以作为重点地研究内容，如克服小波足迹法中图像品质地自阻塞指数特性.随着计算机硬件技术地发展，基于硬件地3D纹理映射体绘制也表现

出巨大地发展潜力，许多新技术、新方法有待于深入研究与开发.此外，在不影

响图像品质地前提下如何将错切形变算法中地某些概念引入到射线投射法或足迹法中以提高绘制速度，以及足迹法地明暗算法等都是值得继续深入研究地问题.1nowfTG4KI目前国内地医学图像三维重建地体绘制技术研究依旧处于起步阶段，面对医学临床已具备地丰富地三维医学体数据和应用需求，尽快尽早地开展此类方面地研究具有广泛地应用价值和深远地发展前景.fjnFLDa5Zo

参考文献

[1]UposonC,FaulhaberT.V-Buffer:

Visiblevolumerendering[J].ComputerGraphics,tfnNhnE6e5

1988,22（4）,59-64

[2]WestoverLA.Interactivevolumerendering[A],Proceedingsofthechapelhillonworkshopvolumevisualization[C].ChapelHill:

DepartmentofComputerScience,UniversityofNorthCarolifornia,1989:

9-16.HbmVN777sL

[3]DongF,GordonJ.Volumerenderingoffinedetailswithinmedicaldata[A].In:

IEETVisualization'

2001[C].SanFrancisco:

IEETComputerSocietyPress,2001:

387-394.V7l4jRB8Hs

[4]孙薇薇,张桦.三种体绘制算法地比较[J].天津理工大学学报,2005,21（4）

[5]张尤赛,陈福民.三维医学图像地体绘制技术综述[J].计算机工程与应用,2002,8:

18-1983lcPA59W9

[6]洪锋,梅炯,李明禄.医学图象三维重建技术综述[J].中国图象图形学报,

2003,8

[7]李燕,谭鸥,段会龙.三维医学图象可视化技术综述[J].中国图象图形学报,

2001,6

[8]陈云芝,张彬,程维春.基于Matlab地人体切片地三维重建及可视化[D]

[9]许元飞.三维数据场可视化中体绘制技术地研究[D].西安:

西安科技大学,2005

[10]彭延军,石教英.体绘制技术在医学可视化中地新发展[J].中国图象图形学报（A版）,2002,12（7）:

1239-1246mZkklkzaaP

[11]秦绪佳.医学图象三维重建及可视化技术研究[D].大连理工大学博士论

文.2001

[12]王文举,侯德文.几种变换域体绘制算法地比较研究[J].计算机技术与发展

[J].2008,18（4）

[13]石玉.直接体绘制技术研究[J].电脑知识与技术.2009,5（10）

[14]TotsukaT,LevoyM.Frequencydomainvolumerendering[A].CoputerGraphics

Proceedings,AnnualConferenceSeries[C],1993:

271-27A8VktR43bpw

[15]缪琳.三维数据场地信息处理与体绘制技术[D].2003

[16]LacrouteP.Analysisofapparallelvolumerenderingvolumerenderingsystermbasedonshear-warpfactorization.IEETTransctionVisualizationandcomputerGraphics,1996,

（2）:

218-231ORjBnOwcEd

[17]WestenbergM,Roerdink.FrequencydomainVolumerenderingbythewaveletX-raytransform[J].IEETTransctionsonImageProcessing,2000,9（7）:

1249-12612MiJTy0dTT

[18]孙宇阳.基于单幅图像地三维重建技术综述[J].北方工业大学学报,2011,23

（1）

版权申明本文部分内容，包括文字、图片、以及设计等在网上搜集整理.版权为个人所有

Thisarticleincludessomeparts,includingtext,pictures,anddesign.Copyrightispersonalownership.gIiSpiue7A

用户可将本文地内容或服务用于个人学习、研究或欣赏，以及其他非商业性或非盈利性用途，但同时应遵守著作权法及其他相关法律地规定，不得侵犯本网站及相关权利人地合法权利.除此以外，将本文任何内容或服务用于其他用途时，须征得本人及相关权利人地书面许可，并支付报酬.uEh0U1Yfmh

Usersmayusethecontentsorservicesofthisarticleforpersonalstudy,researchorappreciation,andothernon-commercialornon-profitpurposes,butatthesametime,theyshallabidebytheprovisionsofcopyrightlawandotherrelevantlaws,andshallnotinfringeuponthelegitimaterightsofthiswebsiteanditsrelevantobligees.Inaddition,whenanycontentorserviceofthisarticleisusedforotherpurposes,writtenpermissionandremunerationshallbe