基于fpga的锥束ct中离散希尔伯特变换的设计与实现-微电子学与固体电子学专业论文.docxVIP

万方数据 万方数据 Classified Index: TN47 U.D.C: 621.3 Dissertation for the Master Degree in Engineering DISCRETE HILBERT TRANSFORM DESIGN AND IMPLEMENTATION FOR THE CONE-BEAM CT WITH FPGA Candidate： Tingliang Ma Supervisor： Prof.Yan ZHANG Academic Degree Applied for： Master of Engineering Specialty： Microelectronics and Physics Solid-State Affiliation： Shenzhen Graduate School Date of Defence： December, 2013 Degree-Conferring-Institution： Harbin Institute of Technology 哈尔滨工业大学工学硕士学位论文 哈尔滨工业大学工学硕士学位论文 摘 要 医用CT（Computed Tomography）图像重建技术中，用螺旋锥束扫描方式重建 技术以其扫面速度快、空间分辨率高的优点得到了越来越多的关注。但是螺旋锥 束CT三维重建理论比较复杂，计算量也非常大，重建速度是这一理论的瓶颈所在。 锥束CT三维重建的数据传输量和运算量巨大，只依靠软件来完成重建任务已经不 能满足医学和工程的需要。由软件和硬件协同完成三维图像重建的任务将会大大 减少CPU的负担，提高图像重建的速率，所以对锥束CT图像重建速度的研究具有 很高的应用研究价值。本文对锥束CT图像重建中离散希尔伯特变换（DHT）部分 做了细致研究，设计了FPGA实现DHT的硬件电路，使锥束CT图像重建算法的DHT 部分由FPGA来实现，其他部分由软件实现，提高重建速度。硬件电路采用并行化 结构，使 1024 点DHT运算可以在 56us完成；数据结构选用单精度浮点数的数据格 式，整个DHT硬件电路的运算精度可以达到 10-5；选用传输速度快的PCIE作为数 据传输总线，DHT接口电路实现了 4Gb/s的传输速度。本文研究内容主要包括四部 分： 首先，比较 DHT 的几种实现方式的硬件实现的速度和面积，根据项目的实际 情况选择用快速傅里叶变换（FFT）实现 DHT；其次，根据按频率抽取基 2FFT 算 法原理设计 DHT 的蝶形流程图，根据蝶形流程图提出了单个运算单元结构和多个 运算单元并行结构的 DHT 硬件电路结构；再次，利用 FPGA 的单精度浮点加法器 和单精度浮点乘法器 IP 核，完成对 DHT 硬件电路的实现，对 DHT 电路性能进行 了分析比较，具有较好的速度优势；最后，根据 PCIE 的传输协议，对 DHT 硬件 电路进行了改进，使 DHT 硬件电路的接口与 PCIE 总线的接口完全契合，在 Xilinx 公司的 XC5VLX50T FPGA 器件的 ML555 开发板上，以 PCIE 为数据总线实现了 锥束 CT 图像重建中 DHT 的硬件加速。并且用 matlab 做了 FPGA 的 DHT 运算结 果的验证界面，验证了 DHT 硬件加速系统 DHT 运算的正确性。 本文设计的 DHT 的 FPGA 硬件电路的性能与其他文献的结构进行了比较，做 一次 DHT 运算的时间比 Altera 的 IP 核计算时间少 58us，具有很好的实时性。本 文中 DHT 硬件结构的最大特点是可以和 PCIE 总线接口对接，实现通过 PCIE 总线 控制 FPGA 的 DHT 运算。整个 DHT 硬件加速系统实现了由 PC 机通过 PCIE 总线 控制 FPGA 实现 DHT 运算的功能，为以后实现软件和硬件共同完成锥束 CT 图像 重建任务打下了坚实的基础。 关键词：DHT；单精度浮点加法器；FFT；PCIE I - Abstract Spiral cone-beam CT (Computed Tomography) scanning surface reconstruction technique has the advantages of high sweep speed and high spatial resolution, which gains more and more attention. But the spiral cone-beam CT reconstruction theory is more complex and the amount of calculation is very large, so