CN109727672A - 患者胸腹部肿瘤呼吸运动预测跟踪方法 - Google Patents
患者胸腹部肿瘤呼吸运动预测跟踪方法 Download PDFInfo
- Publication number
- CN109727672A CN109727672A CN201811624784.XA CN201811624784A CN109727672A CN 109727672 A CN109727672 A CN 109727672A CN 201811624784 A CN201811624784 A CN 201811624784A CN 109727672 A CN109727672 A CN 109727672A
- Authority
- CN
- China
- Prior art keywords
- body surface
- respiratory movement
- prediction
- model
- correlation models
- 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
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
本发明提供一种患者胸腹部肿瘤呼吸运动预测跟踪方法,包括以下步骤:采用光学位置跟踪设备,连续地获取体表呼吸运动;采用X射线立体平面成像定位设备,间断地获取体内肿瘤位置;数据预处理:对体表呼吸运动进行小波分解,得到基线、代表主要呼吸信号的低频部分、以及代表噪声的高频部分;通过呼吸预测算法得到体表呼吸运动的预测结果;以获取的体内肿瘤位置和对应时间点的体表呼吸运动作为输入,采用线性和非线性混合多项式模型,建立体表‑体内呼吸运动相关性模型;复制体表‑体内呼吸运动相关性模型,以预测的连续的体表呼吸运动作为输入,计算连续的体内肿瘤位置;进入相关性模型更新判别逻辑。本发明可用于引导实施精确的跟踪放射治疗。
Description
技术领域
本发明属于医学信号处理领域,特别是涉及一种应用于放射治疗过程的,患者胸腹部肿瘤呼吸运动的预测跟踪方法。
背景技术
在恶性肿瘤中发病率最高的肺部肿瘤和其它胸部、腹部肿瘤,会发生由病人呼吸引起的术中移位。随着技术的发展,采用X射线立体平面成像和光学实时跟踪相结合的跟踪技术,使对患者体内肿瘤呼吸运动的实时跟踪成为可能;且高精度的呼吸预测算法,能补偿数据处理和设备电气机械动作等引入的时间延迟。由此,放射治疗设备可在放射治疗过程中,对运动肿瘤进行跟踪和补偿,实施更精准的放射治疗,可极大地减少放疗射线对肿瘤周边正常组织和关键器官的损伤,改善治疗效果,降低放疗副作用。
呼吸同步追踪***是一种相对完善的***和方法。该技术采用X射线立体平面成像定位设备间断获取患者体内肿瘤位置,采用光学位置跟踪设备实时跟踪患者体表呼吸运动。该技术采用一种线性和非线性混合的多项式模型,能够较好的拟合体内肿瘤-体表呼吸运动之间的运动相关性关系;然后根据体表呼吸运动,得到体内肿瘤的呼吸运动位置。该技术通过在随机或者特定时相上曝光采图,获取体内肿瘤实际位置,验证并更新运动相关性模型。该技术针对设备响应等引入的时间延迟,采用结合了模糊逻辑(FuzzyLogic,FL)、PM和LMS的混合预测算法,呼吸预测精度较好,已有大量临床实践验证。
其他研究主要关注于下一时刻呼吸位置的预测算法,效果较好的有神经网络(NN)、支持向量机(SVM)、自适应神经模糊推理***(ANFIS)等。
呼吸同步追踪***采用的运动相关性模型的验证和更新逻辑,对于呼吸模型变化的响应有滞后性,且对于平稳呼吸过程,又无依据进一步减少验证采图次数,降低患者所受剂量。且其采用的混合预测算法,预测不规则或者异常呼吸的误差较大,仍有较大改进空间。
其他研究的呼吸预测算法,NN算法容易陷入局部最小化或者产生过拟合问题。SVM算法复杂度高,预测结果依赖参数选择。ANFIS算法采用呼吸运动的位置作为模型的输入参数构造模糊集进行预测,预测能力较强,但对不规则信号特别是当呼吸信号的振幅发生突变时,预测误差过大。
发明内容
本发明的目的在于克服现有技术中存在的不足,提供一种患者胸腹部肿瘤呼吸运动预测跟踪方法,能够实时跟踪、预测患者体内肿瘤呼吸运动;(1)为放射治疗门控方法提供输入数据;(2)与实时跟踪放射治疗设备进行通信,引导其实施精确的跟踪放射治疗。本发明采用的技术方案是:
本发明提出一种患者胸腹部肿瘤呼吸运动预测跟踪方法,包括以下步骤:
步骤1,数据采集:采用光学位置跟踪设备,连续地获取体表呼吸运动;采用X射线立体平面成像定位设备,间断地获取体内肿瘤位置;
步骤2,数据预处理:对体表呼吸运动进行小波分解,得到基线、代表主要呼吸信号的低频部分、以及代表噪声的高频部分;
进行小波分解的具体步骤为:
步骤2.1,利用小波函数进行分解和重构的Mallat算法,选取小波基,确定分解层数;同时,为了消除边界效应,在利用小波分解得到基线部分时,边界延拓采用0阶平滑延拓;在利用小波分解得到低频部分和高频部分时,边界延拓采用1阶平滑延拓;
步骤2.2,先对体表呼吸运动的原始信号进行第一次多层分解,将低频分量作为基线部分;
步骤2.3,然后将原信号去除基线的其余部分,再进行一次小波分解,分解得到的低频分量作为低频部分,高频分量作为高频部分。
步骤3,通过呼吸预测算法得到体表呼吸运动的预测结果;经过小波分解后,基线部分较平滑,采用线性算法(LF)进行预测;低频部分是呼吸信号去除基线后的呼吸振幅部分,采用改进的自适应神经模糊推理***(ANFIS)进行预测;高频部分主要是噪声,采用简单移动平均法进行预测;最后综合三部分预测值作为体表呼吸运动的预测结果;
具体地,采用改进的自适应神经模糊推理***(ANFIS)进行预测,步骤如下:
步骤3.1,进行低频信号预测的自适应神经模糊推理***ANFIS采用一种新的预测模型结构,不仅仅采用呼吸运动的位置作为输入参数,而且通过位置信息得到速度(通过邻近的位置信息求得),将呼吸运动的位置和速度作为预测模型的输入参数;
步骤3.2,呼吸运动的位置和速度可以作为预测模型构造中模糊划分的特征,预测模型将位置按大小分成小、中、大三个模糊集;同时又按速度大小分为负的大速度、小速度、正的大速度三个模糊集;从而构造一个结合位置和速度的3×3模糊集;
步骤3.3,依据预测模型输入输出参数构造训练集,应用训练集训练预测模型;采用混合学习算法训练ANFIS的参数,在每一次迭代中,输入信号首先沿网络正向传递到倒数第二层,此时固定前提参数,利用最小二乘调整结论参数,然后信号继续沿网络正向传递,将获得的误差信号沿网络反向传播,利用梯度下降法调节前提参数;
步骤3.4,进行预测时,将最新位置和速度数据输入已训练好的预测模型,得到输入数据对应时间τs个采样点以后的位置预测值;
步骤3.5,在应用ANFIS已训练好的预测模型进行预测时,若当前时刻的呼吸位置或速度超出历史数据位置的幅度范围,首先执行幅度调整,将位置或速度调整到历史数据的幅度范围内,得到幅度调整后的位置或速度;然后将经过幅度调整后的最新输入信号输入已训练好的预测模型,得到输出值;再进行幅度反调整,得到呼吸运动的位置预测值。呼吸信号的位置和速度并不一定同时超出历史数据范围,该ANFIS模型根据信号的不规则情况具体选择需要调整的内容。
步骤4:以获取的体内肿瘤位置和对应时间点的体表呼吸运动作为输入,采用线性和非线性混合多项式模型,建立体表-体内呼吸运动相关性模型;并优化该相关性模型;具体地,
步骤4.1,先采用线性多项式模型,以最小二乘法拟合出多项式最优参数;
步骤4.2,统计相关性模型输出的肿瘤位置与原始体内肿瘤位置之间的偏差;
步骤4.3,若该偏差在阈值内则采用此线性多项式模型;
步骤4.4,否则,更换二阶非线性多项式模型,通过梯度下降法迭代优化获得多项式最优参数。
步骤5:复制体表-体内呼吸运动相关性模型,以预测的连续的体表呼吸运动作为输入,计算连续的体内肿瘤位置;
步骤6:进入相关性模型更新判别逻辑,以小波分解后的体表呼吸运动和相关性模型输出与原始体内肿瘤位置之间的偏差作为输入;进行逻辑判别,并根据判别结果,决定沿用原相关性模型,或者继续获取体内肿瘤位置,更新相关性模型。
具体地,采用的相关性模型更新判别逻辑描述如下:
步骤6.1,通过小波分解后的体表呼吸运动,根据低频部分统计呼吸运动的周期、振幅变化情况,结合呼吸基线变化情况,评估体表呼吸模型变化情况;
步骤6.2,若有一定时间内获取的原始体内肿瘤位置,则计算相关性模型输出与原始体内肿瘤位置之间的偏差;
步骤6.3,进行更新逻辑判别,条件为上述体表呼吸模型变化情况和位置偏差,若满足模型更新条件,则获取体内肿瘤位置,更新相关性模型;否则继续使用原相关性模型。
本发明的优点在于:
(1)通过小波分解将呼吸信号的基线、低频和高频分离出来,用于模型更新和信号预测。
(2)评估呼吸运动的基线、振幅和周期变化情况,结合相关性模型输出与实际肿瘤位置之间的偏差,建立一种合理的运动相关性模型更新逻辑。
(3)采用此相关性模型更新逻辑,及时发现并应对患者呼吸模型的变化,及时曝光采图获取体内肿瘤实际位置,验证并更新运动相关性模型,提高体内肿瘤跟踪精度。
(4)采用此相关性模型更新逻辑,根据呼吸模型变化情况,不断调整更新相关性模型,减少相关性模型输出与实际肿瘤位置之间的偏差过大的情况,可降低重建相关性模型的概率。
(5)采用此相关性模型更新逻辑,且对于平稳呼吸过程,可进一步减少验证采图次数,降低患者所受剂量。
(3)从呼吸信号中提取出相对规则的低频信号成分采用ANFIS进行预测,减小了基线漂移对预测精度的影响,使预测结果更加准确。
(4)在呼吸运动位置的基础上,引入速度作为ANFIS模型结构的输入参数,构造一个结合位置和速度的3×3模糊集,能有效捕获呼吸运动的变化,提升模型对非规则呼吸运动的预测能力。
(5)当信号发生不规则(甚至异常)变化时,对超出历史数据幅度范围的最新输入信号的位置或速度进行相应的幅度调整可提高预测准确度,从而使预测模型对非规则呼吸运动或者随机性异常呼吸的预测能力更强。
(6)将呼吸信号高频部分分离出来通过简单移动平均进行预测,相当于对呼吸信号滤波去噪,减小了噪声对呼吸预测的影响,使预测的呼吸曲线更加平滑。降低高频成分引起放疗呼吸运动补偿中导致机器颤动的概率。
(7)采用小波分解和改进型自适应神经模糊推理相结合的呼吸预测算法,提升呼吸预测进度,对规则和不规则呼吸都有很强的鲁棒性;当呼吸特征改变时能迅速调整适应为新的呼吸模式,并且能够合理处理患者打喷嚏、咳嗽、讲话等干扰信号。
附图说明
图1为本发明的呼吸运动预测跟踪方法的流程图。
图2为本发明的基于小波分解和改进ANFIS的呼吸预测算法的概要流程图。
图3为本发明的采用的小波分解呼吸信号示例图。
图4为本发明提出的以位置和速度作为输入参数的ANFIS预测模型结构图。
图5为本发明的改进ANFIS预测模型中位置和速度的3x3模糊集合示意图。
图6为本发明的呼吸预测幅度调整流程图。
图7为本发明的幅度调整示意图。
图8为本发明的相关性模型更新判别逻辑流程图。
具体实施方式
下面结合具体附图和实施例对本发明作进一步说明。
如图1所示,本发明提出的一种患者胸腹部肿瘤呼吸运动预测跟踪方法,包括以下步骤:
步骤1,数据采集:采用光学位置跟踪设备,连续地获取体表呼吸运动;采用X射线立体平面成像定位设备,间断地获取体内肿瘤位置;
步骤2,数据预处理:对体表呼吸运动进行小波分解,得到基线、代表主要呼吸信号的低频部分、以及代表噪声的高频部分;
利用小波函数进行分解和重构的主流方法Mallat算法,小波基选取5阶Daubechies小波,首先将体表呼吸运动的原信号分解到第8层,其低频分量作为基线部分;然后将原信号去除基线的其余部分,再进行一次小波分解,分解得到的低频分量作为低频部分,高频分量作为高频部分;如图3所示,对图中左侧的原信号进行小波分解,得到右侧的基线、低频和高频信号;其中,基线部分反映了呼吸运动振荡位置的变化;高频部分为噪声;低频部分为去除噪声、去除基线变化的主要呼吸过程。
步骤3,如图2所示,通过呼吸预测算法得到体表呼吸运动的预测结果;经过小波分解后,基线部分较平滑,采用线性算法(LF)进行预测;低频部分是呼吸信号去除基线后的呼吸振幅部分,采用改进的自适应神经模糊推理***(ANFIS)进行预测;高频部分主要是噪声,采用简单移动平均法进行预测;最后综合三部分预测值作为体表呼吸运动的预测结果;
具体地,进行低频信号预测的自适应神经模糊推理***ANFIS采用一种新的模型结构,如图4所示,描述如下:
ANFIS有2个输入和单个输出f。是呼吸运动的位置,是相应时刻的呼吸速度(通过邻近的位置信息求得)。位置和速度可以作为预测模型构造中模糊划分的特征。输入参数来自呼吸信号低频部分输出f是对应时间τs个采样点以后的位置预测值预测模型将位置按大小分成小、中、大三个模糊集,即A1、A2、A3;同时模糊集又按速度大小分为负的大速度、小速度、正的大速度三个模糊集,即B1、B2、B3;位置xlow和速度vlow构成3×3模糊集合,如图5所示;每一个模糊集对应一条模糊规则:
规则1:ifis A1andis B1,then
规则2:ifis A1andis B2,then
……
规则9:ifis A3andis B3,then
其中,fr是规则r对应的函数,{pr,qr,zr}为结论参数(r=1,2,…,9);
ANFIS预测模型结构中每一层的结点具有相同的函数,设第1层的第1个结点的输出为O11;
第1层:有6个结点,将输入变量模糊化,输出对应模糊集的隶属度,例如该层第1个结点的输出为:O1,1是输入变量在模糊集A1中的隶属度,隶属函数选择钟形函数:
{a,b,c}是隶属函数参数,称为前提参数,这些参数值改变,隶属函数也随之改变;
第2层:这一层每个结点的输出是所有输入变量在不同模糊集中隶属度的代数积,例如该层第1个结点的输出为:代表规则1的权重,则规则r的权重为wr=O2,r,r=1,2,…,9;
第3层:将各条规则的权重归一化:r=1,2,…,9;
第4层:这一层每个结点的输出为每条规则的输出,是每条规则函数值与对应归一化权重的乘积:r=1,2,…,9;
第5层:计算所有规则的总输出:r=1,2,…,9。
依据预测模型输入输出参数构造训练集:
上述矩阵共k+1行(k=N-3τs-1),第一列和第二列分别为预测模型的位置输入数据和速度输入数据,第三列为预测模型的输出数据,即输入数据对应时间τs个采样点之后的呼吸位置;应用训练集训练预测模型,采用梯度下降法调节隶属函数参数,采用最小二乘法调节规则函数参数;进行预测时,将最新数据输入已训练好的预测模型,得到输入数据对应时间τs个采样点以后的位置预测值:
采用混合学习算法训练ANFIS的参数,在每一次迭代中,输入信号首先沿网络正向传递到第4层,此时固定前提参数,利用最小二乘调整结论参数,然后信号继续沿网络正向传递,将获得的误差信号沿网络反向传播,利用梯度下降法调节前提参数;
具体地,下面以图6、图7所示的位置幅度调整过程为例进行,说明在应用ANFIS已训练好的预测模型进行预测时,当前时刻的呼吸位置或速度超出历史数据位置的幅度调整过程;
在应用已训练好的预测模型进行预测时,若当前时刻的呼吸位置超出历史数据位置的幅度范围首先执行图7中所示的幅度调整,将位置调整到历史数据位置的幅度范围内,得到幅度调整后的位置:
其中,或者
然后将经过幅度调整后的最新输入信号输入已训练好的预测模型,得到输出值:
最后进行图7中所示的幅度反调整,得到呼吸运动的位置预测值:
呼吸信号的位置和速度并不一定同时超出历史数据范围,该ANFIS模型根据信号的不规则情况具体选择需要调整的内容,上述例子仅调整了位置幅度;
步骤4:以获取的体内肿瘤位置和对应时间点的体表呼吸运动作为输入,采用线性和非线性混合多项式模型,建立体表-体内呼吸运动相关性模型;并优化该相关性模型;
具体地,采用的线性和非线性混合多项式模型描述如下:
步骤4.1,先采用线性一阶多项式模型Pt=aPs+b,Ps为体表呼吸运动,Pt为体内肿瘤位置,a、b为多项式参数;以最小二乘法拟合出多项式最优参数a和b。考虑到每个位置均为3维坐标(x,y,z),采用主成分分析PCA方法对输入数据进行降维处理,将高维特征空间转化成低维特征空间,以减少“维数灾难”;
步骤4.2,统计相关性模型输出与原始体内肿瘤位置之间的偏差;采用偏差的平均值(Mean)和标准偏差(STDEV),评价相关性模型是否符合需求,其中标准偏差是评价线性/非线性的主要指标;
步骤4.3,若标准偏差在阈值内,则采用此线性多项式模型;
步骤4.4,否则,更换二阶非线性多项式模型Pt=aPs 2+bPs+c,Ps为体表呼吸运动,Pt为体内肿瘤位置,a、b、c为多项式参数,通过梯度下降法迭代优化获得多项式最优参数a、b和c。
步骤5:复制体表-体内呼吸运动相关性模型,以预测的连续的体表呼吸运动作为输入,计算连续的体内肿瘤位置。
步骤6:进入相关性模型更新判别逻辑,以小波分解后的体表呼吸运动和相关性模型输出与原始体内肿瘤位置之间的偏差作为输入;进行逻辑判别,并根据判别结果,决定沿用原相关性模型,或者继续获取体内肿瘤位置,更新相关性模型。
具体地,如图8所示,采用的相关性模型更新判别逻辑描述如下:
步骤6.1,通过小波分解后的体表呼吸运动,根据低频部分统计n个历史呼吸运动的周期{Ts1、Ts2…Tsn}、振幅{As1、As2…Asn},评估其变化情况CT&A是否超过阈值;结合呼吸基线变化情况Cbase,评估体表呼吸模型变化情况Cs=CT&A+Cbase。
步骤6.2,若有一定时间t内获取的原始体内肿瘤位置,则获取相关性模型输出与原始体内肿瘤位置之间的偏差ΔP=|Pt-Pm|,Pt为原始体内肿瘤位置,Pm为相关性模型输出的体内肿瘤位置。评估该偏差所反映的相关性模型是否需要更新,为Cm。
步骤6.3,判别相关性模型更新标识C=Cs+Cm。
步骤6.4,满足相关性模型更新条件,则在短时间内获取体内肿瘤位置,更新相关性模型;否则继续使用原相关性模型。
最后所应说明的是,以上具体实施方式仅用以说明本发明的技术方案而非限制,尽管参照实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (6)
1.一种患者胸腹部肿瘤呼吸运动预测跟踪方法,其特征在于,包括以下步骤:
步骤1,数据采集:采用光学位置跟踪设备,连续地获取体表呼吸运动;采用X射线立体平面成像定位设备,间断地获取体内肿瘤位置;
步骤2,数据预处理:对体表呼吸运动进行小波分解,得到基线、代表主要呼吸信号的低频部分、以及代表噪声的高频部分;
步骤3,通过呼吸预测算法得到体表呼吸运动的预测结果;经过小波分解后,基线部分较平滑,采用线性算法进行预测;低频部分是呼吸信号去除基线后的呼吸振幅部分,采用改进的自适应神经模糊推理***ANFIS进行预测;高频部分主要是噪声,采用简单移动平均法进行预测;最后综合三部分预测值作为体表呼吸运动的预测结果;
步骤4:以获取的体内肿瘤位置和对应时间点的体表呼吸运动作为输入,采用线性和非线性混合多项式模型,建立体表-体内呼吸运动相关性模型;并优化该相关性模型;
步骤5:复制体表-体内呼吸运动相关性模型,以预测的连续的体表呼吸运动作为输入,计算连续的体内肿瘤位置;
步骤6:进入相关性模型更新判别逻辑,以小波分解后的体表呼吸运动和相关性模型输出与原始体内肿瘤位置之间的偏差作为输入;进行逻辑判别,并根据判别结果,决定沿用原相关性模型,或者继续获取体内肿瘤位置,更新相关性模型。
2.如权利要求1所述的患者胸腹部肿瘤呼吸运动预测跟踪方法,其特征在于,步骤2中,进行小波分解的具体步骤为:
步骤2.1,利用小波函数进行分解和重构的Mallat算法,选取小波基,确定分解层数;同时,为了消除边界效应,在利用小波分解得到基线部分时,边界延拓采用0阶平滑延拓;在利用小波分解得到低频部分和高频部分时,边界延拓采用1阶平滑延拓;
步骤2.2,先对体表呼吸运动的原始信号进行第一次多层分解,将低频分量作为基线部分;
步骤2.3,然后将原信号去除基线的其余部分,再进行一次小波分解,分解得到的低频分量作为低频部分,高频分量作为高频部分。
3.如权利要求1所述的患者胸腹部肿瘤呼吸运动预测跟踪方法,其特征在于,步骤3中,采用改进的自适应神经模糊推理***ANFIS进行预测,步骤如下:
步骤3.1,进行低频信号预测的自适应神经模糊推理***ANFIS采用一种新的预测模型结构,不仅仅采用呼吸运动的位置作为输入参数,而且通过位置信息得到速度,将呼吸运动的位置和速度作为预测模型的输入参数;
步骤3.2,呼吸运动的位置和速度能够作为预测模型构造中模糊划分的特征,预测模型将位置按大小分成小、中、大三个模糊集;同时又按速度大小分为负的大速度、小速度、正的大速度三个模糊集;从而构造一个结合位置和速度的3×3模糊集;
步骤3.3,依据预测模型输入输出参数构造训练集,应用训练集训练预测模型;采用混合学习算法训练ANFIS的参数,在每一次迭代中,输入信号首先沿网络正向传递到倒数第二层,此时固定前提参数,利用最小二乘调整结论参数,然后信号继续沿网络正向传递,将获得的误差信号沿网络反向传播,利用梯度下降法调节前提参数;
步骤3.4,进行预测时,将最新位置和速度数据输入已训练好的预测模型,得到输入数据对应时间τs个采样点以后的位置预测值;
步骤3.5,在应用ANFIS已训练好的预测模型进行预测时,若当前时刻的呼吸位置或速度超出历史数据位置的幅度范围,首先执行幅度调整,将位置或速度调整到历史数据的幅度范围内,得到幅度调整后的位置或速度;然后将经过幅度调整后的最新输入信号输入已训练好的预测模型,得到输出值;再进行幅度反调整,得到呼吸运动的位置预测值。
4.如权利要求3所述的患者胸腹部肿瘤呼吸运动预测跟踪方法,其特征在于,
ANFIS预测模型中每一个模糊集对应一条模糊规则;ANFIS预测模型结构中每一层的结点具有相同的函数,预测模型具体包括:
第1层:有6个结点,将输入变量模糊化,输出对应模糊集的隶属度;是呼吸运动的位置,是相应时刻的呼吸速度;
第2层:这一层每个结点的输出是所有输入变量在不同模糊集中隶属度的代数积,代表各规则的权重;
第3层:将各条规则的权重归一化;
第4层:这一层每个结点的输出为每条规则的输出,是每条规则函数值与对应归一化权重的乘积;
第5层:计算所有规则的总输出。
5.如权利要求1所述的患者胸腹部肿瘤呼吸运动预测跟踪方法,其特征在于,步骤4具体包括:
步骤4.1,先采用线性多项式模型,以最小二乘法拟合出多项式最优参数;
步骤4.2,统计相关性模型输出的肿瘤位置与原始体内肿瘤位置之间的偏差;
步骤4.3,若该偏差在阈值内则采用此线性多项式模型;
步骤4.4,否则,更换二阶非线性多项式模型,通过梯度下降法迭代优化获得多项式最优参数。
6.如权利要求1所述的患者胸腹部肿瘤呼吸运动预测跟踪方法,其特征在于,步骤6中,采用的相关性模型更新判别逻辑描述如下:
步骤6.1,通过小波分解后的体表呼吸运动,根据低频部分统计呼吸运动的周期、振幅变化情况,结合呼吸基线变化情况,评估体表呼吸模型变化情况;
步骤6.2,若有一定时间内获取的原始体内肿瘤位置,则计算相关性模型输出与原始体内肿瘤位置之间的偏差;
步骤6.3,进行更新逻辑判别,条件为上述体表呼吸模型变化情况和位置偏差,若满足模型更新条件,则获取体内肿瘤位置,更新相关性模型;否则继续使用原相关性模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811624784.XA CN109727672B (zh) | 2018-12-28 | 2018-12-28 | 患者胸腹部肿瘤呼吸运动预测跟踪方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811624784.XA CN109727672B (zh) | 2018-12-28 | 2018-12-28 | 患者胸腹部肿瘤呼吸运动预测跟踪方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109727672A true CN109727672A (zh) | 2019-05-07 |
CN109727672B CN109727672B (zh) | 2023-04-07 |
Family
ID=66296698
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811624784.XA Active CN109727672B (zh) | 2018-12-28 | 2018-12-28 | 患者胸腹部肿瘤呼吸运动预测跟踪方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109727672B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110681074A (zh) * | 2019-10-29 | 2020-01-14 | 苏州大学 | 基于双向gru网络的肿瘤呼吸运动预测方法 |
CN113499091A (zh) * | 2021-08-19 | 2021-10-15 | 四川大学华西医院 | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和*** |
CN113574565A (zh) * | 2021-02-07 | 2021-10-29 | 北京肿瘤医院(北京大学肿瘤医院) | 一种利用光学体表运动信号合成实时图像的*** |
CN114177545A (zh) * | 2022-01-17 | 2022-03-15 | 中国科学院合肥物质科学研究院 | 一种用于放疗中无接触式呼吸节律监测装置和方法 |
CN114927215A (zh) * | 2022-04-27 | 2022-08-19 | 苏州大学 | 基于体表点云数据直接预测肿瘤呼吸运动的方法及*** |
WO2023283866A1 (zh) * | 2021-07-12 | 2023-01-19 | 中国科学院深圳先进技术研究院 | 呼吸运动模型的构建方法和无标记呼吸运动预测方法 |
CN117494058A (zh) * | 2024-01-02 | 2024-02-02 | 中国科学院苏州生物医学工程技术研究所 | 辅助手术机器人穿刺的呼吸运动预测方法、设备及介质 |
WO2024055386A1 (zh) * | 2022-09-14 | 2024-03-21 | 苏州大学 | 一种基于体表显著性分析的呼吸特征提取方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101794357A (zh) * | 2010-03-30 | 2010-08-04 | 江苏瑞尔医疗科技有限公司 | 呼吸同步跟踪***的位置信号的模糊预测方法及预测器 |
CN101843955A (zh) * | 2010-03-30 | 2010-09-29 | 江苏瑞尔医疗科技有限公司 | 呼吸同步跟踪***的位置信号的混合预测方法及预测器 |
CN102306381A (zh) * | 2011-06-02 | 2012-01-04 | 西安电子科技大学 | 基于Beamlet与小波变换的图像融合方法 |
CN102429658A (zh) * | 2011-12-20 | 2012-05-02 | 华南理工大学 | 基于脑电慢皮层电位小波分析的术中运动区功能定位*** |
CN102512162A (zh) * | 2011-12-20 | 2012-06-27 | 华南理工大学 | 基于多模式脑电小波分析的术中运动区功能定位*** |
CN104182989A (zh) * | 2014-07-15 | 2014-12-03 | 江苏大学 | 一种基于压缩传感的粒子滤波视觉跟踪方法 |
CN104285175A (zh) * | 2012-05-02 | 2015-01-14 | 法国国家科学研究中心 | 使用小波分析进行单颗粒定位的方法和装置 |
CN105117593A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 基于小波变换和粒子群改进灰色模型的短期风速预测方法 |
CN106563210A (zh) * | 2016-11-10 | 2017-04-19 | 苏州大学 | 基于ut变换的放疗机器人肿瘤呼吸运动估计及预测***及其使用方法 |
US20170202627A1 (en) * | 2016-01-14 | 2017-07-20 | Auris Surgical Robotics, Inc. | Electromagnetic tracking surgical system and method of controlling the same |
CN107951482A (zh) * | 2017-11-13 | 2018-04-24 | 浙江好络维医疗技术有限公司 | 一种基于小波分解与样条插值的心电信号基线拉平方法 |
CN110660063A (zh) * | 2019-09-19 | 2020-01-07 | 山东省肿瘤防治研究院(山东省肿瘤医院) | 多图像融合的肿瘤三维位置精准定位*** |
-
2018
- 2018-12-28 CN CN201811624784.XA patent/CN109727672B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101794357A (zh) * | 2010-03-30 | 2010-08-04 | 江苏瑞尔医疗科技有限公司 | 呼吸同步跟踪***的位置信号的模糊预测方法及预测器 |
CN101843955A (zh) * | 2010-03-30 | 2010-09-29 | 江苏瑞尔医疗科技有限公司 | 呼吸同步跟踪***的位置信号的混合预测方法及预测器 |
CN102306381A (zh) * | 2011-06-02 | 2012-01-04 | 西安电子科技大学 | 基于Beamlet与小波变换的图像融合方法 |
CN102429658A (zh) * | 2011-12-20 | 2012-05-02 | 华南理工大学 | 基于脑电慢皮层电位小波分析的术中运动区功能定位*** |
CN102512162A (zh) * | 2011-12-20 | 2012-06-27 | 华南理工大学 | 基于多模式脑电小波分析的术中运动区功能定位*** |
CN104285175A (zh) * | 2012-05-02 | 2015-01-14 | 法国国家科学研究中心 | 使用小波分析进行单颗粒定位的方法和装置 |
CN104182989A (zh) * | 2014-07-15 | 2014-12-03 | 江苏大学 | 一种基于压缩传感的粒子滤波视觉跟踪方法 |
CN105117593A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 基于小波变换和粒子群改进灰色模型的短期风速预测方法 |
US20170202627A1 (en) * | 2016-01-14 | 2017-07-20 | Auris Surgical Robotics, Inc. | Electromagnetic tracking surgical system and method of controlling the same |
CN106563210A (zh) * | 2016-11-10 | 2017-04-19 | 苏州大学 | 基于ut变换的放疗机器人肿瘤呼吸运动估计及预测***及其使用方法 |
CN107951482A (zh) * | 2017-11-13 | 2018-04-24 | 浙江好络维医疗技术有限公司 | 一种基于小波分解与样条插值的心电信号基线拉平方法 |
CN110660063A (zh) * | 2019-09-19 | 2020-01-07 | 山东省肿瘤防治研究院(山东省肿瘤医院) | 多图像融合的肿瘤三维位置精准定位*** |
Non-Patent Citations (2)
Title |
---|
史少华: ""基于分离有限状态模型的呼吸预测算法"", 《清华大学学报(自然科学版)》 * |
朱丹等: "小波分解结合自适应神经模糊推理***的呼吸预测研究", 《天津医科大学学报》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110681074B (zh) * | 2019-10-29 | 2021-06-15 | 苏州大学 | 基于双向gru网络的肿瘤呼吸运动预测方法 |
CN110681074A (zh) * | 2019-10-29 | 2020-01-14 | 苏州大学 | 基于双向gru网络的肿瘤呼吸运动预测方法 |
CN113574565A (zh) * | 2021-02-07 | 2021-10-29 | 北京肿瘤医院(北京大学肿瘤医院) | 一种利用光学体表运动信号合成实时图像的*** |
WO2022165812A1 (zh) * | 2021-02-07 | 2022-08-11 | 北京肿瘤医院(北京大学肿瘤医院) | 一种利用光学体表运动信号合成实时图像的*** |
WO2023283866A1 (zh) * | 2021-07-12 | 2023-01-19 | 中国科学院深圳先进技术研究院 | 呼吸运动模型的构建方法和无标记呼吸运动预测方法 |
CN113499091A (zh) * | 2021-08-19 | 2021-10-15 | 四川大学华西医院 | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和*** |
CN113499091B (zh) * | 2021-08-19 | 2023-08-15 | 四川大学华西医院 | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和*** |
CN114177545B (zh) * | 2022-01-17 | 2023-11-07 | 中国科学院合肥物质科学研究院 | 一种用于放疗中无接触式呼吸节律监测装置和方法 |
CN114177545A (zh) * | 2022-01-17 | 2022-03-15 | 中国科学院合肥物质科学研究院 | 一种用于放疗中无接触式呼吸节律监测装置和方法 |
CN114927215A (zh) * | 2022-04-27 | 2022-08-19 | 苏州大学 | 基于体表点云数据直接预测肿瘤呼吸运动的方法及*** |
CN114927215B (zh) * | 2022-04-27 | 2023-08-25 | 苏州大学 | 基于体表点云数据直接预测肿瘤呼吸运动的方法及*** |
WO2024055386A1 (zh) * | 2022-09-14 | 2024-03-21 | 苏州大学 | 一种基于体表显著性分析的呼吸特征提取方法 |
CN117494058A (zh) * | 2024-01-02 | 2024-02-02 | 中国科学院苏州生物医学工程技术研究所 | 辅助手术机器人穿刺的呼吸运动预测方法、设备及介质 |
CN117494058B (zh) * | 2024-01-02 | 2024-04-09 | 中国科学院苏州生物医学工程技术研究所 | 辅助手术机器人穿刺的呼吸运动预测方法、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109727672B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109727672A (zh) | 患者胸腹部肿瘤呼吸运动预测跟踪方法 | |
Kakar et al. | Respiratory motion prediction by using the adaptive neuro fuzzy inference system (ANFIS) | |
US7727160B2 (en) | Ventilator patient synchronization | |
Yu et al. | Rapid prediction of respiratory motion based on bidirectional gated recurrent unit network | |
CN108367161A (zh) | 放疗***、数据处理方法及存储介质 | |
Koçer et al. | Classifying neuromuscular diseases using artificial neural networks with applied Autoregressive and Cepstral analysis | |
Singh et al. | Wavelet transform based pneumonia classification of chest X-ray images | |
Abbod et al. | Survey on the use of smart and adaptive engineering systems in medicine | |
CN113499091B (zh) | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和*** | |
Abdulameer | An improvement of MRI brain images classification using dragonfly algorithm as trainer of artificial neural network | |
CN111815650B (zh) | 一种基于改进模糊c均值的脑磁共振图像分割方法 | |
CN116227606B (zh) | 一种关节角度预测方法、终端设备及介质 | |
Shi et al. | Respiratory prediction based on multi-scale temporal convolutional network for tracking thoracic tumor movement | |
CN101843955B (zh) | 一种用于放射治疗的同步跟踪定位*** | |
Jindal | MobileSOFT: U: A deep learning framework to monitor heart rate during intensive physical exercise | |
Zhu et al. | A kalman filter based approach for markerless pose tracking and assessment | |
Damayanti et al. | Classification of tumor based on magnetic resonance (MR) brain images using wavelet energy feature and neuro-fuzzy model | |
Shuzan et al. | Machine Learning-Based Respiration Rate and Blood Oxygen Saturation Estimation Using Photoplethysmogram Signals. Bioengineering 2023, 10, 167 | |
CN115019388A (zh) | 一种利用单目摄像头拍摄步态视频全自动步态分析方法 | |
Yu et al. | Correlated skin surface and tumor motion modeling for treatment planning in robotic radiosurgery | |
CN110675960B (zh) | 基于时间序列分析的呼吸运动预测建模方法 | |
CN113724304A (zh) | 一种基于深度学习的食管区域图像自动配准方法及*** | |
CN110681074B (zh) | 基于双向gru网络的肿瘤呼吸运动预测方法 | |
Lyu et al. | Skeleton-based sleep posture recognition with BP neural network | |
McGowan et al. | Vocal tract normalization for midsagittal articulatory recovery with analysis-by-synthesis |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |