(19)中华人民共和国国家知识产权局
(12)发明专利说明书
(10)申请公布号 CN 105796121 A
(43)申请公布日 2016.07.27
(21)申请号 CN201610118084.8
(22)申请日 2016.03.02
(71)申请人 中国人民解放军第四军医大学
    地址 710032 陕西省西安市长乐西路169号
(72)发明人 戎军艳 卢虹冰 高鹏 廖琪梅 刘文磊 刘洋 见伟平
(74)专利代理机构 天津市鼎和专利商标代理有限公司
    代理人 李凤
(51)Int.CI
      正则化参数的自适应估计
                                                                  权利要求说明书 说明书 幅图
(54)发明名称
      一种CT和X射线激发荧光双模同步断层成像方法
(57)摘要
      本发明公开了一种CT和X射线激发荧光双模同步断层成像方法，可同步实现小动物X射线激发荧光断层成像和CT成像，提供更准确的成像物体解剖结构和分子影像信息。通过减少CT投影角度数并设置CT投影角度数与XLCT相同或成倍数，保证XLCT投影与CT投影同地点同时刻采集；同时利用压缩感知原理，对CT投影采用图像总变差最小化的少角度重建算法，对XLCT迭代重建引入预条件矩阵和小波变换以降低重建的病态性。本方法的实施可实现CT/XLCT的同步成像，简化双模成像流程，降低实验难度，减少成像时间，减少实验失败率，降低辐射剂量，避免成像过程中物体移动引起的误差，降低图像配准和融合难度。
法律状态
法律状态公告日
法律状态信息
法律状态
权 利 要 求 说 明 书
1.一种CT和X射线激发荧光双模同步断层成像方法，其特征在于，该CT和X射线激发荧            光双模同步断层成像方法通过减少CT投影角度数并设置CT投影角度数与XLCT相同或成倍            数，同时调整单幅投影采集时间，使双模数据同步采集；利用压缩感知原理，对CT投影采用            图像总变差最小化的少角度重建算法，对XLCT重建引入提高重建质量的预条件矩阵和小波            变换以降低光学重建的病态性；通过基于标记点的刚性配准方法对重建得到的CT和XLCT断            层图像进行配准并融合，实现CT和XLCT的同步断层成像。           
2.如权利要求1所述的CT和X射线激发荧光双模同步断层成像方法，其特征在于，所述            CT和X射线激发荧光双模同步断层成像方法包括以下步骤：           
步骤一、设置CT/XLCT同步成像参数：设置CT在360度范围内采集的投影数为M，10≤M≤            30，XLCT荧光采集投影数为J，M＝qJ，其中M、J和q均为正整数，并设置单幅投影采集时间r            秒，CT投影采集所用平板探测器积分时间为s秒，XLCT投影采集所用CCD相机积分时间为h            秒，s<＝r，h<＝r，同时根据成像物体的需求设置X射线源的电压、电流；           
步骤二、CT暗场数据和亮场数据采集：在关闭X射线源状态下，利用平板探测器按照s秒            积分时间采集暗场数据；打开X射线源，利用平板探测器按照s秒积分时间采集亮场数据；           
步骤三、成像物体准备：将成像物体注入荧光探针后置于转台上，固定好；           
步骤四、白光数据采集：在关闭X射线源状态下，打开外界光源，在360度范围内均匀采            集Q幅成像物体的白光图像并保存；           
步骤五、CT/XLCT投影数据同步采集：关闭外界光源，打开X射线源，利用转台控制物体            的旋转，在360度范围内均匀采集M幅CT投影数据，同时均匀采集J幅XLCT投影，采集完成后            两套数据均保存至计算机，并关闭X射线源的成像设备；           
步骤六、CT投影数据预处理并重建：利用亮场、暗场数据及平板探测器自带的坏像素            表，对CT投影数据进行亮场校正，暗场校正以及坏像素校正，生成待重建投影，对待重建投            影利用图像总变差最小化方法进行少角度重建，得到CT断层图像；           
步骤七、XLCT投影数据预处理并重建：对步骤五采集到的XLCT投影数据进行预处理，去            除图像中的孤立亮点和坏点，利用步骤四中采集得到的白光数据构建物体轮廓并剖分，将            各个角度下的经过预处理之后的XLCT少角度投影进行小波变换，得到变换之后的投影并映            射到构建好的轮廓表面，得到表面光强分布矩阵，对系统矩阵和表面光强分布矩阵采用预            条件矩阵处理，建立系统矩阵与表面光强分布矩阵之间的联系，采用迭代重建算法得到            XLCT荧光断层图像；           
步骤八、图像融合：对步骤六和步骤七得到的CT断层图像和XLCT荧光断层图像根据基            于标记点的刚性配准方法进行配准并融合，得到融合后的图像。           
3.如权利要求1所述的CT和X射线激发荧光双模同步断层成像方法，其特征在于，所述            图像总变差最小化的少角度CT重建方法包括以下步骤：           
第一步，采用锥束滤波反投影算法得到初始图像u<sub>0</sub>：对CT待重建的M幅投影采用滤波反            投影算法重建，将重建后图像u<sub>0</sub>作为初始图像；           
第二步，令k＝0，k为迭代次数(k＝0,1,2…K-1)，将u<sub>0</sub>作为迭代初始图像u<sup>(0)</sup>，采用代数            迭代法(algebraicreconstructiontechnique，ART)更新图像：           
<maths><math><mrow><msubsup><mi>u</mi><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow></msubsup><mo>=</mo><msubsup><mi>u</mi><mi>j</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>+</mo><mi>λ</mi><mfrac><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>-</mo><msubsup><mo>Σ</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></msubsup><msubsup><mi>u</mi><mi>n</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow></msubsup><mo>·</mo><msub><mi>m</mi><mrow><mi>i</mi><mi>n</mi></mrow></msub><mo>)</mo></mrow><mrow><msubsup><mo>Σ</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mi>N</
mi></msubsup><msubsup><mi>m</mi><mrow><mi>i</mi><mi>n</mi></mrow><mn>2</mn></msubsup></mrow></mfrac><msub><mi>m</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>;</mo></mrow></math></maths>

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。