CN106960458B - 一种磁共振磁敏感加权成像后处理方法及系统 - Google Patents
一种磁共振磁敏感加权成像后处理方法及系统 Download PDFInfo
- Publication number
- CN106960458B CN106960458B CN201710149852.0A CN201710149852A CN106960458B CN 106960458 B CN106960458 B CN 106960458B CN 201710149852 A CN201710149852 A CN 201710149852A CN 106960458 B CN106960458 B CN 106960458B
- Authority
- CN
- China
- Prior art keywords
- phase
- mask
- image
- local field
- amplitude
- 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
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 87
- 238000012805 post-processing Methods 0.000 title claims abstract description 48
- 230000035945 sensitivity Effects 0.000 title claims abstract description 46
- 238000003384 imaging method Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000008859 change Effects 0.000 claims abstract description 57
- 238000006243 chemical reaction Methods 0.000 claims abstract description 40
- 238000001914 filtration Methods 0.000 claims abstract description 32
- 238000000605 extraction Methods 0.000 claims description 3
- 210000003462 vein Anatomy 0.000 abstract description 6
- 238000012545 processing Methods 0.000 description 7
- 239000000126 substance Substances 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 239000008280 blood Substances 0.000 description 4
- 210000004369 blood Anatomy 0.000 description 4
- 210000004556 brain Anatomy 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- INGWEZCOABYORO-UHFFFAOYSA-N 2-(furan-2-yl)-7-methyl-1h-1,8-naphthyridin-4-one Chemical compound N=1C2=NC(C)=CC=C2C(O)=CC=1C1=CC=CO1 INGWEZCOABYORO-UHFFFAOYSA-N 0.000 description 2
- 108010002255 deoxyhemoglobin Proteins 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000005298 paramagnetic effect Effects 0.000 description 2
- 102000008857 Ferritin Human genes 0.000 description 1
- 238000008416 Ferritin Methods 0.000 description 1
- 108050000784 Ferritin Proteins 0.000 description 1
- 108010061951 Methemoglobin Proteins 0.000 description 1
- 206010028980 Neoplasm Diseases 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 230000036244 malformation Effects 0.000 description 1
- 206010061311 nervous system neoplasm Diseases 0.000 description 1
- 208000015122 neurodegenerative disease Diseases 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000008733 trauma Effects 0.000 description 1
- 230000002792 vascular Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration by the use of more than one image, e.g. averaging, subtraction
-
- G06T5/70—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10088—Magnetic resonance imaging [MRI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30096—Tumor; Lesion
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
Abstract
本发明公开了一种磁共振磁敏感加权成像后处理方法及系统，所述方法利用局部场变化分布生成局部场变化蒙板；再利用低通滤波器对原始K空间数据进行滤波将滤波前后的图像相除提取图像数据，按照正相位的值域转换准则生成相位蒙板；最后将两者按照预设准则结合生成最终蒙板，再结合幅值图像经过N次幂加权得到磁敏感加权图像。本发明通过采用局部场变化分布生成局部场变化蒙板结合图像的相位信息对所述磁共振图像进行后处理，减少了图像内组织交界处的伪影，从而避免了伪影受所选窗宽影响而造成丢失细小静脉信号的问题，提高了磁共振磁敏感加权图像的细节显示。
Description
技术领域
本发明涉及磁敏感加权成像技术领域，特别涉及一种磁共振磁敏感加权成像后处理方法及系统。
背景技术
磁敏感加权成像是由Haccke提出的一种新的磁共振成像方法。该方法利用不同组织间磁敏感性的差异产生图像对比。它能够清晰地显示静脉组织，故最初被称为”高分辨率BOLD静脉成像“（High Resolution BOLD Venography,HRBV)。然而，磁敏感加权成像不仅能够显示静脉，对去氧血红蛋白、含铁血黄素、铁蛋白和铁等顺磁性物质拥有较高的敏感性；除此之外，它还能显示肿瘤的微小血管结构、外伤和卒中后常规磁共振成像不易发现的脑改变。因此，磁敏感加权成像逐渐应用于神经系统肿瘤、脑血管畸形以及脑变性疾病的临床诊断，且具有重要意义的医学检查手段。
目前，磁敏感加权成像技术包括高分辨率的三维梯度回波成像（3D GradientEcho，GRE3D)、在三个方向上完全流动补偿、基于相位蒙板的图像后处理、邻近层面的最小密度投影。梯度回波序列采用短TR长TE的扫描参数获得重T2*加权（Heavy T2* Weight)效应图像。静脉血、去氧血红蛋白、铁等顺磁性物质在图像上能区别于其它组织表现为明显的低信号；而对于静脉血来讲，它和周围的脑实质氢质子进动频率的不同将导致在特定的TE时间的不同相，使静脉血信号与周围信号部分抵消更低。其图像后处理技术也是显示静脉血信号的关键。它是利用含有各种组织间磁敏感性对比信息的图像数据实现图像增强的一种后处理方法。它将图像数据以一定的值域转换规则制成相位蒙板，在利用蒙板对幅值图像进行加权运算，从而进一步提高幅值图像中各组织的对比性。然而，由于组织-空气、组织-骨骼交界处磁感应强度的突变、主磁场的不均匀性使得图像数据产生磁场不均一性伪影，即引入了背景相位伪影。这对交界处乃至整个静脉血管的成像具有不可忽略的影响。在常规的后处理方法中采用高通滤波器滤除包含背景相位伪影的低频信息，但受所选窗宽影响,会丢失有用的信号。同时在求相位过程中不可避免相位卷绕问题，而解相位卷绕是一个相当复杂的过程，且会引入新的背景伪影。为解决背景相位变化伪影,该领域技术人员提出多回波的采集方式,但以损失采集速度为前提;同时也有许多学者提出各种后处理方式,比如空间依赖伪影校正等。但是存在算法复杂，处理速度慢等问题。
因此，现有技术还有待改进和发展。
发明内容
本发明要解决的技术问题在于，针对现有技术的不足，提供一种磁共振磁敏感加权成像的后处理方法及系统，旨在解决现有技术中磁敏感加权成像的后处理方法主要是采用高通滤波器滤除包含背景相位伪影的低频信息，其会受所选窗宽影响,造成丢失有用的信号的问题。
为了解决上述技术问题，本发明所采用的技术方案如下：
一种磁共振磁敏感加权成像后处理方法，其包括：
提取原始图像的相位信息,并根据所述相位信息生成相位蒙板；
利用局部场变化分布生成所述原始图像的局部场变化蒙板；
根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板；
将所述最终蒙板与所述原始图像对应的幅值图像进行N次幂加权得到磁敏感加权图像。
所述磁共振磁敏感加权成像后处理方法，其中，所述提取原始图像的图像数据,并根据所述图像数据生成相位蒙板具体包括：
提取所述原始图像的相位图像以及幅值图像，并对所述第一图像数据进行反傅里叶变换得到幅值为1的K数据空间；
利用低通滤波器对所述K数据空间进行滤波，并对滤波后的K空间数据进行傅里叶变换得到滤波后的第二图像数据，其中，所述相位图像为第一图像数据；
根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板。
所述磁共振磁敏感加权成像后处理方法，其中，所述根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板具体为：
将所述第一图像数据和第二图像数据相除以得到相位；
将所述相位按照预设值域转换准则生成原图像的相位蒙板，其中，所述预设值域转换准则为：若所述相位大于等于0，则相位蒙板为1，反之为0。
所述磁共振磁敏感加权成像后处理方法，其中，所述利用局部场变化分布生成所述原始图像的局部场变化蒙板具体包括：
将所述幅值为1的K数据空间减去低通滤波后的K空间数据得到当前K空间数据；
对所述当前K空间数据进行傅里叶变换得到第三图像数据，其中，所述第三图像数据包括当前K空间数据对应的图像的第一幅值；
将第一幅值减去滤波后的K空间数据对应的图像的第二幅值得到局部场变化，并根据所述局部场变化生成局部场变化蒙板。
所述磁共振磁敏感加权成像后处理方法，其中，所述根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板具体为：
将所述相位蒙板与所述局部场变化蒙板按照预设准则相结合得到所述原始图像的最终蒙板，其中，所述预设准则为：
最终蒙板=1-0.5*局部场变化蒙板*相位蒙板。
一种磁共振磁敏感加权成像后处理系统，其包括：
相位蒙板模块，用于提取原始图像的相位信息,并根据所述相位信息生成相位蒙板；
局部场分布蒙板模块，用于利用局部场变化分布生成所述原始图像的局部场变化蒙板；
最终蒙板模块，用于根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板；
加权模块，用于将所述最终蒙板与所述原始图像对应的幅值图像进行N次幂加权得到磁敏感加权图像。
所述磁共振磁敏感加权成像后处理系统，其中，所述相位蒙板模块具体包括：
提取单元，用于提取所述原始图像的相位图像以及幅值图像，并对所述第一图像数据进行反傅里叶变换得到幅值为1的K数据空间，其中，所述相位图像为第一图像数据；
滤波单元，用于利用低通滤波器对所述K数据空间进行滤波，并对滤波后的K空间数据进行傅里叶变换得到滤波后的第二图像数据；
转换单元，用于根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板。
所述磁共振磁敏感加权成像后处理系统，其中，所述局部场分布蒙板模块具体包括：
计算单元，用于将所述幅值为1的K数据空间减去低通滤波后的K空间数据得到当前K空间数据；
确定单元，用于对所述当前K空间数据进行傅里叶变换得到第三图像数据，其中，所述第三图像数据包括当前K空间数据对应的图像的第一幅值；
生成单元，用于将第一幅值减去滤波后的K空间数据对应的图像的第二幅值得到局部场变化，并根据所述局部场变化生成局部场变化蒙板。
所述磁共振磁敏感加权成像后处理系统，其中，所述转换单元具体包括：
计算子单元，用于将所述第一图像数据和第二图像数据相除以得到相位；
生成子单元，用于将所述相位按照预设值域转换准则生成原图像的相位蒙板，其中，所述预设值域转换准则为：若所述相位大于等于0，则相位蒙板为1，反之为0。
所述磁共振磁敏感加权成像后处理系统，其中，所述最终蒙板模块具体用于：
将所述相位蒙板与所述局部场变化蒙板按照预设准则相结合得到所述原始图像的最终蒙板，其中，所述预设准则为：
最终蒙板=1-0.5*局部场变化蒙板*相位蒙板。
有益效果：与现有技术相比，本发明所提供一种磁共振磁敏感加权成像的后处理方法及系统，所述方法利用局部场变化分布生成局部场变化蒙板；再利用低通滤波器对原始K空间数据进行滤波将滤波前后的图像相除提取图像数据，按照正相位的值域转换准则生成相位蒙板；最后将两者按照预设准则结合生成最终蒙板，再结合幅值图像经过N次幂加权得到磁敏感加权图像。本发明通过采用局部场变化分布生成局部场变化蒙板结合图像的相位信息对所述磁共振图像进行后处理，减少了图像内组织交界处的伪影，从而避免了伪影受所选窗宽影响而造成丢失细小静脉信号的问题，提高了磁共振磁敏感加权图像的细节显示。
附图说明
图1为本发明磁共振磁敏感加权成像的后处理方法较佳实施例的流程图。
图2为本发明实施实例中磁敏感加权成像数据后处理的结果之一。
图3为本发明磁共振磁敏感加权成像的后处理系统的结构原理图。
具体实施方式
本发明提供磁共振磁敏感加权成像的后处理方法及方法，为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
下面结合附图，通过对实施例的描述，对发明内容作进一步说明。
请参见图1，图1为本发明提供磁共振磁敏感加权成像的后处理方法较佳实施例的流程图。所述方法包括：
S100、提取原始图像的相位信息,并根据所述相位信息生成相位蒙板；
S200、利用局部场变化分布生成所述原始图像的局部场变化蒙板；
S300、根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板；
S400、将所述最终蒙板与所述原始图像对应的幅值图像进行N次幂加权得到磁敏感加权图像。
本实施例提供的磁共振磁敏感加权成像的后处理方法，其利用局部场变化分布生成局部场变化蒙板；再利用低通滤波器对原始K空间数据进行滤波将滤波前后的图像相除提取图像数据，按照正相位的值域转换准则生成相位蒙板；最后将两者按照预设准则结合生成最终蒙板，再结合幅值图像经过N次幂加权得到磁敏感加权图像。本发明通过采用局部场变化分布生成局部场变化蒙板结合图像的相位信息对所述磁共振图像进行后处理，减少了图像内组织交界处的伪影，从而避免了伪影受所选窗宽影响而造成丢失细小静脉信号的问题，提高了磁共振磁敏感加权图像的细节显示。
在本实施例中，所述步骤S100具体可以包括：
S101、提取所述原始图像的相位图像以及幅值图像，并对所述第一图像数据进行反傅里叶变换得到幅值为1的K数据空间，其中，所述相位图像为第一图像数据；
S102、利用低通滤波器对所述K数据空间进行滤波，并对滤波后的K空间数据进行傅里叶变换得到滤波后的第二图像数据；
S103、根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板。
具体地，所述第一图像数据和的人图像数据均为复数域上的图像数据。所述对所述K空间数据进行滤波为利用低通滤波器对所述K空间数据进行滤波，并对利用低通滤波器进行滤波后的K空间数据进行傅里叶变换得到滤波后的第二图像数据。
在本实施例中，所述根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板具体可以包括：
S1031、将所述第一图像数据和第二图像数据相除以得到相位；
S1032、将所述相位按照预设值域转换准则生成原图像的相位蒙板，其中，所述预设值域转换准则为：若所述相位大于等于0，则相位蒙板为1，反之为0。
具体地，所述值域转换准则为预预先设置，用于根据相位生成相位蒙板。在本实施例中，所述值域转换准则为用于输出正相位蒙板的值域转换准则，其表达式可以为：
其中, 
为相位蒙板，
，
表示低通滤波，
为原图像的幅值图像，
为相位。
在本实施例的变形实施例中，所述值域转换准则也可以为用于输出负相位蒙板的值域转换准则，其表达式可以为：
当所述相位蒙板为负相位蒙板时，所述值域转换准则为：
其中, 
为相位蒙板，
，
表示低通滤波，
为原图像的幅值图像，
为相位。
在所述步骤S200中，所述局部场变化蒙板为用于反应局部磁敏感变化的局部场变化蒙板。所述局部场变化分布（Local Field Variation Map, LFV）可以通过低通滤波器获取，所述局部变化分布的获取与低通滤波器的关系可以为：
其中，
表示取模运算，
表示低通滤波，
为原图像的幅值图像，
为相位，
表示高通滤波。
进一步，从上述公式(3)可看出，获取的LFV的输出范围在[0,2]之间。从而，在生成局部场变化蒙板时对LFV进行归一化处理，所述归一化处理的公式可以为：
其中，
为局部场变化蒙板，
为局部场变化分布。
在本实施例中，所述利用局部场变化分布生成所述原始图像的局部场变化蒙板具体可以包括：
S201、将所述幅值为1的K数据空间减去低通滤波后的K空间数据得到当前K空间数据；
S202、对所述当前K空间数据进行傅里叶变换得到第三图像数据，其中，所述第三图像数据包括当前K空间数据对应的图像的第一幅值；
S203、将第一幅值减去过滤后的K空间数据对应的图像的第二幅值得到局部场变化，并根据所述局部场变化生成局部场变化蒙板。
进一步，在所述步骤S300中，所述根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板指的是将所述相位蒙板与所述局部场变化蒙板按照预设准则相结合得到所述原始图像的最终蒙板。所述相位蒙板于局部场变化蒙板的结合准则为预先设定的。在本实施例中，所述预设准则为：
其中，
为最终蒙板，
为局部场变化蒙板，
为相位蒙板。
在所述步骤S400中，所述N次幂加权为预先设置的，其可以根据对图像处理的精度或者处理速度的限制而设定所述加权的幂次，如2次幂加权，3次幂加权等。
为了进一步说明本发明提供的磁共振磁敏感加权成像后处理方法，下面给出一个具体例子加以说明。在本实施例中，所采集的片数为36片，层厚1mm,投影层数为9层；图2是本发明实施实例中上述磁敏感加权成像数据后处理的结果之一。在第一组效果图中，圆圈所示的组织交界处的伪影减少。在第二组效果图中，箭头指示的静脉细节显示效果好。
本发明还提供了一种磁共振磁敏感加权成像后处理系统，如图3所示，其包括：
相位蒙板模块100，用于提取原始图像的相位信息,并根据所述相位信息生成相位蒙板；
局部场分布蒙板模块200，用于利用局部场变化分布生成所述原始图像的局部场变化蒙板；
最终蒙板模块300，用于根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板；
加权模块400，用于将所述最终蒙板与所述原始图像对应的幅值图像进行N次幂加权得到磁敏感加权图像。
所述磁共振磁敏感加权成像后处理系统，其中，所述相位蒙板模块具体包括：
提取单元，用于提取所述原始图像的相位图像以及幅值图像，并对所述第一图像数据进行反傅里叶变换得到幅值为1的K数据空间；
滤波单元，用于利用低通滤波器对所述K数据空间进行滤波，并对滤波后的K空间数据进行傅里叶变换得到滤波后的第二图像数据；
转换单元，用于根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板。
所述磁共振磁敏感加权成像后处理系统，其中，所述局部场分布蒙板模块具体包括：
计算单元，用于将所述幅值为1的K数据空间减去低通滤波后的K空间数据得到当前K空间数据；
确定单元，用于对所述当前K空间数据进行傅里叶变换得到第三图像数据，其中，所述第三图像数据包括当前K空间数据对应的图像的第一幅值；
生成单元，用于将第一幅值减去滤波后的K空间数据对应的图像的第二幅值得到局部场变化蒙板。
所述磁共振磁敏感加权成像后处理系统，其中，所述转换单元具体包括：
计算子单元，用于将所述第一图像数据和第二图像数据相除以得到相位；
生成子单元，用于将所述相位按照预设值域转换准则生成原图像的相位蒙板，其中，所述预设值域转换准则为：若所述相位大于等于0，则相位蒙板为1，反之为0。
所述磁共振磁敏感加权成像后处理系统，其中，所述最终蒙板模块具体用于：
将所述相位蒙板与所述局部场变化蒙板按照预设准则相结合得到所述原始图像的最终蒙板，其中，所述预设准则为：
最终蒙板=1-0.5*局部场变化蒙板*相位蒙板。
上述磁共振磁敏感加权成像后处理系统各个单元模块已经在上述方法中进行了详细介绍，这里就不再赘述了。
可以理解的是，对本领域普通技术人员来说，可以根据本发明的技术方案及其发明构思加以等同替换或改变，而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。
Claims (10)
1.一种磁共振磁敏感加权成像后处理方法，其特征在于，其包括：
提取原始图像的相位信息,并根据所述相位信息生成相位蒙板；
利用局部场变化分布生成所述原始图像的局部场变化蒙板；
根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板；
将所述最终蒙板与所述原始图像对应的幅值图像进行N次幂加权得到磁敏感加权图像；
其中，所述提取原始图像的相位信息,并根据所述相位信息生成相位蒙板具体为：
提取所述原始图像的相位图像以及幅值图像，获取相位图像幅值为1的K空间数据；
利用低通滤波器对所述幅值为1的K空间数据进行滤波，并将滤波前后的相位图像相除以提取相位；
将提取到的相位按照预设值域转换准则生成相位蒙板；
其中，所述利用局部场变化分布生成所述原始图像的局部场变化蒙板具体为：
将所述幅值为1的K空间数据减去低通滤波后的K空间数据得到当前K空间数据；
获取当前K空间数据对应的图像的幅值，将当前K空间数据对应的图像的幅值减去滤波后的K空间数据对应的图像的幅值得到局部场变化，并根据所述局部场变化生成局部场变化蒙板。
2.根据权利要求1所述磁共振磁敏感加权成像后处理方法，其特征在于，所述提取原始图像的相位信息,并根据所述相位信息生成相位蒙板具体包括：
提取所述原始图像的相位图像以及幅值图像，并对第一图像数据进行反傅里叶变换得到幅值为1的K空间数据，其中，所述相位图像为第一图像数据；
利用低通滤波器对所述K空间数据进行滤波，并对滤波后的K空间数据进行傅里叶变换得到滤波后的第二图像数据；
根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板。
3.根据权利要求2所述磁共振磁敏感加权成像后处理方法，其特征在于，所述根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板具体为：
将所述第一图像数据和第二图像数据相除以得到相位；
将所述相位按照预设值域转换准则生成原图像的相位蒙板，其中，所述预设值域转换准则为：若所述相位大于等于0，则相位蒙板为1，反之为0。
4.根据权利要求2所述磁共振磁敏感加权成像后处理方法，其特征在于，所述利用局部场变化分布生成所述原始图像的局部场变化蒙板具体包括：
将所述幅值为1的K空间数据减去低通滤波后的K空间数据得到当前K空间数据；
对所述当前K空间数据进行傅里叶变换得到第三图像数据，其中，所述第三图像数据包括当前K空间数据对应的图像的第一幅值；
将第一幅值减去滤波后的K空间数据对应的图像的第二幅值得到局部场变化，并根据所述局部场变化生成局部场变化蒙板。
5.根据权利要求1所述磁共振磁敏感加权成像后处理方法，其特征在于，所述根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板具体为：
将所述相位蒙板与所述局部场变化蒙板按照预设准则相结合得到所述原始图像的最终蒙板，其中，所述预设准则为：
最终蒙板＝1-0.5*局部场变化蒙板*相位蒙板。
6.一种磁共振磁敏感加权成像后处理系统，其特征在于，其包括：
相位蒙板模块，用于提取原始图像的相位信息,并根据所述相位信息生成相位蒙板；
局部场分布蒙板模块，用于利用局部场变化分布生成所述原始图像的局部场变化蒙板；
最终蒙板模块，用于根据所述相位蒙板以及所述局部场变化蒙板得到所述原始图像的最终蒙板；
加权模块，用于将所述最终蒙板与所述原始图像对应的幅值图像进行N次幂加权得到磁敏感加权图像；
其中，所述提取原始图像的相位信息,并根据所述相位信息生成相位蒙板具体为：
提取所述原始图像的相位图像以及幅值图像，获取相位图像幅值为1的K空间数据；
利用低通滤波器对所述幅值为1的K空间数据进行滤波，并将滤波前后的相位图像相除以提取相位；
将提取到的相位按照预设值域转换准则生成相位蒙板；
其中，所述利用局部场变化分布生成所述原始图像的局部场变化蒙板具体为：
将所述幅值为1的K空间数据减去低通滤波后的K空间数据得到当前K空间数据；
获取当前K空间数据对应的图像的幅值，将当前K空间数据对应的图像的幅值减去滤波后的K空间数据对应的图像的幅值得到局部场变化，并根据所述局部场变化生成局部场变化蒙板。
7.根据权利要求6所述磁共振磁敏感加权成像后处理系统，其特征在于，所述相位蒙板模块具体包括：
提取单元，用于提取所述原始图像的相位图像以及幅值图像，并对所述相位图像进行反傅里叶变换得到幅值为1的K空间数据，其中，所述相位图像为第一图像数据；
滤波单元，用于利用低通滤波器对所述K空间数据进行滤波，并对滤波后的K空间数据进行傅里叶变换得到滤波后的第二图像数据；
转换单元，用于根据第一图像数据和第二图像数据确定相位，并将所述相位按照预设值域转换准则转换成相位蒙板。
8.根据权利要求7所述磁共振磁敏感加权成像后处理系统，其特征在于，所述局部场分布蒙板模块具体包括：
计算单元，用于将所述幅值为1的K空间数据减去低通滤波后的K空间数据得到当前K空间数据；
确定单元，用于对所述当前K空间数据进行傅里叶变换得到第三图像数据，其中，所述第三图像数据包括当前K空间数据对应的图像的第一幅值；
生成单元，用于将第一幅值减去滤波后的K空间数据对应的图像的第二幅值得到局部场变化，并根据所述局部场变化生成局部场变化蒙板。
9.根据权利要求7所述磁共振磁敏感加权成像后处理系统，其特征在于，所述转换单元具体包括：
计算子单元，用于将所述第一图像数据和第二图像数据相除以得到相位；
生成子单元，用于将所述相位按照预设值域转换准则生成原图像的相位蒙板，其中，所述预设值域转换准则为：若所述相位大于等于0，则相位蒙板为1，反之为0。
10.根据权利要求6所述磁共振磁敏感加权成像后处理系统，其特征在于，所述最终蒙板模块具体用于：
将所述相位蒙板与所述局部场变化蒙板按照预设准则相结合得到所述原始图像的最终蒙板，其中，所述预设准则为：
最终蒙板＝1-0.5*局部场变化蒙板*相位蒙板。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710149852.0A CN106960458B (zh) | 2017-03-14 | 2017-03-14 | 一种磁共振磁敏感加权成像后处理方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710149852.0A CN106960458B (zh) | 2017-03-14 | 2017-03-14 | 一种磁共振磁敏感加权成像后处理方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106960458A CN106960458A (zh) | 2017-07-18 |
| CN106960458B true CN106960458B (zh) | 2020-08-25 |
Family
ID=59470433
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710149852.0A Active CN106960458B (zh) | 2017-03-14 | 2017-03-14 | 一种磁共振磁敏感加权成像后处理方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106960458B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110634119B (zh) * | 2018-06-05 | 2021-11-02 | 中国科学院深圳先进技术研究院 | 分割磁敏感加权图像中静脉血管的方法、装置和计算设备 |
| CN109934807B (zh) * | 2019-03-06 | 2021-03-23 | 上海东软医疗科技有限公司 | 一种定量磁化率成像的方法及装置 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101143094A (zh) * | 2006-09-13 | 2008-03-19 | 株式会社东芝 | 磁共振图像诊断装置及其控制方法 |
| CN103108586A (zh) * | 2010-04-27 | 2013-05-15 | 磁共振创新公司 | 使用磁敏感加权成像和磁敏感填图(swim)来生成核磁共振图像的方法 |
| CN104126193A (zh) * | 2012-02-14 | 2014-10-29 | 皇家飞利浦有限公司 | 图像分辨率增强 |
| CN104749538A (zh) * | 2015-04-30 | 2015-07-01 | 郑州轻工业学院 | 一种并行磁共振成像相位处理方法 |
| CN105701815A (zh) * | 2016-01-12 | 2016-06-22 | 深圳安科高技术股份有限公司 | 一种磁共振灌注成像后处理方法及系统 |
| CN105997125A (zh) * | 2016-06-15 | 2016-10-12 | 刘丽 | 多针孔单光子spect心肌血流绝对定量方法与用途 |
| CN106339982A (zh) * | 2016-08-24 | 2017-01-18 | 深圳先进技术研究院 | 快速磁共振心脏实时电影成像方法及系统 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8487615B2 (en) * | 2011-04-20 | 2013-07-16 | Imris Inc | Magnetic resonance signal detection using remotely positioned receive coils |
-
2017
- 2017-03-14 CN CN201710149852.0A patent/CN106960458B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101143094A (zh) * | 2006-09-13 | 2008-03-19 | 株式会社东芝 | 磁共振图像诊断装置及其控制方法 |
| CN103108586A (zh) * | 2010-04-27 | 2013-05-15 | 磁共振创新公司 | 使用磁敏感加权成像和磁敏感填图(swim)来生成核磁共振图像的方法 |
| CN104126193A (zh) * | 2012-02-14 | 2014-10-29 | 皇家飞利浦有限公司 | 图像分辨率增强 |
| CN104749538A (zh) * | 2015-04-30 | 2015-07-01 | 郑州轻工业学院 | 一种并行磁共振成像相位处理方法 |
| CN105701815A (zh) * | 2016-01-12 | 2016-06-22 | 深圳安科高技术股份有限公司 | 一种磁共振灌注成像后处理方法及系统 |
| CN105997125A (zh) * | 2016-06-15 | 2016-10-12 | 刘丽 | 多针孔单光子spect心肌血流绝对定量方法与用途 |
| CN106339982A (zh) * | 2016-08-24 | 2017-01-18 | 深圳先进技术研究院 | 快速磁共振心脏实时电影成像方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 基于K空间的MR截断伪影去除算法;王甜甜，等;《北京生物医学工程》;20150831;第375-378页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106960458A (zh) | 2017-07-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP2148212B1 (en) | Magnetic resonance imaging apparatus for contrast enhancement of images | |
| Rauscher et al. | Magnetic susceptibility-weighted MR phase imaging of the human brain | |
| Feng et al. | Catalytic multiecho phase unwrapping scheme (CAMPUS) in multiecho gradient echo imaging: removing phase wraps on a voxel‐by‐voxel basis | |
| Soellinger et al. | 3D cine displacement‐encoded MRI of pulsatile brain motion | |
| JP5394374B2 (ja) | 磁気共鳴イメージング装置及び血管画像取得方法 | |
| CN108778116B (zh) | 磁共振成像装置以及图像处理方法 | |
| US20160012569A1 (en) | Image processing device, magnetic resonance imaging apparatus and image processing method | |
| JP2012070939A (ja) | 磁気共鳴イメージング装置 | |
| JP2015116474A (ja) | 磁気共鳴イメージング装置 | |
| Du et al. | Multi‐echo acquisition of MR angiography and venography of the brain at 3 Tesla | |
| CN110174632B (zh) | 磁共振成像方法、装置、成像设备及磁共振成像系统 | |
| Cao et al. | Optimized multi‐axis spiral projection MR fingerprinting with subspace reconstruction for rapid whole‐brain high‐isotropic‐resolution quantitative imaging | |
| Guo et al. | Multi‐slice myelin water imaging for practical clinical applications at 3.0 T | |
| CN106232002B (zh) | 具有增强的稳定性和速度的电导率重建的ept方法 | |
| JP5417050B2 (ja) | 磁気共鳴イメージング装置 | |
| CN106960458B (zh) | 一种磁共振磁敏感加权成像后处理方法及系统 | |
| WO2014038441A1 (ja) | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 | |
| Kan et al. | Quantitative susceptibility mapping using principles of echo shifting with a train of observations sequence on 1.5 T MRI | |
| Shiraishi et al. | Deep learning–based reconstruction and 3D hybrid profile order technique for MRCP at 3T: evaluation of image quality and acquisition time | |
| Tian et al. | Real‐time water/fat imaging at 0. 55T with spiral out‐in‐out‐in sampling | |
| Wei et al. | Imaging diamagnetic susceptibility of collagen in hepatic fibrosis using susceptibility tensor imaging | |
| US10928475B2 (en) | Dynamic contrast enhanced magnetic resonance imaging with flow encoding | |
| Manduca et al. | Improved inversion of MR elastography images by spatiotemporal directional filtering | |
| Deistung et al. | GUIBOLD: a graphical user interface for image reconstruction and data analysis in susceptibility-weighted MR imaging | |
| Bao et al. | Whole brain susceptibility mapping using harmonic incompatibility removal |
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 |