专利名称：扩散张量表面显象的制作方法
技术领域：
本发明涉及图象处理，更具体地说，涉及一种用于扩散张量MRI表面显象的系统和方法。
背景技术：
扩散张量(DT)图象扫描包含至少六个梯度方向，足以确定脑内的扩散张量。根据扩散张量可以确定扩散各向异性测量值、诸如部分各向异性(FA)。而且，扩散张量的主方向能够用于推断作为束的脑白质连通度。
用于这种束的显象策略是在数据点的子集处绘制扩散椭圆体。由于椭圆体的三维场相互封闭，该显象典型的用于数据的二维切片。此外，为了识别每个椭圆体，仅能绘制稀疏网格上的椭圆体。这种类型的显象可能在视觉上变得混乱，并且传达非常少的信息甚至基本上没有价值。
先前用于显象的技术包括基于纹理的方法(texture base approach)，通过在每个点的矢量场方向使源象模糊来生成图象。线积分卷积(LIC)是一种用于实现这种显象的技术。具有基本上类似视觉效果的另一种方法利用偏微分方程(PDE)的解。
用于基于纹理的DT-MRI显象的绘制技术包括平面绘制和体绘制。在平面绘制中，确定DT-MRI数据的平面切片的图象，并显示为纹理化的四边形(参见图1)。体绘制技术包括确定该体积内每个点的纹理，并通过射线造型法或类似的技术绘制该场。这些技术可能在计算上花费比较大。
因此，存在对大比例解剖信息显象的需求。
发明内容
根据本发明实施例，一种扩散张量显象的计算机实现方法，包括接收扩散加权图象切片数据；根据扩散加权图象切片数据分割扩散张量场，以确定三维三角网；以及确定扩散张量场的部分各向异性场和主扩散方向场。该方法包括根据主扩散方向场确定穿过扩散张量场的分割的表面上每个表面点的流线；根据相应的流线确定表面点的定向纹理强度；通过将定向纹理强度、根据主扩散方向场确定的颜色、以及表面点处部分各向异性场的值进行组合来确定表面点的表面颜色；以及通过绘制具有相应表面颜色的表面点对表面点进行显象。
根据本发明实施例，提供一种机器可读程序存储装置，该装置明确包含可由机器执行的程序指令，以执行扩散张量显象的方法步骤。该方法包括接收扩散加权图象切片数据；根据扩散加权图象切片数据分割扩散张量场，以确定三维三角网；以及确定扩散张量场的部分各向异性场和主扩散方向场。该方法包括根据主扩散方向场确定穿过扩散张量场的分割的表面上每个表面点的流线；根据相应的流线确定表面点的定向纹理强度；通过将定向纹理强度、根据主扩散方向场确定的颜色、以及表面点处部分各向异性场的值进行组合来确定表面点的表面颜色；以及通过绘制具有相应表面颜色的表面点对表面点进行显象。
根据本发明实施例，一种扩散张量显象的计算机实现方法，包括接收扩散加权图象切片数据；根据扩散加权图象切片数据分割扩散张量场，以确定三维三角网；以及确定扩散张量场的部分各向异性场和主扩散方向场。该方法包括根据主扩散方向场确定穿过扩散张量场的分割的表面上每个表面点的流线；根据相应的流线确定表面点的定向纹理强度；通过将定向纹理强度、根据主扩散方向场确定的颜色、以及表面点处部分各向异性场的值中的至少两个进行组合来确定表面点的表面颜色；以及通过绘制具有相应表面颜色的表面点对表面点进行显象。
本发明的优选实施例将在下面结合附图进行更加详细的描述
图1是根据本发明实施例的示范性平面LIC。
图2是根据本发明实施例的方法的流程图。
图3是根据本发明实施例的示范性表面显象。
图4是根据本发明实施例的系统的图。
具体实施例方式
根据本发明的实施例，一种方法包括大比例解剖信息和体素方式(voxel-wise)张量信息的DT图象的显象。体素方式张量信息可包括各向异性或方向相关(directionally dependent)信息及主扩散方向信息。
DT核磁共振成象(DT-MRI)是一种测量生物结构中水扩散速率的方法。规则晶格上每个位置处的DT利用对称的正定矩阵、例如3×3矩阵来描述每个体素内水扩散的方向性的体积平均值。通过扩散在白质纤维束区域内是各向异性的观察结果，DT-MRI使得能够从所确定的张量场推断出白质纤维方向。这些纤维束能够在三维内显象为流线或流管。
流管是流线的三维类似物。通过使流线作为管的中心线确定流管。流管直径用于对关于显象的张量场的附加信息、诸如部分各向异性(FA)值进行编码。
参考图2，该方法包括步骤200，输入图象数据；步骤201，根据输入图象数据分割DT场；步骤202，确定FA场和主扩散方向(PDD)场。该方法包括步骤203，根据PDD场，确定穿过所分割的表面上的每个点(例如如果显象可以看到的每个点)的流线。该方法进一步包括步骤204，根据流线确定线积分卷积(LIC)强度；步骤205，通过对LIC强度、根据PDD场确定的颜色和该点的FA值进行组合来确定每个表面点的表面色；步骤206，通过绘制表面上具有其相应确定的颜色的每个点绘制该表面。
方框200输入图象数据可以是记录到张量场的扩散加权MRI图象切片或CT图象切片。这里，结合MRI图象切片呈现本发明的示范性实施例。
参考方框201，利用产生三维三角网格的测量活动轮廓线(geodesicactive contour)对DT-MRI捕获的DT场进行分割。分割的另一种方式包括例如纤维追踪方法、基于PDD的直接阈值、扩散的标量测量等等。
DT场分割的示例包括将张量转化为特征值/特征向量表示，这是旋转不变量。所转化的张量能够解释为具有沿特征向量取向的主轴、半径等于相应特征值的椭圆体。该椭圆体描述了水分子在固定扩散时间之后的概率分布。利用特征值/特征向量能够确定包括FA的各向异性测量值。
参考方框202，确定FA场。FA场为灰度级图象，例如参见图3中的302，这使白质变暗使其区别于灰质和脑内的脑室或脊柱扫描。FA是DT大小的部分的测量值，其能够归因于各向异性扩散。只要对DT进行了预算(见方框201)，就能够通过得出对应于当前张量的主特征值的特征向量确定PPD。通过对穿过特征向量的场的流线进行显象可描绘纤维束。描绘包括通过RGB配色方案对PPD进行编码，其中，例如，利用绿色对沿前后方向的纤维束进行编码，可用红色对从左向右延伸的纤维束进行编码，垂直于图象平面的纤维由蓝色编码。
在方框203，为三角形网格的分割的表面上的点确定流线。流线沿每个点的DT的主特征向量的方向。
参考方框204，通过LIC将对应于扩散张量的主特征值的三维向量场显象为基于纹理的向量。LIC通过将噪声纹理和与穿过每个象素的流线一致的曲线核(curvilinear kernel)O进行卷积生成强度灰度级值(例如在核内确定流线)。
能够利用LIC技术将对应于扩散张量主特征值的三维向量场显象。LIC技术显象本地向量场的切线。
根据LIC，通过将噪声纹理和与穿过每个象素的流线一致的曲线核进行卷积生成强度值(例如在核内确定流线)，诸如通过I(x0)=∫s0-Ls0+LT(σ(s))k(s0-s)ds]]>这里I(x0)是象素x0处LIC纹理的强度。k是宽度2L的滤波器核(filterkernel)。T是输入噪声纹理，σ是穿过点x0的流线。给定特征向量的离散场，流线σ能够通过数值积分得出。
LIC结果是具有相同流线上相邻象素之间的高相关值和不共享流线的象素的对比值的灰度级纹理。例如，由于扩散场在低FA值不再具有主方向，因此在特定阈值以下的FA值可为积分的停止标准。
由于纤维方向平行于扩散张量的主特征向量，纤维路径(fiber path)可确定为主特征向量场的积分曲线。停止标准基于FA值。例如，当FA下降低于大约0.17时，认为扩散几乎是各向同性的，不再跟踪纤维。
参考方框205，表面颜色由两个或两个以上LIC灰度级强度、根据PDD场确定的颜色以及每个点或象素处FA灰度级值的组合确定，从而根据PDD使DT的可见点(例如表面)的向量场获得颜色。LIC、FA和PDD的组合可以使用乘法、例如将PDD的RGB值、LIC的灰度值和FA的灰度值相乘从而得出绘制包括关于LIC和FA信息的表面的修正颜色值。组合可由PPD和LIC或PDD和FA限定。
在方框206，以颜色绘制出表面点，并可输出到显示器或存储介质。
根据本发明的实施例，DT显象(例如，参见图3)的方法O将测定和显象限制到关注的、解剖的或其他方面的表面。在利用水平集合技术的体积内或在三角网格表面上，均能够以这种方式求解LIC和PDE方法。结果可能为由将三角网格纹理化而绘制的表面或通过射线造型法绘制的体积。因为生成整个扫描体积的子集(例如已确定的表面)的纹理，这种方法可使计算节省。例如，图3表示胼胝体301(哺乳动物脑中连接左右脑半球的白质结构)的段，用1024×1024纹理描绘建立纹理，其中LIC计算为1048576个点，其中白质大约为10个切片厚，并且需要确定512×512纹理O的10个切片(2621440个点)。在图3中，二维FA图象302用作胼胝体301表面的参考。诸如F0的其他图象可用作参考。
应理解，当以灰度级说明图3时，可采用使用LIC的灰度级值和PA值修正的PDD颜色值的颜色绘制图3的扩散张量。
此外，显象使得用户可以更加清楚地看到组织类别之间的分界线，并且看到界面处的扩散。
可以理解的是，本发明可通过不同形式的硬件、软件、固件、专用处理器或其组合实现。在一个实施例中，本发明可通过作为明确包含在程序存储装置中的应用程序软件实现。该应用程序可上传到包含任何适合的体系结构的机器并由其执行。
参考图4，根据本发明的实施例，用于实现扩散张量表面显象方法的计算机系统401可包括中央处理单元(CPU)402、存储器403和输入/输出(I/O)接口404。计算机系统401通常通过I/O接口404与显示器405及诸如鼠标和键盘的各种输入装置406相连。辅助电路可包括诸如高速缓冲存储器、电源、时钟脉冲电路以及通信总线的电路。存储器403可包括随机存取存储器(RAM))、只读存储器(ROM))、磁盘驱动器、磁带驱动器等或其组合。本发明可实现为存储在存储器403内并由CPU 402执行以处理来自信号源408的信号的例程407。同样地，当执行本发明的例程407时，计算机系统401是变为专用计算机系统的通用计算机系统。
计算机平台401还包括操作系统和微指令代码。在此描述的各种处理和功能可以是通过操作系统执行的微指令代码部分或应用程序部分(或其组合)。此外，各种其他外围设备、诸如附加数据存储装置和打印装置可连接到计算机平台。
可以进一步理解的是，因为附图中描述的一些组成系统组件和方法步骤可通过软件实现，所以系统组件(或过程步骤)之间的实际连接可根据对本发明编程的方式而不同。在本文提供的本发明的论述的条件下，本领域技术人员能够想到这些及类似的实现或配置。
已经描述了用于扩散张量表面显象的系统和方法，应注意的是，本领域技术人员能够根据上面的论述作出更改和变化。因此应该理解，对本发明实施例作出的改变在本发明的范围和精神内。
权利要求
1.一种用于扩散张量显象的计算机实现方法，包括接收扩散加权图象切片数据；根据所述扩散加权图象切片数据分割扩散张量场，以确定三维三角网格；确定所述扩散张量场的部分各向异性场和主扩散方向场；根据所述主扩散方向场确定穿过所述扩散张量场的分割的表面上每个表面点的流线；根据相应的流线确定所述表面点的定向纹理强度；通过将所述定向纹理强度、根据所述主扩散方向场确定的颜色、以及所述表面点处部分各向异性场的值进行组合来确定所述表面点的表面颜色；以及通过绘制具有相应表面颜色的所述表面点对所述表面点进行显象。
2.如权利要求1所述的计算机实现方法，其中不是确定点而是确定所述表面点的流线。
3.如权利要求1所述的计算机实现方法，其中根据所述流线的方向以配色方案对所述主扩散方向场编码。
4.如权利要求1所述的计算机实现方法，其中所述定向纹理强度是线积分卷积强度。
5.如权利要求4所述的计算机实现方法，其中所述线积分卷积通过将噪声纹理和与穿过每个象素的所述流线一致的曲线核进行卷积生成强度值。
6.如权利要求5所述的计算机实现方法，其中所述线积分卷积根据以下公式确定I(x0)=∫s0-Ls0+LT(σ(s))k(s0-s)ds]]>这里I(x0)是点x0处线积分卷积纹理的强度，k是宽度2L的滤波器核，T是噪声纹理，σ是穿过所述点x0的流线，其中给定特征向量的离散场，所述流线σ能够通过数值积分得出。
7.如权利要求1所述的计算机实现方法，其中确定所述表面颜色包括将每个表面点的所述定向纹理强度、根据所述主扩散方向场确定的所述颜色以及部分各向异性场的值相乘。
8.一种机器可读程序存储装置，所述装置明确包含可由所述机器执行的程序指令，以执行扩散张量显象的方法步骤，所述方法步骤包括接收扩散加权图象切片数据；根据所述扩散加权图象切片数据分割扩散张量场，以确定三维三角网格；确定所述扩散张量场的部分各向异性场和主扩散方向场；根据所述主扩散方向场确定穿过所述扩散张量场的分割的表面上每个表面点的流线；根据相应的流线确定所述表面点的定向纹理强度；通过将所述定向纹理强度、根据所述主扩散方向场确定的颜色、以及所述表面点处部分各向异性场的值进行组合来确定所述表面点的表面颜色；以及通过绘制具有相应表面颜色的所述表面点对所述表面点进行显象。
9.如权利要求8所述的方法，其中不是确定点而是确定所述表面点的流线。
10.如权利要求8所述的方法，其中根据所述流线的方向以配色方案对所述主扩散方向场编码。
11.如权利要求8所述的计算机实现方法，其中所述定向纹理强度是线积分卷积强度。
12.如权利要求11所述的计算机实现方法，其中所述线积分卷积通过将噪声纹理和与穿过每个象素的所述流线一致的曲线核进行卷积生成强度值。
13.如权利要求12所述的计算机实现方法，其中所述线积分卷积根据以下公式确定I(x0)=∫s0-Ls0+LT(σ(s))k(s0-s)ds]]>这里I(x0)是点x0处线积分卷积纹理的强度，k是宽度2L的滤波器核，T是噪声纹理，σ是穿过所述所述点x0的流线，其中给定特征向量的离散场，所述流线σ能够通过数值积分得出。
14.如权利要求8所述的计算机实现方法，其中确定所述表面颜色包括将每个表面点的所述定向纹理强度、根据所述主扩散方向场确定的所述颜色以及部分各向异性场的值相乘。
15.一种用于扩散张量显象的计算机实现方法，包括接收扩散加权图象切片数据；根据所述扩散加权图象切片数据分割扩散张量场，以确定三维三角网格；确定所述扩散张量场的部分各向异性场和主扩散方向场；根据所述主扩散方向场确定穿过所述扩散张量场的分割的表面上每个表面点的流线；根据相应的流线确定所述表面点的定向纹理强度；通过将所述定向纹理强度、根据所述主扩散方向场确定的颜色、以及所述表面点处部分各向异性场的值中的至少两个进行组合来确定所述表面点的表面颜色；以及通过绘制具有相应表面颜色的所述表面点对所述表面点进行显象。
16.如权利要求15所述的计算机实现方法，其中不是确定点而是确定所述表面点的流线。
17.如权利要求15所述的计算机实现方法，其中所述定向纹理强度是线积分卷积强度。
全文摘要
一种用于扩散张量显象的方法，包括分割扩散张量场(201)；确定部分各向异性场和主扩散方向场(202)；利用主扩散方向场确定穿过分割的表面上每个可见点的流线(203)；根据流线确定定向纹理强度(204)；通过将定向纹理强度、根据主扩散方向场确定的颜色、以及该点处部分各向异性的值进行组合来确定表面颜色(205)；通过绘制具有相应确定颜色的表面上的每个点对表面进行显象(206)。
文档编号A61B5/055GK1818707SQ200610067349
公开日2006年8月16日 申请日期2006年2月9日 优先权日2005年2月9日
发明者T·麦格劳 申请人:西门子共同研究公司