CN1584936A - 一种医学图像的三维分割方法 - Google Patents
一种医学图像的三维分割方法 Download PDFInfo
- Publication number
- CN1584936A CN1584936A CN 200410020673 CN200410020673A CN1584936A CN 1584936 A CN1584936 A CN 1584936A CN 200410020673 CN200410020673 CN 200410020673 CN 200410020673 A CN200410020673 A CN 200410020673A CN 1584936 A CN1584936 A CN 1584936A
- Authority
- CN
- China
- Prior art keywords
- profile
- interpolation
- medical image
- carry out
- automatically
- 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
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000002452 interceptive effect Effects 0.000 claims abstract description 14
- 230000011218 segmentation Effects 0.000 claims description 26
- 238000003708 edge detection Methods 0.000 abstract description 6
- 230000001105 regulatory effect Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 11
- 238000012545 processing Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 210000004185 liver Anatomy 0.000 description 6
- 238000013461 design Methods 0.000 description 5
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 238000003709 image segmentation Methods 0.000 description 4
- 230000007812 deficiency Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000002059 diagnostic imaging Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000002687 intercalation Effects 0.000 description 2
- 238000009830 intercalation Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 210000000056 organ Anatomy 0.000 description 2
- 230000008520 organization Effects 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000003706 image smoothing Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 210000003734 kidney Anatomy 0.000 description 1
- 238000002386 leaching Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Images
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
- Image Analysis (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Image Processing (AREA)
Abstract
一种医学图像的三维分割方法,属于医学成像技术领域;针对现有技术对VOI的提取存在的不足,本发明方法提供一种医学图像的三维分割方法,采用Canny算子进行边缘检测,采用余弦函数进行轮廓插值,利用Canny算子和梯度信息生成图像的全局能量代价函数,作为交互式分割基础,实现对关键层的轮廓提取,进而在关键层之间利用基于余弦函数的插值方法产生中间层轮廓,同时在中间层轮廓上根据该层的局部能量代价自动生成种子点,并依据局部能量代价函数进行轮廓的自动调整,对插值生成的轮廓进行自动的修正,从而得到真正的分割区域;本发明方法能够对各种复杂的医学图像进行准确交互式的分割,具有快速性和可行性,且有极高的应用价值。
Description
技术领域
本发明属于医学成像技术领域,特别涉及一种医学图像的三维分割方法。
背景技术
随着医学成像技术的发展,图像分割在该领域的应用也越来越广泛,医学图像具有复杂性和多样性等特点,并且具有一定的噪音,图像中器官的边缘也会有局部的不清晰,这使得医学图像中VOI(volume of interest,体数据中感兴趣的区域)的提取比较困难,至今在医学图像方面还没有一种比较通用的分割方法。目前该领域的分割方法主要有两大类,一类以计算机为主体的自动的分割方法,另一类是人机结合的交互式分割方法。由于医学图像的复杂性,特别是针对复杂器官的提取,使得自动分割方法的使用受到了局限,并且其分割结果的准确性也不能满足要求;而交互式的分割方式相对而言具有更大的灵活性,因此该分割方法成为医学图像分割方面研究的一个热点。
自动分割的代表方法是Active Contour Models,即活动轮廓法,这类方法需要给出初始的轮廓,然后进行迭代,使轮廓沿能量降低的方向靠近,最后得到一个优化的边界。能量函数包括内外力两方面,如边界曲率和梯度。由于用户无法估计迭代的最后结果,应用活动轮廓法往往需要进行多次的交互工作。特别当目标比较复杂时,初始的轮廓不易确定,而迭代的结果往往不能达到要求。这种方法应用在医学图像中除了以上的缺陷外,时间的开销也是一个瓶颈,不能保证快速的进行三维图像的提取。
人机结合的交互式分割的代表方法是Intelligent Scissors(IS),即智能剪刀法,该方法的基本思想是设定种子点,即设定目标区域边缘任意点,然后计算出图上各点到种子点的最小能量。由于边缘点多为零交叉点,其能量较邻点小,因而使能量图呈“总线”结构,其中边缘为代表低能量的“总线”,从而在交互引导过程中自动勾勒出边缘。对于轮廓比较平滑的目标,这种方法可以获得比较好的效果,但在处理具有尖锐边角的目标时,对尖锐边角往往不能准确定位。
在三维重建的过程中智能剪刀法只对起始层和终止层进行交互式轮廓分割,对于起始层和终止层之间的中间层轮廓采用了线性插值方式生成,进而再进行三维重建。传统插值是一种线性的加权平均,只是沿着切片前进的方向进行简单的以切片间距离为量度的一种线性变换,这种变换会使得插值生成的轮廓在实际轮廓有突变处不能及时的得到改变,使分割的准确性受到了限制。
发明内容
针对现有技术存在的问题,本发明提供一种医学图像的三维分割方法,能够对各种复杂医学图像进行准确交互式的分割。
由于传统的智能剪刀方法在进行分割之前采用了拉普拉斯算子进行预处理,然而拉普拉斯算子并不是一种理想的边缘检测算子,它对图像平滑滤波的程度取决于高斯函数的参数δ,δ值越大,噪声滤波效果越好,但同时也丢失了重要的边缘信息,在平滑相互邻近的两个边缘时还可能把它们连接在一起而只检测到一个边缘。如果选取小δ值,又有可能平滑不完全而留有太多噪声。也就是说抑制噪声和边缘精确定位无法同时得到满足。Canny证明了高斯函数的一阶倒数可以在抗干扰和精确定位之间选择一个折衷的方案,并在此基础上提出了对信噪比与定位之乘积的最优化逼近算子,可以在有效抑制噪声的前提下,最大限度的保证边缘的连续性。因此本发明采用Canny算子(一种边缘检测算子)代替拉普拉斯算子对智能剪刀法进行了改进,使之更加符合医学图像分割的需要。观察图2a,这是一幅典型的CT图像。由于本身组成成分的复杂性,各组织内部的灰度分布很不均匀。由于相邻的组织不同,组织之间的边缘时而清晰,时而模糊。图2b是采用拉普拉斯边缘检测算子进行边缘检测的结果,由于图像的上述特征导致检测结果存在两个问题:第一,组织内部存在大量的伪边缘;第二,边缘存在大量间断。这两个问题会严重影响智能剪刀法交互分割时的分割结果,一方面由于伪边缘的吸引会使动态轮廓线不沿着真正的边缘前进,需要用户增加交互的次数并产生不正确的轮廓;另一方面由于边缘的间断,会导致轮廓线的准确性和光滑性严重下降。图2c是采用Canny算子进行边缘检测的结果。可以看到此时大部分的伪边缘都被去掉了,而检测出来的边缘本身具有很好的连续性,这将大大有利于提高智能剪刀法在分割医学图像时的准确性。
对于传统插值带来的不足,本发明采用基于余弦函数的轮廓插值来取代线性加权平均插值,利用这种插值方法可以明显的改善线性插值存在的不足。余弦插值函数曲线如图3所示,插值函数
其中m=1,2......l-1,l代表总层数。观察基于线性插值图4b和基于余弦插值图5b,可见对于中间插值生成的轮廓,余弦插值可以很好的解决线性插值在轮廓突变处的不足。
为实现上述目的,本发明方法利用Canny算子和梯度信息生成图像的全局能量代价函数,作为交互式分割基础,实现对关键层的轮廓提取,进而在关键层之间利用基于余弦函数的插值方法产生中间层轮廓,同时在中间层轮廓上根据该层的局部能量代价自动生成种子点,并依据局部能量代价函数进行轮廓的自动调整,对插值生成的轮廓进行自动的修正,从而得到真正的分割区域。整个过程将自动分割和交互式的分割相结合,相对于活动轮廓模型更加快速,相对于智能剪刀法更加准确。实践证明本发明提出的方法具有快速性和可行性,且有极高的应用价值。
医学图像的三维分割方法由以下步骤组成:
(一)在关键层利用canny算子和梯度场计算能量代价场,进行交互式分割;
(二)在关键层间进行余弦插值,在各层上产生初始轮廓;
(三)在插值层内根据初始轮廓自动确定种子点,并通过局部轮廓自动调整确定最终完整轮廓。
其中步骤(一)又包括以下步骤:
(1)利用canny算子和邻域梯度进行梯度场的计算,从而得到能量代价场;
(2)选定关键层进行交互式分割。
其中步骤(二)的又包括以下步骤:
(1)利用交互式分割在关键层确定初始轮廓;
(2)利用改进的余弦插值在中间层进行轮廓插值;
(3)对插值轮廓进行细化。
其中步骤(三)又包括以下步骤:
(1)在插值产生的轮廓内部自动确定种子点;
(2)对这些种子点在邻域内自动调整到能量最小的位置;
(3)在每两个种子点之间通过局部能量代价计算轮廓,即轮廓局部自动调整方法。
本发明方法是一种针对复杂医学图像的三维分割方法,对智能剪刀法和活动轮廓法进行了融合改进,实现了快速准确的交互式分割,在医学辅助诊断系统和医学三维重建系统中有着很高的应用价值。
附图说明
图1为医学图像的三维分割方法的流程图;
图2为图像组织的边缘提取图,其中a为原图,b为采用拉普拉斯边缘检测算子进行边缘检测的结果图,c为采用canny边缘检测算子进行边缘检测的结果图;
图3余弦插值函数曲线;
图4是图像经过线性插值处理的结果,其中a为起始层,b为轮廓插层,c为终止层;
图5是图像经过余弦插值处理的结果,其中a为起始层,b为轮廓插层,c为终止层;
图6是为获得肝脏图像而提取的关键层轮廓图;
图7是为获得肝脏图像而进行插值处理的边缘轮廓图;
图8是为获得肝脏图像而进行细化处理的边缘轮廓图;
图9是为获得肝脏图像而进行局部轮廓自动调整的效果图,其中a为调整前效果,b为调整后效果;
图10为应用本发明方法提取的肝脏三维重建结果;
图11是为获得肾脏图像而提取的关键层轮廓图;
图12是为获得肾脏图像而进行插值处理的边缘轮廓图;
图13是为获得肾脏图像而进行细化处理的边缘轮廓图;
图14是为获得肾脏图像而进行局部轮廓自动调整的效果图,其中a为调整前效果,b为调整后效果;
图15为应用本发明方法提取的肾脏三维重建结果。
具体实施方式
结合附图,如图1所示医学图像的三维分割方法的流程图,本发明详细的分割方法,包括以下三个步骤:
(一)在关键层利用canny算子和梯度场计算能量代价场,进行交互式分割;
(二)在关键层间进行余弦插值,在各层上产生初始轮廓;
(三)在插值层内根据初始轮廓自动确定种子点,并通过局部轮廓自动调整确定最终完整轮廓。
本发明方法具体实现过程如下:
步骤(一):
l(p,q)=wZ·fZ(q)+wD·fD(p,q)+wG·fG(q) (2)
其中fZ(q)、fG(q)、fD(p,q)分别表示零交叉点、梯度值、梯度方向,wZ、wD、wG为权值,其参考系数分别为0.4、0.2、0.4,p、q为图像中的像素点。
式(2)中fZ(q)的计算为:
其中IL是对原图像进行Canny变换得到的结果。
式(2)中fG(q)的计算为:
计算梯度时可以采用不同的算子,本方法中令Ix、Iy表示x、y方向的梯度,则梯度G计算为:
为使高梯度产生低能量,令
其中G=G-min(G) (6)
式(2)fD(p,q)的计算为:
因为梯度方向能量是对变化剧烈边缘的一种平滑,且当p、q两点相似时,其梯度方向所占的能量较小,因此这种能量因子的最终影响结果是使相似点归于同一条能量路径。
式(7)中 dp(p,q)=D′(p)·L(p,q) (8)
dq(p,q)=L(p,q)·D′(q) (9)
其中“·”表示矢量点乘。其中矢量
D′(p)=(Iy(p),-Ix(p)) (10)
D′(q)=(Iy(q),-Ix(q)) (11)
Ix(p),Iy(p)、Ix(q),Iy(q)分别表示p、q两点x、y方向的梯度,L(p,q)为p、q两点的单位矢量,D′(p)、D′(q)为梯度矢量顺时针旋转90度后的矢量,利用该函数对关键层进行轮廓的提取,效果如图6。
步骤(二):
在对起始层和终止层的轮廓进行完交互式分割后,中间层采用插值方式生成,在两个相邻的断层图像上抽取轮廓,分别得到两个二值图像,其中“0”像素代表轮廓点。对这两幅二值图像进行距离变换,可以分别得到两个距离图像d1(i,j),d2(i,j),与正规的距离变换不同的是:对于轮廓内的点,它的值是它到轮廓上最近一点距离的负值;而对于轮廓外的点,它的值是它到轮廓上最近一点距离的正值。为节省计算时间,在这里将采用街区距离变换计算最短距离,对于中间断层上的距离图像是通过对两幅距离图像加权平均得到的,即:
dk(i,j)=(l-f)d1(i1,j1)+fd2(i2,j2) (13)
其中,k=1,2......l-1,m=1,2......l-1,l代表总层数,dk(i,j)就是中间断层上的距离图像。
在距离图像中,距离大于零的区域代表轮廓的外部,距离小于零的区域代表轮廓的内部,而轮廓就是过零点构成的,效果如图7。
为了确保插值后的边缘轮廓是单边缘像素宽度,本发明对插值生成轮廓进行了细化处理,以确保在步骤三中局部轮廓自动调整的准确性,效果如图8。
步骤(三):
根据医学图像中前后两层相似性较大的特点,取上一层的种子点为下一层新的初始种子点,然后进行自动调整。
同步骤(一),梯度G表示为:
G′=Max(Gi) (14)
其中Gi为种子点邻域内点的梯度集合
G′取最大值时,{ix,iy}即为种子点的调整位置,然后根据步骤(一)确定最终轮廓,效果如图9。
经过上述三个步骤,对于图层数为123、层厚2mm、像素间距0.6835mm、层间距0.7mm、数据场规模350*350*86的肝脏图像处理后在三维下进行显示的效果图如图10。
图11到图15是利用本发明方法对于图层数为150、层厚2mm、像素间距0.6835mm、层间距0.7mm、数据场规模350*350*105的肾脏图像提取过程中的效果图及处理后的三维显示效果图。
Claims (5)
1.一种医学图像的三维分割方法,其特征在于本发明方法由以下三个步骤组成:
(一)在关键层利用canny算子和梯度场计算能量代价场,进行交互式分割;
(二)在关键层间进行余弦插值,在各层上产生初始轮廓;
(三)在插值层内根据初始轮廓自动确定种子点,并通过局部轮廓自动调整确定最终完整轮廓。
2.根据权利要求1所述的一种医学图像的三维分割方法,其特征在于步骤(一)由以下步骤组成:
(1)利用canny算子和邻域梯度进行梯度场的计算,从而得到能量代价场;
(2)选定关键层进行交互式分割。
3.根据权利要求1所述的一种医学图像的三维分割方法,其特征在于步骤(二)由以下步骤组成:
(1)利用交互式分割在关键层确定初始轮廓;
(2)利用改进的余弦插值在中间层进行轮廓插值;
(3)对插值轮廓进行细化。
4.根据权利要求1所述的一种医学图像的三维分割方法,其特征在于步骤(三)由以下步骤组成:
(1)在插值产生的轮廓内部自动确定种子点;
(2)对这些种子点在邻域内自动调整到能量最小的位置;
(3)在每两个种子点之间通过局部能量代价计算轮廓,即轮廓局部自动调整方法。
5.根据权利要求1所述的一种医学图像的三维分割方法,其特征在于所述的余弦插值的函数为:
其中m=1,2.....l-1,l代表总层数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200410020673XA CN1301494C (zh) | 2004-06-07 | 2004-06-07 | 一种医学图像的三维分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB200410020673XA CN1301494C (zh) | 2004-06-07 | 2004-06-07 | 一种医学图像的三维分割方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1584936A true CN1584936A (zh) | 2005-02-23 |
CN1301494C CN1301494C (zh) | 2007-02-21 |
Family
ID=34600571
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB200410020673XA Expired - Fee Related CN1301494C (zh) | 2004-06-07 | 2004-06-07 | 一种医学图像的三维分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1301494C (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100550051C (zh) * | 2007-04-29 | 2009-10-14 | 威盛电子股份有限公司 | 图像变形的方法 |
CN102289830A (zh) * | 2010-06-07 | 2011-12-21 | 微软公司 | 利用测地仿射性的数据驱动内插 |
CN101766494B (zh) * | 2008-12-31 | 2012-12-12 | 西门子公司 | 使高强度聚焦超声设备的操作参数可再度使用的方法 |
CN103854299A (zh) * | 2014-03-05 | 2014-06-11 | 清华大学 | 一种基于全局能量函数优化的图像层次提取与编辑方法 |
CN104143190B (zh) * | 2014-07-24 | 2017-06-09 | 东软集团股份有限公司 | Ct图像中组织的分割方法及系统 |
CN108122240A (zh) * | 2017-12-19 | 2018-06-05 | 东软集团股份有限公司 | 一种三维图像中的目标区域的确定方法、装置及设备 |
CN108447049A (zh) * | 2018-02-27 | 2018-08-24 | 中国海洋大学 | 一种基于生成式对抗网络的数字化生理生物体分割方法 |
CN110050293A (zh) * | 2016-12-13 | 2019-07-23 | 爱克发医疗保健公司 | 在医学辐射图像中分割3d对象的方法 |
CN110147077A (zh) * | 2019-04-19 | 2019-08-20 | 深圳科瑞技术股份有限公司 | 一种工业机器人在空间直线运行下的余弦插补方法 |
CN111295690A (zh) * | 2017-11-09 | 2020-06-16 | 爱克发医疗保健公司 | 用于在医学辐射图像中分割3d对象的方法 |
CN116912282A (zh) * | 2023-06-05 | 2023-10-20 | 北京纳通医用机器人科技有限公司 | 基于医学图像的三维分割方法、装置、设备及存储介质 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101826209B (zh) * | 2010-04-29 | 2011-12-21 | 电子科技大学 | 一种基于Canny模型的三维医学图像分割方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6331116B1 (en) * | 1996-09-16 | 2001-12-18 | The Research Foundation Of State University Of New York | System and method for performing a three-dimensional virtual segmentation and examination |
US6754376B1 (en) * | 2000-11-22 | 2004-06-22 | General Electric Company | Method for automatic segmentation of medical images |
US7386153B2 (en) * | 2001-11-23 | 2008-06-10 | Infinitt Co., Ltd. | Medical image segmentation apparatus and method thereof |
CN100541541C (zh) * | 2001-12-29 | 2009-09-16 | 田捷 | 基于单层表面跟踪的超大规模医学图像表面重建方法 |
CN1224943C (zh) * | 2002-07-26 | 2005-10-26 | 中国科学院自动化研究所 | 基于活动轮廓模型的多个对象的自动分割方法 |
-
2004
- 2004-06-07 CN CNB200410020673XA patent/CN1301494C/zh not_active Expired - Fee Related
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100550051C (zh) * | 2007-04-29 | 2009-10-14 | 威盛电子股份有限公司 | 图像变形的方法 |
CN101766494B (zh) * | 2008-12-31 | 2012-12-12 | 西门子公司 | 使高强度聚焦超声设备的操作参数可再度使用的方法 |
CN102289830A (zh) * | 2010-06-07 | 2011-12-21 | 微软公司 | 利用测地仿射性的数据驱动内插 |
CN102289830B (zh) * | 2010-06-07 | 2016-06-01 | 微软技术许可有限责任公司 | 利用测地仿射性的数据驱动内插 |
CN103854299A (zh) * | 2014-03-05 | 2014-06-11 | 清华大学 | 一种基于全局能量函数优化的图像层次提取与编辑方法 |
CN103854299B (zh) * | 2014-03-05 | 2017-04-19 | 清华大学 | 一种基于全局能量函数优化的图像层次提取与编辑方法 |
CN104143190B (zh) * | 2014-07-24 | 2017-06-09 | 东软集团股份有限公司 | Ct图像中组织的分割方法及系统 |
CN110050293A (zh) * | 2016-12-13 | 2019-07-23 | 爱克发医疗保健公司 | 在医学辐射图像中分割3d对象的方法 |
CN111295690A (zh) * | 2017-11-09 | 2020-06-16 | 爱克发医疗保健公司 | 用于在医学辐射图像中分割3d对象的方法 |
CN108122240A (zh) * | 2017-12-19 | 2018-06-05 | 东软集团股份有限公司 | 一种三维图像中的目标区域的确定方法、装置及设备 |
CN108122240B (zh) * | 2017-12-19 | 2020-06-12 | 东软集团股份有限公司 | 一种三维图像中的目标区域的确定方法、装置及设备 |
CN108447049A (zh) * | 2018-02-27 | 2018-08-24 | 中国海洋大学 | 一种基于生成式对抗网络的数字化生理生物体分割方法 |
CN110147077A (zh) * | 2019-04-19 | 2019-08-20 | 深圳科瑞技术股份有限公司 | 一种工业机器人在空间直线运行下的余弦插补方法 |
CN110147077B (zh) * | 2019-04-19 | 2021-11-16 | 深圳市科瑞软件技术有限公司 | 一种工业机器人在空间直线运行下的余弦插补方法 |
CN116912282A (zh) * | 2023-06-05 | 2023-10-20 | 北京纳通医用机器人科技有限公司 | 基于医学图像的三维分割方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN1301494C (zh) | 2007-02-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108053417B (zh) | 一种基于混合粗分割特征的3D U-Net网络的肺分割装置 | |
CN104809723B (zh) | 基于超体素和图割算法的三维肝脏ct图像自动分割方法 | |
CN1301494C (zh) | 一种医学图像的三维分割方法 | |
CN106780518B (zh) | 一种基于随机游走和图割的活动轮廓模型的mr图像三维交互分割方法 | |
CN100561518C (zh) | 基于感兴趣区域的自适应医学序列图像插值方法 | |
CN1266655C (zh) | 基于图像轮廓的人肢体三维建模方法 | |
CN102136142B (zh) | 基于自适应三角形网格的非刚性医学图像配准方法 | |
CN102881014B (zh) | 一种基于图割的快速立体匹配方法 | |
CN107067393B (zh) | 一种基于用户交互和形状先验知识的三维医学图像分割方法 | |
CN105741299B (zh) | 一种冠状动脉ct血管造影图像分割方法 | |
CN105957063A (zh) | 基于多尺度加权相似性测度的ct图像肝脏分割方法及系统 | |
CN106485695A (zh) | 基于统计形状模型的医学图像Graph Cut分割方法 | |
CN106485721B (zh) | 从光学相干断层图像获取视网膜结构的方法及其系统 | |
CN107230206A (zh) | 一种基于多模态数据的超体素序列肺部图像的3d肺结节分割方法 | |
CN112140552B (zh) | 一种原位3d打印曲面分层和路径规划方法 | |
CN102426700B (zh) | 基于局部和全局区域信息的水平集sar图像分割方法 | |
CN108053433A (zh) | 一种基于物理对齐和轮廓匹配的多模态颈动脉mri配准方法 | |
CN105389811A (zh) | 一种基于多级阈值分割的多模态医学图像处理方法 | |
CN103295224B (zh) | 一种基于均值漂移和分水岭的乳腺超声图像自动分割方法 | |
CN101303769B (zh) | 基于先验知识测地几何流的二维序列医学图像分割方法 | |
CN105761304B (zh) | 三维脏器模型构造方法和装置 | |
CN101052991A (zh) | 使用距离坐标的特征加权的医学对象轮廓确定 | |
CN105279762A (zh) | 一种口腔软硬组织ct序列与三维网格模型配准方法 | |
CN103714547B (zh) | 一种结合边缘区域和互相关的图像配准方法 | |
CN103544695B (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070221 Termination date: 20110607 |