CN101903912A - 使用多强度重新分配函数进行绘制 - Google Patents

使用多强度重新分配函数进行绘制 Download PDF

Info

Publication number
CN101903912A
CN101903912A CN200880121378XA CN200880121378A CN101903912A CN 101903912 A CN101903912 A CN 101903912A CN 200880121378X A CN200880121378X A CN 200880121378XA CN 200880121378 A CN200880121378 A CN 200880121378A CN 101903912 A CN101903912 A CN 101903912A
Authority
CN
China
Prior art keywords
voxel
value
classification
projection
redefines
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN200880121378XA
Other languages
English (en)
Inventor
H·莱曼
J·威斯
D·格勒
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips Electronics NV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Koninklijke Philips Electronics NV filed Critical Koninklijke Philips Electronics NV
Publication of CN101903912A publication Critical patent/CN101903912A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/08Volume rendering

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Generation (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

本发明涉及一种用于使用射线投射方法对包括多个体素的图像数据集进行可视化的系统(100),所述多个体素中的每个体素属于至少一个分类,所述至少一个分类中的每个分类与用于根据体素的所测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联,系统包括采样单元(120),其用于基于与来自图像像素的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算所述投影射线投射上的所述样本定位处的样本值,其中,使用与所述至少一个体素的至少一个分类相关联的强度重新分配函数,根据所述至少一个体素的所测得体素值计算所述至少一个体素的重定义体素值。

Description

使用多强度重新分配函数进行绘制
技术领域
本发明涉及体积图像数据可视化领域,并且更特别地,涉及使用射线投射将图像数据体积中包括的多个对象进行可视化的领域。
背景技术
诸如,例如,最大强度投影(MIP)、数字重建射线照相(DRR)或直接体积绘制(DVR)的射线投射图像绘制方法通常用于将CT或MRI体积图像可视化。通常,在分段图像数据体积中包括多个对象。例如肋骨的一些对象可以妨碍例如心脏的另一对象的视图。
在此被引用作为参考1的I.Viola、A.Kanitsar和M.E.Groeller的文章“Importance-driven volume rendering”(IEEE Visualization 2004 October10-15,Austin,Texas,USA,139-145页)描述了对分割图像数据中的对象进行可视化的方法。该方法采用修改的DVR,其中,在图像数据中描述的对象的重要性由重要性指数进行限定。在根据图像数据计算的图像中,布置具有较高重要性指数的对象从而使其比具有较低重要性指数的对象更为可视化。这通过将不同水平的稀疏性(sparseness)与不同重要性指数相关联以及通过使用重要性组合达到。在DVR中,传递函数将颜色和不透明性分配给体积数据的体积内的每个样本。沿每条从图像平面到图像数据体积内投射的射线的样本值被组合并计算最终图像。对象的稀疏性由分配给对象样本的不透明性限定。重要对象的样本为不透明的而次重要对象的样本为半透明或全透明的。在参考1中描述的定义稀疏性的水平的方式包括不透明性和/或颜色调制、屏蔽门透明性以及体积稀释。
发明内容
本发明提供了一种使用射线投射的新的、非常直观并且容易施加的方法对图像数据体积中包括的多个对象进行可视化。
在本发明的一方面中,提供了使用射线投射方法,用于对包括多个体素的图像数据集进行可视化的系统,多个体素中的每个体素属于至少一个分类,所述至少一个分类中的每个分类与用于根据体素的所测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联,所述系统包括采样单元,其用于基于与从图像像素投射的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算所述投影射线上的样本定位处的样本值,其中,使用与至少一个体素的至少一个分类相关联的强度重新分配函数,根据至少一个体素的所测得体素值计算所述至少一个体素的重定义体素值。
例如,描述骨骼的第一类体素可以被赋予生成在图像中可见的骨骼的低重定义的灰度值的第一强度重新分配函数。在另一方面,描述血管的第二分类的体素可以被赋予生成在图像中可见的血管的高灰度值的第二强度重新分配函数。例如,图像可以使用MIP计算。因此,检查图像的用户可以在骨骼结构在图像中不可见的情况下看到血管。有利地,图像中要被可视化的感兴趣区域可以由组织类型体素分类或由例如,基于图像数据分割从体素分类获得的其他信息限定。例如,冠状动脉可以具有不同于肺静脉的强度重新分配函数的强度重新分配函数。
本领域技术人员将理解体素值表示体素强度。存在用于表达体素值的不同尺度和单位,所述尺度和单位包括,但不限于,Hounsfield单位(HU)和由[0,255]范围内的整数表示的灰度值。可以使用任意合适的体素强度限定强度重新分配函数。
在系统的实施例中,采样单元包括:
-定位单元,其用于选择从图像像素投射的投影射线上的样本定位;
-体素单元,其用于选择邻近样本定位的至少一个体素;
-重定义单元,其用于使用与至少一个体素的至少一个分类相关联的强度重新分配函数,根据至少一个体素的所测得体素值计算至少一个体素的重定义体素值;以及
-组合单元,其用于基于投影射线上的至少一个体素的重定义体素值,计算在样本定位处的样本值。这些单元表示采样单元的有效实施方式。
在实施例中,系统还包括用于基于用户输入,对强度重新分配函数进行成形的重新分配单元。经改变的强度重新分配函数可以由系统施加给图像数据并且可以实时计算新的图像。所计算的图像可以在显示器上进行显示。用户因此能够交互地重定义分类的体素的灰度值,以达到所观察的图像数据集的最佳可视化。
在系统的实施例中,射线投射方法为最大值强度投影或最小值强度投影。所述最大值强度投影和最小值强度投影为普遍的绘制技术,并且大多数放射学家知道如何使用这一绘制技术解释所绘制的图像。
在系统的实施例中,至少一个分类包括背景分类。所有不能被分类为结构或组织的体素的体素可以被分类为背景分类体素。
在系统的实施例中,至少一个分类包括多个分类。这一实施例有助于处理体素不能被唯一地分类为仅一个分类中的体素的情况。基于体素的多个重定义的灰度值计算从图像像素投射的投影射线上的样本定位处的样本值,其中,使用与体素的不同分类相关联的不同强度重新分配函数,根据体素的所测得体素值计算体素的每个重定义的体素值。
在实施例中,系统还包括用于确定至少一个体素的至少一个分类的分类单元。在实施例中,分类单元可以采用体素分类器。在另一实施例中,分类单元可以采用图像数据分割。
在本发明的另一方面中,提供了一种使用射线投射方法对包括多个体素的图像数据集进行可视化的方法,所述多个体素中的每个体素属于至少一个分类,所述至少一个分类中的每个分类与用于根据体素的所测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联,方法包括采样步骤,其用于基于与从图像像素投射的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算所述投影射线上的样本定位处的样本值,其中,使用与至少一个体素的至少一个分类相关联的强度重新分配函数,根据至少一个体素的所测得体素值计算至少一个体素的重定义体素值。
在本发明的另一方面中,提供了一种要由计算机布置加载的计算机程序产品,所述计算机程序产品包括用于使用射线投射方法,对包括多个体素的图像数据集可视化的指令,所述多个体素中的每个体素属于至少一个分类,所述至少一个分类中的每个分类与用于根据体素的所测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联,所述计算机布置包括处理单元和存储器,在被加载之后,计算机程序产品为所述处理单元提供执行以下任务的能力:基于与从图像像素投射的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算投影射线上的样本定位处的样本值,其中,使用与所述至少一个体素的至少一个分类相关联的强度重新分配函数,根据所述至少一个体素的所测得体素值计算所述至少一个体素的重定义体素值。
在本发明的另一方面,根据本发明的系统包括在图像采集装置中。
在本发明的另一方面,根据本发明的系统包括在工作站中。
本领域技术人员将理解可以以任意被认为是有用的方式组合的两个或更多的上述实施例、实施方式和/或本发明的各方面。
本领域技术人员基于本说明书可以实现对应于所描述的系统的修改和变型的图像采集装置、工作站、方法和/或计算机程序产品的修改和变型。
本领域技术人员将理解本方法可以被应用于由各种采集模态采集的多维图像数据,例如,三维图像或四维图像,所述采集模式诸如,但不限于,标准X-射线成像、计算机断层摄影(CT)、磁共振成像(MRI)、超声(US)、正电子发射断层摄影(PET)、单光子发射计算机断层摄影(SPECT)和核医学(NM)。
附图说明
从下面描述的实施方式和实施例并结合附图,本发明的这些方面和其他方面将变得明显并得以阐明,在附图中:
图1示意性地示出了系统的示例性实施例的方框图;
图2示出了与一类体素相关联的示例性强度重新分配函数;
图3示出了根据本发明的使用应用于CT心脏数据的MIP绘制对心脏进行可视化;
图4示出了根据本发明的使用应用于上述CT心脏数据的MIP绘制对冠状动脉进行可视化;
图5示出了方法的示例性实施方式的流程图;
图6示意性地示出了图像采集装置的示例性实施例;
图7示意性地示出了工作站的示例性实施例。
附图中相同的附图标记用于表示类似的部件。
具体实施方式
图1示意性地示出了用于使用射线投射技术将包括多个体素的图像数据集可视化的系统100的示例性实施例的方框图,多个体素中的每个体素属于至少一个分类,所述至少一个分类中的每个分类与用于根据体素的所测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联,系统包括采样单元120,所述采样单元120用于基于与从图像像素投射的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算所述投影射线上的样本定位处的样本值,其中,使用与至少一个体素的至少一个分类相关联的强度重新分配函数,根据至少一个体素的所测得体素值计算至少一个体素的重定义体素值。
系统100的示例性实施例的采样单元120任选地包括:
-定位单元122,其用于选择从图像像素投射的投影射线上的样本定位;
-体素单元124,其用于选择与样本定位邻近的至少一个体素;
-重定义单元126,其用于使用与至少一个体素的至少一类相关联的强度重新分配函数,根据至少一个体素的所测得体素值计算至少一个体素的重定义体素值;以及
-组合单元128,其用于基于投影射线上的至少一个体素的重定义体素值计算样本定位处的样本值。
系统100的示例性实施例还包括下面单元:
-分类单元110,其用于确定至少一个体素的至少一个分类;
-重新分配单元130,其用于基于用户输入,形成强度重新分配函数;
-图像单元140,其用于基于从所述图像像素投射的投影射线上的样本定位处的样本值计算图像像素的图像像素值;
-控制单元160,其用于控制系统100中的工作流程;
-用户接口165,其用于与系统100的用户交互;以及
-存储器单元170,其用于存储数据。
在系统100的实施例中,存在用于输入数据的三个输入连接器181、182和183。第一输入连接器181布置为接收来自数据存储器件的数据,所述数据存储器件诸如,但不限于,硬盘、磁带、闪存或光盘。第二输入连接器182布置为接收来自用户输入设备的数据,所述用户输入设备诸如,但不限于,鼠标或触摸屏。第三输入连接器183布置为接收来自诸如键盘的用户输入设备的数据。将输入连接器181、182和183连接到输入控制单元180。
在系统100的实施例中,存在用于输出数据的两个输出连接器191和192。第一输出连接器191布置为将数据输出到数据存储器件,所述数据存储器件诸如硬盘、磁带、闪存或光盘。第二输出连接器192布置为将数据输出到显示设备。输出连接器191和192经由输出控制单元190接收相关数据。
本领域技术人员将理解有许多方式将输入设备连接到输入连接器181、182和183,以及将输出设备连接到系统100的输出连接器191、192。这些方式包括,但不限于有线和无线连接,诸如但不限于局域网(LAN)和广域网(WAN)的数字网络、因特网、数字电话网络以及模拟电话网络。
在系统100的实施例中,系统100包括存储器单元170。系统100布置为经由输入连接器181、182和183中的任一个接收来自外部设备的输入数据,并将所接收到的输入数据存储在存储器单元170中。将输入数据加载到存储器单元170允许由系统100的单元对相关数据部分进行快速存取。输入数据可以包括,例如,图像数据集和强度重新分配函数,每个分类具有一个函数的体素分类方案。存储器单元170可以由以下设备实施,所述设备诸如,但不限于,随机存取存储器(RAM)芯片、只读存储器(ROM)芯片和/或硬盘驱动器和硬盘。存储器单元170还被布置为存储输出数据。输出数据可以包括,例如,根据本发明计算的图像。存储器单元170还布置为经由存储器总线175接收来自系统100的单元的数据和/或将数据递送到系统100的单元,所述系统100包括分类单元110、作用单元120、图像单元130、控制单元160和用户接口165。存储器单元170还布置为经由输出连接器191和192的任一个使得外部设备可以获得输出数据。将来自系统100的单元的数据存储在存储器单元170中可以有利地提高系统100的单元的性能并提高将输出数据从系统100的单元到外部设备的传输速率。
可替代地,系统100可以不包括存储器单元170和存储器总线175。由系统100使用的输入数据可以由连接到系统100的单元的诸如外部存储器或处理器的至少一个外部设备提供。类似地,由系统100生成的输出数据可以被提供到连接到系统100的单元的诸如外部存储器或处理器的至少一个外部设备。系统100的单元可以用于经由内部连接或者经由数据总线接收来自彼此的数据。
在系统100的实施例中,系统100包括用于控制系统100中的工作流的控制单元160。控制单元可以布置为从系统100的单元接收控制数据并将控制数据提供给系统100的单元。例如,在计算投影射线上的一个样本定位处的样本值之后,采样单元120可以布置为向控制单元160提供控制数据“所计算的样本值”,并且控制单元160可以布置为向采样单元120提供控制数据“确定投影射线上的下一样本定位处的样本值”。确定下一定位或下一投影射线可以由图像单元130、控制单元160或控制单元120执行。控制函数可以在系统100的任一单元中执行。
在系统100的实施例中,系统100包括用于与系统100的用户通信的用户接口165。用户接口165可以布置为接收强度重新分配函数的用户定义。用户接口可以还提供用于旋转图像数据集以计算对于指定路径有用的不同视图的器件。用户接口还可以向用户提供例如,属于体素分类的体素的柱状图的形式的信息,用于在显示器上显示。任选地,用户接口还可以接收用于选择系统的操作的模式的用户输入,诸如,例如,用于选择图像绘制技术。本领域技术人员将理解更多功能可以有利地在系统100的用户接口165中实施。
在系统100的实施例中,输入数据包括如下图像数据集,其中,每个体素包括体素定位、体素值,即:体素强度,以及体素分类。采样单元120适于计算用于计算投影射线上的给定定位的采样值每个体素的重定义值。使用对应于相应体素的分类的强度重新分配函数根据所测得体素值计算重定义体素值。接下来,采样单元120适于使用所计算的重定义体素值计算在投影射线上的给定定位的样本值。可以使用不同采样技术,所述采样技术包括,但不限于,最近邻近、三线性、高斯或三次样条。
在实施例中,控制单元160布置为确定从图像像素投射的射线,以及用于确定每条射线上的样本定位。相邻样本定位之间的距离可以是相同的。
可以基于诸如,例如,骨、血液、肌肉的组织类型和/或基于由诸如,例如股骨、肋骨、肺、心脏的体素表示的结构定义。可以针对每个分类对强度重新分配函数进行预定义并基于其分类由系统100自动应用到所处理的体素。
根据数据分类潜在的假定,体素可以同时属于若干组织分类。在实施例中,体素的体素值对应于体素体积中所包括的若干组织类型的积累的密度,并且体素分类表示积累的密度有多少属于哪一组织类型。针对每个体素限定分类向量,并且每个向量分量表示体素值中的多少属于与分类向量内的所述分量的位置对应的组织类型。因此,能够计算不同组织类型的作用的总和。特定组织类型对所测得体素值的作用为所测得体素值与对应于特定组织类型的向量分量的乘积。这一测得的作用用于使用对应于特定组织类型的分类的强度重新分配函数,计算对应于特定组织类型的分类的重定义作用。
在系统100的实施例中,图像数据集包括多个体素,其中,每个体素包括体素坐标、体素值和组织类型向量,每个向量分量描述组织类型t的权重ct。针对每个组织类型t,提供强度重新分配函数ft。使用标准射线投射处理计算像素强度,即:对于图像平面中的每个像素p,射线被在观察方向投射,并且在投射射线上的每个样本定位i,例如通过邻近体素的最近邻或三线性内插从体积数据获取样本值si和组织类型向量
Figure BPA00001160461800081
基于样本值和组织类型向量,每个组织类型t在每个样本定位i的作用可以被写为:
w i t = s i · c i t
权重
Figure BPA00001160461800083
满足下述条件:
Σ t c i t ≤ 1
因此,在样本定位i处未被分类的强度的权重由下式给出
c i = 1 - Σ t c i t
这一权重ci被称为背景分类权重。背景分类对在每个样本定位i处的样本值的作用为:
wi=si·ci
背景分类被分配有背景强度重新分配函数f。在样本定位处的样本值vi由下式给出:
v i = max ( max t f t ( w i ) , f ( w i ) )
可替代地,样本定位处的样本值vi可以被定义为所有作用的总和:
v i = Σ t f t ( w i t ) + f ( w i )
在最大值强度投影中,像素值vp为从所述像素投射的射线上的所有采样值vi的最大值:
v p = max i ( v i )
可替代地,体素值可以被计算为从所述像素投射的射线上的所有采样值vi的平均:
v p = 1 N Σ i v i
其中,N表示样本的数量。
在实施例中,系统100包括用于基于用户输入对强度重新分配函数成形的重新分配单元130。图2示出了与一类体素相关联的示例性强度重新分配函数。在用于在显示器上显示的窗20中,用户可以通过画出强度重新分配函数的曲线21定义强度重新分配函数。以Hounsfield单位(HU)表达所测得体素值。将重定义体素值表达为灰度级值。曲线可以被实现为例如,由用户放置的若干控制点22控制的折线或贝塞尔曲线。有利地,可以实时使用用户定义的强度重新分配函数以计算用于在例如,图像窗口的显示器上显示的图像。用户可以基于来自所显示的图像的反馈交互地继续调节强度重新分配函数。任选地,窗口20可以布置为显示与所成形的强度重新分配函数对应的分类的体素的体素值的柱状图。任选地,重新分配单元130可以包括在用户接口165中。
在系统100的实施例中,系统的图像单元140使用最大值强度投影(MIP)技术用于图像绘制。对于来自图像平面内的像素的每个射线投影,像素值为沿该射线的最大样本值。图3示出了根据本发明使用应用于CT心脏数据的MIP技术对心脏进行可视化。左侧的图像31示出了CT心脏数据集的标准MIP。右侧的图像32示出了相同数据集的相同视图,其中,应用三个不同线性强度重新分配函数,两个分别用于突出左心室的腔和心肌,而一个用于对提供有助于对图像进行取向的参照结构的一些组织进行可视化。图4示出了根据本发明的使用应用于上述CT心脏数据的MIP绘制对冠状动脉进行可视化。左侧的图像41示出了上述CT心脏数据集的另一标准MIP视图。右侧的图像42示出了相同数据集的相同视图,其中,通过施加适当的强度重新分配函数抑制七个不同区域,每个区域对不同组织类型进行分类。因此,图像生成冠状动脉中的一个的自由视图,不能得到针对其的分割数据。
本领域技术人员将理解任意适合的技术可以被用于计算体素值,所述适合技术诸如,但不限于,用于生成直接体积绘制的技术、最靠近的导管投影或数字重建放射照片。
本领域技术人员将理解系统100可以为用于在其工作的许多方面辅助医生的有用工具。
本领域技术人员还将理解系统100的其他实施例也是可以的。尤其地,可以重定义系统的单元,并且重新分配其功能。尽管所描述的实施例应用于医学图像,但是系统的与医学应用无关的其他应用也是可能的。
可以使用处理器实施系统100的单元。通常地,其功能在软件程序产品的控制下执行。在执行期间,软件程序产品通常被加载到如RAM的存储器中,并从这里执行。程序可以从背景存储器中加载,所述背景存储器诸如ROM、硬盘或磁存储器和/或光学存储器,或者可以经由如因特网的网络进行加载。任选地,专用集成电路可以提供所描述的功能性。
图5示出了使用射线投射方法,对包括多个体素的图像数据集进行可视化的方法500的示例性实施方式的流程图。多个体素中的每个体素数据属于至少一个分类。所述至少一个分类中的每个分类与用于根据体素的所测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联。方法从用于基于图像数据集确定初始射线和初始样本定位的初始化步骤502开始。在初始化步骤502之后,方法继续到采样步骤520,所述采样步骤520用于基于与从图像像素投射的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算所述投影射线上的样本定位处的样本值,其中,使用与至少一个体素的至少一个分类相关联的强度重新分配函数,根据至少一个体素的所测得体素值计算所述至少一个体素的重定义体素值。这在后续子步骤中执行。在定位步骤522,选择从图像像素投射的投影射线上的样本定位。在体素步骤524,从图像数据集中选择邻近样本定位的至少一个体素。在重定义步骤526,使用与至少一个体素的至少一个分类相关联的强度重新分配函数根据至少一个体素的所测得体素值计算所述至少一个体素的重定义体素值。在组合步骤528,基于投影射线上的至少一个体素的重定义体素值计算样本定位处的样本值。在采样步骤520之后,方法继续道定位更新步骤532用于更新射线上的样本定位。在定位更新步骤532之后,方法500继续到采样步骤520,或者,如果已经计算射线上的所有样本定位处的样本值,方法500继续到图像步骤540用于基于从所述图像像素投射的投影射线上的样本定位处的样本值计算图像像素的图像像素值。在图像步骤540之后,方法500继续到射线更新步骤534,所述更新步骤534用于选择下一像素和从该像素投射的下一射线,以及该射线上的初始样本定位。在射线更新步骤534之后,方法继续到采样步骤520,或者,如果图像平面中的所有像素的像素值都已经被计算,方法结束。
本领域技术人员在不脱离有本发明提出的概念的情况下,可以使用线程模型、多处理器系统或多个处理改变一些步骤的顺序或者连续地执行一些步骤。任选地,本发明的方法的两个或更多步骤可以被合并到一个步骤。任选地,本发明的方法的步骤可以被分为多个步骤。
图6示意性地示出了采用系统100的图像采集装置600的示例性实施例,所述图像采集系统600包括经由内连接与系统100连接的CT图像采集单元、输入连接器601、输出连接器602。这一布置有利地提高了图像采集装置600的性能,为所述图像采集装置600提供系统100的有利性能。
图7示意性地示出了工作站700的示例性实施例。工作站包括系统总线701。处理器710、存储器720、磁盘输入/输出(I/O)适配器730以及用户接口(UI)740被可操作地连接到系统总线701。磁盘存储器设备731可操作地耦合到磁盘I/O适配器730。键盘741、鼠标742和显示器743被可操作地耦合到UI 740。将被实施为计算机程序的本发明的系统100存储在磁盘存储器设备731中。工作站700布置为将程序和输入数据加载到存储器720并在处理器710上执行该程序。用户可以使用键盘741和/或鼠标742向工作站700输入信息。工作站布置为向显示设备743和/或磁盘731输出信息。本领域技术人员将理解存在本领域已知的工作站700的若干其他实施例,并且本实施例用于示出本发明,而不应该被理解为将本发明限制于该特定实施例。
应该注意上述实施例示出而非限制本发明,并且本领域技术人员在不脱离所附权利要求的范围的情况下能够设计可替代实施例。在权利要求中,置于括号之间的任何附图标记都不应该被理解为限制该权利要求。词语“包括”不排除权利要求或说明书中未列出的元件或步骤的存在。元件前的词语“一”或“一个”不排除存在多个这样的元件。本发明可以借助于包括若干不同元件的硬件实施和借助于编程计算机实施。在列举若干单元的系统权利要求中,若干这些单元可以由一个并且是相同的硬件或软件示例。使用词语第一、第二、第三等不表示任何顺序。应该将这些词语理解为名称。

Claims (10)

1.一种用于使用射线投射方法对包括多个体素的图像数据集进行可视化的系统(100),所述多个体素中的每个体素属于至少一个分类,所述至少一个分类中的每个分类与用于根据体素的测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联,所述系统包括采样单元(120),其用于基于与从图像像素投射的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算所述投影射线上的所述样本定位处的样本值,其中,使用与所述至少一个体素的至少一个分类相关联的强度重新分配函数,根据所述至少一个体素的测得体素值计算所述至少一个体素的所述重定义体素值。
2.根据权利要求1所述的系统(100),其中,所述采样单元(120)包括:
-定位单元(122),其用于选择从所述图像像素投射的所述投影射线上的所述样本定位;
-体素单元(124),其用于选择邻近所述样本定位的所述至少一个体素;
-重定义单元(126),其用于使用与所述至少一个体素的至少一个分类相关联的强度重新分配函数,根据所述至少一个体素的测得体素值计算所述至少一个体素的重定义体素值;以及
-组合单元(128),其用于基于所述投影射线上的所述至少一个体素的重定义体素值,计算在所述样本定位处的所述样本值。
3.根据权利要求1所述的系统(100),还包括用于基于用户输入对所述强度重新分配函数进行成形的重新分配单元(130)。
4.根据权利要求1所述的系统(100),其中,所述射线投射方法为最大强度投影或最小强度投影。
5.根据权利要求1所述的系统(100),其中,所述多个分类包括背景分类。
6.根据权利要求1所述的系统,其中,所述至少一个分类包括多个分类。
7.一种使用射线投射方法对包括多个体素的图像数据集进行可视化的方法(500),所述多个体素中的每个体素属于至少一个分类,所述至少一个分类中的每个分类与用于根据体素的测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联,所述方法包括采样步骤(520),其用于基于与从图像像素投射的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算所述投影射线上的所述样本定位处的样本值,其中,使用与所述至少一个体素的至少一个分类相关联的强度重新分配函数,根据所述至少一个体素的测得体素值计算所述至少一个体素的所述重定义体素值。
8.一种包括如权利要求1所述的系统(100)的图像采集装置(600)。
9.一种包括如权利要求1所述的系统(100)的工作站(700)。
10.一种要由计算机布置加载的计算机程序产品,包括用于使用射线投射方法对包括多个体素的图像数据集进行可视化的指令,所述多个体素中的每个体素属于至少一个分类,所述至少一个分类中的每个分类与用于根据体素的测得体素值计算所述体素的重定义体素值的强度重新分配函数相关联,所述计算机布置包括处理单元和存储器,在被加载之后,所述计算机程序产品为所述处理单元提供执行以下任务的能力:基于与从图像像素投射的投影射线上的样本定位邻近的至少一个体素的重定义体素值,计算所述投影射线上的所述样本定位处的样本值,其中,使用与所述至少一个体素的至少一个分类相关联的强度重新分配函数,根据所述至少一个体素的测得体素值计算所述至少一个体素的所述重定义体素值。
CN200880121378XA 2007-12-20 2008-12-16 使用多强度重新分配函数进行绘制 Pending CN101903912A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP07123817 2007-12-20
EP07123817.4 2007-12-20
PCT/IB2008/055326 WO2009083861A1 (en) 2007-12-20 2008-12-16 Rendering using multiple intensity redistribution functions

Publications (1)

Publication Number Publication Date
CN101903912A true CN101903912A (zh) 2010-12-01

Family

ID=40451194

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200880121378XA Pending CN101903912A (zh) 2007-12-20 2008-12-16 使用多强度重新分配函数进行绘制

Country Status (4)

Country Link
US (1) US20100265252A1 (zh)
EP (1) EP2223286A1 (zh)
CN (1) CN101903912A (zh)
WO (1) WO2009083861A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103021019A (zh) * 2013-01-10 2013-04-03 广东工业大学 一种基于膝关节ct图像的高逼真度模型体绘制方法

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663803B (zh) * 2012-04-13 2014-11-26 北京工业大学 一种基于RayCasting改进算法的仿真投影DRR生成方法
CN104658028B (zh) * 2013-11-18 2019-01-22 清华大学 在三维图像中快速标记目标物的方法和装置
US10593099B2 (en) * 2017-11-14 2020-03-17 Siemens Healthcare Gmbh Transfer function determination in medical imaging

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2005055148A1 (en) * 2003-11-29 2005-06-16 Vital Images, Inc. Segmented volume rendering using a programmable graphics pipeline
US20050135555A1 (en) * 2003-12-23 2005-06-23 Claus Bernhard Erich H. Method and system for simultaneously viewing rendered volumes

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103021019A (zh) * 2013-01-10 2013-04-03 广东工业大学 一种基于膝关节ct图像的高逼真度模型体绘制方法

Also Published As

Publication number Publication date
WO2009083861A1 (en) 2009-07-09
EP2223286A1 (en) 2010-09-01
US20100265252A1 (en) 2010-10-21

Similar Documents

Publication Publication Date Title
Kurgan et al. Knowledge discovery approach to automated cardiac SPECT diagnosis
US8355553B2 (en) Systems, apparatus and processes for automated medical image segmentation using a statistical model
Shareef et al. Segmentation of medical images using LEGION
US8144949B2 (en) Method for segmentation of lesions
CN101036165B (zh) 用于树模型显像以检测肺栓塞的系统和方法
US8335359B2 (en) Systems, apparatus and processes for automated medical image segmentation
US9427173B2 (en) Determining mechanical force on aneurysms from a fluid dynamic model driven by vessel blood flow information
US8116542B2 (en) Determining hazard of an aneurysm by change determination
US8229188B2 (en) Systems, methods and apparatus automatic segmentation of liver in multiphase contrast-enhanced medical images
US20030099390A1 (en) Lung field segmentation from CT thoracic images
CN101903912A (zh) 使用多强度重新分配函数进行绘制
CN107833229A (zh) 信息处理方法、装置及系统
Tory et al. 4D space-time techniques: a medical imaging case study
CN101802877B (zh) 路径近似绘制
CN101601071B (zh) 应力水平心脏功能分析结果的可视化的系统和方法
Schlachter et al. Principles of visualization in radiation oncology
CN101410870A (zh) 乳腺mri上的自动心脏带探测
Suganthi et al. Interactive Visualization for Understanding and Analyzing Medical Data
Beyer Gpu-based multi-volume rendering of complex data in neuroscience and neurosurgery
Medina et al. Accuracy of connected confidence left ventricle segmentation in 3-D multi-slice computerized tomography images
Nguyen et al. Visualizing veterinary medical data sets with Jetstream
Banerjee et al. Brain Tumor Detection And Analysis From MRI-A Review
Preimb et al. Principles of Visualization in Radiation Oncology
Richardson Image enhancement of cancerous tissue in mammography images
Raspe GPU-assisted diagnosis and visualization of medical volume data

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20101201