CN102411780B - 一种基于配准的ct图像全心脏自动分割系统 - Google Patents
一种基于配准的ct图像全心脏自动分割系统 Download PDFInfo
- Publication number
- CN102411780B CN102411780B CN201110264063.4A CN201110264063A CN102411780B CN 102411780 B CN102411780 B CN 102411780B CN 201110264063 A CN201110264063 A CN 201110264063A CN 102411780 B CN102411780 B CN 102411780B
- Authority
- CN
- China
- Prior art keywords
- image
- registration
- module
- transformation
- parameter
- 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
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种基于配准的CT图像全心脏自动分割系统,该系统利用不同个体间心脏CT的相似性,把配准应用于分割,实现全心脏的自动分割。该系统包括如下模块:输入模块、分割模块、粗配准模块、精配准模块、变换模块和输出模块。该系统能实现全心脏的自动分割并能同时分割出多个心脏腔室。
Description
技术领域
本发明涉及医学图像分割系统,具体涉及一种基于配准的CT图像全心脏自动分割系统。
背景技术
医学图像分割就是对感兴趣区域进行提取的一个过程,是进行医学图像三维可视化前的关键步骤。心脏分割能把心脏各部位分割出来,从而为临床诊断提供大量的解剖信息,并能为心脏手术的术前计划以及心脏介入手术提供依据。对分割后的各腔室进行体积计算还能为心功能评价提供量化的信息。
目前,医学图像的分割可分为手动分割、半自动分割和自动分割。手动分割就是医生利用临床知识在医学图像上勾画出感兴趣区域的轮廓,精度最高。但效率较低,分割结果不可重复,而且分割结果受到分割者的经验制约。半自动分割是利用计算机通过人机交互进行分割,分割结果在一定程度仍受到分割者经验知识的制约。自动分割是由计算机完成医学图像分割的全过程,其分割结果能够重现。传统的自动分割完全脱离医生的临床经验,因此应用于临床时有一定的局限性。而本发明把医学专家的临床经验知识与自动分割过程相结合,能克服传统自动分割的局限性,因此对临床诊断有重要的意义。
发明内容
本发明利用不同个体间心脏CT的相似性,把配准应用于分割中,把医学专家的临床经验知识与自动分割相结合,实现了全心脏及心脏各腔室的自动分割。
本发明的全心脏自动分割系统通过以下技术方案实现:
本发明提出一种基于配准的CT图像全心脏自动分割系统,包括:
输入模块100,用于输入心脏CT序列图像模板和待分割的心脏CT图像,并将心脏CT模板图像传送至分割模块200,将输入心脏CT序列图像模板和待分割的心脏CT图像传送至粗配准模块300;
分割模块200,用于在CT序列图像模板上分割出心脏的各个部位,为各部位标记上标记符,制作出Label图像,并传送至变换模块500;该Label图像能应用于所有分割实例中,无需重复制作;
粗配准模块300,用于基于仿射变换的粗配准,对模板作出平移、缩放、旋转变换,并将粗配准结果传送给精配准模块400,将粗配准变换参数传送至变换模块500;
精配准模块400,用于基于B样条变换的精配准,并将精配准变换参数传送至变换模块500;
变换模块500,根据粗配准和精配准变换参数对Label图像作出变换,并将变换结果传送至输出模块600,变换结果即自动分割结果;
输出模块600,用于输出自动分割结果,分割结果中不同的标记符即对应于心脏中不同的部位。
所述粗配准模块300依照步骤(3.1)~(3.8)进行处理:
(3.1)将输入心脏CT序列图像模板作为浮动图像,将待分割的心脏CT图像作为参考图像,把浮动图像和参考图像分解为若干个不同的分辨率层;
(3.2)在当前分辨率层上对浮动图像进行仿射变换,并记录变换参数;
(3.3)对变换后的浮动图像进行插值;
(3.4)对参考图像与插值后的浮动图像进行相似性度量;
(3.5)利用优化算法搜寻下一次变换参数;
(3.6)重复步骤(3.2)~(3.5),直至优化算法的迭代次数到达最大值;
(3.7)用当前辨率层的参数结果对下一分辨率层的变换参数进行初始化;
(3.8)在所有分辨率层上重复步骤(3.2)~(3.7),得到基于仿射变换的粗配准结果和粗配准变换参数。
所述精配准模块400依照步骤(4.1)~(4.8)进行处理:
(4.1)将基于仿射变换的粗配准的结果分解为若干个不同的分辨率层,并将其作为新的浮动图像;
(4.2)在当前分辨率层上对新的浮动图像进行B样条变换,并记录变换参数;
(4.3)对变换后的浮动图像进行插值;
(4.4)对参考图像与插值后的浮动图像进行相似性度量;
(4.5)通过优化算法搜寻下一次变换的参数;
(4.6)重复步骤(4.2)~(4.5),直至优化算法的迭代次数到达最大值;
(4.7)用当前辨率层的参数结果对下一分辨率层的变换参数进行初始化;
(4.8)在所有分辨率层上重复步骤(4.2)~(4.7),得到精配准变换参数。
所述粗配准模块300和精配准模块400的多分辨率分解采用金字塔算法;所述插值方法采用B样条插值法;所述相似性度量采用互信息度量法;所述优化算法采用自适应随机梯度下降法。
本发明与现有技术相比具有如下优势:
(1)把医学专家的临床经验知识与自动分割过程相结合,在医学专家的临床经验知识的指导下实现全自动分割。
(2)通过一次自动分割即可同时分割并标记心脏各腔室。
附图说明
图1为基于配准的CT图像全心脏自动分割系统的结构图。
图2为基于配准的CT图像全心脏自动分割系统配准方法的流程图。
图3为几种Sigmoid函数的函数曲线。
图4为Label自动分割结果效果图,图中不同值对应于心脏不同部分。
图5为全心脏的分割结果效果图。
图6为主动脉的分割结果效果图。
图7为左心房的分割结果效果图。
图8为左心室腔的分割结果效果图。
图9为左心室心肌的分割结果效果图。
图10为右心房的分割结果效果图。
图11为右心室的分割结果效果图。
具体实施方式
为了更好地理解本发明的技术方案,以下结合附图和实施例进行详细的描述,但本发明的实施和保护范围不限于此。
实施例
本实施例的序列模板图像和待分割的心脏图像分辨率为512*512。
如图1所示,本发明提出一种基于配准的CT图像全心脏自动分割系统,其工作流程如下:
第一步,输入模块100输入心脏CT序列图像模板和待分割的心脏CT图像。
第二步,输入模块100将心脏CT模板图像传送至分割模块200,在CT序列图像模板上分割出心脏的各个部位,为各部位标记上标记符,制作出Label图像,并传送至变换模块500;该Label图像能应用于所有分割实例中,无需重复制作。
第三步,输入模块100将心脏CT序列图像模板和待分割的心脏CT图像传送至粗配准模块300,进行粗配准处理,具体步骤如下:
(3.1)将输入心脏CT序列图像模板作为浮动图像,将待分割的心脏CT图像作为参考图像,把浮动图像和参考图像分解为6不同的分辨率层。
(3.2)在当前分辨率层上对浮动图像进行仿射变换,并记录变换参数;
所述仿射变换可以分解为平移变换、旋转变换以及缩放变换。其数学表达式为:
T(x)=RSx+t (1)
其中x ∈Ω为变换前的点,Ω为图像域,x=[xyz]T,R为旋转变换矩阵,S为缩放变换矩阵,t为平移变换矩阵。
平移变换矩阵t的表达式为:
其中tx,ty,tz分别为沿x,y,z轴方向平移的距离。
缩放变换矩阵S的表达式为:
其中Sx,Sy,Sz分别为x,y,z轴方向上的缩放比例。
旋转变换又可分解为绕x,y,z轴旋转,绕各轴旋转的旋转变换矩阵分别为:
其中α,β,θ分别为绕x,y,z轴旋转的角度。
整体旋转是绕x,y,z轴旋转的串联,可表示为:
R=RxRyRz (7)
把(2)~(7)式代入(1)式可得:
由(8)式可知,仿射变换有12个参数,即μ=[a11,a12,a13,a21,a22,a23,a31,a32,a33,tx,ty,tz],所述仿射变换根据这12个参数对浮动图像进行变换。
首次仿射变换时,把[tx,ty,tz,Sx,Sy,Sz,α,β,θ]设置为[0,0,0,1,1,1,0,0,0],并通过(1)~(8)式求取a11~a33,然后根据(8)式进行仿射变换;此后的各次迭代中这12个参数由优化过程求取,并根据(8)式进行仿射变换。
(3.3)对变换后的浮动图像进行B样条插值。
在图像中只能直接获取整数点上的值。当一个变换将一个点从一个空间映射到另一个空间时,一般情况下会被映射到非整数点区域上。此时,需要通过插值来计算该点的值。
本发明中的插值方法采用B样条插值法。假设图像f(x)由采样点集fi=f(xi)描述,其中xi∈Ω是有整数间距的点,Ω为图像域。利用B样条基函数实现插值,非整数点上的像素值利用下式插值求得:
式中,x为图像域中任意实值像素位置,x=[xyz]T;xi为整数点位置的坐标矢量,xi=[xiyizi]T;ci为B样条的系数,利用递归滤波器计算;β3(x)为可分离的三次B样条卷积核,β3(x)=β3(x)β3(y)β3(z)。其中:
(3.4)对参考图像与插值后的浮动图像进行相似性度量。
本发明采用互信息度量法,它是基于信息理论的。互信息通常用于描述两个系统间的相似性,它可通过熵进行计算。熵是对两个系统间的混乱程度的度量,互信息越大熵越小,两图像间的相似性越大。
假设参考图像和浮动图像的像素的值可看成两个随机变量F和M,则变量F和M的熵H(F)、H(M)及其联合熵H(F,M)可表示为:
其中f∈F,m∈M,PF(f)和PM(m)分别是F和M完全独立时的概率分布,PFM(f,m)是F和M的联合概率分布。
利用下式计算参考图像和浮动图像间的互信息:
I(F,M)=H(F)+H(M)-H(F,M) (14)
(3.5)利用优化算法搜寻下一次变换参数。
本发明所用优化算法为自适应随机梯度下降法。假设x为图中的坐标,可表示为x=[x,y,z]T;若配准中的参数化变换公式为T(x,μ),ΩF为参考图像的图像域,采用下式作为代价函数:
其中:
Ψ(u)=u (17)
ξ(u,v)=(u-v)2 (18)
因此得出代价函数的梯度函数:
然后利用以下迭代优化策略寻找最佳参数:
μk+1=μk-γ(tk)g(μk),k=0,1,Λ,K (20)
tk+1=[tk+f(-g(μk)Tg(μk-1))]+ (21)
γ(tk)=a/(tk+1+A)α (22)
式(20)中g(μk)为代价函数C在μk处的梯度,γ为步长;式(21)中[x]+=max(x,0);其中t0=0;f为Sigmoid函数,其中几种Sigmoid函数的曲线如图3所示。式(22)中a>0,A≥1,0<α≤1;式(22)中tk与当前梯度g(μk)和前一梯度g(μk-1)的内积有关,如果两个连续迭代步骤的梯度同向,则内积为正,-g(μk)Tg(μk-1)为负,如图3所示,当-g(μk)Tg(μk-1)为负时,f(-g(μk)Tg(μk-1))为负,因此tk+1相对于tk减少。从式(22)可见,当tk+1减少,则步长γ(tk+1)增大。自适应随机梯度下降法正是利用这种方法使得迭代步长得到自动调节。
在式(22)中取A=20,α=1。利用下式计算a:
其中ε为逼近代价函数梯度g时产生的误差,E为对g作出N次估计时的期望,aMAX、E||g||2、E||εk||2利用下式计算:
式(25)中δ为参考图像和浮动图像各方向像素间距的平均值,Ji为空间变换的雅可比矩阵,利用式(28)进行计算:
式(25)~(27)中tr(X)为矩阵X的迹,C利用下式进行计算:
式(25)~(27)中σ1、σ2利用下式计算:
式(30)(31)中μn是从μk附近随机抽取的估计参数,其服从正态分布:
其中I为单位矩阵,σ3利用下式计算:
(33)式中首先求取分母的最大值,当分母最大时即能求得σ3。
式(30)(31)中N为对梯度g作出估计的次数,估计的次数N满足下列各式:
联立(34)~(36)得:
若(37)式右边小于2,则N=2,否则N利用下式计算:
其中K=1.5。
把各参数代入(25)式中,并求的最大值,则能求得aMAX,再通过式(24)求参数a。
式(23)中fMAX取1,fMIN及ω利用下式求取:
至此,则可利用式(20)~(23)求取各次迭代的变换参数。
(3.6)重复步骤(3.2)~(3.5),直至优化算法的迭代次数到达250次。
(3.7)用当前辨率层的参数结果对下一分辨率层的变换参数进行初始化。
(3.8)在所有分辨率层上重复步骤(3.2)~(3.7),得到基于仿射变换的粗配准结果和粗配准变换参数,粗配准结果传送至精配准模块400,粗配准变换参数传送至变换模块500。
第四步,粗配准模块300将粗配准结果传送至精配准模块400,进行精配准处理,具体步骤如下:
(4.1)将基于仿射变换的粗配准的结果分解为6个不同的分辨率层,并将其作为新的浮动图像。
(4.2)在当前分辨率层上对新的浮动图像进行B样条变换,并记录变换参数。
本发明利用B样条变换进行精配准,它在形变中采用了B样条基函数。由于改变某一控制点的位置时只会影响有限的一段B样条曲线,因此改变局部控制点时不会影响整体的形变,只对局部形变产生影响,所以能利用它进行局部精配准。
假设有三维体数据V={(x,y,z)|0≤x<X,0≤y<Y,0≤z<Z},Φ是一组网格控制点φi,j,k的集合,其大小为nx×ny×nz,它们在三个方向上的间距分别为δx,δy,δz。利用三次B样条的三维张量积表示形变模型:
其中:
(41)式中φi+l,j+m,k+n为网格点(i+l,j+m,k+n)对应的坐标,它就是B样条变换的参数μ,把初始参数设置为网格点间距为16mm的网格。
B样条的基函数为:
B0(u)=(1-u)3/6
B1(u)=(3u3-6u2+4)/6 (44)
B2(u)=(1-3u3+3u2+3u)/6
B3(u)=u3/6
(4.3)利用如步骤(3.3)所述的方法对变换后的浮动图像进行B样条插值。
(4.4)利用如步骤(3.4)所述的方法对参考图像与插值后的浮动图像进行相似性度量。
(4.5)利用如步骤(3.5)所述的方法通过优化算法搜寻下一次迭代的变换参数。
(4.6)重复步骤(4.2)~(4.5),直至优化算法的迭代次数到达500次。
(4.7)用当前辨率层的参数结果对下一分辨率层的变换参数进行初始化。
(4.8)在所有分辨率层上重复步骤(4.2)~(4.7),得到精配准变换参数,传送至变换模块500。至此完成了整个配准处理过程。
第五步,变换模块500根据上述步骤得到的变换参数对Label图像进行变换。变换参数记录了配准过程中的浮动图像的变换过程,因此把它应用于Label图像即能重现浮动图像的变换过程。将变换结果传送至输出模块600,变换结果即自动分割结果。
第六步,输出模块600输出自动分割结果,分割结果中不同的标记符则对应于心脏中不同的部位。
至此,即实现整个基于配准的CT图像全心脏自动分割,并得到如图4至图11所示的分割结果。
Claims (3)
1.一种基于配准的CT图像全心脏自动分割系统,其特征在于包括:
输入模块(100),用于输入三维的心脏CT序列图像模板和三维的待分割的心脏CT图像,并将心脏CT模板图像传送至分割模块(200),将输入的心脏CT序列图像模板和待分割的心脏CT图像传送至粗配准模块(300);
分割模块(200),用于在CT序列图像模板上分割出心脏的各个部位,为各部位标记上标记符,制作出Label图像,并传送至变换模块(500);该Label图像能应用于所有分割实例中,无需重复制作;
粗配准模块(300),用于基于仿射变换的粗配准,对模板作出平移、缩放、旋转变换,并将粗配准结果传送给精配准模块(400),将粗配准变换参数传送至变换模块(500);
精配准模块(400),用于基于B样条变换的精配准,并将精配准变换参数传送至变换模块(500);
变换模块(500),根据粗配准和精配准变换参数对Label图像作出变换,并将变换结果传送至输出模块(600),变换结果即自动分割结果;
输出模块(600)用于输出三维的自动分割结果,分割结果中不同的标记符即对应于心脏中不同的部位;
所述粗配准模块(300)依照步骤(3.1)~(3.8)进行处理:
(3.1)将输入的心脏CT序列图像模板作为浮动图像,将待分割的心脏CT图像作为参考图像,把浮动图像和参考图像分解为若干个不同的分辨率层;
(3.2)在当前分辨率层上对浮动图像进行仿射变换,并记录变换参数;
(3.3)对变换后的浮动图像进行插值;
(3.4)对参考图像与插值后的浮动图像进行相似性度量;
(3.5)利用优化算法搜寻下一次变换参数;
(3.6)重复步骤(3.2)~(3.5),直至优化算法的迭代次数到达最大值;
(3.7)用当前辨率层的参数结果对下一分辨率层的变换参数进行初始化;
(3.8)在所有分辨率层上重复步骤(3.2)~(3.7),得到基于仿射变换的粗配准结果和粗配准变换参数;
所述精配准模块(400)依照步骤(4.1)~(4.8)进行处理:
(4.1)将基于仿射变换的粗配准的结果分解为若干个不同的分辨率层,并将其作为新的浮动图像;
(4.2)在当前分辨率层上对新的浮动图像进行B样条变换,并记录变换参数;
(4.3)对变换后的浮动图像进行插值;
(4.4)对参考图像与插值后的浮动图像进行相似性度量;
(4.5)通过优化算法搜寻下一次变换的参数;
(4.6)重复步骤(4.2)~(4.5),直至优化算法的迭代次数到达最大值;
(4.7)用当前辨率层的参数结果对下一分辨率层的变换参数进行初始化;
(4.8)在所有分辨率层上重复步骤(4.2)~(4.7),得到精配准变换参数。
2.根据权利要求1所述的自动分割系统,其特征在于,对于所述粗配准模块(300)和精配准模块(400),其中多分辨率分解采用金字塔算法;所述插值方法采用B样条插值法;所述相似性度量采用互信息度量法;所述优化算法采用自适应随机梯度下降法。
3.根据权利要求2所述的自动分割系统,其特征在于将浮动图像和参考图像分解分解为6个不同的分辨率层。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110264063.4A CN102411780B (zh) | 2011-09-07 | 2011-09-07 | 一种基于配准的ct图像全心脏自动分割系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110264063.4A CN102411780B (zh) | 2011-09-07 | 2011-09-07 | 一种基于配准的ct图像全心脏自动分割系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102411780A CN102411780A (zh) | 2012-04-11 |
CN102411780B true CN102411780B (zh) | 2014-07-02 |
Family
ID=45913846
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110264063.4A Expired - Fee Related CN102411780B (zh) | 2011-09-07 | 2011-09-07 | 一种基于配准的ct图像全心脏自动分割系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102411780B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103854276B (zh) * | 2012-12-04 | 2018-02-09 | 东芝医疗系统株式会社 | 图像配准及分割装置和方法,以及医学图像设备 |
CN104240287B (zh) * | 2013-06-08 | 2017-10-20 | 北京思创贯宇科技开发有限公司 | 一种利用ct图像生成冠脉全景图的方法及系统 |
US9305358B2 (en) * | 2013-07-01 | 2016-04-05 | Kabushiki Kaisha Toshiba | Medical image processing |
CN104021547B (zh) * | 2014-05-17 | 2017-02-15 | 清华大学深圳研究生院 | 肺部 ct 的三维配准方法 |
CN105243657B (zh) * | 2015-09-08 | 2018-01-19 | 首都医科大学附属北京安贞医院 | 基于增强弹性形变的carto电解剖图与ct图像配准方法和装置 |
US20170337682A1 (en) | 2016-05-18 | 2017-11-23 | Siemens Healthcare Gmbh | Method and System for Image Registration Using an Intelligent Artificial Agent |
CN106991694B (zh) * | 2017-03-17 | 2019-10-11 | 西安电子科技大学 | 基于显著区域面积匹配的心脏ct与超声图像配准方法 |
CN107492104B (zh) * | 2017-07-24 | 2020-09-15 | 广东顺德中山大学卡内基梅隆大学国际联合研究院 | 低磁场磁共振图像的皮层下结构子区域的自动分割方法 |
CN108549906A (zh) * | 2018-04-10 | 2018-09-18 | 北京全域医疗技术有限公司 | 放疗勾靶图像配准方法及装置 |
CN109509193B (zh) * | 2018-10-19 | 2020-08-25 | 中南民族大学 | 一种基于高精度配准的肝脏ct图谱分割方法及系统 |
CN109614991A (zh) * | 2018-11-19 | 2019-04-12 | 成都信息工程大学 | 一种基于Attention的多尺度扩张性心肌的分割分类方法 |
CN111402221B (zh) * | 2020-03-11 | 2023-02-17 | 联想(北京)有限公司 | 一种图像处理方法、装置及电子设备 |
CN112990444B (zh) * | 2021-05-13 | 2021-09-24 | 电子科技大学 | 一种混合式神经网络训练方法、系统、设备及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630256A (zh) * | 2009-08-27 | 2010-01-20 | 浙江大学 | 多模医学影像信息处理系统及滤波器的实现方法 |
-
2011
- 2011-09-07 CN CN201110264063.4A patent/CN102411780B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101630256A (zh) * | 2009-08-27 | 2010-01-20 | 浙江大学 | 多模医学影像信息处理系统及滤波器的实现方法 |
Non-Patent Citations (4)
Title |
---|
医学图像配准算法研究;张红颖;《中国博士学位论文全文数据库 信息科技辑》;20090415(第4期);第2-5章 * |
基于配准的颅脑CT图像分割的研究;徐峰;《中国优秀硕士学位论文全文数据库 信息科技辑》;20110115(第1期);第4章 * |
张红颖.医学图像配准算法研究.《中国博士学位论文全文数据库 信息科技辑》.2009,(第4期),第2-5章. |
徐峰.基于配准的颅脑CT图像分割的研究.《中国优秀硕士学位论文全文数据库 信息科技辑》.2011,(第1期),第4章. |
Also Published As
Publication number | Publication date |
---|---|
CN102411780A (zh) | 2012-04-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102411780B (zh) | 一种基于配准的ct图像全心脏自动分割系统 | |
Toga | Brain warping | |
CN102880866B (zh) | 一种人脸特征提取方法 | |
Seifert et al. | Hierarchical parsing and semantic navigation of full body CT data | |
JP2022518446A (ja) | 深層学習に基づく医用画像検出方法及び装置、電子機器及びコンピュータプログラム | |
US9785858B2 (en) | Method and system for hierarchical parsing and semantic navigation of full body computed tomography data | |
CN102208109B (zh) | X射线图像和激光图像的异源图像配准方法 | |
Chen et al. | Parametric shape representation by a deformable NURBS model for cardiac functional measurements | |
CN104091337A (zh) | 一种基于PCA及微分同胚Demons的变形医学图像配准方法 | |
CN104881680A (zh) | 一种基于二维特征和三维特征融合的阿尔茨海默病及轻度认知功能障碍识别方法 | |
US9508140B2 (en) | Quantifying curvature of biological structures from imaging data | |
CN109325971A (zh) | 图像配准方法及装置 | |
CN104050666A (zh) | 基于分割的脑部mr图像配准方法 | |
CN104240287A (zh) | 一种利用ct图像生成冠脉全景图的方法及系统 | |
Xu et al. | Supervised action classifier: Approaching landmark detection as image partitioning | |
CN102663762B (zh) | 医学图像中对称器官的分割方法 | |
CN107507189A (zh) | 基于随机森林与统计模型的小鼠ct图像肾脏分割方法 | |
CN104545999B (zh) | 一种超声图像膀胱容积测量方法及装置 | |
CN111260702B (zh) | 激光三维点云与ct三维点云配准方法 | |
Xue et al. | Lung 4D CT image registration based on high-order markov random field | |
CN113902779B (zh) | 一种基于张量投票方法的点云配准方法 | |
CN113706514B (zh) | 基于模板图像的病灶定位方法、装置、设备及存储介质 | |
CN116763295B (zh) | 牲畜体尺测量方法、电子设备及存储介质 | |
CN112734814B (zh) | 三维颅面锥形束ct图像配准方法 | |
CN106251359A (zh) | 基于Clifford代数几何相对不变量的3D图像配准方法 |
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: 20140702 Termination date: 20210907 |
|
CF01 | Termination of patent right due to non-payment of annual fee |