CN101373541A - 医学图像体绘制方法及装置 - Google Patents

医学图像体绘制方法及装置 Download PDF

Info

Publication number
CN101373541A
CN101373541A CNA2008101705432A CN200810170543A CN101373541A CN 101373541 A CN101373541 A CN 101373541A CN A2008101705432 A CNA2008101705432 A CN A2008101705432A CN 200810170543 A CN200810170543 A CN 200810170543A CN 101373541 A CN101373541 A CN 101373541A
Authority
CN
China
Prior art keywords
volume data
texture
plane
image
behalf
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.)
Granted
Application number
CNA2008101705432A
Other languages
English (en)
Other versions
CN101373541B (zh
Inventor
赵大哲
杨金柱
栗伟
冯朝路
耿欢
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.)
Neusoft Corp
Original Assignee
Neusoft Corp
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 Neusoft Corp filed Critical Neusoft Corp
Priority to CN2008101705432A priority Critical patent/CN101373541B/zh
Publication of CN101373541A publication Critical patent/CN101373541A/zh
Application granted granted Critical
Publication of CN101373541B publication Critical patent/CN101373541B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Image Generation (AREA)

Abstract

本发明公开了一种医学图像体绘制方法及装置,所述方法包括:将体数据载入显存,并生成对应不同组织的体数据纹理;建立与各组织对应的传输函数;利用所述传输函数和体数据纹理分别绘制对应各组织的图像。利用本发明,可以提高绘制速度与绘制质量。

Description

医学图像体绘制方法及装置
技术领域
本发明涉及体绘制技术,具体涉及一种医学图像体绘制方法及装置。
背景技术
体绘制又称直接体绘制,是近年来发展迅速的一种三维数据场可视化方法。体绘制不生成中间图元,直接由三维数据场产生屏幕的二维图像。首先通过对离散的三维数据场的重构得到初始的连续三维数据场,然后在该连续的三维数据场中进行重采样,得到采样点的函数值。对于新的采样点,按照设定的阻光度和颜色模型进行分类和颜色赋值,最后根据体绘制的光学模型,通过图像合成得到绘制结果图像。
体绘制方法从提出以后,经过研究者20多年的完善与改进,从理论、方法及实现上都已经形成比较完善的体系,但是仍然存在一些实际的问题,主要有以下几点:
(1)对于大规模的体数据(如5123),在普通的PC机上实现体绘制算法,虽然经过优化,其绘制速度仍然较慢,不能满足实时绘制的要求;
(2)没有较好的数据分类方法,传输函数的设计直接影响绘制效果;
(3)无法解决对多于三个以上的组织标定。
发明内容
本发明提供一种医学图像体绘制方法及装置,以提高绘制速度与绘制质量。
本发明提供一种医学图像体绘制方法,包括:
将体数据载入显存,并生成对应不同组织的体数据纹理;
建立与各组织对应的传输函数;
利用所述传输函数和体数据纹理分别绘制对应各组织的图像。
可选地,在将体数据载入显存之前,对所述体数据进行预处理。
所述预处理的过程包括:对所述体数据进行增强、和/或滤波、和/或插值处理。
所述对体数据进行插值处理具体为:利用三线性插值算法对所述体数据进行插值处理。
优选地,所述方法还包括:将所述传输函数通过采样转化为二维映射表,并将所述二维映射表作为纹理载入显存,所述纹理的横轴表示数据的标量值,纵轴表示对应的标注号;
所述利用所述传输函数和体数据纹理分别绘制对应各组织的图像包括:
利用对应不同组织的体数据纹理查找所述二维映射表,获得所述体数据的光学属性;
根据所述体数据的光学属性分别绘制对应各组织的图像。
优选地,所述方法还包括:
利用CPU生成代理几何体;
根据所述代理几何体产生纹理坐标;
根据所述纹理坐标将对应各组织的图像显示到屏幕上。
优选地,所述利用CPU生成代理几何体包括:
利用观察模型矩阵将垂直于视线方向的平面切割形成的体数据包围盒顶点形成观察坐标;
通过计算垂直于视线方向的平面的平面方程得到一组互相平行的平面;
按照从前到后的顺序,计算每个平面与所述体数据包围盒的各条棱的交点;
根据所述观察坐标以及得到的交点,将每个平面形成的多边形转化为三角形网格。
本发明提供一种医学图像体绘制装置,包括:
体数据纹理生成单元,用于将体数据载入显存,并生成对应不同组织的体数据纹理;
传输函数建立单元,用于建立与各组织对应的传输函数;
绘制单元,用于利用所述传输函数和体数据纹理分别绘制对应各组织的图像。
可选地,还包括:
预处理单元,用于在将体数据载入显存之前,对所述体数据进行预处理。
可选地,所述预处理单元的预处理的过程包括:对所述体数据进行增强、和/或滤波、和/或插值处理。
优选地,还包括:
映射单元,用于将所述传输函数通过采样转化为二维映射表,并将所述二维映射表作为纹理载入显存,所述纹理的横轴表示数据的标量值,纵轴表示对应的标注号;
所述绘制单元包括:
光学属性获取子单元,用于利用对应不同组织的体数据纹理查找所述二维映射表,获得所述体数据的光学属性;
绘制子单元,用于根据所述体数据的光学属性分别绘制对应各组织的图像。
优选地,还包括:
代理几何体生成单元,用于利用CPU生成代理几何体;
纹理坐标生成单元,用于根据所述代理几何体产生纹理坐标;
显示单元,用于根据所述纹理坐标将对应各组织的图像显示到屏幕上。
优选地,所述代理几何体生成单元包括:
观察坐标生成子单元,用于利用观察模型矩阵将垂直于视线方向的平面切割形成的体数据包围盒顶点形成观察坐标;
平面获取子单元,用于通过计算垂直于视线方向的平面的平面方程得到一组互相平行的平面;
交点计算子单元,用于按照从前到后的顺序,计算每个平面与所述体数据包围盒的各条棱的交点;
转化子单元,用于根据所述观察坐标以及得到的交点,将每个平面形成的多边形转化为三角形网格。
本发明医学图像体绘制方法及装置,基于计算机图形处理器(GraphicProcessing Unit,GPU)进行体绘制,将体数据载入显存,并生成对应不同组织的体数据纹理,建立与各组织对应的传输函数,利用所述传输函数和体数据纹理分别绘制对应各组织的图像,对不同规模的数据集,均能有效地提高绘制速度与绘制质量。
进一步地,将对应于各组织的传输函数转化为二维映射表,并将所述二维映射表作为纹理载入显存,解决了医学影像多组织的标定问题,便于同时观察这些组织的空间位置关系以及形状结构特点。
进一步地,利用CPU生成代理几何体的方法产生纹理坐标,减小了GPU运算资源的耗费,使CPU与GPU的资源得到了合理的利用。
附图说明
图1是本发明实施例医学图像体绘制方法的实现流程图;
图2是本发明实施例中与传输函数对应的二维映射表的结构示意图;
图3是本发明实施例医学图像体绘制方法中基于CPU生成代理几何体的实现流程图;
图4是本发明实施例医学图像体绘制装置的一种结构示意图;
图5是本发明实施例医学图像体绘制装置的另一种结构示意图;
图6是本发明实施例医学图像体绘制装置的另一种结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明实施例的方案,下面结合附图和实施方式对本发明实施例作进一步的详细说明。
GPU是一种图像处理的处理器,它使图像处理计算减少了对CPU的依赖,在图像处理时它的工作效率远高于CPU,并代替部分原本属于CPU的工作。因此,本发明实施例医学图像体绘制方法,基于GPU进行体绘制,将体数据载入显存,并生成对应不同组织的体数据纹理,建立与各组织对应的传输函数,利用所述传输函数和体数据纹理分别绘制对应各组织的图像,对不同规模的数据集,均能有效地提高绘制速度与绘制质量。
如图1所示,是本发明实施例医学图像体绘制方法的一种实现流程图,主要包括以下步骤:
步骤101,将体数据载入显存。
在将体数据载入显存时,需要首先将这些体数据载入CPU,然后再载入显存;
在某些情况下,比如使用直接采集的体数据成像质量不佳时,在将所述体数据载入CPU后,还需要对载入显存的体数据进行预处理,所述预处理的过程包括:对所述体数据进行增强、和/或滤波处理等一些操作,然后再将处理后的体数据载入显存。所述增强处理是指增强图像中的有用信息,它可以是一个失真的过程,其目的是要增强视觉效果,将原来不清晰的图像变得清晰或强调某些感兴趣的特征,抑制不感兴趣的特征,以改善图像质量、丰富信息量、加强图像判读和识别效果。所述滤波处理是指只让一定频率信号通过,可去掉图像中的噪声,增强边缘信号,使模糊的图像变得清晰。
由于体数据通常会有不同的大小和类型,载入的数据多数都不能满足各向比例相同的特性,为了完成体绘制,还可能需要对这些数据进行插值处理,还原体数据原有的比例。所述插值处理是指应用几个相邻已知数据点,通过三线性算法计算出一个插值点,这样可以增加图像分辨率,提高成像质量。
具体地,可以使用三线性插值算法对所述体数据进行插值处理。该处理可以在CPU上进行,也就是说,在将体数据载入显存前进行,也可以在GPU上进行,即将体数据载入显存后进行,这样,可以提高处理效率。
除此之外,在将体数据生成纹理载入内存之前,有时还需要选择计算梯度、减少采样数据等操作。不同的计算梯度的方法对成像质量和运算时间会有影响,针对不同的数据源(不同人体、不同组织)可以选择不同的计算梯度方法。
步骤102,将载入显存的体数据生成对应不同组织的体数据纹理。
载入显存的数据以纹理形式存在,通常,纹理包括一维纹理、二维纹理和三维纹理。对于体数据,包括三个维度,以体现数据的三维空间性质。
在本发明实施例中,所述体数据可能对应了不同的组织,因此,可以将这些体数据按照不同的组织进行分类,生成对应不同组织的体数据纹理。
步骤103,建立与各组织对应的传输函数。
为了能够进一步提高绘制的医学图像的质量和效率,并能够明显地从绘制的图像上区分不同的组织,本发明实施例中,针对不同组织各自的特点,对每个组织建立一个与其对应的传输函数。
步骤104,利用所述传输函数和体数据纹理分别绘制对应各组织的图像。
本发明医学图像体绘制方法,对不同规模的数据集,均能有效地提高绘制速度与绘制质量。
如果在绘制时,对每个采样点(即体数据)都要根据传输函数计算映射的光学属性,无疑会耗费大量的计算时间,因此,为了进一步提高绘制效率,还可以预先将各传输函数转化为一个二维映射表,并将所述二维映射表作为纹理载入显存,所述纹理的横轴表示数据的标量值,纵轴表示对应的标注号。
这样,在进行绘制时,就可以利用对应不同组织的体数据纹理查找相应的二维映射表,获得所述体数据的光学属性,比如颜色、不透明度等,然后,根据所述体数据的光学属性分别绘对应各组织的图像。
查表的过程,可以在GPU上完成。将所述二维映射表作为纹理载入显存中,通过在片元着色器中查找纹理的方法即可取得需要的映射值,获得所述体数据的光学属性。
在本发明实施例中,使用一维传输函数,一维传输函数本身可以只使用一维纹理,使用数据标量值作为自变量,对应结果是RGBA(red,green,blue,alpha)的光学属性。所述RGBA代表像素点三基色的成份构成,其中alpha用来归一化,还可以做不透明度来应用。
在设计传输函数时,考虑了体数据的标注问题,对于医学可视化系统,通常需要对检测出的不同组织、器官加以标注,以不同的颜色或者绘制效果显示出来,这就需要满足对于同一个数据值,其对应的光学属性可能不相同。
在本发明实施例中,将标注与传输函数相结合,如图2所示,将传输函数转化为一个二维映射表,并将所述二维映射表作为纹理载入显存,纹理的x轴表示数据值,y轴表示对应的标注值。比如,本实施例设计可以标注8个不同的组织,每个组织对应一个传输函数。将这个二维表生成二维纹理载入显存,颜色和透明度都需要转化到0~255之间。在对图像进行渲染时,在片元着色器中完成查找,将采样点的数据值和标注值作为纹理坐标,查询传输函数二维纹理,即可获得采样点的光学属性。
可见,本发明实施例的方法将对应于各组织的传输函数转化为二维映射表,并将所述二维映射表作为纹理载入显存。从而可以针对不同的组织、器官,使用户可以自行选择应用标定的颜色、阻光度等,解决了医学影像多组织的标定问题,可以更好地观察不同组织的空间位置,结构等信息。
在将对应各组织的图像显示到屏幕上时,需要代理几何体生成纹理坐标,根据所述纹理坐标将纹理切片绘制到相应的位置,即将对应各组织的图像显示到屏幕上,恢复体数据表示的原有结构。具体可以采用现有技术中通过图形API(Application Program Interface,应用程序编程接口)自动生成纹理坐标的方法。但这种方法需要耗费GPU的运算资源。
在绘制阶段,在按照一定顺序画出所述纹理切片之前,需要适当的设置渲染状态,包括禁用光照和剔除、建立透明混合状态等。为了在不透明的几何体中进行绘制,还需要激活深度检测(即绘制过程中程序会以不同数据点的深度值来决定遮挡问题),禁用写入深度缓冲的功能。因为要绘制的物体都有0-1的透明度,为防止透明物体的深度值比不透明物体的小,而发生深度缓冲被重写,不透明物体信息丢失的现象,需要将深度缓冲设置为只读模式,这样可以保证它的值不被修改。另外,还需要根据用户配置绑定不同的着色程序,按照一定的顺序绘制代理几何体。光栅化过程中执行片元着色程序,在GPU上完成颜色映射、图像合成。如果开启光照,还需要实时地计算光照效果,输出最终的像素点颜色。
上述绘制阶段的这些操作与现有技术中的实现过程相类似,在此不再一一详细描述。
为了进一步降低对GPU资源的消耗,,使CPU与GPU的资源得到了合理的利用。在本发明实施例中,可以在CPU上生成代理几何体。
所述代理几何体由多边形组构成,而多边形则是垂直于视线方向把体切成薄片形成的。
对于任意的视线方向,需要生成与之垂直的一系列纹理多边形。具体地,首先求出一系列垂直于视线方向的平面切割形成的体数据包围盒所形成的截面,这些截面是一些多边形,即切片多边形,这些多边形组就是代理几何体。
如图3所示,是本发明实施例医学图像体绘制方法中基于CPU生成代理几何体的实现流程图,主要包括以下步骤:
步骤301,利用观察模型矩阵将垂直于视线方向的平面切割形成的体数据包围盒顶点形成观察坐标。
所述观察模型矩阵是指将物体由物体坐标系转换到世界坐标系,并决定它的位置及朝向的矩阵。
步骤302,通过计算垂直于视线方向的平面的平面方程得到一组互相平行的平面。
假定视点为vp(x,y,z),视线的负方向为vec(x,y,z),垂直于视线方向的平面的平面方程为Ax+By+Cz+D=0,则有
A = vec . x B = vec . y C = vec . z
取D的初始值为零,即初始切片平面通过体数据的原点。通过调整D的值,即可得到一组互相平行的平面。
步骤303,计算每个平面与所述体数据包围盒的各条棱的交点。
在计算每个平面与所述体数据包围盒的各条棱的交点时,可以按照一定的平面顺序来计算,比如从前到后的顺序、或者从后到前的顺序。
体数据与垂直于视线方向的平面的切片形成的多边形,可以通过所述平面与体数据包围盒的各条棱的交点来确定。
设置数组pts[8]存储长方体的各顶点坐标,整型的二维数组edges[12][2]存储每条棱的端点在pts[8]中的索引。
对于每条棱,判断其与已知平面是否相交,用函数表示已知平面的平面方程为:F(x,y,z)=Ax+By+Cz+D。
根据点与平面的关系,若点pts在平面上方,则F(pts)>0;若点在平面下方,则F(pts)<0;若点在平面上,则F(pts)=0。
设棱的两个端点分别为pt1、pt2,判断该棱与平面是否有交点,有以下几种情况:
若F(pt1)·F(pt2)>0,则两点在平面同一侧,棱与平面无交点;
若F(pt1)·F(pt2)<0,则两点在平面两侧,棱与平面相交;
若F(pt1)=0,F(pt2)=0,则两点都在平面上,平面经过该棱;
若F(pt1)=0,F(pt2)≠0,则点pt1在平面上;
若F(pt1)≠0,F(pt2)=0,则点pt2在平面上。
若棱与平面相交,需要求出它们的交点。
t = F ( pt 2 ) F ( pt 2 ) - F ( pt 1 ) , 则交点edgeInt的坐标计算如下:
edgeInt . x = t &times; pt 1 . x + ( 1 - t ) &times; pt 2 . x edgeInt . y = t &times; pt 1 . y + ( 1 - t ) &times; pt 2 . y edgeInt . z = t &times; pt 1 . z + ( 1 - t ) &times; pt 2 . z
步骤304,根据所述观察坐标以及得到的交点,将每个平面形成的多边形转化为三角形网格。
为了提高渲染效率,需要将非三角形的切面转化为由三角形组成的切面。在将所述多边形转化为三角网格时,可以对一个面上的相交点(即所述平面与体数据包围盒的各条棱的交点)求平均值,利用该平均值计算代理多边形(也就是前面提到的每个平面形成的多边形)的中心点,将相交点投影到X-Y平面上,使用该代理多边形的第一个顶点或者X轴作参照,计算他们围绕中心点的角度,对该代理多边形顶点按顺时针或者逆时针排序。从而可以将每个平面形成的代理多边形转化为以所述中心点为公共点的三角形扇面或者三角形带。
在实际应用中,可以在系统初始化时将体数据及一些参数设置载入显存,生成三维纹理,并进行体绘制。当用户输入的体数据发生变化时,例如观察或者渲染参数改变时,则需要更新三维纹理,重新进行体绘制。
需要说明的是,在初始化和用户输入发生变化时,都需要重新计算代理几何体,然后将所述代理几何体存储在顶点数组(即存储三维空间点坐标信息的数组)之中。对于用户输入改变,如视角变换、传输函数变更和渲染模式改变,还有具体环境下的配置参数,都需要重新计算代理几何体。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,所述的程序可以存储于一计算机可读取存储介质中,所述的存储介质,如:ROM/RAM、磁碟、光盘等。
相应于本发明实施例医学图像体绘制方法,本发明实施例还提供了一种医学图像体绘制装置。
如图4所示,是本发明实施例医学图像体绘制装置的一种结构示意图:
该装置包括:体数据纹理生成单元401、传输函数建立单元402、绘制单元403,还可进一步包括:预处理单元404。其中:
体数据纹理生成单元401,用于将体数据载入显存,并生成对应不同组织的体数据纹理;传输函数建立单元402,用于建立与各组织对应的传输函数;绘制单元403,用于利用所述传输函数和体数据纹理分别绘制对应各组织的图像;预处理单元404,用于在将体数据载入显存之前,对所述体数据进行预处理,所述预处理单元404的预处理的过程包括:对所述体数据进行增强、和/或滤波、和/或插值处理,所述插值处理可以使用三线性插值算法对所述体数据进行插值处理。该处理可以在CPU上进行,也就是说,在将体数据载入显存前进行,也可以在GPU上进行,即将体数据载入显存后进行,这样,可以提高处理效率。
本发明实施例医学图像体绘制装置可以适用于不同规模的数据集,均能有效地提高绘制速度与绘制质量。
参照图5,是本发明实施例医学图像体绘制装置的另一种结构示意图:
与图4所示实施例的区别在于,在该实施例中,除了包括上述各功能单元之外,还包括:映射单元501,用于将所述传输函数通过采样转化为二维映射表,并将所述二维映射表作为纹理载入显存,所述纹理的横轴表示数据的标量值,纵轴表示对应的标注号。
这样,在该实施例中,绘制单元503可以进一步包括:光学属性获取子单元531和绘制子单元532。其中:
光学属性获取子单元531,用于利用对应不同组织的体数据纹理查找所述二维映射表,获得所述体数据的光学属性;绘制子单元532,用于根据所述体数据的光学属性分别绘对应各组织的图像。
所述传输函数通过采样转化为二维映射表的过程,以及绘制过程的具体实现可参照前面本发明实施例医学图像体绘制方法中的描述,在此不再赘述。
利用本发明实施例医学图像体绘制装置,不仅可以有效地提高图像绘制速度与绘制质量,而且解决了医学影像多组织的标定问题使用户可以针对不同的组织,自行选择应用标定的颜色、阻光度,可以更好地观察各组织的空间位置、结构等信息。
参照图6,是本发明实施例医学图像体绘制装置的另一种结构示意图:
与图5所示实施例的区别在于,在该实施例中,除了包括上述各功能单元之外,还包括:代理几何体生成单元601、纹理坐标生成单元602和显示单元603。其中:
代理几何体生成单元601,用于利用CPU生成代理几何体;纹理坐标生成单元602,用于根据所述代理几何体产生纹理坐标;显示单元603,用于根据所述纹理坐标将对应各组织的图像显示到屏幕上。
代理几何体生成单元601包括:
观察坐标生成子单元611,用于利用观察模型矩阵将垂直于视线方向的平面切割形成的体数据包围盒顶点形成观察坐标;
平面获取子单元612,用于通过计算垂直于视线方向的平面的平面方程得到一组互相平行的平面;
交点计算子单元613,用于按照从前到后的顺序,计算每个平面与所述体数据包围盒的各条棱的交点;
转化子单元614,用于根据所述观察坐标以及得到的交点,将每个平面形成的多边形转化为三角形网格。
利用CPU生成代理几何体的详细过程可参照对图3所示流程的描述,在此不再赘述。
本发明实施例医学图像体绘制装置,不仅可以有效地提高图像绘制速度与绘制质量,而且利用CPU生成代理几何体的方法产生纹理坐标,减小了GPU运算资源的耗费,使CPU与GPU的资源得到了合理的利用。
本发明实施例进行了详细介绍,本文中应用了具体实施方式对本发明进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及设备;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (13)

1.一种医学图像体绘制方法,其特征在于,包括:
将体数据载入显存,并生成对应不同组织的体数据纹理;
建立与各组织对应的传输函数;
利用所述传输函数和体数据纹理分别绘制对应各组织的图像。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
在将体数据载入显存之前,对所述体数据进行预处理。
3.根据权利要求2所述的方法,其特征在于,所述预处理的过程包括:
对所述体数据进行增强、和/或滤波、和/或插值处理。
4.根据权利要求3所述的方法,其特征在于,所述对体数据进行插值处理具体为:
利用三线性插值算法对所述体数据进行插值处理。
5.根据权利要求1所述的方法,其特征在于,所述方法还包括:将所述传输函数通过采样转化为二维映射表,并将所述二维映射表作为纹理载入显存,所述纹理的横轴表示数据的标量值,纵轴表示对应的标注号;
所述利用所述传输函数和体数据纹理分别绘制对应各组织的图像包括:
利用对应不同组织的体数据纹理查找所述二维映射表,获得所述体数据的光学属性;
根据所述体数据的光学属性分别绘制对应各组织的图像。
6.根据权利要求5所述的方法,其特征在于,还包括:
利用CPU生成代理几何体;
根据所述代理几何体产生纹理坐标;
根据所述纹理坐标将对应各组织的图像显示到屏幕上。
7.根据权利要求6所述的方法,其特征在于,所述利用CPU生成代理几何体包括:
利用观察模型矩阵将垂直于视线方向的平面切割形成的体数据包围盒顶点形成观察坐标;
通过计算垂直于视线方向的平面的平面方程得到一组互相平行的平面;
按照从前到后的顺序,计算每个平面与所述体数据包围盒的各条棱的交点;
根据所述观察坐标以及得到的交点,将每个平面形成的多边形转化为三角形网格。
8.一种医学图像体绘制装置,其特征在于,包括:
体数据纹理生成单元,用于将体数据载入显存,并生成对应不同组织的体数据纹理;
传输函数建立单元,用于建立与各组织对应的传输函数;
绘制单元,用于利用所述传输函数和体数据纹理分别绘制对应各组织的图像。
9.根据权利要求8所述的装置,其特征在于,还包括:
预处理单元,用于在将体数据载入显存之前,对所述体数据进行预处理。
10.根据权利要求9所述的装置,其特征在于,所述预处理单元的预处理的过程包括:对所述体数据进行增强、和/或滤波、和/或插值处理。
11.根据权利要求8所述的装置,其特征在于,还包括:
映射单元,用于将所述传输函数通过采样转化为二维映射表,并将所述二维映射表作为纹理载入显存,所述纹理的横轴表示数据的标量值,纵轴表示对应的标注号;
所述绘制单元包括:
光学属性获取子单元,用于利用对应不同组织的体数据纹理查找所述二维映射表,获得所述体数据的光学属性;
绘制子单元,用于根据所述体数据的光学属性分别绘制对应各组织的图像。
12.根据权利要求11所述的装置,其特征在于,还包括:
代理几何体生成单元,用于利用CPU生成代理几何体;
纹理坐标生成单元,用于根据所述代理几何体产生纹理坐标;
显示单元,用于根据所述纹理坐标将对应各组织的图像显示到屏幕上。
13.根据权利要求12所述的装置,其特征在于,所述代理几何体生成单元包括:
观察坐标生成子单元,用于利用观察模型矩阵将垂直于视线方向的平面切割形成的体数据包围盒顶点形成观察坐标;
平面获取子单元,用于通过计算垂直于视线方向的平面的平面方程得到一组互相平行的平面;
交点计算子单元,用于按照从前到后的顺序,计算每个平面与所述体数据包围盒的各条棱的交点;
转化子单元,用于根据所述观察坐标以及得到的交点,将每个平面形成的多边形转化为三角形网格。
CN2008101705432A 2008-10-17 2008-10-17 医学图像体绘制方法及装置 Expired - Fee Related CN101373541B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2008101705432A CN101373541B (zh) 2008-10-17 2008-10-17 医学图像体绘制方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008101705432A CN101373541B (zh) 2008-10-17 2008-10-17 医学图像体绘制方法及装置

Publications (2)

Publication Number Publication Date
CN101373541A true CN101373541A (zh) 2009-02-25
CN101373541B CN101373541B (zh) 2010-09-15

Family

ID=40447692

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008101705432A Expired - Fee Related CN101373541B (zh) 2008-10-17 2008-10-17 医学图像体绘制方法及装置

Country Status (1)

Country Link
CN (1) CN101373541B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103559733A (zh) * 2013-10-09 2014-02-05 浙江大学 支持三维体数据内部视点漫游的球面体绘制方法
CN103745496A (zh) * 2013-12-11 2014-04-23 南方医科大学 基于二维图像为交互界面的转换函数的直接体绘制方法
CN105719333A (zh) * 2014-12-02 2016-06-29 上海联影医疗科技有限公司 三维图像数据处理方法及装置
CN106530386A (zh) * 2016-11-25 2017-03-22 上海联影医疗科技有限公司 医学图像的体绘制方法及其系统
CN106934841A (zh) * 2017-03-07 2017-07-07 上海联影医疗科技有限公司 生成彩色医学影像的方法及医学系统
CN107280703A (zh) * 2016-07-22 2017-10-24 珠海医凯电子科技有限公司 基于gpu平台的实时3d超声扫描变换方法
US10748280B2 (en) 2015-12-31 2020-08-18 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for image processing
US11094067B2 (en) 2014-12-02 2021-08-17 Shanghai United Imaging Healthcare Co., Ltd. Method and system for image processing
CN113436305A (zh) * 2021-06-23 2021-09-24 北京理工大学 一种由直接体绘制引导的高质量几何建模方法
CN114255249A (zh) * 2021-12-22 2022-03-29 沈阳东软智能医疗科技研究院有限公司 一种医疗影像纹理特征的处理方法、装置及设备

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103559733A (zh) * 2013-10-09 2014-02-05 浙江大学 支持三维体数据内部视点漫游的球面体绘制方法
CN103559733B (zh) * 2013-10-09 2017-02-08 浙江大学 支持三维体数据内部视点漫游的球面体绘制方法
CN103745496A (zh) * 2013-12-11 2014-04-23 南方医科大学 基于二维图像为交互界面的转换函数的直接体绘制方法
CN103745496B (zh) * 2013-12-11 2017-01-18 南方医科大学 基于二维图像为交互界面的转换函数的直接体绘制方法
CN105719333A (zh) * 2014-12-02 2016-06-29 上海联影医疗科技有限公司 三维图像数据处理方法及装置
US11094067B2 (en) 2014-12-02 2021-08-17 Shanghai United Imaging Healthcare Co., Ltd. Method and system for image processing
US10748280B2 (en) 2015-12-31 2020-08-18 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for image processing
US11769249B2 (en) 2015-12-31 2023-09-26 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for image processing
CN107280703A (zh) * 2016-07-22 2017-10-24 珠海医凯电子科技有限公司 基于gpu平台的实时3d超声扫描变换方法
CN106530386B (zh) * 2016-11-25 2020-08-04 上海联影医疗科技有限公司 医学图像的体绘制方法及其系统
CN106530386A (zh) * 2016-11-25 2017-03-22 上海联影医疗科技有限公司 医学图像的体绘制方法及其系统
CN106934841A (zh) * 2017-03-07 2017-07-07 上海联影医疗科技有限公司 生成彩色医学影像的方法及医学系统
CN113436305A (zh) * 2021-06-23 2021-09-24 北京理工大学 一种由直接体绘制引导的高质量几何建模方法
CN114255249A (zh) * 2021-12-22 2022-03-29 沈阳东软智能医疗科技研究院有限公司 一种医疗影像纹理特征的处理方法、装置及设备

Also Published As

Publication number Publication date
CN101373541B (zh) 2010-09-15

Similar Documents

Publication Publication Date Title
CN101373541B (zh) 医学图像体绘制方法及装置
CN112270756B (zh) 一种应用于bim模型文件的数据渲染方法
Naylor Interactive solid geometry via partitioning trees
Heckbert Graphics gems
Weiskopf et al. Volume clipping via per-fragment operations in texture-based volume visualization
US5861891A (en) Method, system, and computer program for visually approximating scattered data
Lu et al. Illustrative interactive stipple rendering
CN103500463B (zh) 一种gpu上多层形状特征融合的可视化方法
Grimm et al. Flexible direct multi-volume rendering in dynamic scenes
Livnat et al. Interactive point-based isosurface extraction
CN105549929A (zh) 用于改善图形性能的方法、装置和计算机程序产品
CN101763649A (zh) 一种增强模型轮廓的表面点绘制方法
US20240212863A1 (en) Unified anisotropic volume and surface mesh storage
Kämpe et al. Fast, memory-efficient construction of voxelized shadows
Suarez et al. Real-time 3D rendering with hatching
CN1776747A (zh) 医学图像中基于gpu硬件加速的体绘制方法
Krishnamurthy et al. Optimized GPU evaluation of arbitrary degree NURBS curves and surfaces
Nah et al. MobiRT: an implementation of OpenGL ES-based CPU-GPU hybrid ray tracer for mobile devices
Andersen et al. Hybrid fur rendering: combining volumetric fur with explicit hair strands
Crawfis Real-time slicing of data space
US20030160776A1 (en) Geometric folding for cone-tree data compression
CN105301652A (zh) 一种三维地震数据固定轴向切片动态判断体绘制方法
Ahrens et al. A modular extensible visualization system architecture for culled prioritized data streaming
Adams et al. Boolean operations on surfel-bounded solids using programmable graphics hardware
Lu et al. Multivariate visualization for atmospheric pollution

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100915

Termination date: 20201017

CF01 Termination of patent right due to non-payment of annual fee