CN114264997B - 梯度灵敏度校准方法、装置及磁共振设备 - Google Patents
梯度灵敏度校准方法、装置及磁共振设备 Download PDFInfo
- Publication number
- CN114264997B CN114264997B CN202111525705.1A CN202111525705A CN114264997B CN 114264997 B CN114264997 B CN 114264997B CN 202111525705 A CN202111525705 A CN 202111525705A CN 114264997 B CN114264997 B CN 114264997B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- ellipsoid
- fitting
- target axis
- resonance system
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 230000035945 sensitivity Effects 0.000 title claims abstract description 119
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000004590 computer program Methods 0.000 claims description 27
- 238000003384 imaging method Methods 0.000 claims description 9
- 230000000694 effects Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 239000012528 membrane Substances 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R35/00—Testing or calibrating of apparatus covered by the other groups of this subclass
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Signal Processing (AREA)
- High Energy & Nuclear Physics (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请涉及一种梯度灵敏度校准方法、装置及磁共振设备。所述方法包括：获取球形模体的三维图像；基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸；根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对所述磁共振系统的梯度灵敏度进行校准。采用本方法能够确保对磁共振系统的梯度灵敏度的校准准确度。
Description
技术领域
本申请涉及磁共振技术领域，特别是涉及一种梯度灵敏度校准方法、装置及磁共振设备。
背景技术
在磁共振系统中，为了确保采集的磁共振图像的尺寸与实物尺寸保持一致，需要对磁共振系统的梯度灵敏度进行校准。
传统技术中，大多使用球形膜体，利用球形模体中心点上冠状面、矢状面、横断面三个层面尺寸的一致性对磁共振系统的梯度灵敏度进行校准。
然而，传统的磁共振系统的梯度灵敏度的校准方法，存在校准准确度较低的问题。
发明内容
基于此，有必要针对上述技术问题，提供一种能够提高磁共振系统的梯度灵敏度的校准准确度的梯度灵敏度校准方法、装置及磁共振设备。
第一方面，本申请提供了一种梯度灵敏度校准方法，所述方法包括：
获取球形模体的三维图像；
基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸；
根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对所述磁共振系统的梯度灵敏度进行校准。
在其中一个实施例中，所述根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准，包括：
获取所述拟合尺寸与所述球形模体的目标轴的实际尺寸的偏差值；
根据所述偏差值，对所述磁共振系统的梯度灵敏度因子进行调整，以对所述磁共振系统的梯度灵敏度进行校准。
在其中一个实施例中，所述根据所述偏差值，对所述磁共振系统的梯度灵敏度因子进行调整，以对所述磁共振系统的梯度灵敏度进行校准，包括：
确定所述偏差值是否满足预设的偏差范围；
若否，则根据所述偏差值对所述磁共振系统的梯度灵敏度因子进行调整，以对所述磁共振系统的梯度灵敏度进行校准
在其中一个实施例中，所述目标轴包括以磁体成像区域的中心位置为原点，以空间正交的三轴分别作为三轴方向建立的坐标系所对应的三个轴。
在其中一个实施例中，所述获取球形模体的三维图像，包括：
使用多片层扫描协议对所述球形模体进行扫描，得到所述球形模体的三维图像。
在其中一个实施例中，所述方法应用于悬臂床磁共振系统。
第二方面，本申请还提供了一种梯度灵敏度校准装置，所述装置包括：
第一获取模块，用于获取球形模体的三维图像；
第二获取模块，用于基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸；
校准模块，用于根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对所述磁共振系统的梯度灵敏度进行校准。
第三方面，本申请还提供了一种磁共振设备，所述计算机设备包括存储器和处理器，所述存储器存储有计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
获取球形模体的三维图像；
基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸；
根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对所述磁共振系统的梯度灵敏度进行校准。
第四方面，本申请还提供了一种计算机可读存储介质，所述计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
获取球形模体的三维图像；
基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸；
根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对所述磁共振系统的梯度灵敏度进行校准。
第五方面，本申请还提供了一种计算机程序产品，所述计算机程序产品，包括计算机程序，该计算机程序被处理器执行时实现以下步骤：
获取球形模体的三维图像；
基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸；
根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对所述磁共振系统的梯度灵敏度进行校准。
上述梯度灵敏度校准方法、装置及磁共振设备，通过获取球形模体的三维图像，基于球形模体的三维图像进行椭球拟合，能够得到椭球的目标轴的拟合尺寸，从而可以根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准，该过程降低了校准时对球模摆位的要求，通过拟合的椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸就能够对磁共振系统的梯度灵敏度进行校准，确保了对磁共振系统的梯度灵敏度的校准准确度。
附图说明
图1为一个实施例中梯度灵敏度校准方法的应用环境图；
图2为一个实施例中梯度灵敏度校准方法的流程示意图；
图3为另一个实施例中梯度灵敏度校准方法的流程示意图；
图4为另一个实施例中梯度灵敏度校准方法的流程示意图
图5为一个实施例中梯度灵敏度校准装置的结构框图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
本申请实施例提供的梯度灵敏度校准方法，可以适用于如图1所示的磁共振设备。该磁共振设备包括通过系统总线连接的处理器、存储器，该存储器中存储有计算机程序，处理器执行该计算机程序时可以执行下述方法实施例的步骤。可选的，该磁共振设备还可以包括网络接口、显示屏和输入装置。其中，该磁共振设备的处理器用于提供计算和控制能力。该磁共振设备的存储器包括非易失性存储介质、内存储器，该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该磁共振设备的网络接口用于与外部的终端通过网络连接通信。
在一个实施例中，如图2所示，提供了一种梯度灵敏度校准方法，以该方法应用于图1中的计算机设备为例进行说明，包括以下步骤：
S201，获取球形模体的三维图像。
其中，球形模体为用于对磁共振系统进行检测的模体，该球形模型能够满足对磁共振的轴向面、冠状面、矢状面图像性能的检测，磁共振的检测功能主要包括空间均匀性、扫描层厚/层间距、准直定位系统验证、空间分辨力、几何畸变率(空间线性)、信噪比(SNR)、低对比灵敏度(低对比分辨力)、T1和T2驰豫时间值等测量内容。
具体地，球形模体的三维图像能够用来表征球形模体的三维结构特征，磁共振设备可以通过对球形模体进行三维扫描，得到球形模体的三维图像。或者，磁共振设备也可以从存储图像的数据库中获取预先对球形模体进行扫描得到的球形模体的三维图像，本实施中对获取球形模体的三维图像的具体时间和步骤不加以限制，只要能够通过对球形模体进行磁共振扫描成像得到球形模体的三维图像即可。
S202，基于三维图像进行椭球拟合，得到椭球的目标轴的拟合尺寸。
可选的，磁共振设备可以利用预设的拟合算法对球形模体的三维图像进行椭球拟合，基于球形模体的三维图像拟合出一个椭球，或者，磁共振设备也可以对球形模体的三维图像进行重建调整，拟合出一个椭球，从而磁共振设备基于拟合的椭球，得到椭球的目标轴的拟合尺寸。可选的，本实施例中的目标轴包括以磁体成像区域的中心位置为原点，以空间正交的三轴分别作为三轴方向建立的坐标系所对应的三个轴，也就是磁共振设备得到的椭球的目标轴的拟合尺寸为椭球的空间正交的三轴的拟合尺寸。可以理解的是，以磁体成像区域的中心位置为原点，以空间正交的三轴分别作为三轴方向建立的坐标系中每点的坐标是可以获取的，因此，磁共振设备可以通过该坐标系中每点的坐标，得到椭球的最长轴的拟合尺寸，再基于得到的椭球的最长轴的拟合尺寸和坐标系中各点的坐标，得到椭球的上述目标轴的拟合尺寸。
S203，根据拟合尺寸和球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准。
可选的，磁共振设备可以获取椭球的目标轴的拟合尺寸和上述球形模体的目标轴的实际尺寸的相似度，根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的相似度，对磁共振系统的梯度灵敏度进行校准。例如，若椭球的目标轴的拟合尺寸和上述球形模体的目标轴的实际尺寸的相似度值较大，则说明椭球的目标轴的拟合尺寸和上述球形模体的目标轴的实际尺寸较为相似，磁共振设备可以对磁共振系统的梯度灵敏度进行细微地调整或者保持当前的磁共振系统的梯度灵敏度不变；若椭球的目标轴的拟合尺寸和上述球形模体的目标轴的实际尺寸的相似度值较小，则说明椭球的目标轴的拟合尺寸和上述球形模体的目标轴的实际尺寸相差较大，磁共振设备可以对磁共振系统的梯度灵敏度进行大幅度调整，直至椭球的目标轴的拟合尺寸和上述球形模体的目标轴的实际尺寸相近。可以理解的是，在理想情况下，椭球的目标轴的拟合尺寸为球形模体的目标轴的实际尺寸，也就是说在理想情况下，椭球的目标轴的拟合尺寸等于球形模体的目标轴的实际尺寸。
在一些实施例中，本申请实施例提供的梯度灵敏度校准方法应用于悬臂床磁共振系统。例如超高场动物磁共振成像系统。动物磁共振成像系统的扫描腔内径比较小而且长度较长，一般采用固定在扫描腔外的可伸缩的悬臂床将待扫描动物送入扫描腔的磁体成像区域的中心位置。相比于一般的临床采用的滑轨床磁共振系统，通常，动物磁共振成像系统采用的悬臂床磁共振系统的悬臂末端下垂造成模体中心位置不易确定，采用本申请的实施例的梯度灵敏度校准方法，降低了校准时对球模摆位的要求，避免了悬臂床磁共振系统的悬臂床摆动的影响，操作简单，准确度高。
上述梯度灵敏度校准方法中，通过获取球形模体的三维图像，基于球形模体的三维图像进行椭球拟合，能够得到椭球的目标轴的拟合尺寸，从而可以根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准，该过程降低了校准时对球模摆位的要求，通过拟合的椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸就能够对磁共振系统的梯度灵敏度进行校准，确保了对磁共振系统的梯度灵敏度的校准准确度。
在上述根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准的场景中，在一个实施例中，如图3所示，上述S203，包括：
S301，获取拟合尺寸与球形模体的目标轴的实际尺寸的偏差值。
可选的，磁共振设备可以根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的差值，得到椭球的目标轴的拟合尺寸与球形模体的目标轴的实际尺寸的偏差值，与可以根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的比值，得到椭球的目标轴的拟合尺寸与球形模体的目标轴的实际尺寸的偏差值。可选的，若椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的差值越大，则说明椭球的目标轴的拟合尺寸与球形模体的目标轴的实际尺寸的偏差值越大；若椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的差值越小，则说明椭球的目标轴的拟合尺寸与球形模体的目标轴的实际尺寸的偏差值越小。可选的，若椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的比值越靠近1，则说明椭球的目标轴的拟合尺寸与球形模体的目标轴的实际尺寸的偏差值越小；若椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的比值越远离1，则说明椭球的目标轴的拟合尺寸与球形模体的目标轴的实际尺寸的偏差值越大。
S302，根据偏差值，对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
可选的，磁共振设备可以根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值的大小对磁共振系统的梯度灵敏度因子进行调整，例如，若对磁共振系统的梯度灵敏度因子进行调整后，椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值减小，则磁共振设备可以按照当前的调整方向继续对磁共振系统的梯度灵敏度因子进行调整，否则，磁共振设备可以调整当前的调整方向，以与之前调整方向相反的调整方向对磁共振系统的梯度灵敏度因子进行调整。例如，若获取的椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值为0.8，对磁共振系统的梯度灵敏度因子进行调整后，椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值为0.5，则磁共振设备可以继续以当前的调整方向继续对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。又例如，若获取的椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值为0.8，对磁共振系统的梯度灵敏度因子进行调整后，椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值为0.9，则磁共振设备可以调整当前的调整方向，以与之前调整方向相反的调整方向对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
本实施例中，磁共振设备通过获取椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值，能够根据该偏差值，准确地对磁共振系统的梯度灵敏度因子进行调整，由于对磁共振系统的梯度灵敏度因子调整的准确度得到了提高从而也能够准确地对磁共振系统的梯度灵敏度进行准确地校准，提高了对磁共振系统的梯度灵敏度校准的准确度。
在上述根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值，对磁共振系统的梯度灵敏度因子进行调整的场景中，可以根据该偏差值和预设的偏差范围，对磁共振系统的梯度灵敏度因子进行调整，在上述实施例的基础上，在一个实施例中，如图4所示，上述S302，包括：
S401，确定偏差值是否满足预设的偏差范围。
其中，预设的偏差范围可以表征磁共振系统的梯度灵敏度因子已调整至正常的范围值内。可选的，磁共振设备可以将上述偏差值和预设的偏差范围内的值进行比较，确定上述偏差值是否在预设的偏差范围内。可以理解的是，若该偏差值超出了上述预设的偏差范围内的值，则表示该偏差值不在上述预设的偏差范围内；若该偏差值在上述预设的偏差范围内，则表示该偏差值在上述预设的偏差范围内。
S402，若偏差值不满足预设的偏差范围，则根据偏差值对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
具体地，若磁共振设备确定上述偏差值不在预设的偏差范围内，则磁共振设备可以根据该偏差值对磁共振系统的梯度灵敏度因子进行调整，以使偏差值落入上述预设的偏差范围内。可选的，磁共振设备可以根据上述偏差值与预设的偏差范围内的值的差值，对磁共振系统的梯度灵敏度因子进行调整，例如，若上述偏差值与上述预设的偏差范围内的差值较大，则磁共振设备可以以当前磁共振系统的梯度灵敏度因子为参考，向该梯度灵敏度因子的相反方向进行调整，以对磁共振系统的梯度灵敏度进行校准。
本实施例中，磁共振设备可以准确地确定出椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值是否在预设的偏差范围内，进而可以基于准确度较高的确定结果，根据椭球的目标轴的拟合尺寸和球形模体的目标轴的实际尺寸的偏差值对磁共振系统的梯度灵敏度因子进行准确地调整，以提高对磁共振系统的梯度灵敏度的校准准确度。
在上述获取球形模体的三维图像的场景中，在上述实施例的基础上，在一个实施例中，上述S201，包括：使用多片层扫描协议对球形模体进行扫描，得到球形模体的三维图像。
具体地，多片层扫描协议能够使磁共振设备对球形模体进行多层扫描，从多个层次得到球形模体的三维图像，确保了得到的球形模体的三维图像的准确度。需要说明的是，本实施例中的球形模体的三维图像为重建后的图像，图像中表征的是球形模体的三维立体信息。
本实施例中，磁共振设备使用多片层扫描协议能够对球形模体进行准确地扫描，这样得到的球形模体的三维图像是从多个片层进行扫描得到的，从而提高了得到的球形模体的三维图像。
应该理解的是，虽然如上的各实施例所涉及的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，如上的各实施例所涉及的流程图中的至少一部分步骤可以包括多个步骤或者多个阶段，这些步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤中的步骤或者阶段的至少一部分轮流或者交替地执行。
基于同样的发明构思，本申请实施例还提供了一种用于实现上述所涉及的梯度灵敏度校准方法的梯度灵敏度校准装置。该装置所提供的解决问题的实现方案与上述方法中所记载的实现方案相似，故下面所提供的一个或多个梯度灵敏度校准装置实施例中的具体限定可以参见上文中对于梯度灵敏度校准方法的限定，在此不再赘述。
在一个实施例中，如图5所示，提供了一种梯度灵敏度校准装置，包括：第一获取模块、第二获取模块和校准模块，其中：
第一获取模块，用于获取球形模体的三维图像；
第二获取模块，用于基于三维图像进行椭球拟合，得到椭球的目标轴的拟合尺寸；
校准模块，用于根据拟合尺寸和球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准。
可选的，上述目标轴包括以磁体成像区域的中心位置为原点，以空间正交的三轴分别作为三轴方向建立的坐标系所对应的三个轴。
可选的，本实施例提供的梯度灵敏度校准装置应用于悬臂床磁共振系统。
本实施例提供的梯度灵敏度校准装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
在上述实施例的基础上，可选的，上述校准模块，包括：第一获取单元和校准单元，其中：
第一获取单元，用于获取拟合尺寸与球形模体的目标轴的实际尺寸的偏差值。
校准单元，用于根据偏差值，对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
本实施例提供的梯度灵敏度校准装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
在上述实施例的基础上，可选的，上述校准单元，用于确定偏差值是否满足预设的偏差范围；若否，则根据偏差值对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
本实施例提供的梯度灵敏度校准装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
在上述实施例的基础上，可选的，上述第一获取模块，包括：第二获取单元，其中：
第二获取单元，用于使用多片层扫描协议对球形模体进行扫描，得到球形模体的三维图像。
本实施例提供的梯度灵敏度校准装置，可以执行上述方法实施例，其实现原理和技术效果类似，在此不再赘述。
上述梯度灵敏度校准装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
本领域技术人员可以理解，图1中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种磁共振设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
获取球形模体的三维图像；
基于三维图像进行椭球拟合，得到椭球的目标轴的拟合尺寸；
根据拟合尺寸和球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准。
可选的，目标轴包括以磁体成像区域的中心位置为原点，以空间正交的三轴分别作为三轴方向建立的坐标系所对应的三个轴。
可选的，该处理器执行的计算机程序应用于悬臂床磁共振系统。
上述实施例提供的磁共振设备，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
获取拟合尺寸与球形模体的目标轴的实际尺寸的偏差值；
根据偏差值，对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
上述实施例提供的磁共振设备，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
确定偏差值是否满足预设的偏差范围；
若否，则根据偏差值对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
上述实施例提供的磁共振设备，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
使用多片层扫描协议对球形模体进行扫描，得到球形模体的三维图像。
上述实施例提供的磁共振设备，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取球形模体的三维图像；
基于三维图像进行椭球拟合，得到椭球的目标轴的拟合尺寸；
根据拟合尺寸和球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准。
可选的，目标轴包括以磁体成像区域的中心位置为原点，以空间正交的三轴分别作为三轴方向建立的坐标系所对应的三个轴。
可选的，该处理器执行的计算机程序应用于悬臂床磁共振系统。
上述实施例提供的计算机可读存储介质，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
获取拟合尺寸与球形模体的目标轴的实际尺寸的偏差值；
根据偏差值，对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
上述实施例提供的计算机可读存储介质，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
确定偏差值是否满足预设的偏差范围；
若否，则根据偏差值对磁共振系统的梯度灵敏度因子进行调整，以对磁共振系统的梯度灵敏度进行校准。
上述实施例提供的计算机可读存储介质，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
使用多片层扫描协议对球形模体进行扫描，得到球形模体的三维图像。
上述实施例提供的计算机可读存储介质，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
需要说明的是，本申请所涉及的用户信息(包括但不限于用户设备信息、用户个人信息等)和数据(包括但不限于用于分析的数据、存储的数据、展示的数据等)，均为经用户授权或者经过各方充分授权的信息和数据。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(Read-OnlyMemory，ROM)、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(ReRAM)、磁变存储器(Magnetoresistive Random Access Memory，MRAM)、铁电存储器(Ferroelectric Random Access Memory，FRAM)、相变存储器(Phase Change Memory，PCM)、石墨烯存储器等。易失性存储器可包括随机存取存储器(Random Access Memory，RAM)或外部高速缓冲存储器等。作为说明而非局限，RAM可以是多种形式，比如静态随机存取存储器(Static Random Access Memory，SRAM)或动态随机存取存储器(Dynamic RandomAccess Memory，DRAM)等。本申请所提供的各实施例中所涉及的数据库可包括关系型数据库和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数据库等，不限于此。本申请所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等，不限于此。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请的保护范围应以所附权利要求为准。
Claims (9)
1.一种梯度灵敏度校准方法，其特征在于，所述梯度灵敏度校准方法应用于悬臂床磁共振系统，所述方法包括：
获取球形模体的三维图像；
基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸；
根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准，包括：根据所述拟合尺寸和所述实际尺寸的相似度或偏差值，对所述磁共振系统的梯度灵敏度进行校准；
所述基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸，包括：
利用预设的拟合算法对所述球形模体的三维图像进行椭球拟合拟合出一个椭球，或者，对所述球形模体的三维图像进行重建调整拟合出一个椭球，基于拟合的椭球，得到所述椭球的目标轴的拟合尺寸。
2.根据权利要求1所述的方法，其特征在于，所述根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准，包括：
获取所述拟合尺寸与所述球形模体的目标轴的实际尺寸的偏差值；
根据所述偏差值，对所述磁共振系统的梯度灵敏度因子进行调整，以对所述磁共振系统的梯度灵敏度进行校准。
3.根据权利要求2所述的方法，其特征在于，所述根据所述偏差值，对所述磁共振系统的梯度灵敏度因子进行调整，以对所述磁共振系统的梯度灵敏度进行校准，包括：
确定所述偏差值是否满足预设的偏差范围；
若否，则根据所述偏差值对所述磁共振系统的梯度灵敏度因子进行调整，以对所述磁共振系统的梯度灵敏度进行校准。
4.根据权利要求1所述的方法，其特征在于，所述目标轴包括以磁体成像区域的中心位置为原点，以空间正交的三轴分别作为三轴方向建立的坐标系所对应的三个轴。
5.根据权利要求1所述的方法，其特征在于，所述获取球形模体的三维图像，包括：
使用多片层扫描协议对所述球形模体进行扫描，得到所述球形模体的三维图像。
6.一种梯度灵敏度校准装置，其特征在于，所述梯度灵敏度校准装置应用于悬臂床磁共振系统，所述装置包括：
第一获取模块，用于获取球形模体的三维图像；
第二获取模块，用于基于所述三维图像进行椭球拟合，得到所述椭球的目标轴的拟合尺寸；
校准模块，用于根据所述拟合尺寸和所述球形模体的目标轴的实际尺寸，对磁共振系统的梯度灵敏度进行校准，包括：根据所述拟合尺寸和所述实际尺寸的相似度或偏差值，对所述磁共振系统的梯度灵敏度进行校准；
所述第二获取模块，还用于利用预设的拟合算法对所述球形模体的三维图像进行椭球拟合拟合出一个椭球，或者，对所述球形模体的三维图像进行重建调整拟合出一个椭球，基于拟合的椭球，得到所述椭球的目标轴的拟合尺寸。
7.一种磁共振设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至5中任一项所述的方法的步骤。
8.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至5中任一项所述的方法的步骤。
9.一种计算机程序产品，包括计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至5中任一项所述的方法的步骤。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111525705.1A CN114264997B (zh) | 2021-12-14 | 2021-12-14 | 梯度灵敏度校准方法、装置及磁共振设备 |
| PCT/CN2022/100376 WO2023109064A1 (zh) | 2021-12-14 | 2022-06-22 | 梯度灵敏度校正方法和系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111525705.1A CN114264997B (zh) | 2021-12-14 | 2021-12-14 | 梯度灵敏度校准方法、装置及磁共振设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114264997A CN114264997A (zh) | 2022-04-01 |
| CN114264997B true CN114264997B (zh) | 2024-03-22 |
Family
ID=80826969
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111525705.1A Active CN114264997B (zh) | 2021-12-14 | 2021-12-14 | 梯度灵敏度校准方法、装置及磁共振设备 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN114264997B (zh) |
| WO (1) | WO2023109064A1 (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114264997B (zh) * | 2021-12-14 | 2024-03-22 | 武汉联影生命科学仪器有限公司 | 梯度灵敏度校准方法、装置及磁共振设备 |
| CN117607770B (zh) * | 2024-01-22 | 2024-04-16 | 中国科学技术大学先进技术研究院 | 磁共振图像重建方法、装置、电子设备及存储介质 |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014025886A1 (en) * | 2012-08-09 | 2014-02-13 | Hologic, Inc. | System and method of overlaying images of different modalities |
| CN103761519A (zh) * | 2013-12-20 | 2014-04-30 | 哈尔滨工业大学深圳研究生院 | 一种基于自适应校准的非接触式视线追踪方法 |
| CN106102585A (zh) * | 2015-02-16 | 2016-11-09 | 深圳迈瑞生物医疗电子股份有限公司 | 三维成像数据的显示处理方法和三维超声成像方法及系统 |
| CN106483475A (zh) * | 2015-09-01 | 2017-03-08 | 意法半导体股份有限公司 | 用于磁场感测设备的校准方法、对应系统、装置和计算机程序产品 |
| WO2017063386A1 (zh) * | 2015-10-13 | 2017-04-20 | 上海华测导航技术股份有限公司 | 一种姿态测量系统的精度校准方法 |
| CN106896337A (zh) * | 2015-12-18 | 2017-06-27 | 意法半导体公司 | 用于磁传感器校准的方法 |
| CN107024671A (zh) * | 2016-02-01 | 2017-08-08 | 中国科学院上海高等研究院 | 电子设备的磁力计校准进度可视化方法及系统 |
| CN108037318A (zh) * | 2017-12-07 | 2018-05-15 | 智灵飞（北京）科技有限公司 | 一种基于椭球拟合的无人机加速度计校准方法 |
| CN108195399A (zh) * | 2017-12-22 | 2018-06-22 | 苏州捷研芯纳米科技有限公司 | 用于动态校准磁场传感器的方法和系统 |
| CN108872909A (zh) * | 2018-06-06 | 2018-11-23 | 上海东软医疗科技有限公司 | 磁共振系统的梯度线圈位置校准方法和装置 |
| CN112964278A (zh) * | 2021-03-25 | 2021-06-15 | 北京三快在线科技有限公司 | 确定磁力计校准参数的方法、装置、电子设备及存储介质 |
Family Cites Families (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5545995A (en) * | 1994-11-29 | 1996-08-13 | General Electric Company | Measurement of geometric distortion in MRI images |
| AU2002953540A0 (en) * | 2002-12-24 | 2003-01-16 | The University Of Queensland | Correction of non-linear gradients effects on magnetic resonance imaging |
| JP2008229277A (ja) * | 2007-03-23 | 2008-10-02 | Ge Medical Systems Global Technology Co Llc | 磁気共鳴イメージング装置、磁気共鳴イメージング方法、および、感度分布計測装置 |
| CN101046506B (zh) * | 2007-03-28 | 2011-07-27 | 新奥博为技术有限公司 | 一种基于球谐函数的磁共振图像梯度变形校正方法 |
| US8098068B2 (en) * | 2007-07-20 | 2012-01-17 | Medical College Of Georgia Research Institute, Inc. | Systems and methods for parametric mapping for correction of noise-based systematic bias of DTI metrics, using a DTI mapping phantom |
| JP2009284958A (ja) * | 2008-05-27 | 2009-12-10 | Hitachi Medical Corp | 磁気共鳴イメージング装置及びファントム |
| EP2478383A1 (en) * | 2009-09-14 | 2012-07-25 | INSERM (Institut National de la Santé et de la Recherche Médicale) | Calibration phantom and method for measuring and correcting geometric distorsions in medical images |
| CN104142485B (zh) * | 2013-05-08 | 2017-09-15 | 意法半导体(中国)投资有限公司 | 校准磁传感器的方法和装置 |
| CN103632367B (zh) * | 2013-11-27 | 2016-04-13 | 中国科学技术大学 | 一种基于多组织区域拟合的mri线圈灵敏度估计方法 |
| FR3077877B1 (fr) * | 2018-02-14 | 2020-03-13 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Procede de calibration d'un capteur tri-axe avec selection d'une methode de calibration en fonction de la distribution spatiale des mesures |
| CN111060042A (zh) * | 2019-12-30 | 2020-04-24 | 中国工程物理研究院应用电子学研究所 | 一种基于球模体进行工业ct几何尺寸测量不确定度的方法 |
| CN113782174B (zh) * | 2021-11-12 | 2022-02-18 | 极限人工智能有限公司 | Ct机的校验方法、装置、电子设备及可读存储介质 |
| CN114264997B (zh) * | 2021-12-14 | 2024-03-22 | 武汉联影生命科学仪器有限公司 | 梯度灵敏度校准方法、装置及磁共振设备 |
-
2021
- 2021-12-14 CN CN202111525705.1A patent/CN114264997B/zh active Active
-
2022
- 2022-06-22 WO PCT/CN2022/100376 patent/WO2023109064A1/zh unknown
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2014025886A1 (en) * | 2012-08-09 | 2014-02-13 | Hologic, Inc. | System and method of overlaying images of different modalities |
| CN103761519A (zh) * | 2013-12-20 | 2014-04-30 | 哈尔滨工业大学深圳研究生院 | 一种基于自适应校准的非接触式视线追踪方法 |
| CN106102585A (zh) * | 2015-02-16 | 2016-11-09 | 深圳迈瑞生物医疗电子股份有限公司 | 三维成像数据的显示处理方法和三维超声成像方法及系统 |
| CN106483475A (zh) * | 2015-09-01 | 2017-03-08 | 意法半导体股份有限公司 | 用于磁场感测设备的校准方法、对应系统、装置和计算机程序产品 |
| WO2017063386A1 (zh) * | 2015-10-13 | 2017-04-20 | 上海华测导航技术股份有限公司 | 一种姿态测量系统的精度校准方法 |
| CN106896337A (zh) * | 2015-12-18 | 2017-06-27 | 意法半导体公司 | 用于磁传感器校准的方法 |
| CN107024671A (zh) * | 2016-02-01 | 2017-08-08 | 中国科学院上海高等研究院 | 电子设备的磁力计校准进度可视化方法及系统 |
| CN108037318A (zh) * | 2017-12-07 | 2018-05-15 | 智灵飞（北京）科技有限公司 | 一种基于椭球拟合的无人机加速度计校准方法 |
| CN108195399A (zh) * | 2017-12-22 | 2018-06-22 | 苏州捷研芯纳米科技有限公司 | 用于动态校准磁场传感器的方法和系统 |
| CN108872909A (zh) * | 2018-06-06 | 2018-11-23 | 上海东软医疗科技有限公司 | 磁共振系统的梯度线圈位置校准方法和装置 |
| CN112964278A (zh) * | 2021-03-25 | 2021-06-15 | 北京三快在线科技有限公司 | 确定磁力计校准参数的方法、装置、电子设备及存储介质 |
Non-Patent Citations (2)
| Title |
|---|
| Fangzhou Hu等.An Improved Method for the Magnetometer Calibration Based on Ellipsoid Fitting.《2019 12th International Congress on Image and Signal Processing, BioMedical Engineering and Informatics (CISP-BMEI)》.2020,全文. * |
| 刘轶凡等.基于牛顿迭代和椭球拟合的磁力和惯性传感器校准方法.《仪器仪表学报》.2020,第41卷(第8期),全文. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114264997A (zh) | 2022-04-01 |
| WO2023109064A1 (zh) | 2023-06-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114264997B (zh) | 梯度灵敏度校准方法、装置及磁共振设备 | |
| US20130279786A1 (en) | Rapid parallel reconstruction for arbitrary k-space trajectories | |
| CN107392876B (zh) | 获取归一化校正因子的方法、装置及医学成像方法 | |
| JP2011501200A (ja) | 陽電子射出断層撮影方法及びｐｅｔスキャナ | |
| Zahneisen et al. | Reverse retrospective motion correction | |
| WO2023174333A9 (zh) | 磁共振梯度校正补偿因子的确定方法、磁共振梯度校正方法和装置 | |
| CN114469151A (zh) | 数据校正方法、装置、计算机设备、存储介质和程序产品 | |
| Zahneisen et al. | Fast noniterative calibration of an external motion tracking device | |
| CN107479015A (zh) | 磁共振校准扫描序列配置和图像获取的方法及系统 | |
| US10866298B2 (en) | Low rank and spatial regularization model for magnetic resonance fingerprinting | |
| CN116256681A (zh) | 磁共振成像方法、装置、设备和存储介质 | |
| CN112006713B (zh) | 动态重建、扫描对象定位、计数率显示方法和pet扫描设备 | |
| US20230200653A1 (en) | Diffusion kurtosis imaging method, computer device and storage medium | |
| US10317495B2 (en) | Method and computer use of an MR overview dataset for reconstruction of MR raw data | |
| CN116819412B (zh) | 磁共振系统的校正方法、装置和磁共振系统 | |
| CN116264964A (zh) | 弥散张量成像方法、装置、磁共振设备和存储介质 | |
| CN116973825A (zh) | 磁共振成像系统匀场方法、装置、系统和存储介质 | |
| CN113384260B (zh) | 加速因子调节方法、磁共振成像扫描方法、装置和设备 | |
| CN117765101A (zh) | 磁共振图像重建方法和装置 | |
| CN117392242B (zh) | 成像系统标定方法、装置、计算机设备和存储介质 | |
| CN114299182A (zh) | 扩散峰度系数的确定方法和装置、参数图像的生成方法 | |
| CN117092570A (zh) | 磁共振成像方法、装置和设备 | |
| CN115994874B (zh) | 玻片图像处理方法、装置、玻片、计算机设备和存储介质 | |
| CN116008884A (zh) | 梯度幅值修正方法、装置、设备、存储介质和程序产品 | |
| CN112967297B (zh) | 磁共振氧十七代谢成像方法、装置、存储介质及终端设备 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |