CN104408734A - 联合图像分割与形变配准技术的自适应目标区域转换方法 - Google Patents
联合图像分割与形变配准技术的自适应目标区域转换方法 Download PDFInfo
- Publication number
- CN104408734A CN104408734A CN201410764128.5A CN201410764128A CN104408734A CN 104408734 A CN104408734 A CN 104408734A CN 201410764128 A CN201410764128 A CN 201410764128A CN 104408734 A CN104408734 A CN 104408734A
- Authority
- CN
- China
- Prior art keywords
- phi
- image
- function
- deformation
- sigma
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/77—Determining position or orientation of objects or cameras using statistical methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
-
- 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/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- 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
- G06T2207/30056—Liver; Hepatic
-
- 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
- G06T2207/30096—Tumor; Lesion
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Image Processing (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种联合图像分割与形变配准技术的自适应目标区域转换方法,本方法自动分割危及器官及胸腹部肿瘤区域,采用形变配准算法,将CBCT图像与计划CT图像进行多模态配准,通过计算出的形变域将计划CT上的信息自动转移到CBCT上。本发明在计划CT与CBCT图像刚体配准后,为补偿不能恢复的形变范围,采用基于正交小波基函数的形变配准算法,通过计算出的形变域,实现危及器官与肿瘤的轮廓以及放疗计划从计划CT到CBCT的转移,精确定位靶区位置。
Description
技术领域
本发明涉及一种联合图像分割与形变配准技术的自适应目标区域转换方法。
背景技术
临床精确放射治疗的目标是在最大可能地提高肿瘤靶区射线吸收剂量的同时尽可能地降低危及器官的射线接收剂量。针对胸腹部肿瘤,目前临床上通常的做法是先CT定位,再按制定的放疗计划实施逐次照射。但在放疗过程中,由于器官运动和变形或各种客观因素使肿瘤和危及器官(主要是肝脏)偏离射野,导致肿瘤欠剂量照射而危及器官却过剂量。
利用近些年发展起来的图像引导放疗(Image guided radiation therapy,IGRT)技术,可对肿瘤靶区和危及器官进行实时监控,根据器官位置的变化调整治疗条件达到使照射野跟踪靶区的目的。临床上主要采用基于CBCT(Cone beam CT)的图像引导放疗技术,但目前该技术仅限用于调整摆位误差,即将计划CT与放疗前获取的KV级CBCT进行基于骨形的刚体配准,再通过计算放疗肿瘤靶区的位置偏移量来指导摆位;同时,目前的放疗计划也仅仅基于静态的计划CT,并未做到真正意义上的“自适应”。此外,计划CT上的危及器官及肿瘤也均由医生手工勾画,费时费力,具有主观性。
发明内容
本发明为了解决上述问题,提出了一种联合图像分割与形变配准技术的自适应目标区域转换方法,本方法自动分割危及器官及胸腹部肿瘤区域,采用形变配准算法,将CBCT图像与计划CT图像进行多模态配准,通过计算出的形变域将计划CT上的信息自动转移到CBCT上。
为了实现上述目的,本发明采用如下技术方案:
一种联合图像分割与形变配准技术的自适应目标区域转换方法,包括以下步骤:
(1)采用基于图谱的自动分割方法,对采集的N组CT图像去除伪影、噪声成分,采用梯度矢量流(Gradient vector flow,GVF)进行胸肋骨自动提取;
(2)在N组图像中随机选取一组作为参考图像,其余作为目标图像,将得到的肝脏训练集和胸肋骨分别进行肝脏图谱和胸肋骨图谱的创建;
(3)输入患者当前计划CT图像,采用GVF进行胸肋骨自动提取,将胸肋骨图谱作为参考图像,患者计划CT图像作为目标图像,以胸肋骨为特征,采用仿射配准算法与胸肋骨图谱进行刚体配准,得到患者计划CT图像中肝脏区域的初始值;
(4)使用基于正交小波基函数的形变配准算法对得到的肝脏区域与肝脏图谱进行形变配准,使肝脏边界更为精准,通过基于图谱的分割方法在患者计划CT上自动分割出肝脏区域;
(5)结合水平集模型与似然函数模型在计划CT上自动分割胸腹部肿瘤区域;
(6)以胸肋骨为特征,通过仿射变换将最新CBCT图像与计划CT图像进行刚体配准,为补偿刚体配准时不能恢复的形变范围,使用基于正交小波基函数的形变配准算法,计算出CBCT与计划CT的形变域,将计划CT中的危及器官和肿瘤的轮廓以及等剂量线信息自动转移到CBCT上。
所述步骤(2)中,具体方法为:
在N组满足条件的CT中随机选取一组作为参考图像,将其与其余的目标图像采用仿射变换进行配准,所有的图像对齐后,统计空间中每个像素存在的概率;N组图像数据配准后,得到一个概率图谱;采用多次循环迭代,以提高所建图谱的可靠性;
建立概率图谱后,根据具体应用环境选取合适的阈值,按照灰度值随图谱概率增大而增大的准则,把图谱从概率空间转化到灰度空间。
所述步骤(5)中,具体方法为:
采用似然函数与水平集函数相结合的能量函数,在患者计划CT图像上自动获取胸腹部肿瘤信息;似然能量函数Ep(φ)由感兴趣区域能量函数eT(x)和背景区域能量函数eB(x)两部分构成;用似然函数的有限高斯模型估计感兴趣肿瘤区域密度分布,似然函数的高斯混合模型估计胸腹部背景区域密度分布,通过以上似然函数的估计来划分图像概率密度;水平集函数用于获得边缘轮廓信息,防止肿瘤边界泄漏。
所述步骤(5)中,肿瘤区域的密度分布采用有限高斯模型估计,具体的表达式为:
其中,ΩT为肿瘤区域,p(X)为肿瘤区域像素的概率,μT为高斯函数均值,σT为标准差,u0(x)为患者CT图像;
背景区域密度分布采用高斯混合模型估计,具体的表达式为:
其中,ΩB为背景区域,p(X)为背景区域像素的概率,和分别为高斯模型的均值和标准差,K为背景区域的分类个数,为分类的混合系数;
似然能量函数的具体表达式为:
其中,e1(x)=-eT(x),
距离正则化水平集能量函数的具体表达式为:
其中:μ、λ、α分别是正则化项、长度项、面积
项的权重系数,
当似然能量函数与水平集能量函数最小时,水平集函数演化曲线收敛到感兴趣区域边界,水平集函数演化表达式为:
其中,div(·)为散度算子。
所述步骤(6)中,具体方法为:
将已导入加速器的计划CT作为目标图像,CBCT图像作为参考图像,采用GVF对CBCT图像进行胸肋骨自动提取,以胸肋骨为特征,利用仿射配准算法与患者计划CT进行刚体配准;为补偿刚体配准时不能恢复的形变范围,使用基于正交小波基函数的形变配准算法,获得形变域,将形变域三维小波分解后,使小波系数的初始值为0;采用非线性最小二乘LM(Levenberg-Marquadt,LM)优化算法,寻找合适的小波系数,根据所需精度控制在某个尺度上收敛;由估计到的小波系数得到参考图像与目标图像之间的形变域;将计划CT中的危及器官和肿瘤的轮廓以及等剂量线信息自动转移到CBCT上,进而实现剂量体积直方图(Dosevolume histograms,DVH)分析。
所述步骤(6)中,在形变配准时,使用Navier偏微分方程描述发生形变的各向同性物质的平衡状态,通过Navier偏微分方程建立形变能量函数,将形变域建模为小波系数的函数,其具体表达式为:
其中,(x,y,z)表示参考图像的空间坐标系,(x′,y′,z′)表示目标图像的空间坐标系,u表示形变域,c表示小波系数,u是c的函数。
所述Navier偏微分方程的外力具体表达式为:
H(X)是图像X的香农熵;
Navier偏微分内力方程具体表达式为:
其中,λ和μ是Lame系数;
形变能量函数具体表达式为:
E(c)=inter(c)+w*exter(c)
其中,w是加权常数,本发明w为常数1,inter(c)表示内力约束函数,exter(c)表示外力约束函数。
本发明的有益效果为:
(1)在患者计划CT上自动分割危及器官与胸腹部肿瘤,避免手工勾画引入的人为误差,保证剂量使用的准确性;
(2)在计划CT与CBCT图像刚体配准后,为补偿不能恢复的形变范围,采用基于正交小波基函数的形变配准算法,通过计算出的形变域,实现危及器官与肿瘤的轮廓以及放疗计划从计划CT到CBCT图像的自适应转移,精确定位靶区位置。
附图说明
图1为本发明方法流程图;
图2为创建图谱流程图;
图3为基于小波基函数形变配准流程图。
具体实施方式:
下面结合附图与实施例对本发明作进一步说明。
对于胸腹部肿瘤,进行放射治疗时主要的危及器官为肝脏,通过50组CT图像建立正常肝脏图谱,采用基于图谱的分割算法实现患者计划CT中肝脏组织区域的自动分割;同时通过水平集模型与似然函数模型相结合的方法在计划CT上自动提取胸腹部肿瘤区域,达到自动分割计划CT的目的。之后使用基于正交小波基函数的形变配准方法,将计划CT中各种信息通过计算出的形变域自动转移到CBCT上。
具体步骤如图1所示:
第一步:为实现胸腹部肿瘤患者计划CT中肝脏组织的自动分割,本发明采用基于图谱的自动分割方法,根据患者的年龄、性别、身高、体重、病种、治疗状况、种族和区域等因素对100组增强型计划CT图像进行详细地分类筛选,选取其中50组条件相近,年龄在20-40岁之间且肝脏组织正常的CT图像,这些CT图像中的肝脏区域均被多位物理师和医学专家详细勾画与核对,可确保图谱的准确性。从中随机选择其中一组作为参考图像,其余49组作为目标图像,并去除伪影、噪声成分;采用梯度矢量流GVF进行胸肋骨自动提取。
第二步:对以上50组数据的肝脏训练集和胸肋骨分别进行肝脏图谱和胸肋骨图谱的创建。将随机选取的一组参考图像与其余的49组目标图像采用仿射变换进行配准,所有图像对齐后,统计空间中每个像素存在的概率。为了提高所建图谱的可靠性,本发明采用两次循环进行迭代,以减少由于采用特定初始参考图像而引入的误差。
使用上述同样的方法建立胸肋骨图谱。
建立概率图谱后,根据具体应用环境,选择合适的阈值,按照灰度值随图谱概率增大而增大的准则,将图谱从概率空间转化到灰度空间。选取阈值时要避免选取过高排除连带边缘组织,而选取过低又会包含其他组织的情况。
第三步:输入患者计划CT图像,采用GVF进行胸肋骨自动提取。以胸肋骨为特征,利用仿射配准算法与胸肋骨图谱进行刚体配准,在计划CT图像上得到肝脏区域的初始值。
第四步:为得到计划CT中更为精准的肝脏区域边界,本发明将计划CT作为目标图像,肝脏图谱作为参考图像,使用基于正交小波基函数的形变配准算法,实现患者计划CT肝脏组织与肝脏图谱的形变配准;根据以上几步所述,通过基于图谱的分割方法便可在计划CT上自动获取肝脏位置。基于正交小波基函数的形变配准算法与第七步相同,其具体详情参照步骤(7)。
第五步:构建似然函数与水平集函数相结合的能量函数,用来自动获取计划CT图像中胸腹部肿瘤区域。
本发明将计划CT图像用各个相互独立的区域组成,这些区域分为感兴趣区域ΩT和背景区域ΩB,感兴趣区域即计划CT图像中的肿瘤区域。采用的分割模型由似然能量函数和水平集能量函数两部分构成,具体表达式如下:
E(φ)=Ep(φ)+Eε(φ) (1)
其中,Ep(φ)为似然能量函数,用以估计各像素的概率,进而区分该像素是属于肿瘤区域还是背景区域,似然能量函数包括感兴趣区域能量函数和背景区域能量函数两部分。Ee(φ)为局部空间边缘能量函数,用以防止分割时边界泄漏。
依据胸腹部CT图像全局密度分布多峰性特点,可用多个高斯分布的叠加拟合实现肿瘤分割。本发明使用有限高斯模型估计肿瘤区域的密度分布,具体表达式为:
其中,p(X)为肿瘤区域像素的概率,μT为高斯函数均值,σT为标准差,u0(x)为患者CT图像。然后使用高斯混合模型估计背景区域密度分布,具体表达式为:
其中,p(X)为背景区域像素的概率,和
分别为高斯模型的均值和标准差,K为背景区域的分类个数,为分类的混合系数。
结合水平集方法,用水平集函数的正负表示目标函数和背景函数,目标函数ΩT={x,φ(x)≥0},背景函数ΩB={x,φ(x)<0}。通过水平集函数使用阶梯函数H(φ)(HeavisideFunction)以区分ΩT和ΩB,因此可得出似然能量函数表达式:
其中,e1(x)=-eT(x),
接下来本发明采用距离正则化水平集模型以获得肿瘤边缘轮廓信息,水平集能量函数的
具体表达式为:
其中:为加权长度项,用以驱动零水平集
曲线向目标边界演化;为加权面积项,用以加快演化速度。为正则化
项,用以纠正水平集函数与符号距离函数的误差,保证水平集函数在整个曲线演化过程中最
大程度地接近符号距离函数,避免在演化过程中重新初始化。μ、λ、α分别是正则化项、
长度项、面积项的权重系数。
通过以上所述即可得出本发明中似然能量函数与水平集能量函数相结合的分割模型能量函数,当该能量函数最小时,水平集演化曲线收敛到感兴趣区域边界。水平集函数演化表达式如下:
其中,div(·)为散度算子。
第六步:通过治疗计划***TPS(Treatment planning system,TPS)在已自动分割危及器官和肿瘤区域的计划CT上制定放疗计划,设计照射野。
第七步:患者开始接受治疗时,将已导入加速器的计划CT作为目标图像,治疗日CBCT图像作为参考图像,采用GVF对CBCT图像进行胸肋骨自动提取,以胸肋骨为特征,通过仿射变换与患者计划CT进行刚体配准,为补偿刚体配准时不能恢复的形变范围,本发明使用与上述第四步自动分割方法中相同的形变配准算法,即基于正交小波基函数的形变配准算法,计算形变域,将计划CT中各种信息自动转移到CBCT图像上,其结果可对离线状态下的放疗计划进行修正或更新,避免人工修订的误差,在CBCT中自动产生危及器官和肿瘤靶区的轮廓线,并自动生成相应等剂量线,进行剂量体积直方图分析,从而有效指导自适应放疗。基于正交小波基函数的形变配准算法的具体流程如图3所示。
计划CT与CBCT刚体配准后,对其进行小波分解,获得形变域。本发明采用接近三次正交的不完全正交三阶样条小波来表示形变域,以减少极小化能量函数过程的复杂度,将形变域三维小波分解后,小波系数的初始值设为0。采用非线性最小二乘LM(Levenberg-Marquadt,LM)优化算法递归的极小化公式,寻找合适的小波系数,以小波系数的尺度和子带作为循环参数进行递归,根据所需精度控制在某个尺度上收敛。最终由估计到的小波系数得出参考图像与目标图像之间的形变域。
小波系数的函数通过形变域建模,参考图像、目标图像、形变域,小波系数之间的数学表达如下:
其中,(x,y,z)表示参考图像的空间坐标系,(x′,y′,z′)表示目标图像的空间坐标系。形变域u=(u1,u2,u3),c表示小波系数,形变域u为小波系数c的函数。作为目标参数,它通过极小化能量函数获得。
对多尺度小波系数进行估计来建立形变能量函数,函数由两部分组成,其中一项是内力形变能量函数,另外一项是外力约束能量函数。通过内力产生各种自由形变,通过外力约束达到平衡状态。在形变配准中,使用Navier偏微分方程描述发生形变的各向同性物质的平衡状态,具体的Navier偏微分方程表达如下:
其中,θ是三维膨胀系数,表达式为:
其中,X=(x1,x2,x3)T是形变组织的三维空间坐标系,F=(F1,F2,F3)T是作用于形变组织各处的外力,μ和λ均为形变常量,其取值取决于形变组织的自身性质。u=(u1,u2,u3)T是配准算法所需恢复目标对象与参考对象之间的形变域。在发生形变的组织边缘,外力近似为零,对应医学图像的器官轮廓边界;而对于发生形变的其它部位,均需选择合适的内力和外力模型,并用适当的方法表示形变域。对于给定形变域的受力和表示方法,物体的形变域完全由形变常数和外力决定。
Navier偏微分方程前两项是造成组织发生形变的内力部分,将(9)代入(8)中前两项,得到的内力方程如下:
其中,λ和μ是Lame系数,其取值由弹性形变物质的自身性质决定。本发明设置λ=0、μ=1,这种情况下的医学图像性质满足各向同性物质属性,且不影响算法的精度。
本发明为多模态配准,由于归一化互信息与待配准图像数据对的重叠程度无关,鲁棒性较好,因此采用基于归一化互信息的相似性测度,具体表达式为:
其中,H(X)是图像X的香农熵。通过极小化模拟平衡状态的能量函数来估计目标函数中的小波系数,因此,外力用归一化互信息的倒数表示:
exter(c)=1/NMI(X,X(u)) (12)
由此便得到用于构成Navier偏微分方程来解决弹性形变问题的内力和外力部分,其中形变域用正交小波基函数表示,通过内力和外力的平衡最终得出形变域的能量函数,能量函数具体的表达式为:
E(c)=inter(c)+w*exter(c) (13)
其中,w是加权常数,本发明w为常数1。因此,得到的形变域及其能量函数均为关于小波系数的函数。
为了减少被估计小波参数的数量,把小波系数根据尺度和子带的不同进行分类,而形变域便由不同尺度和不同子带的小波系数来描述。本发明采用递进策略寻找小波系数,估计大尺度上的小波系数后,再估计小尺度上的小波系数。不同尺度和不同子带小波系数的估计次序为:1-2-3-4-22-33-44,具体如下表所示:
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (8)
1.一种联合图像分割与形变配准技术的自适应目标区域转换方法,其特征是:包括以下步骤:
(1)采用基于图谱的自动分割方法,对采集的N组CT图像去除伪影、噪声成分,采用梯度矢量流进行胸肋骨自动提取;
(2)在N组图像中随机选取一组作为参考图像,其余作为目标图像,将得到的肝脏训练集和胸肋骨分别进行肝脏图谱和胸肋骨图谱的创建;
(3)输入患者当前计划CT图像,采用GVF进行胸肋骨自动提取,将胸肋骨图谱作为参考图像,患者计划CT图像作为目标图像,以胸肋骨为特征,采用仿射配准算法与胸肋骨图谱进行刚体配准,得到患者计划CT图像中肝脏区域的初始值;
(4)使用基于正交小波基函数的形变配准算法对得到的肝脏区域与肝脏图谱进行形变配准,使肝脏边界更为精准,通过基于图谱的分割方法在患者计划CT上自动分割出肝脏区域;
(5)结合水平集模型与似然函数模型在计划CT上自动分割胸腹部肿瘤区域;
(6)以胸肋骨为特征,通过仿射变换将最新CBCT图像与计划CT图像进行刚体配准,为补偿刚体配准时不能恢复的形变范围,使用基于正交小波基函数的形变配准算法,计算出CBCT与计划CT的形变域,将计划CT中的危及器官和肿瘤的轮廓以及等剂量线信息自动转移到CBCT上。
2.如权利要求1所述的自适应目标区域转换方法,其特征是:所述步骤(2)中,具体方法为:
在N组满足条件的CT中随机选取一组作为参考图像,将其与其余的目标图像采用仿射变换进行配准,所有的图像对齐后,统计空间每个像素存在的概率;N组图像数据配准后,得到一个概率图谱;采用多次循环迭代,以提高所建图谱的可靠性;
建立概率图谱后,根据具体应用环境选取合适的阈值,按照灰度值随图谱概率增大而增大的准则,把图谱从概率空间转化到灰度空间。
3.如权利要求1所述的自适应目标区域转换方法,其特征是:所述步骤(5)中,其具体方法为:
采用似然函数与水平集函数相结合的能量函数,在患者计划CT图像上自动获取胸腹部肿瘤信息;似然能量函数Ep(φ)由感兴趣区域能量函数eT(x)和背景区域能量函数eB(x)两部分构成;用似然函数的有限高斯模型估计感兴趣肿瘤区域密度分布,似然函数的高斯混合模型估计胸腹部背景区域密度分布,通过以上似然函数的估计来划分图像概率密度;水平集函数用于获得边缘轮廓信息,防止肿瘤边界泄漏。
4.如权利要求3所述的自适应目标区域转换方法,其特征是:所述步骤(5)中,所述肿瘤区域的密度分布采用有限高斯模型估计,具体的表达式为:
其中,ΩT为肿瘤区域,p(X)为肿瘤区域像素的
概率,μT为高斯函数均值,σT为标准差,u0(x)为患者CT图像;
背景区域密度分布采用高斯混合模型估计,具体的表达式为:
其中,ΩB为背景区域,p(X)为背景区域像素
的概率,和分别为高斯模型的均值和标准差,K为背景区域的分类个数,为分类
的混合系数;
似然函数具体表达式为:
其中,
距离正则化水平集模型能量函数的具体表达式为:
Ee(φ)=μRP(φ)+λLg(φ)+αAg(φ)=
其中:μ、λ、α分别是正则化项、长度项、面
积项的权重系数,
当似然能量函数与水平集能量函数最小时,水平集函数演化曲线收敛到感兴趣区域边界,水平集函数演化表达式为:
其中,div(·)为散度算子。
5.如权利要求1所述的自适应目标区域转换方法,其特征是:所述步骤(6)中,其具体方法为:
将已导入加速器的计划CT作为目标图像,CBCT图像作为参考图像,采用GVF对CBCT图像进行胸肋骨自动提取,以胸肋骨为特征,利用仿射配准算法与患者计划CT进行刚体配准;为补偿刚体配准时不能恢复的形变范围,使用基于正交小波基函数的形变配准算法,获得形变域,将形变域三维小波分解后,使小波系数的初始值为0;采用非线性最小二乘LM优化算法,寻找合适的小波系数,根据所需精度控制在某个尺度上收敛;由估计到的小波系数得到参考图像与目标图像之间的形变域;将计划CT中的危及器官和肿瘤的轮廓以及等剂量线信息自动转移到CBCT上,进而实现剂量体积直方图分析。
6.如权利要求1所述的自适应目标区域转换方法,其特征是:所述步骤(6)中,在形变配准中,使用Navier偏微分方程描述发生形变的各向同性物质的平衡状态,通过Navier偏微分方程来建立形变能量函数,将形变域建模为小波系数的函数,其具体表达式为:
其中,(x,y,z)表示参考图像的空间坐标系,(x′,y′,z′)表示目标图像的空间坐标系,u表示形变域,c表示小波系数,u是c的函数。
7.如权利要求6所述的自适应目标区域转换方法,其特征是:所述Navier偏微分方程的外力具体表达式为:
H(X)是图像X的香农熵。
8.如权利要求6所述的自适应目标区域转换方法,其特征是:所述Navier偏微分方程的内力方程具体表达式为:
其中,λ和μ是Lame系数;
形变能量函数具体表达式为:
E(c)=inter(c)+w*exter(c)
其中,w是加权常数,本发明w为常数1,inter(c)表示内力约束函数,exter(c)表示外力约束函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410764128.5A CN104408734B (zh) | 2014-12-11 | 2014-12-11 | 联合图像分割与形变配准技术的自适应目标区域转换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410764128.5A CN104408734B (zh) | 2014-12-11 | 2014-12-11 | 联合图像分割与形变配准技术的自适应目标区域转换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104408734A true CN104408734A (zh) | 2015-03-11 |
CN104408734B CN104408734B (zh) | 2017-10-27 |
Family
ID=52646363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410764128.5A Active CN104408734B (zh) | 2014-12-11 | 2014-12-11 | 联合图像分割与形变配准技术的自适应目标区域转换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104408734B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106683042A (zh) * | 2017-01-05 | 2017-05-17 | 南京觅踪电子科技有限公司 | 一种基于非线性最小二乘的多模态医学图像融合方法以及*** |
CN107924565A (zh) * | 2015-06-24 | 2018-04-17 | 光线搜索实验室公司 | 处理图像数据的***和方法 |
CN108460809A (zh) * | 2017-02-22 | 2018-08-28 | 西门子保健有限责任公司 | 用于***癌检测和分类的深度卷积编码器-解码器 |
WO2019041262A1 (en) * | 2017-08-31 | 2019-03-07 | Shenzhen United Imaging Healthcare Co., Ltd. | SYSTEM AND METHOD FOR IMAGE SEGMENTATION |
CN109615642A (zh) * | 2018-11-05 | 2019-04-12 | 北京全域医疗技术有限公司 | 一种放疗计划中的危及器官自动勾画方法及装置 |
CN110363750A (zh) * | 2019-06-28 | 2019-10-22 | 福建师范大学 | 一种基于多模态数据融合的牙根管形态自动提取方法 |
CN110570426A (zh) * | 2018-06-06 | 2019-12-13 | 国际商业机器公司 | 使用深度学习的图像联合配准和分割 |
CN111344737A (zh) * | 2017-09-07 | 2020-06-26 | 医科达有限公司 | 自适应放射疗法*** |
CN112150569A (zh) * | 2020-09-24 | 2020-12-29 | 北京羽医甘蓝信息技术有限公司 | 将cbct图像生成为ct图像的方法、装置及终端设备 |
WO2021104061A1 (zh) * | 2019-11-26 | 2021-06-03 | 中国科学院深圳先进技术研究院 | 一种肝脏变形预测方法、***及电子设备 |
CN113706409A (zh) * | 2021-08-18 | 2021-11-26 | 苏州雷泰医疗科技有限公司 | 基于人工智能的cbct图像增强方法、装置和存储介质 |
WO2021238732A1 (en) * | 2020-05-23 | 2021-12-02 | Ping An Technology (Shenzhen) Co., Ltd. | Device and method for alignment of multi-modal clinical images using joint synthesis, segmentation, and registration |
CN114022524A (zh) * | 2021-11-04 | 2022-02-08 | 上海体素信息科技有限公司 | 基于图谱配准的ct颈肋定位方法及*** |
CN116188377A (zh) * | 2022-12-27 | 2023-05-30 | 中山大学肿瘤防治中心(中山大学附属肿瘤医院、中山大学肿瘤研究所) | 一种锥形束ct图像配准结果评价方法及*** |
CN116365400A (zh) * | 2023-03-28 | 2023-06-30 | 上海新联合电气有限公司 | 一种全电操控无人值守智能配电*** |
CN117152543A (zh) * | 2023-10-30 | 2023-12-01 | 山东浪潮科学研究院有限公司 | 一种图像分类方法、装置、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2088556A1 (en) * | 2008-02-05 | 2009-08-12 | Siemens Aktiengesellschaft | Method and system for registration to align at least one preoperative image to at least one intraoperativley acquired image |
CN101984462A (zh) * | 2010-09-08 | 2011-03-09 | 天津科技大学 | 电阻抗断层功能图像与解剖图像的合成方法 |
-
2014
- 2014-12-11 CN CN201410764128.5A patent/CN104408734B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2088556A1 (en) * | 2008-02-05 | 2009-08-12 | Siemens Aktiengesellschaft | Method and system for registration to align at least one preoperative image to at least one intraoperativley acquired image |
CN101984462A (zh) * | 2010-09-08 | 2011-03-09 | 天津科技大学 | 电阻抗断层功能图像与解剖图像的合成方法 |
Non-Patent Citations (3)
Title |
---|
CHANGYANG LI等: "A Likelihood and Local Constraint Level Set Model for Liver Tumor Segmentation from CT Volumes", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 * |
YAOQIN XIE等: "Feature-based rectal contour propagation from planning CT to cone beam CT", 《MEDICAL PHYSICS》 * |
李登旺: "医学图像配准和分割技术研究及在图像引导放射治疗***中的应用", 《万方学位论文全文数据库》 * |
Cited By (32)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107924565A (zh) * | 2015-06-24 | 2018-04-17 | 光线搜索实验室公司 | 处理图像数据的***和方法 |
CN107924565B (zh) * | 2015-06-24 | 2021-09-14 | 光线搜索实验室公司 | 处理图像数据的***和方法 |
CN106683042A (zh) * | 2017-01-05 | 2017-05-17 | 南京觅踪电子科技有限公司 | 一种基于非线性最小二乘的多模态医学图像融合方法以及*** |
CN106683042B (zh) * | 2017-01-05 | 2020-08-11 | 卡本(深圳)医疗科技有限公司 | 一种基于非线性最小二乘的多模态医学图像融合方法以及*** |
CN108460809A (zh) * | 2017-02-22 | 2018-08-28 | 西门子保健有限责任公司 | 用于***癌检测和分类的深度卷积编码器-解码器 |
US10789709B2 (en) | 2017-08-31 | 2020-09-29 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image segmentation |
WO2019041262A1 (en) * | 2017-08-31 | 2019-03-07 | Shenzhen United Imaging Healthcare Co., Ltd. | SYSTEM AND METHOD FOR IMAGE SEGMENTATION |
US11508068B2 (en) | 2017-08-31 | 2022-11-22 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image segmentation |
US11715206B2 (en) | 2017-08-31 | 2023-08-01 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image segmentation |
US10621724B2 (en) | 2017-08-31 | 2020-04-14 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for image segmentation |
CN113577577A (zh) * | 2017-09-07 | 2021-11-02 | 医科达有限公司 | 自适应放射疗法*** |
CN111344737A (zh) * | 2017-09-07 | 2020-06-26 | 医科达有限公司 | 自适应放射疗法*** |
CN113577577B (zh) * | 2017-09-07 | 2024-07-19 | 医科达有限公司 | 自适应放射疗法*** |
US11318327B2 (en) | 2017-09-07 | 2022-05-03 | Elekta, Inc. | Adaptive radiotherapy system |
CN110570426A (zh) * | 2018-06-06 | 2019-12-13 | 国际商业机器公司 | 使用深度学习的图像联合配准和分割 |
CN110570426B (zh) * | 2018-06-06 | 2023-04-25 | 国际商业机器公司 | 使用深度学习的图像联合配准和分割 |
CN109615642B (zh) * | 2018-11-05 | 2021-04-06 | 北京全域医疗技术集团有限公司 | 一种放疗计划中的危及器官自动勾画方法及装置 |
CN109615642A (zh) * | 2018-11-05 | 2019-04-12 | 北京全域医疗技术有限公司 | 一种放疗计划中的危及器官自动勾画方法及装置 |
CN110363750A (zh) * | 2019-06-28 | 2019-10-22 | 福建师范大学 | 一种基于多模态数据融合的牙根管形态自动提取方法 |
CN110363750B (zh) * | 2019-06-28 | 2023-05-09 | 福建师范大学 | 一种基于多模态数据融合的牙根管形态自动提取方法 |
WO2021104061A1 (zh) * | 2019-11-26 | 2021-06-03 | 中国科学院深圳先进技术研究院 | 一种肝脏变形预测方法、***及电子设备 |
WO2021238732A1 (en) * | 2020-05-23 | 2021-12-02 | Ping An Technology (Shenzhen) Co., Ltd. | Device and method for alignment of multi-modal clinical images using joint synthesis, segmentation, and registration |
CN112150569A (zh) * | 2020-09-24 | 2020-12-29 | 北京羽医甘蓝信息技术有限公司 | 将cbct图像生成为ct图像的方法、装置及终端设备 |
CN113706409A (zh) * | 2021-08-18 | 2021-11-26 | 苏州雷泰医疗科技有限公司 | 基于人工智能的cbct图像增强方法、装置和存储介质 |
CN114022524A (zh) * | 2021-11-04 | 2022-02-08 | 上海体素信息科技有限公司 | 基于图谱配准的ct颈肋定位方法及*** |
CN114022524B (zh) * | 2021-11-04 | 2024-07-23 | 上海体素信息科技有限公司 | 基于图谱配准的ct颈肋定位方法及*** |
CN116188377A (zh) * | 2022-12-27 | 2023-05-30 | 中山大学肿瘤防治中心(中山大学附属肿瘤医院、中山大学肿瘤研究所) | 一种锥形束ct图像配准结果评价方法及*** |
CN116188377B (zh) * | 2022-12-27 | 2023-10-27 | 中山大学肿瘤防治中心(中山大学附属肿瘤医院、中山大学肿瘤研究所) | 一种锥形束ct图像配准结果评价方法及*** |
CN116365400A (zh) * | 2023-03-28 | 2023-06-30 | 上海新联合电气有限公司 | 一种全电操控无人值守智能配电*** |
CN116365400B (zh) * | 2023-03-28 | 2024-02-13 | 上海新联合电气有限公司 | 一种全电操控无人值守智能配电*** |
CN117152543A (zh) * | 2023-10-30 | 2023-12-01 | 山东浪潮科学研究院有限公司 | 一种图像分类方法、装置、设备及存储介质 |
CN117152543B (zh) * | 2023-10-30 | 2024-06-07 | 山东浪潮科学研究院有限公司 | 一种图像分类方法、装置、设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN104408734B (zh) | 2017-10-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104408734A (zh) | 联合图像分割与形变配准技术的自适应目标区域转换方法 | |
CN105031833B (zh) | 放射治疗装置的剂量验证*** | |
Anders et al. | Performance of an atlas-based autosegmentation software for delineation of target volumes for radiotherapy of breast and anorectal cancer | |
Zhang et al. | Respiratory liver motion estimation and its effect on scanned proton beam therapy | |
US9282944B2 (en) | C-arm pose estimation using intensity-based registration of imaging modalities | |
Schubert et al. | A comprehensive assessment by tumor site of patient setup using daily MVCT imaging from more than 3,800 helical tomotherapy treatments | |
Hostettler et al. | A real-time predictive simulation of abdominal viscera positions during quiet free breathing | |
Nassef et al. | Quantification of dose uncertainties in cumulated dose estimation compared to planned dose in prostate IMRT | |
Xing et al. | Computational challenges for image-guided radiation therapy: framework and current research | |
CN104117151A (zh) | 一种在线自适应放疗计划优化方法 | |
McGurk et al. | Extension of the NCAT phantom for the investigation of intra-fraction respiratory motion in IMRT using 4D Monte Carlo | |
Zhang et al. | A kernel-based method for markerless tumor tracking in kV fluoroscopic images | |
Huang et al. | 2D ultrasound imaging based intra-fraction respiratory motion tracking for abdominal radiation therapy using machine learning | |
Shao et al. | Real-time liver tumor localization via a single x-ray projection using deep graph neural network-assisted biomechanical modeling | |
CN110706779B (zh) | 一种放射治疗计划精确目标函数的自动生成方法 | |
Geneser et al. | Quantifying variability in radiation dose due to respiratory-induced tumor motion | |
Shao et al. | Automatic liver tumor localization using deep learning‐based liver boundary motion estimation and biomechanical modeling (DL‐Bio) | |
Kim et al. | Prostate localization on daily cone-beam computed tomography images: accuracy assessment of similarity metrics | |
CN113499091A (zh) | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和*** | |
Shi et al. | Fast shading correction for cone-beam CT via partitioned tissue classification | |
Hooshangnejad et al. | DAART: a deep learning platform for deeply accelerated adaptive radiation therapy for lung cancer | |
Camps et al. | Various approaches for pseudo-CT scan creation based on ultrasound to ultrasound deformable image registration between different treatment time points for radiotherapy treatment plan adaptation in prostate cancer patients | |
Wang et al. | An individualized strategy to estimate the effect of deformable registration uncertainty on accumulated dose in the upper abdomen | |
Götz et al. | On the use of particle filters for electromagnetic tracking in high dose rate brachytherapy | |
Bueno et al. | A geodesic deformable model for automatic segmentation of image sequences applied to radiation therapy |
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 |