CN109727672B - 患者胸腹部肿瘤呼吸运动预测跟踪方法 - Google Patents
患者胸腹部肿瘤呼吸运动预测跟踪方法 Download PDFInfo
- Publication number
- CN109727672B CN109727672B CN201811624784.XA CN201811624784A CN109727672B CN 109727672 B CN109727672 B CN 109727672B CN 201811624784 A CN201811624784 A CN 201811624784A CN 109727672 B CN109727672 B CN 109727672B
- Authority
- CN
- China
- Prior art keywords
- respiratory
- prediction
- body surface
- model
- tumor
- 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
- 230000000241 respiratory effect Effects 0.000 title claims abstract description 114
- 206010028980 Neoplasm Diseases 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000001727 in vivo Methods 0.000 claims abstract description 36
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 31
- 230000029058 respiratory gaseous exchange Effects 0.000 claims abstract description 30
- 210000001015 abdomen Anatomy 0.000 claims abstract description 7
- 238000003384 imaging method Methods 0.000 claims abstract description 6
- 230000003287 optical effect Effects 0.000 claims abstract description 6
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 230000008859 change Effects 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 12
- 230000003044 adaptive effect Effects 0.000 claims description 10
- 238000012549 training Methods 0.000 claims description 10
- 230000001537 neural effect Effects 0.000 claims description 8
- 238000011478 gradient descent method Methods 0.000 claims description 7
- 210000000038 chest Anatomy 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 238000010276 construction Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 210000000115 thoracic cavity Anatomy 0.000 claims description 2
- 230000003187 abdominal effect Effects 0.000 claims 1
- 238000010606 normalization Methods 0.000 claims 1
- 230000005855 radiation Effects 0.000 abstract description 2
- 238000001959 radiotherapy Methods 0.000 description 11
- 230000008569 process Effects 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 230000001788 irregular Effects 0.000 description 6
- 238000012795 verification Methods 0.000 description 4
- 238000013528 artificial neural network Methods 0.000 description 3
- 238000012706 support-vector machine Methods 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 208000006678 Abdominal Neoplasms Diseases 0.000 description 1
- 206010011224 Cough Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 208000037841 lung tumor Diseases 0.000 description 1
- 230000036210 malignancy Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 208000025426 neoplasm of thorax Diseases 0.000 description 1
- 210000000056 organ Anatomy 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 206010041232 sneezing Diseases 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 201000003957 thoracic cancer Diseases 0.000 description 1
Images
Landscapes
- Radiation-Therapy Devices (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明提供一种患者胸腹部肿瘤呼吸运动预测跟踪方法，包括以下步骤：采用光学位置跟踪设备，连续地获取体表呼吸运动；采用X射线立体平面成像定位设备，间断地获取体内肿瘤位置；数据预处理：对体表呼吸运动进行小波分解，得到基线、代表主要呼吸信号的低频部分、以及代表噪声的高频部分；通过呼吸预测算法得到体表呼吸运动的预测结果；以获取的体内肿瘤位置和对应时间点的体表呼吸运动作为输入，采用线性和非线性混合多项式模型，建立体表‑体内呼吸运动相关性模型；复制体表‑体内呼吸运动相关性模型，以预测的连续的体表呼吸运动作为输入，计算连续的体内肿瘤位置；进入相关性模型更新判别逻辑。本发明可用于引导实施精确的跟踪放射治疗。
Description
技术领域
本发明属于医学信号处理领域，特别是涉及一种应用于放射治疗过程的，患者胸腹部肿瘤呼吸运动的预测跟踪方法。
背景技术
在恶性肿瘤中发病率最高的肺部肿瘤和其它胸部、腹部肿瘤，会发生由病人呼吸引起的术中移位。随着技术的发展，采用X射线立体平面成像和光学实时跟踪相结合的跟踪技术，使对患者体内肿瘤呼吸运动的实时跟踪成为可能；且高精度的呼吸预测算法，能补偿数据处理和设备电气机械动作等引入的时间延迟。由此，放射治疗设备可在放射治疗过程中，对运动肿瘤进行跟踪和补偿，实施更精准的放射治疗，可极大地减少放疗射线对肿瘤周边正常组织和关键器官的损伤，改善治疗效果，降低放疗副作用。
其他研究主要关注于下一时刻呼吸位置的预测算法，效果较好的有神经网络(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：if
is A1and
is B1,then
规则2：if
is A1and
is B2,then
……
规则9：if
is A3and
is 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 (4)
1.一种患者胸腹部肿瘤呼吸运动预测跟踪方法，其特征在于，包括以下步骤：
步骤1，数据采集：采用光学位置跟踪设备，连续地获取体表呼吸运动；采用X射线立体平面成像定位设备，间断地获取体内肿瘤位置；
步骤2，数据预处理：对体表呼吸运动进行小波分解，得到基线、代表主要呼吸信号的低频部分、以及代表噪声的高频部分；
步骤3，通过呼吸预测算法得到体表呼吸运动的预测结果；经过小波分解后，基线部分较平滑，采用线性算法进行预测；低频部分是呼吸信号去除基线后的呼吸振幅部分，采用改进的自适应神经模糊推理系统ANFIS进行预测；高频部分主要是噪声，采用简单移动平均法进行预测；最后综合三部分预测值作为体表呼吸运动的预测结果；
步骤4：以获取的体内肿瘤位置和对应时间点的体表呼吸运动作为输入，采用线性和非线性混合多项式模型，建立体表-体内呼吸运动相关性模型；并优化该相关性模型；
步骤5：复制体表-体内呼吸运动相关性模型，以预测的连续的体表呼吸运动作为输入，计算连续的体内肿瘤位置；
步骤6：进入相关性模型更新判别逻辑，以小波分解后的体表呼吸运动和相关性模型输出与原始体内肿瘤位置之间的偏差作为输入；进行逻辑判别，并根据判别结果，决定沿用原相关性模型，或者继续获取体内肿瘤位置，更新相关性模型；
步骤4中，采用的线性和非线性混合多项式模型描述如下：
步骤4.1，先采用线性一阶多项式模型Pt＝aPs+b，Ps为体表呼吸运动，Pt为体内肿瘤位置，a、b为多项式参数；以最小二乘法拟合出多项式最优参数a和b；
步骤4.2，统计相关性模型输出与原始体内肿瘤位置之间的偏差；采用偏差的平均值(Mean)和标准偏差(STDEV)，评价相关性模型是否符合需求，其中标准偏差是评价线性/非线性的主要指标；
步骤4.3，若标准偏差在阈值内，则采用此线性多项式模型；
步骤4.4，否则，更换二阶非线性多项式模型Pt＝aPs 2+bPs+c，Ps为体表呼吸运动，Pt为体内肿瘤位置，a、b、c为多项式参数，通过梯度下降法迭代优化获得多项式最优参数a、b和c；
步骤6中，采用的相关性模型更新判别逻辑描述如下：
步骤6.1，通过小波分解后的体表呼吸运动，根据低频部分统计呼吸运动的周期、振幅变化情况，结合呼吸基线变化情况，评估体表呼吸模型变化情况；
步骤6.2，若有一定时间内获取的原始体内肿瘤位置，则计算相关性模型输出与原始体内肿瘤位置之间的偏差；
步骤6.3，进行更新逻辑判别，条件为上述体表呼吸模型变化情况和位置偏差，若满足模型更新条件，则获取体内肿瘤位置，更新相关性模型；否则继续使用原相关性模型。
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层：计算所有规则的总输出。
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 CN109727672A (zh) | 2019-05-07 |
| CN109727672B true 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) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110681074B (zh) * | 2019-10-29 | 2021-06-15 | 苏州大学 | 基于双向gru网络的肿瘤呼吸运动预测方法 |
| WO2022165812A1 (zh) * | 2021-02-07 | 2022-08-11 | 北京肿瘤医院(北京大学肿瘤医院) | 一种利用光学体表运动信号合成实时图像的系统 |
| CN113674393B (zh) * | 2021-07-12 | 2023-09-26 | 中国科学院深圳先进技术研究院 | 呼吸运动模型的构建方法和无标记呼吸运动预测方法 |
| CN113499091B (zh) * | 2021-08-19 | 2023-08-15 | 四川大学华西医院 | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统 |
| CN114177545B (zh) * | 2022-01-17 | 2023-11-07 | 中国科学院合肥物质科学研究院 | 一种用于放疗中无接触式呼吸节律监测装置和方法 |
| CN114927215B (zh) * | 2022-04-27 | 2023-08-25 | 苏州大学 | 基于体表点云数据直接预测肿瘤呼吸运动的方法及系统 |
| CN115187608B (zh) * | 2022-09-14 | 2023-02-03 | 苏州大学 | 一种基于体表显著性分析的呼吸特征提取方法 |
| CN117494058B (zh) * | 2024-01-02 | 2024-04-09 | 中国科学院苏州生物医学工程技术研究所 | 辅助手术机器人穿刺的呼吸运动预测方法、设备及介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101794357A (zh) * | 2010-03-30 | 2010-08-04 | 江苏瑞尔医疗科技有限公司 | 呼吸同步跟踪系统的位置信号的模糊预测方法及预测器 |
| CN102512162A (zh) * | 2011-12-20 | 2012-06-27 | 华南理工大学 | 基于多模式脑电小波分析的术中运动区功能定位系统 |
| CN104285175A (zh) * | 2012-05-02 | 2015-01-14 | 法国国家科学研究中心 | 使用小波分析进行单颗粒定位的方法和装置 |
| CN110660063A (zh) * | 2019-09-19 | 2020-01-07 | 山东省肿瘤防治研究院（山东省肿瘤医院） | 多图像融合的肿瘤三维位置精准定位系统 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101843955B (zh) * | 2010-03-30 | 2015-07-15 | 江苏瑞尔医疗科技有限公司 | 一种用于放射治疗的同步跟踪定位系统 |
| CN102306381B (zh) * | 2011-06-02 | 2013-02-27 | 西安电子科技大学 | 基于Beamlet与小波变换的图像融合方法 |
| CN102429658A (zh) * | 2011-12-20 | 2012-05-02 | 华南理工大学 | 基于脑电慢皮层电位小波分析的术中运动区功能定位系统 |
| CN104182989A (zh) * | 2014-07-15 | 2014-12-03 | 江苏大学 | 一种基于压缩传感的粒子滤波视觉跟踪方法 |
| CN105117593A (zh) * | 2015-08-18 | 2015-12-02 | 河海大学 | 基于小波变换和粒子群改进灰色模型的短期风速预测方法 |
| US10932861B2 (en) * | 2016-01-14 | 2021-03-02 | Auris Health, Inc. | Electromagnetic tracking surgical system and method of controlling the same |
| CN106563210B (zh) * | 2016-11-10 | 2020-07-10 | 苏州大学 | 基于ut变换的放疗机器人肿瘤呼吸运动估计及预测方法 |
| CN107951482A (zh) * | 2017-11-13 | 2018-04-24 | 浙江好络维医疗技术有限公司 | 一种基于小波分解与样条插值的心电信号基线拉平方法 |
-
2018
- 2018-12-28 CN CN201811624784.XA patent/CN109727672B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101794357A (zh) * | 2010-03-30 | 2010-08-04 | 江苏瑞尔医疗科技有限公司 | 呼吸同步跟踪系统的位置信号的模糊预测方法及预测器 |
| CN102512162A (zh) * | 2011-12-20 | 2012-06-27 | 华南理工大学 | 基于多模式脑电小波分析的术中运动区功能定位系统 |
| CN104285175A (zh) * | 2012-05-02 | 2015-01-14 | 法国国家科学研究中心 | 使用小波分析进行单颗粒定位的方法和装置 |
| CN110660063A (zh) * | 2019-09-19 | 2020-01-07 | 山东省肿瘤防治研究院（山东省肿瘤医院） | 多图像融合的肿瘤三维位置精准定位系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109727672A (zh) | 2019-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109727672B (zh) | 患者胸腹部肿瘤呼吸运动预测跟踪方法 | |
| Wang et al. | A feasibility of respiration prediction based on deep Bi-LSTM for real-time tumor tracking | |
| Ruan et al. | Real-time prediction of respiratory motion based on local regression methods | |
| CN107273925B (zh) | 基于局部感受野和半监督深度自编码的肺实质ct影像处理装置 | |
| Viswanath et al. | Design and analysis performance of kidney stone detection from ultrasound image by level set segmentation and ANN classification | |
| Yu et al. | Rapid prediction of respiratory motion based on bidirectional gated recurrent unit network | |
| CN110378881B (zh) | 一种基于深度学习的肿瘤定位系统 | |
| CN108159576B (zh) | 一种放疗中人体胸腹表面区域呼吸运动预测方法 | |
| CN106485721B (zh) | 从光学相干断层图像获取视网膜结构的方法及其系统 | |
| Gupta et al. | Gammatonegram based triple classification of lung sounds using deep convolutional neural network with transfer learning | |
| Goodband et al. | A comparison of neural network approaches for on‐line prediction in IGRT | |
| CN114266729A (zh) | 一种基于机器学习的胸部肿瘤放疗后放射性肺炎预测方法和系统 | |
| Koçer et al. | Classifying neuromuscular diseases using artificial neural networks with applied Autoregressive and Cepstral analysis | |
| CN112155549B (zh) | 基于深度卷积脉冲神经网络的adhd疾病诊断辅助决策系统 | |
| CN114247061A (zh) | 肿瘤动态跟踪控制方法、装置及放射治疗设备 | |
| CN113499091B (zh) | 一种患者体表和体内肿瘤运动相关性及肿瘤内部动度的预测方法和系统 | |
| CN117918789A (zh) | 一种基于机器学习的全身麻醉药效预测方法及系统 | |
| Dogan et al. | Application of substitution box of present cipher for automated detection of snoring sounds | |
| CN116935009B (zh) | 基于历史数据分析进行预测的手术导航系统 | |
| Zhu et al. | A kalman filter based approach for markerless pose tracking and assessment | |
| CN109191425B (zh) | 基于多层神经网络模型医学影像分析方法 | |
| Peng et al. | Prostate segmentation of ultrasound images based on interpretable-guided mathematical Model | |
| CN110675960B (zh) | 基于时间序列分析的呼吸运动预测建模方法 | |
| Yu et al. | Correlated skin surface and tumor motion modeling for treatment planning in robotic radiosurgery | |
| CN113724304A (zh) | 一种基于深度学习的食管区域图像自动配准方法及系统 |
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 |