发明名称 呼吸同步跟踪系统的位置信号的模糊预测方法及预测器
摘要 呼吸同步跟踪系统在治疗射线投放时,可以监测病人的呼吸,并对肿瘤目标运动作出补偿。在治疗射线投放系统的指令和响应之间存在时间延迟,这种延迟可能会导致意外的、甚至危险的振动。为了提高同步跟踪性能,可以用信号预测技术来弥补时间延迟。本发明提供了一种新的同步跟踪系统的位置信号的模糊预测方法以及预测器,该预测器能够处理随时间变化的幅度和频率问题,比较好的解决了放射治疗中同步跟踪的问题。
申请公布号 CN101794357A 申请公布日期 2010.08.04
申请号 CN201010134869.7 申请日期 2010.03.30
申请人 江苏瑞尔医疗科技有限公司 发明人 盛晔
分类号 G06F19/00(2006.01)I;A61N5/00(2006.01)I 主分类号 G06F19/00(2006.01)I
代理机构 苏州华博知识产权代理有限公司 32232 代理人 傅靖
主权项 1.一种呼吸同步跟踪系统的位置信号的模糊预测方法,其特征在于包括以下步骤:(1)在病人体表确定一体表位置P<sub>0</sub>并计算出相对应的体内病灶的目标位置参数P(x,y,z);(2)采集该位置P(x,y,z)在一段时间内的位置参数(x,y,z)存入计算机;(3)根据P(x,y,z)的位置参数计算出相应的速度参数(dx,dy,dz)并存入计算机:<maths num="0001"><![CDATA[<math><mrow><msub><mi>dx</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>x</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>x</mi><mi>i</mi></msub></mrow><mrow><msub><mi>t</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>t</mi><mi>i</mi></msub></mrow></mfrac></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>dy</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>y</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>y</mi><mi>i</mi></msub></mrow><mrow><msub><mi>t</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>t</mi><mi>i</mi></msub></mrow></mfrac></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><msub><mi>dz</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>z</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>z</mi><mi>i</mi></msub></mrow><mrow><msub><mi>t</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>-</mo><msub><mi>t</mi><mi>i</mi></msub></mrow></mfrac><mo>;</mo></mrow></math>]]></maths>(4)根据计算机内存储的位置参数(x,y,z),找出最大和最小的位置(x<sub>max</sub>,y<sub>max</sub>,z<sub>max</sub>)和(x<sub>min</sub>,y<sub>min</sub>,z<sub>min</sub>);(5)根据计算机内存储的速度参数(dx,dy,dz),找出最大和最小的速度(dx<sub>max</sub>,dy<sub>max</sub>,dz<sub>max</sub>)和(dx<sub>min</sub>,dy<sub>min</sub>,dz<sub>min</sub>);(6)制定模糊规则:根据最大和最小的位置参数,把x模糊划分出n1个子集,这n1个子集的每一个子集再被根据最大和最小的速度参数模糊划分为n2个子集,得到M=n1*n2个子模型,每个子模型都对应一个规则,该规则用线性方程<maths num="0004"><![CDATA[<math><mrow><msub><mover><mi>x</mi><mo>^</mo></mover><mi>m</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>a</mi><mrow><mi>m</mi><mo>,</mo><mn>0</mn></mrow></msub><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>+</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>+</mo><msub><mi>a</mi><mrow><mi>m</mi><mo>,</mo><mi>r</mi></mrow></msub><mi>x</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>r</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>c</mi><mi>m</mi></msub></mrow></math>]]></maths>表示,<img file="FSA00000065570600015.GIF" wi="48" he="52" />是第m规则的模式输出,m的范围是在1~M,而M是规则的总数,p是预测步骤,r是预测阶数,a<sub>m</sub>和c<sub>m</sub>是在第m规则的参数;对于y,z作同样处理;(7)根据各个子模型所对应的计算机内已经存储的已知的位置参数(x,y,z)、速度参数(dx,dy,dz),采用梯度下降算法,训练出对应子模型的参数a<sub>m</sub>和c<sub>m</sub>;(8)采集下一时间点体表目标的位置参数,并计算出该位置参数对应的体内病灶的位置参数,根据已经制定的模糊规则,把计算出的参数映射到相应的子模型里,得到输出<img file="FSA00000065570600021.GIF" wi="664" he="53" />对相应映射的子模型的输出结果加权平均,即为所需要的预测结果,<maths num="0005"><![CDATA[<math><mrow><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mover><mi>x</mi><mo>^</mo></mover><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub></mrow></mfrac></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><mover><mi>y</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mover><mi>y</mi><mo>^</mo></mover><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub></mrow></mfrac></mrow></math>]]></maths><maths num="0007"><![CDATA[<math><mrow><mover><mi>z</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub><mo>&CenterDot;</mo><msub><mover><mi>z</mi><mo>^</mo></mover><mi>i</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mi>p</mi><mo>)</mo></mrow></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>&mu;</mi><mi>i</mi></msub></mrow></mfrac></mrow></math>]]></maths>μ<sub>i</sub>是子模型的隶属度。
地址 214000 江苏省无锡市锡山经济开发区芙蓉中三路99号科创园