CN113900056A - 流速测量方法、装置及存储介质 - Google Patents
流速测量方法、装置及存储介质 Download PDFInfo
- Publication number
- CN113900056A CN113900056A CN202111209652.2A CN202111209652A CN113900056A CN 113900056 A CN113900056 A CN 113900056A CN 202111209652 A CN202111209652 A CN 202111209652A CN 113900056 A CN113900056 A CN 113900056A
- Authority
- CN
- China
- Prior art keywords
- flow velocity
- fluid
- flow rate
- phase
- contrast image
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000012530 fluid Substances 0.000 claims abstract description 83
- 230000008859 change Effects 0.000 claims abstract description 35
- 238000000691 measurement method Methods 0.000 claims abstract description 26
- 238000004590 computer program Methods 0.000 claims description 15
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims description 11
- 230000002708 enhancing effect Effects 0.000 claims description 7
- 238000009738 saturating Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 abstract description 24
- 238000003384 imaging method Methods 0.000 abstract description 11
- 238000000338 in vitro Methods 0.000 abstract description 5
- 230000001965 increasing effect Effects 0.000 abstract description 5
- 239000000523 sample Substances 0.000 description 37
- 238000002595 magnetic resonance imaging Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 8
- 238000004891 communication Methods 0.000 description 6
- 230000005284 excitation Effects 0.000 description 6
- 230000005415 magnetization Effects 0.000 description 6
- 239000012472 biological sample Substances 0.000 description 5
- 238000009792 diffusion process Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 4
- 229910052739 hydrogen Inorganic materials 0.000 description 4
- 239000001257 hydrogen Substances 0.000 description 4
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 4
- 230000007170 pathology Effects 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000004907 flux Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- IXPNQXFRVYWDDI-UHFFFAOYSA-N 1-methyl-2,4-dioxo-1,3-diazinane-5-carboximidamide Chemical compound CN1CC(C(N)=N)C(=O)NC1=O IXPNQXFRVYWDDI-UHFFFAOYSA-N 0.000 description 1
- 229920001661 Chitosan Polymers 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 125000002057 carboxymethyl group Chemical group [H]OC(=O)C([H])([H])[*] 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 238000013399 early diagnosis Methods 0.000 description 1
- 210000003722 extracellular fluid Anatomy 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 239000000499 gel Substances 0.000 description 1
- 239000000017 hydrogel Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001575 pathological effect Effects 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 239000002356 single layer Substances 0.000 description 1
- 235000010413 sodium alginate Nutrition 0.000 description 1
- 229940005550 sodium alginate Drugs 0.000 description 1
- 239000000661 sodium alginate Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Heart & Thoracic Surgery (AREA)
- Veterinary Medicine (AREA)
- High Energy & Nuclear Physics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Signal Processing (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请提供一种流速测量方法、装置及存储介质，其中方法包括：在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；基于所述相位对比图像，确定所述流体的流速大小。本申请实施例提供的流速测量方法、装置及存储介质，通过在离体的新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像，根据相位变化确定流体的流速大小，通过增加流速编码梯度之间的间隔，实现了高空间分辨率和成像效率的新鲜样品流速的定量测量。
Description
技术领域
本申请涉及磁共振成像技术领域，尤其涉及一种流速测量方法、装置及存储介质。
背景技术
磁共振成像(Magnetic Resonance Imaging，MRI)中除了幅值图外还可以得到相位分布图(Phase Map，PC)，相位分布图除用于检查主磁场的均匀性外，还可以用于观察病灶区磁化率变化，测量流体的流速。
小型低场核磁仪一般检测非生物样品和固定后的生物样品的流速，通过长时间重复扫描来提高信号及灵敏度。
但是，这种测量流速的方案并不适用于离体的新鲜样品的流速的测量，使用小型低场核磁仪检测离体的新鲜样品的流速时，分辨率和灵敏度低，扫描时间长，测量结果具有显著的误差。
发明内容
本申请提供一种流速测量方法、装置及存储介质，用以解决现有技术中无法精确测量离体的新鲜样品的流速的技术问题。
本申请提供一种流速测量方法，包括：
在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；
基于所述相位对比图像，确定所述流体的流速大小。
根据本申请提供的一种流速测量方法，所述基于所述相位对比图像，确定所述流体的流速大小，包括：
基于所述相位对比图像进行解卷绕，获得解卷绕相位图像；
基于所述解卷绕相位图像，确定流体的流速大小。
根据本申请提供的一种流速测量方法，所述在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像，包括：
在所述流速方向上施加第一对流速编码梯度，得到第一相位变化；
在所述流速方向上施加第二对流速编码梯度，得到第二相位变化；所述第一对流速编码梯度的脉冲顺序和所述第二对流速编码梯度的脉冲顺序之间呈镜像对称；
根据所述第一相位变化和所述第二相位变化确定所述流速方向的相位对比图像。
根据本申请提供的一种流速测量方法，所述第一对流速编码梯度为大小相等、持续时间相同但方向相反的双极梯度；所述第二对流速编码梯度为大小相等，持续时间相同但方向相反的双极梯度。
根据本申请提供的一种流速测量方法，所述方法还包括：
增强所述新鲜样品中所述流体的磁共振信号。
根据本申请提供的一种流速测量方法，所述增强所述新鲜样品中所述流体的磁共振信号，包括：
通过辐照激发所述流体的电子自旋以及电子质子交换，饱和所述流体的自旋磁矩；
关闭辐照并发射射频脉冲，激发所述流体的磁共振信号。
本申请还提供一种流速测量装置，包括：
获取模块，用于在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；
确定模块，用于基于所述相位对比图像，确定所述流体的流速大小。
本申请还提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现如上述任一种所述流速测量方法的步骤。
本申请还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现如上述任一种所述流速测量方法的步骤。
本申请还提供一种计算机程序产品，包括计算机程序，所述计算机程序被处理器执行时实现如上述任一种所述流速测量方法的步骤。
本申请实施例提供的流速测量方法、装置及存储介质，通过在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像，根据相位变化确定流体的流速大小，通过增加流速编码梯度之间的间隔，实现了高空间分辨率和成像效率的新鲜样品流速的定量测量。
附图说明
为了更清楚地说明本申请或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本申请实施例提供的流速测量方法的流程示意图；
图2是本申请实施例提供的相位分布序列图；
图3是本申请实施例提供的流速测量系统的结构示意图；
图4是本申请实施例提供的水模型流动仿体的结构示意图；
图5是本申请实施例提供的流速测量装置的结构示意图；
图6是本申请实施例提供的电子设备的结构示意图。
具体实施方式
为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请中的附图，对本申请中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
以磁共振成像为代表的现代医学影像技术由单一的结构成像发展到结构-定量成像。随着磁共振成像技术的发展进步，获取和分析更加丰富的流体信息，应用于疾病的早期诊断、精准诊断、预后等。磁共振成像技术影像学模态和病理学的密切结合和深度融合得到医学领域的广泛认同和大力推进。
对于离体的新鲜样品进行磁共振成像，可以将影像学和病理学连接起来，为传统病理学检测补充流体行为方面的关键信息，弥补脱水固定生物样品后流体信息的缺失。
磁共振成像得到的相位分布图可以用于观察病灶区磁化率变化，测量流体的流速。相位的周期是360°或者2π弧度，一般相位差值
的表达式为：
其中，
表示相位值，n表示整数，Z表示整数集合。
小型低场核磁仪一般检测非生物样品和固定后的生物样品的流速，通过长时间重复扫描来提高信号及灵敏度。但是，这种测量流速的方案并不适用于离体的新鲜样品的流速的测量。
离体的新鲜样品的流速测量范围不同于传统医用磁共振仪的测量范围，其相位信息的处理需要更加的精准。直接使用原始相位计算流速将导致显著的误差。
为解决现有技术中存在的上述问题，本申请提供一种流速测量方法、装置及存储介质。
图1是本申请实施例提供的流速测量方法的流程示意图，如图1所示，本申请实施例提供一种流速测量方法，包括：
步骤101、在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像。
具体来说，磁共振图像的空间信息是用梯度磁场形成的，人为地使磁场强度在空间上形成有规律的分布。本申请采用梯度回波序列，在离体的新鲜样品中流体的流速方向上施加流速编码梯度，可以获取流速方向的相位对比图像。
新鲜样品中的流体可以为水、血液、组织液等。离体的新鲜样品相较于非生物样品和处理后的生物样品，在流速测量过程中需要对平面内的流动更加敏感。通过施加流速编码梯度，利用样品中流体的流动所导致的宏观横向磁化矢量的相位变化来抑制组织背景、突出血管信号等。
在新鲜样品中流体的流速方向上施加一对双极的流速编码梯度，双极梯度由一对大小相等，持续时间相同但方向相反的梯度脉冲组成；经过这次双极梯度脉冲激发后，静止的质子受到大小相等但方向相反的脉冲作用，相位变化为零，而运动的质子在梯度场中的位置发生变化，产生了相位变化。如此将静止与运动的质子的相位差异显示出来。
然后再施加一对双极的流速编码梯度，与第一对双极的流速编码梯度的脉冲顺序之间呈镜像对称，再次产生相位变化。将两次相位变化作差，得到一个相位差值Δφ。
磁共振成像系统采集经过上述编码的磁共振信号，对采集到的磁共振信号进行解码，即进行图像重构，得到流体流速方向截面的磁共振图像。根据磁共振图像上提取相位信号的变化可以得到相位对比图像。
步骤102、基于所述相位对比图像，确定所述流体的流速大小。
具体来说，通过相位对比图像中的相位信息，可以对样品中流体的流速进行测量。
流动质子群在施加双极梯度场期间积聚的相位变化与其流速有关，流速越快，相位变化越明显。通过施加两对脉冲顺序呈镜像对称的双极的流速编码梯度后，得到两次相位变化。
将两次相位变化作差，得到的相位差值Δφ的表达式为：
Δφ＝2γ·m·Δ·v
其中，γ表示磁旋比，m表示设置的流速编码梯度场的面积，m＝∫Gdt，G表示流速编码梯度场的幅值，t表示时间，Δ表示流速编码梯度之间的间隔，v表示流体的流速。
由此可以得到流体的流速大小v的表达式为：
相位对比图像中相位值的取值区间为(-π，π]，相位测量流速的灵敏度的表达式为：
根据该表达式可以得到，测量的流速的最小值由流速编码梯度场的面积m和流速编码梯度场之间的间隔Δ这两个参量的乘积决定。
同时，磁共振信号强度还会受到流速编码梯度场的影响，导致信号衰减。磁共振信号强度S的表达式为：
S＝e-b·DS0
其中，e表示常量，b表示扩散因子，D表示扩散系数，S0表示b的取值为0时采集到的磁共振信号强度。扩散因子b的表达式为：
其中，δ表示磁场施加时间，对于梯形编码梯度场，其面积m＝Gδ。在本申请实施例中，通过在信号采集过程中增大流速编码梯度之间的间隔Δ，在相同的相位测量流速的灵敏度下，使扩散因子b的取值最小化，从而达到减轻信号损失，提高流速测量的检测精度的效果。
本申请实施例提供的流速测量方法，通过在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像，根据相位变化确定流体的流速大小，通过增加流速编码梯度之间的间隔，实现了高空间分辨率和成像效率的新鲜样品流速的定量测量。
可选地，所述基于所述相位对比图像，确定所述流体的流速大小，包括：
基于所述相位对比图像进行解卷绕，获得解卷绕相位图像；
基于所述解卷绕相位图像，确定流体的流速大小。
具体来说，在一个相位像的一个周期内，相位灰度值随着角度线性变化，当相位角达到周期内的最大阈值时，灰度值急剧变化至最小灰度值，此时产生相位数据溢出，发生卷绕(Wrapping Around)现象。
在本申请实施例中，相位对比图像中相位值的取值区间为(-π，π]，当流速太大时，相对应的相位大于π，就会发生卷绕，在相位对比图像中反映为一个负的相位。
将相位对比图像在(-π，π]边界跳变产生的卷绕修改为连续变化的灰度图像，获得解卷绕相位图像，并根据流体的流速v的表达式确定新鲜样品中流体的流速。流速v的表达式为：
其中，Δφ表示相位差值，γ表示磁旋比，m表示设置的流速编码梯度场的面积，Δ表示流速编码梯度之间的间隔。
本申请实施例提供的流速测量方法，通过对原始的相位对比图像进行解卷绕，使得相位在π处不发生跳变，得到反应真实的相位变化的解卷绕相位图像，进而使得流速测量更加精准。
可选地，所述在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像，包括：
在所述流速方向上施加第一对流速编码梯度，得到第一相位变化；
在所述流速方向上施加第二对流速编码梯度，得到第二相位变化；所述第一对流速编码梯度的脉冲顺序和所述第二对流速编码梯度的脉冲顺序之间呈镜像对称；
根据所述第一相位变化和所述第二相位变化确定所述流速方向的相位对比图像。
具体来说，磁共振相位对比成像是基于流体内质子相位变化这一原理进行成像，在本申请实施例中，采用梯度回波序列，在流体的流速方向上施加第一对流速编码梯度，第一对流速编码梯度由一对大小相等，持续时间相同但方向相反的双极梯度脉冲组成。
经过第一对流速编码梯度脉冲激发后，静止的质子受到大小相等、方向相反的脉冲作用，相位变化为零；而运动的质子在梯度场中的位置发生了变化，产生了第一相位变化。静止与运动的质子的相位差异显现出来。
然后再在流体的流速方向上施加第二对流速编码梯度，得到第二相位变化。第二对流速编码梯度由一对大小相等，持续时间相同但方向相反的双极梯度脉冲组成，与第一对流速编码梯度的脉冲顺序呈镜像对称。通过这一脉冲序列，得到两个相位变化，将这两个相位变化作差，得到相位差值Δφ。
在磁共振成像的坐标系中，顺主磁场方向的为Z轴或纵轴，垂直于主磁场方向的平面为XY平面或者横平面。在本申请实施例中，将磁共振信号从横平面激发到纵平面后，使用一个180°脉冲，使得纵向磁化倾斜180°，在纵平面对信号进行回聚，之后再激发到横平面进行成像。
图2是本申请实施例提供的相位分布序列图，如图2所示，水平方向表示时间，竖直方向表示梯度幅值，第一条轴线表示射频发射轴，第二条轴线表示选层梯度轴，第三条轴线表示扩散编码梯度轴，矩形框内表示信号采集结果，Δ表示流速编码梯度之间的间隔，δ表示磁场施加时间。
本申请实施例将流速编码梯度融合得到纵平面回聚回波序列中。以第三个脉冲为中心，在第二个和第四个脉冲两侧使用面积形状相同的梯度场进行相位编码，得到第一相位变化；第二次通过虚线部分标记，使用与第一次镜像对称的梯度场进行相位编码，得到第二相位变化，两次相位变化作差得到相位差值。
本申请实施例提供的流速测量方法，通过在新鲜样品中流体的流速方向上施加两对镜像对称的流速编码梯度，获取流速方向的相位对比图像，根据相位变化确定流体的流速大小，通过增加流速编码梯度之间的间隔，实现了高空间分辨率和成像效率的新鲜样品流速的定量测量。
可选地，所述第一对流速编码梯度为大小相等、持续时间相同但方向相反的双极梯度；所述第二对流速编码梯度为大小相等，持续时间相同但方向相反的双极梯度。
可选地，所述方法还包括：
增强所述新鲜样品中所述流体的磁共振信号。
可选地，所述增强所述新鲜样品中所述流体的磁共振信号，包括：
通过辐照激发所述流体的电子自旋以及电子质子交换，饱和所述流体的自旋磁矩；
关闭辐照并发射射频脉冲，激发所述流体的磁共振信号。
具体来说，生物组织中含有大量原子核，原子核由中子和质子组成，其中以氢核最多。氢原子核的质子具有自旋的特性，带正电，产生磁矩，绕自旋轴旋转。将氢原子放置在外加的恒定磁场中，其自旋轴将按照外加磁场磁力线的方向有序排列。在这种状态下，用特定频率的射频脉冲进行激发，氢原子核吸收一定能量而共振，即发生磁共振现象。
停止发射射频脉冲后，被激发的氢原子核把吸收的能量逐步释放出来，其相位和能级都恢复到被激发前的状态，这个恢复过程被称为弛豫过程，恢复到原来的平衡状态所需要的时间称为弛豫时间，分为纵向弛豫时间和横向弛豫时间。
磁共振成像需要通过射频线圈来接收磁共振信号，通过外加射频脉冲，使磁化矢量发生偏转，磁化矢量在与此处垂直的横向平面内进动，使穿过线圈的磁通量发生变化，磁通量的变化在线圈中感应出电流。线圈内流动的电流大小或者说信号强度，与磁化矢量横向分量的大小正相关，与氢原子的密度以及弛豫时间有关。
采集到的磁共振信号为复数信号，得到的磁共振图像也为复数，其幅度、相位又氢原子密度、弛豫时间和扫描参数确定。在采集磁共振信号之前，需要增强新鲜样品中流体的磁共振信号。具体操作方式为：
通过辐照激发流体的电子自旋以及电子质子自旋交换，饱和流体的自旋磁矩；关闭辐照并发射射频脉冲，激发流体的磁共振信号。
图3是本申请实施例提供的流速测量系统的结构示意图，如图3所示，本申请实施例提供一种流速测量系统，包括：
增强激发模块301，用于通过辐照激发电子自旋以及电子质子自旋交换，饱和流体的自旋磁矩；关闭辐照并发射射频脉冲，激发流体到的磁共振信号。
控制采集模块302，用于所激发的流体施加外部的梯度磁场，采集原始磁共振图像。
存储模块303，用于提取原始磁共振图像的相位对比图像并进行存储。
处理模块304，用于读取存储的相位对比图像，将在(-π，π]边界跳变产生的卷绕修改为连续变化的灰度图像，获得解卷绕相位图像，并计算流体的流速大小v。
图4是本申请实施例提供的水模型流动仿体的结构示意图，如图4所示，在离心管与离心管管盖之间平行映射打孔，孔间插入软胶质导水管，导水管与离心管间填充氧化海藻酸钠和羧甲基壳聚糖在常温常压下形成水凝胶，从导水孔注入设定流速的纯水，从出水孔流出。
通过对该水模型流动仿体进行磁共振相位对比成像，得到水流速度的实际测量值和理论值的偏差均在0.2mm/s内，空间分辨率0.1mm×0.1mm，单样品的单层流速测量扫描时间为2s。
通过本申请实施例提供的流速测量系统，可以在0.5T低场核磁共振成像仪上对新鲜样品进行流速的定量测量，流速测量精度可以实现0.2mm/s，空间分辨率可以达到0.1mm×0.1mm，能够实现高空间分辨率和成像效率的新鲜样品流速定量测量，有助于病理的精准诊断，以及疾病的早期发现，对发展新的病理检测手段提供技术支撑。
图5是本申请实施例提供的流速测量装置的结构示意图，如图5所示，本申请实施例提供一种流速测量装置，包括：
获取模块501，用于在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；
确定模块502，用于基于所述相位对比图像，确定所述流体的流速大小。
可选地，所述确定模块包括：
第一获取子模块，用于基于所述相位对比图像进行解卷绕，获得解卷绕相位图像；
第一确定子模块，用于基于所述解卷绕相位图像，确定流体的流速大小。
可选地，所述获取模块包括：
第二获取子模块，用于在所述流速方向上施加第一对流速编码梯度，得到第一相位变化；
第三获取子模块，用于在所述流速方向上施加第二对流速编码梯度，得到第二相位变化；所述第一对流速编码梯度的脉冲顺序和所述第二对流速编码梯度的脉冲顺序之间呈镜像对称；
第二确定子模块，用于根据所述第一相位变化和所述第二相位变化确定所述流速方向的相位对比图像。
可选地，所述第一对流速编码梯度为大小相等、持续时间相同但方向相反的双极梯度；所述第二对流速编码梯度为大小相等，持续时间相同但方向相反的双极梯度。
可选地，所述装置还包括：
增强模块，用于增强所述新鲜样品中所述流体的磁共振信号。
可选地，所述增强模块包括：
第一激发子模块，用于通过辐照激发所述流体的电子自旋以及电子质子交换，饱和所述流体的自旋磁矩；
第二激发子模块，用于关闭辐照并发射射频脉冲，激发所述流体的磁共振信号。
在此需要说明的是，本申请实施例提供的上述装置，能够实现上述方法实施例所实现的所有方法步骤，且能够达到相同的技术效果，在此不再对本实施例中与方法实施例相同的部分及有益效果进行具体赘述。
图6是本申请实施例提供的电子设备的结构示意图，如图6所示，该电子设备可以包括：处理器(processor)601、通信接口(Communications Interface)602、存储器(memory)603和通信总线604，其中，处理器601，通信接口602，存储器603通过通信总线604完成相互间的通信。处理器601可以调用存储器603中的逻辑指令，以执行流速测量方法，该方法包括：
在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；
基于所述相位对比图像，确定所述流体的流速大小。
此外，上述的存储器603中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本申请的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本申请各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
另一方面，本申请还提供一种计算机程序产品，所述计算机程序产品包括计算机程序，计算机程序可存储在非暂态计算机可读存储介质上，所述计算机程序被处理器执行时，计算机能够执行上述各方法所提供的流速测量方法，该方法包括：
在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；
基于所述相位对比图像，确定所述流体的流速大小。
又一方面，本申请还提供一种非暂态计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现以执行上述各方法提供的流速测量方法，该方法包括：
在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；
基于所述相位对比图像，确定所述流体的流速大小。
以上所描述的装置实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到各实施方式可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的方法。
最后应说明的是：以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (10)
1.一种流速测量方法，其特征在于，包括：
在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；
基于所述相位对比图像，确定所述流体的流速大小。
2.根据权利要求1所述的流速测量方法，其特征在于，所述基于所述相位对比图像，确定所述流体的流速大小，包括：
基于所述相位对比图像进行解卷绕，获得解卷绕相位图像；
基于所述解卷绕相位图像，确定流体的流速大小。
3.根据权利要求1所述的流速测量方法，其特征在于，所述在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像，包括：
在所述流速方向上施加第一对流速编码梯度，得到第一相位变化；
在所述流速方向上施加第二对流速编码梯度，得到第二相位变化；所述第一对流速编码梯度的脉冲顺序和所述第二对流速编码梯度的脉冲顺序之间呈镜像对称；
根据所述第一相位变化和所述第二相位变化确定所述流速方向的相位对比图像。
4.根据权利要求3所述的流速测量方法，其特征在于，所述第一对流速编码梯度为大小相等、持续时间相同但方向相反的双极梯度；所述第二对流速编码梯度为大小相等，持续时间相同但方向相反的双极梯度。
5.根据权利要求1所述的流速测量方法，其特征在于，所述方法还包括：
增强所述新鲜样品中所述流体的磁共振信号。
6.根据权利要求5所述的流速测量方法，其特征在于，所述增强所述新鲜样品中所述流体的磁共振信号，包括：
通过辐照激发所述流体的电子自旋以及电子质子交换，饱和所述流体的自旋磁矩；
关闭辐照并发射射频脉冲，激发所述流体的磁共振信号。
7.一种流速测量装置，其特征在于，包括：
获取模块，用于在新鲜样品中流体的流速方向上施加流速编码梯度，获取流速方向的相位对比图像；
确定模块，用于基于所述相位对比图像，确定所述流体的流速大小。
8.一种电子设备，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至6中的任一项所述流速测量方法的步骤。
9.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6中的任一项所述流速测量方法的步骤。
10.一种计算机程序产品，包括计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6中的任一项所述流速测量方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111209652.2A CN113900056A (zh) | 2021-10-18 | 2021-10-18 | 流速测量方法、装置及存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111209652.2A CN113900056A (zh) | 2021-10-18 | 2021-10-18 | 流速测量方法、装置及存储介质 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN113900056A true CN113900056A (zh) | 2022-01-07 |
Family
ID=79192425
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111209652.2A Pending CN113900056A (zh) | 2021-10-18 | 2021-10-18 | 流速测量方法、装置及存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113900056A (zh) |
Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5034692A (en) * | 1989-02-24 | 1991-07-23 | Siemens Aktiengesellschaft | Magnetic resonance imaging method for acquiring flux-compensated, T2 -w |
| US5111145A (en) * | 1988-07-01 | 1992-05-05 | Instrumentarium Corp. | Method and apparatus for studying the properties of a material |
| WO1995017684A1 (en) * | 1993-12-22 | 1995-06-29 | Ardenkjaer Larsen Jan Henrik | Method and apparatus for measuring weak magnetic fields |
| US20110080471A1 (en) * | 2009-10-06 | 2011-04-07 | Iowa State University Research Foundation, Inc. | Hybrid method for 3D shape measurement |
| KR20110075859A (ko) * | 2009-12-29 | 2011-07-06 | (주)에이아이랩 | Ｍｒｉ 시스템의 위상오차 보정 방법 |
| CN102488497A (zh) * | 2011-12-12 | 2012-06-13 | 中国科学院深圳先进技术研究院 | 磁共振温度测量方法及系统 |
| CN103110420A (zh) * | 2013-01-31 | 2013-05-22 | 深圳先进技术研究院 | 磁共振血管成像方法和系统 |
| CN106923813A (zh) * | 2015-12-30 | 2017-07-07 | 中国科学院深圳先进技术研究院 | 一种磁共振心脏血流速度的测量方法及装置 |
| CN108426909A (zh) * | 2018-02-11 | 2018-08-21 | 中国石油大学(北京) | 核磁共振弛豫时间和流速同步快速测量的方法及装置 |
| CN110558986A (zh) * | 2019-08-27 | 2019-12-13 | 苏州朗润医疗系统有限公司 | 适用于磁共振成像的相位对比成像方法 |
| CN110992435A (zh) * | 2019-11-06 | 2020-04-10 | 上海东软医疗科技有限公司 | 图像重建方法及设备、成像数据的处理方法及装置 |
| CN111265206A (zh) * | 2020-01-22 | 2020-06-12 | 上海联影医疗科技有限公司 | 一种磁共振血管成像方法、装置及设备 |
| CN113313778A (zh) * | 2021-05-13 | 2021-08-27 | 中国科学院深圳先进技术研究院 | 磁共振图像的重建方法、计算机设备及存储介质 |
-
2021
- 2021-10-18 CN CN202111209652.2A patent/CN113900056A/zh active Pending
Patent Citations (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5111145A (en) * | 1988-07-01 | 1992-05-05 | Instrumentarium Corp. | Method and apparatus for studying the properties of a material |
| US5034692A (en) * | 1989-02-24 | 1991-07-23 | Siemens Aktiengesellschaft | Magnetic resonance imaging method for acquiring flux-compensated, T2 -w |
| WO1995017684A1 (en) * | 1993-12-22 | 1995-06-29 | Ardenkjaer Larsen Jan Henrik | Method and apparatus for measuring weak magnetic fields |
| US20110080471A1 (en) * | 2009-10-06 | 2011-04-07 | Iowa State University Research Foundation, Inc. | Hybrid method for 3D shape measurement |
| KR20110075859A (ko) * | 2009-12-29 | 2011-07-06 | (주)에이아이랩 | Ｍｒｉ 시스템의 위상오차 보정 방법 |
| CN102488497A (zh) * | 2011-12-12 | 2012-06-13 | 中国科学院深圳先进技术研究院 | 磁共振温度测量方法及系统 |
| CN103110420A (zh) * | 2013-01-31 | 2013-05-22 | 深圳先进技术研究院 | 磁共振血管成像方法和系统 |
| CN106923813A (zh) * | 2015-12-30 | 2017-07-07 | 中国科学院深圳先进技术研究院 | 一种磁共振心脏血流速度的测量方法及装置 |
| CN108426909A (zh) * | 2018-02-11 | 2018-08-21 | 中国石油大学(北京) | 核磁共振弛豫时间和流速同步快速测量的方法及装置 |
| CN110558986A (zh) * | 2019-08-27 | 2019-12-13 | 苏州朗润医疗系统有限公司 | 适用于磁共振成像的相位对比成像方法 |
| CN110992435A (zh) * | 2019-11-06 | 2020-04-10 | 上海东软医疗科技有限公司 | 图像重建方法及设备、成像数据的处理方法及装置 |
| CN111265206A (zh) * | 2020-01-22 | 2020-06-12 | 上海联影医疗科技有限公司 | 一种磁共振血管成像方法、装置及设备 |
| CN113313778A (zh) * | 2021-05-13 | 2021-08-27 | 中国科学院深圳先进技术研究院 | 磁共振图像的重建方法、计算机设备及存储介质 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Buonincontri et al. | MR fingerprinting with simultaneous B1 estimation | |
| US7230424B1 (en) | Magnetic resonance imaging | |
| CN102612657B (zh) | 针对磁共振成像的图像强度校正 | |
| US9766313B2 (en) | MR imaging using apt contrast enhancement and sampling at multiple echo times | |
| RU2523687C2 (ru) | Мр-томография, использующая параллельное получение сигнала | |
| JP4203259B2 (ja) | 核スピントモグラフィ装置および異なる化学シフトを有する２つのスピン集団の像を提供する装置 | |
| US20140070805A1 (en) | Mr imaging with b1 mapping | |
| JP5848713B2 (ja) | 磁気共鳴イメージング装置及びコントラスト強調画像取得方法 | |
| JP6356809B2 (ja) | 水／脂肪分離を有するゼロエコー時間ｍｒイメージング | |
| US10120053B2 (en) | Magnetic resonance imaging method and apparatus | |
| RU2538421C2 (ru) | Картирование градиента восприимчивости | |
| US10203387B2 (en) | MR imaging with enhanced susceptibility contrast | |
| van der Heide et al. | High‐resolution in vivo MR‐STAT using a matrix‐free and parallelized reconstruction algorithm | |
| WO2014059237A1 (en) | Systems and methods for susceptibility tensor imaging in the p-space | |
| JP3847554B2 (ja) | 磁気共鳴イメージング装置 | |
| He et al. | Phantom evaluation of electrical conductivity mapping by MRI: Comparison to vector network analyzer measurements and spatial resolution assessment | |
| JP2019536990A (ja) | 重合式フーリエ変換に基づく分光学とイメージング法 | |
| CN113900056A (zh) | 流速测量方法、装置及存储介质 | |
| KR101458557B1 (ko) | 상이한 숙임각을 갖는 자기 공명 영상 시스템에서 주자장 정보 및 라디오 펄스 관련 정보를 획득하기 위한 방법 및 장치 | |
| JP2005111059A (ja) | 磁気共鳴イメージング装置 | |
| US11815582B2 (en) | Dual echo steady state MR imaging using bipolar diffusion gradients | |
| US11733334B2 (en) | Image quality assessment of magnetic resonance images using a phantom | |
| KR101475932B1 (ko) | 하이브리드 자기 공명 영상 처리 장치 및 방법 | |
| WO2012030548A1 (en) | Traveling wave mri in high bo field strengths | |
| Smith et al. | Design and use of variable flip angle schedules in transient balanced SSFP subtractive imaging |
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 |