CN102169578B - 一种基于有限元模型的非刚性医学图像配准方法 - Google Patents
一种基于有限元模型的非刚性医学图像配准方法 Download PDFInfo
- Publication number
- CN102169578B CN102169578B CN 201110071997 CN201110071997A CN102169578B CN 102169578 B CN102169578 B CN 102169578B CN 201110071997 CN201110071997 CN 201110071997 CN 201110071997 A CN201110071997 A CN 201110071997A CN 102169578 B CN102169578 B CN 102169578B
- Authority
- CN
- China
- Prior art keywords
- image
- registration
- images
- reference picture
- finite element
- 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
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000008569 process Effects 0.000 claims abstract description 20
- 238000007667 floating Methods 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 5
- 238000005381 potential energy Methods 0.000 claims description 23
- 239000004744 fabric Substances 0.000 claims description 12
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 238000009499 grossing Methods 0.000 claims description 9
- 230000008878 coupling Effects 0.000 claims description 5
- 238000010168 coupling process Methods 0.000 claims description 5
- 238000005859 coupling reaction Methods 0.000 claims description 5
- 238000003384 imaging method Methods 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- RKTYLMNFRDHKIL-UHFFFAOYSA-N copper;5,10,15,20-tetraphenylporphyrin-22,24-diide Chemical compound [Cu+2].C1=CC(C(=C2C=CC([N-]2)=C(C=2C=CC=CC=2)C=2C=CC(N=2)=C(C=2C=CC=CC=2)C2=CC=C3[N-]2)C=2C=CC=CC=2)=NC1=C3C1=CC=CC=C1 RKTYLMNFRDHKIL-UHFFFAOYSA-N 0.000 claims description 2
- 230000004927 fusion Effects 0.000 abstract description 3
- 238000002474 experimental method Methods 0.000 abstract description 2
- 230000002349 favourable effect Effects 0.000 abstract 1
- 238000004422 calculation algorithm Methods 0.000 description 8
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003709 image segmentation Methods 0.000 description 2
- 238000012216 screening Methods 0.000 description 2
- 238000011524 similarity measure Methods 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 229920000535 Tan II Polymers 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000005587 bubbling Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000002405 diagnostic procedure Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 230000016507 interphase Effects 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Image Processing (AREA)
- Image Analysis (AREA)
Abstract
本发明属于医学图像处理领域,涉及一种基于有限元模型非刚性医学图像配准方法。本发明所涉及配准方法在配准过程中采用分级配准的思想,先对浮动图像和参考图像进行全局配准,再对全局配准图像和参考图像进行有限元建模,最后计算与该模型相关的能量项,包内部能量项和驱动模型发生形变的外部能量项。在内部能量项和外部能量项的共同控制下模型的整体总能量达到最小时,模型收敛于平衡位置,实现图像局部精确配准,最终实现非刚性图像的配准。实验表明本方法能够能很好的模拟非刚性图像的局部形变,从而提高配准的精确度及速度、增强配准的鲁棒性,为医学图像融合奠定一个良好的基础。
Description
技术领域
本发明涉及一种图像处理技术领域的方法,具体是一种基于有限元模型的非刚性医学图像配准方法。
背景技术
医学图像配准是医学图像融合的基础,也是医学图像分析一项重要的技术,其主要通过寻找一种空间变换,使得两幅医学图像上的对应点达到空间位置或者解剖点上的一致性。在医学诊断过程中,由于存在不同模式图像表现不同性质的物理机制、患者的移动、成像参数的变化、不同成像设备间空间分辨率不相同等现实问题,所以单单凭借医生手动将两张相同或不同模式的图像在空间上作配准会受到很多局限,且常带有较大的主观性,其可靠性往往不高,不可避免地会产生误差。特别是在定向放射外科和心脏手术可视化等应用领域,对于图像配准的精度要求很高,使得医学图像配准成为一项必要而有相当困难的任务,诸如需要进一步解决刚性变换或者放射变换等全局配准方法通常无法应付局部变换的缺点。为此本文提出了一种基于有限元模型的非刚性医学图像配准方法。
医学图像配准方法的步骤主要包括:确定待配准图像与目标图像的一种空间变换;确定经过空间变换后的图像与目标图像的相似性测度;寻找一种参数优化策略使两幅图像的相似性测度达到最大。现有配准方法主要分为两大类:基于图像角点医学图像配准方法与基于医学图像灰度信息的配准方法。两者技术缺陷主要为:
基于图像特征的医学图像配准方法的技术缺陷为:它需要对图像进行分割来提取图像的特征,由于非刚性组织的结构很复杂,有些分界面不是很明显,通常需要人工预选定特征,这样会费时也费力而且配准的精度受分割精度影响,一般很难自动完成,使得配准时间过程、速度慢、配准精确不高。
基于图像灰度信息的医学图像配准方法的技术缺陷为:它不需要对图像进行分割处理,直接对整幅图像进行运算,会造成配准的速度较慢、配准时间长、鲁棒性差。
非刚性医学图像相对刚性医学图像来讲,由于其形变量较大较复杂,现有的配准算法不能很好的模拟其形变并克服配准时消耗时间多、精确度低、鲁棒性不强等缺点。
发明内容
本发明的目的在于提供一种能够提高配准的精确度及速度,增强配准鲁棒性的基于有限元模型的非刚性医学图像配准方法。
技术解决方案:
本发明一种基于有限元模型非刚性医学图像的配准方法,方法步骤如下:1)首先通过图像的一阶矩寻找浮动图像和参考图像这两幅待配准图像的质心;再次,通过图像的二阶中心矩寻找两幅待配准图像的主轴和坐标系的夹角;最后,通过平移和旋转使得两幅待配准图像的质心和主轴对齐,得到旋转分量,实现两幅图像的全局配准。
2)然后对全局配准图像及参考图像这两幅待配准图像进行网格划分,首先,对全局配准图像和参考图像分别进行高斯平滑滤波,消除图像成像过程中生成的噪声;然后分别对平滑处理后的两幅待配准图像利用Harris算子提取图像的角点,对其角点进行排序,其横坐标中最小值、最大值,取其纵坐标中最小值、最大值并记录下来,这样就形成四个点,连结四个角点形成一个四边形,使得四边形包含所有的角点。本发明将四边形进行划分成为横纵坐标相等的区域,将形成规则的四边形,接着按照正对角线连接四边形单元的对角线形成三角形网格,对三角形网格进行编号1、2、......、M,M为大于1的整数。
3)建立数学模型:将待配准图像的形变过程看作弹性的伸缩过程,利用线性弹性模型中的最小势能原理约束每个匹配后的形变,外力为两幅图像的灰度差值,当外力为零时有限元单元的势能到最小,图像完成精确配准。
4)建立有限元方程:在对步骤3)进行求解的过程中。由于对整幅图像进行建模,使得方程求解过程比较复杂,所以将浮动图像和参考图像进行离散化,将图像感兴趣的区域划分成规则形的三角形网格,将整幅图像的形变看作为三角形单元形变的总和。如步骤2)所述。将步骤2)全局配准图像和参考图像这两幅待配准图像网格单元建立有限元方程,首先,对两幅待配准图像的网格单元进行匹配,将匹配后的网格单元,进行编号1、2、......、N,N为大于1的整数:然后图像的总势能看作为三角形势能的总和,建立有限元方程:
其中,σt表示弹性体的应力,ε表示弹性体的应变,∫ΩσtεdΩ表示弹性***的形变势能,Ω表示二维图像区域,∫ΩFiudΩ表示弹性***的外力势能,Fi表示作用于非刚性医学图像的外力,u表示图像的形变位移,m表示三角形单元的数目,i表示第i个三角形;
5)求解方程:求解步骤4)中的有限元方程,在每个单元内选择至少一个节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,通过插值函数计算出各个单元内场函数的近似解;利用三角形单元的灰度差值的平方和作为外力,当总势能最小时,外力为零,来求解方程,得到全局配准图像的偏移量,实现图像的局部精确配准。
本发明将待配准的浮动图像和参考图像的质心和主轴对齐,实现两幅待配准图像的全局配准;然后根据有限分析中的变分原理将待配准图像划分成为不同的矩形区域,利用线性弹性模型最小势能原理约束图像的局部形变,实现非刚性医学图像局部精确配准,本发明能够能很好的模拟非刚性图像的局部形变,从而提高配准的精确度及速度、增强配准的鲁棒性,为医学图像融合奠定一个良好的基础。
附图说明
图1为本发明非刚性医学图像配准算法的流程图;
图2为本发明非刚性医学图像配准算法的三角形单元;
图3为本发明所用的浮动图像;
图4为本发明所用的参考图像;
图5为本发明全局配准图像;
图6为本发明配准算法的配准结果图;
具体实施方式
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1所示,本实施例包含以下步骤:
第一步:对两幅待配准图像进行全局配准
首先通过一阶矩寻找浮动图像和参考图像这两幅图像的质心,然后通过二阶中心矩寻找浮动图像和参考图像的主轴和坐标系的夹角,根据两幅图像的质心的位置差作为平移分量,然后计算各自的主轴方向并根据其差值将两幅图像旋转对齐,得到旋转分量,实现两幅待配准图像全局配准。具体如下步:
通过经典力学计算方法计算图像的质心如公式(1)所示:
其中:xi、yi代表像素点,pi代表像素的灰度值,n代表图像总共的像素点个数,Xcenter、Ycenter为图像的重心坐标。所述的图像的主轴通过以下方式获得:
其中: x,y代表图像中任一点的坐标,p(x,y)代表该点的灰度值,m取值范围为[0,2],n取值范围为[0,2],α为图像主轴和参考坐标系的夹角。所述的平移和旋转是指:根据两幅图像的质心的位置差作为平移分量,然后计算各自的主轴方向并根据其差 值将两幅图像旋转对齐,得到旋转分量,实现图像全局配准。实验证明,主轴质心法具有自动、快速、易实现等特点。
第二步:对全局配准图像和参考图像建立三角形网格
首先,将步骤1)中的全局配准图像和参考图像分别进行高斯平滑滤波,消除图像成像过程中生成的噪声;然后分别对平滑处理后的两幅待配准图像利用Harris算子提取全局配准图像和参考图像的角点,其次,利用冒泡算法对所有的角点按照横坐标由小到大进行排序,取其横坐标的最小值xmin和最大值xmax,并记录下来。同理对角点按照纵坐标由小到大进行排序,取其纵坐标的最小值ymin和最大值ymax,记录下来,由横纵坐标的最小和最大值,形成四个二维空间的点ABCD,它们坐标分别为A点(xmin,ymin),B点(xmax,ymin),C点(xmax,ymax),D点(xmin,ymax),连接四点形成四边形ABCD。再次,将四边形ABCD的两个邻边AB(xmax-xmin)、BC(ymax-ymin,)进行等距离的平分,将形成四边形网格。然后,分别正对角线的连接四边形的对角线,将形成三角形网格。最后,将三角形网格进行编号分别为1、2、......、M,M为大于1的整数。
其角点的提取过程如下:
2.1)将所处理的矩形区域窗w向任意方向移动位移(x,y),得到对应灰度值的改变量为:
其中:E(u,v)|(x,y)表示在点(x,y)处移动一个(u,v)小窗口所发生的灰度值的变化情况,I(u,v)表示图像像素点(x,y)的灰度值,(u,v)表示图像的移动变量,w(x,y)为高斯平滑因子。I(x+u,y+v)-I(x,y)表示图像灰度差值, 表示图像横坐标的梯度值, 表示图像纵坐标的梯度值, 为高斯窗口,对图像窗口进行高斯平滑,目的是提高抗噪能力。
所述的卷积是指:采用水平差分算子 对图像中的每个点可以用水平算子进行卷积,以及采用垂直差分算子 对图像中的每个点可以用垂直算子进行卷积。
2.3)对矩阵M的四个元素进行高斯平滑滤波,得到更新后的矩阵M, 以计算图像中每个像素的对应角点量R,
2.4)对图像进行分块得到若干个m*n大小的图像块,对图像块中的角点量R值进行由大到小的排序,经阈值筛选后得到图像的角点。
所述的分块是指:按照长度和宽度方向正整数方式进行图像划分。
所述的阈值筛选是指:取图像块中所有角点量的中值或平均值为阈值,当该角点量大于阈值且为某邻域的局部极大值时,该角点为图像的角点。
第三步:对全局配准图像和参考图像建立数学模型
本发明将非刚性医学图像的形变过程看作二维弹性薄板的伸缩过程,当作用于弹性体的内力与外力相等时,物体停止形变。因此我们将非刚性医学图像的形变过程看作有一个外力迫使它发生形变的过程,当作用于物体的外力为零时,图像形变结束。具体步骤如下:
假设I1(x,y)表示为浮动图像的灰度函数,I2(x,y)表示参考图像的灰度函数,[x,y]∈R2表示像素空间坐标,h表示从浮动图像到目标图像的空间变换,则任意一点的空间变换可以由如下的映射关系来表示: 其中u1(x,y),u2(x,y)均为关于坐标(x,y)的函数,记((x,y)-u(x,y))表示形变后的源图像,u(x,y)为形变位移场。配准的目的是期望形变的源图像和目标图像之间的相似度尽可能达到最大。本文中相似度用弹性***总势能函数来表示,如公式(5)
其中,σt表示弹性体的应力,ε表示弹性体的应变,∫ΩσtεdΩ表示弹性***的形变势能,Ω表示二维图像区域,∫ΩFiudΩ表示弹性***的外力势能,Fi表示作用于非刚性医学图像的外力,u表示图像的形变位移。
当两幅待配准图像的对应解剖点在空间位置上达到一致时,则***的总势能达到最小。
其次,我们将作用于待配准图像的外力看作是浮动图像和参考图像的灰度差值,如公式6所示表示两幅图像的灰度差值。
Fi=I1((x,y)+u(x,y))-I2(x,y)..........................................(6)
则得到:E总=∫ΩσtεdΩ+∫Ω(I1((x,y)+u(x,y))-I2(x,y))udΩ......(7)
当作用于***外力为零时***的总势能达最小,图像形变结束,即图像实现局部精确配准。
第四步:对全局配准图像和参考图像建立有限元方程
在对步骤3)进行求解的过程中。由于对整幅图像进行建模,使得方程求解过程比较复杂,所以将浮动图像和参考图像进行离散化,将图像感兴趣的区域划分成规则形的三角形网格,将整幅图像的形变看作为三角形单元形变的总和。如步骤2)所述。将步骤2)全局配准图像和参考图像这两幅待配准图像网格单元建立有限元方程,首先,对两幅待配准图像的网格单元进行匹配,将匹配后的网格单元,进行编号1、2、......、N,N为大于1的整数:然后图像的总势能看作为三角形势能的总和,建立有限元方程。
根据有限元分析中的变分原理和加权余量法将连续的求解图像区域离散为一组有限且按一定方式相互联结在一起单元的组合体。在每个单元内选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,通过插值函数计算出各个单元内场函数的近似解,从而实现两幅图像的配准。
如步骤3)所述当***的总势能达到最小时,非刚性医学图像的配准将结束。根据第二步我们可以对非刚性医学图像建立三角形网格,这样***的总位等于有限三角单元的势能之和,如公式8所示:
其中,m表示三角形单元的数目,i表示第几个三角形如图2所示。假设非刚性医学图像三角形单元区域的形变u(x,y)量看作较小且光滑,形变后的灰度信息可以表示为如公式9所 示:
由此可知非刚性医学图像总势能如公式10所示:
第五步:求解有限元方程
如步骤2)所述我们将非刚性医学图像划分为不同的三角形单元,形成三角形网格。
三角形单元可以由三角形单元线性积分近似,则三角形单元内任意点位移可以由三角形单元的三个节点位移和三角形单元形函数加权和得到如公式11所示:
其中ui表示三角形单元内的任意点的位移, 表示三角形单元的形函数,S为三角形单元的面积,l,m,n为三角形单元的三节点其坐标为l(xl,yl),m(xm,ym),n(xn,yn),则:ael=xmyn-xnym,bel=ym-yn,cel=xn-xm。公式6中的其余系数可由循环置换角标l,m,n的坐标而得到。
如步骤4)所述当***的总势能达到最小时,非刚性医学图像的配准将结束。则有: 将非刚性医学图像配准问题看作是
求解在各个单元节点的最佳位移,依次迭代可求的整个非刚性医学图像的形变,实现非刚性医学图像的配准。
实验结果分析
表1.1全局配准参数
质心X轴的坐标 | 222.485 |
质心Y轴的坐标 | 254.008 |
水平偏移量 | -0.859382 |
垂直偏移量 | 0.861449 |
角度偏移量 | 0.0120708 |
[0075] 由表1.1可知图像全局配准是水平、垂直偏移量、角度偏移量。
表1.2几种配准算法配准结果对比表
由表1.2可知本发明基于有限元模型的非刚性医学图像配准算法在配准的精确度有所提高,克服了以往配准算法精确度不高的缺点。
Claims (1)
1.一种基于有限元模型非刚性医学图像的配准方法,其特征在于:包含如下步骤:
1)全局配准:首先通过图像的一阶矩寻找浮动图像和参考图像的质心;再次,通过浮动图像和参考图像的二阶中心矩寻找所述浮动图像和参考图像的主轴和坐标系的夹角;最后,通过平移和旋转使得浮动图像和参考图像的质心和主轴对齐,得到旋转分量,实现浮动图像和参考图像的全局配准;
2)构建浮动图像和参考图像的网格:首先,将步骤1)中的全局配准图像和参考图像分别进行高斯平滑滤波,消除图像成像过程中生成的噪声;然后分别对平滑处理后的两幅待配准图像利用Harris算子提取图像的角点,对其角点进行排序,取角点的横坐标中最小值、最大值及纵坐标中最小值、最大值并记录下来,这样就形成四个点,连结四个点形成一个四边形,使得四边形包含所有的角点;等分四边形,使得四边形的两组对边平分成等距离,形成四边形网格,将等分后的四边形按照相同的方向连接对角线形成三角形网格,对三角形网格进行编号1、2、……、M,M为大于1的整数;
3)建立数学模型:将全局配准图像和参考图像的形变过程看作弹性的伸缩过程,利用线性弹性模型中的最小势能约束每个匹配后的形变,外力为全局配准图像和参考图像两幅图像的灰度差值,当外力为零时有限元单元的势能到最小,图像完成精确配准;
4)建立有限元方程:在对步骤3)进行求解的过程中,将全局配准图像和参考图像进行离散化,将图像感兴趣的区域划分成规则形的三角形网格,将整幅图像的形变看作为三角形单元形变的总和,对步骤2)中全局配准图像和参考图像这两幅待配准图像网格单元建立有限元方程,首先,对两幅待配准图像的网格单元进行匹配,将匹配后的网格单元,进行编号1、2、……、N,N为大于1的整数,然后图像的总势能看作为三角形势能的总和,建立有限元方程,
5)求解方程:求解步骤4)中的有限元方程,在每个单元内选择至少一个节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,通过插值函数计算出各个单元内场函数的近似解;利用三角形单元的灰度差值的平方和作为外力,当总势能最小时,外力为零,来求解方程,得到全局配准图像的偏移量,实现图像的局部精确配准。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110071997 CN102169578B (zh) | 2011-03-16 | 2011-03-16 | 一种基于有限元模型的非刚性医学图像配准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110071997 CN102169578B (zh) | 2011-03-16 | 2011-03-16 | 一种基于有限元模型的非刚性医学图像配准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102169578A CN102169578A (zh) | 2011-08-31 |
CN102169578B true CN102169578B (zh) | 2013-02-13 |
Family
ID=44490730
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110071997 Expired - Fee Related CN102169578B (zh) | 2011-03-16 | 2011-03-16 | 一种基于有限元模型的非刚性医学图像配准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102169578B (zh) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103065322B (zh) * | 2013-01-10 | 2015-03-25 | 合肥超安医疗科技有限公司 | 一种基于双x射线成像的2d和3d医学图像配准方法 |
CN103077522B (zh) * | 2013-01-22 | 2015-06-10 | 南方医科大学 | 一种改进的局部形变明显图像的配准方法 |
CN104680486A (zh) * | 2013-11-29 | 2015-06-03 | 上海联影医疗科技有限公司 | 一种非刚体的配准方法 |
KR102233966B1 (ko) * | 2014-05-12 | 2021-03-31 | 삼성전자주식회사 | 의료 영상 정합 방법 및 그 장치 |
CN104318568B (zh) * | 2014-10-24 | 2017-07-28 | 武汉华目信息技术有限责任公司 | 一种图像配准的方法和*** |
CN106780575A (zh) * | 2016-11-23 | 2017-05-31 | 山东大学 | 基于图像特征和Demons的非刚性图像配准方法 |
CN107038716B (zh) * | 2017-03-29 | 2019-06-04 | 大连外国语大学 | 一种基于椭球模型的非刚性医学图像抗旋转配准方法 |
CN107977465A (zh) * | 2017-12-21 | 2018-05-01 | 福建省晨曦信息科技股份有限公司 | 一种工程造价显示方法及终端 |
CN108399412A (zh) * | 2018-03-07 | 2018-08-14 | 鞍钢集团矿业有限公司 | 一种x型角点亚像素提取方法 |
CN108492306A (zh) * | 2018-03-07 | 2018-09-04 | 鞍钢集团矿业有限公司 | 一种基于图像轮廓的x型角点提取方法 |
CN109461140A (zh) * | 2018-09-29 | 2019-03-12 | 沈阳东软医疗***有限公司 | 图像处理方法及装置、设备和存储介质 |
CN110021068B (zh) * | 2019-02-14 | 2020-10-23 | 清华大学 | 一种基于有限元的三维医学图像重构方法 |
CN110197503B (zh) * | 2019-05-14 | 2021-07-30 | 北方夜视技术股份有限公司 | 基于增强型仿射变换的非刚性点集配准方法 |
CN112150404B (zh) * | 2019-06-27 | 2023-10-27 | 上海交通大学 | 基于联合显著图的全局到局部非刚性图像配准方法及装置 |
CN112863641A (zh) * | 2019-11-12 | 2021-05-28 | 西安大医集团股份有限公司 | 放射治疗***及其放射源的偏移量确定方法、装置 |
CN113779823B (zh) * | 2021-08-20 | 2024-01-02 | 桂林电子科技大学 | 基于微观组织图像三角形网格划分的有限元建模方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1426021A (zh) * | 2002-12-19 | 2003-06-25 | 上海交通大学 | 遥感图像的非线性配准方法 |
CN101238993A (zh) * | 2008-02-01 | 2008-08-13 | 哈尔滨工业大学 | 基于整数提升小波多分辨分析的医学超声图像配准方法 |
CN101493936A (zh) * | 2008-05-30 | 2009-07-29 | 内蒙古科技大学 | 一种基于图像边缘的多分辨非刚性头部医学图像配准方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6909794B2 (en) * | 2000-11-22 | 2005-06-21 | R2 Technology, Inc. | Automated registration of 3-D medical scans of similar anatomical structures |
-
2011
- 2011-03-16 CN CN 201110071997 patent/CN102169578B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1426021A (zh) * | 2002-12-19 | 2003-06-25 | 上海交通大学 | 遥感图像的非线性配准方法 |
CN101238993A (zh) * | 2008-02-01 | 2008-08-13 | 哈尔滨工业大学 | 基于整数提升小波多分辨分析的医学超声图像配准方法 |
CN101493936A (zh) * | 2008-05-30 | 2009-07-29 | 内蒙古科技大学 | 一种基于图像边缘的多分辨非刚性头部医学图像配准方法 |
Non-Patent Citations (5)
Title |
---|
Yurui Xie, et al..A New Non-rigid Image Registration Algorithm Using the Finite-Element Method.《2010 Second International Workshop on Education Technology and Computer Science(ETCS)》.2010,第1卷206-210. * |
吕晓琪,等..一种基于角点特征的医学图像配准方法.《内蒙古科技大学学报》.2010,第29卷(第1期),49-52. * |
吕晓琪,等。.一种基于角点特征的医学图像配准方法.《内蒙古科技大学学报》.2010,第29卷(第1期),49-52. |
龚永义,等..基于改进的弹簧质子模型的医学图像配准.《计算机学报》.2008,第31卷(第7期),1224-1233. * |
龚永义,等。.基于改进的弹簧质子模型的医学图像配准.《计算机学报》.2008,第31卷(第7期),1224-1233. |
Also Published As
Publication number | Publication date |
---|---|
CN102169578A (zh) | 2011-08-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102169578B (zh) | 一种基于有限元模型的非刚性医学图像配准方法 | |
CN105427308B (zh) | 一种稀疏和稠密特征匹配结合的图像配准方法 | |
CN102722890B (zh) | 基于光流场模型的非刚性心脏图像分级配准方法 | |
Sorzano et al. | Elastic registration of biological images using vector-spline regularization | |
Shi et al. | Temporal sparse free-form deformations | |
CN105787901A (zh) | 一种用于太阳高分辨图像序列中相邻两帧间的多尺度速度场测量方法 | |
CN108389226A (zh) | 一种基于卷积神经网络和双目视差的无监督深度预测方法 | |
US8218909B2 (en) | System and method for geodesic image matching using edge points interpolation | |
CN110363802B (zh) | 基于自动分割和骨盆对齐的***图像配准***及方法 | |
CN106204561A (zh) | 基于混合模型的***多模态图像非刚性配准方法 | |
CN101814192A (zh) | 真实感3d人脸重建的方法 | |
CN104021547A (zh) | 肺部 ct 的三维配准方法 | |
Agudo et al. | Real-time 3D reconstruction of non-rigid shapes with a single moving camera | |
CN113012208B (zh) | 多视角遥感图像配准方法及*** | |
CN105488541A (zh) | 增强现实***中基于机器学习的自然特征点识别方法 | |
Malti et al. | A linear least-squares solution to elastic shape-from-template | |
CN113570658A (zh) | 基于深度卷积网络的单目视频深度估计方法 | |
CN113516693A (zh) | 一种快速通用的图像配准方法 | |
CN106952292A (zh) | 基于6自由度场景流聚类的3d运动目标检测方法 | |
Agudo et al. | FEM models to code non-rigid EKF monocular SLAM | |
Cofaru et al. | Evaluation of digital image correlation techniques using realistic ground truth speckle images | |
Cai et al. | A target-free video structural motion estimation method based on multi-path optimization | |
CN104091364A (zh) | 单幅图像超分辨重建方法 | |
CN100446034C (zh) | 基于有限采样全局优化的图像弹性配准方法 | |
CN111696167A (zh) | 自范例学习引导的单张影像超分辨率重构方法 |
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 | ||
C53 | Correction of patent of invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: Lv Xiaoqi Inventor after: Ma Hongli Inventor after: Yu Dahua Inventor after: Zhang Baohua Inventor before: Lv Xiaoqi Inventor before: Ma Hongli Inventor before: Zhang Baohua |
|
COR | Change of bibliographic data |
Free format text: CORRECT: INVENTOR; FROM: LV XIAOQI MA HONGLI ZHANG BAOHUA TO: LV XIAOQI MA HONGLI YU DAHUA ZHANG BAOHUA |
|
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: 20130213 |