一种mri超范围编码成像方法
【技术领域】
[0001]本发明涉及MRI图像处理技术领域，尤其是一种MRI超范围编码成像方法。
【背景技术】
[0002]磁共振成像(MRI)已发展成为一种非常强大的医学诊断手段，它特别适应于脑组织、神经系统以及人体软组织等部位的诊断检测，能够在清晰的解剖影像背景上明确显示出病变影像。MRI已经被广泛应用于医学临床诊断，成为重要的疾病诊断手段之一。中国发明专利申请CN 103300859A、CN 104161517A等现有技术公开了不同的核磁共振成像技术，但是，对于核磁共振对于视野小于被检部位时的图像伪影问题，不能彻底解决，仍然有明显的伪影。
【发明内容】
[0003]本发明要解决的技术问题是提供一种MRI超范围编码成像方法，能够解决现有技术的不足，降低了伪影的影响。
[0004]为解决上述技术问题，本发明所采取的技术方案如下。
[0005]一种MRI超范围编码成像方法，步骤如下:
[0006]A、获取视野范围内的第一图像，以及第一图像外侧的第二图像、第三图像…第η图像，其中n ^ 4，第二图像?第η图像均与第一图像部分重叠；
[0007]B、将第一图像作为中心图像进行编码成像，形成目标图像，其编码参数根据相关性分成若干组，每组参数的数量相同，其数量等于包含参数最多的参数组的参数数量，其余参数组空缺的参数使用本组参数的加权平均值代替，将所有参数组合成第一编码矩阵，以目标图像为目标值，按照上述编码方法对剩余图像进行依次编码，得到第二编码矩阵、第三编码矩阵…第η编码矩阵，在形成编码矩阵时保证每个编码矩阵均为非奇异矩阵；
[0008]C、依次选择第二图像?第η图像为中心图像，重复步骤B的过程，每个中心图像均得到其对应的n-Ι个编码矩阵；
[0009]D、将与第一图像相关的编码矩阵归为甲类编码矩阵，其余编码矩阵归为乙类编码矩阵，将甲类编码矩阵内的编码矩阵分为若干个矩阵对，每个矩阵对内的两个图像互为中心图像和目标图像，将矩阵对内未以第一图像作为中心图像的编码矩阵求逆，然后与另一个以第一图像作为中心图像的编码矩阵进行融合，形成若干个融合矩阵，逐一求取乙类编码矩阵的特征向量，将特征向量与其相关的融合矩阵进行二次融合；
[0010]其中特征向量与融合矩阵相关的定义为:求取特征向量的乙类编码矩阵所涉及的两个图像中任意一个图像与融合矩阵所涉及的非第一图像相同，即定义特征向量与融合矩阵相关；
[0011]Ε、使用步骤D中得到的二次融合矩阵对第二图像?第η图像进行重新编码，编码的结果与以第一图像作为中心图像进行编码形成的目标图像进行比较，若比较差异小于等于阈值，则将η个图像进行合并，得到最终编码结果，若比较差异大于阈值，则以最终编码结果作为新的第一图像，重复步骤B?步骤E，直至比较差异小于阈值。
[0012]作为本发明的一种优选技术方案，步骤D中，两个编码矩阵的融合采用K-means聚类算法融合，特征向量与其相关的融合矩阵的融合采用特征向量与融合矩阵相乘的方式进行。
[0013]作为本发明的一种优选技术方案，步骤E中，编码后形成的图像进行比较的方法为，选取相等数量的参考像素点，两个图像中像素点灰度值相减，其结果中的最大值与两个图像像素点灰度值的平均值之比即为差异值。
[0014]为本发明的一种优选技术方案，所述阈值为5%?r10A。
[0015]为本发明的一种优选技术方案，选取参考像素点时，参考像素点在图像中不同位置的分布密度与步骤A中不同图像间的重叠次数成反比。
[0016]采用上述技术方案所带来的有益效果在于:本发明通过不同图像间的互相编码处理，得出多组相互关联的编码参数矩阵，通过对编码矩阵进行融合处理，得到较佳编码参数，可以逐渐消除编码过程中不同图像带来的伪影。
【附图说明】
[0017]图1是未采用本发明成像方法得到的腰椎矢状位图像。
[0018]图2是使用了本发明提供的成像方法得到的腰椎矢状位图像。
【具体实施方式】
[0019]由于图像选取的数量与循环计算的循环次数成反比，以腰椎处的核磁共振成像为例，这个视野范围选取6张图像进行编码处理可以以最少的计算量得到合格的图像。
[0020]一种MRI超范围编码成像方法，步骤如下:
[0021]A、获取视野范围内的第一图像，以及第一图像外侧的第二图像、第三图像…第六图像，第二图像?第六图像均与第一图像部分重叠；
[0022]B、将第一图像作为中心图像进行编码成像，形成目标图像，其编码参数根据相关性分成若干组，每组参数的数量相同，其数量等于包含参数最多的参数组的参数数量，其余参数组空缺的参数使用本组参数的加权平均值代替，将所有参数组合成第一编码矩阵，以目标图像为目标值，按照上述编码方法对剩余图像进行依次编码，得到第二编码矩阵、第三编码矩阵…第六编码矩阵，在形成编码矩阵时保证每个编码矩阵均为非奇异矩阵；在计算本组参数的加权平均值时，每个参数的加权量与此参数对图像进行编码时所进行的映射路径的长度成反比；
[0023]C、依次选择第二图像?第六图像为中心图像，重复步骤B的过程，每个中心图像均得到其对应的五个编码矩阵；
[0024]D、将与第一图像相关的编码矩阵归为甲类编码矩阵，其余编码矩阵归为乙类编码矩阵，将甲类编码矩阵内的编码矩阵分为若干个矩阵对，每个矩阵对内的两个图像互为中心图像和目标图像，将矩阵对内未以第一图像作为中心图像的编码矩阵求逆，然后与另一个以第一图像作为中心图像的编码矩阵进行融合，形成若干个融合矩阵，逐一求取乙类编码矩阵的特征向量，将特征向量与其相关的融合矩阵进行二次融合；
[0025]其中特征向量与融合矩阵相关的定义为:求取特征向量的乙类编码矩阵所涉及的两个图像中任意一个图像与融合矩阵所涉及的非第一图像相同，即定义特征向量与融合矩阵相关；
[0026]两个编码矩阵的融合采用K-means聚类算法融合，若融合后的编码矩阵与任意一个融合前的编码矩阵的欧氏距离大于融合后的编码矩阵与另一个融合前的编码矩阵的欧氏距离的2倍，则使用融合后的编码矩阵再次与和其欧氏距离较近的编码矩阵进行融合，直至融合后的编码矩阵与任意一个融合前的编码矩阵的欧氏距离小于等于融合后的编码矩阵与另一个融合前的编码矩阵的欧氏距离的2倍；特征向量与其相关的融合矩阵的融合采用特征向量与融合矩阵相乘的方式进行，融合后再乘以一个融合因子，融合因子的大小与特征向量所对应的特征值成正比；
[0027]E、使用步骤D中得到的二次融合矩阵对第二图像?第六图像进行重新编码，编码的结果与以第一图像作为中心图像进行编码形成的目标图像进行比较，若比较差异小于等于阈值，则将六个图像进行合并，得到最终编码结果，若比较差异大于阈值，则以最终编码结果作为新的第一图像，重复步骤B?步骤E，直至比较差异小于阈值，编码后形成的图像进行比较的方法为，选取相等数量的参考像素点，两个图像中像素点灰度值相减，其结果中的最大值与两个图像像素点灰度值的平均值之比即为差异值，阈值为6.5%，选取参考像素点时，参考像素点在图像中不同位置的分布密度与步骤A中不同图像间的重叠次数成反比。
[0028]参考图1和图2，可以明显看出，本发明所提供的成像方法可以明显减低伪影的影响，提高图像清晰度。
[0029]上述描述仅作为本发明可实施的技术方案提出，不作为对其技术方案本身的单一限制条件。
【主权项】
1.一种MRI超范围编码成像方法，其特征在于步骤如下: A、获取视野范围内的第一图像，以及第一图像外侧的第二图像、第三图像…第η图像，其中n ^ 4，第二图像?第η图像均与第一图像部分重叠； B、将第一图像作为中心图像进行编码成像，形成目标图像，其编码参数根据相关性分成若干组，每组参数的数量相同，其数量等于包含参数最多的参数组的参数数量，其余参数组空缺的参数使用本组参数的加权平均值代替，将所有参数组合成第一编码矩阵，以目标图像为目标值，按照上述编码方法对剩余图像进行依次编码，得到第二编码矩阵、第三编码矩阵…第η编码矩阵，在形成编码矩阵时保证每个编码矩阵均为非奇异矩阵； C、依次选择第二图像?第η图像为中心图像，重复步骤B的过程，每个中心图像均得到其对应的n-Ι个编码矩阵； D、将与第一图像相关的编码矩阵归为甲类编码矩阵，其余编码矩阵归为乙类编码矩阵，将甲类编码矩阵内的编码矩阵分为若干个矩阵对，每个矩阵对内的两个图像互为中心图像和目标图像，将矩阵对内未以第一图像作为中心图像的编码矩阵求逆，然后与另一个以第一图像作为中心图像的编码矩阵进行融合，形成若干个融合矩阵，逐一求取乙类编码矩阵的特征向量，将特征向量与其相关的融合矩阵进行二次融合； 其中特征向量与融合矩阵相关的定义为:求取特征向量的乙类编码矩阵所涉及的两个图像中任意一个图像与融合矩阵所涉及的非第一图像相同，即定义特征向量与融合矩阵相关; Ε、使用步骤D中得到的二次融合矩阵对第二图像?第η图像进行重新编码，编码的结果与以第一图像作为中心图像进行编码形成的目标图像进行比较，若比较差异小于等于阈值，则将η个图像进行合并，得到最终编码结果，若比较差异大于阈值，则以最终编码结果作为新的第一图像，重复步骤B?步骤Ε，直至比较差异小于阈值。
2.根据权利要求1所述的MRI超范围编码成像方法，其特征在于:步骤D中，两个编码矩阵的融合采用K-means聚类算法融合，特征向量与其相关的融合矩阵的融合采用特征向量与融合矩阵相乘的方式进行。
3.根据权利要求1所述的MRI超范围编码成像方法，其特征在于:步骤E中，编码后形成的图像进行比较的方法为，选取相等数量的参考像素点，两个图像中像素点灰度值相减，其结果中的最大值与两个图像像素点灰度值的平均值之比即为差异值。
4.根据权利要求3所述的MRI超范围编码成像方法，其特征在于:所述阈值为5%?7%。
5.根据权利要求3所述的MRI超范围编码成像方法，其特征在于:选取参考像素点时，参考像素点在图像中不同位置的分布密度与步骤A中不同图像间的重叠次数成反比。
【专利摘要】本发明公开了一种MRI超范围编码成像方法，包括以下步骤：A、获取若干个图像；B、以第一图像作为中心图像进行编码成像，形成目标图像，并依此目标图像为目标值，对剩余图像进行依次编码；C、依次选取其它图像为中心图像，进行上述步骤B的编码过程；D、对编码矩阵进行融合；E、进行循环编码比较，最终得到最终图像。本发明能够解决现有技术的不足，降低了伪影的影响。
【IPC分类】G06T5-50, A61B5-055, G06T9-00, G06T5-00
【公开号】CN104766351
【申请号】CN201510197277
【发明人】陈亮, 王天峰
【申请人】惠仁望都医疗设备科技有限公司
【公开日】2015年7月8日
【申请日】2015年4月24日