CN100498839C - 一种多模态医学体数据三维可视化方法 - Google Patents
一种多模态医学体数据三维可视化方法 Download PDFInfo
- Publication number
- CN100498839C CN100498839C CNB2006100497486A CN200610049748A CN100498839C CN 100498839 C CN100498839 C CN 100498839C CN B2006100497486 A CNB2006100497486 A CN B2006100497486A CN 200610049748 A CN200610049748 A CN 200610049748A CN 100498839 C CN100498839 C CN 100498839C
- Authority
- CN
- China
- Prior art keywords
- registration
- data
- coordinate
- image
- modal
- 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.)
- Expired - Fee Related
Links
Landscapes
- Image Generation (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明涉及一种多模态医学体数据的三维可视化方法。现有技术中的配准多基于二维-二维层面或二维-三维层面,同时不能完整的实现医学体数据配准和融合显示过程。本发明主要包括基于多分辨率的标准化互信息法配准和单点多模态直接体绘制法融合显示。其中基于多分辨率的标准化互信息法配准包括坐标变换、定位判据、多分辨率优化三个步骤;单点多模态直接体绘制法融合显示包括传递函数定义、光照模型的计算、图像合成、多模态数据的单点多模表示、融合显示五个步骤。采用本发明方法可以完整的实现医学体数据配准和融合显示过程,是真正三维-三维意义上的多模态可视化。
Description
技术领域
本发明属于生物医学工程领域,涉及一种多模态医学体数据的三维可视化方法。
背景技术
医学图像数据可分为两大类:一类是X线计算机断层(CT)和核磁共振(MR)等解剖图像,适合于显示人体的解剖细节,具有较高的空间分辨率;另一类是正电子发射断层(PET)、单光子发射断层(SPECT)等功能图像,适合于反映人体组织的功能改变和代谢异常等信息,其空间分辨率差。由于成像原理不同所造成的图像信息局限性,使得单独使用一种模态数据的效果并不理想。通过多模态可视化可以把多个扫描源或不同时间获取的同一病人或多个病人间的数据进行适当综合,将功能信息、新陈代谢及解剖信息等融合在一起为医生提供与解剖位置相关的功能变异,实现解剖图像与功能图像的信息互补,有利于发现新的信息,弥补单模态图像信息不完整、部分信息不准确或不确定引起的缺陷,使临床诊断和治疗、放疗的定位和计划设计、外科手术和疗效评估等更加准确完善。多模态可视化的关键在于多模态数据间的配准和融合显示。
配准是为了多模态体数据间空间上的对齐,主要有基于图像特征匹配和基于图像灰度信息相似性两类方法,目前已申请的配准方法如:刚体变换下基于轮廓多源图像配准方法(200310108731.X)主要针对多源二维图像采用轮廓法实现图像配准,属于基于图像特征的方法;一种基于互信息敏感区域的多模态医学图像配准方法(200310122434.0)提出一种基于互信息敏感区域的图像配准方法,该法在互信息敏感区域的选取上需要经验知识,因此并不是全自动的方法;图像配准方法(200410015656.7)提出基于快速相关法和最大化互信息法来进行图像配准,在该法中坐标变换是基于刚体变换的,即没有考虑弹性缩放问题,而互信法存在重叠区域易变性问题。
目前已申请的有关医学体数据显示方法的发明有:分布式心脏图像多维重建与交互可视化方法(200510020318.7)提出一种采用分布式结构实现旋转扫描二维图像快速重建三维数据场的交互显示的方法;具有封闭轮廓多平面重定格式的体绘制数据的可视化(200510065505.7)提出一种用于可视化与查看矢量相关的体绘制图像的方法;基于点重建的超大规模医学影像三维可视化方法(200310121173.0)提出一种基于点重建的超大规模医学影像三维可视化方法;脑外科手术导航系统中解剖图谱的三维可视化应用方法(03116689.X)提出一种脑外科手术导航系统中解剖图谱的三维可视化应用方法。以上这些发明都是针对单模态医学体数据的显示方法,没有解决多种模态的医学体数据的融合显示问题。容积-容积融合的可视化(03805385.3)提出利用多个容积的方法来融合显示图像的方法,该发明在三维体数据的融合显示采用提取二维层面,其实质是基于二维-二维的融合显示方式,没有解决真正基于三维立体的多模态数据融合显示问题。
以上申请的发明专利表明,在医学多模态体数据可视化领域,配准与融合显示的研究工作都是各自进行的,即专门从事配准研究或对已配准的多模态数据进行可视化,且配准多基于二维-二维层面或二维-三维层面。
发明内容
本发明针对现有技术的不足,从信息论的角度出发,提供一种基于多分辨率的标准化互信息的方法来对三维-三维间的医学体数据进行高效、全自动的弹性配准,并提供一种新的基于单点多模态的直接体绘制法来进行多模态医学体数据的融合显示的方法,从而完整的实现医学体数据配准和融合显示过程,是真正三维-三维意义上的多模态可视化。
本发明主要包括两大步骤:基于多分辨率的标准化互信息法配准和单点多模态直接体绘制法融合显示。
基于多分辨率的标准化互信息法配准包括坐标变换、定位判据、多分辨率优化三个步骤:
(1)坐标变换方法
设要进行多模态可视化的两个体数据分别称为参考体和浮动体,对两个待配准体数据分别定义一个物空间坐标系,其中X轴为沿行扫描方向,Y轴为沿列扫描方向,Z轴则沿从颅顶到颅底的方向。传统的刚体变换包括平移、旋转矩阵,分别以T、R表示,本发明通过增加比例因子S,解决弹性配准问题,假设A是参考体到浮动体的坐标变换函数,则A可表示成:
A=T(tx,ty,tz)*R(φx,φy,φz)*S(sx,sy,sz)
tx、ty、tz分别表示x、y、z轴方向的平移,φx、φy、φz分别表示绕x、y、z轴的旋转角度,sx、sy、sz分别为x、y、z轴方向的缩放比例,利用坐标变换函数A将从浮动体的网格点提取的样本进行坐标变换,如果坐标变换后的坐标未映射到参考体的网格点上,这时通过三次线性插值法得到整数坐标。
(2)定位判据方法
在坐标变换过程中寻找参数:
使之满足下式:
此时,两个体数据对应点间的标准化互信息NIFR最大。通过计算两个体数据重叠部分体素的联合灰度直方图来估计浮动体的熵HF、参考体的熵HR及两者的联合熵HFR,从而计算当前配准位置的标准化互信息NIFR。
标准化互信息NIFR计算公式为:
NIFR=(HF+HR)/HFR (2)
p(f)、p(r)和p(f,r)分别为浮动体和参考体的概率密度函数及它们的联合概率密度函数,nF、nR是参考体和浮动体数据灰度级别的个数,一般都标准值为255。
(3)多分辨率优化方法
由于计算配准判据的时间与采样点数成比例,通过欠采样浮动图像在低分辨率下优化可极大的提高速度。将图像配准过程分解为由粗到细:先将原始图像进行子采样,降低其分辨率,按照(1)和(2)的坐标变换和定位判据的方法在低分辨率下进行配准,得到一个最优的变换数据;然后以此变换参数为初始点,按照(1)和(2)的坐标变换和定位判据的方法进行更高分辨率的配准;由于前次配准的结果已经接近于最终结果,因此,后一次配准可以减小变换参数的搜索范围,从而达到减少迭代次数提高配准速度的效果。另外引入多分辨率搜索方法可以避免陷入局部极值。
配准后的多模态数据采用单点多模态直接体绘制法进行融合显示,包括传递函数定义、光照模型的计算、图像合成、多模态数据的单点多模表示、融合显示五个步骤。
(1)传递函数定义的方法
采集到的医学体数据多为灰度数据,直接体绘制中对每个体素赋予颜色值和不透明值,不同的人体组织器官具有不同的灰度分布,本发明中的传递函数则根据灰度直方图中不同的灰度分布赋予颜色值和不透明度值。
(2)光照模型的计算方法
在直接体绘制中,环境光和周围的点光源都会对最终的显示效果产生影响,在医学体数据的可视化中,光照的逼真问题并不是最重要的因素,为了提高绘制速度,本发明将Phone模型加以简化,采用单个点光源,简化公式为:
I=kaIa+fIl[kd(NL)+ks(NH)] (6)
I为局部光亮度;ka、kd、ks分别为环境、漫反射及镜面反射系数;Ia为环境光亮度分量;I1、f为点光源光强及衰减因子;N为单位法向量;L为点光源的单位入射向量;H为点光源的单位入射向量与单位视线向量的平均值。
采用基于灰度梯度的中心差分法计算中间体素的单位法向量N,对于边界点处的法向量计算则采用向前或向后差分法。中心差分法可表示成:
Grad_x=[f(xi+1,yj,zk)-f(xi-1,yj,zk)]/2Δx
Grad_y=[f(xi,yj+1,zk)-f(xi,yj-1,zk)]/2Δy (7)
Grad_z=[f(xi,yj,zk+1)-f(xi,yj,zk-1)]/2Δz
Grad_x表示x轴方向的灰度梯度,Grad_y表示y轴方向的灰度梯度,Grad_z表示z轴方向的灰度梯度。
(3)图像合成的方法
本发明采用由前向后的图像合成方法,在合成过程中,不透明度β必然逐步增大,当β值趋近于1时,说明该像素点的图像已接近于完全不透明,后面的体元不会再对该像素点的图像有所贡献,因而可以不必再计算了,从而可以省去无效的计算,提高计算速度。对于多个(2)中光照模型的计算结果由前向后进行图像合成,具体公式可表示成:
Coutβout=Cinβin+Cnowβnow(1-βin) (9)
βout=βin+βnow(1-βin) (10)
Cout、βout为经过第i个体元后的颜色值和不透明度值,Cnow、βnow为第i个体元的颜色值和不透明度值,Cin、βin为进入第i个体元的颜色值和不透明度值。
(4)多模态数据的单点多模表示的方法
为了在一个像素点上同时显示多模态体数据信息,我们采用视觉颜色HLS模型表示(3)中图像合成的结果。人眼的颜色知觉主要包括三个要素:即色调(Hue)、饱和度(Saturation)和亮度(Luminance),色调指光的颜色,不同波长的光呈现不同的颜色,具有不同的色调。饱和度指颜色的深浅或浓淡程度,其深浅与颜色中加入白色的比例有关。亮度就是人眼能感觉到的光的明暗程度,光波的能量越大,亮度就越大。
视觉颜色模型以人眼的视觉特征为基础,用色调、饱和度和亮度来表示各个多模态数据的颜色属性。如:HLS模型中的亮度表示参考体属性,而色调和饱和度则表示浮动体的颜色属性。通常以高空间分辨率的模态数据来表示HLS模型中的亮度信息。
(5)融合显示的方法
将HLS模型转化为RGB模型,并加以显示。
采用本发明方法可以完整的实现医学体数据配准和融合显示过程,是真正三维-三维意义上的多模态可视化,同时本发明具有以下特点:
1、用标准化互信息作为多模态医学体数据的配准测度,不需要对两种模态数据强度间关系的性质作任何假设,也不需要对图像作预分割或任何预处理,因此用这种方法来进行多模态体数据的配准,能提高自动化程度和配准精度;同时在判据中考虑了坐标变换中的弹性缩放,适合任何不同模态体数据的配准;且标准化互信息法可解决互信息法难以克服的重叠区域易变性问题。
2、通过多分辨率优化配准过程可极大的提高配准的速度,同时避免陷入局部极值。
3、单点单模态的表示仅仅提供多个单模态可视化结果的数学叠加,表示效果受限,单点多模态直接体绘制法可使每个像素点同时表示多种模态数据属性,能更有效的表现多模态融合信息。
具体实施方式
我们对CT和PET两种模态数据的三维可视化进行实例说明。令CT体数据为参考体,PET体数据为浮动体,各自的坐标为物空间坐标。首先我们设置多分辨率的粗配准为8个像素长度单位,对CT和PET数据进行坐标变换,并计算坐标变换后的标准化互信息,找出标准化互信息最大时的坐标参数,将该最优点作为起始点,缩小坐标变换的范围,同时将坐标变换单位长度减半,即提高分辨率,继续前面的坐标变换及标准化互信息判据过程,直至最优,在新的最优基础上,重复前面的步骤,由粗配准到精细配准步步逼近,直至单位长度为1,得到此时的坐标变换参数即为配准参数。
得到配准参数后,从像空间平面发出射线,确定采样点在射线上的分布间隔,采样点像空间坐标变换至CT物空间坐标,读取CT的相应灰度值,同时利用配准参数得到对应PET物空间中点的灰度值;CT和PET各自对采样点灰度样本分类,赋颜色值和不透明度值,计算光照效应,利用公式(9)(10)由前向后合成图像,把空间分辨率较高的CT颜色值映射到HSL模型中的亮度值,而PET颜色值映射到HSL中的色调和饱和度,最后将HSL模型转换为RGB模型得到融合显示结果。
Claims (1)
1、一种多模态医学体数据三维可视化方法,包括配准和融合显示,其特征在于配准采用基于多分辨率的标准化互信息法,具体包括坐标变换、定位判据、多分辨率优化,其中
(1)坐标变换方法
设要进行多模态可视化的两个体数据分别称为参考体和浮动体,对两个待配准体数据分别定义一个物空间坐标系,其中X轴为沿行扫描方向,Y轴为沿列扫描方向,Z轴则沿从颅顶到颅底的方向,传统的刚体变换包括平移、旋转矩阵,分别以T、R表示,通过增加比例因子S,解决弹性配准问题,假设A是参考体到浮动体的坐标变换函数,则A可表示成:
A=T(tx,ty,tz)*R(φx,φy,φz)*S(sx,sy,sz)
tx、ty、tz分别表示x、y、z轴方向的平移,分别表示绕x、y、z轴的旋转角度,sx、sy、sz分别为x、y、z轴方向的缩放比例,利用坐标变换函数A将从浮动体的网格点提取的样本进行坐标变换,如果坐标变换后的坐标未映射到参考体的网格点上,通过三次线性插值法得到整数坐标;
(2)定位判据方法
在坐标变换过程中寻找参数:
使之满足下式:
此时,两个体数据对应点间的标准化互信息NIFR最大,通过计算两个体数据重叠部分体素的联合灰度直方图来估计浮动体的熵HF、参考体的熵HR及两者的联合熵HFR,从而计算当前配准位置的标准化互信息NIFR,计算公式为:
NIFR=(HF+HR)/HFR
p(f)、p(r)和p(f,r)分别为浮动体和参考体的概率密度函数及它们的联合概率密度函数,nF、nR是参考体和浮动体数据灰度级别的个数,采用标准值255;
(3)多分辨率优化方法
将图像配准过程分解为由粗到细:先将原始图像进行子采样,降低其分辨率,在低分辨率下按照(1)和(2)的坐标变换和定位判据的方法进行配准,得到一个最优的变换数据;然后以此变换参数为初始点,按照(1)和(2)的坐标变换和定位判据的方法进行更高分辨率的配准;在前次配准的基础上减小变换参数的搜索范围;
配准后的多模态数据采用单点多模态直接体绘制法进行融合显示,具体包括传递函数定义、光照模型的计算、图像合成、多模态数据的单点多模表示、融合显示,其中
(1)传递函数定义的方法
将采集到的医学体数据根据人体不同组织和器官的分类对每个体素赋予颜色值和不透明值,其中分类采用灰度直方图的方法;
(2)光照模型的计算方法
对Phone模型加以简化,采用单个点光源,公式为:
I=kaIa+fIl[kd(NL)+ks(NH)]
I为局部光亮度;ka、kd、ks分别为环境、漫反射及镜面反射系数;Ia为环境光亮度分量;I1、f为点光源光强及衰减因子;N为单位法向量;L为点光源的单位入射向量;H为点光源的单位入射向量与单位视线向量的平均值,
采用基于灰度梯度的中心差分法计算中间体素的单位法向量N,表示成:
Grad_x表示x轴方向的灰度梯度,Grad_y表示y轴方向的灰度梯度,Grad_z表示z轴方向的灰度梯度;
(3)图像合成的方法
对于多个(2)中光照模型的计算结果由前向后进行图像合成,具体公式表示成:
Coutβout=Cinβin+Cnowβnow(1-βin)
βout=βin+βnow(1-βin)
Cout、βout为经过第i个体元后的颜色值和不透明度值,Cnow、βnow为第i个体元的颜色值和不透明度值,Cin、βin为进入第i个体元的颜色值和不透明度值;
(4)多模态数据的单点多模表示的方法
采用视觉颜色HLS模型表示(3)中图像合成的结果,该模型根据人眼的颜色知觉三个要素,即色调、饱和度和亮度来表示各个多模态数据的颜色属性;
(5)融合显示的方法
将HLS模型转化为RGB模型,并加以显示。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100497486A CN100498839C (zh) | 2006-03-08 | 2006-03-08 | 一种多模态医学体数据三维可视化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2006100497486A CN100498839C (zh) | 2006-03-08 | 2006-03-08 | 一种多模态医学体数据三维可视化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1818974A CN1818974A (zh) | 2006-08-16 |
CN100498839C true CN100498839C (zh) | 2009-06-10 |
Family
ID=36918973
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100497486A Expired - Fee Related CN100498839C (zh) | 2006-03-08 | 2006-03-08 | 一种多模态医学体数据三维可视化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100498839C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101814193A (zh) * | 2010-03-09 | 2010-08-25 | 哈尔滨工业大学 | 基于gpu加速的实时三维心脏数据体绘制方法 |
Families Citing this family (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101515374B (zh) * | 2008-02-20 | 2010-12-01 | 中国科学院自动化研究所 | 基于图像的个性化真实感虚拟人物造型方法 |
DE102009006147A1 (de) * | 2008-12-23 | 2010-06-24 | Siemens Aktiengesellschaft | Modellgenerator für kardiologische Erkrankungen |
CN101814191B (zh) * | 2009-02-25 | 2011-08-24 | 中国科学院自动化研究所 | 基于二维传递函数的三维图像可视化方法 |
CN101794460A (zh) * | 2010-03-09 | 2010-08-04 | 哈尔滨工业大学 | 基于光线投射体绘制算法的人体心脏三维解剖组织结构模型可视化方法 |
CN101887590B (zh) * | 2010-06-25 | 2011-12-07 | 哈尔滨工业大学 | 一种数字图像可视化组织展示的方法 |
CN102385748B (zh) * | 2010-08-31 | 2013-12-25 | 上海微创医疗器械(集团)有限公司 | 一种图像配准方法 |
US9710730B2 (en) | 2011-02-11 | 2017-07-18 | Microsoft Technology Licensing, Llc | Image registration |
KR20130015146A (ko) * | 2011-08-02 | 2013-02-13 | 삼성전자주식회사 | 의료 영상 처리 방법 및 장치, 영상 유도를 이용한 로봇 수술 시스템 |
CN102945127B (zh) * | 2012-11-05 | 2016-04-27 | 深圳市旭东数字医学影像技术有限公司 | 体绘制显示的交互方法及其系统 |
CN104103083A (zh) * | 2013-04-03 | 2014-10-15 | 株式会社东芝 | 图像处理装置和方法以及医学成像设备 |
CN105144175A (zh) * | 2013-04-24 | 2015-12-09 | 皇家飞利浦有限公司 | 图像可视化 |
CN103268605B (zh) * | 2013-05-15 | 2015-12-23 | 北京理工大学 | 基于高阶统计形变模型的脑图谱与脑图像配准方法 |
CN103295259B (zh) * | 2013-05-31 | 2016-05-18 | 浙江工业大学 | 一种自适应采样的最小梯度夹角预积分光照方法 |
CN103413274A (zh) * | 2013-07-25 | 2013-11-27 | 沈阳东软医疗系统有限公司 | 一种图像补偿方法及装置 |
CN105793897B (zh) * | 2013-12-04 | 2019-12-10 | 皇家飞利浦有限公司 | 图像数据处理 |
CN106296792A (zh) * | 2015-05-12 | 2017-01-04 | 周建龙 | 一种拓扑驱动的体视化光照定义方法 |
CN106296621B (zh) * | 2015-05-22 | 2019-08-23 | 腾讯科技(深圳)有限公司 | 图像处理方法和装置 |
CN104933751B (zh) * | 2015-07-20 | 2017-10-20 | 上海交通大学医学院附属瑞金医院 | 基于局部直方图的心血管冠脉增强的体绘制方法及系统 |
US10943342B2 (en) * | 2016-11-30 | 2021-03-09 | Capsovision Inc. | Method and apparatus for image stitching of images captured using a capsule camera |
CN106652012A (zh) * | 2016-12-09 | 2017-05-10 | 江西中科九峰智慧医疗科技有限公司 | 医学图像体绘制的方法及系统 |
CN107644427B (zh) * | 2017-09-08 | 2021-06-25 | 康达洲际医疗器械有限公司 | 一种基于msew-ca结合灰度直方图的多模态图像分割方法 |
CN107993236A (zh) * | 2017-11-27 | 2018-05-04 | 上海交通大学 | 一种多模态图像处理的方法和平台 |
CN108742678B (zh) * | 2018-06-01 | 2022-02-18 | 妙智科技(深圳)有限公司 | 图像配准方法、装置及计算机可读存储介质 |
US20200103309A1 (en) * | 2018-09-28 | 2020-04-02 | Computational Systems, Inc. | Historical Playback of Waveform Data |
CN110599501B (zh) * | 2019-09-05 | 2022-02-11 | 北京理工大学 | 一种胃肠结构的真实尺度三维重建及可视化方法 |
CN111145336B (zh) * | 2019-12-12 | 2023-05-30 | 东软医疗系统股份有限公司 | 图像绘制方法及装置 |
CN111743564B (zh) * | 2020-07-15 | 2024-05-03 | 北京永新医疗设备有限公司 | 一种临床快速采集引导配准的方法及系统 |
CN113450294A (zh) * | 2021-06-07 | 2021-09-28 | 刘星宇 | 多模态医学图像配准融合方法、装置及电子设备 |
CN113256500B (zh) * | 2021-07-02 | 2021-10-01 | 北京大学第三医院(北京大学第三临床医学院) | 一种多模态图像合成的深度学习神经网络模型系统 |
-
2006
- 2006-03-08 CN CNB2006100497486A patent/CN100498839C/zh not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101814193A (zh) * | 2010-03-09 | 2010-08-25 | 哈尔滨工业大学 | 基于gpu加速的实时三维心脏数据体绘制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1818974A (zh) | 2006-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100498839C (zh) | 一种多模态医学体数据三维可视化方法 | |
Liu et al. | CT synthesis from MRI using multi-cycle GAN for head-and-neck radiation therapy | |
Gao et al. | Pycortex: an interactive surface visualizer for fMRI | |
US7813535B2 (en) | System and method for fused PET-CT visualization for heart unfolding | |
US7439974B2 (en) | System and method for fast 3-dimensional data fusion | |
Tian et al. | Medical image processing and analysis | |
CN103077557B (zh) | 一种自适应分层次胸部大数据显示的实现方法 | |
CN103106685B (zh) | 一种基于gpu的腹部脏器三维可视化方法 | |
CN103295234B (zh) | 基于形变表面模型的医学图像分割系统及方法 | |
CN109191510B (zh) | 一种病理切片的3d重建方法及其装置 | |
AU2012350363A1 (en) | Method and apparatus for the assessment of medical images | |
CN107146262B (zh) | 一种oct图像的三维可视化方法和系统 | |
Philipp et al. | Three-dimensional volume rendering of multidetector-row CT data: applicable for emergency radiology | |
TW201123076A (en) | Three-dimensional display method of medical images | |
CN102509296B (zh) | 基于最大相似性区域合并的胃部ct图像交互式分割方法 | |
Sendra-Balcells et al. | Domain generalization in deep learning for contrast-enhanced imaging | |
CN102426694B (zh) | 一种基于Alpha通道位图技术的图像融合方法 | |
Dai et al. | The application of multi-modality medical image fusion based method to cerebral infarction | |
Levin et al. | Techniques for efficient, real-time, 3D visualization of multi-modality cardiac data using consumer graphics hardware | |
Woo et al. | Super-resolution reconstruction for tongue MR images | |
Adeshina et al. | Locating abnormalities in brain blood vessels using parallel computing architecture | |
Teßmann et al. | GPU accelerated normalized mutual information and B-Spline transformation. | |
CN103778658A (zh) | 一种快速展示体数据特征的可视化方法 | |
Lei | A brain nuclear magnetic resonance image segmentation algorithm based on non-rigid registration | |
Xu et al. | Seamless warping of diffusion tensor fields |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20090610 Termination date: 20120308 |