专利名称：：一种锥束ct三维图像的并行重建方法
技术领域：
：本发明属于CT系统图像重建领域，涉及锥束CT系统并行化三维图像重建的方法。
背景技术：
：CT(ComputedTomography)技术利用射线(一般为X射线)透射被检测物体，并在探测器上获取一组投影图像，再配合相应的重建算法得到物体的切片图像。目前研究和应用的CT可分为二维CT和三维CT两大类，锥束CT属于三维CT。与传统二维CT相比，锥束CT扫描速度快、射线利用率高，能获得均匀、高精度的空间分辨率。锥束CT的重建是三维图像重建，即需要重建得到连续的很多层切片图像，涉及巨大的计算量。在锥束CT实际应用中，在保证良好的重建质量的基础上，提高重建速度成为必须突破的瓶颈。目前在商用锥束CT系统中，应用最广泛的重建算法是FDK滤波反投影算法，该算法是一种近似重建，计算量相对较小，在小锥角的情况下可以获得质量良好的重建图像。即便如此，FDK算法的计算量依然较大，主要集中在反投影过程，计算复杂度为0(N4)，其中N为投影数据的尺寸。例如，用360幅5122投影图像重建5123图像，在普通PC机上通常需要计算数小时才能完成。在现有的文献中所提到的提高图像重建速度的方法，主要可以分为以下三大类(1)重建算法的改进，如几何变量分离法、增量法、图像细分法等；(2)采用特殊的硬件来加速，如采用专用的DSP模块，对反投影计算过程进行硬件化实施；(3)并行计算技术，如采用阵列处理器或分布式计算机对图像重建过程进行并行化。毛海鹏、张定华、梁亮等人在《系统仿真学报》(2004，16(ll):2486-2489)的文章"一种基于PC的快速三维图像重建方法"中提出了一种Z线优先重建算法，它是对传统FDK算法进行改进并通过数据并行计算来实现快速三维图像重建的算法，其基本原理是根据FDK反投影过程中反投影点的p值与z无关，重建图像空间的每条Z线可映射为滤波投影平面上的一条u向亚像素线，标记为L，L可由两条相邻的u线进行线性插值得到，最后根据L线与Z线的对应关系得到重建图像；在此基础上通过有效地组织和划分重建数据，使得对重建数据的内存访问非常连续，采用单指令多数据(SingleInstructionMultipleData,SIMD)技术进行数据并行处理。该方法的计算效率很高，但只实现了单线程的计算，而且不能通过简单增加线程的方法进行并行计算，其原因在于对于同一个重建图像空间，简单多线程计算可能出现多个线程同时操作一个重建体素的情况，此时可能出现意想不到的结果；由于Z线优先重建算法本身已经进行了SIMD并行处理，而一般并行计算的线程通信与数据保护技术很难应用于SIMD计算过程。随着CPU技术的发展，目前双核CPU和四核CPU已经普及，而六核CPU、八核CPU也将很快得到广泛应用。因此，研究并实现Z线优先重建算法在多核计算机中的并行计算方法，是一个具有应用基础和实际价值的项目。
发明内容为了克服现有技术不能通过简单增加线程的方法进行并行计算的不足，本发明提供一种锥束CT三维图像的并行重建方法，使得Z线优先重建算法能够直接在多核计算机中实行并行计算。本发明解决其技术问题所采用的技术方案包括以下步骤(1)对试件进行圆周锥束CT扫描，采集一组用于重建的投影图像；(2)在投影图像中选取一个边长为￡个象素的正方形区域，将该组所有投影图像按该正方形区域裁剪为边长为￡个象素的正方形图像，并将所有裁剪后的投影图像进行公知的对数运算，此时投影图像的象素灰度为单精度浮点型；(3)对该组投影图像进行FDK算法中的滤波处理，滤波函数采用公知的S-L滤波器、R-L滤波器或SL-W滤波器；(4)从计算机中一次分配P大小的单精度浮点型内存空间，此内存空间称为重建空间，是一个边长为￡个体素的立方体；(5)对重建空间进行分割，得到若干个大小基本相同的子重建空间；(6)对各个子重建空间同时采用SIMD加速的Z线优先重建算法进行重建，当所有子重建空间完成其重建计算时，重建结束；(7)根据需要将重建空间存储为公知的X向、Y向或Z向的序列切片图像。在上述方法第5步中，为减少不必要的重建计算，同时达到避免线程通信和实现负载平衡的目的，进行重建空向分割的具体步骤如下(图2):1)选取立方体重建空间的内切圆柱为感兴趣重建空间ROI，即每一层切片的重建区域只局限于切片图像的内切圆中，这样可将重建体素的数目减少至原始数目的7l/4;2)设置2S子重建空间的个数^计算机中所有CPU的总核心数；3)对感兴趣重建空间ROI按子重建空间个数进行扇形等分，对位于等分线上的重建空间体素，将其归入包含该体素较多部分的子重建空间，这样得到的各子重建空间的大小基本相同，且各子重建空间之间无交义重叠部分。在上述方法第6步中，所有子重建空间同时采用SIMD加速的Z线优先重建算法进行重建的具体步骤如下1)分别为每一个子重建空间设置一个整数标志位，O表示空闲或已完成重建计算，1表示正在进行重建计算，并设置所有标志位的初始值为0;2)判断是否还有未读入的滤波处理后的投影图像，若是，则继续进行下一步，若否，则全部重建结束；3)读入下一幅滤波处理后的投影图像；4)设置所有标志位的值为1，同时从该滤波图像中读取所需的数据进行各子重建空间的重建，由于对滤波图像仅进行读操作，因此不会出现子重建空间重建任务间的排斥和滤波图像数据被改写的情况。另外，由于操作系统的原因，实际中不可能所有子重建空间同时完成各自的重建计算，因此当某个子重建空间完成其重建计算时，将其标志位设置为0;5)判断直到所有标志位的值均为0，转第2步。本发明的有益效果是在保持Z线优先重建算法原有的重建图像质量基础上，通过对重建空间进行简便实用的分割，同时达到避免线程通信和实现负载平衡的目的，解决了Z线优先重建算法不能直接在多核计算机中实行并行计算的问题。下面结合附图和实施例对本发明进一步说明。图l是本发明方法流程图。图2是三维图像重建中的Z-L线映射关系及重建空间按线程分割方式。具体实施例方式在IntelXeon2.5GHz四核处理器、4G内存的计算机上，实施本发明方法的步骤如下(1)对一阶梯铝柱零件进行圆周锥束CT扫描，采集一组360幅768x%0的投影图像用于重建-,(2)在投影图像中选取一个边长为512个象素的正方形区域，该正方形区域的中心位于整幅投影图像的中心，将该组所有投影图像按该正方形区域裁剪为边长为512个象素的正方形图像，并将所有裁剪后的投影图像进行公知的对数运算，此时投影图像的象素灰度为单精度浮点型；(3)对该组投影图像进行FDK算法中的滤波处理，滤波函数采用公知的S-L滤波器；(4)从计算机中一次分配5123大小的单精度浮点型内存空间，此内存空间称为重建空间，是一个边长为512个体素的立方体；(5)对重建空间进行分割，得到若干个大小基本相同的子重建空间，为减少不必要的重建计算，同时达到避免线程通信和实现负载平衡的目的，进行重建空间分割的具体步骤为1)选取立方体重建空间的内切圆柱为感兴趣重建空间ROI，即每一层切片的重建区域只局限于切片图像的内切圆中，这样可将重建体素的数目减少至原始数目的兀/4;2)设置子重建空间的个数二计算机中所有CPU的总核心数=4;3)对感兴趣重建空间ROI按子重建空间个数进行扇形4等分，对位于等分线上的重建空间体素，将其归入包含该体素较多部分的子重建空间，这样得到的各子重建空间的大小基本相同，且各子重建空间之间无交叉重叠部分。(6)对各个子重建空间同时采用SIMD加速的Z线优先重建算法进行重建，当所有子重建空间完成其重建计算时，重建结束，具体步骤如下1)分别为每一个子重建空间设置一个整数标志位，O表示空闲或已完成重建计算，1表示正在进行重建计算，并设置所有标志位的初始值为0;2)判断是否还有未读入的滤波处理后的投影图像，若是，则继续进行下一步，若否，则全部重建结束；3)读入下一幅滤波处理后的投影图像；4)设置所有标志位的值为1，同时从该滤波图像中读取所需的数据进行各子重建空间的重建，由于对滤波图像仅进行读操作，因此不会出现子重建空间重建任务间的排斥和滤波图像数据被改写的情况。另外，由于操作系统的原因，实际中不可能所有子重建空间同时完成各自的重建计算，因此当某个子重建空间完成其重建计算时，将其标志位设置为0;5)判断直到所有标志位的值均为0，转第2步。(7)根据需要将重建空间存储为公知的Z向序列切片图像。针对上述裁剪后的投影图像，表1给出了FDK算法、Z线优先算法和本发明方法在锥束CT三维图像重建的反投影部分的计算速度比较，可见本发明方法很好地实现了Z线优先算法的并行加速计算，而且由于没有改变Z线优先算法本身的计算策略，所以重建图像质量与Z线优先算法保持一致。表l三维图像重建的反投影部分的计算速度比较<table>tableseeoriginaldocumentpage8</column></row><table>权利要求1、一种锥束CT三维图像的并行重建方法，其特征在于包括下述步骤(1)对试件进行圆周锥束CT扫描，采集一组用于重建的投影图像；(2)在投影图像中选取一个边长为E个象素的正方形区域，将该组所有投影图像按该正方形区域裁剪为边长为E个象素的正方形图像，并将所有裁剪后的投影图像进行对数运算，此时投影图像的象素灰度为单精度浮点型；(3)对该组投影图像进行FDK算法中的滤波处理，滤波函数采用S-L滤波器、R-L滤波器或SL-W滤波器；(4)从计算机中一次分配E3大小的单精度浮点型内存空间，此内存空间称为重建空间，是一个边长为E个体素的立方体；(5)对重建空间进行分割，得到若干个大小基本相同的子重建空间；(6)对各个子重建空间同时采用SIMD加速的Z线优先重建算法进行重建，当所有子重建空间完成其重建计算时，重建结束；(7)根据需要将重建空间存储为公知的X向、Y向或Z向的序列切片图像。2、根据权利要求1所述的一种锥束CT三维图像的并行重建方法，其特征在于所述的第(5)步中，进行重建空间分割的具体步骤如下1)选取立方体重建空间的内切圆柱为感兴趣重建空间ROI，即每一层切片的重建区域只局限于切片图像的内切圆中；2)设置2^子重建空间的个数^计算机中所有CPU的总核心数；3)对感兴趣重建空间ROI按子重建空间个数进行扇形等分，对位于等分线上的重建空间体素，将其归入包含该体素较多部分的子重建空间。3、根据权利要求1所述的一种锥束CT三维图像的并行重建方法，其特征在于所述的第(6)步中，所有子重建空间同时采用SIMD加速的Z线优先重建算法进行重建的具体步骤如下1)分别为每一个子重建空间设置一个整数标志位，0表示空闲或已完成重建计算，l表示正在进行重建计算，并设置所有标志位的初始值为0;2)判断是否还有未读入的滤波处理后的投影图像，若是，则继续进行下一步，若否，则全部重建结束；3)读入下一幅滤波处理后的投影图像；4)设置所有标志位的值为1，同时从该滤波图像中读取所需的数据进行各子重建空间的重建，当某个子重建空间完成其重建计算时，将其标志位设置为0;5)判断直到所有标志位的值均为0，转第2)步。全文摘要本发明公开了一种锥束CT三维图像的并行重建方法，对试件进行圆周锥束CT扫描，采集一组用于重建的投影图像；将该组所有投影图像按该正方形区域裁剪为边长为E个象素的正方形图像，并将投影图像进行对数运算；对该组投影图像进行滤波处理；从计算机中一次分配E<sup>3</sup>大小的单精度浮点型内存空间作为重建空间；对重建空间分割得若干个大小基本相同的子重建空间；对各子重建空间采用SIMD加速的Z线优先重建算法进行重建；根据需要将重建空间存储为公知的X向、Y向或Z向的序列切片图像。本发明解决了Z线优先重建算法不能直接在多核计算机中实行并行计算的问题。文档编号G06T1/60GK101520899SQ20091002191公开日2009年9月2日申请日期2009年4月8日优先权日2009年4月8日发明者昆卜,张定华,李明君,黄魁东申请人:西北工业大学