[发明专利]基于非线性扩散模型的有限角CT图像重建方法及装置

说明书

本申请提供了基于非线性扩散模型的有限角CT图像重建方法及装置，步骤如下：利用有限角CT扫描数据集，通过与扫描几何参数集G相关的图像重建算子对候选图像进行更新后得到图像；根据图像中像素点的像素值，确定图像的边缘信息；对图像的边缘点进行赋值后得到边缘指示符函数；根据边缘指示符函数，利用二维非线性扩散模型对图像进行更新后得到图像；判断，或是否达到迭代次数上限，若满足，则终止迭代，输出图像。采用上述方法，以提高重建得到的CT图像的准确性。

技术领域

本发明涉及CT图像重建技术领域，具体而言，涉及基于非线性扩散模型的有限角CT图像重建方法及装置。

背景技术

X射线计算机断层成像 (Computed Tomography, 简称CT) 在医学、工业、安检、考古等领域有广泛的应用。受 X 射线高吸收、被扫描物体的结构特点、设备限制等很多因素的影响，通常无法做到全角度的扫描，导致采集不到全角度的数据，所以人们非常重视有限角 CT 图像重建问题。

有限角 CT 图像重建问题，其过程是通过利用有限角度扫描的投影数据，对被扫描物体的断层图像进行重建；发明人在研究中发现，在使用传统的 CT 图像重建算法重建图像时，由于所获取的投影数据为非全角度的投影数据，所以重建得到的图像会出现伪影，这些伪影会造成图像无法准确体现出被扫描物体的结构，即，降低了重建得到的CT图像的准确性。

发明内容

有鉴于此，本发明的目的在于提供基于非线性扩散模型的有限角CT图像重建方法及装置，以提高重建得到的CT图像的准确性。

第一方面，本申请实施例提供了一种基于非线性扩散模型的有限角CT图像重建方法，所述方法包括：

步骤1，输入已知数据集：有限角CT扫描获得的数据集p，CT扫描几何参数集G；

步骤2，初始化：初始候选图像，确定迭代次数上限/和迭代终止阈值/；

步骤3，迭代处理候选图像，/=0，1，2，···，/，利用有限角CT扫描数据集，通过与扫描几何参数集G相关的图像重建算子/对候选图像/进行更新后得到图像，其中，/表示扫描几何参数集G相关的图像重建算子；

步骤4，根据图像中的每个像素点在预设八个方向上相邻的八个像素点的像素值，确定出图像/的边缘信息；

步骤5，根据预设的赋值规则对图像的边缘信息中的边缘点进行赋值后得到边缘指示符函数；

步骤6，根据所述边缘指示符函数，利用预设的二维非线性扩散模型对图像进行更新后得到图像/；

步骤7，判断相邻两次迭代图像间的差别是否小于所述迭代终止阈值，即，或是否达到迭代次数上限/，若不满足，转至步骤3，若满足，则终止迭代，输出图像/。

可选地，所述步骤3中，在利用有限角CT扫描数据集，通过与扫描几何参数集G相关的图像重建算子/对候选图像/进行更新后得到图像/前，所述方法还包括：

将候选图像的投影数据进行随机化处理；

在利用有限角CT扫描数据集，通过与扫描几何参数集G相关的图像重建算子/对候选图像/进行更新后得到图像/后，所述方法还包括：

将图像中的负像素赋值为0。

可选地，所述步骤4中，所述根据图像中的每个像素点在预设八个方向上相邻的八个像素点的像素值，确定出图像/的边缘信息，包括：

对于图像中的每个像素点，根据该像素点在预设八个方向上相邻的八个像素点的像素值确定出该像素点的梯度值；

判断该像素点的梯度值是否超过预设的标准梯度值；

若该像素点的梯度值超过所述标准梯度值，则将该像素点标记为边缘点；