基于DSPCT图像重建加速方法实现.pdfVIP

研讨会·论文集 全【日射线数字成像与cT新技术 基于DSP的CT图像重建加速方法实现 王超，李建新，王大会，闫镔凋，李磊 (国家数字交换系统工程技术研究中心，郑州450002) 摘要：随着cT技术的发展，重建速度已成为CT系统的重要指标之一，图像重建加速的研究也 成为CT研究的热点之一。DSP由J：采用哈佛结构，』￡自．专门的硬件乘法器以及J“泛采用流水线 操作等特点，已，I：始应用于图像重建加速。本文针对近年来开始应川的■维锥束cT的图像重建 算法进行了DSP加速方法研究。首先介绍了DSP应Hj于图像垂建的研究现状，其次，通过对三 维重建的典型算法一FDK的算粒分解分析，结合DSP的硬件结构的特点，提}l；了一种基’J二DSP 的CT图像币=建的综合加速技术。该方法通过C的优化、汇编优化和编译优化的综合使用，充分 利川了硬件结构和软件流水技术，采JlJTMS320C6455芯片实现了CT图像重建的加速计算，实验 结果表明，该方法取得了较岛加速比。 关键词：cT图像；三维重建；FDK：DSP：软件流水线 X射线计算机断层成像技术(Computed 令人瞩同，该技术最初应用于医学，极大推动了临床医学的发展，随后应用于工业领域， 在工业无损检测与无损评价方面也显示出了独特的优点…。作为cT的核心技术之一，同前 i维图像重建算法理论已有较好的丛础，实际的应用难点主要集中在计算规模与计算效率 上。工业CT图像■维重建的数据晕大，计算帚大，数据处理需要消耗大量时问，成为整个 CT系统的瓶颈制约幢1。囚此，CT图像重建速度成为衡量CT系统的重要指标之一。 针对CT图像三维重建的实时性要求，重建加速技术的研究主要集中在算法的改进，机 和数据分开的哈佛结构，具有专门的硬件乘法器以及广泛采用流水线操作，基于DSP的cT 图像重建加速引起研究人员的关注。 5 速度达到0．000 与汇编混合编程，以及将FFT函数拆成反序、核心计算等几个小函数，在DSP上比CPU计 算速度提高10倍。2004年陈雪松等…设计了一种丛于单指令集，多数据流并行处理结构的 可扩展的多数字信号处理器并行计算系统模犁，提出了DSP数最的选择依据，使用多DSP Mvoxel／s。 模为6002的cT图像重建，计算速度达到1．3170 基金项目：陶家863 省肇础与前沿技术研究计划项Il(072300450240)． 下超等：摹于DSP的CT图像重建加速方法实现 13 为了进一步加速锥束cT图像三维蕈建，本文基于DSP给出一种综合优化的加速方法。 使用C的优化、汇编优化和编译选项优化相结合的方法，充分发挥硬件结构和软件流水技 加速效能。 1 CT图像重建算法简介 在各种三维锥柬CT图像霞建算法中，FDK算法是一种基于圆轨道扫描的l-tj-维投影进 行三维蕈建的近似莺建算法，fll于该算法计算效率高，成为实际锥束重建中应用最为广泛 的算法b1。本文丰要讨论FDK算法。FDK算法的几何关系如图l所示。 罔1 FDK几何关系示意罔 ofFDK Fig．1diagram geometry 其中。D为射线源到旋转中心的距离，D+Y’为射线源到待重建点的距离在旋转后的 J，’轴上的投影。g(“，l，，0)为投影图像，f(x，Y，z)为重建后图像。 FDK算法的公式可以如下表示： ㈩ 脚国=三r(南]Ⅷ一卜们，赢斋g(u。-u)妇 其中，1——兰一为对投影数据q(u，，，秒)进行修正的加权囚子，适当地修正体素√D2+U2+V2 到源点的距离和角度幕。g(u)为常用的Ramp滤波函数，对不同投影角度的投影数据进行 水平方向上的一维滤波。外层积分的上下限为[O．2冗]，代