CN103337065A - 小鼠三维ct图像的非刚性配准方法 - Google Patents
小鼠三维ct图像的非刚性配准方法 Download PDFInfo
- Publication number
- CN103337065A CN103337065A CN201310193801XA CN201310193801A CN103337065A CN 103337065 A CN103337065 A CN 103337065A CN 201310193801X A CN201310193801X A CN 201310193801XA CN 201310193801 A CN201310193801 A CN 201310193801A CN 103337065 A CN103337065 A CN 103337065A
- Authority
- CN
- China
- Prior art keywords
- image
- mouse
- registration
- unique point
- target image
- 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
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
小鼠三维CT图像的非刚性配准方法,所述方法包括:A对小鼠三维CT图像的源图像的方向位置进行校正;B用阈值分割算法分割小鼠图像,得到源图像和目标图像的小鼠骨骼,分别提取骨骼切片连通域的质心作为特征点,得到源图像特征点与目标图像特征点;C对步骤B中源图像特征点与目标图像特征点进行匹配,由匹配后的特征点计算薄板样条变换的矩阵,并采用分层变换的方式对小鼠源图像进行薄板样条变换;D对薄板样条变换后的小鼠源图像与目标图像进行精细配准。本发明提取图像特征点时减少了特征点的数量并且采用分层变换的方式,克服了配准时间过长的缺点,在配准图像时基本能够实现自动配准,配准三维小鼠图像更加方便。
Description
技术领域
本发明涉及医学图像处理领域,更进一步涉及医学图像的配准方法。本发明可以用于小鼠三维CT图像的配准。
背景技术
医学图像配准在医学图像分析和临床实验中都有着重要的作用,可以用于图像分割、图像融合、手术导航等技术。在医学实验中,小鼠常作为研究的对象,在分析和处理小鼠三维图像前,图像配准是很有必要的。但是小鼠的身体结构比较复杂,不同的图像之间经常有不规则的非刚性形变,而且软组织与骨骼的形变程度是不一样的,在配准的时候需要将软组织与骨骼区别对待,否则不能得到很好的配准结果。传统的刚性配准方法并不能很好的完成这种配准任务,必须采用非刚性配准方法,但是这样会造成配准时间太长,实用性不强等问题。
Xia Li等人在“Automatic nonrigid registration of whole body CT mice images.Medical Physics35(4)20081507–1520.”中提出一种小鼠整体三维CT图像的配准方法。该方法首先利用基于仿射变换的配准方法对小鼠在图像中的相对位置进行校正,然后利用点配准算法对小鼠图像进行粗配准,最后利用一种非刚性配准算法对小鼠图像进行精细配准。该方法配准后的结果有非常高的配准精度,并且能够实现不同小鼠个体图像的配准。该方法的不足之处有二:其一是配准耗费的时间太长,没有太大的实用性;其二是不能用于小鼠身体形变太大的情况。
Martin Baiker等人在“Atlas-based whole-body segmentation of mice from low-contrast Micro-CT data.”中提出用小鼠三维CT图像与小鼠图谱配准的方法。该方法首先对小鼠图谱中的骨骼按照小鼠的解剖学结构手工分割每个骨头关节,在配准是先将小鼠CT图像的骨骼与图谱的骨骼按每个关节进行配准,然后配准小鼠CT图像与图谱的皮肤表面,得到的变换应用于整个小鼠的软组织。该方法的配准结果有很高的配准精度,并且能够应用于小鼠器官的精确分割。该方法的不足之处:小鼠图谱骨骼的手工分割一般难以实现,配准的过程复杂。
发明内容
本发明的目的在于克服上述已有技术的不足,提出一种三维CT图像非刚性配准方法,用于实现小鼠三维CT图像的配准。
为实现上述目的,本发明1.小鼠三维CT图像的非刚性配准方法,所述方法包括:
A对小鼠三维CT图像的源图像的方向位置等进行校正;
B用阈值分割算法分割小鼠图像,得到源图像和目标图像的小鼠骨骼,分别提取骨骼切片连通域的质心作为特征点,得到源图像特征点与目标图像特征点;
C对步骤B中源图像特征点与目标图像特征点进行匹配,由匹配后的特征点计算薄板样条变换的矩阵,并采用分层变换的方式对小鼠源图像进行薄板样条变换。
D最后对薄板样条变换后的小鼠源图像与目标图像进行精细配准。
在上述技术方案的基础上,所述步骤C中匹配的算法采用TPS-RPM点匹配算法
在上述技术方案的基础上,所述步骤C中匹配的算法采用ICP点匹配算法。
本发明还提供小鼠三维CT图像的非刚性配准方法,包括如下步骤:
1)源图像刚性校准步骤,其包括:
对源图像进行刚性变换,使源图像中小鼠的头尾方向,小鼠仰卧姿势和小鼠的主轴线与目标图像一致;
(2)特征点提取步骤,其包括:
2a)采用阈值分割方法分割源图像和目标图像,得到源骨骼图像和目标骨骼图像;
2b)计算目标骨骼图像中每个与小鼠的主轴线垂直切面中连通域的面积,计算面积大于设定阈值的连通区域的质心,得到目标图像特征点;
2c)计算源骨骼图像中每个与小鼠主轴线垂直切面中连通域的面积,计算面积大于设定阈值的连通区域的质心,得到源图像特征点;
(3)特征点匹配步骤,其包括:
匹配源图像特征点与目标图像特征点,得到匹配特征点;
(4)薄板样条变换步骤,其包括:
4a)对源图像进行采样间隔为4*4*4的下采样,得到下采样图像;
4b)利用源图像特征点和匹配特征点,对下采样图像进行薄板样条变换,得到配准采样图像;
4c)计算配准采样图像在薄板样条变换时每个坐标的偏移量,得到薄板样条偏移场;
4d)采用三次线性插值方法对薄板样条偏移场进行采样间隔为4*4*4的上采样插 值,得到配准偏移场;
4e)利用配准偏移场对源图像进行变换,得到点配准图像;
(5)精细配准步骤,其包括:
对点配准图像和目标图像进行精细配准,得到灰度配准图像,配准过程结束。
在上述技术方案的基础上,步骤(3)中匹配源图像特征点与目标图像特征点采用TPS-RPM点匹配算法,其按如下步骤进行:
第一步:选取模拟退火算法的退火系数,源图像特征点作为待匹配特征点;
第二步:利用下式计算待匹配特征点与目标图像特征点的对应矩阵:
其中,M表示待匹配特征点与目标图像特征点的匹配矩阵;t表示退火系数;exp表示指数运算;x表示目标图像特征点的坐标;y表示待匹配特征点的坐标;T表示矩阵的转置运算;
第三步:计算待匹配特征点与目标图像特征点之间薄板样条变换的核矩阵,非刚性变换矩阵和刚性变换矩阵;
第四步:按照下式重新计算待匹配特征点的坐标:
X=x*d+K*c
其中,X表示重新计算后的待匹配特征点的坐标;x表示重新计算前的待匹配特征点的坐标;*表示矩阵相乘运算;d表示薄板样条变换的刚性变换矩阵;+表示矩阵相加运算;K表示薄板样条变换的核矩阵;c表示薄板样条变换的非刚性变换矩阵;
第五步:判断退火系数是否低于设定的阈值;如果是,待匹配特征点即为匹配特征点,结束匹配过程,否则按下式计算退火系数:
t=0.93*t0
其中,t表示退火系数;t0表示前一次迭代的退火系数;并返回第二步。
在上述技术方案的基础上,所述步骤(5)利用B样条非刚性配准算法对点配准图像和目标图像进行精细配准,其具体步骤如下:
第一步:将点配准图像作为灰度预配准图像;
第二步:对点配准图像和目标图像进行随机采样,每幅图像得到2000-5000个采 样点,计算两幅图像采样点之间的归一化互信息;
第三步:利用灰度预配准图像和目标图像采样点之间的归一化互信息,采用自适应梯度下降优化算法对灰度预配准图像与目标图像之间B样条变换参数进行梯度优化,得到最优B样条变换参数;
第四步:利用最优B样条变换参数对灰度预配准图像进行三次均匀B样条变换,采用一阶B样条差值方法对变换后的图像进行插值,将插值后的图像作为灰度预配准图像;
第五步:判断梯度优化时步长是否小于人工设定的最小步长,如果是,灰度预配准图像即是灰度配准图像,非刚性配准结束;如果否,继续返回第二步。
与现有技术相比较,本发明在提取图像特征点时通过采样减少了源图像的特征点的数量,在薄板样条变换时采用分层变换的方式,克服了现有技术中配准时间过长的缺点,使得本发明在保证配准精度的前提下,大大缩减了配准的时间,使用本发明配准三维图像更加快捷;
且在配准图像时基本能够实现自动配准,克服了现有技术中预处理复杂的缺点,使本发明易于实现,使用本发明配准三维图像更加方便。
附图说明
图1为本发明小鼠三维CT图像的非刚性配准方法流程图;
图2本发明Micro-CT扫描的目标图像;
图3为本发明Micro-CT扫描的需要配准的源图像;
图4为本发明目标图像阈值分割的小鼠骨骼图像;
图5为本发明源图像阈值分割的小鼠骨骼图像;
图6为本发明处理后的配准图像;
图7为本发明为配准结果图像阈值分割的小鼠骨骼图像。
具体实施方式
下面结合图1对本发明的具体实施步骤做进一步的详细描述。
步骤1.对源图像刚性校准
在刚性校准时,首先校正图像中小鼠的头部和尾部方向一致,然后校正小鼠仰卧姿势一致,最后小鼠的中轴线一致对齐。
步骤2.提取源图像和目标图像的特征点
(1)选取阈值,采用阈值分割方法提取源图像和目标图像中小鼠的骨骼,得到源骨骼图像和目标骨骼图像,在本实施例中,阈值的大小是105;
(2)计算目标骨骼图像中每个与小鼠的主轴线垂直切面连通域的面积,计算面积大于10的连通区域的质心,得到目标图像特征点;
(3)计算源骨骼图像中每个与小鼠的主轴线垂直切面连通域的面积,计算面积大于10的连通区域的质心,得到源图像特征点。
步骤3.匹配源图像和目标图像的特征点
将源图像特征点与目标图像特征点进行匹配,得到匹配特征点;其算法可采用TPS-RPM点匹配算法和ICP点匹配算法,其中,采用TPS-RPM点匹配算其具体步骤如下:
第一步:选取模拟退火算法的退火系数,源图像特征点作为待匹配特征点;
第二步:利用下式计算待匹配特征点与目标图像特征点的对应矩阵:
其中,M表示待匹配特征点与目标图像特征点的匹配矩阵;t表示退火系数;exp表示指数运算;x表示目标图像特征点的坐标;y表示待匹配特征点的坐标;T表示矩阵的转置运算;
第三步:计算待匹配特征点与目标图像特征点之间薄板样条变换的核矩阵,非刚性变换矩阵和刚性变换矩阵;
第四步:按照下式重新计算待匹配特征点的坐标:
X=x*d+K*c
其中,X表示重新计算后的待匹配特征点的坐标;x表示重新计算前的待匹配特征点的坐标;*表示矩阵相乘运算;d表示薄板样条变换的刚性变换矩阵;+表示矩阵相加运算;K表示薄板样条变换的核矩阵;c表示薄板样条变换的非刚性变换矩阵;
第五步:判断退火系数是否低于设定的阈值;如果是,待匹配特征点即为匹配特征点,结束匹配过程,否则按下式计算退火系数:
t=0.93*t0
其中,t表示退火系数;t0表示前一次迭代的退火系数;并返回第二步。
步骤4.对源图像进行薄板样条变换
第一步:对源图像进行采样间隔为4*4*4的下采样,得到下采样图像;
第二步:利用源图像特征点和匹配特征点,对下采样图像进行薄板样条变换,得到配准采样图像;
第三步:计算下采样图像在薄板样条变换时每个坐标的偏移量,得到薄板样条偏移场;
第四步:采用三线性插值对薄板样条偏移场进行采样间隔为4*4*4的上采样,得到配准偏移场;
第五步:利用配准偏移场对源图像进行变换,得到点配准图像。
步骤5.对薄板样条变换后的源图像进行精细配准
利用改进的B样条非刚性配准算法对点配准图像和目标图像进行精细配准。其具体步骤如下:
第一步:将点配准图像作为灰度预配准图像;
第二步:对点配准图像和目标图像进行随机采样,每幅图像得到2000-5000个采样点,计算两幅图像采样点之间的归一化互信息;
第三步:利用灰度预配准图像和目标图像采样点之间的归一化互信息,采用自适应梯度下降优化算法对灰度预配准图像与目标图像之间B样条变换参数进行梯度优化,得到最优B样条变换参数;
第四步:利用最优B样条变换参数对灰度预配准图像进行三次均匀B样条变换,采用一阶B样条差值方法对变换后的图像进行插值,将插值后的图像作为灰度预配准图像;
第五步:判断梯度优化时步长是否小于设定的最小步长,如果是,灰度预配准图像即是灰度配准图像,非刚性配准结束;如果否,继续返回第二步。
下面结合附图2至附图7对本发明的配准结果做进一步的描述。
1.仿真条件
在CPU为Intel corei5、3.1GHZ、内存8G、Window7操作系统上使用MATLAB2009a,VS2010,3DMed4.1进行仿真
2.仿真内容
同一只小鼠由于在Micro-CT扫描时的时间不同,在扫描时的小鼠的姿势出现了变化,致使小鼠的骨骼和内脏器官出现了不可预知的变形。通过本发明对出现变形的图像与另一幅图像进行配准,使两幅图像能够共同分析和处理。
图2为Micro-CT扫描的目标图像;图3为Micro-CT扫描的需要配准的源图像。图4为目标图像阈值分割的小鼠骨骼;图5为源图像阈值分割的小鼠骨骼。图6为本发明处理后的配准图像;图7为配准图像阈值分割的小鼠骨骼。
3.仿真结果
根据图7可以看出本发明对小鼠发生姿势形变时的配准结果是很准确的,下面我们通过计算配准前后源图像和目标图像小鼠骨骼之间的相关系数和Dice系数来定量分析本发明的结果:
配准前 | 配准后 | |
相关系数 | 0.179 | 0.958 |
Dice系数 | 0.172 | 0.893 |
本发明的另一个优点是配准三维图像的时间大大的缩减,在仿真时所用的时间是:
点匹配 | 薄板样条变换 | B样条非刚性配准 | 总时间 |
40.6s | 206.33s | 72s | 318.93s |
以上仿真结果可以看出,采用三维CT图像非刚性配准方法,能够很好的对三维CT小鼠图像非刚性变形进行配准,在保证一定配准精度的前提下,缩短了配准的时间。
Claims (6)
1.小鼠三维CT图像的非刚性配准方法,其特征在于,所述方法包括:
A对小鼠三维CT图像的源图像的方向位置进行校正;
B用阈值分割算法分割小鼠图像,得到源图像和目标图像的小鼠骨骼,分别提取骨骼切片连通域的质心作为特征点,得到源图像特征点与目标图像特征点;
C对步骤B中源图像特征点与目标图像特征点进行匹配,由匹配后的特征点计算薄板样条变换的矩阵,并采用分层变换的方式对小鼠源图像进行薄板样条变换;
D对薄板样条变换后的小鼠源图像与目标图像进行精细配准。
2.如权利要求1所述的小鼠三维CT图像的非刚性配准方法,其特征在于:所述步骤C中匹配的算法采用TPS-RPM点匹配算法。
3.如权利要求1所述的小鼠三维CT图像的非刚性配准方法,其特征在于:所述步骤C中匹配的算法采用ICP点匹配算法。
4.小鼠三维CT图像的非刚性配准方法,其特征在于,包括如下步骤:
(1)源图像刚性校准步骤,其包括:
对源图像进行刚性变换,使源图像中小鼠的头尾方向,小鼠仰卧姿势和小鼠的主轴线与目标图像一致;
(2)特征点提取步骤,其包括:
2a)采用阈值分割方法分割源图像和目标图像,得到源骨骼图像和目标骨骼图像;
2b)计算目标骨骼图像中每个与小鼠的主轴线垂直切面中连通域的面积,提取面积大于设定阈值的连通区域的质心,得到目标图像特征点;
2c)计算源骨骼图像中每个与小鼠主轴线垂直切面中连通域的面积,提取面积大于设定阈值的连通区域的质心,得到源图像特征点;
(3)特征点匹配步骤,其包括:
匹配源图像特征点与目标图像特征点,得到匹配特征点;
(4)薄板样条变换步骤,其包括:
4a)对源图像进行采样间隔为4*4*4的下采样,得到下采样图像;
4b)利用源图像特征点和匹配特征点,对下采样图像进行薄板样条变换,得到配准采样图像;
4c)计算配准采样图像在薄板样条变换时每个坐标的偏移量,得到薄板样条偏移场;
4d)采用三次线性插值方法对薄板样条偏移场进行采样间隔为4*4*4的上采样插值,得到配准偏移场;
4e)利用配准偏移场对源图像进行变换,得到点配准图像;
(5)精细配准步骤,其包括:
对点配准图像和目标图像进行精细配准,得到灰度配准图像,配准过程结束。
5.如权利要求4中所述的小鼠三维CT图像的非刚性配准方法,其特征在于:步骤(3)中匹配源图像特征点与目标图像特征点采用TPS-RPM点匹配算法,其按如下步骤进行:
第一步:选取模拟退火算法的退火系数,将源图像特征点作为待匹配特征点;
第二步:利用下式计算待匹配特征点与目标图像特征点的对应矩阵:
其中,M表示待匹配特征点与目标图像特征点的匹配矩阵;t表示退火系数;exp表示指数运算;x表示目标图像特征点的坐标;y表示待匹配特征点的坐标;T表示矩阵的转置运算;
第三步:计算待匹配特征点与目标图像特征点之间薄板样条变换的核矩阵,非刚性变换矩阵和刚性变换矩阵;
第四步:按照下式重新计算待匹配特征点的坐标:
X=x*d+K*c
其中,X表示重新计算后的待匹配特征点的坐标;x表示重新计算前的待匹配特征点的坐标;*表示矩阵相乘运算;d表示薄板样条变换的刚性变换矩阵;+表示矩阵相加运算;K表示薄板样条变换的核矩阵;c表示薄板样条变换的非刚性变换矩阵;
第五步:判断退火系数是否低于设定的阈值;如果是,待匹配特征点即为匹配特征点,结束匹配过程,否则按下式计算退火系数:
t=0.93*t0
其中,t表示退火系数;t0表示前一次迭代的退火系数;并返回第二步。
6.如权利要求4中所述的小鼠三维CT图像的非刚性配准方法,其特征在于:所述步骤(5)利用B样条非刚性配准算法对点配准图像和目标图像进行精细配准,其具体步骤如下:
第一步:将点配准图像作为灰度预配准图像;
第二步:对点配准图像和目标图像进行随机采样,每幅图像得到2000-5000个采样点,计算两幅图像采样点之间的归一化互信息;
第三步:利用灰度预配准图像和目标图像采样点之间的归一化互信息,采用自适应梯度下降优化算法对灰度预配准图像与目标图像之间B样条变换参数进行梯度优化,得到最优B样条变换参数;
第四步:利用最优B样条变换参数对灰度预配准图像进行三次均匀B样条变换,采用一阶B样条差值方法对变换后的图像进行插值,将插值后的图像作为灰度预配准图像;
第五步:判断梯度优化时步长是否小于设定的最小步长,如果是,灰度预配准图像即是灰度配准图像,非刚性配准结束;如果否,继续返回第二步。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310193801.XA CN103337065B (zh) | 2013-05-22 | 2013-05-22 | 小鼠三维ct图像的非刚性配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310193801.XA CN103337065B (zh) | 2013-05-22 | 2013-05-22 | 小鼠三维ct图像的非刚性配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103337065A true CN103337065A (zh) | 2013-10-02 |
CN103337065B CN103337065B (zh) | 2016-03-02 |
Family
ID=49245213
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310193801.XA Expired - Fee Related CN103337065B (zh) | 2013-05-22 | 2013-05-22 | 小鼠三维ct图像的非刚性配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103337065B (zh) |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103632371A (zh) * | 2013-12-06 | 2014-03-12 | 河海大学常州校区 | 基于兼容性网格分割的骨骼参数计算方法 |
CN103745467A (zh) * | 2014-01-07 | 2014-04-23 | 东南大学 | 一种基于数字体积相关法的三维图像配准方法 |
CN105118021A (zh) * | 2015-07-28 | 2015-12-02 | 深圳市汇川技术股份有限公司 | 基于特征点的图像配准方法及系统 |
CN105139382A (zh) * | 2015-08-03 | 2015-12-09 | 华北电力大学(保定) | 一种冠状动脉内超声图像序列的弹性配准方法 |
CN106353354A (zh) * | 2015-07-14 | 2017-01-25 | 艾斯拜克特成像有限公司 | 在高温高压下样品和过程的核磁共振成像的装置和方法 |
CN106447613A (zh) * | 2016-09-27 | 2017-02-22 | 西安蒜泥电子科技有限责任公司 | 一种基于图像局部配准的全景去虚影方法及系统 |
CN106934821A (zh) * | 2017-03-13 | 2017-07-07 | 中国科学院合肥物质科学研究院 | 一种基于icp算法和b样条的锥形束ct和ct图像配准方法 |
CN108154525A (zh) * | 2017-11-21 | 2018-06-12 | 四川大学 | 一种基于特征匹配的骨骼碎片拼接方法 |
CN109754396A (zh) * | 2018-12-29 | 2019-05-14 | 上海联影智能医疗科技有限公司 | 图像的配准方法、装置、计算机设备和存储介质 |
CN109934861A (zh) * | 2019-01-22 | 2019-06-25 | 广东工业大学 | 一种头颈部多模态医学图像自动配准方法 |
CN110473241A (zh) * | 2019-07-31 | 2019-11-19 | 上海联影智能医疗科技有限公司 | 图像配准方法、存储介质及计算机设备 |
WO2020024135A1 (zh) * | 2018-08-01 | 2020-02-06 | 深圳市大疆创新科技有限公司 | 图像配准的方法、装置、计算机系统和可移动设备 |
CN111487612A (zh) * | 2020-04-19 | 2020-08-04 | 中国人民解放军海军航空大学 | 基于cpd的异地配置雷达/esm航迹抗差关联方法 |
CN111602173A (zh) * | 2017-10-23 | 2020-08-28 | 布莱恩欧米克斯有限公司 | 断层扫描数据分析 |
CN114404039A (zh) * | 2021-12-30 | 2022-04-29 | 华科精准(北京)医疗科技有限公司 | 三维模型的组织漂移校正方法、装置、电子设备及存储介质 |
CN114742869A (zh) * | 2022-06-15 | 2022-07-12 | 西安交通大学医学院第一附属医院 | 基于图形识别的脑部神经外科配准方法及电子设备 |
CN115100661A (zh) * | 2022-08-24 | 2022-09-23 | 上海安维尔信息科技股份有限公司 | 一种集装箱箱号提取方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1772827A2 (en) * | 2005-10-05 | 2007-04-11 | Siemens Medical Solutions USA, Inc. | Method and apparatus for the rigid registration of 3D ear impression shapes with skeletons |
CN102100565A (zh) * | 2009-12-17 | 2011-06-22 | 通用电气公司 | 使用非刚性配准校正门控pet图像中运动的系统和方法 |
-
2013
- 2013-05-22 CN CN201310193801.XA patent/CN103337065B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1772827A2 (en) * | 2005-10-05 | 2007-04-11 | Siemens Medical Solutions USA, Inc. | Method and apparatus for the rigid registration of 3D ear impression shapes with skeletons |
CN102100565A (zh) * | 2009-12-17 | 2011-06-22 | 通用电气公司 | 使用非刚性配准校正门控pet图像中运动的系统和方法 |
Non-Patent Citations (3)
Title |
---|
DI XIAO ET AL.: "Non-rigid registration of small animal skeletons from micro-CT using 3D shape context", 《MEDICAL IMAGING 2009:IMAGE PROCESSING》, vol. 7259, 27 March 2009 (2009-03-27), pages 1 - 8 * |
MARTIN BAIKER ET AL.: "Atlas-based whole-body segmentation of mice from low-contrast Micro-CT data", 《MEDICAL IMAGE ANALYSIS》, vol. 14, no. 6, 31 December 2010 (2010-12-31), pages 723 - 737, XP027141535 * |
XIA LI ET AL.: "Automatic nonrigid registration of whole body CT mice images", 《MEDICAL PHYSICS》, vol. 35, no. 4, 30 April 2008 (2008-04-30), pages 1507 - 1520, XP012116009, DOI: doi:10.1118/1.2889758 * |
Cited By (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103632371A (zh) * | 2013-12-06 | 2014-03-12 | 河海大学常州校区 | 基于兼容性网格分割的骨骼参数计算方法 |
CN103745467A (zh) * | 2014-01-07 | 2014-04-23 | 东南大学 | 一种基于数字体积相关法的三维图像配准方法 |
CN103745467B (zh) * | 2014-01-07 | 2016-08-17 | 东南大学 | 一种基于数字体积相关法的三维图像配准方法 |
CN106353354A (zh) * | 2015-07-14 | 2017-01-25 | 艾斯拜克特成像有限公司 | 在高温高压下样品和过程的核磁共振成像的装置和方法 |
CN105118021A (zh) * | 2015-07-28 | 2015-12-02 | 深圳市汇川技术股份有限公司 | 基于特征点的图像配准方法及系统 |
CN105118021B (zh) * | 2015-07-28 | 2019-04-19 | 深圳市汇川技术股份有限公司 | 基于特征点的图像配准方法及系统 |
CN105139382A (zh) * | 2015-08-03 | 2015-12-09 | 华北电力大学(保定) | 一种冠状动脉内超声图像序列的弹性配准方法 |
CN106447613A (zh) * | 2016-09-27 | 2017-02-22 | 西安蒜泥电子科技有限责任公司 | 一种基于图像局部配准的全景去虚影方法及系统 |
CN106934821B (zh) * | 2017-03-13 | 2020-06-23 | 中国科学院合肥物质科学研究院 | 一种基于icp算法和b样条的锥形束ct和ct图像配准方法 |
CN106934821A (zh) * | 2017-03-13 | 2017-07-07 | 中国科学院合肥物质科学研究院 | 一种基于icp算法和b样条的锥形束ct和ct图像配准方法 |
US10902596B2 (en) | 2017-10-23 | 2021-01-26 | Brainomix Limited | Tomographic data analysis |
CN111602173A (zh) * | 2017-10-23 | 2020-08-28 | 布莱恩欧米克斯有限公司 | 断层扫描数据分析 |
CN108154525B (zh) * | 2017-11-21 | 2022-06-07 | 四川大学 | 一种基于特征匹配的骨骼碎片拼接方法 |
CN108154525A (zh) * | 2017-11-21 | 2018-06-12 | 四川大学 | 一种基于特征匹配的骨骼碎片拼接方法 |
CN110945566A (zh) * | 2018-08-01 | 2020-03-31 | 深圳市大疆创新科技有限公司 | 图像配准的方法、装置、计算机系统和可移动设备 |
WO2020024135A1 (zh) * | 2018-08-01 | 2020-02-06 | 深圳市大疆创新科技有限公司 | 图像配准的方法、装置、计算机系统和可移动设备 |
CN109754396A (zh) * | 2018-12-29 | 2019-05-14 | 上海联影智能医疗科技有限公司 | 图像的配准方法、装置、计算机设备和存储介质 |
CN109934861A (zh) * | 2019-01-22 | 2019-06-25 | 广东工业大学 | 一种头颈部多模态医学图像自动配准方法 |
CN109934861B (zh) * | 2019-01-22 | 2022-10-18 | 广东工业大学 | 一种头颈部多模态医学图像自动配准方法 |
CN110473241B (zh) * | 2019-07-31 | 2022-03-01 | 上海联影智能医疗科技有限公司 | 图像配准方法、存储介质及计算机设备 |
CN110473241A (zh) * | 2019-07-31 | 2019-11-19 | 上海联影智能医疗科技有限公司 | 图像配准方法、存储介质及计算机设备 |
CN111487612A (zh) * | 2020-04-19 | 2020-08-04 | 中国人民解放军海军航空大学 | 基于cpd的异地配置雷达/esm航迹抗差关联方法 |
CN114404039A (zh) * | 2021-12-30 | 2022-04-29 | 华科精准(北京)医疗科技有限公司 | 三维模型的组织漂移校正方法、装置、电子设备及存储介质 |
CN114742869A (zh) * | 2022-06-15 | 2022-07-12 | 西安交通大学医学院第一附属医院 | 基于图形识别的脑部神经外科配准方法及电子设备 |
CN114742869B (zh) * | 2022-06-15 | 2022-08-16 | 西安交通大学医学院第一附属医院 | 基于图形识别的脑部神经外科配准方法及电子设备 |
CN115100661A (zh) * | 2022-08-24 | 2022-09-23 | 上海安维尔信息科技股份有限公司 | 一种集装箱箱号提取方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103337065B (zh) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103337065A (zh) | 小鼠三维ct图像的非刚性配准方法 | |
US11373303B2 (en) | Systems and methods for ultrasound imaging | |
Huang et al. | Anatomical prior based vertebra modelling for reappearance of human spines | |
CN102722890B (zh) | 基于光流场模型的非刚性心脏图像分级配准方法 | |
CN106447764B (zh) | 一种三维可视化的人体骨盆参数自动测量方法 | |
MXPA02001035A (es) | Sistema y metodo de fusion/alineacion de imagen automatizada. | |
CN108830852B (zh) | 三维超声肿瘤辅助测量系统及方法 | |
Hacihaliloglu et al. | Automatic adaptive parameterization in local phase feature-based bone segmentation in ultrasound | |
Zhang et al. | 3-D reconstruction of the spine from biplanar radiographs based on contour matching using the hough transform | |
CN107708550A (zh) | 用于检测和测量解剖异常的分段回声结构的表面建模 | |
CN103544695B (zh) | 一种高效的基于博弈框架的医学图像分割方法 | |
US20200121279A1 (en) | Method and device for measuring spinal column curvature | |
CN102663738A (zh) | 三维图像配准方法及系统 | |
CN104867104A (zh) | 基于xct图像非刚度配准的目标鼠解剖结构图谱获取方法 | |
CN102360495A (zh) | 基于平均密度投影和平移高斯模型的肺结节分割方法 | |
CN106780715A (zh) | 一种骨盆ct三维坐标体系的建立方法 | |
CN109994188A (zh) | 基于ndi的神经外科导航配准测试方法和系统 | |
Gutiérrez-Becker et al. | Automatic segmentation of the fetal cerebellum on ultrasound volumes, using a 3D statistical shape model | |
CN108670301B (zh) | 一种基于超声影像的脊柱横突定位方法 | |
CN111627011A (zh) | 基于ct值椎弓根钉置钉钉道选择方法、装置和设备 | |
Hacihaliloglu et al. | Statistical shape model to 3D ultrasound registration for spine interventions using enhanced local phase features | |
CN109692041A (zh) | 一种针对覆盖软骨的骨表面配准方法 | |
EP2742483B1 (en) | Image processing method | |
Alam et al. | Evaluation of medical image registration techniques based on nature and domain of the transformation | |
CN108670302B (zh) | 一种基于2.5维超声宽景成像的脊柱三维结构再现方法 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160302 Termination date: 20210522 |