CN105160693B - 一种用于医学图像的曲率测量方法和装置 - Google Patents
一种用于医学图像的曲率测量方法和装置 Download PDFInfo
- Publication number
- CN105160693B CN105160693B CN201510631492.9A CN201510631492A CN105160693B CN 105160693 B CN105160693 B CN 105160693B CN 201510631492 A CN201510631492 A CN 201510631492A CN 105160693 B CN105160693 B CN 105160693B
- Authority
- CN
- China
- Prior art keywords
- point
- interval
- obtaining
- derivative
- abscissa
- 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.)
- Active
Links
- 238000000691 measurement method Methods 0.000 title abstract description 5
- 230000014509 gene expression Effects 0.000 claims abstract description 97
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 37
- 238000000034 method Methods 0.000 claims abstract description 36
- 238000005452 bending Methods 0.000 claims description 52
- 238000009795 derivation Methods 0.000 claims description 14
- 230000010354 integration Effects 0.000 claims description 10
- 230000001575 pathological effect Effects 0.000 abstract description 13
- 230000006870 function Effects 0.000 description 56
- 238000005259 measurement Methods 0.000 description 7
- 238000002604 ultrasonography Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000000747 cardiac effect Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000008602 contraction Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Quality & Reliability (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Geometry (AREA)
- Image Analysis (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明实施例提供了一种用于医学图像的曲率测量方法和装置，所述方法包括：获得医学图像中待测点所在的轨迹；将所述轨迹划分为至少两个区间；根据三弯矩算法获得Mi和Mi+1的值，Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数；根据Mi和Mi+1的值分别获得S′(x)和S″(x)的表达式；根据S′(x)和S″(x)的表达式分别获得所述轨迹在所述待测点的一阶导数和二阶导数；根据所述一阶导数和所述二阶导数获得所述待测点的曲率。可见本发明能够更加精确地测量出医学图像中所述待测点的曲率，从而更好地反映病理特征。
Description
技术领域
本发明涉及图像处理领域，尤其是涉及一种用于医学图像的曲率测量方法和装置。
背景技术
医学图像指的是用作医疗用途的图像，例如，超声图像以及核磁共振图像等都是常见的医学图像。对于一些医学图像来说，该医学图像的曲率能够反映出一定的病理特征。例如，心脏超声图像的曲率，能够反映心脏的舒张和收缩情况。
目前在测量医学图像的曲率时，通常采用近似理想圆的方式。具体的测量过程是：为了测量图1中形状A上一点的曲率，以形状A的几何中心为圆心得到理想圆B，其中形状A与理想圆B的面积相同。在得到理想圆B后，将理想圆B的半径rideal作为该点的曲率半径，根据曲率半径和曲率的倒数关系获得该点的近似曲率。
显然，通过这种方式测量出的形状A上每一点的近似曲率都相同，测量精度较低，在很多情况下并不能较好地反映出医学图像中的病理特征。因此，如何能够更加精确地测量医学图像中的曲率，从而更好地反映病理特征，是目前亟待解决的问题。
发明内容
本发明解决的技术问题在于提供一种用于医学图像的曲率测量方法和装置，以实现能够更加精确地测量出医学图像中的曲率，从而更好地反映病理特征。
为此，本发明解决技术问题的技术方案是：
本发明提供了一种用于医学图像的曲率测量方法，包括：
获得医学图像中待测点所在的轨迹；
将所述轨迹划分为至少两个区间；
根据三弯矩算法获得Mi和Mi+1的值，Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点；
根据Mi和Mi+1的值分别获得S′(x)和S″(x)的表达式，S′(x)表示对S(x)一阶求导后的函数，S″(x)表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数；
根据S′(x)和S″(x)的表达式分别获得所述轨迹在所述待测点的一阶导数和二阶导数；
根据所述一阶导数和所述二阶导数获得所述待测点的曲率。
可选的，所述至少两个区间共包括n+1个区间点；
所述根据三弯矩算法获得Mi和Mi+1的值，包括：
获得三弯矩方程组，所述三弯矩方程组的未知数为所述轨迹在所述n+1个区间点中各个区间点的二阶导数；
对所述三弯矩方程组求解得到Mi和Mi+1的值。
可选的，所述n+1个区间点具体为{Ek|0≤k≤n}，其中区间点E0和区间点En分别为所述轨迹的起始点和终止点；
所述获得三弯矩方程组，包括：
获得S1(x)的表达式，S1(x)为横坐标x位于区间[xj,xj+1]时的所述轨迹的三次样条插值函数，S1(x)的表达式中包括未知数Mj和Mj+1，1≤j≤n-1，xj为区间点Ej的横坐标，xj+1为区间点Ej+1的横坐标，区间点Ej和区间点Ej+1为相邻的区间点，Mj为所述轨迹在区间点Ej的二阶导数，Mj+1为所述轨迹在区间点Ej+1的二阶导数；
获得S2(x)的表达式，S2(x)为横坐标x位于区间[xj-1,xj]时的所述轨迹的三次样条插值函数，S2(x)的表达式中包括未知数Mj-1和Mj，xj-1为区间点Ej-1的横坐标，区间点Ej-1和区间点Ej为相邻的区间点，Mj-1为所述轨迹在区间点Ej-1的二阶导数；
根据S1(x)和S2(x)在横坐标x＝xj时的关系，获得n-1个线性方程；
获得横坐标x＝x0和横坐标x＝xn时的两个边界方程，x0为区间点E0的横坐标，xn为区间点En的横坐标；
将所述n-1个线性方程和所述两个边界方程，作为所述三弯矩方程组。
可选的，所述n+1个区间点具体为{Ek|0≤k≤n}，所述三弯矩方程组具体为：
其中，λ0＝1,μn＝1， hj＝xj+1-xj，hj-1＝xj-xj-1，1≤j≤n-1，xj-1为区间点Ej-1的横坐标，yj-1为区间点Ej-1的纵坐标，xj为区间点Ej的横坐标，yj为区间点Ej的纵坐标，xj+1为区间点Ej+1的横坐标，yj+1为区间点Ej+1的纵坐标，Mk为所述轨迹在区间点Ek的二阶导数，0≤k≤n。
可选的，所述根据Mi和Mi+1的值分别获得S′(x)和S″(x)的表达式，包括：
根据Mi和Mi+1的值获得S″(x)的表达式；
对S″(x)两次积分获得S(x)的表达式，S(x)的表达式包括积分常数w1和w2；
根据S(xi)＝yi和S(xi+1)＝yi+1获得积分常数w1和w2的值，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标；
对S(x)的表达式一阶求导获得S′(x)的表达式。
可选的，，S′(x)和S″(x)的表达式分别为：
其中，hi＝xi+1-xi，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标。
可选的，所述根据所述一阶导数和所述二阶导数获得所述待测点的曲率，包括：
根据获得所述待测点的曲率K；
其中，S′(xt)为所述一阶导数，S″(xt)为所述二阶导数。
可选的，获得所述轨迹在所述待测点的一阶导数和二阶导数之前，所述方法还包括：获得用户在所述轨迹上选中的点，将所述选中的点作为所述待测点；
获得所述待测点的曲率之后，所述方法还包括：向用户显示所述待测点的曲率。
本发明提供了一种用于医学图像的曲率测量装置，包括：
第一获得单元，用于获得医学图像中待测点所在的轨迹；
划分单元，用于将所述轨迹划分为至少两个区间；
第二获得单元，用于根据三弯矩算法获得Mi和Mi+1的值，Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点；
第三获得单元，用于根据Mi和Mi+1的值分别获得S′(x)和S″(x)的表达式，S′(x)表示对S(x)一阶求导后的函数，S″(x)表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数；
第四获得单元，用于根据S′(x)和S″(x)的表达式分别获得所述轨迹在所述待测点的一阶导数和二阶导数；
第五获得单元，用于根据所述一阶导数和所述二阶导数获得所述待测点的曲率。
可选的，所述至少两个区间共包括n+1个区间点；
所述第二获得单元包括：
第一获得子单元，用于获得三弯矩方程组，所述三弯矩方程组的未知数为所述轨迹在所述n+1个区间点中各个区间点的二阶导数；
第二获得子单元，用于对所述三弯矩方程组求解得到Mi和Mi+1的值。
可选的，所述n+1个区间点具体为{Ek|0≤k≤n}，其中区间点E0和区间点En分别为所述轨迹的起始点和终止点；
所述第一获得子单元具体用于：
获得S1(x)的表达式，S1(x)为横坐标x位于区间[xj,xj+1]时的所述轨迹的三次样条插值函数，S1(x)的表达式中包括未知数Mj和Mj+1，1≤j≤n-1，xj为区间点Ej的横坐标，xj+1为区间点Ej+1的横坐标，区间点Ej和区间点Ej+1为相邻的区间点，Mj为所述轨迹在区间点Ej的二阶导数，Mj+1为所述轨迹在区间点Ej+1的二阶导数；
获得S2(x)的表达式，S2(x)为横坐标x位于区间[xj-1,xj]时的所述轨迹的三次样条插值函数，S2(x)的表达式中包括未知数Mj-1和Mj，xj-1为区间点Ej-1的横坐标，区间点Ej-1和区间点Ej为相邻的区间点，Mj-1为所述轨迹在区间点Ej-1的二阶导数；
根据S1(x)和S2(x)在横坐标x＝xj时的关系，获得n-1个线性方程；
以及，获得横坐标x＝x0和横坐标x＝xn时的两个边界方程，将所述n-1个线性方程和所述两个边界方程，作为所述三弯矩方程组，x0为区间点E0的横坐标，xn为区间点En的横坐标。
可选的，所述n+1个区间点具体为{Ek|0≤k≤n}，所述三弯矩方程组具体为：
其中，λ0＝1,μn＝1， hj＝xj+1-xj，hj-1＝xj-xj-1，1≤j≤n-1，xj-1为区间点Ej-1的横坐标，yj-1为区间点Ej-1的纵坐标，xj为区间点Ej的横坐标，yj为区间点Ej的纵坐标，xj+1为区间点Ej+1的横坐标，yj+1为区间点Ej+1的纵坐标，Mk为所述轨迹在区间点Ek的二阶导数，0≤k≤n。
可选的，所述第三获得单元具体用于：
根据Mi和Mi+1的值获得S″(x)的表达式；
对S″(x)两次积分获得S(x)的表达式，S(x)的表达式包括积分常数w1和w2；
根据S(xi)＝yi和S(xi+1)＝yi+1获得积分常数w1和w2的值，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标；
以及，对S(x)的表达式一阶求导获得S′(x)的表达式。
可选的，S′(x)和S″(x)的表达式分别为：
其中，hi＝xi+1-xi，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标。
可选的，所述第五获得单元具体用于，根据获得所述待测点的曲率K；
其中，S′(xt)为所述一阶导数，S″(xt)为所述二阶导数。
可选的，还包括：
第六获得单元，用于在所述第四获得单元获得所述轨迹在所述待测点的一阶导数和二阶导数之前，获得用户在所述轨迹上选中的点，将所述选中的点作为所述待测点；
显示单元，用于在所述第五获得单元获得所述待测点的曲率之后，向用户显示所述待测点的曲率。
通过上述技术方案可知，在本发明实施例中，不再采用近似理想圆的方式测量近似曲率，而是利用三弯矩算法获得待测点的一阶导数和二阶导数，从而能够根据所述一阶导数和所述二阶导数，更加精确地测量出所述待测点的曲率。具体的测量过程包括：获得所述待测点所在的轨迹，将所述轨迹划分为至少两个区间，从而能够根据三弯矩算法获得Mi和Mi+1的值，其中Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点。之后根据Mi和Mi+1的值获得S′(x)和S″(x)的表达式，其中S′(x)表示对S(x)一阶求导后的函数，S″(x)表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数，并根据S′(x)和S″(x)的表达式分别获得所述轨迹在所述待测点的一阶导数和二阶导数，从而能够根据所述一阶导数和所述二阶导数，更加精确地测量出所述待测点的曲率。可见本发明通过这种方式测量的不再是近似曲率，能够更加精确地测量出医学图像中所述待测点的曲率，从而更好地反映病理特征。
附图说明
为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，还可以根据这些附图获得其它的附图。
图1为通过近似理想圆的方式测量近似曲率时的示意图；
图2为本发明实施例提供的一种方法实施例的流程示意图；
图3为一种心脏超声图像的示意图；
图4为本发明实施例提供的三弯矩方程组的一种获得方式的流程示意图；
图5为本发明实施例提供的步骤S204的一种实现方式的流程示意图；
图6为本发明实施例提供的一种装置实施例的结构示意图。
具体实施方式
目前在测量医学图像的曲率时，通常采用近似理想圆的方式。这种测量方式只能测量出近似曲率，而近似曲率在很多情况下并不能较好地反映出医学图像中的病理特征。因此，如何能够更加精确地测量出医学图像中一个点的曲率，更好地反映病理特征，是目前亟待解决的问题。
为了能够计算出医学图像中一个待测点的曲率，需要拟合出所述待测点所在的轨迹的表达式，之后根据拟合出的所述轨迹的表达式获得所述待测点的曲率。在拟合所述轨迹的表达式时，可以采用三次样条插值算法。下面说明通过三次样条插值算法拟合所述轨迹的表达式的一种常规思路。
根据三次样条插值算法的要求，将所述轨迹划分为n个区间，n≥2。所述n个区间共有n+1个区间点，其中，沿着所述轨迹将所述n+1个区间点顺序编号为E0，E1，……，En-1，En。区间点E0和区间点En分别为所述轨迹的起始点和终止点，区间点Ek的横坐标为xk，区间点Ek的纵坐标为yk，0≤k≤n，且a＝x0<x1<...<xn＝b。在本发明中，区间点指的是区间的边界点。
当所述待测点的横坐标x位于区间[a,b]时，如果函数S(x)满足条件1和条件2，则称S(x)为所述轨迹在区间[a,b]上的三次样条插值函数。
条件1：S(x)具有连续的一阶导数和连续的二阶导数。
条件2：S(x)在每一个区间[xm,xm+1]上为三次多项式。
其中，0≤m≤n-1，区间点Em的横坐标为xm，纵坐标为ym。
根据条件2，可以令：S(x)＝S1(x)＝am+bmx+cmx2+dmx3 (公式1)
其中，x∈[xm,xm+1]。
公式1中共包括4n个系数am、bm、cm和dm，因此为了获得该4n个系数的值，共需要4n个独立条件。下面提供三类独立条件。
第一类独立条件为插值条件：
根据n+1个区间点中每个区间点的横坐标和纵坐标的值，有：S(xm)＝ym，其中m从0依次取值到n-1，因此共提供了n+1个独立条件。
第二类独立条件为连续性条件：
令m＝j，公式1可表示为：
S(x)＝S1(x)＝aj+bjx+cjx2+djx3
其中，x∈[xj,xj+1]，1≤j≤n-1。
根据条件2，令S(x)＝S2(x)＝aj-1+bj-1x+cj-1x2+dj-1x3
其中，x∈[xj-1,xj]。
根据S1(x)和S2(x)在x＝xj时的连续性，提供了3n-3个独立条件：
其中，S1'(xj)表示对S1(xj)一阶求导后的函数，S2'(xj)表示对S2(xj)一阶求导后的函数，S1"(xj)表示对S1(xj)二阶求导后的函数，S2"(xj)表示对S2(xj)二阶求导后的函数。
所述第一类独立条件和所述第二类独立条件共提供了4n-2个独立条件，为获得所述4n个系数的值，还需要两个边界条件作为独立条件。例如可以使用固定边界条件，即令起始点和终止点的一阶导数为零：因此有S'(x0)＝S'(xn)＝0
因此，根据所述第一类独立条件、所述第二类独立条件和所述两个边界条件，能够计算出4n个系数am、bm、cm和dm的值，将计算出的4n个系数的值代入到公式1中，即可得到所述轨迹在每一个区间[xm,xm+1]的表达式，从而根据所述待测点所在的区间的表达式，获得所述待测点的曲率。
然而在采用上述方式获得所述轨迹在每一个区间[xm,xm+1]的表达式时，由于要获得4n个系数的值，因此需要对4n个独立方程进行求解，计算量比较大。
为了解决上述问题，本发明提供了一种用于医学图像的曲率测量方法和装置，不仅能够更加精确地计算出医学图像中待测点的曲率，从而更好地反映病理特征，而且具有计算量更小的特点。本发明的主要思想包括但不局限于：利用三弯矩算法获得待测点的一阶导数和二阶导数，从而根据所述一阶导数和所述二阶导数，直接测量出所述待测点的曲率。
请参阅图2，本发明提供了用于医学图像的曲率测量方法的一种方法实施例，本实施例的所述方法包括：
S201：获得医学图像中待测点所在的轨迹。
所述待测点所在的轨迹一般为所述医学图像中的边界轨迹，可以通过人工方式标出，也可以由设备根据所述医学图像的灰度值自动获得。例如在图3所示的心脏超声图像中，可以通过人工方式标记出一段边界轨迹S，作为所述待测点所在的轨迹。
其中，可以先确定出所述待测点，再根据所述待测点获得所述轨迹；也可以先获得所述轨迹，再从所述轨迹上确定出所述待测点，例如在图3中，在获得所述边界轨迹S之后，可以将用户在所述边界轨迹S上选中的一点作为所述待测点。其中，用户可以通过单击、双击、或者光标停留时间大于预设时间等操作选中所述待测点。
S202：将所述轨迹划分为至少两个区间。
由于本发明采用的是三次样条插值算法，因此根据该算法的要求，需要将所述轨迹划分为n个区间，n大于或等于2。其中，可以将所述轨迹划分为n个等间距的区间。例如在图3所示的心脏超声图像中，将边界轨迹S划分成10个等间距的区间。该10个区间共包括11个区间点，沿着所述轨迹将所述n+1个区间点顺序编号为E0，E1，……，E9，E10，区间点E0和区间点E10分别为边界轨迹S的起始点和终止点。
S203：根据三弯矩算法获得Mi和Mi+1的值，Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点。
在本发明实施例中，采用三次样条插值算法获得所述轨迹的表达式，即所述轨迹的三次样条插值函数。其中，根据三次样条插值算法的要求，将所述轨迹划分为n个区间，n≥2。所述n个区间共有n+1个区间点{Ek|0≤k≤n},Mk为所述轨迹在区间点Ek的二阶导数。所述轨迹的三次样条插值函数可以有多种表达式，将用每个区间点Ek的二阶导数值Mk表示所述轨迹的三次样条插值函数的算法称为三弯矩算法。通过三弯矩算法能够构造包括Mk的三弯矩方程组，从而获得Mk的值，根据所述待测点所在的区间从Mk的值中获得Mi和Mi+1的值。
S204：根据Mi和Mi+1的值获得S′(x)和S″(x)的表达式。
其中，S′(x)的表达式以Mi和Mi+1为系数，以横坐标x为变量，表示对S(x)一阶求导后的函数。S″(x)的表达式以Mi和Mi+1为系数，以横坐标x为变量，表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数。
本发明实施例中不再通过公式1，而是通过Mi和Mi+1表示S(x)、S′(x)和S″(x)表达式，当获得Mi和Mi+1的值后，就能够根据Mi和Mi+1的值获得S(x)、S′(x)和S″(x)的表达式。
S205：根据S′(x)和S″(x)的表达式分别获得所述轨迹在所述待测点的一阶导数和二阶导数。
由于S′(x)表示对S(x)一阶求导后的函数，S″(x)表示对S(x)二阶求导后的函数，因此S′(x)在所述待测点处的函数值即为所述待测点的一阶导数，S″(x)在所述待测点处的函数值即为所述待测点的二阶导数，因此，令x＝xt时，S′(x)和S″(x)的表达式的值即分别为所述待测点的一阶导数S′(xt)和二阶导数S″(xt)。
S206：根据所述一阶导数和所述二阶导数获得所述待测点的曲率。
例如，可以根据获得所述待测点的曲率K。
其中，S′(xt)为所述轨迹在所述待测点的一阶导数，S″(xt)为所述轨迹在所述待测点的二阶导数。
通过上述技术方案可知，在本发明实施例中，不再采用近似理想圆的方式测量近似曲率，而是利用三弯矩算法获得待测点的一阶导数和二阶导数，从而能够根据所述一阶导数和所述二阶导数，更加精确地测量出所述待测点的曲率。具体的测量过程包括：获得所述待测点所在的轨迹，将所述轨迹划分为至少两个区间，从而能够根据三弯矩算法获得Mi和Mi+1的值，其中Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点。之后根据Mi和Mi+1的值获得S′(x)和S″(x)的表达式，其中S′(x)表示对S(x)一阶求导后的函数，S″(x)表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数，也就是说本发明中通过Mi和Mi+1作为系数表示S′(x)的表达式和S″(x)的表达式，当获得Mi和Mi+1的值后，就能够根据Mi和Mi+1的值获得S′(x)的表达式和S″(x)的表达式，之后根据S′(x)的表达式和S″(x)的表达式获得所述轨迹在所述待测点的一阶导数和二阶导数，从而能够根据所述一阶导数和所述二阶导数，更加精确地测量出所述待测点的曲率。可见本发明实施例通过这种方式测量的不再是近似曲率，能够更加精确地确定出医学图像中所述待测点的曲率，从而更好地反映病理特征。其中，所述待测点可以为所述轨迹上的任一点，因此本发明实施例中可以测量出所述轨迹上任一点的曲率。
同时，本发明实施例中不再通过公式1构造S(x)的常规表达式，因此无需通过对4n个方程求解以得到4n个系数am、bm、cm和dm，即无需计算出所述轨迹的函数表达式。只需根据三弯矩算法计算出Mi和Mi+1，并且根据Mi和Mi+1获得所述待测点的一阶导数和二阶导数，从而根据所述一阶导数和所述二阶导数直接计算出所述待测点的曲率。其中，由于根据三弯矩算法计算Mi和Mi+1时，一般只需对n+1个方程求解，相比于对4n个方程求解，计算量大幅度减少。尤其是n的值越大时，本发明对计算量的改善就越明显。
并且，在对4n个方程求解以得到4n个系数am、bm、cm和dm时，并不能保证一定有解，因此可能无法测量出曲率。而本发明利用三弯矩算法计算Mi和Mi+1属于一种数值解法，因此能够保证一定有解，从而一定能够测量出所述待测点的曲率。
本发明实施例所具有的精确地更高、计算量更小以及保证一定有解的特点，使得本发明实施例非常适合工程应用。
下面分别举例说明S203和S204的具体实现方式。
所述至少两个区间共包括n+1个区间点，S203可以包括：获得三弯矩方程组，所述三弯矩方程组的未知数为所述轨迹在所述n+1个区间点中各个区间点的二阶导数；对所述三弯矩方程组求解得到Mi和Mi+1的值。
下面对所述获得三弯矩方程组这一步骤进行具体说明。
如图4所示，所述获得三弯矩方程组可以包括S2031至S2035。
本实施例中的所述n+1个区间点具体为{Ek|0≤k≤n}，其中区间点E0和区间点En分别为所述轨迹的起始点和终止点。
S2031：获得S1(x)的表达式，S1(x)为横坐标x位于区间[xj,xj+1]时的所述轨迹的三次样条插值函数，S1(x)的表达式中包括未知数Mj和Mj+1，1≤j≤n-1，xj为区间点Ej的横坐标，xj+1为区间点Ej+1的横坐标，区间点Ej和区间点Ej+1为所述n+1个区间点中相邻的区间点，Mj为所述轨迹在区间点Ej的二阶导数，Mj+1为所述轨迹在区间点Ej+1的二阶导数。其中，xj＜xj+1。
下面说明S1(x)的表达式的一种获取方式。
由于S1″(x)在区间[xj,xj+1]上为一次函数，因此由两点拉格朗日插值可以表示为：
(公式2)
其中，S1″(x)表示对S1(x)二阶求导后的函数，hj＝xj+1-xj。
对公式2进行积分，可以得到：
(公式3)
其中，S1′(x)表示对S1(x)一阶求导后的函数，w1为积分常数。
对公式3积分，可以得到：
(公式4)
其中，w2为积分常数。
根据S(xj)＝yj和S(xj+1)＝yj+1，可以确定出积分常数w1和w2分别为：
(公式5)
其中，yj为区间点Ej的纵坐标，yj+1为区间点Ej+1的纵坐标。
将公式5代入公式4和公式3，可以得出S1(x)和S1′(x)的表达式：
S2032：获得S2(x)的表达式，S2(x)为横坐标x位于区间[xj-1,xj]时的所述轨迹的三次样条插值函数，S2(x)的表达式中包括未知数Mj-1和Mj，xj-1为区间点Ej-1的横坐标，区间点Ej-1和区间点Ej为相邻的区间点，Mj-1为所述轨迹在区间点Ej-1的二阶导数。其中，xj-1＜xj。
与获取S1(x)的表达式的过程类似，可以获得S2(x)和S2′(x)的表达式：
其中，S2′(x)表示对S2(x)一阶求导后的函数，yj-1为区间点Ej-1的纵坐标，hj-1＝xj-xj-1。
S2033：根据S1(x)和S2(x)在横坐标x＝xj时的关系，获得n-1个线性方程。
为了获得Mk的值，需要使用横坐标x＝xj时的导数连续条件：
S1′(xj)＝S2′(xj) (公式6)
根据公式6、以及S1′(x)和S2′(x)的表达式，可以得到：
(公式7)
对公式7整理后得到：
μjMj-1+2Mj+λjMj+1＝dj (公式8)
其中，
其中，j从1取到n-1，即可得到n-1个线性方程：
μ1M0+2M1+λ1M2＝d1
μ2M1+2M2+λ2M3＝d2
……
μn-1Mn-2+2Mn-1+λn-1Mn＝dn-1
S2034：获得横坐标x＝x0和横坐标x＝xn时的两个边界方程，x0为区间点E0的横坐标，xn为区间点En的横坐标。
所述n-1个线性方程共含有n+1个未知数Mk，因此为了确定n+1个未知数Mk的值，还需横坐标x＝x0和横坐标x＝xn时的两个边界方程。例如可以采用固定边界条件：S1′(x0)＝0和S2′(xn)＝0。
其中，根据S1′(x)的表达式有：
(公式9)
将j＝0代入公式9，并根据固定边界条件：S1′(x0)＝0可得到：
(公式10)
类似地，根据S2′(x)的表达式和固定边界条件：S2′(xn)＝0可得到：
(公式11)
S2035：将所述n-1个线性方程和所述两个边界方程，作为所述三弯矩方程组，所述三弯矩方程组的未知数为所述轨迹在所述n+1个区间点中各个区间点的二阶导数。
令λ0＝1，μn＝1，
根据所述n-1个线性方程、以及根据固定边界条件得到的公式10和公式11，可以得到三弯矩方程组：
其中，λ0＝1,μn＝1， hj＝xj+1-xj，hj-1＝xj-xj-1，1≤j≤n-1，xj-1为区间点Ej-1的横坐标，yj-1为区间点Ej-1的纵坐标，xj为区间点Ej的横坐标，yj为区间点Ej的纵坐标，xj+1为区间点Ej+1的横坐标，yj+1为区间点Ej+1的纵坐标，Mk为所述轨迹在区间点Ek的二阶导数。
对三弯矩方程组求解可以得到n+1个未知数Mk的值，从而得到Mi和Mi+1的值。
如图5所示，S204可以包括S2041至S2044。
S2041：根据Mi和Mi+1的值获得S″(x)的表达式，S″(x)表示对S(x)二阶求导后的函数。
令S(x)＝S1(x)，j＝i，根据公式2可得：
其中，所述待测点的横坐标x∈[xi,xi+1]，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，hi＝xi+1-xi。
S2042：对S″(x)两次积分获得S(x)的表达式，S(x)的表达式包括积分常数w1和w2。
对S”(x)两次积分后可获得：
S2043：根据S(xi)＝yi和S(xi+1)＝yi+1获得积分常数w1和w2的值，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标。
具体地，有：
得到积分常数w1和w2的值后，即可得到S(x)的表达式：
S2044：对S(x)的表达式一阶求导获得S′(x)的表达式。
在本发明实施例中，n+1个区间点E0至En沿着所述轨迹，按照顺序进行排布，根据三次样条插值函数的要求，需要满足n+1个区间点的横坐标满足：a＝x0<x1<...<xn＝b。因此本发明实施例中，可以在获取轨迹时即保证满足上述条件，或者，如果所述n+1个区间点的横坐标的不满足上述条件，则可以将所述轨迹划分成若干个子轨迹，并且保证每个子轨迹上的区间点满足上述条件，根据所述待测点所在的子轨迹进行曲率测量。
本发明实施例中在获得所述待测点的曲率后，还可以进一步显示所述曲率，以使得根据所述曲率确定所述医学图像能够表示的病理特征。具体地，S205中的获得所述轨迹在所述待测点的一阶导数和二阶导数之前，所述方法还包括：获得用户在所述轨迹上选中的点，将所述选中的点作为所述待测点；S206中的获得所述待测点的曲率之后，所述方法还包括：向用户显示所述待测点的曲率。
本发明实施例中，当所述待测点发生变化时，还可以实时地更新所述待测点的曲率。例如，用户可以将光标在图3所示的边界轨迹S上进行移动，当光标移动到任一点时，可以将该点作为所述待测点，根据本发明实施例测量出该点的曲率并进行显示。
对应上述方法实施例，本发明还提供了一种用于医学图像的曲率测量装置，具体通过下面的实施例进行说明。
请参阅图6，本发明提供了用于医学图像的曲率测量装置的一种装置实施例，本实施例的所述装置包括：第一获得单元601、划分单元602、第二获得单元603、第三获得单元604、第四获得单元605和第五获得单元606。
第一获得单元601，用于获得医学图像中待测点所在的轨迹。
所述待测点所在的轨迹一般为所述医学图像中的边界轨迹，可以通过人工方式标出，也可以由设备根据所述医学图像的灰度值自动获得。例如在图3所示的心脏超声图像中，可以通过人工方式标记出一段边界轨迹S，作为所述待测点所在的轨迹。
其中，可以先确定出所述待测点，再根据所述待测点获得所述轨迹；也可以先获得所述轨迹，再从所述轨迹上确定出所述待测点，例如在图3中，在获得所述边界轨迹S之后，可以将用户在所述边界轨迹S上选中的一点作为所述待测点。其中，用户可以通过单击、双击、或者光标停留时间大于预设时间等操作选中所述待测点。
划分单元602，用于将所述轨迹划分为至少两个区间。
由于本发明采用的是三次样条插值算法，因此根据该算法的要求，需要将所述轨迹划分为n个区间，n大于或等于2。其中，可以将所述轨迹划分为n个等间距的区间。例如在图3所示的心脏超声图像中，将边界轨迹S划分成10个等间距的区间。该10个区间共包括11个区间点，沿着所述轨迹将所述n+1个区间点顺序编号为E0，E1，……，E9，E10，区间点E0和区间点E10分别为边界轨迹S的起始点和终止点。
第二获得单元603，用于根据三弯矩算法获得Mi和Mi+1的值，Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点。
在本发明实施例中，采用三次样条插值算法获得所述轨迹的表达式，即所述轨迹的三次样条插值函数。其中，根据三次样条插值算法的要求，将所述轨迹划分为n个区间，n≥2。所述n个区间共有n+1个区间点{Ek|0≤k≤n},Mk为所述轨迹在区间点Ek的二阶导数。通过三弯矩算法能够构造包括Mk的三弯矩方程组，从而获得Mk的值，根据所述待测点所在的区间从Mk的值中获得Mi和Mi+1的值。
第三获得单元604，用于根据Mi和Mi+1的值获得S′(x)和S″(x)的表达式。
其中，S′(x)的表达式以Mi和Mi+1为系数，以横坐标x为变量，表示对S(x)一阶求导后的函数。S″(x)的表达式以Mi和Mi+1为系数，以横坐标x为变量，表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数。
本发明实施例中不再通过公式1，而是通过Mi和Mi+1表示S(x)、S′(x)和S″(x)表达式，当获得Mi和Mi+1的值后，就能够根据Mi和Mi+1的值获得S(x)、S′(x)和S″(x)的表达式。
第四获得单元605，用于根据S′(x)和S″(x)的表达式分别获得所述轨迹在所述待测点的一阶导数和二阶导数。
由于S′(x)表示对S(x)一阶求导后的函数，S″(x)表示对S(x)二阶求导后的函数，因此S′(x)在所述待测点处的函数值即为所述待测点的一阶导数，S″(x)在所述待测点处的函数值即为所述待测点的二阶导数，因此，令x＝xt时，S′(x)和S″(x)的表达式的值即分别为所述待测点的一阶导数S′(xt)和二阶导数S″(xt)。
第五获得单元606，用于根据所述一阶导数和所述二阶导数获得所述待测点的曲率。
例如，第五获得单元606可以具体用于根据获得所述待测点的曲率K。
其中，S′(xt)为所述轨迹在所述待测点的一阶导数，S″(xt)为所述轨迹在所述待测点的二阶导数。
通过上述技术方案可知，在本发明实施例中，不再采用近似理想圆的方式测量近似曲率，而是利用三弯矩算法获得待测点的一阶导数和二阶导数，从而能够根据所述一阶导数和所述二阶导数，更加精确地测量出所述待测点的曲率。具体地，由第一获得单元601获得所述待测点所在的轨迹，划分单元602将所述轨迹划分为至少两个区间，从而第二获得单元603根据三弯矩算法获得Mi和Mi+1的值，其中Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点。第三获得单元604根据Mi和Mi+1的值获得S′(x)和S″(x)的表达式，其中S′(x)表示对S(x)一阶求导后的函数，S″(x)表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数，也就是说本发明中通过Mi和Mi+1作为系数表示S′(x)的表达式和S″(x)的表达式，当获得Mi和Mi+1的值后，就能够根据Mi和Mi+1的值获得S′(x)的表达式和S″(x)的表达式，第四获得单元605根据S′(x)的表达式和S″(x)的表达式获得所述轨迹在所述待测点的一阶导数和二阶导数，第六获得单元606能够根据所述一阶导数和所述二阶导数，更加精确地测量出所述待测点的曲率。可见本发明实施例通过这种方式测量的不再是近似曲率，能够更加精确地确定出医学图像中所述待测点的曲率，从而更好地反映病理特征。其中，所述待测点可以为所述轨迹上的任一点，因此本发明实施例中可以测量出所述轨迹上任一点的曲率。
同时，本发明实施例中不再通过公式1构造S(x)的常规表达式，因此无需通过对4n个方程求解以得到4n个系数am、bm、cm和dm，即无需计算出所述轨迹的函数表达式。只需根据三弯矩算法计算出Mi和Mi+1，并且根据Mi和Mi+1获得所述待测点的一阶导数和二阶导数，从而根据所述一阶导数和所述二阶导数直接计算出所述待测点的曲率。其中，由于根据三弯矩算法计算Mi和Mi+1时，一般只需对n+1个方程求解，相比于对4n个方程求解，计算量大幅度减少。尤其是n的值越大时，本发明对计算量的改善就越明显。
并且，在对4n个方程求解以得到4n个系数am、bm、cm和dm时，并不能保证一定有解，因此可能无法测量出曲率。而本发明利用三弯矩算法计算Mi和Mi+1属于一种数值解法，因此能够保证一定有解，从而一定能够测量出所述待测点的曲率。
本发明实施例所具有的精确地更高、计算量更小以及保证一定有解的特点，使得本发明实施例非常适合工程应用。
下面分别举例说明第二获得单元603和第三获得单元604的具体实现方式。
所述至少两个区间共包括n+1个区间点；所述第二获得单元包括：第一获得子单元和第二获得子单元，所述第一获得子单元用于获得三弯矩方程组，所述三弯矩方程组的未知数为所述轨迹在所述n+1个区间点中各个区间点的二阶导数；所述第二获得子单元用于对所述三弯矩方程组求解得到Mi和Mi+1的值。
下面对所述第一获得子单元进行具体说明。
本实施例中的所述n+1个区间点具体为{Ek|0≤k≤n}，其中区间点E0和区间点En分别为所述轨迹的起始点和终止点。
所述第一获得子单元可以具体用于：
获得S1(x)的表达式，S1(x)为横坐标x位于区间[xj,xj+1]时的所述轨迹的三次样条插值函数，S1(x)的表达式中包括未知数Mj和Mj+1，1≤j≤n-1，xj为区间点Ej的横坐标，xj+1为区间点Ej+1的横坐标，区间点Ej和区间点Ej+1为相邻的区间点，Mj为所述轨迹在区间点Ej的二阶导数，Mj+1为所述轨迹在区间点Ej+1的二阶导数；
获得S2(x)的表达式，S2(x)为横坐标x位于区间[xj-1,xj]时的所述轨迹的三次样条插值函数，S2(x)的表达式中包括未知数Mj-1和Mj，xj-1为区间点Ej-1的横坐标，区间点Ej-1和区间点Ej为相邻的区间点，Mj-1为所述轨迹在区间点Ej-1的二阶导数；
根据S1(x)和S2(x)在横坐标x＝xj时的关系，获得n-1个线性方程；
以及，获得横坐标x＝x0和横坐标x＝xn时的两个边界方程，将所述n-1个线性方程和所述两个边界方程，作为所述三弯矩方程组，x0为区间点E0的横坐标，xn为区间点En的横坐标。
其中，所述n+1个区间点具体为{Ek|0≤k≤n}，所述三弯矩方程组可以具体为：
其中，λ0＝1,μn＝1， hj＝xj+1-xj，hj-1＝xj-xj-1，1≤j≤n-1，xj-1为区间点Ej-1的横坐标，yj-1为区间点Ej-1的纵坐标，xj为区间点Ej的横坐标，yj为区间点Ej的纵坐标，xj+1为区间点Ej+1的横坐标，yj+1为区间点Ej+1的纵坐标，Mk为所述轨迹在区间点Ek的二阶导数，0≤k≤n。
第三获得单元604可以具体用于：
根据Mi和Mi+1的值获得S″(x)的表达式；
对S″(x)两次积分获得S(x)的表达式，S(x)的表达式包括积分常数w1和w2；
根据S(xi)＝yi和S(xi+1)＝yi+1获得积分常数w1和w2的值，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标；
以及，对S(x)的表达式一阶求导获得S′(x)的表达式。
其中，S′(x)和S″(x)的表达式可以分别为：
其中，hi＝xi+1-xi。
在本发明实施例中，n+1个区间点E0至En沿着所述轨迹，按照顺序进行排布，根据三次样条插值函数的要求，需要满足n+1个区间点的横坐标满足：a＝x0<x1<...<xn＝b。因此本发明实施例中，可以使得第一获得单元601在获取轨迹时即保证满足上述条件，或者，如果所述n+1个区间点的横坐标的不满足上述条件，则可以由将所述轨迹划分成若干个子轨迹，并且保证每个子轨迹上的区间点满足上述条件，根据所述待测点所在的子轨迹进行曲率测量。
本发明实施例中在获得所述待测点的曲率后，还可以进一步显示所述曲率，以使得根据所述曲率确定所述医学图像能够表示的病理特征。具体地，所述曲率测量装置还包括：第六获得单元，用于在第四获得单元605获得所述轨迹在所述待测点的一阶导数和二阶导数之前，获得用户在所述轨迹上选中的点，将所述选中的点作为所述待测点；显示单元，用于在第五获得单元606获得所述待测点的曲率之后，向用户显示所述待测点的曲率。
本发明实施例中，当所述待测点发生变化时，还可以实时地更新所述待测点的曲率。例如，用户可以将光标在图3所示的边界轨迹S上进行移动，当光标移动到任一点时，可以将该点作为所述待测点，根据本发明实施例测量出该点的曲率并进行显示。
所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的装置和单元的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。
在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个***，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本申请各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-OnlyMemory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述，以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (16)
1.一种用于医学图像的曲率测量方法，其特征在于，包括：
获得医学图像中待测点所在的轨迹；
将所述轨迹划分为至少两个区间；
根据三弯矩算法获得Mi和Mi+1的值，Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点；
根据Mi和Mi+1的值分别获得S′(x)和S″(x)的表达式，S′(x)的表达式以Mi和Mi+1为系数，S′(x)表示对S(x)一阶求导后的函数，S″(x)的表达式以Mi和Mi+1为系数，S″(x)表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数；
根据S′(x)和S″(x)的表达式分别获得所述轨迹在所述待测点的一阶导数和二阶导数；
根据所述一阶导数和所述二阶导数获得所述待测点的曲率。
2.根据权利要求1所述的方法，其特征在于，所述至少两个区间共包括n+1个区间点；
所述根据三弯矩算法获得Mi和Mi+1的值，包括：
获得三弯矩方程组，所述三弯矩方程组的未知数为所述轨迹在所述n+1个区间点中各个区间点的二阶导数；
对所述三弯矩方程组求解得到Mi和Mi+1的值。
3.根据权利要求2所述的方法，其特征在于，所述n+1个区间点具体为{Ek|0≤k≤n}，其中区间点E0和区间点En分别为所述轨迹的起始点和终止点；
所述获得三弯矩方程组，包括：
获得S1(x)的表达式，S1(x)为横坐标x位于区间[xj,xj+1]时的所述轨迹的三次样条插值函数，S1(x)的表达式中包括未知数Mj和Mj+1，1≤j≤n-1，xj为区间点Ej的横坐标，xj+1为区间点Ej+1的横坐标，区间点Ej和区间点Ej+1为相邻的区间点，Mj为所述轨迹在区间点Ej的二阶导数，Mj+1为所述轨迹在区间点Ej+1的二阶导数；
获得S2(x)的表达式，S2(x)为横坐标x位于区间[xj-1,xj]时的所述轨迹的三次样条插值函数，S2(x)的表达式中包括未知数Mj-1和Mj，xj-1为区间点Ej-1的横坐标，区间点Ej-1和区间点Ej为相邻的区间点，Mj-1为所述轨迹在区间点Ej-1的二阶导数；
根据S1(x)和S2(x)在横坐标x＝xj时的关系，获得n-1个线性方程；
获得横坐标x＝x0和横坐标x＝xn时的两个边界方程，x0为区间点E0的横坐标，xn为区间点En的横坐标；
将所述n-1个线性方程和所述两个边界方程，作为所述三弯矩方程组。
4.根据权利要求2或3所述的方法，其特征在于，所述n+1个区间点具体为{Ek|0≤k≤n}，所述三弯矩方程组具体为：
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中， hj＝xj+1-xj，hj-1＝xj-xj-1，1≤j≤n-1，xj-1为区间点Ej-1的横坐标，yj-1为区间点Ej-1的纵坐标，xj为区间点Ej的横坐标，yj为区间点Ej的纵坐标，xj+1为区间点Ej+1的横坐标，yj+1为区间点Ej+1的纵坐标，Mk为所述轨迹在区间点Ek的二阶导数，0≤k≤n。
5.根据权利要求1所述的方法，其特征在于，所述根据Mi和Mi+1的值分别获得S′(x)和S″(x)的表达式，包括：
根据Mi和Mi+1的值获得S″(x)的表达式；
对S″(x)两次积分获得S(x)的表达式，S(x)的表达式包括积分常数w1和w2；
根据S(xi)＝yi和S(xi+1)＝yi+1获得积分常数w1和w2的值，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标；
对S(x)的表达式一阶求导获得S′(x)的表达式。
6.根据权利要求1或5所述的方法，其特征在于，S′(x)和S″(x)的表达式分别为：
<mrow>
<msup>
<mi>S</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>M</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
</mrow>
<mn>6</mn>
</mfrac>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mo>;</mo>
</mrow>
<mrow>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mfrac>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mfrac>
<msub>
<mi>M</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>;</mo>
</mrow>
其中，hi＝xi+1-xi，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标。
7.根据权利要求1所述的方法，其特征在于，所述根据所述一阶导数和所述二阶导数获得所述待测点的曲率，包括：
根据获得所述待测点的曲率K；
其中，S′(xt)为所述一阶导数，S″(xt)为所述二阶导数。
8.根据权利要求1所述的方法，其特征在于，
获得所述轨迹在所述待测点的一阶导数和二阶导数之前，所述方法还包括：获得用户在所述轨迹上选中的点，将所述选中的点作为所述待测点；
获得所述待测点的曲率之后，所述方法还包括：向用户显示所述待测点的曲率。
9.一种用于医学图像的曲率测量装置，其特征在于，包括：
第一获得单元，用于获得医学图像中待测点所在的轨迹；
划分单元，用于将所述轨迹划分为至少两个区间；
第二获得单元，用于根据三弯矩算法获得Mi和Mi+1的值，Mi为所述轨迹在第一区间点的二阶导数，Mi+1为所述轨迹在第二区间点的二阶导数，所述第一区间点和所述第二区间点为所述待测点所在区间的两个区间点；
第三获得单元，用于根据Mi和Mi+1的值分别获得S′(x)和S″(x)的表达式，S′(x)的表达式以Mi和Mi+1为系数，S′(x)表示对S(x)一阶求导后的函数，S″(x)的表达式以Mi和Mi+1为系数，S″(x)表示对S(x)二阶求导后的函数，S(x)为所述轨迹在所述待测点所在区间的三次样条插值函数；
第四获得单元，用于根据S′(x)和S″(x)的表达式分别获得所述轨迹在所述待测点的一阶导数和二阶导数；
第五获得单元，用于根据所述一阶导数和所述二阶导数获得所述待测点的曲率。
10.根据权利要求9所述的装置，其特征在于，所述至少两个区间共包括n+1个区间点；
所述第二获得单元包括：
第一获得子单元，用于获得三弯矩方程组，所述三弯矩方程组的未知数为所述轨迹在所述n+1个区间点中各个区间点的二阶导数；
第二获得子单元，用于对所述三弯矩方程组求解得到Mi和Mi+1的值。
11.根据权利要求10所述的装置，其特征在于，所述n+1个区间点具体为{Ek|0≤k≤n}，其中区间点E0和区间点En分别为所述轨迹的起始点和终止点；
所述第一获得子单元具体用于：
获得S1(x)的表达式，S1(x)为横坐标x位于区间[xj,xj+1]时的所述轨迹的三次样条插值函数，S1(x)的表达式中包括未知数Mj和Mj+1，1≤j≤n-1，xj为区间点Ej的横坐标，xj+1为区间点Ej+1的横坐标，区间点Ej和区间点Ej+1为相邻的区间点，Mj为所述轨迹在区间点Ej的二阶导数，Mj+1为所述轨迹在区间点Ej+1的二阶导数；
获得S2(x)的表达式，S2(x)为横坐标x位于区间[xj-1,xj]时的所述轨迹的三次样条插值函数，S2(x)的表达式中包括未知数Mj-1和Mj，xj-1为区间点Ej-1的横坐标，区间点Ej-1和区间点Ej为相邻的区间点，Mj-1为所述轨迹在区间点Ej-1的二阶导数；
根据S1(x)和S2(x)在横坐标x＝xj时的关系，获得n-1个线性方程；
以及，获得横坐标x＝x0和横坐标x＝xn时的两个边界方程，将所述n-1个线性方程和所述两个边界方程，作为所述三弯矩方程组，x0为区间点E0的横坐标，xn为区间点En的横坐标。
12.根据权利要求10或11所述的装置，其特征在于，所述n+1个区间点具体为{Ek|0≤k≤n}，所述三弯矩方程组具体为：
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>0</mn>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>&mu;</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
<mtd>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>&mu;</mi>
<mn>1</mn>
</msub>
</mtd>
<mtd>
<mn>2</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>M</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mn>0</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>...</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mrow>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mi>n</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中， hj＝xj+1-xj，hj-1＝xj-xj-1，1≤j≤n-1，xj-1为区间点Ej-1的横坐标，yj-1为区间点Ej-1的纵坐标，xj为区间点Ej的横坐标，yj为区间点Ej的纵坐标，xj+1为区间点Ej+1的横坐标，yj+1为区间点Ej+1的纵坐标，Mk为所述轨迹在区间点Ek的二阶导数，0≤k≤n。
13.根据权利要求9所述的装置，其特征在于，所述第三获得单元具体用于：
根据Mi和Mi+1的值获得S″(x)的表达式；
对S″(x)两次积分获得S(x)的表达式，S(x)的表达式包括积分常数w1和w2；
根据S(xi)＝yi和S(xi+1)＝yi+1获得积分常数w1和w2的值，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标；
以及，对S(x)的表达式一阶求导获得S′(x)的表达式。
14.根据权利要求9或13所述的装置，其特征在于，S′(x)和S″(x)的表达式分别为：
<mrow>
<msup>
<mi>S</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mrow>
</mfrac>
<msub>
<mi>M</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>y</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>M</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
</mrow>
<mn>6</mn>
</mfrac>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
<mo>;</mo>
</mrow>
<mrow>
<msup>
<mi>S</mi>
<mrow>
<mo>&prime;</mo>
<mo>&prime;</mo>
</mrow>
</msup>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<mi>x</mi>
</mrow>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mfrac>
<msub>
<mi>M</mi>
<mi>i</mi>
</msub>
<mo>+</mo>
<mfrac>
<mrow>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>i</mi>
</msub>
</mrow>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mfrac>
<msub>
<mi>M</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>;</mo>
</mrow>
其中，hi＝xi+1-xi，xi为所述第一区间点的横坐标，xi+1为所述第二区间点的横坐标，yi为所述第一区间点的纵坐标，yi+1为所述第二区间点的纵坐标。
15.根据权利要求9所述的装置，其特征在于，所述第五获得单元具体用于，根据获得所述待测点的曲率K；
其中，S′(xt)为所述一阶导数，S″(xt)为所述二阶导数。
16.根据权利要求9所述的装置，其特征在于，还包括：
第六获得单元，用于在所述第四获得单元获得所述轨迹在所述待测点的一阶导数和二阶导数之前，获得用户在所述轨迹上选中的点，将所述选中的点作为所述待测点；
显示单元，用于在所述第五获得单元获得所述待测点的曲率之后，向用户显示所述待测点的曲率。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510631492.9A CN105160693B (zh) | 2015-09-28 | 2015-09-28 | 一种用于医学图像的曲率测量方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510631492.9A CN105160693B (zh) | 2015-09-28 | 2015-09-28 | 一种用于医学图像的曲率测量方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105160693A CN105160693A (zh) | 2015-12-16 |
| CN105160693B true CN105160693B (zh) | 2018-04-27 |
Family
ID=54801535
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510631492.9A Active CN105160693B (zh) | 2015-09-28 | 2015-09-28 | 一种用于医学图像的曲率测量方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105160693B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110500202B (zh) * | 2019-08-19 | 2022-03-15 | 西安航天动力测控技术研究所 | 一种自动寻找固体火箭发动机燃烧终止时间的方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103693042A (zh) * | 2013-12-03 | 2014-04-02 | 重庆交通大学 | 一种基于前视轨迹曲率的山区复杂道路汽车行驶速度预测方法 |
| CN103839251A (zh) * | 2012-11-26 | 2014-06-04 | 三星电子株式会社 | 用于医学图像配准的方法和设备 |
| US8878929B2 (en) * | 2009-05-27 | 2014-11-04 | Koh Young Technology Inc. | Three dimensional shape measurement apparatus and method |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005185560A (ja) * | 2003-12-25 | 2005-07-14 | Konica Minolta Medical & Graphic Inc | 医用画像処理装置及び医用画像処理システム |
-
2015
- 2015-09-28 CN CN201510631492.9A patent/CN105160693B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8878929B2 (en) * | 2009-05-27 | 2014-11-04 | Koh Young Technology Inc. | Three dimensional shape measurement apparatus and method |
| CN103839251A (zh) * | 2012-11-26 | 2014-06-04 | 三星电子株式会社 | 用于医学图像配准的方法和设备 |
| CN103693042A (zh) * | 2013-12-03 | 2014-04-02 | 重庆交通大学 | 一种基于前视轨迹曲率的山区复杂道路汽车行驶速度预测方法 |
Non-Patent Citations (5)
| Title |
|---|
| A review of methods for quantitative evaluation of spinal curvature;Tomaz Vrtovec等;《Eur Spine J(2009)》;20090227;第18卷;第593-607页 * |
| Biomedical Image Processing of Magnetic Resonance Imaging of the Pathological Human Brain: An Intensity-Curvature Based Approach;Carlo Ciulla等;《ICT Innovations 2014》;20141231;第56-65页 * |
| B样条曲面等曲率线的高阶抽取方法;王世亮等;《第六届全国几何设计与计算学术会议论文集》;20130719;第1-5页 * |
| 人体尺寸测量数据的局部特征提取方法;孙晓鹏等;《计算机辅助设计与图形学学报》;20060228;第18卷(第2期);第194-199页 * |
| 基于伪Zenike矩的离线地图匹配算法;刘智等;《计算机工程与应用》;20141231;第50卷(第10期);第23-26页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105160693A (zh) | 2015-12-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20200098481A1 (en) | Methods and systems for predicting sensitivity of blood flow calculations to changes in anatomical geometry | |
| Takizawa et al. | Aorta modeling with the element-based zero-stress state and isogeometric discretization | |
| JP6386752B2 (ja) | 医用画像処理装置、医用画像診断装置及び画像処理方法 | |
| CN104699931B (zh) | 一种基于人脸的神经网络血压预测方法及手机 | |
| Nickisch et al. | Learning patient-specific lumped models for interactive coronary blood flow simulations | |
| US20230310085A1 (en) | Systems and methods for estimation of blood flow using response surface and reduced order modeling | |
| CN110838359A (zh) | 基于对话机器人的分诊方法、装置、存储介质及机器人 | |
| CN112711831A (zh) | 血管仿真分析方法、设备、装置、计算机设备和存储介质 | |
| WO2022001026A1 (en) | Systems and methods for determining blood vessel parameters | |
| CN110119540A (zh) | 一种用于生存风险分析的多输出梯度提升树建模方法 | |
| CN105160693B (zh) | 一种用于医学图像的曲率测量方法和装置 | |
| Fonken et al. | Ultrasound-based fluid-structure interaction modeling of abdominal aortic aneurysms incorporating pre-stress | |
| CN113870215A (zh) | 中线提取方法及装置 | |
| Pajaziti et al. | Shape-driven deep neural networks for fast acquisition of aortic 3D pressure and velocity flow fields | |
| CN110099502A (zh) | 一种x光成像设备的自适应控制方法及*** | |
| CN112163271B (zh) | 大气数据传感***的大气参数解算方法 | |
| Huang et al. | Aero-engine blade profile reconstruction based on adaptive step size bat algorithm and visualization of machining error | |
| CN115861471A (zh) | 一种血管曲面重建方法、血管曲面重建装置、设备及介质 | |
| CN111640102A (zh) | 动脉瘤检测方法、装置以及计算机设备 | |
| CN106251359A (zh) | 基于Clifford代数几何相对不变量的3D图像配准方法 | |
| EP3438989A1 (en) | Method and apparatus for predicting fluid flow through a subject conduit | |
| CN103971339B (zh) | 一种基于参数法的核磁共振图像分割方法及设备 | |
| JP7346419B2 (ja) | エンティティ間のコンテキスト的類似性の学習および適用 | |
| CN112419308A (zh) | 斑块评估方法、装置、电子设备及存储介质 | |
| Lim et al. | Empirical quantification of optic nerve strain due to horizontal duction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP03 | Change of name, title or address |
Address after: 110179 No. 177-1 Innovation Road, Hunnan District, Shenyang City, Liaoning Province Patentee after: DongSoft Medical System Co., Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Patentee before: Dongruan Medical Systems Co., Ltd., Shenyang |
|
| CP03 | Change of name, title or address |