CN105931262B - 一种肿瘤位置预测装置 - Google Patents

一种肿瘤位置预测装置 Download PDF

Info

Publication number
CN105931262B
CN105931262B CN201610551244.8A CN201610551244A CN105931262B CN 105931262 B CN105931262 B CN 105931262B CN 201610551244 A CN201610551244 A CN 201610551244A CN 105931262 B CN105931262 B CN 105931262B
Authority
CN
China
Prior art keywords
data
module
knub position
sample
input
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.)
Expired - Fee Related
Application number
CN201610551244.8A
Other languages
English (en)
Other versions
CN105931262A (zh
Inventor
吴宜灿
胡丽琴
龙鹏程
贾婧
曹瑞芬
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hefei Institutes of Physical Science of CAS
Original Assignee
Hefei Institutes of Physical Science of CAS
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hefei Institutes of Physical Science of CAS filed Critical Hefei Institutes of Physical Science of CAS
Priority to CN201610551244.8A priority Critical patent/CN105931262B/zh
Publication of CN105931262A publication Critical patent/CN105931262A/zh
Application granted granted Critical
Publication of CN105931262B publication Critical patent/CN105931262B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/20Special algorithmic details
    • G06T2207/20081Training; Learning

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Radiation-Therapy Devices (AREA)

Abstract

本发明涉及一种肿瘤位置预测装置,该装置包括数据集的采集和预处理模块,样本训练集生成模块和数据学习及预测模块,通过数据集的采集和预处理模块获得呼吸运动数据和肿瘤位置历史数据,通过样本训练集生成模块和数据学习及预测模块生成肿瘤位置预测函数,从而使用该预测函数对病人新采集的肿瘤位置信息进行快速并准确地预测。用于在精准放射治疗过程中提高病人肿瘤靶区的实时跟踪定位的精度,在治疗过程中对病人位置进行纠正,从而最大限度的减少靶区和器官运动带来的误差影响。

Description

一种肿瘤位置预测装置
技术领域
本发明涉及放射治疗中病人肿瘤靶区的定位跟踪,具体指在精准放射治疗过程中对肿瘤位置的预测,涉及一种快速精确的肿瘤位置预测装置。
背景技术
随着放疗技术的发展,肿瘤放疗已经进入到“精准定位,精准计划,精准实施,精准验证”的阶段,精准放疗技术使得在肿瘤治疗上实现高精度,高剂量,高疗效和低损伤的现代放疗模式逐步成为现实。然而临床放疗的精度往往受到人体呼吸运动的影响,导致肿瘤靶区发生不同程度的移动甚至错位,严重影响了治疗效果。为了能够减小呼吸运动带来的影响,可以通过对肿瘤的实时跟踪来进行补偿。然而,跟踪信号发出与治疗设备的响应之间有着一定的时间延迟,无法真正实现实时的放疗跟踪,延迟时间一般为几百毫秒,主要包括获取数据的时间、计算处理的时间、传输数据的时间以及存在的机械延迟等。为了补偿系统的延迟时间,消除延迟对精准照射的影响,可以采用对肿瘤靶区进行位置预测的方法来实现。在放疗中的对肿瘤的位置预测需要有良好的实时性,这对预测模型提出了很高的要求。传统的预测方法有:线性估计、卡尔曼滤波器、线性外推、神经网络等。这些方法中线性估计和线性外推法对预测的准确性较差;对卡尔曼滤波器来说,预测模型的初始参数难以确定,而且不适合进行长时间预测;神经网络的预测精度较高,在实际应用中,神经网络是用有限的训练数据来训练,对于单隐层前馈神经网络,其N个隐层结点和激励函数都通过N个样本集来学习,需要一定的训练数据的时间,在实时预测调整和较长时间预测中效果也受到一定限制。
对基于单隐层神经网络的学习算法,其学习速度比传统的单隐层神经网络学习算法如Back-Propagation(BP)快上千倍同时获得更好的泛化性能。传统的梯度下降法的会遇见局部最小值、不适当的步长以及过拟合的问题,为避免这些问题往往引入复杂的停止方法或增加迭代次数。基于单隐层神经网络的学习算法学习阶段可以在几秒内完成,泛化性能在大部分情况下都优于基于梯度下降法的学习算法,对于单隐层和多个单层神经网络都适用,可以逼近任意连续函数,因此适用范围更广。与其他学习算法相比的优越之处在于,该算法以解析的方法获得单隐层前馈神经网络唯一的权值解,同时能够保证此解具有最小训练误差和最小范数输出权值,因此该学习算法具有易于实现,获得最小的训练误差,最小的权值范数,良好的泛化性能,运行速度特别快等特点。目前放疗过程中对肿瘤位置的预测多数直接是对已有的肿瘤位置数据进行回顾性分析,有效性和实时性还有待检验。从医生和物理师的角度出发,在整个治疗过程中,可以说如果精准定位不能实现或是出现较大偏差,那么在精准放疗计划的制定和实施就毫无意义。
发明内容
本发明技术解决问题:克服现有技术的不足,提供一种在放射治疗过程中快速精确的肿瘤位置预测装置,快速有效解决肿瘤靶区位置预测问题,用于在精准放射治疗过程中提高病人肿瘤靶区的实时跟踪定位的精度,使预测结果能够满足放疗计划制定的要求,而又使得计算时间能够可以在接受的范围,从而最大限度的减少靶区和器官运动带来的误差影响。
为了实现上述目的,本发明所采用的技术方案是:
一种快速精确的肿瘤位置预测装置,包括数据集的采集和预处理模块,样本训练集生成模块和数据学习及预测模块;数据集的采集和预处理模块将采集和预处理后的数据通过计算机连接传送至样本训练集生成模块,再将样本训练集生成模块得到的数据传送至数据学习及预测模块;
数据集采集和预处理模块,包括红外空间位置探测单元和分析与处理模块;红外空间位置探测单元用于提供红外定位小球作为标记点的位置信息,进而反馈给分析与处理模块,用于实时显示和自动量化分析,将当前获取的位置信息和分析与处理模块中保存的原始治疗计划中的位置信息进行比较和配准,从而得到校正后的空间位置信息,作为经过数据采集和预处理后的数据并传递给样本训练集生成模块作为输入数据;
所述红外空间位置探测单元通过记录红外线发射和接收器得到红外定位小球标记点的位置信息,采用红外线与红外定位小球标记点相连,接收检测区域的红外线,红外空间位置探测单元反馈红外定位小球标记点的位置和方向的变化信息,将转换的数据和信息输出给分析与处理模块;
样本训练集生成模块,接收数据集采集和预处理模块经过数据采集和预处理后的数据作为样本训练集生成的输入,基于这些输入数据,从中选取肿瘤患者正常呼吸的N(N>5)个周期,基于体内骨性标记点或植入标记点的位置信息,获得肿瘤位置信息作为训练用历史数据,构建N个样本训练集,作为肿瘤位置预测函数的输入{(x1,y1),(x2,y2),L,(xN,yN)},其中x1表示样本集1中肿瘤位置输入,y1表示样本集1中肿瘤位置输出,x2表示样本集2中肿瘤位置输入,y2表示样本集2中肿瘤位置输出,xN表示样本集N中肿瘤位置输入,yN表示样本集N中肿瘤位置输出,N个样本训练集肿瘤位置输入输出数据作为样本训练集生成模块得到的数据,传递给数据学习及预测模块;
数据学习及预测模块,接收由样本训练集生成模块得到的数据作为输入,用来模拟肿瘤位置输入与肿瘤位置输出之间的预测函数关系y=f(x),在数据学习过程中,y是经过预测得到的肿瘤位置输出,x是肿瘤位置输入,在此基础上构建单隐层神经网络,设置隐层节点数为K和激活函数g(x),激活函数可选择为s型函数、基本周期函数、正弦函数、余弦函数中的一种,进而在区间[-1,1]上,随机选择输入层到隐层的权值w和隐层的偏置值b,通过将采集的肿瘤位置的历史数据样本信息输入计算隐层输出矩阵H,计算输出层权值其中Y=[y1,...,yN]T为肿瘤位置样本输出组成的矩阵,进而通过求得的输出层权值β对肿瘤位置进行预测。
本发明与现有技术相比的优点在于:
(1)本发明装置包括了数据集的采集和预处理模块,样本训练集生成模块和数据学习及预测模块,可以实现快速实时无创的肿瘤位置预测。其中数据集的采集和预处理模块采用了红外空间位置探测单元进行肿瘤位置信息采集,可以实现实时无创的呼吸运动数据采集和处理;对于样本训练集生成模块,接收数据集采集和预处理模块经过数据采集和预处理后的数据作为样本训练集生成的输入,基于体内骨性标记点或植入标记点的位置信息,获得肿瘤位置信息作为训练用历史数据,构建样本训练集,作为肿瘤位置预测函数的输入;同时数据学习及预测模块采用基于单隐层神经网络的学习回归模型来模拟肿瘤位置输入与肿瘤位置输出之间的函数关系,输入与输出间的关系由历史的肿瘤位置数据来描述,能够实现随机选择输入层到隐层的权值和隐层的偏置值,通过将采集的肿瘤位置的历史数据样本信息输入计算输出层权值,进而通过求得的输出层权值对肿瘤位置进行预测。在历史数据充足的情况下,通过充分利用历史数据中的信息,准确地预测放疗过程中肿瘤位置情况,提高了放疗过程中肿瘤位置预测的处理速度和精度。
(2)本发明不仅能够准确地建立肿瘤位置预测的模型,实现了可以快速高效进行肿瘤靶区的位置预测,满足了实际的用户需求。
附图说明
图1是本发明装置的总体结构及装置中各个模块之间的连接示意图;
图2是具体实施方式中数据学习及预测模块的单隐层神经网络结构示意图。
具体实施方式
下面结合附图对本发明的具体实施方式作进一步说明:
如图1所示,为本发明装置的总体结构及模块之间连接示意图;从图1可知,本发明装置包括数据集的采集和预处理模块,样本训练集生成模块和数据学习及预测模块;数据集的采集和预处理模块包含红外空间位置探测单元,采集和预处理后的数据通过计算机连接传送至样本训练集生成模块,进而将样本训练集生成模块得到的数据传送至数据学习及预测模块;
数据集采集和预处理模块,包括红外空间位置探测单元和分析与处理模块;红外空间位置探测单元用于提供红外定位小球作为标记点的位置信息,进而反馈给分析与处理模块,用于实时显示和自动量化分析,将当前获取的位置信息和分析与处理模块中保存的原始治疗计划中的位置信息进行比较和配准,从而得到校正后的空间位置信息,作为经过数据采集和预处理后的数据并传递给样本训练集生成模块作为输入数据;
红外空间位置探测单元通过记录红外线发射和接收器得到红外定位小球标记点的位置信息,采用红外线与红外定位小球标记点相连,接收检测区域的红外线,红外空间位置探测单元反馈红外定位小球标记点的位置和方向的变化信息,将转换的数据和信息输出给分析与处理模块;
在样本训练集生成模块中采集肿瘤患者正常呼吸的N(N>5)个周期,如选取N=10,基于体内的骨性标记点或植入标记点的位置信息,获得肿瘤位置的信息作为训练用的历史数据,采集肿瘤位置的历史数据信息,构建10个训练样本{(x1,y1),(x2,y2),L,(x10,y10)},构成一个训练集,作为预测函数的输入;其中输入xi=[xi1,xi2,...,xi10]∈R10,输入中的每个分量表征肿瘤位置属性;
数据学习及预测模块中的预测函数选择通过单隐层前馈神经网络来近似肿瘤位置输入与肿瘤位置输出之间的函数关系y=f(x),输出yi=[yi1,yi2,...,yi10]∈R10是经过预测得到的肿瘤位置结果,构建单隐层神经网络,网络结构如图2所示,m=10,设置隐层节点数为K和激活函数g(x),激活函数可选择为s型函数、基本周期函数、正弦函数、余弦函数中的一种,例如选择最常用的非线性变换函数Sigmoid函数曲线,函数定义如下:
进而在区间[-1,1]上,随机选择输入层到隐层的权值wj和隐层的偏置值bj
根据随机选择的输入层到隐层的权值wj和隐层的偏置值bj,通过将数据采集和预处理模块采集的肿瘤位置的历史数据样本信息输入计算隐层输出矩阵H:
通过计算得到的隐层输出矩阵H,进而计算输出层权值其中Y=[y1,...,y10]T为肿瘤位置样本输出组成的矩阵;
通过求得的输出层权值β,在数据学习及预测模块采用如下预测函数对肿瘤位置进行预测:
其中,f(x)代表肿瘤位置输入与肿瘤位置输出之间的函数关系,x代表肿瘤位置输入,β代表输出层权值,K为隐层节点数,g(x)为激活函数,wj是随机选择的输入层到隐层的权值,bj是隐层的偏置值。
将预测值和样本训练集生成模块肿瘤位置样本信息进行归一化均方根误差比较,如果满足精度要求则输出预测值,得到预测函数;如果精度不满足要求则重新选择样本集进行训练,直至得到满足误差精度的预测函数。
提供以上实施例仅仅是为了描述本发明的目的,而并非要限制本发明的范围。本发明的范围由所附权利要求限定。不脱离本发明的精神和原理而做出的各种等同替换和修改,均应涵盖在本发明的范围之内。

Claims (3)

1.一种肿瘤位置预测装置,其特征在于:包括数据集的采集和预处理模块,样本训练集生成模块和数据学习及预测模块;数据集的采集和预处理模块将采集和预处理后的数据通过计算机连接传送至样本训练集生成模块,再将样本训练集生成模块得到的数据传送至数据学习及预测模块;
数据集采集和预处理模块,包括红外空间位置探测单元和分析与处理模块;红外空间位置探测单元提供红外定位小球作为标记点的位置信息,反馈给分析与处理模块,用于实时显示和自动量化分析,将当前获取的位置信息和分析与处理模块中保存的原始治疗计划中的位置信息进行比较和配准,从而得到校正后的空间位置信息,作为经过数据采集和预处理后的数据并传递给样本训练集生成模块作为输入数据;
样本训练集生成模块,接收经过数据采集和预处理后的数据作为样本训练集生成的输入,基于这些输入数据从中选取肿瘤患者正常呼吸的N个周期,基于体内骨性标记点或植入标记点的位置信息,获得肿瘤位置信息作为训练用历史数据,构建N个样本训练集,作为肿瘤位置预测函数的输入{(x1,y1),(x2,y2),L,(xN,yN)},其中x1表示样本集1中肿瘤位置输入,y1表示样本集1中肿瘤位置输出,x2表示样本集2中肿瘤位置输入,y2表示样本集2中肿瘤位置输出,xN表示样本集N中肿瘤位置输入,yN表示样本集N中肿瘤位置输出,N个样本训练集肿瘤位置输入输出数据作为样本训练集生成模块得到的数据,传递给数据学习及预测模块;
数据学习及预测模块,接收由样本训练集生成模块得到的数据作为输入,用来模拟肿瘤位置输入与肿瘤位置输出之间的预测函数关系y=f(x);在数据学习过程中,y是经过预测得到的肿瘤位置输出,x是肿瘤位置输入,在此基础上构建单隐层神经网络,设置隐层节点数为K和激活函数g(x),激活函数可选择为s型函数、基本周期函数、正弦函数、余弦函数中的一种,进而在区间[-1,1]上,随机选择输入层到隐层的权值w和隐层的偏置值b,通过将采集的肿瘤位置的历史数据样本信息输入计算隐层输出矩阵H,计算输出层权值其中Y=[y1,...,yN]T为肿瘤位置样本输出组成的矩阵,进而通过求得的输出层权值β对肿瘤位置进行预测。
2.根据权利要求1所述的肿瘤位置预测装置,其特征在于:所述红外空间位置探测单元通过记录红外线发射和接收器得到红外定位小球标记点的位置信息,采用红外线与红外定位小球标记点相连,接收检测区域的红外线,红外空间位置探测单元反馈红外定位小球标记点的位置和方向的变化信息,将转换的数据和信息输出给分析与处理模块。
3.根据权利要求1所述的肿瘤位置预测装置,其特征在于:所述N>5个周期。
CN201610551244.8A 2016-07-13 2016-07-13 一种肿瘤位置预测装置 Expired - Fee Related CN105931262B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610551244.8A CN105931262B (zh) 2016-07-13 2016-07-13 一种肿瘤位置预测装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610551244.8A CN105931262B (zh) 2016-07-13 2016-07-13 一种肿瘤位置预测装置

Publications (2)

Publication Number Publication Date
CN105931262A CN105931262A (zh) 2016-09-07
CN105931262B true CN105931262B (zh) 2017-02-22

Family

ID=56827967

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610551244.8A Expired - Fee Related CN105931262B (zh) 2016-07-13 2016-07-13 一种肿瘤位置预测装置

Country Status (1)

Country Link
CN (1) CN105931262B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106563210B (zh) * 2016-11-10 2020-07-10 苏州大学 基于ut变换的放疗机器人肿瘤呼吸运动估计及预测方法
CN107126192B (zh) * 2017-04-18 2020-03-06 四川省肿瘤医院 一种肿瘤位置实时监测系统及其监测方法
CN106963383A (zh) * 2017-04-21 2017-07-21 南京大学 一种基于呼吸状态空间重构的体内组织呼吸运动估计方法
CN108037520B (zh) * 2017-12-27 2021-01-22 中国人民解放军战略支援部队信息工程大学 阵列幅相误差条件下基于神经网络的直接定位偏差修正方法
CN110378881B (zh) * 2019-07-05 2020-03-24 北京航空航天大学 一种基于深度学习的肿瘤定位系统
CN110681074B (zh) * 2019-10-29 2021-06-15 苏州大学 基于双向gru网络的肿瘤呼吸运动预测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101053531A (zh) * 2007-05-17 2007-10-17 上海交通大学 基于多模式增敏成像融合的早期肿瘤定位跟踪方法
CN102592283A (zh) * 2012-02-21 2012-07-18 华南理工大学 一种脑部肿瘤头皮定位图像处理方法
CN103598915A (zh) * 2013-11-13 2014-02-26 浙江德尚韵兴图像科技有限公司 一种肿瘤射频消融技术中精确定位三球覆盖肿瘤的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104338238A (zh) * 2013-07-24 2015-02-11 江苏富科思科技有限公司 一种快速精确定位肿瘤的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101053531A (zh) * 2007-05-17 2007-10-17 上海交通大学 基于多模式增敏成像融合的早期肿瘤定位跟踪方法
CN102592283A (zh) * 2012-02-21 2012-07-18 华南理工大学 一种脑部肿瘤头皮定位图像处理方法
CN103598915A (zh) * 2013-11-13 2014-02-26 浙江德尚韵兴图像科技有限公司 一种肿瘤射频消融技术中精确定位三球覆盖肿瘤的方法

Also Published As

Publication number Publication date
CN105931262A (zh) 2016-09-07

Similar Documents

Publication Publication Date Title
CN105931262B (zh) 一种肿瘤位置预测装置
Chan et al. Integration of AI and machine learning in radiotherapy QA
Nguyen et al. Incorporating human and learned domain knowledge into training deep neural networks: a differentiable dose‐volume histogram and adversarial inspired framework for generating Pareto optimal dose distributions in radiation therapy
Wang et al. A feasibility of respiration prediction based on deep Bi-LSTM for real-time tumor tracking
Brand et al. Predicting pregnancy status from mid-infrared spectroscopy in dairy cow milk using deep learning
CN110211664B (zh) 一种基于机器学习自动设计放射治疗方案的系统
CN108717866A (zh) 一种预测放疗计划剂量分布的方法、装置、设备及存储介质
CN112581475B (zh) 一种预测放疗计划的伽马通过率的方法
CN101623198A (zh) 动态肿瘤实时跟踪方法
CN110841205A (zh) 一种肿瘤患者精准剂量验证方法、装置及设备
Ren et al. An optimal neural network and concrete strength modeling
CN101628154A (zh) 基于预测的图像引导跟踪方法
Garijo et al. Computational evaluation of different numerical tools for the prediction of proximal femur loads from bone morphology
KR102639558B1 (ko) 관심영역별 골 성숙 분포를 이용한 성장 분석 예측 장치 및 방법
WO2023184598A1 (zh) 基于人工智能的心脏模拟器数据矫正系统和方法
CN111345823A (zh) 一种远程运动康复方法、装置和计算机可读存储介质
Mashayekhi et al. Site‐agnostic 3D dose distribution prediction with deep learning neural networks
CN109785962A (zh) 一种基于学习模型智能预测放疗中射野参数的方法和系统
Fan et al. Verification of the machine delivery parameters of a treatment plan via deep learning
CN115271033B (zh) 基于联邦知识蒸馏医学图像处理模型构建及其处理方法
Fan et al. Independent verification of brachytherapy treatment plan by using deep learning inference modeling
Reijtenbagh et al. Multi-center dosimetric predictions to improve plan quality for brachytherapy for cervical cancer treatment
Yip et al. Evaluating performance of a user‐trained MR lung tumor autocontouring algorithm in the context of intra‐and interobserver variations
Jaberi et al. Artificial neural network based gynaecological image-guided adaptive brachytherapy treatment planning correction of intra-fractional organs at risk dose variation
Zarr et al. Foundations of human spatial problem solving

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170222

CF01 Termination of patent right due to non-payment of annual fee