一种放射影像引导下诊治肺肿瘤装置的制作方法

技术特征：

1.一种放射影像引导下诊治肺肿瘤装置，其特征在于，所述放射影像引导下诊治肺肿瘤装置由CT重建设备和CT检查床组成；

所述CT重建设备的一侧放置有CT检查床；

所述CT重建设备由探测器、光管、360°旋转载物平台运动系统、探测器Y向直线运动系统、探测器X向直线运动系统、探测器Z向直线运动系统、光管Z向直线运动系统和被测物组成；

所述探测器Y向直线运动系统、探测器X向直线运动系统、探测器Z向直线运动系统设置在扁平构件的CT重建设备的上方；探测器Y向直线运动系统上架设探测器X向直线运动系统、探测器X向直线运动系统上架架设探测器Z向直线运动系统，探测器固定在探测器Z向直线运动系统上；所述360°旋转载物平台运动系统位于整体结构的中间，所述光管Z向直线运动系统位于扁平构件的CT重建设备下方，光管固定在光管Z向直线运动系统；

所述CT检查床由床体、头枕、头枕导轨、头枕伸缩杆、手臂固定器、手臂固定器伸缩杆、电机、脚轮、脚轮锁止器、床腿、床腿高度调节器、头枕伸缩杆控制键、手臂固定器伸缩杆控制键、床腿高度调节器控制键、电源开关和电源接口，所述的头枕设置在头枕导轨上面；所述的头枕导轨固定在床体上面；所述的头枕伸缩杆一端固定于床体上面，另一端与头枕相连；所述的手臂固定器伸缩杆固定于床体上面；所述的手臂固定器设置在手臂固定器伸缩杆上；所述的电机设置在床体下面中间部位；所述的床腿设置在床体下面四个角内里面；所述的脚轮设置在床腿下面；所述的脚轮锁止器设置于脚轮上面；所述的床腿高度调节器设置于床腿中间；所述的头枕伸缩杆控制键设置在床体侧面，通过导线分别与头枕伸缩杆和电机相连；所述的手臂固定器伸缩杆控制键设置在床体侧面，通过导线分别与手臂固定器伸缩杆和电机相连；所述的床腿高度调节器控制键设置在床体侧面，通过导线分别与床腿高度调节器和电机相连；所述的电源开关在床体侧面；所述的电源接口设置在床体侧面。

2.如权利要求1所述的放射影像引导下诊治肺肿瘤装置，其特征在于，所述探测器与计算机连接，所述计算机设置有图像处理重建模块，所述图像处理重建模块的图像处理方法包括：

第一步，采用锥束滤波反投影算法得到初始图像u₀：对XCT待重建的M幅投影采用滤波反投影算法重建，对圆弧轨迹扫描的锥束投影数据进行半精确重建；圆弧轨迹是短扫描轨迹，即π加扇角扫描轨迹；也是超短扫描轨迹，即小于短扫描的轨迹包括：权重投影数据的偏导数；沿水平与非水平方向的一维希尔伯特变换；圆弧轨迹的反投影；将重建后图像u₀作为初始图像；

第二步，令k＝0，k为迭代次数(k＝0，1，2...K-1)，将u₀作为迭代初始图像u⁽⁰⁾，采用代数迭代法更新图像：

$u_j^{(k+1)}=u_j^{\left(k\right)}+\mathrm{\λ}\frac{(p_i-{\mathrm{\Σ}}_{n=1}^N{u}_n^{\left(k\right)}\·m_{in})}{{\mathrm{\Σ}}_{n=1}^N{m}_{in}^2}{m}_{ij};$

其中，k为迭代次数(k＝0，1，2...K)，即为第k÷1次迭代后更新的图像，为更新前图像，j代表图像角标，p_i为第i个探测器单元的投影值，m_m代表图像中第n个像素值对探测器单元i的贡献大小，λ为松弛因子，取值范围0～2之间；

第三步，对得到的图像进行非负值约束：对上步中所有小于0的值赋值为零，即如果则令

第四步，根据梯度下降法对u^(k+1)求图像总变差最小化后的图像：

$u^{(k+1,l+1)}=u^{(k+1,l)}-a\left|\right|u^{(k+1,0)}-u^{\left(k\right)}\left|\right|*d{\hat{u}}^{(k+1,l)};$

其中l为梯度下降法的迭代次数(l＝0，1...L-1)，a为步长因子；

第五步，令k＝k+1，u^(k)＝u^(k，L+1)，重复步第二步-第四步，直到满足迭代次数要求。

3.如权利要求1所述的放射影像引导下诊治肺肿瘤装置，其特征在于，所述CT重建设备进一步设置有CT扫描成像检测优化模块，所述CT扫描成像检测优化的方法包括：

前一帧图像通过单帧图像方法进行行人检测，提取检测出的矩形框的坐标并保存；扫描窗口的下边缘坐标小于等于上区域的下边界时用最小窗口遍历匹配；扫描窗口下边缘坐标大于上区域的下边界小于等于中上区域的下边界时用较小窗口遍历匹配；扫描窗口下边缘坐标大于中上区域的下边界小于等于中下区域的下边界时用较大窗口遍历匹配；扫描窗口下边缘坐标大于中下区域下边界时用较大窗口遍历匹配；单帧图像采用划分区域遍历式方法，把图像从上到下依次分为上、中上、中下、下四个区域，并用不同比例的扫描窗口依次对每一区域进行匹配；根据前一帧图像检测的结果对下一帧图像进行跟踪和识别，通过改变前一帧矩形框的横纵坐标值使矩形框放大一定的比例作为当前帧的感兴趣区域，当前帧图像只需对感兴趣区域通过的单帧图像的方法进行扫描窗口遍历匹配，并且根据扫描窗口的下边界区域判断窗口是否会增大。

4.如权利要求1所述的放射影像引导下诊治肺肿瘤装置，其特征在于，所述电机的电子限位控制方法步骤为：

步骤一，调整电机与传感器的相序，若电机的相序与传感器的相序相匹配，电机输出轴受电机控制器控制约转半圈，正转和反转能自动停止；

若电机的相序与传感器的相序不匹配，则电机不会自停，电机只正转不能反转；

步骤二，在显示控制模块上选定电机运行的位置A和位置B，在显示控制模块上按下设置键，设置灯闪烁，将电机输出轴运转至A位置，输出轴向B位置的方向运转，通过高位键和低位键调控电机的上、下位置，按下设置键，设置灯停止闪烁，电机A位置和B位置均设置完成；

步骤三，将电机的输出轴运转至B位置，设置灯闪烁，输出轴向A位置方向运转，通过低位键和高位键调控电机的上、下位置，设置灯停止闪烁；按下设置键，电机B位置和A位置均设置完成；

步骤四，通过电机控制器内的记忆模块将此调整好的相序记忆；

步骤五，设置好电机输出轴的输出路径后，当电机控制器发出正转指令后由电机输出轴带动减速机正转运行，同时传感器发出正转脉冲，通过电机控制器内的信号接收模块接收由传感器发送的信号后，将此信号发送至显示控制模块，显示控制模块检测若此信号是与记忆模块内的相序相符，则向电机发出停机指令由电机输出轴带动减速机停机，减速机停止运行，反之，若显示控制模块检测此信号与记忆模块内的相序不相符，则不发出停机指令，减速机正常运行；

步骤六，设置好电机输出轴的输出路径后，当电机控制器发出反转指令后由电机输出轴带动减速机反转运行，同时传感器发出反转脉冲，通过电机控制器内的信号接收模块接收由传感器发送的信号后，将此信号发送至显示控制模块，显示控制模块检测若此信号是与记忆模块内的相序相符，则向电机发出停机指令由电机输出轴带动减速机停机，减速机停止运行，反之，若显示控制模块检测此信号与记忆模块内的相序不相符，则不发出停机指令，减速机正常运行。

5.如权利要求4所述的放射影像引导下诊治肺肿瘤装置，其特征在于，所述电机控制器设置有混合矩阵列向量和跳频频率估计模块，所述混合矩阵列向量和跳频频率估计模块的信号处理方法包括：

第一步，在p(p＝0，1，2，...P-1)时刻，对表示的频率值进行聚类，得到的聚类中心个数表示p时刻存在的载频个数，个聚类中心则表示载频的大小，分别用表示；

第二步，对每一采样时刻p(p＝0，1，2，...P-1)，利用聚类算法对进行聚类，同样可得到个聚类中心，用表示；

第三步，对所有求均值并取整，得到源信号个数的估计即

$\hat{N}=round\left(\frac{1}{p}\underset{p=0}{\overset{P-1}{\Σ}}{\hat{N}}_p\right);$

第四步，找出的时刻，用p_h表示，对每一段连续取值的p_h求中值，用表示第l段相连p_h的中值，则表示第l个频率跳变时刻的估计；

第五步，根据第二步中估计得到的以及第四步中估计得到的频率跳变时刻估计出每一跳对应的个混合矩阵列向量具体公式为：

${\hat{a}}_n\left(l\right)=\left\{\begin{array}{cc}\frac{1}{{\stackrel{\‾}{p}}_h\left(1\right)}\·\underset{p=1,p\≠p_h}{\overset{{\stackrel{\‾}{p}}_h\left(1\right)}{\Σ}}{b}_{n,p}^0& l=1,\\ \frac{1}{{\stackrel{\‾}{p}}_h\left(l\right)-{\stackrel{\‾}{p}}_h(l-1)}\·\underset{p={\stackrel{\‾}{p}}_h(l-1)+1,p\≠p_h}{\overset{{\stackrel{\‾}{p}}_h\left(l\right)}{\Σ}}{b}_{n,p}^0& l\>1,\end{array}\right.,n=1,2,\mathrm{...}\hat{N}$

这里表示第l跳对应的个混合矩阵列向量估计值；

第六步，估计每一跳对应的载频频率，用表示第l跳对应的个频率估计值，计算公式如下：

${\hat{f}}_{o,n}\left(l\right)=\left\{\begin{array}{cc}\frac{1}{{\stackrel{\‾}{p}}_h\left(1\right)}\·\underset{p=1,p\≠p_h}{\overset{{\stackrel{\‾}{p}}_h\left(1\right)}{\Σ}}{f}_o^n\left(p\right)& l=1,\\ \frac{1}{{\stackrel{\‾}{p}}_h\left(l\right)-{\stackrel{\‾}{p}}_h(l-1)}\·\underset{p={\stackrel{\‾}{p}}_h(l-1)+1,p\≠p_h}{\overset{{\stackrel{\‾}{p}}_h\left(l\right)}{\Σ}}{f}_o^n\left(p\right)& l\>1,\end{array}\right.,n=1,2,\mathrm{...}\hat{N}.$