[发明专利]一种基于稀疏表示和非局部相似的核磁共振图像重建方法

说明书

技术领域

本发明属于图像处理技术领域，涉及一种核磁共振图像重建方法，具体是一种基于稀疏表示和非局部相似的核磁共振图像重建方法。

背景技术

磁共振成像是一种重要的医学成像技术，在临床中有着广泛的应用。近年来的研究表明，压缩感知理论能够利用相比其它方法更少的采样数据较好重建图像。该理论指出在重建磁共振图像的问题中，如果图像在某个确定的变换域中能够稀疏的表示，那么利用核磁共振图像对应的频域下采样数据能够以很高的概率完整重建原始图像。一种重建图像的方法是将图像的重建问题转化为凸优化问题。已有的采用全变差约束和稀疏域l¹范数约束的凸优化重建方法，对纹理结构复杂的图像重建效果不佳，图像的细节很难重建。通过利用非局部相似性先验信息对图像的稀疏表示系数进行约束，可以降低重建图像中的块效应，改善图像中纹理结构的重建效果。进而降低所需的采样数据量，具有实际意义。

发明内容

本发明的目的是提供一种充分利用图像非局部相似性的核磁共振图像重建方法，使得能够有效重建图像中的各种结构，提高重建质量。

为了实现上述目的，本发明提供的技术方案如下：

首先，对磁共振图像对应的傅里叶变换系数，采用变密度随机下采样的方式采样，由采样得到的数据通过傅里叶反变换得到用于重建的初始参考图像；其次，对初始参考图像进行分块，依据图像子块中的边缘和结构对图像子块进行分类，并得到每类图像子块对应的字典和图像子块的稀疏表示系数；最后，利用非局部相似性先验信息引入针对图像稀疏系数的约束项，结合图像在小波域的稀疏性，通过混合正则项求解模型对图像进行迭代重建。

具体实现方案包括以下步骤：

步骤(1)获取用于重建的初始参考图像，具体是：

对磁共振图像对应的傅里叶变换系数，采用变密度随机下采样的方式采样，即对傅里叶系数对应图像低频信息的部分较多的采样，对傅里叶系数对应图像高频信息的部分较少采样；对得到的采样数据矩阵缺失的部分补零值，然后用二维傅里叶反变换得到用于重建的初始参考图像x⁽⁰⁾。

步骤(2)对参考图像分块并对图像子块分类，求得每类对应的字典和图像子块对应类字典的稀疏系数，具体是：

将输入图像x＝x⁽⁰⁾进行分块，即利用提取图像子块的矩阵P_i通过式x_i＝P_ix从x中得到大小为n×n的图像子块x_i，总的分块数为Π，其中对图像子块进行高通滤波，然后使用K均值算法对图像子块进行分类，得到具有相似结构特征的K类图像子块；利用主成分分析法(PCA)求得每一类的字典，由K类图像子块对应的K个子字典构成整个图像的字典；对于每一图像子块x_i选择该块所属类k对应的子字典则图像子块x_i对应的稀疏系数可以通过求解l¹范数最小问题得到，如式(1)所示：

其中常数λ表示正则化参数，由此，得到对应的稀疏表示为

步骤(3)对图像子块稀疏系数进行约束，迭代重建磁共振图像，具体是：

对于每个图像子块x_i，以x_i的中心为中心、大小为S×S的范围内搜索其非局部相似的块，得到x_i所有非局部相似图像子块的集合表示为C_i，Q(C_i)为该集合中的元素个数；利用图像的非局部相似性对重建图像进行估计，令x_i,j为x_i的非局部相似的图像子块，其稀疏系数表示为则重建图像子块x_i对应的稀疏系数的无偏估计表示为可通过式(2)计算：

其中，ω_i,j是权重，可以通过式(3)计算，其中h为权重控制常数；