超级流水线三维成像CT系统.pdfVIP

簟10届垒暖计算机在现代科学技术镊墟应用拳术会设论盘量 yrI_脚佣co忡^ppII盟l胁hM喇咖S岫Te由·咖吖 P喊啪dn∞啊腑10廿T曲№S 超级流水线三维成像cT系统 刘 亚 强 (清华大学工程物理粟，北京，lo0084) 在Cr理论界，三堆Cr田像重建在算法的复杂度和计算量方面一直存在着巨大的困难．本文提出将三 雏Cr运算分为多级。采用局部同络进行并行分布i算．莨田像重建速度大为墨膏． 美■调。三堆CT田像重建 I，i n?一，f 一’}‘¨” 1黜 ‘丑：‘ 。l“¨ 7。l 三维CT图像重建的计算量远远超过了二维CT图像重建．为了缩短图像重建时间。人们 提出各种方法．算法改进是最经济也是最有效的办法．但现阶段如果投有理论上的突破．各种 办法也只是修修补补．使用专用硬件可以大大缩短图像重建时间．但其算法固定，不能适应 CT广泛应用的要求．并行运算在近十年有了巨大的发展，但由于其通用性的考虑。在图像重 建方面并没有独到之处，因而工业CT中没有得到广泛的应用．由于网络技术的飞速发展．特 别是局域网据传输率的大幅度提高．使效据的传输时同不再成为工作的主要矛后，使阿络并行 分布运算成为一种经济可行的办法．本文将三维cT运算分为多级．采用局郁阿络进行并行分 布运算，使图像重建速度大为提高． 2超级流水线三维CT系统 L．A．Feldkmp等人在80年代中提出的卷积反投影算法∞是最早进入实用的三维CT算 法tB．D．smith指出了三维cT精确置建的充分必要条件，并在文献[3]提出了相应的重建算 法．后来提出的TTCR算法采用一组空间不变的滤波器的线性组合来逼近空间变化的点扩散 函数．总之。绝大部分的三维cT解析重建算法主要包古滤波和反投影两个过程．文献[1]给出 了圆轨迹扫描的三维图像重建方式。文献[2]对其进行了分析并提出了相应的改进算法，总的 来说。实现锥柬卷积反投影重建需要经过三个过程．以传统的圆轨迹锥柬扫描(如图1)重建为 倒，其方法如下； 1)数据加权： 丸(y，z)皇P．(y，z)———生一 ∥万干7邪 2)卷积滤波t f r+∞，+∞ 。 凡(y，z)=Idy，I羽也(y，一y)玩(z’一z)^(r．z，) 其中t风ct，=嘉≯n(詈)一(争n罢)2，丑c幻=去(罢)． 100 3)加权反投影· m．n‘ ，6)=1万÷!-_示P．(y，z)∞ J●～口1-r’；J =．_ ：．o 其中，z(≯)掌孑：南d·y(声)昌孑：；-；}歹d· 对上述三个步骤的计算量进行了分析比较，投影数据的加权并不需要太大的计算量t卷积 的计算量比较大。但由于其二维卷积可分为两次一维卷积计，同时文献[4]表明·不进行z方 向的卷积对重建效果没有太大的影响，因此只稠下一次一维卷积运算．如果采用硬件卷积器卷 积运算还将减少一到两个数量级I加权反投影的运算量很大t文献[5]提出了锥柬三维CT快 速图像重建的坐标地址映射反投影方法．可以将其运算量减步一个敦量级．即使采用上述方法 减少计算量，三维CT的图像重煮时间仍然是很长的．仅以反投髟为饲t重建256x256×100 5．这仅仅是反投影的时问·不包括致据加权、 的图像，在‘．7MFIDPS的SGI工作站需要424 · Z 卷积以及数据的增密和插值运算． ／ 对于常甩的圆轨迹和黑旋线轨迹扫描的 ／，一 芗 三维圈像重建，其算法都可以分为三个步骤． 7