一种图像重建方法技术

一种图像重建方法，包括：读取原始数据；将原始数据降维以获得一维向量；对一维向量取自相关，以获得自相关矩阵；对自相关矩阵进行特征分解，以获得自相关矩阵的特征值和特征向量；将特征向量按照其对应的特征值大小进行排序；对前N个较小的特征值对应的特征向量构成的子空间取自相关，以获得目标子空间；定义相位编码方向上的谱搜索的格点位置；对于每个相位编码方向执行以下谱搜索步骤：定义对应于相位编码方向的阵列流形；获得对应于相位编码方向的阵列流形与目标子空间的相关系数或者获得阵列流形与目标子空间的逆矩阵的相关系数，然后对相关系数进行负相关变换以获得谱；从谱中提取元素以获得重建出的图像的一行；输出重建出的图像。建出的图像。建出的图像。

【技术实现步骤摘要】
一种图像重建方法


[0001]本专利技术涉及一种图像重建方法，具体地，涉及一种磁共振图像重建方法。

技术介绍

[0002]磁共振成像在医学领域中有着广泛的应用，通过磁共振成像技术能够对脑部、脊柱和心血管系统等人体器官进行成像。磁共振成像技术通过外部射频信号使得处于磁场中的待测样品发生磁共振现象，然后通过线圈接收磁共振信号，并将磁共振信号通过图像重建方法转化为图像数据，以获得待测样品的图像。
[0003]在现有技术中，绝大部分上述图像重建方法是基于傅里叶变换进行的，但是基于傅里叶变换的图像重建方法存在一些如下所述的局限性。一是当进行高分辨成像时，随着成像资源(主要是相位编码)的增加，所能换得的分辨率增加是迅速衰减的，即浪费成像资源。二是磁共振机器报告的名义分辨率(实际含义为每像素代表多少空间距离)很高，但实际上由于部分容积效应、瑞利极限等原因，导致看不清物体，即真实分辨能力不好。换句话说，由于信噪比的迅速衰减，傅里叶方法本身具有的旁瓣叠加的干扰问题，导致即使用特别多的相位编码数量所获得的图像也是模糊的图像。三是在待测样品的性质(例如，磁敏感性、质子密度)在不同区域剧烈变化的情况下，傅里叶变换方法将导致严重的伪影现象。上述局限性使得基于傅里叶变换的图像重建方法所获得的图像质量仍有待改进。
[0004]为此，本专利技术提出了一种图像重建方法，改进了上述现有技术中存在的问题。

技术实现思路

[0005]根据本专利技术的第一方面，提供了一种图像重建方法，包括：读取原始数据；将原始数据降维，以获得一维向量；对一维向量取自相关，以获得自相关矩阵；对自相关矩阵进行特征分解，以获得自相关矩阵的特征值和与相应特征值对应的特征向量，特征值和特征向量的数量均为M；将特征值按照大小进行排序，获得前N个较小的特征值对应的特征向量构成的子空间，其中N小于M；对前N个较小的特征值对应的特征向量构成的子空间取自相关，以获得目标子空间；定义相位编码方向上的谱搜索的格点位置；对于每个相位编码方向执行谱搜索步骤，谱搜索步骤包括：定义对应于相位编码方向的阵列流形；获得对应于相位编码方向的阵列流形与目标子空间的相关系数，然后对相关系数进行负相关变换以获得谱；从谱中以预定规则提取元素以获得重建出的图像的一行，预定规则使得谱的每一行和每一列均有且仅有一个元素被提取；输出重建出的图像。
[0006]根据该方案，通过谱搜索方法而非传统的傅里叶变换方法获得了图像，与传统的傅里叶方法相比，克服了一些傅里叶方法的缺点，改进了图像质量。
[0007]在一些方案中，定义对应于相位编码的阵列流形的步骤包括：定义相位编码方向上的原始数据的格点位置和频率编码方向上的原始数据的格点位置，其中相位编码方向上的原始数据的格点位置和频率编码方向上的原始数据的格点位置正交；定义频率编码方向上的谱搜索的格点位置；对于每个频率编码方向上的谱搜索格点位置执行以下步骤：根据
相位编码方向上的谱搜索的格点位置、频率编码方向上的谱搜索的格点位置、相位编码方向上的原始数据的格点位置和频率编码方向上的原始数据的格点位置获得对应于相位编码和频率编码时的导向矢量，以获得对应于所述相位编码的阵列流形的一列。
[0008]在一些方案中，通过将对应于相位编码的阵列流形与目标子空间通过做内积以获得相关系数。
[0009]在一些方案中，可以设定信噪比阈值，如果信噪比低于信噪比阈值，则将对应于相位编码的阵列流形与目标子空间的逆矩阵做内积以获得相关系数，如果信噪比高于信噪比阈值，则将对应于相位编码的阵列流形与目标子空间做内积以获得相关系数。
[0010]在一些方案中，通过对相关系数取倒数以获得谱。
[0011]在一些方案中，通过对相关系数取倒数再乘以阵列流形的模以获得谱。
[0012]在一些方案中，通过对相关系数取倒数再乘以阵列流形的转置与阵列流形的乘积以获得谱。
[0013]在一些方案中，从谱中提取对角元素以获得重建出的图像的一行。
[0014]在一些方案中，图像重建方法还包括对重建出的图像进行归一化，以获得归一化图像。
[0015]在一些方案中，N比M小1。
[0016]在一些方案中，可以设定谱搜索的范围，以在特定的谱搜索范围内基于谱搜索范围确定相位编码方向上的原始数据的格点位置和频率编码方向上的原始数据的格点位置，以执行谱搜索步骤。
[0017]在一些方案中，图像重建方法是磁共振图像重建方法。
[0018]在一些方案中，原始数据是k空间数据。
[0019]在一些方案中，对应于相位编码和频率编码时的导向矢量Avec由下式定义：
[0020]Avec＝exp(2πi*k*(f1(Yx*Yer(vid))+f2(Yy*Hor(hid))))
[0021]其中，f1和f2分别表示一种变换，vid和hid分别表示第vid个相位编码方向上的谱搜索格点位置和第hid个频率编码方向上的谱搜索格点位置，k是波数，表示原子进动时的相位随距离变化的速率，其与共振频率的关系为其中c是光速。
[0022]在一些方案中，对应于相位编码和频率编码时的导向矢量Avec由下式定义：
[0023]Avec＝exp(2πi*k*(Yx*Ver(vid)+Yy*Hor(hid)))
[0024]根据本专利技术的第二方面，提供了一种计算机可读存储介质，其特征在于，其上存储有计算机可读的指令，当利用计算机执行所述指令时执行根据本专利技术的第一方面所述的图像重建方法。
附图说明
[0025]图1显示了根据本专利技术的实施例的图像重建方法的流程图；
[0026]图2显示了根据本专利技术的实施例的图像重建方法的程序代码示意图；
[0027]图3显示了根据本专利技术的实施例的图像重建方法获得的图像的分辨率与传统的傅里叶方法获得的图像的分辨率的对比；
[0028]图4显示了根据本专利技术的实施例的图像重建方法获得的图像的真实分辨能力与传
统的傅里叶方法获得的图像的真实分辨能力的可视化对比；
[0029]图5和图6显示了根据本专利技术的实施例的图像重建方法获得的图像的真实分辨能力与传统的傅里叶方法获得的图像的真实分辨能力的定量对比；
[0030]图7显示了根据本专利技术的实施例的图像重建方法的目标子空间与传统的傅里叶方法的k空间的对比；
[0031]图8显示了根据本专利技术的实施例的图像重建方法获得的重建信号与傅里叶变换方法获得的重建信号的对比。
具体实施方式
[0032]为了使得本专利技术的技术方案的目的、方案和优点更加清楚，下文中将结合本专利技术的具体实施例的附图，对本专利技术实施例的技术方案进行清楚、完整的描述。除非另有说明，否则本文所使用的术语具有本领域通常的含义。
[0033]图1显示了根据本专利技术的实施例的图像重建方法100的流程图，图2显示了根据本专利技术的实施例的图像重建方法100的程序代码示意图，需要注意的是，图2所示的程序代码仅是示例性的而非限制性的，可以采用任何其他合本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种图像重建方法，其特征在于，包括：读取原始数据(Y)；将所述原始数据(Y)降维，以获得一维向量(F)；对所述一维向量(F)取自相关，以获得自相关矩阵(R)；对所述自相关矩阵进行特征分解，以获得所述自相关矩阵(R)的特征值(Lam)和与相应特征值对应的特征向量(Vec)，所述特征值和特征向量的数量均为M；将所述特征值(Lam)按照大小进行排序，获得前N个较小的特征值对应的特征向量构成的子空间(SVN)，其中N小于M；对所述前N个较小的特征值对应的特征向量构成的子空间(SVN)取自相关，以获得目标子空间(SUN)；定义相位编码方向上的谱搜索的格点位置(Hor)；对于每个所述相位编码方向执行谱搜索步骤，所述谱搜索步骤包括：定义对应于所述相位编码方向的阵列流形(A)；获得对应于所述相位编码方向的阵列流形(A)与所述目标子空间(SUN)的相关系数，然后对所述相关系数进行负相关变换以获得谱(Spec)；从所述谱(Spec)中以预定规则提取元素以获得重建出的图像(X)的一行，所述预定规则使得所述谱(Spec)的每一行和每一列均有且仅有一个元素被提取；输出所述重建出的图像(X)。2.根据权利要求1所述的图像重建方法，其特征在于，定义对应于所述相位编码的阵列流形(A)的步骤包括：定义相位编码方向上的原始数据的格点位置(Yx)和频率编码方向上的原始数据的格点位置(Yy)，其中所述相位编码方向上的原始数据的格点位置(Yx)和所述频率编码方向上的原始数据的格点位置(Yy)正交；定义频率编码方向上的谱搜索的格点位置(Ver)；对于每个所述频率编码方向上的谱搜索格点位置(Ver)执行以下步骤：根据所述相位编码方向上的谱搜索的格点位置(Hor)、所述频率编码方向上的谱搜索的格点位置(Ver)、所述相位编码方向上的原始数据的格点位置(Yx)和所述频率编码方向上的原始数据的格点位置(Yy)获得对应于所述相位编码和所述频率编码时的导向矢量(Avec)，以获得对应于所述相位编码的阵列流形(A)的一列。3.根据权利要求1所述的图像重建方法，其特征在于，通过将对应于所述相位编码的阵列流形(A)与所述目标子空间(SUN)通过做内积以获得所述相关系数。4.根据权利要求3所述的图像重建方法，其特征在于，还包括设定信噪比阈值，如果信噪比低于所述信噪比阈值，则将对应于所述相位编码的阵列流形(A)与所述目标子空间(SUN)的逆矩阵做内积以...

【专利技术属性】
技术研发人员：孙东标，张紫豪，卓彦，王帆，
申请(专利权)人：中国科学院生物物理研究所，
类型：发明
国别省市：