CN110101379B - 流速测量方法、装置、磁共振成像设备 - Google Patents
流速测量方法、装置、磁共振成像设备 Download PDFInfo
- Publication number
- CN110101379B CN110101379B CN201910329011.7A CN201910329011A CN110101379B CN 110101379 B CN110101379 B CN 110101379B CN 201910329011 A CN201910329011 A CN 201910329011A CN 110101379 B CN110101379 B CN 110101379B
- Authority
- CN
- China
- Prior art keywords
- flow rate
- flow velocity
- pulse sequence
- target
- included angle
- 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
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000002595 magnetic resonance imaging Methods 0.000 title claims abstract description 8
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 claims abstract description 115
- 239000007788 liquid Substances 0.000 claims abstract description 77
- 230000005284 excitation Effects 0.000 claims abstract description 55
- 238000005259 measurement Methods 0.000 claims abstract description 24
- 238000013507 mapping Methods 0.000 claims description 63
- 238000000691 measurement method Methods 0.000 claims description 21
- 238000004590 computer program Methods 0.000 claims description 17
- 238000001514 detection method Methods 0.000 claims 1
- 230000000747 cardiac effect Effects 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 7
- 238000011002 quantification Methods 0.000 description 5
- 230000017531 blood circulation Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 108091026890 Coding region Proteins 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012806 monitoring device Methods 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording for evaluating the cardiovascular system, e.g. pulse, heart rate, blood pressure or blood flow
- A61B5/026—Measuring blood flow
- A61B5/0265—Measuring blood flow using electromagnetic means, e.g. electromagnetic flowmeter
-
- 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
-
- 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
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Molecular Biology (AREA)
- Animal Behavior & Ethology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- High Energy & Nuclear Physics (AREA)
- Surgery (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Radiology & Medical Imaging (AREA)
- Electromagnetism (AREA)
- Hematology (AREA)
- Cardiology (AREA)
- Physiology (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本申请涉及一种流速测量方法、装置、磁共振成像设备，计算机设备先确定包括流动液体的目标区域，并采用第一脉冲序列在预设层面上激发该目标区域，得到流动液体的第一流速，以及采用第二脉冲序列在与预设层面存在预设夹角的方向上激发该目标区域，得到流动液体的第二流速，然后根据该第一流速和第二流速确定流动液体的目标流速，该方法中，第二脉冲序列的激发层面方向与第一脉冲序列的激发层面方向不同，且该第二脉冲序列的激发层面方向是预设的，使得获取的第一流速、第二流速、第一流速和第二流速之间的夹角均为已知值，这样通过多个已知量就可以确定出目标流速，以该目标流速作为流速值测量结果，大大提高了流速值测量结果的精确性。
Description
技术领域
本申请涉及磁共振技术领域，特别是涉及一种流速测量方法、装置、磁共振成像设备。
背景技术
流动定量(flow quantification，FQ)表示基于gre_fq序列，利用ECG分段触发扫描，在一个心动周期内可显示十余甚至数十幅不同心脏时相的图像，从而获得一个心动周期内的血流速度变化。
在进行2D流动定量进行流速测量时，扫描层面方向S必须垂直于血流方向，才能得到正确的流速值V，但在实际应用中，扫描层面往往会因为多种因素导致实际扫描层面的方向S1与理想扫描层面方向S存在一定的误差角度θ，使得实际的流速值测量结果为V1＝V*cos(θ)。
因此，流动定量进行流速测量的实际应用中，存在流速值测量结果不精确的技术问题。
发明内容
基于此，有必要针对上述流动定量进行流速测量的实际应用中，存在流速值测量结果不精确的技术问题，提供一种流速测量方法、装置、磁共振成像设备。
第一方面，本申请实施例提供一种流速测量方法，该方法包括：
获取目标区域；其中，目标区域内包含流动液体；
在预设层面采用第一脉冲序列激发目标区域，获取流动液体的第一流速；其中，预设层面方向与流动液体的流动方向成垂直关系；
在与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，获取流动液体的第二流速；其中，第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向不同；
根据第一流速和第二流速确定流动液体的目标流速。
在其中一个实施例中，第一脉冲序列的激发层面方向与预设层面方向之间存在任意夹角。
在其中一个实施例中，上述根据第一流速和第二流速获取流动液体的目标流速，包括：
获取第一流速与目标流速的第一映射关系、第二流速与目标流速之间的第二映射关系；
根据第一映射关系和第二映射关系，确定流动液体的目标流速。
在其中一个实施例中，上述获取第一流速与目标流速的第一映射关系、第二流速与目标流速之间的第二映射关系，包括：
根据第一夹角的角度值，确定第一流速与目标流速的第一映射关系；其中，第一夹角表示第一脉冲序列的激发层面方向与预设层面方向之间的夹角；
根据第一夹角的角度值和第二夹角角度值，确定第二流速与目标流速之间的第二映射关系；其中，第二夹角表示第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向之间的夹角。
在其中一个实施例中，上述第一映射关系包括：第一流速等于目标流速与第一夹角角度值的余弦值的乘积。
在其中一个实施例中，上述第二映射关系包括：第二流速等于目标流速与第一夹角角度值和第二夹角角度值之和的余弦值的乘积。
在其中一个实施例中，上述在预设层面采用第一脉冲序列激发目标区域，获取流动液体的第一流速，包括：
分别从X、Y、Z方向使用第一脉冲序列在预设层面激发所述目标区域，得到第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量；
根据第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量确定第一流速。
在其中一个实施例中，上述与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，获取流动液体的第二流速，包括：
分别从X、Y、Z方向与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，得到第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量；
根据第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量确定第二流速。
第二方面，本申请实施例提供一种流速测量装置，所述装置包括：
目标区域获取模块，用于获取目标区域；其中，目标区域内包含流动液体；
第一流速获取模块，用于在预设层面采用第一脉冲序列激发目标区域，获取流动液体的第一流速；其中，预设层面方向与流动液体的流动方向成垂直关系；
第二流速获取模块，用于在与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，获取流动液体的第二流速；其中，第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向不同；
目标流速确定模块，用于根据第一流速和第二流速确定流动液体的目标流速。
第三方面，本申请实施例提供一种磁共振成像设备，包括存储器和处理器以及扫描设备，存储器存储有计算机程序，处理器控制上述扫描设备执行计算机程序时实现上述第一方面实施例提供的任一项方法的步骤。
第四方面，本申请实施例提供一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述第一方面实施例提供的任一项方法的步骤。
本申请实施例提供的一种流速测量方法、装置、磁共振成像设备，计算机设备先确定包括流动液体的目标区域，并采用第一脉冲序列在预设层面上激发该目标区域，得到流动液体的第一流速，以及采用第二脉冲序列在与预设层面存在预设夹角的方向上激发该目标区域，得到流动液体的第二流速，然后根据该第一流速和第二流速确定流动液体的目标流速，由于该方法中，第二脉冲序列的激发层面方向与第一脉冲序列的激发层面方向不同，且该第二脉冲序列的激发层面方向是预设的，使得获取的第一流速、第二流速、第一流速和第二流速之间的夹角均为已知值，这样通过多个已知量就可以确定出目标流速，以该目标流速作为流速值测量结果，大大提高了流速值测量结果的精确性。
附图说明
图1为一个实施例提供的一种流速测量方法的应用环境图；
图2为一个实施例提供的一种流速测量方法的流程示意图；
图2a为一个实施例提供的一种流速测量参数方向示意图；
图3为一个实施例提供的一种流速测量方法的流程示意图；
图4为一个实施例提供的一种流速测量方法的流程示意图；
图4a为一个实施例所示的所用的gre_fq序列的示意图；
图5为一个实施例提供的一种流速测量方法的流程示意图；
图6为一个实施例提供的一种流速测量方法的流程示意图；
图7为一个实施例提供的一种流速测量装置的结构框图；
图8为一个实施例提供的一种流速测量装置的结构框图；
图9为一个实施例提供的一种流速测量装置的结构框图；
图10为一个实施例提供的一种流速测量装置的结构框图；
图11为一个实施例提供的一种流速测量装置的结构框图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。
本申请提供的一种流速测量方法，可以应用于如图1所示的应用环境中，该计算机设备包括通过系统总线连接的处理器、存储器、网络接口和数据库。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统、计算机程序和数据库。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的数据库用于存储流速测量数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种流速测量方法。
本申请的实施例提供一种流速测量方法、装置、计算机设备和存储介质，旨在解决流动定量进行流速测量的实际应用中，存在流速值测量结果不精确的技术问题。下面将通过实施例并结合附图具体地对本申请的技术方案以及本申请的技术方案如何解决上述技术问题进行详细说明。下面这几个具体的实施例可以相互结合，对于相同或相似的概念或过程可能在某些实施例中不再赘述。需要说明的是，本申请提供的一种流速测量方法，图2-图6的执行主体为计算机设备，其中，其执行主体还可以是流速测量装置，其中该装置可以通过软件、硬件或者软硬件结合的方式实现成为流速测量的部分或者全部。
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。
在一个实施例中，图2提供了一种流速测量方法，本实施例涉及的是计算机设备先获取目标区域，以及目标区域中流动液体的第一流速和第二流速，然后根据该第一流速和第二流速确定目标流速的具体过程，如图2所示，所述方法包括：
S101，获取目标区域；其中，目标区域内包含流动液体。
其中，目标区域表示的是待进行磁共振扫描的区域，即待进行脉冲序列激发的区域，其中该目标区域中包含流动液体，其中该流动液体可以是血液、脑脊液等，本实施例对此不做限定。则在实际应用中，计算机设备获取目标区域的方式可以是根据用户输入的区域标识确定该目标区域，或者是根据待扫描要求分析结果确定该目标区域等，还可以是其他方式，本实施例对此不做限定。
S102，在预设层面采用第一脉冲序列激发目标区域，获取流动液体的第一流速；其中，预设层面方向与流动液体的流动方向成垂直关系。
本步骤中，基于上述S101步骤中获取的目标区域，计算机设备在预设层面采用第一脉冲序列激发目标区域，获取流动液体的第一流速，其中，预设层面方向与流动液体的流动方向成垂直关系。其中，预设层面为预先设定的目标区域中的一个层面，则在实际应用中，计算机设备采用第一脉冲序列激发在该预设层面上激发该目标区域，对应得到该流动液体的第一流速，该第一脉冲序列例如可以是gre_fq等，本实施例对此不做限定。可选地，第一脉冲序列的激发层面方向与预设层面之间存在任意夹角。其中，任意夹角表示的是该第一脉冲序列的激发层面方向与预设层面之间的夹角为未知的，且其具体角度值为任意值，也就是说该第一脉冲序列的方向与预设层面之间存在夹角是由于客观因素引起的，不是设定的。需要说明的是，如图2a所示，预设层面的方向与流动液体的流动方向成垂直关系，相应地，可以得到第一脉冲序列的激发层面方向与第一流速方向成垂直关系，以及后续涉及的第二脉冲序列的激发层面方向与第二流速方向也成垂直关系，具体可参见图2a所示。
S103，在与所述预设层面存在预设夹角的方向采用第二脉冲序列激发所述目标区域，获取流动液体的第二流速；其中，第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向不同。
本步骤中，基于上述S102步骤中获取的流动液体的第一流速，计算机设备在与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，获取流动液体的第二流速，其中，该第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向不同，其中，第二脉冲序列的激发层面方向为根据预设层面方向设定的已知方向，利用该方法可以使第一脉冲序列的激发层面方向和第二脉冲序列的激发层面方向之间的夹角为已知值。则在实际应用中，计算机设备在与所述预设层面存在预设夹角的方向采用第二脉冲序列激发所述目标区域，对应得到该流动液体的第二流速。基于上述对第一脉冲序列层面与第二脉冲序列层面之间的夹角Φ是已知的关系，再结合图2a所示，可以得到第一流速V1和第二流速V2之间的夹角也为已知值Φ。
S104，根据第一流速和第二流速确定流动液体的目标流速。
基于上述S103步骤中的第一流速和S104步骤中的第二流速，计算机设备确定流动液体的目标流速，示例地，计算机设备根据根据第一流速和第二流速确定流动液体的目标流速的方式可以是将该第一流速和第二流速作为输入数据，输入进预先训练好的神经网络模型中得到该目标流速，也可以是采用预设算法根据该第一流速和第二流速之间的已知夹角角度值确定该目标流速等，当然还可以是其他方式，本实施例对此不做限定。
本实施例提供的流速测量方法，计算机设备先确定包括流动液体的目标区域，并采用第一脉冲序列在预设层面上激发该目标区域，得到流动液体的第一流速，以及采用第二脉冲序列在与预设层面存在预设夹角的方向上激发该目标区域，得到流动液体的第二流速，然后根据该第一流速和第二流速确定流动液体的目标流速，由于该方法中，第二脉冲序列的激发层面方向与第一脉冲序列的激发层面方向不同，且该第二脉冲序列的激发层面方向是预设的，使得获取的第一流速、第二流速、第一流速和第二流速之间的夹角均为已知值，这样通过多个已知量就可以确定出目标流速，以该目标流速作为流速值测量结果，大大提高了流速值测量结果的精确性。
对于计算机设备根据第一流速和第二流速确定目标流速的具体过程，本申请实施例将通过以下几个实施例进行详细说明，则在以上实施例的基础上，本申请实施例还提供了一种流速测量方法，其涉及的是计算机设备根据第一流速与目标流速之间的第一映射关系和第二流速与目标流速之间的第二映射关系，确定流动液体的目标流速的具体过程，如图3所示，上述S104步骤包括：
S201，获取第一流速与目标流速的第一映射关系、第二流速与目标流速之间的第二映射关系。
本实施例中，计算机设备获取第一流速和目标流速的第一映射关系、第二流速与目标流速之间的第二映射关系，其获取方式可以是根据第一流速方向和目标流速方向之间的关系确定该第一映射关系，根据第二流速方向与目标流速方向之间的关系确定该第二映射关系，本实施例对此不做限定。下面提供一种计算机设备获取第一映射关系和第二映射关系的可实现方式。
可选地，如图4所示，该S201步骤的一种可实现方式包括：
S301，根据第一夹角的角度值，确定第一流速与目标流速的第一映射关系；其中，第一夹角表示第一脉冲序列的激发层面方向与预设层面方向之间的夹角。
本步骤中，计算机设备根据第一夹角的角度值确定第一映射关系，其中，该第一夹角表示第一脉冲序列的激发层面方向与预设层面方向之间的夹角，可选地，第一映射关系包括：第一流速等于目标流速与第一夹角角度值的余弦值的乘积。
示例地，以上述图2a为例，设第一夹角为θ，第一流速为V1,目标流速为V,则该第一映射关系可以表示为V1＝Vcosθ。
另外本申请实施例提供了一种通过激发gre_fq序列来确定流速的具体方式，具体地，利用心电监控设备ECG分段触发扫描，然后在一个心动周期内可显示十余甚至数十幅不同心脏时相的图像，这样就可以获得一个心动周期内的血流速度变化。如图4a所示为本申请实施例所用的gre_fq序列示意图，该图中是利用ECG获取受检者的心脏运动曲线，在每次心脏运动的R波触发流动补偿和流动编码。其中，图中的SS表示层选梯度编码；RF表示射频脉冲；PE表示相位梯度编码；RO表示频率读出梯度编码。本申请实施例中，在每次心脏运动周期的R波触发由三个梯度方向的梯度脉冲组成的流动补偿和流动编码序列。需要说明的是gre_fq序列的原理是利用流动编码对流动信号产生一额外的相位，以获得流速信息。示例性的，利用速度(流动)编码梯度将横向磁化矢量的速度编码到磁共振信号的相位信息中。可得相位与流速关系如下三个公式所示，其中，下述公式中γ表示磁旋比，G表示梯度，x表示位置，υ表示运动速度，φ表示相位。
S302，根据第一夹角的角度值和第二夹角角度值，确定第二流速与目标流速之间的第二映射关系；其中，第二夹角表示第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向之间的夹角。
本步骤中，计算机设备根据第一夹角的角度值和第二夹角角度值确定第二映射关系，其中，第二夹角表示第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向之间的夹角，根据上述实施例可知该第二脉冲序列的激发层面方向与第一脉冲序列的激发层面方向之间第二夹角角度值为已知值，基于上述实施例中脉冲激发层面与液体流速之间的关系，则可选地，第二映射关系包括：第二流速等于目标流速与第一夹角角度值和预设的第二夹角角度值之和的余弦值的乘积。
示例地，继续以上述图2a为例，设第一夹角为θ，第二夹角为Φ，第一流速为V1,第二流速为V2,目标流速为V,则该第一映射关系可以表示为V2＝V cos(θ+Φ)。
S202，根据第一映射关系和第二映射关系，确定流动液体的目标流速。
基于上述S201步骤中的第一映射关系和第二映射关系，计算机设备确定流动液体的目标流速，示例地，以上述第一映射关系为V1＝Vcosθ，第二映射关系为V2＝V cos(θ+Φ)为例，则计算机设备确定的流动液体的目标流速V可以表示为
该表达式中，V1、V2以及Φ均为已知值，这样就可以根据该表达式确定出目标流速v的值。
本实施例提供的流速测量方法，由于计算机设备先确定了第一流速与目标流速之间的第一映射关系，确定了第二流速与目标流速之间的第二映射关系，然后根据该第一映射关系和第二映射关系确定出目标流速的表达式，由于目标流速的表达式中除了目标流速，其余值均是已知值，就可以准确的确定出目标流速，这样以该目标流速作为流速值测量结果，大大提高了流速值测量结果的精确性。
考虑到本申请实施例提供的流速测量方法在3D环境中应用时，计算机设备获取第一流速和第二流速，可以是通过三维空间的X、Y、Z方向进行合成，则如图5所示，本申请实施例提供了一种流速测量方法，其涉及的是计算机设备从X、Y、Z方向合成第一流速的具体过程，上述S102步骤包括：
S401，分别从X、Y、Z方向使用第一脉冲序列在预设层面激发所述目标区域，得到第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量。
本实施例中，计算机设备从X、Y、Z方向使用第一脉冲序列在预设层面激发所述目标区域，相应的得到第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量。
S402，根据第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量确定第一流速。
基于上述S401步骤中得到的第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量，计算机设备确定第一流速，示例地，设第一流速X方向分量是V1X、第一流速Y方向分量V1Y和第一流速Z方向分量V1Z,则
同样，在另外一个实施例中，如图6所示，本申请实施例提供了一种流速测量方法，其涉及的是计算机设备从X、Y、Z方向合成第二流速的具体过程，上述S103步骤包括：
S501，分别从X、Y、Z方向在与所述预设层面存在预设夹角的方向采用第二脉冲序列激发所述目标区域，得到第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量。
本实施例中，计算机设备从X、Y、Z方向使用第二脉冲序列在预设层面激发所述目标区域，相应的得到第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量。
S502，根据第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量确定第二流速。
基于上述S501步骤中得到的第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量，计算机设备确定第二流速，示例地，设第二流速X方向分量是V2X、第一流速Y方向分量V2Y和第一流速Z方向分量V2Z,则
本实施例提供的流速测量方法，计算机设备先确定出第一流速或者第二流速的XYZ三个方向上分量，然后将该三个分量合成为第一流速或者第二流速，这样，通过分量合成流速，大大提高第一流速或者第二流速的精确度。
应该理解的是，虽然图2-6的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2-6中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
在一个实施例中，如图7所示，提供了一种流速测量装置，该装置包括：目标区域获取模块10、第一流速获取模块11、第二流速获取模块12和目标流速确定模块13，其中，
目标区域获取模块10，用于获取目标区域；其中，目标区域内包含流动液体；
第一流速获取模块11，用于在预设层面采用第一脉冲序列激发目标区域，获取流动液体的第一流速；其中，预设层面方向与流动液体的流动方向成垂直关系；
第二流速获取模块12，用于在与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，获取流动液体的第二流速；其中，第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向不同；
目标流速确定模块13，用于根据第一流速和第二流速确定流动液体的目标流速。
上述实施例提供的一种流速测量装置，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在其中一个实施例中，第一脉冲序列的激发层面方向与预设层面方向之间存在任意夹角。
在一个实施例中，如图8所示，提供了一种流速测量装置，上述目标流速确定模块13包括：映射关系获取单元131和目标流速确定单元132，其中，
映射关系获取单元131，用于获取第一流速与目标流速的第一映射关系、第二流速与目标流速之间的第二映射关系；
目标流速确定单元132，用于根据第一映射关系和第二映射关系，确定流动液体的目标流速。
上述实施例提供的一种流速测量装置，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，如图9所示，提供了一种流速测量装置，上述映射关系获取单元131包括：第一映射关系确定子单元1311和第一映射关系确定子单元1312，其中，
第一映射关系确定子单元1311，用于根据第一夹角的角度值，确定第一流速与目标流速的第一映射关系；其中，第一夹角表示第一脉冲序列的激发层面方向与预设层面方向之间的夹角；
第一映射关系确定子单元1312，用于根据第一夹角的角度值和第二夹角角度值，确定第二流速与目标流速之间的第二映射关系；其中，第二夹角表示第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向之间的夹角。
上述实施例提供的一种流速测量装置，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，上述第一映射关系包括：第一流速等于目标流速与第一夹角角度值的余弦值的乘积。
在一个实施例中，上述第二映射关系包括：第二流速等于目标流速与第一夹角角度值和第二夹角角度值之和的余弦值的乘积。
在一个实施例中，如图10所示，提供了一种流速测量装置，上述第一流速获取模块11，包括：第一分量确定单元111和第一流速确定单元112，其中，
第一分量确定单元111，用于分别从X、Y、Z方向使用第一脉冲序列在预设层面激发所述目标区域，得到第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量；
第一流速确定单元112，用于根据第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量确定第一流速。
上述实施例提供的一种流速测量装置，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，如图11所示，提供了一种流速测量装置，上述第二流速获取模块12，包括：第二分量确定单元121和第二流速确定单元122，其中，
第二分量确定单元121，用于分别从X、Y、Z方向与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，得到第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量；
第二流速确定单元122，用于根据第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量确定第二流速。
上述实施例提供的一种流速测量装置，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
关于流速测量装置的具体限定可以参见上文中对于流速测量方法的限定，在此不再赘述。上述流速测量装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如上述图1所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种流速测量方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，上述图1中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器和扫描设备，存储器中存储有计算机程序，该处理器控制该扫描设备执行计算机程序时实现以下步骤：
获取目标区域；其中，目标区域内包含流动液体；
在预设层面采用第一脉冲序列激发目标区域，获取流动液体的第一流速；其中，预设层面方向与流动液体的流动方向成垂直关系；
在与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，获取流动液体的第二流速；其中，第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向不同；
根据第一流速和第二流速确定流动液体的目标流速。
其中，扫描设备包括磁体、梯度线圈、射频发射线圈、射频接受线圈等，则示例性地，在预设层面采用第一脉冲序列激发目标区域具体为：处理器控制扫描设备中的梯度线圈在检测对象心脏周期的R波发射第一脉冲序列，该第一脉冲序列包括流动补偿和流动编码梯度；示例地，在预设层面采用第二脉冲序列激发目标区域具体为：处理器控制扫描设备中的梯度线圈在检测对象心脏周期的R波发射第二脉冲序列，该第二脉冲序列包括流动补偿和流动编码梯度，第一脉冲序列相关的施加梯度脉冲序列与第二脉冲序列相关的施加梯度脉冲序列强度相同，不同之处在于施加的方向。更具体地，第一脉冲序列的RF脉冲对应的选层梯度和第二脉冲序列的RF脉冲对应的选层梯度不同。
上述实施例提供的一种计算机设备，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
获取目标区域；其中，目标区域内包含流动液体；
在预设层面采用第一脉冲序列激发目标区域，获取流动液体的第一流速；其中，预设层面方向与流动液体的流动方向成垂直关系；
在与预设层面存在预设夹角的方向采用第二脉冲序列激发目标区域，获取流动液体的第二流速；其中，第一脉冲序列的激发层面方向与第二脉冲序列的激发层面方向不同；
根据第一流速和第二流速确定流动液体的目标流速。
上述实施例提供的一种计算机可读存储介质，其实现原理和技术效果与上述方法实施例类似，在此不再赘述。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种流速测量方法，其特征在于，所述方法包括：
获取目标区域，所述目标区域内包含流动液体；
在预设层面采用第一脉冲序列激发所述目标区域，获取所述流动液体的第一流速；所述预设层面方向与所述流动液体的流动方向成垂直关系；所述第一脉冲序列的激发层面方向与所述预设层面方向之间存在任意角度值的第一夹角；
在与所述预设层面存在预设夹角的方向采用第二脉冲序列激发所述目标区域，获取所述流动液体的第二流速；所述第一脉冲序列的激发层面方向与所述第二脉冲序列的激发层面方向不同；所述第一脉冲序列的激发层面方向与所述第二脉冲序列的激发层面方向之间的夹角为第二夹角；
根据所述第一夹角的角度值，确定所述第一流速与所述流动液体的目标流速之间的第一映射关系，根据所述第二夹角的角度值，确定所述第二流速与所述目标流速之间的第二映射关系，并根据所述第一映射关系和所述第二映射关系确定所述流动液体的目标流速。
2.根据权利要求1所述的方法，其特征在于，所述第一映射关系包括：所述第一流速等于所述目标流速与所述第一夹角角度值的余弦值的乘积。
3.根据权利要求2所述的方法，其特征在于，所述第二映射关系包括：所述第二流速等于所述目标流速与所述第一夹角角度值和所述预设的第二夹角角度值之和的余弦值的乘积。
4.根据权利要求1所述的方法，其特征在于，所述在预设层面采用第一脉冲序列激发所述目标区域，获取所述流动液体的第一流速，包括：
分别从X、Y、Z方向使用所述第一脉冲序列在所述预设层面激发所述目标区域，得到第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量；
根据所述第一流速X方向分量、所述第一流速Y方向分量和所述第一流速Z方向分量确定所述第一流速。
5.根据权利要求1所述的方法，其特征在于，所述在与所述预设层面存在预设夹角的方向采用第二脉冲序列激发所述目标区域，获取所述流动液体的第二流速，包括：
分别从X、Y、Z方向在与所述预设层面存在预设夹角的方向采用第二脉冲序列激发所述目标区域，得到第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量；
根据所述第二流速X方向分量、所述第二流速Y方向分量和所述第二流速Z方向分量确定所述第二流速。
6.一种流速测量装置，其特征在于，所述装置包括：
目标区域获取模块，用于获取目标区域;所述目标区域内包含流动液体；
第一流速获取模块，用于在预设层面采用第一脉冲序列激发所述目标区域，获取所述流动液体的第一流速；所述预设层面方向与所述流动液体的流动方向成垂直关系；所述第一脉冲序列的激发层面方向与所述预设层面方向之间存在任意角度值的第一夹角；
第二流速获取模块，用于在与所述预设层面存在预设夹角的方向采用第二脉冲序列激发所述目标区域，获取所述流动液体的第二流速；所述第一脉冲序列的激发层面方向与所述第二脉冲序列的激发层面方向不同；所述第一脉冲序列的激发层面方向与所述第二脉冲序列的激发层面方向之间的夹角为第二夹角；
目标流速确定模块，所述目标流速确定模块包括：
第一映射关系确定子单元，用于根据所述第一夹角的角度值，确定所述第一流速与所述流动液体的目标流速之间的第一映射关系；
第二映射关系确定子单元，用于根据所述第二夹角的角度值，确定所述第二流速与所述目标流速之间的第二映射关系；
目标流速确定单元，用于根据所述第一映射关系和所述第二映射关系确定所述流动液体的目标流速。
7.根据权利要求6所述装置，其特征在于，所述第一流速获取模块包括：第一分量确定单元和第一流速确定单元，其中，
所述第一分量确定单元，用于分别从X、Y、Z方向使用所述第一脉冲序列在所述预设层面激发所述目标区域，得到第一流速X方向分量、第一流速Y方向分量和第一流速Z方向分量；
所述第一流速确定单元，用于根据所述第一流速X方向分量、所述第一流速Y方向分量和所述第一流速Z方向分量确定所述第一流速。
8.根据权利要求6所述装置，其特征在于，所述第二流速获取模块包括：第二分量确定单元和第二流速确定单元，其中，
所述第二分量确定单元，用于分别从X、Y、Z方向与预设层面存在预设夹角的方向采用所述第二脉冲序列激发所述目标区域，得到第二流速X方向分量、第二流速Y方向分量和第二流速Z方向分量；
所述第二流速确定单元，用于根据所述第二流速X方向分量、所述第二流速Y方向分量和所述第二流速Z方向分量确定所述第二流速。
9.一种磁共振成像设备，包括存储器和处理器以及扫描设备，所述存储器存储有计算机程序，其特征在于，所述处理器控制所述扫描设备执行所述计算机程序时实现权利要求1至5中任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行所述计算机程序时实现权利要求1至5中任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910329011.7A CN110101379B (zh) | 2019-04-23 | 2019-04-23 | 流速测量方法、装置、磁共振成像设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910329011.7A CN110101379B (zh) | 2019-04-23 | 2019-04-23 | 流速测量方法、装置、磁共振成像设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110101379A CN110101379A (zh) | 2019-08-09 |
| CN110101379B true CN110101379B (zh) | 2022-03-01 |
Family
ID=67486313
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910329011.7A Active CN110101379B (zh) | 2019-04-23 | 2019-04-23 | 流速测量方法、装置、磁共振成像设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110101379B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5891032A (en) * | 1997-04-10 | 1999-04-06 | Elscint Ltd | Fat free TOF angiography |
| CN107536609A (zh) * | 2016-06-28 | 2018-01-05 | 上海联影医疗科技有限公司 | 磁共振成像设备及磁共振成像的扫描方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9700220B2 (en) * | 2006-04-25 | 2017-07-11 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
| JP5665040B2 (ja) * | 2009-09-10 | 2015-02-04 | 学校法人上智学院 | 変位計測方法及び装置、並びに、超音波診断装置 |
| DE102011077197B4 (de) * | 2011-06-08 | 2013-05-16 | Siemens Aktiengesellschaft | Verzeichnungskorrektur bei einer Magnetresonanz-Bildgebung |
| US9075121B2 (en) * | 2011-07-15 | 2015-07-07 | Wisconsin Alumni Research Foundation | System and method for rotating angle velocity encoding, phase contrast magnetic resonance imaging |
-
2019
- 2019-04-23 CN CN201910329011.7A patent/CN110101379B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5891032A (en) * | 1997-04-10 | 1999-04-06 | Elscint Ltd | Fat free TOF angiography |
| CN107536609A (zh) * | 2016-06-28 | 2018-01-05 | 上海联影医疗科技有限公司 | 磁共振成像设备及磁共振成像的扫描方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110101379A (zh) | 2019-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8191359B2 (en) | Motion estimation using hidden markov model processing in MRI and other applications | |
| JP5198859B2 (ja) | 位相マッピングと、位相基準として用いる基準媒体が関係するｍｒｉ温度測定 | |
| CN103800010B (zh) | 用利用多回波磁振成像的多步自适应拟合法的脂和铁量化 | |
| CN105467340B (zh) | 具有多个子系统的医学成像检查设备的运行 | |
| US20120262172A1 (en) | Motion correction in accelerated t1-weighted magnetic resonance imaging | |
| JP6320357B2 (ja) | 複数のサブシステムを含む医療用撮像検査デバイスの動作 | |
| US10338183B2 (en) | Magnetic resonance imaging with different frequency encoding patterns | |
| CN103885016B (zh) | 相位差图的确定 | |
| US8290566B2 (en) | Magnetic resonance imaging apparatus and image generating method | |
| CN110101379B (zh) | 流速测量方法、装置、磁共振成像设备 | |
| CN109620228B (zh) | 磁共振成像中脂肪零点偏差校正方法及磁共振成像方法 | |
| CN113495242B (zh) | 相位误差探测方法、装置、磁共振系统及其成像方法 | |
| US10054653B2 (en) | Magnetic resonance method and apparatus for quantitative simultaneous multi-slice assessment of tissue displacement, deformation, and related biomarker parameters | |
| US9588209B2 (en) | Method of multislice MR elastography with multiband acquisition | |
| CN111009020A (zh) | 图像重建方法、装置、计算机设备和存储介质 | |
| CN110542872A (zh) | 磁共振成像方法及设备、相位信息的获取方法及装置 | |
| JP2006304818A (ja) | 画像合成方法及び磁気共鳴イメージング装置 | |
| CN114675221A (zh) | 磁共振梯度校正补偿因子的确定方法、校正方法和装置 | |
| CN112652049B (zh) | 心脏的扫描方法、装置、计算机设备和存储介质 | |
| CN115222640A (zh) | 磁共振图像的校正方法、装置、计算机设备和存储介质 | |
| CN115439335A (zh) | 图像校正方法、装置、计算机设备和存储介质 | |
| CN114463451B (zh) | 图像重建方法、装置、计算机设备和存储介质 | |
| CN114089242B (zh) | 激发脉冲的设计方法、装置、计算机设备和存储介质 | |
| CN113625209B (zh) | 磁共振系统频率漂移量的确定方法、装置和计算机设备 | |
| CN117414124B (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 | ||
| 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 |