CN104107045B - 磁共振成像方法和装置 - Google Patents
磁共振成像方法和装置 Download PDFInfo
- Publication number
- CN104107045B CN104107045B CN201410306978.0A CN201410306978A CN104107045B CN 104107045 B CN104107045 B CN 104107045B CN 201410306978 A CN201410306978 A CN 201410306978A CN 104107045 B CN104107045 B CN 104107045B
- Authority
- CN
- China
- Prior art keywords
- image data
- value
- image
- values
- real part
- 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
- 238000003384 imaging method Methods 0.000 title abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 68
- 238000012545 processing Methods 0.000 claims abstract description 59
- 238000001914 filtration Methods 0.000 claims abstract description 30
- 230000005291 magnetic effect Effects 0.000 claims abstract description 30
- 230000011218 segmentation Effects 0.000 claims abstract description 15
- 238000010606 normalization Methods 0.000 claims description 37
- 230000008569 process Effects 0.000 claims description 29
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 20
- 230000009466 transformation Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000006870 function Effects 0.000 description 36
- 210000001519 tissue Anatomy 0.000 description 9
- 210000004204 blood vessel Anatomy 0.000 description 6
- 210000003462 vein Anatomy 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000003672 processing method Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000005298 paramagnetic effect Effects 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 2
- 210000001175 cerebrospinal fluid Anatomy 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 230000004069 differentiation Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 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 1
- 208000003174 Brain Neoplasms Diseases 0.000 description 1
- 206010008111 Cerebral haemorrhage Diseases 0.000 description 1
- 108010064719 Oxyhemoglobins Proteins 0.000 description 1
- 108010001267 Protein Subunits Proteins 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 108010002255 deoxyhemoglobin Proteins 0.000 description 1
- 230000005292 diamagnetic effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明实施例公开了一种磁共振成像方法。该方法包括:对原始图像进行复数高通滤波,获得滤波后的图像;从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理;将处理后的数值作为权重值与幅值图像的灰度值相乘,最终得到磁共振图像。根据本发明实施例,可以在保证图像清晰度的同时,提供一种快速便捷的磁共振成像方法。解决现有技术中磁敏感加权图像处理过程过于复杂的问题。本发明实施例还提供了一种磁共振成像装置。
Description
技术领域
本发明涉及医学图像处理领域,特别是涉及磁共振成像方法和装置。
背景技术
磁敏感加权成像(Susceptibility Weighted Imaging,SWI)技术是近年来新发展起来的一种磁共振对比增强成像技术。不同于以往的质子密度、T1或T2加权成像,SWI根据不同组织间的磁敏感性(磁敏感性可以量化为磁化率)差异提供图像对比增强。
研究发现,由于静脉血的主要成份为顺磁性的去氧血红蛋白,动脉血的主要成份则是抗磁性的氧合血红蛋白,因此,静脉血管和动脉血管之间存在磁敏感性差异,这种差异会最终导致两种血管的信号强度不同,从而为静脉血管能够独立于动脉血管清晰成像提供了可能。在临床应用上,SWI技术可应用于脑肿瘤、脑出血或其它与静脉血管有关的病灶研究中。
目前,通过现有的磁共振成像设备还无法直接得到SWI图像,而是需要对基于T2*加权梯度回波序列获得的原始图像(即,实际检测到的图像)进行复杂的图像处理才能得到SWI图像。
在现有技术中,提供了一种获得SWI图像的图像处理方法,其处理流程为:先根据原始图像数据生成相位图像和幅值图像(或者称为磁距图像);然后对相位图像进行高通滤波,以去除由于背景磁场的不均匀造成的低频扰动,得到滤波后的相位图像;再用滤波后的相位图像计算相位图像掩模(或相位蒙片);最后将相位图像掩模应用到幅值图像n次,得到最终的SWI图像。
在实现本发明的过程中,本发明的发明人发现现有技术中至少存在如下问题:由于现有的这种图像处理方法需要计算相位图像,并根据相位图像计算相位蒙片n次幂,而幂指数n的取值由图像信噪比决定(即,在不同的图像信噪比下,n的取值也会不同),因此,在图像处理过程中还需要先确定图像信噪比,从而增加了整个图像处理过程的复杂性。并且,n次幂的计算本身也使整个图像处理过程相对比较繁琐。
发明内容
为了解决上述技术问题,本发明实施例提供了磁共振成像方法和装置,以在保证图像清晰度的同时,解决现有技术中图像处理过程过于复杂的问题。
本发明实施例公开了如下技术方案:
一种磁共振成像方法,包括:
对原始图像进行复数高通滤波,获得滤波后的图像;
从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理;
将处理后的数值作为权重值与幅值图像的灰度值相乘,最终得到磁共振图像。
优选的,所述方法还包括:
对获得的各个层面所述磁共振图像进行最小信号强度投影。
优选的,所述复数高通滤波为同态高通滤波。
优选的,所述从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理,包括:
当所述图像数据中的虚部数值大于或等于0时,将所述图像数据中的实部数值归一化为1;当虚部数值小于0时,对所述图像数据中的实部数值归一化为位于0和1之间的值;
或者,
当所述图像数据中的虚部数值大于或等于0时,将所述图像数据中的实部数值归一化为位于0和1之间的值;当虚部数值小于0时,将所述图像数据中的实部数值归一化为1;
或者,
直接将所述图像数据中的实部数值归一化为位于0和1之间的值。
进一步优选的,所述将所述图像数据中的实部数值归一化为位于0和1之间的值具体为:
按照将所述图像数据中的实部数值归一化为位于0 和1之间的值;
其中,W(r)为处理后的数值,R(ρHF(r))为实部数值。
一种磁共振成像装置,包括:
高通滤波单元,用于对原始图像进行复数高通滤波,获得滤波后的图像;
归一化处理单元,用于从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理;
加权计算单元,用于将处理后的数值作为权重值与幅值图像的灰度值相乘,得到磁共振图像。
优选的,所述装置还包括:
投影单元,用于对获得的各个层面所述磁共振图像进行最小信号强度投影。
优选的,所述复数高通滤波为同态高通滤波。
优选的,所述归一化处理单元包括:
第一设置子单元,用于当所述图像数据中的虚部数值大于或等于0,将所述图像数据中的实部数值归一化为1;
第一变换子单元,用于当所述图像数据中的虚部数值小于0,将所述图像数据中的实部数值归一化为位于0和1之间的值;
或者,
所述归一化处理单元包括:
第二设置子单元,用于当所述图像数据中的虚部数值小于0,将所述图像数据中的实部数值归一化为1;
第二变换子单元,用于当所述图像数据中的虚部数值大于或等于0,将所述图像数据中的实部数值归一化为位于0和1之间的值;
或者,
所述归一化处理单元包括:
第三变换子单元,直接将所述图像数据的实部数值归一化为位于0和1之间的值。
进一步优选的,所述第一变换子单元、第二变换子单元或者第三变换子 单元具体用于,
按照将所述图像数据中的实部数值归一化为位于0和1之间的值;
其中,W(r)为处理后的数值,R(ρHF(r))为实部数值。
由上述实施例可以看出,与现有技术相比,本发明技术方案的优点在于:
将原始图像经滤波的图像数据进行不同的归一化处理,利用处理后的权重函数与现有技术中相位蒙片的幂指数函数(即,(φMASK(r))n)之间的相似性,由处理后的权重函数替换现有技术中相位蒙片的幂指数函数。由于该处理后的权重函数不随图像信噪比的变化而变化,稳定性较好。因此,当采用该处理后的权重函数替换相位蒙片的幂指数函数,在图像处理过程中不需要计算相位图像,不需要通过图像信噪比确定幂次,从而简化整个图像处理过程。并且,无需进行n次幂的计算,也使整个图像处理过程相对简单容易。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例一提供的一种磁共振成像方法的流程图;
图2为本发明的权重曲线与现有技术的幂指数曲线的一种相似对比图;
图3为本发明的权重曲线与现有技术的幂指数曲线的另一种相似对比图;
图4为本发明的权重曲线与现有技术的幂指数曲线的另一种相似对比图;
图5为本发明实施例二提供的一种磁共振成像方法的流程图;
图6为本发明实施例三提供的一种磁共振成像装置的结构框图;
图7为本发明实施例三提供的另一种磁共振成像装置的结构框图;
图8为本发明中归一化处理单元的一种结构框图;
图9(a)-(b)分别为原始图像与本发明处理后的图像。
具体实施方式
为了去除背景磁场不均匀造成的低频相位干扰,并进一步增强组织间的磁敏感对比度,从而更加清晰地显示解剖结构,需要对原始图像进行一系列的后处理。
在现有技术中,根据原始图像的图像数据(该图像数据为复数)可以分别获得相位图像和幅值图像。其中,
相位图像为:
φ(r)=arctan(I(r)/R(r))
幅值图像为:
I(r)和R(r)分别为图像数据的虚部和实部,r为相位值。
常规的磁共振成像(Magnetic Resonance Imaging,MRI)只利用了单一的幅值信息,SWI技术则利用了一直被忽略的相位信息,并经过一系列的图像处理将相位图像应用到幅值图像中,形成独特的图像对比增强。该处理过程包括:
1、先对相位图像进行高通滤波,去除由于背景磁场不均匀造成的低频扰动。
例如,首先对相位图像进行低通滤波,然后在复数域中,用原始图像除以低通滤波后的k空间数据,最终得到的就是高通滤波后的相位图像。
2、再对滤波后的相位图像进行归一化处理,生成相位蒙片φMASK(r)。
3、最后根据相位蒙片和幅值图像,生成SWI图像:
ρSWI(r)=(φMASK(r))n×ρm(r),n决定了权重的大小,一般n取3~6可以得到信噪比较高的图像。
本发明的发明人在研究中发现,将原始图像经复数高通滤波处理后,如果将处理后的图像的图像数据中的实部数值按照对虚部数值的分段,分别进行不同的归一化处理,处理后所得到的权重函数与上述相位蒙片的幂指数函数(即,(φMASK(r))n)非常近似,并且,该权重函数不随图像信噪比的变化而变 化,稳定性较好。因此,当采用该权重函数替换相位蒙片的幂指数函数,在图像处理过程中就不需要计算相位图像及相位蒙片,并确定图像信噪比,从而简化整个图像处理过程。并且,无需进行n次幂的计算,也使整个图像处理过程相对简单容易。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图对本发明实施例进行详细描述。
实施例一
请参阅图1,其为本发明实施例一提供的一种磁共振成像方法的流程图,该方法包括以下步骤:
步骤101:对原始图像进行复数高通滤波,获得滤波后的图像。
其中,为了消除由于背景磁场的不均匀造成的低频扰动,获得更加清晰的图像,需要先对原始图像进行高通滤波。由于原始图像中的图像数据为复数,因此,该高通滤波为复数形式的高通滤波。
在本发明的一个优选实施方式中,可以采用同态高通滤波实现复数形式的高通滤波。
例如,原始图像为ρ(r),对原始图像进行同态高通滤波后的图像为:
ρHF(r)=exp(FFT-1(Han(r)×FFT(ln(ρ(r)))))
其中,ρHF(r)为高通滤波后的图像的图像数据,Han(r)为Hanning函数,
步骤102:从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理。
由于在相位图像中,顺磁性物质(如静脉血管)的相位值表现为明显的负值,而大部分脑实质及脑脊液等的相位值通常为正值或较小的负值,因此,在现有技术中,需要对相位图像计算相位掩模(或相位蒙片),并将相位掩模的n次幂与幅度图像幅值相乘,突出小结构的可见性。
由于经高通滤波图像数据中的虚部数值包含相位值的正弦函数信息,实 部包含相位值的余弦函数信息,因此,当顺磁性物质(如静脉血管)的相位值表现为明显的负值,该部分组织的图像数据中的虚部数值也会相应地表现为明显的负值,而当大部分脑实质及脑脊液等的相位值通常为正值或较小的负值时,该部分组织的图像数据的虚部数值也会相应地表现为正值。也就是说,不同组织的图像数据的虚部数值是不同的。因此,可以利用不同组织的图像数据的虚部数值是不同的特性,对获得的图像数据中的虚部数值进行分段,而分段后,位于同一个分段内的虚部数值即为同一组织的图像数据的虚部数值。然后对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理。即,针对不同的组织的图像数据的实部数值,采取差异化的归一处理方式,以便最终达到不同的组织差异化显示的效果。
后续将详细描述如何进行不同的归一化处理。
步骤103:将处理后的数值作为权重值与幅值图像的灰度值相乘,最终得到磁共振图像。
显然,最终获得的磁共振图像可以突出显示原始图像中的小结构,如,静脉血管。
在现有技术中,通常有三种对相位图像进行掩模处理,生成相位蒙片的方法。
第一种方法:当相位值位于[0,π]之间时,将相位值设置为1;当相位值位于[-π,0)之间时,将相位值变换为0和1之间的值,如,φ(r)为相位值。
第二种方法:当相位值位于[-π,0)之间时,将相位值设置为1;当相位值位于[0,π]之间时,将相位值变换为0和1之间的值。
第三种方法:当相位值位于[-π,π]之间时,将相位值变换为0和1之间的值。
对应于上述三种方法,上述步骤102中的不同的归一化处理方式可以为:
第一种方法、当所述图像数据中的虚部数值大于或等于0时,将所述图 像数据中的实部数值归一化为1;当虚部数值小于0时,将所述图像数据中的实部数值归一化为位于0和1之间的值;
第二种方法、当所述图像数据中的虚部数值大于或等于0时,将所述图像数据中的实部数值归一化为位于0和1之间的值;当虚部数值小于0时,将所述图像数据中的实部数值归一化为1。
第三种方法、直接将所述图像数据中的实部数值归一化为位于0和1之间的值。
在本发明的一个优选实施方式中,将所述图像数据中的实部数值归一化为位于0和1之间的值具体为:
按照将所述图像数据中的实部数值归一化为位于0和1之间的值;
其中,W(r)为归一化处理后的数值,R(ρHF(r))为实部数值。
例如,以第一种方法为例,分段归一化处理后的权重数值为:
其中,ρHF(r)为高通滤波后的图像的图像数据,I(ρHF(r))为该图像数据的虚部数值,R(ρHF(r))为该图像数据的实部数值。
最后将分段归一化处理后的W(r)作为权重值,与幅值图像ρM(r)进行点对点的乘积,最终得到突出显示小结构的磁共振图像为:
当采用该公式进行变换处理时,针对第一种归一化处理方法,本发明的权重曲线与现有技术的幂指数曲线的相似对比图如图2所示,针对第二种归一化处理方法,本发明的权重曲线与现有技术的幂指数曲线的相似对比图如图3所示,针对第三种归一化处理方法,本发明的权重曲线与现有技术的幂指数曲线的相似对比图如图4所示。其中,线1为相位蒙片的1次幂函数,线2为由归一化处理后的虚部数值构成的正弦函数,线3为相位蒙片的4次 幂函数。
由上述实施例可以看出,与现有技术相比,本发明技术方案的优点在于:
将滤波后的图像中的图像数据的虚部数值作为该图像数据的分段依据,利用实部数值对图像数据进行不同的归一化处理,利用处理后的权重函数与现有技术中相位蒙片的幂指数函数(即,(φMASK(r))n)之间的相似性,由处理后的权重替换现有技术中相位蒙片的幂指数函数。由于该处理后的权重函数不随图像信噪比的变化而变化,稳定性较好。因此,当采用该处理后的权重函数替换相位蒙片的幂指数函数,在图像处理过程中就不需要计算相位图像和相位蒙片,不需要通过图像信噪比确定幂次,从而简化整个图像处理过程。并且,无需进行n次幂的计算,也使整个图像处理过程相对简单容易。
实施例二
通过实施例一中磁共振成像方法,可以获得各个层面的磁共振图像。为了进一步使分散在各个层面的静脉血管的磁敏感信号连续化,最终显示出连续的静脉血管结构,本实施例二在实施例一的基础上,对获得的各个层面的磁共振图像进行最小信号强度投影。
请参阅图5,其为本发明实施例二提供的一种磁共振成像方法的流程图,该方法包括以下步骤:
步骤501:对原始图像进行复数高通滤波,获得滤波后的图像。
步骤502:从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理。
步骤503:将处理后的数值作为权重值与幅值图像的灰度值相乘,最终得到磁共振图像。
步骤504:对获得的各个层面所述磁共振图像进行最小信号强度投影。
由上述实施例可以看出,与现有技术相比,本发明技术方案的优点在于:
将滤波后的图像中的图像数据的虚部数值作为该图像数据的分段依据,利用实部数值对图像数据进行不同的归一化处理,利用处理后的权重函数与现有技术中相位蒙片的幂指数函数(即,(φMASK(r))n)之间的相似性,由处理后 的权重替换现有技术中相位蒙片的幂指数函数。由于该处理后的权重函数不随图像信噪比的变化而变化,稳定性较好。因此,当采用该处理后的权重函数替换相位蒙片的幂指数函数,在图像处理过程中就不需要计算相位图像和相位蒙片,不需要通过图像信噪比确定幂次,从而简化整个图像处理过程。并且,无需进行n次幂的计算,也使整个图像处理过程相对简单容易。
实施例三
与上述一种磁共振成像方法相对应,本发明实施例还提供了一种磁共振成像装置。请参阅图6,其为本发明实施例三提供的一种磁共振成像装置的结构框图,该装置包括:高通滤波单元601、归一化处理单元602和加权计算单元603。下面结合该装置的工作原理进一步介绍其内部结构以及连接关系。
高通滤波单元601,用于对原始图像进行复数高通滤波,获得滤波后的图像。
归一化处理单元602,用于从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理。
加权计算单元603,用于将处理后的数值作为权重值与幅值图像的灰度值相乘,得到磁共振图像。
在本发明的一个优选实施方式中,如图7所示,该装置还包括:
投影单元604,用于对获得的各个层面所述磁共振图像进行最小信号强度投影。
在本发明的另一个优选实施方式中,所述复数高通滤波为同态高通滤波。
在本发明的另一个优选实施方式中,如图8所示,归一化处理单元602包括:
第一设置子单元6021,用于当所述图像数据中的虚部数值大于或等于0,将所述图像数据中的实部数值归一化为1。
第一变换子单元6022,用于当所述图像数据中的虚部数值小于0,将所述图像数据中的实部数值归一化为位于0和1之间的值。
作为一种替换方案,归一化处理单元602包括:
第二设置子单元,用于当所述图像数据中的虚部数值小于0,将所述图像数据中的实部数值归一化为1。
第二变换子单元,用于当所述图像数据中的虚部数值大于或等于0,将所述图像数据中的实部数值归一化为位于0和1之间的值。
或者,作为另一种替换方案,归一化处理单元602包括:
第三变换子单元,直接将所述图像数据的实部数值归一化为位于0和1之间的值。
在本发明的另一优选实施方式中,所述第一变换子单元、第二变换子单元或者第三变换子单元具体用于,
按照将所述图像数据中的实部数值归一化为位于0和1之间的值;
其中,W(r)为处理后的数值,R(ρHF(r))为实部数值。
由上述实施例可以看出,与现有技术相比,本发明技术方案的优点在于:
将滤波后的图像中的图像数据的虚部数值进行分段,利用实部数值的不同的归一化处理得到加权图像数据,利用处理后的权重函数与现有技术中相位蒙片的幂指数函数(即,(φMASK(r))n)之间的相似性,由处理后的权重替换现有技术中相位蒙片的幂指数函数。由于该处理后的权重函数不随图像信噪比的变化而变化,稳定性较好。因此,当采用该处理后的权重函数替换相位蒙片的幂指数函数,在图像处理过程中就不需要计算相位图像和相位蒙片,不需要通过图像信噪比确定幂次,从而简化整个图像处理过程。并且,无需进行n次幂的计算,也使整个图像处理过程相对简单容易。
通过健康受试者的实验数据,验证本发明的技术方案可以得到突出显示小结构的效果。结果实例如图9所示,其中图9(a)为原始图像;图9(b)为经本方法处理的图像,突出显示了静脉。
所述领域的技术人员可以清楚地了解到,为了描述的方便和简洁,上述描述的系统、装置和单元的具体工作过程,可以参考前述方法实施例中的对 应过程,在此不再赘述。
在本发明所提供的几个实施例中,应该理解到,所揭露的系统、装置和方法,可以通过其它的方式实现。例如,以上所描述到的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性、机械或其它的形式。
所述作为分离部件说明的单元可以是或者也可以是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现,可以采用软件功能单元的形式实现。
需要说明的是,本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random Access Memory,RAM)等。
以上对本发明所提供的磁共振成像方法和装置进行了详细介绍,本文中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种磁共振成像方法,其特征在于,包括:
对原始图像进行复数高通滤波,获得滤波后的图像;
从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理;
将处理后的数值作为权重值与幅值图像的灰度值相乘,最终得到磁共振图像。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对获得的各个层面所述磁共振图像进行最小信号强度投影。
3.根据权利要求1所述的方法,其特征在于,所述复数高通滤波为同态高通滤波。
4.根据权利要求1所述的方法,其特征在于,所述从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理,包括:
当所述图像数据中的虚部数值大于或等于0时,将所述图像数据中的实部数值归一化为1;当虚部数值小于0时,对所述图像数据中的实部数值归一化为位于0和1之间的值;
或者,
当所述图像数据中的虚部数值大于或等于0时,将所述图像数据中的实部数值归一化为位于0和1之间的值;当虚部数值小于0时,将所述图像数据中的实部数值归一化为1;
或者,
直接将所述图像数据中的实部数值归一化为位于0和1之间的值。
5.根据权利要求4所述的方法,其特征在于,所述将所述图像数据中的实部数值归一化为位于0和1之间的值具体为:
按照将所述图像数据中的实部数值归一化为位于0和1之间的值;
其中,W(r)为处理后的数值,R(ρHF(r))为实部数值。
6.一种磁共振成像装置,其特征在于,包括:
高通滤波单元,用于对原始图像进行复数高通滤波,获得滤波后的图像;
归一化处理单元,用于从滤波后的图像中获得图像数据,对所述图像数据中的虚部数值进行分段,并对各分段范围内的所述图像数据中的实部数值进行不同的归一化处理;
加权计算单元,用于将处理后的数值作为权重值与幅值图像的灰度值相乘,得到磁共振图像。
7.根据权利要求6所述的装置,其特征在于,所述装置还包括:
投影单元,用于对获得的各个层面所述磁共振图像进行最小信号强度投影。
8.根据权利要求6所述的装置,其特征在于,所述复数高通滤波为同态高通滤波。
9.根据权利要求6所述的装置,其特征在于,
所述归一化处理单元包括:
第一设置子单元,用于当所述图像数据中的虚部数值大于或等于0,将所述图像数据中的实部数值归一化为1;
第一变换子单元,用于当所述图像数据中的虚部数值小于0,将所述图像数据中的实部数值归一化为位于0和1之间的值;
或者,
所述归一化处理单元包括:
第二设置子单元,用于当所述图像数据中的虚部数值小于0,将所述图像数据中的实部数值归一化为1;
第二变换子单元,用于当所述图像数据中的虚部数值大于或等于0,将所述图像数据中的实部数值归一化为位于0和1之间的值;
或者,
所述归一化处理单元包括:
第三变换子单元,直接将所述图像数据的实部数值归一化为位于0和1之间的值。
10.根据权利要求9所述的装置,其特征在于,所述第一变换子单元、第二变换子单元或者第三变换子单元具体用于,
按照将所述图像数据中的实部数值归一化为位于0和1之间的值;
其中,W(r)为处理后的数值,R(ρHF(r))为实部数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410306978.0A CN104107045B (zh) | 2014-06-30 | 2014-06-30 | 磁共振成像方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410306978.0A CN104107045B (zh) | 2014-06-30 | 2014-06-30 | 磁共振成像方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104107045A CN104107045A (zh) | 2014-10-22 |
CN104107045B true CN104107045B (zh) | 2016-08-17 |
Family
ID=51704172
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410306978.0A Active CN104107045B (zh) | 2014-06-30 | 2014-06-30 | 磁共振成像方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104107045B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104490393B (zh) * | 2014-12-17 | 2017-04-12 | 中国科学院深圳先进技术研究院 | 一种磁共振大脑血氧水平测量方法 |
CN105261051B (zh) * | 2015-09-25 | 2018-10-02 | 沈阳东软医疗系统有限公司 | 一种获取图像掩码的方法及装置 |
CN105785296A (zh) * | 2016-03-10 | 2016-07-20 | 哈尔滨医科大学 | 一种基于1H-19F-23Na原子核的多通道三维核磁共振成像方法 |
CN105842641A (zh) * | 2016-03-10 | 2016-08-10 | 哈尔滨医科大学 | 一种基于1h-19f-31p原子核的多通道三维核磁共振成像方法 |
CN108784695A (zh) * | 2018-06-12 | 2018-11-13 | 暨南大学附属第医院(广州华侨医院) | 磁敏感加权成像用于慢性大脑中动脉狭窄或闭塞的诊断系统 |
CN112767259B (zh) * | 2020-12-29 | 2024-12-24 | 上海联影智能医疗科技有限公司 | 图像处理方法、装置、计算机设备和存储介质 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2148212B1 (en) * | 2008-07-24 | 2019-08-21 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus for contrast enhancement of images |
US8452065B2 (en) * | 2009-03-11 | 2013-05-28 | Fred S. Azar | System for detecting malignant lymph nodes using an MR imaging device |
US9448289B2 (en) * | 2010-11-23 | 2016-09-20 | Cornell University | Background field removal method for MRI using projection onto dipole fields |
-
2014
- 2014-06-30 CN CN201410306978.0A patent/CN104107045B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104107045A (zh) | 2014-10-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104107045B (zh) | 磁共振成像方法和装置 | |
Tamada et al. | Motion artifact reduction using a convolutional neural network for dynamic contrast enhanced MR imaging of the liver | |
West et al. | BOLD hemodynamic response function changes significantly with healthy aging | |
Shimizu et al. | Wavelet-based multifractal analysis of fMRI time series | |
Zarinabad et al. | Voxel‐wise quantification of myocardial perfusion by cardiac magnetic resonance. Feasibility and methods comparison | |
WO2022183988A1 (en) | Systems and methods for magnetic resonance image reconstruction with denoising | |
JPWO2010073923A1 (ja) | 位相差強調画像化法(PhaseDifferenceEnhancedImaging;PADRE)、機能画像作成法、位相差強調画像化プログラム、位相差強調画像化装置、機能画像作成装置および磁気共鳴画像化(MagneticResonanceImaging;MRI)装置 | |
Delbany et al. | One‐millimeter isotropic breast diffusion‐weighted imaging: Evaluation of a superresolution strategy in terms of signal‐to‐noise ratio, sharpness and apparent diffusion coefficient | |
Westin et al. | Three‐dimensional adaptive filtering in magnetic resonance angiography | |
Barnhill et al. | Fast robust dejitter and interslice discontinuity removal in MRI phase acquisitions: application to magnetic resonance elastography | |
Li et al. | Quantitative assessment of susceptibility‐weighted imaging processing methods | |
Tabelow et al. | Local estimation of the noise level in MRI using structural adaptation | |
CN105849804A (zh) | 过滤噪声的计算高效方法 | |
CN106232002B (zh) | 具有增强的稳定性和速度的电导率重建的ept方法 | |
US8275181B2 (en) | Method for tracking of contrast enhancement pattern for pharmacokinetic and parametric analysis in fast-enhancing tissues using high-resolution MRI | |
Bouhrara et al. | Spatially adaptive unsupervised multispectral nonlocal filtering for improved cerebral blood flow mapping using arterial spin labeling magnetic resonance imaging | |
Gulban et al. | Improving a probabilistic cytoarchitectonic atlas of auditory cortex using a novel method for inter-individual alignment | |
James et al. | Impact of sampling rate on statistical significance for single subject fMRI connectivity analysis | |
Huang et al. | A kernel-based image denoising method for improving parametric image generation | |
Lindquist et al. | Spatial smoothing in fMRI using prolate spheroidal wave functions | |
JP6991728B2 (ja) | 画像処理装置、磁気共鳴イメージング装置及び画像処理方法 | |
CN105261051A (zh) | 一种获取图像掩码的方法及装置 | |
CN115984131A (zh) | 双维度图像边缘增强方法及应用 | |
CN109389651B (zh) | 磁共振化学位移编码成像方法与装置 | |
Datta et al. | A deep framework for enhancement of diagnostic information in CSMRI reconstruction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20170321 Address after: 200241 Shanghai, Minhang District purple Road, No. 119, room 117, room 1000 Patentee after: Shanghai Neusoft Medical Technology Co., Ltd. Address before: Hunnan New Century Road 110179 Shenyang city of Liaoning Province, No. 16 Patentee before: Dongruan Medical Systems Co., Ltd., Shenyang |