[发明专利]图像处理方法

说明书

技术领域

本发明实施例一般涉及图像处理领域，更具体地说，涉及图像处理方法。 

背景技术

断层扫描成像的核心是将扫描物体时采集的一维投影数据经算法重建出目标物体的二维图像。如图1所示，扫描仪在扫描目标物体时，记录关于该物体成像信息关于不同视角(θ)及不同距离(d)的线性积分值。一个投影即是在同一视角方向θ₀([0，π])而不同距离d所测得的一组线性积分值。视角θ与距离d之间的数学关系式为： 

d＝xcosθ+ysinθ            (1) 

其中，x、y分别代表笛卡尔坐标系中横、纵坐标值。 

如图2所示，把所有视角θ的投影汇总到一起，则形成一张由不同频率的正弦曲线组成的正弦图(sinogram)。 

在断层扫描显像中，图像重建(image reconstruction)是指从正弦图中重建出被扫描的目标物体图像的过程。滤波反投影算法(FilteredBackprojection，FBP)是一种最常用的解析重建方法，在Avinash Kak &Malcolm Slaney(1988)，Principles of Computerized Tomographic Imaging，IEEEPress，ISBN 0-87942-198-3中有具体的论述，该文件通过引用方式而全文包括在本文中。滤波反投影算法的重建过程可大致分为四个步骤：1)将每个时域(或笛卡尔坐标系)的投影经傅里叶变换(Fourier transform)至频域；2)将转换的频域数据乘以一个斜坡滤波器(Ramp filter)；3)将乘得的数据经过反傅里叶变换回时域。由于数据与滤波器在频域内相乘，即等价于在时域内进行卷积计算，此时在时域的投影，即为滤波后的投影(Filteredprojection)；4)将滤波后的每个投影做反投影(back projection)，并将反投影相加，得到目标物体的图像f(x，y)。由于滤波反投影算法中涉及傅里叶变换，在重建图像里往往会出现条状伪影(streak artifacts)。 

条状伪影限制了滤波反投影算法在临床的一些应用。例如，在计算机断层扫描(CT)显影时，若扫描物内含有金属移植物，如髋关节假体或牙齿填充物，则会导致严重的条状伪影，对于临床诊断带来干扰及误导。对于核医学影像，放射性示踪剂(radioactive tracer)通常会在某些个别组织器官(如膀胱)累积，从而在该区域产生压倒性的高强信号。这种与常规信号强烈对比所形成的条状伪影，破坏周围腹部的重建成像，大大增加了诊断前列腺癌或宫颈癌的难度，也给信号定量分析带来困难。 

之前克服条状伪影的研究大多集中在两个方面：或是在相邻的两个投影数据中进行线性或更高阶的插值；或是使用及改进诸如期望最大化等迭代重建算法。两种方法均带来一定效果。但是，在大部分的插值算法得到的重建图像中，包围在高强信号区域的图像仍然存有强烈的难以修复的扭曲变形；而在迭代算法中，因算法自身存在非负约束(non-negative constraint)，得到的重建图像存有很高的偏离率(high bias)及较低的定量鲁棒性。因此，在一些对于生理学定量分析有很高要求的核医学应用中，迭代重建算法并不适用。 

发明内容

本发明实施例一方面的目的是降低条状伪影的干扰，提高重建图像质量。 

一方面，本发明实施例提出了一种图像处理方法，所述方法包括： 

(a)采用滤波反投影算法对断层成像扫描获得的目标物体的L帧动态正弦图y(t)进行图像重建，得到L帧动态图像x′(t)，其中t∈[0，L-1]； 

(b)利用归类算法将所述L帧动态图像x′(t)的信号强度归类为预设的K类，得到每类的平均强度曲线其中K＞1，k∈[1，K]； 

(c)根据所述将所述动态正弦图y(t)中的单帧正弦图y(t₀)按照信号强度分解成K幅次正弦图y′_p(t₀)，其中t₀∈[0，L-1]、p∈[1，K]，所述y′_p(t₀)表示所述动态正弦图y(t)中的单帧正弦图y(t₀)上信号强度为的区域；