基于大滑环锥形束双源能谱CT的网络学习控制方法及设备与流程

基于大滑环锥形束双源能谱ct的网络学习控制方法及设备
技术领域
1.本发明属于加速器放射治疗技术领域，具体涉及一种基于大滑环锥形束双源能谱ct的网络学习控制方法。


背景技术：

2.从1971年hounsfield发明ct出现至今，ct升级一直侧重在图像质量的优化以及检查成功率的提升方面，转速、覆盖宽度、空间分辨率、采集通道等性能指标已发展到接近极限，导致ct成像问题愈发明显：ct扫描获得的是关于组织密度的差异图像，对于物质密度相近的组织无法区分，导致难于发现微小病灶和隐匿性病灶。尤其是进行肿瘤学检查时，ct更多的是基于病灶的密度变化来判断病灶的形态、大小、边缘，增强程度等，这更多是基于病变的解剖结构，而解剖结构变化往往意味着疾病发展到了中、后期，错过了诊疗的最佳时机。
3.多能谱ct的出现解决了上述问题，其突破传统ct的等密度病灶局限，形成直接反应病灶组织成分的图像，能够带来更丰富、精准、定量的多参数成像信息，对于肿瘤早期病灶发现、冠心病、大血管造影、肌骨成像等多种疾病影像学精准诊断具有独特价值。
4.随着现代计算机技术和医学影像技术的发展，肿瘤放射治疗已步入精确放疗时代，而精确放疗的基础是患者三维图像信息的精确获取，目前ct是三/四维图像信息获取的主要途径，它不仅能提供肿瘤靶区的信息，也能提供肿瘤周围正常组织的信息。同时放疗计划系统通过重建的ct图像提供的电子密度信息计算校准照射剂量，从而达到精确的肿瘤靶区剂量并减少正常组织的受照剂量。但由于受常规ct机架孔径限制而无法使所有的患者在理想的体位上进行模拟位，因此大孔径ct模拟定位系统在放疗领域逐渐崭露头角。大孔径ct模拟定位系统有别于常规诊断用ct，扫描孔径一般要求最好≥75cm。这是因为放疗模拟扫描，根据不同部位，会有一些特殊体位的摆位及需要结合一些放疗体架。很多小型医院，放疗科刚开始开展，病人不算多，可能科室有1-2台加速器，却没有适合放疗科使用的专业ct模拟机。
5.大孔径ct模拟定位系统包括：多套x光源和多套影像探测器。由于1s影像探测器需要均匀采多张图，因而对ct滑环旋转速度要求非常稳定精确，对多套影像探测器的触发时序也需要非常的同步和精确，对这大滑环ct的控制来说是一大挑战，也是本领域的攻克难题。


技术实现要素：

6.为了解决上述技术问题，本发明提出了一种基于大滑环锥形束双源能谱ct的网络学习控制方法及设备。
7.为了达到上述目的，本发明的技术方案如下：
8.一方面，本发明公开了一种基于大滑环锥形束双源能谱ct的网络学习控制方法，用于控制多套x光源以及多套影像探测器进行双源能谱ct扫描，包括以下内容：
9.采用基于ethercat网络学习的pid自适应控制方法控制大滑环ct旋转运动的速度，利用bp神经网络调节pid控制器的参数；
10.控制多套x光源同步出束；
11.采用一个控制信号同时控制多套影像探测器同步触发采图。
12.在上述技术方案的基础上，还可做如下改进：
13.作为优选的方案，采用基于ethercat网络学习的pid自适应控制方法控制大滑环ct旋转运动的速度，具体包括以下步骤：
14.s1.1：选定bp神经网络的结构；
15.s1.2：实时测量大滑环ct的位置；
16.s1.3：计算大滑环ct的实时速度，并计算速度误差和加速度；
17.s1.4：前m个时刻的速度误差和加速度作为运动控制器bp神经网络的输入，同时也作为智能驱动器的pid控制器输入；
18.s1.5：计算bp神经网络的各层神经元的输入和输出，输出层的输出即为pid控制器的三个可调参数k
p
、ki、kd；
19.s1.6：计算pid控制器的控制输出；
20.s1.7：计算修正bp神经网络输出层的加权系数；
21.s1.8：计算修正bp神经网络隐含层的加权系数；
22.s1.9：重复步骤s1.2～s1.9直至修正完所有pid控制器的采样点。
23.作为优选的方案，pid控制器根据定值r(t)与实际输出值y(t)构成控制偏差：
24.e(t)＝r(t)-y(t)；
25.连续情况下pid控制器的控制输出为：
[0026][0027]
式中，k
p
为比例系数，ti为积分时间常数，td为微分时间常数；
[0028]
离散情况下pid控制器的控制输出为：
[0029][0030]
δu(k)＝k
p
(e(k)-e(k-1)) kie(k) kd(e(k)-2e(k-1) e(k-2))；(2)
[0031]
式中，ki＝k
p
*t/ti，kd＝k
p
*t/td；
[0032]
t为采样周期，k为采样序号。
[0033]
作为优选的方案，bp神经网络的结构具体为：
[0034]
bp神经网络的输入层输入为：
[0035][0036]
式中，m为输入层的节点数；
[0037]
bp神经网络的隐含层输出为：
[0038][0039]
式中，w
ij(2)
为隐含层权系数；w
im(2)
为阈值，w
im(2)
＝θi，θi代表隐含层第i个神经元的阈值；
[0040]
f[
·
]为活化函数，f[
·
]＝tanh(x)；
[0041]
q为隐含层节点数；
[0042]
bp神经网络的输出层的输入输出为：
[0043][0044]
式中，w
li(3)
为输出层加权系数；w
lq(3)
为阈值，w
lq(3)
＝θ
l
，θ
l
代表输出层第l个神经元的阈值；
[0045]
g[
·
]为活化函数，满足g[
·
]＝(1/2)[1 tanh(x)]；
[0046]
上标(1)，(2)，(3)分别表示输入层、隐含层、输出层；
[0047]
设性能指标函数为：
[0048][0049]
式中，r代表期望速度，y代表实际速度，z代表速度误差；
[0050]
根据最速下降法修正网络的加权系数，即按j对加权系数的负梯度方向搜索调整，并附加一个搜索快速收敛全局极小的惯性项，得到下式：
[0051][0052]
式中：η为学习速率，α为惯性系数。
[0053]
作为优选的方案，选定bp神经网络的结构，具体包括以下内容：选定bp神经网络输入层的节点数m和隐含层节点数q，并给出bp神经网络隐含层加权系数的初值w
ij(2)
(0)、输出层加权系数的初值w
li(3)
(0)，选定学习速率η和惯性系数α,k＝1。
[0054]
作为优选的方案，bp神经网络输出层的加权系数的计算公式为：
[0055][0056]
式中，g'[
·
]＝g(x)[1-g(x)]；
[0057]
s1.7中，根据式(8)计算修正bp神经网络输出层的加权系数。
[0058]
作为优选的方案，bp神经网络隐含层的加权系数的计算公式为：
[0059][0060]
式中，f'[
·
]＝[1-f2(x)]/2；
[0061]
s1.8中，根据式(9)计算修正bp神经网络隐含层的加权系数。
[0062]
作为优选的方案，多套x光源的开始时间同步，且结束时间同步。
[0063]
作为优选的方案，多套影像探测器同步触发采图具体包括以下内容：
[0064]
大滑环ct开始旋转，发送影像触发脉冲触发多套影像探测器同步采图；
[0065]
根据大滑环ct实时速度v、当前大滑环ct的位置编码器值y和前几个时刻大滑环ct的位置编码器值yi，利用三次样条预测算法预测下一时刻的位置编码器值yj，同时预测下一时刻的旋转速度vj，则下一个影像触发脉冲的时间为360/(t*p)*x/vj；
[0066]
其中：t为大滑环ct旋转一周的时间，p为1s影像探测器均匀采图的张数，x为大滑环ct旋转1
°
的编码器值。
[0067]
另一方面，本发明还公开一种基于大滑环锥形束双源能谱ct的网络学习控制设备，包括：
[0068]
一个或多个处理器；
[0069]
存储器；
[0070]
以及一个或多个程序，其中一个或多个程序存储在存储器中并被配置为由一个或多个处理器执行，一个或多个程序包括上述任一种基于大滑环锥形束双源能谱ct的网络学习控制方法的指令。
[0071]
本发明一种基于大滑环锥形束双源能谱ct的网络学习控制方法及设备，采用基于ethercat网络学习的pid自适应控制方法控制大滑环ct旋转运动的速度，利用bp神经网络调节pid控制器的参数。pid控制器参数的自适应调节能够兼顾速度控制的快速性、稳定性和鲁棒性，采用一个控制信号同时控制多套影像探测器的触发。
附图说明
[0072]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0073]
图1为本发明实施例提供的大滑环锥形束双源能谱ct的结构示意图。
[0074]
图2为本发明实施例提供的基于ethercat网络学习的pid自适应控制方法流程图。
[0075]
图3为本发明实施例提供的基于网络学习的控制系统结构图。
[0076]
图4为本发明实施例提供的bp神经网络(nn)结构图。
[0077]
图5为本发明实施例提供的大滑环ct旋转过程中影像板的触发控制结构图。
具体实施方式
[0078]
下面结合附图详细说明本发明的优选实施方式。
[0079]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0080]“包括”元件的表述是“开放式”表述，该“开放式”表述仅仅是指存在对应的部件或步骤，不应当解释为排除附加的部件或步骤。
[0081]
为了达到本发明的目的，一种基于大滑环锥形束双源能谱ct的网络学习控制方法及设备的其中一些实施例中，基于大滑环锥形束双源能谱ct的网络学习控制方法用于控制如图1所示两套x光源以及四套影像探测器进行双源能谱ct扫描，大滑环ct旋转按照一周4s或6s的速度进行匀速旋转。
[0082]
网络学习控制方法包括以下内容：
[0083]
采用基于ethercat网络学习的pid自适应控制方法控制大滑环ct旋转运动的速度，利用bp神经网络调节pid控制器的参数；
[0084]
控制两套x光源同步出束；
[0085]
采用一个控制信号同时控制四套影像探测器同步触发采图。
[0086]
x光源可以但不限于为x球管，影像探测器可以但不限于为探测板。
[0087]
为了进一步地优化本发明的实施效果，在另外一些实施方式中，其余特征技术相同，不同之处在于，如图2所示，采用基于ethercat网络学习的pid自适应控制方法控制大滑环ct旋转运动的速度，具体包括以下步骤：
[0088]
s1.1：选定bp神经网络的结构；
[0089]
s1.2：实时测量大滑环ct的位置；
[0090]
s1.3：计算大滑环ct的实时速度，并计算速度误差和加速度；
[0091]
s1.4：前m个时刻的速度误差和加速度作为运动控制器bp神经网络的输入，同时也作为智能驱动器的pid控制器输入；
[0092]
s1.5：计算bp神经网络的各层神经元的输入和输出，输出层的输出即为pid控制器的三个可调参数k
p
、ki、kd；
[0093]
s1.6：计算pid控制器的控制输出；
[0094]
s1.7：计算修正bp神经网络输出层的加权系数；
[0095]
s1.8：计算修正bp神经网络隐含层的加权系数；
[0096]
s1.9：重复步骤s1.2～s1.9直至修正完所有pid控制器的采样点。
[0097]
大滑环ct旋转运动的速度控制算法采用基于ethercat网络学习的pid自适应控制方法，pid控制参数的优化在dsp fpga的运动控制器完成，pid控制器在具有ethercat网络的智能运动驱动器完成，如图3所示。pid控制具有结构简单、稳定性能好、可靠性高等优点，尤其适用于可建立精确数学模型的确定性控制系统。但是，在实际的应用中，大多数工业过程都不同程度地存在非线性、参数时变性和模型不确定性，因而，一般的pid控制无法实现对这样的过程的精确控制。
[0098]
基于网络学习的pid自适应控制方法对数学模型的依赖性弱，不需要建立过程的
精确数学模型。网络学习算法采用bp神经网络学习算法，可以说是目前最成功的神经网络学习算法。bp神经网络就是一个“万能的模型 误差修正函数”，每次根据训练得到的结果与预想结果进行误差分析，进而修改权值和阈值，一步一步得到能输出和预想结果一致的模型。
[0099]
进一步，在上述实施例的基础上，pid控制器是一种线性控制器，根据定值r(t)与实际输出值y(t)构成控制偏差：
[0100]
e(t)＝r(t)-y(t)；
[0101]
连续情况下pid控制器的控制输出为：
[0102][0103]
式中，k
p
为比例系数，ti为积分时间常数，td为微分时间常数；
[0104]
离散情况下pid控制器的控制输出为：
[0105][0106]
δu(k)＝k
p
(e(k)-e(k-1)) kie(k) kd(e(k)-2e(k-1) e(k-2))；(2)
[0107]
式中，ki＝k
p
*t/ti，kd＝k
p
*t/td；
[0108]
t为采样周期，k为采样序号。
[0109]
网络学习环节采用bp神经网络控制器，控制回路采用普通的pid控制器，其系统结构如图3所示。
[0110]
控制系统由二部分组成：
[0111]
(1)传统pid控制器：直接对被控对象过程进行闭环控制，并且三个参数k
p
、ki、kd为在线整定式；
[0112]
(2)神经网络(nn)：经过神经网络的学习，以加权系数的形式表现出来，规则的生成就转化为加权系数初值的确定和修改。根据系统的运行状态，调节pid控制器的参数，以达到性能指标的最优化。即输出层神经元的输出状态对应于pid控制器的三个可调参数k
p
、ki、kd，通过神经网络的自身学习、加权系数调整，从而使其稳定状态对应于误差平方和最小的最优控制律下的pid控制器参数。
[0113]
本发明采用三层bp神经网络，其网络结构如图4所示，有m个输入节点、q个隐含节点、3个输出节点。输入节点对应经预测处理后的系统状态变量e{i}。输出节点对应pid控制器的三个可调参数kp、ki、kd。由于k
p
、ki、kd不能为负值，所以输出层神经元的活化函数取非负的sigmoid，而隐含层神经元的活化函数可取正负对称的sigmoid函数。
[0114]
进一步，在上述实施例的基础上，bp神经网络的结构具体为：
[0115]
bp神经网络的输入层输入为：
[0116][0117]
式中，m为输入层的节点数，输入变量的个数m取决于被控系统的复杂程度；
[0118]
bp神经网络的隐含层输出为：
[0119][0120]
式中，w
ij(2)
为隐含层权系数；w
im(2)
为阈值，w
im(2)
＝θi，θi代表隐含层第i个神经元的阈值，也就是最大或最小值；
[0121]
f[
·
]为活化函数，f[
·
]＝tanh(x)；
[0122]
q为隐含层节点数；
[0123]
bp神经网络的输出层的输入输出为：
[0124][0125]
式中，w
li(3)
为输出层加权系数；w
lq(3)
为阈值，w
1q(3)
＝θ
l
，θ
l
代表输出层第l个神经元的阈值，也就是最大或最小值；
[0126]
g[
·
]为活化函数，满足g[
·
]＝(1/2)[1 tanh(x)]；
[0127]
上标(1)，(2)，(3)分别表示输入层、隐含层、输出层；
[0128]
设性能指标函数为：
[0129][0130]
式中，r代表期望速度，y代表实际速度，z代表速度误差；
[0131]
根据最速下降法修正网络的加权系数，即按j对加权系数的负梯度方向搜索调整，并附加一个搜索快速收敛全局极小的惯性项，得到下式：
[0132][0133]
式中：η为学习速率，α为惯性系数。
[0134]
进一步，在上述实施例的基础上，选定bp神经网络的结构，具体包括以下内容：选定bp神经网络输入层的节点数m和隐含层节点数q，并给出bp神经网络隐含层加权系数的初值w
ij(2)
(0)、输出层加权系数的初值w
li(3)
(0)，选定学习速率η和惯性系数α,k＝1。
[0135]
进一步，在上述实施例的基础上，bp神经网络输出层的加权系数的计算公式为：
[0136][0137]
式中，g'[
·
]＝g(x)[1-g(x)]；
[0138]
s1.7中，根据式(8)计算修正bp神经网络输出层的加权系数。
[0139]
进一步，在上述实施例的基础上，bp神经网络隐含层的加权系数的计算公式为：
[0140][0141]
式中，f
′
[
·
]＝[1-f2(x)]/2；
[0142]
s1.8中，根据式(9)计算修正bp神经网络隐含层的加权系数。
[0143]
具体地，网络学习控制方法的计算步骤可归纳如下：
[0144]
1)选定bp神经网络的结构，即选定bp神经网络输入层的节点数m和隐含层节点数q，并给出bp神经网络隐含层加权系数的初值w
ij(2)
(0)、输出层加权系数的初值w
li(3)
(0)，选定学习速率η和惯性系数α，k＝1；
[0145]
2)根据已知的e1，e2，
…
，em、以及前一段时间的速度误差和其微分，共同作为bp神经网络的输入，e1，e2，
…
，em为实际时刻的速度误差，为其对应微分；
[0146]
3)根据式(3)-(5)前向计算bp神经网络的各层神经元的输入和输出，输出层的输出即为pid控制器的三个可调参数k
p
(n)、ki(n)、kd(n)；
[0147]
4)根据式u(n)＝u(n-1) k
p
[e(n)-e(n-1)] kie(n) kd[e(n)-2e(n-1) e(n-2)]计算pid控制器的控制输出u(n)，参与控制和计算；
[0148]
5)由式(8)，计算修正输出层的加权系数w
li(3)
(n)；
[0149]
6)由式(9)，计算修正隐含层的加权系数w
ij(2)
(n)；
[0150]
7)置k＝k 1，ei＝e
i 1
返回步骤2)。
[0151]
上述步骤中：k＝1，2，...n。
[0152]
本发明采用了基于ethercat网络学习的pid自适应控制方法，pid控制参数的优化在dsp fpga的运动控制器完成，pid控制器在具有ethercat网络的智能运动驱动器完成。在驱动器完成基本控制，保证了控制的稳定性和可靠性，pid控制器参数优化在运动控制器完成，保证了控制性能的优化、精准和抗干扰的鲁棒性。
[0153]
为了进一步地优化本发明的实施效果，在另外一些实施方式中，其余特征技术相同，不同之处在于，两套x光源的开始时间同步，且结束时间同步。
[0154]
二套x球管a，b进行同步连续出束，不是脉冲出束，所以只要保证开始时间和结束时间同步就可以。
[0155]
为了进一步地优化本发明的实施效果，在另外一些实施方式中，其余特征技术相同，不同之处在于，四个影像板a1，a2，b1，b2触发需要同步、等间隔快速触发，控制结构如图5所示。
[0156]
四个影像板同步触发采图具体包括以下内容：
[0157]
大滑环ct开始旋转，发送影像触发脉冲触发多套影像探测器同步采图；
[0158]
根据大滑环ct实时速度v、当前大滑环ct的位置编码器值y和前几个时刻大滑环ct的位置编码器值yi，利用三次样条预测算法预测下一时刻的位置编码器值yj，同时预测下一时刻的旋转速度vj，则下一个影像触发脉冲的时间为360/(t*p)*x/vj；
[0159]
其中：t为大滑环ct旋转一周的时间，p为1s影像探测器均匀采图的张数，x为大滑
环ct旋转1
°
的编码器值。
[0160]
具体地，以大滑环ct旋转一周为4s为例，p＝80,则t*p＝320，可以采集320张图，需要320个影像触发脉冲，影像板触发控制流程如下：
[0161]
(1)大滑环ct开始旋转的角度是-10
°
，停止旋转的角度是370
°
，保证在0
°
至360
°
之间匀速运动；
[0162]
(2)大滑环ct开始旋转时，x球管a、b开始连续出束，ct停止旋转时，x球管a、b停止出束，保证在0
°
至360
°
之间出束正常；
[0163]
(3)根据现在实时的速度v，现在的位置编码器值y和前几个时刻的位置编码器值yi，利用三次样条预测算法预测下一时刻的位置编码器值yj，从而可以预测出下一时刻的旋转速度vj，下一个影像触发脉冲的时间就是360/320*x/vj；
[0164]
(4)dsp控制器有高精度的定时器，根据预估时间360/320*x/vj进行精确定时触发。
[0165]
由于有了旋转速度预测，可以保证在旋转过程中，速度有微小波动情况下的影像板等角度精确触发。若没有预测，则理论上每间隔360/320*x进行脉冲触发。
[0166]
此外，本发明实施例还公开一种基于大滑环锥形束双源能谱ct的网络学习控制设备，包括：
[0167]
一个或多个处理器；
[0168]
存储器；
[0169]
以及一个或多个程序，其中一个或多个程序存储在存储器中并被配置为由一个或多个处理器执行，一个或多个程序包括上述任一实施例公开的基于大滑环锥形束双源能谱ct的网络学习控制方法的指令。
[0170]
本发明一种基于大滑环锥形束双源能谱ct的网络学习控制方法及设备具有以下有益效果：
[0171]
第一，利用两套x光源以及四套影像探测器进行双源能谱ct扫描，每套光源增加了一套影像探测器，增加采图的fov，能够诊断更多的病例，能够适应更多的人群。
[0172]
第二，采用基于ethercat网络学习的pid自适应控制方法控制大滑环ct旋转运动的速度，利用bp神经网络调节pid控制器的参数。pid控制器参数的自适应调节能够兼顾速度控制的快速性、稳定性和鲁棒性。本发明大滑环双源能谱ct的旋转速度稳定、可靠且精准。
[0173]
第三，采用一个控制信号同时控制多套影像探测器的触发，保证了影像触发的同步性。触发信号的时序控制，利用dsp控制器有高精度的定时器和对下一时刻速度的预测，预测算法采用三次样条预测算法进行误差平方和最小的预测，保证在旋转过程中，速度有微小波动情况下的影像板等角度精确触发。
[0174]
应当理解，这里描述的各种技术可结合硬件或软件，或者它们的组合一起实现。从而，本发明的方法和设备，或者本发明的方法和设备的某些方面或部分可采取嵌入有形媒介，例如软盘、cd-rom、硬盘驱动器或者其它任意机器可读的存储介质中的程序代码(即指令)的形式，其中当程序被载入诸如计算机之类的机器，并被该机器执行时，该机器变成实践本发明的设备。
[0175]
上述实施例只为说明本发明的技术构思及特点，其目的在于让本领域普通技术人
员能够了解本发明的内容并加以实施，并不能以此限制本发明的保护范围，凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围内。