CN109917314A - 磁共振场图确定方法、装置、电子设备及存储介质 - Google Patents

磁共振场图确定方法、装置、电子设备及存储介质 Download PDF

Info

Publication number
CN109917314A
CN109917314A CN201910357668.4A CN201910357668A CN109917314A CN 109917314 A CN109917314 A CN 109917314A CN 201910357668 A CN201910357668 A CN 201910357668A CN 109917314 A CN109917314 A CN 109917314A
Authority
CN
China
Prior art keywords
resonance frequency
reference picture
item
objective function
time
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.)
Granted
Application number
CN201910357668.4A
Other languages
English (en)
Other versions
CN109917314B (zh
Inventor
郑远
叶永泉
丁彧
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai United Imaging Healthcare Co Ltd
Original Assignee
Shanghai United Imaging Healthcare Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai United Imaging Healthcare Co Ltd filed Critical Shanghai United Imaging Healthcare Co Ltd
Priority to CN201910357668.4A priority Critical patent/CN109917314B/zh
Publication of CN109917314A publication Critical patent/CN109917314A/zh
Priority to US16/750,012 priority patent/US11067654B2/en
Application granted granted Critical
Publication of CN109917314B publication Critical patent/CN109917314B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • 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/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/243—Spatial mapping of the polarizing magnetic field
    • 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
    • G01R33/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
    • 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
    • G01R33/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56536—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to magnetic susceptibility variations

Landscapes

  • Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (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)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明实施例公开了一种磁共振场图确定方法、装置、电子设备及存储介质。该方法包括：获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；构建以偏共振频率为自变量的目标函数；其中，所述目标函数包括所述参考图像在各相邻回波时间的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项；根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。本发明实施例的技术方案，避免磁共振场图确定过程中对图像空间连续性以及种子点选取的依赖，提高了求解过程的稳定性，提高了磁共振场图的信噪比。

Description

磁共振场图确定方法、装置、电子设备及存储介质
技术领域
本发明实施例涉及医疗设备技术领域，尤其涉及一种磁共振场图确定方法、装置、电子设备及存储介质。
背景技术
由于磁共振成像能够在无损伤、无电离辐射的情况下得到被检测对象内部的高对比度清晰图像，因此其在诸多领域特别是医学诊断中得到了广泛的应用。
磁共振图像中每个像素点均为复数，其幅值和相位都包含了被检测对象与磁场相互作用的信息。在某些特定的场景下，例如磁化率加权成像、定量磁化率成像、温度成像或者磁共振系统的匀场校准过程中，均需要基于各像素点的相位得到场图信息。
场图的计算通常依赖于不同回波时间采集的图像的相位差，但是由于磁共振图像中每个像素点的相位差由两个回波时间的复数信号相除取幅角得到，因此相位差限定在一个2π周期内(通常为(-π,π])，当相位差超过该取值范围的上下限时，则会以2π周期进行增减。因此，将会造成图像空间相位的跳变，称为相位卷褶。而具有相位卷褶的图像不能直接用于场图的确定。
现有技术通常对两个回波时间对应的相位差图进行去卷褶处理，然后通过去卷褶后的相位差图进行磁共振场图的确定。然而，相位去卷褶算法依赖于真实相位的空间连续性，在空间连续性不好的位置，例如磁共振图像中的眼睛、头皮等，去卷褶算法的效果较差。另外，去卷褶后的相位的仅相对于预先选取的种子点无相位卷褶，而种子点本身的相位与真实相位可能还存在n(n为整数)个2π周期的差值，因此也无法保证去卷褶后绝对相位的准确性，进而影响确定的磁共振场图的信噪比。
发明内容
本发明提供一种磁共振场图确定方法、装置、电子设备及存储介质，以避免磁共振场图确定过程中对图像空间连续性以及种子点选取的依赖，同时提高场图的信噪比。
第一方面，本发明实施例提供了一种磁共振场图确定方法，包括：
获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；
构建以偏共振频率为自变量的目标函数；其中，所述目标函数包括所述参考图像在各相邻回波时间的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项；
根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。
第二方面，本发明实施例还提供了一种磁共振场图确定装置，包括：
参考图像获取模块，用于获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；
目标函数构建模块，用于构建以偏共振频率为自变量的目标函数；其中，所述目标函数包括所述参考图像在各相邻回波时间的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项；
磁共振场图确定模块，用于根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。
第三方面，本发明实施例还提供了一种电子设备，包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如第一方面实施例所提供的一种磁共振场图确定方法。
第四方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如第一方面实施例所提供的一种磁共振场图确定方法。
本发明实施例通过获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中至少两个相邻回波时间的差值不相等；以偏共振频率为自变量，构建包括参考图像在各相邻回波时间的相位偏差构建的相位偏差项和偏共振频率在至少一个变换域下的稀疏性参数构建的系数项约束项的目标函数；根据目标函数确定参考图像对应的目标偏共振频率，并根据目标偏共振频率确定磁共振场图。上述技术方案通过利用多个回波时间的参考图像的信息，避免磁共振场图确定过程中对图像空间连续性以及种子点选取的依赖，提高了求解过程的稳定性，并通过包括时间域下的相位偏差项和至少一个变换域下的稀疏性约束项的目标函数，进行目标偏共振频率的确定，进而确定磁共振场图，兼顾了时间域和变换域下的数据特征，进一步提高了磁共振场图的信噪比。
附图说明
图1A是本发明实施例一中的一种磁共振场图确定方法的流程图；
图1B是本发明实施例一中的一种回波采集序列图；
图2是本发明实施例二中的一种磁共振场图确定方法的流程图；
图3A是本发明实施例三中的一种磁共振场图确定方法的流程图；
图3B是本发明实施例三提供的采用现有技术所确定的磁共振场图；
图3C是本发明实施例三提供的采用一种目标函数确定的磁共振场图；
图3D是本发明实施例三提供的采用另一种目标函数确定的磁共振场图；
图4是本发明实施例四中的一种磁共振场图确定装置的结构图；
图5是本发明实施例五中的一种电子设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1A是本发明实施例一中的一种磁共振场图确定方法的流程图。本发明实施例适用于对磁共振图像进行场图计算的情况，该方法由磁共振场图确定装置执行，该装置通过软件和/或硬件实现，并具体配置于电子设备中。例如可以是磁共振设备。
如图1A所示的一种磁共振场图确定方法，包括：
S110、获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等。
其中，参考图像为采用磁共振设备扫描被检测对象后所得到的磁共振图像；回波时间为施加激发脉冲开始时到回波完成之间的时间间隔。
可选的，可以在电子设备本地、与电子设备相关联的其他存储设备或云端中预先存储参考图像；相应的，可以从电子设备本地、与电子设备相关联的其他存储设备或云端中直接读取至少三个不同回波时间下的扫描同一被检测对象所得到的磁共振图像作为参考图像。
或者可选的，直接控制磁共振设备对被检测对象进行扫描，并根据得到的扫描数据重建后得到的磁共振图像；获取至少三个不同回波时间下对应的磁共振图像作为参考图像。
其中，相邻回波时间可以是连续的两个回波时间，也可以是两个不连续的回波时间；可以是同一扫描过程对应的回波时间，还可以是不同扫描过程对应的回波时间，本发明实施例对此不做任何限定。
其中，所获取的参考图像所满足的条件为：至少两个相邻回波时间的差值不相等，也即各相邻回波时间的差值不完全相等。例如，当第k个回波时间TEk和相邻的第k+1的回波时间TEk+1的差值可以是：dTEk,k+1＝TEk+1-TEk。仅需保证至少存在一个dTEk,k+1≠dTEm,m+1即可，其中k≠m，m为回波时间的序号。
图1B为同一被检测对象在同一扫描过程中进行扫描成像时所对应的回波采集序列图，并以三个不同回波时间为例对所选取的参考图像所满足的条件进行说明。其中TE1、TE2以及TE3对应不同的回波时间，记相邻回波时间TE2与TE1的差值dTE1,2，记相邻回波时间TE3与TE2的差值为dTE3,2，那么TE1、TE2以及TE3需要满足的基本条件为：dTE2,3-dTE1,2＝(TE3-TE2)-(TE2-TE1)≠0，也即dTE2,3≠dTE1,2。
S120、构建以偏共振频率为自变量的目标函数。
其中，所述目标函数包括所述参考图像在各相邻回波时间的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项。
其中，参考图像中各像素点对应的像素值均为复数，并通过复数中的幅值和相位值，共同反映被检测对象与磁场相互作用后的相关信息。
可选的，构建包括各参考图像的相位偏差的相位偏差项，可以是：根据参考图像中各像素点对应的相位值，确定各参考图像的相位偏差；根据各参考图像的相位偏差，构建相位偏差项。其中，各参考图像的相位偏差用于表征不同回波时间对应的参考图像由于磁共振设备的机械误差以及外部环境干扰等因素导致的实际相位值与理想相位值之间的差距。因此，相位偏差项反映了参考图像在时间域上的变化程度。
可选的，构建包括至少一个变换域下的稀疏性参数的稀疏性约束项，可以是：对偏共振频率在至少一个变换域下进行处理得到至少一个稀疏性参数，并根据所得到的至少一个稀疏性参数，构建稀疏性约束项。其中，稀疏性约束项用于表征偏共振频率在至少一个变换域下的稀疏性，反映了偏共振频率在变换域上的变化程度。
可选的，构建以目标偏共振频率为自变量的目标函数，可以是：对所述相位偏差项与所述稀疏性约束项进行加权求和，得到以所述偏共振频率为自变量的目标函数。其中，相位偏差项和稀疏性约束项的权重可以由技术人员根据需要或经验值进行设定，还可以根据大量试验加以确定。
S130、根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。
具体的，在目标函数的函数值满足设定约束条件时，确定满足设定约束条件时所对应的偏共振频率为目标偏共振频率，并根据参考图像中各像素点对应的目标偏共振频率确定磁共振场图。具体的，根据各目标偏共振频率确定磁共振场图可通过现有技术中的任何方式加以实现，再次不再赘述。其中，设定约束条件可以由技术人员根据需要或经验值进行确定。
本发明实施例通过获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中至少两个相邻回波时间的差值不相等；以偏共振频率为自变量，构建包括各参考图像的相位偏差构建的相位偏差项和偏共振频率在至少一个变换域下的稀疏性参数构建的系数项约束项的目标函数；根据目标函数确定参考图像对应的目标偏共振频率，并根据目标偏共振频率确定磁共振场图。上述技术方案通过利用多个回波时间的参考图像的信息，避免磁共振场图确定过程中对图像空间连续性以及种子点选取的依赖，提高了求解过程的稳定性，并通过包括时间域下的相位偏差项和至少一个变换域下的稀疏性约束项的目标函数，进行目标偏共振频率的确定，进而确定磁共振场图，兼顾了时间域和变换域下的数据特征，进一步提高了磁共振场图的信噪比。
需要说明的是，本发明实施例的技术方案可应用于磁共振设备执行多通道数据采集的情况。
示例性地，在磁共振设备进行多通道数据采集得到各通道对应的磁共振图像时，可以通过通道合并的方式加以实现。具体的，在获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像时，具体采用以下步骤进行参考图像的获取：获取不同通道在至少三个不同回波时间下，扫描同一被检测对象所得到的初始图像；将同一回波时间下各通道对应的初始图像合并，得到所述参考图像。
示例性地，在磁共振设备进行多通道数据采集得到各通道对应的磁共振图像时，还可以通过将同一相邻回波时间中各通道对应的相位偏差合并的方式加以实现。具体的，针对同一通道所采集的至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，分别确定各相邻回波时间对应的初始相位偏差；根据各通道对应的同一相邻回波时间对应的参考图像的像素值，确定初始相位偏差的权重值；针对同一相邻回波时间，根据各通道对应的权重值对初始相位偏差进行加权，得到最终的相位偏差。
实施例二
图2是本发明实施例二中的一种磁共振场图确定方法的流程图。本发明实施例在上述各实施例的技术方案的基础上进行了细分优化。
进一步地，将操作“根据所述目标函数确定所述参考图像对应的目标偏共振频率”细化为“根据初始偏共振频率，确定目标函数的函数值最小时所对应的偏共振频率；将函数值最小时所对应的偏共振频率作为目标偏共振频率”，以完善目标偏共振频率的确定机制，同时加快迭代收敛速度，提高算法鲁棒性。
如图2所示的一种磁共振场图确定方法，包括：
S210、获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等。
S220、构建以偏共振频率为自变量的目标函数。
其中，所述目标函数包括所述参考图像在各相邻回波时间的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项。
S230、根据初始偏共振频率，确定目标函数的函数值最小时所对应的偏共振频率。
具体的，根据初始偏共振频率对偏共振频率的频率取值范围进行初步限定；确定初始偏共振频率所限定的频率取值范围内目标函数的函数值最小时所对应的偏共振频率，后续将函数值最小时所对应的偏共振频率确定为目标偏共振频率，用以作为磁共振场图的确定依据。
其中，初始偏共振频率可以由技术人员根据经验值进行设定，还可以采用以下方法加以确定：获取至少三个相邻的初始回波时间下的初始参考图像；根据各初始回波时间的回波时间差值的累计值，以及各初始参考图像的相位差的累计值，确定所述初始偏共振频率。
具体的，根据公式确定初始偏共振频率fstart；
其中，i和j对应初始回波时间的序号；aij为任意实数；TEi、TEj分别为第i和第j个初始回波时间；φi、φj分别为第i和第j个初始回波时间对应的初始参考图像的相位矩阵；Test为预先估计的可避免卷褶的时间差；fstart为初始偏共振频率矩阵。
以三个相邻回波时间为例，示例性地，可以根据以下公式确定初始偏共振频率fstart：
S240、将函数值最小时所对应的偏共振频率作为目标偏共振频率。
S250、根据所述目标偏共振频率确定磁共振场图。
本发明实施例通过将目标偏共振频率的确定步骤细化为根据初始偏共振频率确定目标函数的函数值最小时所对应的偏共振频率，并将函数值最小时对应的偏共振频率作为目标偏共振频率，完善了目标偏共振频率的确定机制，同时通过初始偏共振频率的引入，加快了迭代收敛速度，提高了算法的鲁棒性。
实施例三
图3A是本发明实施例三中的一种磁共振场图确定方法的流程图。本发明实施例在上述各实施例的技术方案的基础上进行了优化改进。
进一步地，将操作“构建以偏共振频率为自变量的目标函数”细化为“根据所述偏共振频率、各相邻回波时间对应的参考图像以及所述相邻回波时间的时间差值，分别确定各相邻回波时间对应的相位偏差，并根据各所述相位偏差构建所述相位偏差项；确定所述偏共振频率在至少一个变换域下的稀疏性参数，并根据各所述稀疏性参数构建所述稀疏性约束项；根据所述相位偏差项和所述稀疏性约束项，得到以所述偏共振频率为自变量的目标函数”，以完善目标函数的构建机制。
如图3A所示的一种磁共振场图确定方法，包括：
S310、获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等。
S320、根据所述偏共振频率、各相邻回波时间对应的参考图像以及所述相邻回波时间的时间差值，分别确定各相邻回波时间对应的相位偏差，并根据各所述相位偏差构建所述相位偏差项。
示例性地，根据所述偏共振频率、各相邻回波时间对应的参考图像以及所述相邻回波时间的时间差值，分别确定各相邻回波时间对应的相位偏差，可以是：根据相邻回波时间的时间差值，对相邻回波时间对应的参考图像的相位差进行估计，得到以偏共振频率为自变量的估计相位差；分别根据所述相邻回波时间的各参考图像的当前像素值，确定实际相位差；根据所述估计相位差和所述实际相位差，分别确定各相邻回波时间的相位偏差。
其中，估计相位差可以理解为在理想状态下，也即忽略磁共振设备的机械误差以及外部环境干扰的情况下，得到的相邻回波时间对应的相位的理想差值。其中，实际相位差可以理解为磁共振设备对待检测对象进行图像扫描后，得到的相邻回波时间对应的相位的实际差值，也即在未忽略磁共振设备的机械误差以及外部环境干扰时的实际扫描环境下，所得到的相邻回波时间对应的相位的实际差值。相应的，相位偏差即为相邻回波时间所获取的参考图像的实际相位差与理想相位差之间的偏差，用以反映参考图像在时间域上的变化程度。
示例性地，根据所述估计相位差和所述实际相位差，分别确定各相邻回波时间的相位偏差，可以是：根据各相邻回波时间对应的估计相位差和实际相位差之间的距离，分别确定各相邻回波时间的相位偏差。其中，距离可以是欧几里得距离、或者马哈拉诺比斯距离等。示例性地，确定各相邻回波时间对应的所述估计相位差与所述实际相位差的欧几里得范数，将得到的结果作为所述相位偏差。
在本发明实施例的一个可选实施方式中，根据各所述相位偏差构建所述相位偏差项，可以是对各相位偏差进行加权求和，得到相位偏差项。
可选的，对相位偏差进行加权时的权重可以由技术人员根据经验值设定，并预先存储在电子设备本地、与电子设备关联的存储设备或云端中；相应的，在对相位偏差进行加权之前，从电子设备本地、与电子设备关联的存储设备或云端中进行获取。
示例性地，可以采用以下公式，分别确定各相邻回波时间对应的相位偏差，并根据各所述相位偏差构建所述相位偏差项：
其中，f为自变量偏共振频率矩阵；k为回波时间的序数，其中k＝1,2,…,n-1，其中n为最大回波时间序数，且n≥3；TEk为第k个回波时间；φk为第k个回波时间对应的参考图像的相位矩阵；dTEk,k+1为第k+1和k两个相邻回波时间的差值；dφk,k+1为第k+1和k两个相邻回波时间时相位矩阵之间的相位差矩阵；G1(f)为相位偏差项。
可以理解的是，为了有效权衡不同相邻回波时间对应的各参考图像在时间域上的变化程度，在在根据所述估计相位差和所述实际相位差，分别确定各相邻回波时间的相位偏差之前，还可以对不同相邻回波时间对应的实际相位差和估计相位差进行加权。
可选的，根据相邻回波时间对应的参考图像的像素值，确定各相邻回波时间对应的相位差权重；根据所述相位差权重，分别对所述估计相位差和所述实际相位差加权。
示例性地，可以根据公式确定各相邻回波时间对应的相位差权重；其中，wk为第k个回波时间对应的相位差权重矩阵，Ak为第k个回波时间对应的参考图像的幅值矩阵。
相应的，根据所述估计相位差和所述实际相位差，分别确定各相邻回波时间的相位偏差，可以根据公式进行确定。
S330、确定所述偏共振频率在至少一个变换域下的稀疏性参数，并根据各所述稀疏性参数构建所述稀疏性约束项。
示例性地，确定所述偏共振频率在至少一个变换域下的稀疏性参数，可以是基于变分法、小波变换、离散傅里叶变换、离散余弦变换、以及有限差分变换中的至少一种对所述偏共振频率进行处理；根据至少一个处理结果，确定稀疏性参数。其中，稀疏性参数用于表征偏共振频率在相应变换域内的稀疏性，反映偏共振频率在变换域中的频率变化程度。
示例性地，根据至少一个处理结果，确定稀疏性参数，可以是根据至少一个处理结果中所包含元素的绝对值之和，确定稀疏性参数。
在本发明实施例的一个可选实施方式中，根据各所述稀疏性参数构建所述稀疏性约束项，可以是对各所述稀疏性参数加权求和，得到所述稀疏性约束项。
可选的，对稀疏性参数加权时的权重可以由技术人员根据经验值设定或根据大量的重复性试验加以确定，并预先存储在电子设备本地、与电子设备关联的存储设备或云端中；相应的，在对稀疏性参数加权之前，从电子设备本地、与电子设备关联的存储设备或云端中进行获取。
示例性地，可以采用以下公式，确定所述偏共振频率在至少一个变换域下的稀疏性参数，并根据各所述稀疏性参数构建所述稀疏性约束项：
其中，j为采用的变换域算法的序数，其中j＝1,2,…,m，其中m为使用的变换域算法的总数量，且m≥1；λj为各变换域对应的稀疏性参数的权重，func()为所采用的变换域算法，可以是变分法、小波变换、离散傅里叶变换、离散余弦变换、或者有限差分变换等；G2(f)为稀疏性约束项。
S340、根据所述相位偏差项和所述稀疏性约束项，得到以所述偏共振频率为自变量的目标函数。
示例性地，将相位偏差项和稀疏性约束项加权求和，得到以偏共振频率为自变量的目标函数。其中，相位偏差项和稀疏性约束项的权重可以由技术人员根据经验值或根据大量试验确定。
具体的，根据公式G(f)＝G1(f)+G2(f)，得到以偏共振频率为自变量的目标函数。
S350、根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。
图3B、图3C和图3D分别为采用现有技术对相位差图进行去卷褶处理所确定的磁共振场图、仅根据相位偏差项所构建的目标函数确定的磁共振场图以及采用本发明实施例的技术方案所得到的磁共振场图。经比对可知，根据相位偏差项所构建的目标函数确定的磁共振场图噪音较低，具有较好的信噪比。然而由于仅考虑时间域上的变化情况，导致图像的局部结构的空间连续性相对较差，在某些信噪比较低的点出现场图异常。而采用本发明实施例的技术方案通过引入变换域对应的稀疏性约束项，利用场图本身的连续性避免了异常点的出现，并进一步提高了磁共振场图的信噪比。
本发明实施例通过将目标函数的构建步骤细化为根据偏共振频率和各相邻回波时间对应的参考图像以及相邻回波时间的时间差值，分别确定各相邻回波时间对应的相位偏差，并根据各相位偏差构建相位偏差项；确定偏共振频率在至少一个变换域下的稀疏性参数，并根据各稀疏性参数构建稀疏性约束项；根据相位偏差项和稀疏性约束项，得到以偏共振频率为自变量的目标函数。采用上述技术方案完善了目标函数的构建步骤，并通过时间域下参考图像的相位变化情况以及至少一个变换域下偏共振频率的变化情况，构建目标函数，使得在根据目标函数确定目标偏共振频率时，能够权衡时间域和变化域多个维度，避免了单独考虑时间域时在磁共振场图中出现异常点的情况，有效的保留了场图的连续性，并进一步提高了磁共振场图的信噪比。
实施例四
图4是本发明实施例四中的一种磁共振场图确定装置的结构图。本发明实施例适用于对磁共振图像进行场图计算的情况，该装置通过软件和/或硬件实现，并具体配置于电子设备中。例如可以是磁共振设备。
如图4所示的一种磁共振场图确定装置，包括：参考图像获取模块410，目标函数构建模块420以及磁共振场图确定模块430。
其中，参考图像获取模块410，用于获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；
目标函数构建模块420，用于构建以偏共振频率为自变量的目标函数；其中，所述目标函数包括所述参考图像在各相邻回波时间的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项；
磁共振场图确定模块430，用于根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。
本发明实施例通过参考图像获取模块获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中至少两个相邻回波时间的差值不相等；通过目标函数构建模块以偏共振频率为自变量，构建包括参考图像在各相邻回波时间的相位偏差构建的相位偏差项和偏共振频率在至少一个变换域下的稀疏性参数构建的系数项约束项的目标函数；通过磁共振场图确定模块根据目标函数确定参考图像对应的目标偏共振频率，并根据目标偏共振频率确定磁共振场图。上述技术方案通过利用多个回波时间的参考图像的信息，避免磁共振场图确定过程中对图像空间连续性以及种子点选取的依赖，提高了求解过程的稳定性，并通过包括时间域下的相位偏差项和至少一个变换域下的稀疏性约束项的目标函数，进行目标偏共振频率的确定，进而确定磁共振场图，兼顾了时间域和变换域下的数据特征，进一步提高了磁共振场图的信噪比。
进一步地，目标函数构建模块420，包括：
相位偏差项构建单元，用于根据所述偏共振频率、各相邻回波时间对应的参考图像以及所述相邻回波时间的时间差值，分别确定各相邻回波时间对应的相位偏差，并根据各所述相位偏差构建所述相位偏差项；
稀疏性约束项构建单元，用于确定所述偏共振频率在至少一个变换域下的稀疏性参数，并根据各所述稀疏性参数构建所述稀疏性约束项；
目标函数得到单元，用于根据所述相位偏差项和所述稀疏性约束项，得到以所述偏共振频率为自变量的目标函数。
进一步地，相位偏差项构建单元，在执行根据所述偏共振频率、各相邻回波时间对应的参考图像以及所述相邻回波时间的时间差值，分别确定各相邻回波时间对应的相位偏差时，具体包括：
估计相位差确定子单元，用于根据相邻回波时间的时间差值，对相邻回波时间对应的参考图像的相位差进行估计，得到以偏共振频率为自变量的估计相位差；
实际相位差确定子单元，用于分别根据所述相邻回波时间的各参考图像的当前像素值，确定实际相位差；
相位偏差确定子单元，用于根据所述估计相位差和所述实际相位差，分别确定各相邻回波时间的相位偏差。
进一步地，相位偏差确定子单元，具体用于：
确定各相邻回波时间对应的所述估计相位差与所述实际相位差的欧几里得范数，将得到的结果作为所述相位偏差。
进一步地，稀疏性约束项构建单元，在执行确定所述偏共振频率在至少一个变换域下的稀疏性参数时，包括：
偏共振频率处理子单元，用于基于变分法、小波变换、离散傅里叶变换、离散余弦变换、以及有限差分变换中的至少一种对所述偏共振频率进行处理；
稀疏性参数确定子单元，用于根据至少一个处理结果，确定稀疏性参数。
进一步地，稀疏性参数确定子单元，具体用于：
根据至少一个处理结果中所包含元素的绝对值之和，确定稀疏性参数。
进一步地，相位偏差项构建单元，，还包括：
权重确定子单元，用于根据相邻回波时间对应的参考图像的像素值，确定各相邻回波时间对应的相位差权重；
相位差加权子单元，用于根据所述相位差权重，分别对所述估计相位差和所述实际相位差加权。
进一步地，相位偏差项构建单元，在执行根据各所述相位偏差构建所述相位偏差项时，包括：
相位偏差项得到子单元，用于对各所述相位偏差加权求和，得到所述相位偏差项。
进一步地，稀疏性约束项构建单元，在执行根据各所述稀疏性参数构建所述稀疏性约束项时，包括：
稀疏性约束项得到子单元，用于对各所述稀疏性参数加权求和，得到所述稀疏性约束项。
进一步地，目标函数得到单元，包括：
目标函数得到子单元，用于对所述相位偏差项与所述稀疏性约束项进行加权求和，得到以所述偏共振频率为自变量的目标函数。
进一步地，磁共振场图确定模块430，在执行根据所述目标函数确定所述参考图像对应的目标偏共振频率时，具体用于：
根据初始偏共振频率，确定目标函数的函数值最小时所对应的偏共振频率；
将函数值最小时所对应的偏共振频率作为目标偏共振频率。
进一步地，该装置还包括，初始偏共振频率确定模块，具体用于：
在确定目标函数的函数值最小时所对应的偏共振频率之前，获取三个相邻的初始回波时间下的初始参考图像；
根据各初始回波时间的回波时间差值的累计值，以及各初始参考图像的相位差的累计值，确定所述初始偏共振频率。
上述磁共振场图确定装置可执行本发明任意实施例所提供的磁共振场图确定方法，具备执行磁共振场图确定方法相应的功能模块和有益效果。
实施例五
图5是本发明实施例五中的一种电子设备的结构示意图。该设备包括处理器510以及存储装置520。
一个或多个处理器510；
存储装置520，用于存储一个或多个程序。
该电子设备，还包括输入装置530，用于获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；
输出装置540，用于显示参考图像和/或磁共振场图。
图5中以一个处理器510为例，该设备中的处理器510和存储装置520通过总线或其他方式连接，输入装置530和输出装置550也可以通过总线或其他方式进行连接，图5中以通过总线连接为例。
在本实施例中，电子设备中的处理器510可以从输入装置530中获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；还可以构建以偏共振频率为自变量的目标函数；其中，所述目标函数包括各参考图像的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项；还可以根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图；还可以控制输出装置540进行参考图像和/或磁共振场图的显示。
该电子设备中的存储装置520作为一种计算机可读存储介质，可用于存储一个或多个程序，所述程序可以是软件程序、计算机可执行程序以及模块，如本发明实施例中磁共振场图确定方法对应的程序指令/模块(例如，附图4所示的参考图像获取模块410，目标函数构建模块420以及场图确定模块430)。处理器510通过运行存储在存储装置520中的软件程序、指令以及模块，从而执行电子设备的各种功能应用以及数据处理，即实现上述方法实施例中的磁共振场图确定方法。
存储装置520可包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储数据等(如上述实施例中的参考图像、目标偏共振频率以及磁共振场图等)。此外，存储装置520可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储装置520可进一步包括相对于处理器510远程设置的存储器，这些远程存储器可以通过网络连接至服务器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
此外，本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，该程序被磁共振场图确定装置执行时实现本发明实施提供的磁共振场图确定方法，该方法包括：获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；构建以偏共振频率为自变量的目标函数；其中，所述目标函数包括各参考图像的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项；根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。
通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(Read-Only Memory,ROM)、随机存取存储器(RandomAccess Memory,RAM)、闪存(FLASH)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的磁共振场图确定方法。
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。

Claims (14)

1.一种磁共振场图确定方法，其特征在于，包括：
获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；
构建以偏共振频率为自变量的目标函数；其中，所述目标函数包括所述参考图像在各相邻回波时间的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项；
根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。
2.根据权利要求1所述的方法，其特征在于，构建以偏共振频率为自变量的目标函数，包括：
根据所述偏共振频率、各相邻回波时间对应的参考图像以及所述相邻回波时间的时间差值，分别确定各相邻回波时间对应的相位偏差，并根据各所述相位偏差构建所述相位偏差项；
确定所述偏共振频率在至少一个变换域下的稀疏性参数，并根据各所述稀疏性参数构建所述稀疏性约束项；
根据所述相位偏差项和所述稀疏性约束项，得到以所述偏共振频率为自变量的目标函数。
3.根据权利要求2所述的方法，其特征在于，根据所述偏共振频率、各相邻回波时间对应的参考图像以及所述相邻回波时间的时间差值，分别确定各相邻回波时间对应的相位偏差，包括：
根据相邻回波时间的时间差值，对相邻回波时间对应的参考图像的相位差进行估计，得到以偏共振频率为自变量的估计相位差；
分别根据所述相邻回波时间的各参考图像的当前像素值，确定实际相位差；
根据所述估计相位差和所述实际相位差，分别确定各相邻回波时间的相位偏差。
4.根据权利要求3所述的方法，其特征在于，根据所述估计相位差和所述实际相位差，分别确定各相邻回波时间的相位偏差，包括：
确定各相邻回波时间对应的所述估计相位差与所述实际相位差的欧几里得范数，将得到的结果作为所述相位偏差。
5.根据权利要求2所述的方法，其特征在于，确定所述偏共振频率在至少一个变换域下的稀疏性参数，包括：
基于变分法、小波变换、离散傅里叶变换、离散余弦变换、以及有限差分变换中的至少一种对所述偏共振频率进行处理；
根据至少一个处理结果，确定稀疏性参数。
6.根据权利要求5所述的方法，其特征在于，根据至少一个处理结果，确定稀疏性参数，包括：
根据至少一个处理结果中所包含元素的绝对值之和，确定稀疏性参数。
7.根据权利要求3所述的方法，其特征在于，在根据所述估计相位差和所述实际相位差，分别确定各相邻回波时间的相位偏差之前，还包括：
根据相邻回波时间对应的参考图像的像素值，确定各相邻回波时间对应的相位差权重；
根据所述相位差权重，分别对所述估计相位差和所述实际相位差加权。
8.根据权利要求2所述的方法，其特征在于，根据各所述相位偏差构建所述相位偏差项，包括：
对各所述相位偏差加权求和，得到所述相位偏差项；
和/或，
根据各所述稀疏性参数构建所述稀疏性约束项，包括：
对各所述稀疏性参数加权求和，得到所述稀疏性约束项。
9.根据权利要求2所述的方法，其特征在于，根据所述相位偏差项和所述稀疏性约束项，得到以所述偏共振频率为自变量的目标函数，包括：
对所述相位偏差项与所述稀疏性约束项进行加权求和，得到以所述偏共振频率为自变量的目标函数。
10.根据权利要求1-9任一项所述的方法，其特征在于，根据所述目标函数确定所述参考图像对应的目标偏共振频率，包括：
根据初始偏共振频率，确定目标函数的函数值最小时所对应的偏共振频率；
将函数值最小时所对应的偏共振频率作为目标偏共振频率。
11.根据权利要求10所述的方法，其特征在于，在确定目标函数的函数值最小时所对应的偏共振频率之前，还包括：
获取至少三个相邻的初始回波时间下的初始参考图像；
根据各初始回波时间的回波时间差值的累计值，以及各初始参考图像的相位差的累计值，确定所述初始偏共振频率。
12.一种磁共振场图确定装置，其特征在于，包括：
参考图像获取模块，用于获取至少三个不同回波时间下扫描同一被检测对象所得到的参考图像，其中，至少两个相邻回波时间的差值不相等；
目标函数构建模块，用于构建以偏共振频率为自变量的目标函数；其中，所述目标函数包括所述参考图像在各相邻回波时间的相位偏差构建的相位偏差项和所述偏共振频率在至少一个变换域下的稀疏性参数构建的稀疏性约束项；
磁共振场图确定模块，用于根据所述目标函数确定所述参考图像对应的目标偏共振频率，并根据所述目标偏共振频率确定磁共振场图。
13.一种电子设备，其特征在于，包括：
一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如权利要求1-11任一项所述的一种磁共振场图确定方法。
14.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现如权利要求1-11任一项所述的一种磁共振场图确定方法。
CN201910357668.4A 2019-04-29 2019-04-29 磁共振场图确定方法、装置、电子设备及存储介质 Active CN109917314B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201910357668.4A CN109917314B (zh) 2019-04-29 2019-04-29 磁共振场图确定方法、装置、电子设备及存储介质
US16/750,012 US11067654B2 (en) 2019-04-29 2020-01-23 Systems and methods for determining field map

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910357668.4A CN109917314B (zh) 2019-04-29 2019-04-29 磁共振场图确定方法、装置、电子设备及存储介质

Publications (2)

Publication Number Publication Date
CN109917314A true CN109917314A (zh) 2019-06-21
CN109917314B CN109917314B (zh) 2021-05-07

Family

ID=66978965

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910357668.4A Active CN109917314B (zh) 2019-04-29 2019-04-29 磁共振场图确定方法、装置、电子设备及存储介质

Country Status (2)

Country Link
US (1) US11067654B2 (zh)
CN (1) CN109917314B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11435424B1 (en) * 2021-08-27 2022-09-06 Spintech, Inc. Systems and methods of constrained reconstruction of images with white noise

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116029963A (zh) * 2021-10-26 2023-04-28 深圳先进技术研究院 血管反应性的评价方法、装置、电子设备及存储介质

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101034149A (zh) * 2007-04-19 2007-09-12 新奥博为技术有限公司 一种基于预扫描测量校正磁共振信号相位偏移的方法
JP2007301159A (ja) * 2006-05-11 2007-11-22 Utsunomiya Univ 磁気共鳴画像化装置、磁気共鳴画像化方法、演算処理プログラム及びそれを記録した情報記録媒体
CN1530072B (zh) * 2003-03-12 2010-05-05 Ge医药系统环球科技公司 核磁共振成像设备和核磁共振成像方法
CN105392421A (zh) * 2013-08-07 2016-03-09 株式会社日立制作所 磁共振摄影装置以及水脂肪分离方法
US20160202335A1 (en) * 2012-04-20 2016-07-14 University Of Virginia Patent Foundation Systems and methods for reduced off-resonance blurring in spiral imaging
CN105809661A (zh) * 2014-12-30 2016-07-27 中国科学院深圳先进技术研究院 磁共振成像的图像水脂分离方法和系统
CN107037386A (zh) * 2016-11-01 2017-08-11 上海联影医疗科技有限公司 一种平面回波成像方法以及系统
CN107358588A (zh) * 2017-07-12 2017-11-17 上海联影医疗科技有限公司 相位去卷褶方法、磁共振成像扫描方法及磁共振成像设备

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1530072B (zh) * 2003-03-12 2010-05-05 Ge医药系统环球科技公司 核磁共振成像设备和核磁共振成像方法
JP2007301159A (ja) * 2006-05-11 2007-11-22 Utsunomiya Univ 磁気共鳴画像化装置、磁気共鳴画像化方法、演算処理プログラム及びそれを記録した情報記録媒体
CN101034149A (zh) * 2007-04-19 2007-09-12 新奥博为技术有限公司 一种基于预扫描测量校正磁共振信号相位偏移的方法
US20160202335A1 (en) * 2012-04-20 2016-07-14 University Of Virginia Patent Foundation Systems and methods for reduced off-resonance blurring in spiral imaging
CN105392421A (zh) * 2013-08-07 2016-03-09 株式会社日立制作所 磁共振摄影装置以及水脂肪分离方法
CN105809661A (zh) * 2014-12-30 2016-07-27 中国科学院深圳先进技术研究院 磁共振成像的图像水脂分离方法和系统
CN107037386A (zh) * 2016-11-01 2017-08-11 上海联影医疗科技有限公司 一种平面回波成像方法以及系统
CN107358588A (zh) * 2017-07-12 2017-11-17 上海联影医疗科技有限公司 相位去卷褶方法、磁共振成像扫描方法及磁共振成像设备

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11435424B1 (en) * 2021-08-27 2022-09-06 Spintech, Inc. Systems and methods of constrained reconstruction of images with white noise

Also Published As

Publication number Publication date
CN109917314B (zh) 2021-05-07
US11067654B2 (en) 2021-07-20
US20200341091A1 (en) 2020-10-29

Similar Documents

Publication Publication Date Title
Li et al. Joint image filtering with deep convolutional networks
US10915990B2 (en) Systems and methods for denoising medical images with deep learning network
US20210295474A1 (en) Methods and system for selective removal of streak artifacts and noise from images using deep neural networks
CN104155623B (zh) 用于自动确定组织种类的磁场反转时间的方法和系统
US20080107319A1 (en) Practical Image Reconstruction for Magnetic Resonance Imaging
CN101933045A (zh) 短轴后期增强心脏mri的自动3d分割
CN112396560A (zh) 用于使用深度神经网络对医疗图像进行去模糊的系统和方法
US8314792B2 (en) Prediction of cardiac shape by a motion model
CN109917314A (zh) 磁共振场图确定方法、装置、电子设备及存储介质
WO2020119581A1 (zh) 磁共振参数成像方法、装置、设备及存储介质
CN105518477A (zh) 使用空间自适应正则化以用于图像重建的mri
WO2023077942A1 (zh) 超声图像的处理方法和超声设备
CN105849804A (zh) 过滤噪声的计算高效方法
Tan et al. Pixelwise Estimation of Signal‐Dependent Image Noise Using Deep Residual Learning
US11426094B2 (en) Methods for iterative reconstruction of medical images using primal-dual optimization with stochastic dual variable updating
Yancheng et al. RED-MAM: A residual encoder-decoder network based on multi-attention fusion for ultrasound image denoising
WO2023027804A1 (en) Denoising depth data of low-signal pixels
Hanika et al. Edge‐Optimized À‐Trous Wavelets for Local Contrast Enhancement with Robust Denoising
CN103784164A (zh) 超声信号的预处理方法及系统
WO2024093083A1 (zh) 一种基于变分自编码器的磁共振加权图像合成方法和装置
JP2021089493A (ja) 情報処理装置およびその学習方法
CN109839608B (zh) 磁共振场图确定方法、装置、电子设备及存储介质
Haq et al. Block-based compressed sensing of MR images using multi-rate deep learning approach
Chilukuri et al. Analysing Of Image Quality Computation Models Through Convolutional Neural Network
CN108629798A (zh) 基于gpu的图像快速配准方法

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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant after: Shanghai Lianying Medical Technology Co., Ltd

Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258

Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd.

GR01 Patent grant
GR01 Patent grant