基于小波变换的多尺度锥束CT图像快速三维重建方法技术

一种基于小波变换的多尺度锥束CT图像三维重建快速算法，包括以下步骤：(1)对试件进行圆周锥束CT扫描，采集一组用于重建的投影图像；(2)根据对重建数据分辨率要求的不同，对采集到的投影图像进行相应尺度的小波变换，得到各尺度小波分解系数；(3)选择相应尺度的小波系数进行FDK重建；如果需要对小波一次分解后的所有分量分别进行FDK重建，则继续；(4)将得到的重建数据分别沿着y轴方向取断层图像，进行相应的小波逆变换；(5)将小波逆变换后得到的三维数据沿着z轴方向取断层图像，进而进行线性插值，得到最终重建图像。本发明专利技术能够在一定的精度内快速实现基于FDK算法的三维重建，在此过程中能够针对不同的要求来呈现不同分辨率的三维数据。

【技术实现步骤摘要】

本专利技术涉及一种用于计算机图像处理的基于小波变换的多尺度锥束CT图像的三维重建方法，属于计算机图像处理

技术介绍
X光计算机断层成像技术(X-ray Computer Tomography,XCT)，是利用X射线对物体进行不同角度的扫描而获得物体断层信息的成像技术，该技术涵盖了物理学、数学、计算机视觉、计算机图形学、数字信号处理等多方面的学科知识，并且它被广泛应用在很多科学领域，这使得人们观察物体内部结构的能力极大地增强。尤其是在医学临床诊断中，它被用来当做一种获取人体内部结构信息的最佳手段。目前主流的商用CT机是扇束CT扫描机，也就是说这种CT机的X射线源每扫描一周只能获取被检测对象某一个断层的数据，这样导致X射线的使用效率低、扫描速度慢及重建速度慢。与扇束CT机相比，锥束CT由于采用了锥束X射线源和面阵探测器来采集被检测对象不同角度的投影数据，大大提高了X射线利用率、扫描速度和CT重建图像的分辨率，在医学临床检测和工业无损检测上的应用越来越广泛。一般地，锥束XCT重建算法分为两类，一类是迭代法，另一类是解析法(如FDK算法)。迭代法的计算量很大、效率低，难以适应现代商业实时性的要求。与迭代法相比，解析算法的重建效率高，需要的数据存储空间小。一般地解析法分为近似重建法和精确重建法，其中FDK算法是从扇束滤波反投影算法推导出来的锥束XCT近似重建算法，在X射线源锥角较小的情况下重建出来的断层图像质量非常好。由于其简单而且易实现的特点，使FDK算法的研究具有实际应用价值。锥束CT完成一次三维扫描的数据量达几百GB甚至TB，重建算法复杂，运算量大，耗时较多，难以满足实际应用中实时性的要求，特别是在医学应用领域。其次，针对不同的应用和显示载体的限制，对投影数据重建后的三维数据的分辨率的要求是不同的，对一些重建后分辨率要求不高的情况，不需要耗时太多来重建高分辨率图像。
技术实现思路
针对锥束CT三维重建算法存在的耗时多以及不同的应用对重建后的三维数据的分辨率要求不同等问题，本专利技术提出在保证精度前提下快速实现三维重建的基于小波变换的多尺度锥束CT图像三维重建快速算法。本专利技术的基于小波变换的多尺度锥束CT图像三维重建快速算法，包括以下步骤：(1)对试件进行圆周锥束CT扫描，采集一组用于重建的投影图像；(2)根据对重建数据分辨率的要求，对采集到的投影图像进行相应尺度的小波变换，得到各尺度小波分解系数；(3)选择相应尺度的小波系数进行FDK重建；已知得到分辨率为2N×2N的投影图像，对投影图像进行小波一次分解，进而对分解后的低频分量进行FDK重建，最终得到分辨率为(N/2)×(N/2)×(N/2)的重建数据；或者对小波一次分解后的所有分量分别进行FDK重建，得到四组重建数据，继续执行步骤(4)和(5)；以一次小波分解为例，由步骤(1)得到一组分辨率为2N×2N的二维投影图像，对投影图像进行小波一次分解，得到四组分辨率为N×N的小波分解系数，若对重建数据分辨率要求较低，可仅对低频分量进行FDK重建，得到分辨率为(N/2)×(N/2)×(N/2)的三维重建数据；若对重建数据分辨率要求较高，可对小波分解后的所有分量分别进行FDK重建，得到四组重建数据，继续执行步骤(4)和(5)。(4)将得到的重建数据分别沿着y轴方向取断层图像，进行对应的小波逆变换；(5)将小波逆变换后得到的三维数据沿着z轴方向取断层图像，进行线性插值，得到最终重建图像。本专利技术的方法能够在一定的精度内快速实现基于FDK算法的三维重建，在此过程中能够针对不同的要求来呈现不同分辨率的三维数据。附图说明图1是一幅图像在二维频域被分解为四个子带的框图。具体实施方式以下详细说明本专利技术基于小波变换的多尺度锥束CT图像三维重建快速算法的具体过程。(1)对试件进行圆周锥束CT扫描，采集一组用于重建的投影图像；锥束扫描及投影数据的参数如下表所示：参数名称参数射线源到试件中心距离(mm)DSO射线源到探测器距离(mm)DSD探测器长度(mm)2L探测器分辨率2N×2N重建物体分辨率N×N×N重建物体边长长度(mm)L采集投影间隔(°)M采集投影个数360/m按照上述参数进行锥束扫描，会得到360/m张，分辨率为2N×2N的二维投影图像Pn(u,v)，n＝1,2,…,360/m，u,v分别标记投影图像像素的行标和列标，它们的零点都以平面探测器的中心像素为基准。(2)根据对重建数据分辨率的要求，对采集到的投影图像进行相应尺度的二维离散小波变换，得到各尺度小波分解系数；二维离散小波变换：二维小波变换由一维小波变换扩展而来，由一维尺度函数和相应的小波函数ψ，可得到一个二维尺度函数和三个二维小波ψH(x,y)，ψV(x,y),ψD(x,y)：ψD(x,y)＝ψ(x)ψ(y)，例如，Haar尺度函数定义为:小波函数定义为:那么，尺寸为M×N的函数f(x,y)的离散小波变换是： W ψ i ( j , m , n ) = 1 M N Σ x = 0 M - 1 Σ y = 0 N - 1 f ( x , y ) ψ j , m , n i ( x , y ) , i = { H , V , D本文档来自技高网...

【技术保护点】
一种基于小波变换的多尺度锥束CT图像三维重建快速算法，其特征是：包括以下步骤：(1)对试件进行圆周锥束CT扫描，采集一组用于重建的投影图像；(2)根据对重建数据分辨率要求的不同，对采集到的投影图像进行相应尺度的小波变换，得到各尺度小波分解系数；(3)选择相应尺度的小波系数进行FDK重建；已知得到分辨率为2N×2N的投影图像，对投影图像进行小波一次分解，进而对分解后的低频分量进行FDK重建，最终得到分辨率为(N/2)×(N/2)×(N/2)的重建数据；或者对小波一次分解后的所有分量分别进行FDK重建，得到四组重建数据，继续执行步骤(4)和(5)；(4)将得到的重建数据分别沿着y轴方向取断层图像，进行相应的小波逆变换；(5)将小波逆变换后得到的三维数据沿着z轴方向取断层图像，进而进行线性插值，得到最终重建图像。

【技术特征摘要】
1.一种基于小波变换的多尺度锥束CT图像三维重建快速算法，其特征是：包括以下步骤：(1)对试件进行圆周锥束CT扫描，采集一组用于重建的投影图像；(2)根据对重建数据分辨率要求的不同，对采集到的投影图像进行相应尺度的小波变换，得到各尺度小波分解系数；(3)选择相应尺度的小波系数进行FDK重建；已知得到分辨率为2N×2N的投影图像，对投影图像进行小波一次分解，进而...

【专利技术属性】
技术研发人员：韩民，成旭，
申请(专利权)人：山东大学，
类型：发明
国别省市：山东;37