基于PCA的磁共振成像参数虚拟仿真方法技术

本发明专利技术公开了一种基于PCA的磁共振成像参数虚拟仿真方法，包括以下步骤：S1:确定所需的磁共振成像参数种类,且计算所有参数种类的全组合；S2:对全组合中的每一条参数组合，在真实磁共振扫描仪上按参数设定进行扫描成像或进行有限元计算仿真成像，获得成像的三维体数据的数据库；S3:对数据库利用PCA算法进行降维处理，得到一维压缩向量；S4:根据成像参数设定找到临近参数组合，查询其对应压缩向量，进行线性插值，将其重新折叠为三维体数据，获得虚拟仿真成像。本发明专利技术使用PCA算法对磁共振多参数组合预成像数据库进行降维，降维后只需存储低维的压缩向量与少量主成分图像，在仿真过程中对压缩向量进行插值，随后进行图像恢复，具有低存储量和低计算量。低存储量和低计算量。低存储量和低计算量。

【技术实现步骤摘要】
基于PCA的磁共振成像参数虚拟仿真方法


[0001]本专利技术涉及磁共振成像领域，尤其涉及一种基于PCA的磁共振成像参数虚拟仿真方法。

技术介绍

[0002]磁共振成像系统作为高技术含量和高附加值的大型医学成像设备，具有原理复杂、价格高昂的特点，并且在医院临床应用中机时紧张。因此，在高校的课程教学和医院的临床培训中存在着进行磁共振成像参数虚拟仿真的需要。其具体的内容为，对于一个虚拟病例或体模，使用者可以调整与成像相关的各种参数，仿真系统根据参数的设置给出与真实磁共振扫描在相同参数下结果相近的仿真成像结果。在此过程中，使用者并不关心真实成像过程的物理原理，而是为了熟悉不同的参数设置对成像结果有什么样的影响。在以往的实践中，磁共振成像虚拟仿真主要包括以下2种实现方法。
[0003]一种方法是基于数字化体模与理想化物理模型的实时计算仿真。该方法对磁共振成像过程的物理原理进行了理想化简，再根据使用者给出的成像参数设置以及数字化体模中与磁共振成像有关的物理量的数值分布，通过有限元方法按一定的时间步长与空间分辨率进行成像过程的计算推演。待推演结束后，即获得了模拟的磁共振原始信号。随后按与真实成像过程相同的图像重建算法处理原始信号后就能得到仿真的成像结果。这种仿真方法只需要存储数字化体模与计算程序，占用的存储空间较小，但是其计算步骤量非常大，有很多三角与指数运算，并且出于精度考虑需要使用双精度浮点数据类型。为了满足使用者对实时性的需要，通常使用高性能专业计算加速卡(一般GPU由于双精度计算能力不足无法满足这一方法)或者数十核的计算服务器来进行计算，硬件成本高昂。
[0004]另一种方法是基于多参数组合预成像数据库的插值仿真方法。该方法首先根据每一个成像参数的可调整范围取有限个数的离散取值，再将所有参数的离散取值生成全组合。对全组合中的每一条参数组合，进行实际的成像或前一种方法所述的计算仿真，获得成像结果并保存，构建一个预成像的数据库。实际使用时，根据使用者给出的成像参数确定其在全组合中的邻近参数组合。从数据库中取出邻近参数组合对应的预成像结果进行线性插值，得到最终仿真结果。这一方法不需要进行成像过程的实时计算仿真，不需要昂贵的计算平台，但是预成像数据库的容量非常大。首先，全组合的数量非常大。其次磁共振的成像结果为3维高精度体数据，单个图像的容量就可达数十兆。例如，在满足基本需求的6个参数可调，每个参数在数据库中有5个离散取值的情况下，整个数据库的总容量将达到700GB以上，这对数据发布和传输产生了很大压力。另一方面，虽然这一方法的计算量相对于实时仿真有所降低，但对许多邻近参数组合对应的成像结果的三维体数据进行插值仍然具有可观的运算量。

技术实现思路

[0005]本专利技术提出使用PCA算法对磁共振多参数组合预成像数据库进行降维，降维后只
需存储低维的压缩向量与少量主成分图像，具有低存储量；提出在仿真过程中只需对压缩向量进行插值，随后进行图像恢复，具有低计算量，可适用于在WEB平台或移动设备平台等存储和计算受限的情境下进行磁共振成像参数虚拟仿真。
[0006]为实现上述目的，本专利技术提供一种基于PCA的磁共振成像参数虚拟仿真方法，包括以下步骤：
[0007]S1:确定所需的磁共振成像参数种类,且计算所有参数种类的全组合；
[0008]S2:对全组合中的每一条参数组合，在真实磁共振扫描仪上按参数设定进行扫描成像或进行有限元计算仿真成像，获得成像的三维体数据的数据库；
[0009]S3:对数据库利用PCA算法进行降维处理，得到一维压缩向量；
[0010]S4:根据成像参数设定找到临近参数组合，查询其对应压缩向量，进行线性插值，将其重新折叠为三维体数据，获得虚拟仿真成像。
[0011]作为优选，在步骤S1中，根据所需的磁共振成像参数的种类，在每个成像参数可调整范围中确定有限个数的离散数值。
[0012]作为优选，在步骤S2中，对全组合中的每一条参数组合，获得成像结果后进行保存，构建一个预成像的数据库，这个数据库中保存磁共振成像参数组合和对应的成像三维体数据。
[0013]作为优选，对数据库利用PCA算法进行降维时，首先将每一个三维体数据展开为一维压缩向量，则数据库中的原本所有三维体数据构成了第一矩阵，计算第一矩阵的行平均向量，最终得到归一化的第二矩阵。
[0014]作为优选，获得第二矩阵后代入公式计算得到协方差矩阵，得到协方差矩阵的特征值和特征向量，取对应特征值最大的前若干个特征向量按照排列组成变换矩阵，计算降维后的结果矩阵和主成分矩阵。
[0015]作为优选，在步骤S4中，根据使用者给出的成像参数设定找到若干个参数向上和向下最接近的两个离散取值，进而确定全组合中与使用者参数设定最接近的邻近参数组合以及对应的压缩向量，进行线性插值，然后按照相关公式进行重建，将其重新折叠为三维体数据，即得到虚拟仿真成像的结果；
[0016]本专利技术的有益效果是：与现有技术相比，本专利技术提供的基于PCA的磁共振成像参数虚拟仿真方法，通过在预成像数据库中查找邻近参数的图像进行插值来得到仿真结果，不需要进行成像过程的实时计算仿真，不需要昂贵的计算平台，使用PCA算法对预成像数据进行降维，把大量高容量的三维体数据转化为低维的压缩向量和少量主成分图像；仿真成像过程中对压缩向量进行插值，再进行数据重建，而不是对三维体数据进行插值，进一步降低了计算量。
附图说明
[0017]图1为本专利技术的步骤流程图。
具体实施方式
[0018]为了更清楚地表述本专利技术，下面结合附图和实施例对本专利技术作进一步地描述，本
技术人员可以理解，除非特意声明，这里使用的单数形式“一”、“一个”、“所述”和
“
该”也可包括复数形式。应该进一步理解的是，本申请的说明书中使用的措辞“包括”是指存在所述特征、整数、步骤、操作、元件和/或组件，但是并不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、组件和/或它们的组。应该理解，这里使用的措辞“和/或”包括一个或更多个相关联的列出项的全部或任一单元和全部组合。
[0019]本
技术人员可以理解，除非另外定义，这里使用的所有术语(包括技术术语和科学术语)，具有与本申请所属领域中的普通技术人员的一般理解相同的意义。还应该理解的是，诸如通用字典中定义的那些术语，应该被理解为具有与现有技术的上下文中的意义一致的意义，并且除非像这里一样被特定定义，否则不会用理想化或过于正式的含义来解释。
[0020]请参阅图1，本专利技术公开了一种基于PCA的磁共振成像参数虚拟仿真方法，包括以下步骤：S1:确定所需的磁共振成像参数种类,且计算所有参数种类的全组合；S2:对全组合中的每一条参数组合，在真实磁共振扫描仪上按参数设定进行扫描成像，获得具有成像的三维体数据的数据库；S3:对数据库利用PCA算法进行降维处理，得到一维压缩向量；S4:根据成像参数设定找到临近参数组合，查询其对应压缩向量，进行线性插值，将其重新折叠为三维本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于PCA的磁共振成像参数虚拟仿真方法，其特征在于，包括以下步骤：S1:确定所需的磁共振成像参数种类,且计算所有参数种类的全组合；S2:对全组合中的每一条参数组合，在真实磁共振扫描仪上按参数设定进行扫描成像或进行有限元计算仿真成像，获得成像的三维体数据的数据库；S3:对数据库利用PCA算法进行降维处理，得到一维压缩向量；S4:根据成像参数设定找到临近参数组合，查询其对应压缩向量，进行线性插值，将其重新折叠为三维体数据，获得虚拟仿真成像。2.根据权利要求1所述的基于PCA的磁共振成像参数虚拟仿真方法，其特征在于，在步骤S1中，根据所需的磁共振成像参数的种类，在每个成像参数可调整范围中确定有限个数的离散数值。3.根据权利要求1所述的基于PCA的磁共振成像参数虚拟仿真方法，其特征在于，在步骤S2中，对全组合中的每一条参数组合，获得成像结果后进行保存，构建一个预成像的数据库，这个数据库中保存磁共振成像参数组合和对应的成像三维体...

【专利技术属性】
技术研发人员：孙钰，梁嘉炜，沈靖菲，符谦益，杨明，
申请(专利权)人：南京伯睿生命科学研究院有限公司，
类型：发明
国别省市：