CN1357302A - 图像处理方法和设备、记录介质和成像设备 - Google Patents
图像处理方法和设备、记录介质和成像设备 Download PDFInfo
- Publication number
- CN1357302A CN1357302A CN01143519A CN01143519A CN1357302A CN 1357302 A CN1357302 A CN 1357302A CN 01143519 A CN01143519 A CN 01143519A CN 01143519 A CN01143519 A CN 01143519A CN 1357302 A CN1357302 A CN 1357302A
- Authority
- CN
- China
- Prior art keywords
- pixel value
- interest
- variance
- pixels
- pixel
- 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
Links
- 238000003672 processing method Methods 0.000 title claims description 25
- 238000003384 imaging method Methods 0.000 claims description 29
- 230000005764 inhibitory process Effects 0.000 claims description 6
- 230000002708 enhancing effect Effects 0.000 abstract 1
- 230000005291 magnetic effect Effects 0.000 description 50
- 230000006870 function Effects 0.000 description 15
- 238000000034 method Methods 0.000 description 14
- 230000017531 blood circulation Effects 0.000 description 13
- 230000003068 static effect Effects 0.000 description 9
- 230000000694 effects Effects 0.000 description 8
- 238000013480 data collection Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 239000004744 fabric Substances 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 238000005481 NMR spectroscopy Methods 0.000 description 2
- 210000003050 axon Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 210000000038 chest Anatomy 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000002600 positron emission tomography Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000000007 visual 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/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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
Landscapes
- Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Heart & Thoracic Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- High Energy & Nuclear Physics (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Radiology & Medical Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Image Processing (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Image Analysis (AREA)
- Facsimile Image Signal Circuits (AREA)
Abstract
本发明涉及图像处理方法和装置以及记录介质和成像装置。为了增大图像中涉及物质结构的像素与其他像素之间的差异,把构成图像的每个像素定义为感兴趣的像素,在所述感兴趣的像素所属的局部区域中确定像素值的方差(508);并且当所述像素值的方差明显大于噪声方差时,维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值(510,512)。
Description
技术领域 本发明涉及图像处理方法和设备、记录介质和成像设备,具体地说,本发明涉及用于调整构成图像的像素值的图像处理方法和设备、带有使计算机实现这种图像处理功能的程序的记录介质、以及包括这种图像处理设备的成像设备。
背景技术 在磁共振成像(MRI)设备中,把要成像的对象载入磁系统的内部空间中,即这样的空间:其中产生静态磁场;施加梯度磁场和高频磁场而导致对象内部的自旋、从而产生磁共振信号;并且根据所接收的信号产生层析图像。
在身体内部移动的自旋、诸如血流内的自旋与不移动的自旋、如组织中的自旋之间,梯度磁场和高频磁场对自旋的影响是不同的。利用这种差异,可以捕获身体内部移动的自旋的图像,即例如血流图像。
在捕获血流图像时,采用了飞行时间(TOF)技术、相位对比(PC)技术等等。
通过利用这些技术之一来捕获关于三维区域的多切片血流层析图像、并且在切片厚度方向上对多切片血流层析图像进行最大强度投影(MIP),获得三维区域中的血流投影图像。
当如上所述获得血流的投影图像时,模糊的血流无法被投影,因为噪声使其变模糊。另外,当各切片之间、图像的平均信号强度不同时,在平均信号强度小的图像中的血流图像无法被投影,因为平均信号强度大的图像中的噪声使其变得不明显。
发明内容 因此,本发明的一个目的是提供一种提高图像中有关物质结构的像素与其他像素之间差异的图像处理方法和设备、记录有使计算机实现这种图像处理功能的程序的介质、以及包括这种图像处理设备的成像设备。
(1)在解决上述问题的一个方面中,本发明是一种图像处理方法,其特征在于:把构成某个图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;并且当所述确定的方差明显大于噪声方差时,维持所述感兴趣的像素的像素值,否则就抑制所述感兴趣的像素的像素值。
(2)在解决上述问题的另一个方面中,本发明是一种图像处理设备,其特征在于包括:方差计算装置,用于把构成某个图像的每个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中的像素值的方差;以及像素值调整装置,用于当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值。
(3)在解决上述问题的又一方面中,本发明是记录介质,其特征在于以计算机可读的方式记录了使计算机实现以下功能的程序:把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;并且当所述确定的方差明显大于噪声方差时,维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值。
(4)在解决上述问题的再一个方面中,本发明是用于根据从对象收集的信号产生图像的成像设备,其特征在于包括:方差计算装置,用于把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;以及像素值调整装置,用于当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值。
按照本发明、在以上(1)-(4)中描述的各方面中,把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;并且当所确定的方差明显大于噪声方差时,维持感兴趣的像素的像素值,否则抑制感兴趣的像素的像素值;因此,在图像中有关物质结构的像素与其他像素之间的差异得以增大。
(5)在解决上述问题的又一方面中,本发明是图像处理方法,其特征在于:把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;并且当所述确定的方差明显大于噪声方差时,增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值。
(6)在解决上述问题的再一方面中,本发明是图像处理设备,其特征在于包括:方差计算装置,用于把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;以及像素值调整装置,用于当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值、否则维持所述感兴趣的像素的像素值。
(7)在解决上述问题的又一方面中,本发明是记录介质,其特征在于以计算机可读的方式记录了使计算机实现以下功能的程序:把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;并且当所述确定的方差明显大于噪声方差时,增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值。
(8)在解决上述问题的再一方面中,本发明是用来根据从对象收集的信号产生图像的成像设备,其特征在于包括:方差计算装置,用于把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;以及像素值调整装置,用于当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值、否则维持所述感兴趣的像素的像素值。
按照本发明、在以上(5)-(8)中描述的各方面中,把构成图像的每个像素定义为感兴趣的像素,在感兴趣的像素所属的局部区域中确定像素值的方差,并且当确定的方差明显大于噪声方差时,增大感兴趣的像素的像素值,否则维持感兴趣的像素的像素值;因此,能增大涉及图像中物质结构的像素与其他像素之间的差异。
另外,在以上(1)-(8)中描述的各方面中,图像可以是血流图像,由此增大涉及血流图像的像素与其他像素之间的差异。
另外,按照本发明,在以上(4)-(8)中描述的各方面中,信号可以是磁共振信号,由此实现用于磁共振图像的发明。
(9)在解决上述问题的又一方面中,本发明是图像处理方法,其特征在于:把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属局部区域中的像素值的方差;进行像素值调整,包括当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;并且对受到所述像素值调整的多切片图像进行最大强度投影。
(10)在解决上述问题的再一个方面中,本发明是图像处理设备,其特征在于包括:方差计算装置,用于把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;像素值调整装置,用于当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及最大强度投影装置,用于对受到所述像素值调整的多切片图像进行最大强度投影。
(11)在解决上述问题的再一个方面中,本发明是记录介质,其特征在于以计算机可读的方式记录使计算机实现以下功能的程序:利把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;进行像素值调整,包括当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及对受到所述像素值调整的多切片图像进行最大强度投影。
(12)在解决上述问题的又一个方面中,本发明是用来根据从对象收集的信号产生图像的成像设备,其特征在于包括:方差计算装置,用于把构成多切片图像的每个像素定义为感兴趣的像素、确定所述感兴趣的像素所属局部区域中的像素值的方差;像素值调整装置,用于当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及最大强度投影装置,用于对受到所述像素值调整的多切片图像进行最大强度投影。
按照本发明在上面(9)-(12)中描述的各方面中,把构成多切片图像的每个像素定义为感兴趣的像素,在所述感兴趣的像素所属局部区域中确定像素值的方差,并且当确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;因此,可增大涉及图像中物质结构的像素与其他像素之间的差异。另外,由于对受到这种差异增大的多切片图像进行最大强度投影,可获得图像中具有微弱信号强度并且涉及物质结构的像素的投影图像。
(13)在解决上述问题的再一个方面中,本发明是图像处理方法,其特征在于:把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属局部区域中的像素值的方差;执行像素值调整,包括当所述确定的方差明显大于噪声方差时,增大所述感兴趣像素的像素值,否则维持所述感兴趣像素的像素值;并对受到所述像素值调整的多切片图像进行最大强度投影。
(14)在解决上述问题的又一个方面中,本发明是图像处理设备,其特征在于包括:方差计算装置,用于把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣像素所属的局部区域中的像素值的方差;像素值调整装置,用于当所述确定的方差明显大于噪声方差时、增大所述感兴趣像素的像素值,否则维持所述感兴趣像素的像素值;以及最大强度投影装置,用于对受到所述像素值调整的多切片图像进行最大强度投影。
(15)在解决上述问题的再一个方面中,本发明是记录介质,其特征在于以计算机可读的方式记录有使计算机实现以下功能的程序:把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣像素所属的局部区域中像素值的方差;进行像素值调整,包括当所述确定的方差明显大于噪声方差时、增大所述感兴趣像素的像素值,否则维持所述感兴趣像素的像素值;以及对受到所述像素值调整的多切片图像进行最大强度投影。
(16)在解决上述问题的再一个方面中,本发明是用来根据从对象收集的信号产生图像的成像设备,其特征在于包括:方差计算装置,用于把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣像素所属于的局部区域中的像素值的方差;像素值调整装置,用于当所述确定的方差明显大于噪声方差时、增大所述感兴趣像素的像素值,否则维持所述感兴趣像素的像素值;以及最大强度投影装置,用来对受到所述像素值调整的多切片图像进行最大强度投影。
按照本发明、在上面(13)-(16)中描述的各方面中,把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣像素所属的局部区域中的像素值的方差,并且当确定的方差明显大于噪声方差时,增大感兴趣像素的像素值,否则维持感兴趣像素的像素值;因此,可增大涉及图像中物质结构的像素和其他像素之间的差异。另外,由于对受到这种差异增大的多切片图像进行最大强度投影,因此可获得图像中的具有微弱信号强度并涉及物质结构的像素的投影图像。
另外,按照本发明在以上(9)-(16)中描述的各方面中,图像可以是血流图像,从而增大涉及血流图像的像素与其他像素之间的差异;而且通过对受到这样的差异增大的多切片图像进行最大强度投影,可获得微弱血流的投影图像。
此外,按照本发明在以上(12)和(16)中描述的各方面中,所述信号可以是磁共振信号,从而实现用于磁共振图像的发明。
(17)在解决上述问题的再一个方面中,本发明是图像处理方法,其特征在于:把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣像素所属局部区域中的像素值的方差;在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值方差;进行像素值调整,包括当所述相加后的方差明显大于噪声方差时,维持所述感兴趣像素的像素值,否则抑制所述感兴趣像素的像素值;以及对受到所述像素值调整的多切片图像进行最大强度投影。
(18)在解决上述问题的再一个方面中,本发明是图像处理设备,其特征在于包括:方差计算装置,用来把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣像素所属的局部区域中的像素值的方差;加法装置,用于在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值的方差;像素值调整装置,用于当所述相加后的方差明显大于噪声方差时、维持所述感兴趣像素的像素值,否则抑制所述感兴趣像素的像素值;以及最大强度投影装置,用于对受到所述像素值调整的多切片图像进行最大强度投影。
(19)在解决上述问题的再一个方面中,本发明是记录介质,其特征在于以计算机可读的方式记录有使计算机实现以下功能的程序:把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属局部区域中的像素值的方差;在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值的方差;进行像素值调整,包括当所述相加后的方差明显大于噪声方差时,维持所述感兴趣像素的像素值,否则抑制所述感兴趣像素的像素值;以及对受到所述像素值调整的多切片图像进行最大强度投影。
(20)在解决上述问题的再一个方面中,本发明是用来根据从对象收集的信号产生图像的成像设备,其特征在于包括:方差计算装置,用于把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣像素所属的局部区域中的像素值的方差;加法装置,用于在所述确定的方差上加上相邻切片图像中相应的感兴趣像素所属局部区域中的像素值的方差;像素值调整装置,用于当所述相加后的方差明显大于噪声方差时、维持所述感兴趣像素的像素值,否则抑制所述感兴趣像素的像素值;以及最大强度投影装置,用于对受到所述像素值调整的多切片图像进行最大强度投影。
按照本发明在上面(17)-(20)描述的各方面中,把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣像素所属的局部区域中的像素值的方差;在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中像素值的方差,并且当所述相加的方差明显大于噪声方差时、维持感兴趣的像素的像素值,否则抑制感兴趣像素的像素值;因此,可进一步增大涉及图像中物质结构的像素与其他像素之间的差异。另外,由于对受到这种差异增强的多切片图像进行最大强度投影,所以可获得在图像中具有微弱信号强度并涉及物质结构的像素的更佳投影图像。
(21)在解决上述问题的再一个方面中,本发明是图像处理方法,其特征在于:把构成多切片图像的每一个像素定义为感兴趣的像素,确定所述感兴趣像素所属的局部区域中的像素值的方差;在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值的方差;进行像素值调整,包括当所述相加后的方差明显大于噪声方差时,增大所述感兴趣像素的像素值,否则维持所述感兴趣像素的像素值;以及对受到所述像素值调整的多切片图像进行最大强度投影。
(22)在解决上述问题的再一个方面中,本发明是图像处理设备,其特征在于包括:方差计算装置,用于把构成多切片图像的每一个像素定义为感兴趣的像素,确定所述感兴趣像素所属的局部区域中的像素值的方差;加法装置,用于在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值的方差;像素值调整装置,用于当所述相加后的方差明显大于噪声方差时、增大所述感兴趣像素的像素值,否则维持所述感兴趣像素的像素值;以及最大强度投影装置,用于对受到所述像素值调整的多切片图像进行最大强度投影。
(23)在解决上述问题的再一个方面中,本发明是记录介质,其特征在于以计算机可读方式记录有使计算机实现以下功能的程序:把构成多切片图像的每一个像素定义为感兴趣的像素,确定所述感兴趣的像素所属局部区域中的像素值的方差;在所述确定的方差上加上相邻切片的图像中相应感兴趣像素所属的局部区域中像素值的方差;进行像素值调整,包括当所述相加后的方差明显大于噪声方差时、增大所述感兴趣像素的像素值,否则维持所述感兴趣像素的像素值;以及对受到所述像素值调整的多切片图像进行最大强度投影。
(24)在解决上述问题的又一个方面中,本发明是用来根据从对象收集的信号产生图像的成像设备,其特征在于包括:方差计算装置,用来把构成多切片图像的每一个像素定义为感兴趣的像素、确定所述感兴趣像素所属的局部区域中像素值的方差;加法装置,用来在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属局部区域中的像素值的方差;像素值调整装置,用于当所述相加后的方差明显大于噪声方差时、增大所述感兴趣像素的像素值,否则维持所述感兴趣像素的像素值;以及最大强度投影装置,用来对受到所述像素值调整的多切片图像进行最大强度投影。
按照本发明在上面(21)-(24)描述的各方面中,把构成多切片图像的每个像素定义为感兴趣的像素、确定所述感兴趣的像素所属局部区域中的像素值的方差,在所述确定的方差上加上相邻切片的图像中相应感兴趣的像素所属局部区域中的像素值的方差,并且当相加后的值明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣像素的像素值;因此,可以进一步增强涉及图像中物质结构的像素与其他像素之间的差异。另外,由于对受到这种差异增强的多切片图像进行最大强度投影,所以可获得图像中具有微弱信号强度并涉及物质结构的像素的更佳投影图像。
另外,按照本发明在以上(17)-(24)中描述的各方面中,图像可以是血流图像,从而进一步增强涉及血流图像的像素与其他像素之间的差异;而且通过对受到这样的差异增强的多切片图像进行最大强度投影,可获得微弱血流的更佳投影图像。
此外,按照本发明在如以上(20)和(24)中描述的各方面中,所述信号可以是磁共振信号,从而实现用于磁共振图像的发明。
在本发明中在以上(1)-(24)中描述的各方面中,最好对在整个图像上定义的多个局部区域中的每个、确定像素值的剩余平方和,确定剩余平方和的直方图,然后根据给出直方图峰值的剩余平方和计算噪声方差,其中可直接根据图像获得噪声方差。
因此,本发明可以提供增大图像中与物质结构相关的像素与其他像素之间差异的图像处理方法和设备、记录了可使计算机完成这样的图像处理功能的程序的介质、以及包括这样的图像处理设备的成像设备。
从如附图中举例说明的本发明最佳实施例的以下描述,本发明的其他目标和优点将显而易见。
附图说明
图1是按照本发明的一个实施例的设备的框图;
图2是按照本发明的一个实施例的设备的框图;
图3是图1或2中所示设备的工作流程图;
图4是多切片图像的概念示意图;
图5是图3中所示流程图的一部分的详细流程图;
图6是图5中所示流程图的一部分的详细流程图;
图7是直方图的概念示意图;
图8是直方图的概念示意图;
图9是说明感兴趣的像素与局部区域之间关系的示意图;
图10是作为图5中所示流程图的一部分而插入的程序的流程图;
图11是表示像素值调整效果的图像曲线的实例;
图12是表示像素值调整效果的图像曲线的实例;
图13是图3中所示流程图的一部分的详细流程图;
图14是表示像素值调整效果的图像曲线的实例;
图15是表示像素值调整效果的图像曲线的实例;
图16是最大强度投影的概念图。
具体实施方式 现参照附图详细描述本发明的实施例。应该指出,本发明不限于这些实施例。图1表示成像设备或磁共振成像(MRI)设备的框图,这是本发明的一个实施例。该设备的配置代表按照本发明的设备的实施例。该设备的工作代表按照本发明的方法的实施例。
如图1中所示,本设备具有磁系统100。磁系统100具有主磁场线圈部分102、梯度线圈部分106和RF(射频)线圈部分108。这些线圈部分一般为圆筒形并且同轴设置。要成像的对象300躺在托架500上,并且通过运载装置(未示出)运进和运出磁系统100的一般为圆筒形的内部空间(膛)。
主磁场线圈部分102在磁系统100的内部空间中产生静态磁场。静态磁场的方向一般与对象300的体轴的方向平行。就是说,产生“水平”磁场。主磁场线圈部分102是用例如超导线圈制成的。容易理解,主磁场线圈部分102不限于超导线圈,而是可以用普通导电线圈等制成。
梯度线圈部分106产生梯度磁场,以便把梯度赋予静态磁场强度。要产生的梯度磁场为如下三种:切片梯度磁场、读出梯度磁场和相位编码梯度磁场。梯度线圈部分106具有未示出的三种梯度线圈,对应于这三种梯度磁场。
RF线圈部分108产生高频磁场、用于在静态磁场空间中激发对象300内的自旋。在下文中,产生高频磁场有时称为发射RF激发信号。RF线圈部分108还接收电磁波、即由激发的自旋产生的磁共振信号。
RF线圈部分108具有未示出的发射和接收线圈。对于发射和接收线圈,可以采用同一线圈或分开的专用线圈。
梯度线圈部分106与梯度驱动部分130连接。梯度驱动部分130向梯度线圈部分106提供驱动信号、以便产生梯度磁场。梯度驱动部分130具有对应于梯度线圈部分106中三种梯度线圈的三种未示出的驱动电路。
RF线圈部分108与RF驱动部分140连接。RF驱动部分140向RF线圈部分108提供驱动信号、以便发出RF激发信号,以此激发对象300内的自旋。
RF线圈部分108与数据收集部分150连接,数据收集部分150收集由RF线圈部分108接收的接收信号,并且收集这些信号作为图像数据。
梯度驱动部分130、RF驱动部分140和数据收集部分150与控制部分160连接。控制部分160控制梯度驱动部分130、RF驱动部分140和数据收集部分150以实现成像。
把数据收集部分150的输出连接到数据处理部分170。数据处理部分170是例如用计算机构成的。数据处理部分170具有未示出的存储器。存储器存储数据处理部分170所用的程序和几种类型的数据。通过数据处理部分170执行存储器中存储的程序来实现本设备的功能。
数据处理部分170把从数据收集部分150收集的图像数据存入存储器。在存储器中形成数据空间。该数据空间构成二维傅里叶空间。二维傅里叶空间有时称为k空间。数据处理部分170对二维傅里叶空间中的数据进行二维傅里叶逆变换,从而产生(重构)对象300的图像。
由二维傅里叶逆变换重构的图像具有复数的像素值。该复数的绝对值用来构造绝对值图像。复数的实部可用来构造实部图像。复数的虚部可用来构造虚部图像。虚部和实部都可以是正值和负值。这样的图像有时称为正负图像。
数据处理部分170具有为确定关于重构图像的像素值的方差而进行图像处理的功能。数据处理部分170还具有为确定关于重构图像的噪声方差而进行图像处理的功能。数据处理部分170还具有为调整关于重构图像的像素值而进行图像处理的功能。数据处理部分170还具有进行图像处理、以便对受到像素值调整的图像进行最大强度投影(MIP)的功能。下面将描述数据处理部分170的这些图像处理功能。
数据处理部分170是本发明的图像处理设备的一个实施例。该设备的配置代表按照本发明的设备的实施例。该设备的工作代表按照本发明的方法的实施例。
数据处理部分170连接到控制部分160。数据处理部分170高于控制部分160并且控制后者。数据处理部分170与显示部分180和操作部分190相连。显示部分180包括图形显示器等。操作部分190包括键盘等,并且设有定位设备。
显示部分180显示从数据处理部分170输出的重构图像和几类信息。操作部分190由操作者操作,并且该部分190向数据处理部分170输入几种命令、信息等等。操作者经由显示部分180和操作部分190交互地操作本设备。
图2是作为本发明的一个实施例的另一种类型的MRI设备的框图。该设备的配置代表按照本发明的设备的实施例。
图2中所示设备具有与图1中所示设备不同类型的磁系统100’。由于除磁系统100’以外,该设备具有类似于图1所示设备的配置,所以类似部分由类似的参考标号标明并且省略其说明。
磁系统100’具有主磁场磁体部分102’、梯度线圈部分106’和RF线圈部分108’。主磁场磁体部分102’和线圈部分均由跨过一定空间彼此面对的一对组件组成。这些部分一般具有盘状外形并且布置成具有共同的中心轴。对象300躺在托架500上,并且由未示出的运载装置运入和运出磁系统100’的内部空间(膛)。
主磁场磁体部分102’在磁系统100’的内部空间中产生静态磁场。静态磁场的方向一般与对象300的体轴方向垂直。也就是说,产生“垂直磁场”。主磁场磁体部分102’是用例如永磁体制成的。容易认识到,主磁场磁体部分102’不限于永磁体,而是可用超导或普通导电的电磁体等来制造。
梯度线圈部分106’产生梯度磁场,以便把梯度赋予静态磁场强度。要产生的梯度磁场是下列三种:切片梯度磁场、读出梯度磁场和相位编码梯度磁场。梯度线圈部分106’具有与这三种梯度磁场对应的三种未示出的梯度线圈。
RF线圈部分108’发出RF激发信号、以便在静态磁场空间中激发对象300内的自旋。RF线圈部分108’还接收由激发的自旋产生的磁共振信号。RF线圈部分108’具有未示出的发射和接收线圈。对于发射和接收线圈,可以使用同一线圈或分开的专用线圈。
图3表示本设备的工作流程图。图1和2中所示的设备都以同样的方式工作。如图3中所示,在步骤302进行血流成像。为了血流成像,采用飞行时间(TOF)技术、相位对比(PC)技术等。此外,在多切片中进行成像。这样,如图4中从概念上说明的,就对象300的三维区域来捕获多切片血流层析图像S1、S2、S3、…、Sm。
接着,在步骤304,对血流层析图像S1、S2、S3、…、Sm进行像素值调整。像素值调整是由数据处理部分170的数据处理功能来实现的。以下把血流层析图像简称为图像。
图5表示像素值调整的详细流程图。如所示,在步骤502进行切片选择。于是,选择图像S1、S2、S3、…、Sm中的一个,例如,图像S1。
接着,在步骤504,进行噪声方差的计算。在步骤504计算噪声方差的数据处理部分170是本发明的噪声方差计算装置的一个实施例。图6表示噪声方差计算的详细流程图。如图所示,在步骤602在图像中定义一个局部区域。局部区域是供下一个步骤中计算用的像素值所属的区域。例如,图像中心处的局部区域被定义为第一区域。
作为局部区域,采用N×N像素矩阵。例如,N为9。应该指出,矩阵大小不限于此值,任何适当的数值均可。此外,像素矩阵不限于方阵,以某个像素为中心的任何适当区域均可。在下文中局部区域有时简称为区域。
接着,在步骤604,确定属于该区域的像素值的剩余平方和S,具体为: 式中:Pi是像素值,而
Pi是以Pi为中心的N×N区域中的像素值的平均值。另外,k是例如81。
接着,在步骤606,判断是否对所有局部区域已完成上述处理,如果不是,则在步骤608移动局部区域。这样,例如,选择相邻的N×N区域作为新的局部区域。
在新的局部区域上执行步骤604的处理,以便确定像素值的剩余平方和。此后,以类似方式对该图像中的每个局部区域确定像素值的剩余平方和。
这样得出的剩余平方和具有χ2分布,而且其平均值为k·σ2。若k大,则χ2分布趋近于高斯分布,且其峰值的位置大约在k·σ2处。
接着,在步骤610,产生剩余平方和S的直方图。
图7表示图像是绝对值图像时剩余平方和S的直方图的概念。如图所示,直方图包括三条分布曲线a、b和c。
分布曲线a是高斯分布曲线,是均匀结构部分中的噪声造成的。分布曲线b是瑞利分布曲线,是不包含对象300的FOV(视场)部分中的噪声、即背景中的噪声造成的。因为图像是绝对值图像,所以背景中的噪声造成的分布曲线与高斯分布不一致,而是符合瑞利分布。分布曲线c是对象的精细结构造成的,并且呈现出不同于其他两条曲线的不确定分布。
在步骤612,对这样的直方图进行峰值位置检测。于是,针对高斯分布曲线a检测峰值位置s1,并且针对瑞利分布曲线b检测峰值位置s2。
因为直方图实际上具有离散的值,最好在峰值检测之前在步骤612执行拟合到某个函数,以此可以高精度地检测峰值位置。拟合中用的函数分别是例如高斯分布函数和瑞利分布函数。但是,不限于这些函数,任何其他适当的函数均可。
接着,在步骤614,计算噪声方差。噪声方差的计算是基于峰值位置s1或s2来进行的。
因为s1、s2和σ分别具有以下关系:
S1=k·σ2 (2)
以及
S2=(2-π/2)k·σ2 (3)从这些关系中确定σ值。不论是由方程式(2)还是由方程式(3)确定σ,σ值都相同。确定的σ值作为噪声方差Vn存储在存储器中。
在分布曲线c的某些情况下,高斯分布曲线a的峰值位置s1可能无法精确检测。在这种情况下,σ值是根据瑞利分布曲线b的峰值位置s2确定的。另外,对于具有较大比例的背景部分区域的图像,瑞利分布曲线b对于高精度地确定噪声方差更为合适。
虽然以上描述是针对绝对值图像的情况而作的,但是当要处理的图像是正负像、即实部图像或虚部图像时,背景部分中的噪声具有以零为中心的正值和负值。
因此,在步骤610产生的直方图变成如图8中举例示出的那样,而它不再具有瑞利分布。在这种情况下,在步骤614根据高斯分布曲线a的峰值位置s1确定噪声方差。
于是,可以直接根据实际捕获的图像求得噪声方差的值。若噪声方差是以前已知的,则可使用该方差并且省略所述计算步骤。
这样确定噪声方差Vn之后,在图5的流程图中,在步骤506在图像中定义感兴趣的像素。第一个感兴趣的像素是例如在图像中心的某个像素。
接着,在步骤508,计算包含该感兴趣像素的局部区域中的像素值的方差Vi。包含该感兴趣像素的局部区域是、例如、如图9中所示的以感兴趣的像素i为中心的5×5矩阵。应该指出,矩阵的大小不限于此值,任何适当的值均可。另外,像素矩阵不限于方形矩阵,任何以某像素为中心的适当区域均可。下文中有时把局部区域简称为区域。在步骤508计算像素值的方差Vi的数据处理部分170是本发明的方差计算装置的实施例。
以下的方程式用来计算像素值的方差Vi:
式中k=25。
接着,在步骤510,判断像素值的方差Vi是否明显大于噪声方差Vn。利用下列公式做出该判断: 式中:
γ:阈值。
对于阈值γ,可以采用任何大于1的适当值。
若包含感兴趣像素的局部区域中的像素值的方差并非明显大于噪声方差,则局部区域中的图像或许没有突出的结构,像素值的方差可能来源于噪声。
因此,在这种情况下,在步骤512抑制所述感兴趣像素的像素值。例如,通过把该像素值乘上某个系数α来达到像素值的抑制。系数α的值是小于1是正数,例如0.8。于是,感兴趣的像素的像素值被减小为例如初始值的0.8倍。但是,系数α的值不限于0.8,任何适当的数值均可。另外,可以例如通过从像素值中减去某个预定的值来达到像素值的抑制。应该指出,所述的常数值不得超过像素值的最小值。
若包含感兴趣的像素的局部区域中的像素值方差Vi明显大于噪声方差Vn,则局部区域中的图像可能具有特定的结构、如边缘,而像素值的方差大概来源于该图像的结构。在这种情况下,对该像素值不进行特殊操作。于是,所述感兴趣的像素的像素值维持其原值。执行这种像素值调整的数据处理部分170是本发明的像素值调整装置的实施例。
接着,在步骤514,判断是否对于所有感兴趣的像素均已完成上述处理,若未完成,则在步骤516把感兴趣的像素移到例如相邻的一个,并且执行从步骤508开始的处理。此后,重复同样的处理,直至为图像S1中每一个像素调整了像素值。
然后,在步骤518,判断是否对所有切片均已完成上述处理,若未完成,则在步骤520转移切片,在该切片的图像上进行同样的处理。此后,重复同样的处理,从而对所有图像S1-Sm中的像素都进行了像素值调整。
在步骤508和510之间,可以加上图10的流程图中所示的步骤。具体地说,在步骤702,对相邻切片中包含相应感兴趣的像素的局部区域计算像素值的方差Vi’。
术语“相邻切片”意味着与在步骤508已确定像素值方差Vi的切片相邻的一个或多个切片。对于这类切片,例如,可采用贴近前面或后面的一个或几个切片。
在步骤704,把像素值的方差Vi’加到Vi上,并且把相加后的值定义为供下一个步骤510判断使用的像素值方差Vi。在加法中可以对Vi’加上适当的权重。在步骤702中计算像素值方差的数据处理部分170是本发明的像素值方差计算装置的一个实施例。在步骤704中把像素值方差相加的数据处理部分170是本发明的相加装置的一个实施例。
于是,在通过上述处理获得的像素值方差Vi中,反映了跨多个切片的结构。因此,例如,若血流图像在穿过各切片的方向上存在,该图像呈现为一个图像上的一点,则可获得准确反映这种结构的像素值方差,根据该方差可以得到更准确的像素值调整。
图11把这种像素值调整的效果表示成像素值分布的变化。图11中符号B表示像素值调整之前的分布,在背景噪声上存在明显的血流图像b1和模糊的血流图像b2。
作为上述像素值调整的结果,这样的分布具有被抑制了0.8倍的背景噪声的像素值,而同时维持血流图像b1和b2的像素值不变,结果得到如图11中A处所示的分布。在分布A中,原始图像中模糊的血流图像b2呈现出与背景噪声之间放大的差异,从而变得明显,更不用说血流图像b1了。因而,可增强在原始图像中模糊的血流图像b2的可辨认性。
图12表示像素值调整的另一种效果。图12中的符号P和Q表示不同切片的两幅图像的分布,分布Q的背景噪声电平大于分布P中明显的血流图像b1的信号强度。
通过上述像素值调整,这样的分布具有被抑制了0.8倍的背景噪声的像素值,因而如图12中Q’所示,可以获得相对于血流图像b1和b2的信号强度、噪声电平减小的分布。因而,血流图像b1和b2与切片Q的图像的噪声电平之间的差异也变明显了,并且两幅图像都可辨认。
图13表示像素值调整的另一种程序的流程图。在图13中,与图5中所示的类似的步骤标以类似的参考标号并省略其说明。图5和13中所示程序之间的差异在于步骤510中判断之后的像素值处理。
具体地讲,若在包含感兴趣的像素的局部区域中、像素值方差明显大于噪声方差,则在步骤512’增大感兴趣的像素的像素值。像素值的增大是通过例如把像素值乘以系数β来达到的。系数β的值是大于1的正数、例如1.2。于是,感兴趣的像素的像素值增大为原值的例如1.2倍。应该指出,系数β的值不限于1.2,它可以是任何适当的值。另外,不是乘以某个系数,而是例如通过给像素值加上某个预定值,亦可得到像素值的增大。
若在包含感兴趣的像素的局部区域中、像素值的方差Vi并非明显大于噪声方差Vn,则对该像素值不作特殊处理。因而,感兴趣的像素的像素值维持其原值。执行这种像素值调整的数据处理部分170是本发明的像素值调整装置的一个实施例。
图14通过像素值分布的变化来表示这种像素值调整的效果。如图所示,在像素值调整之前的分布中血流图像b1和b2具有由于上述像素值调整而增大的像素值,如血流图像b1’和b2’所示。于是,增大了与背景噪声的差异并且可辨认性增强了。
图15表示像素值调整的另一种效果。图15中的符号P和Q表示不同切片的两幅图像的分布。即使当分布Q的背景噪声电平大于分布P中明显的血流图像b1的信号强度时,通过上述像素值调整,分布P中的血流图像b1和b2的像素值增大了例如1.2倍,结果得出血流图像b1’和b2’。这样,相对于切片Q的图像的噪声电平,血流图像b1’和b2’也可辨认。
对于如上所述的像素调整之后的多切片图像,在图3的流程图中,在步骤306进行最大强度投影。在步骤306执行最大强度投影的数据处理部分170是本发明的最大强度投影装置的一个实施例。
以计算机可读的方式在记录介质上记录使计算机实现上述功能的程序。对于所述记录介质,采用例如磁记录介质、光记录介质、磁光记录介质以及任何其他适当类型的记录介质中的任何一种。记录介质可以是半导体存储介质。在本说明中存储介质与记录介质同义。
图16表示最大强度投影的概念示意图。如图所示,沿着穿过多切片图像S1-Sm的视线E提取像素值的最大值,并且用所提取的值作为投影图像R的像素值。用等于投影图像R中的像素个数的若干条视线作为视线E。
按照如上所述的像素值调整,因为对于每幅图像S1-Sm,已经增大了血流图像与噪声之间的差异,即使模糊的血流图像也能明显地呈现出来而不被噪声所掩盖。因此,即使血流图像模糊,也能获得具有清晰的血流图像的MIP图像。在步骤308,在显示部分180上显示这样的MIP图像。
已经就磁共振成像设备中数据处理部分执行图像处理的实例做了上述描述;但是,容易理解,所述图像处理可由与磁共振成像设备分开的数据处理设备、诸如EWS(工程工作站)或PC(个人计算机)来执行。
此外,尽管所述成像设备是作为MRI设备进行描述的,但是成像设备不限于此,而是可以为任何其他类型的成像设备,诸如X射线CT(计算层析成像)设备、X射线成像设备、PET(正电子放射断层照相法)或γ摄像机。
另外,尽管参照处理医学图像的实例做了描述,但是,要处理的对象不限于医学图像,而是可以对各种各样的图像、如由光学仪器捕获的数字图像进行图像处理。
尽管已参照上述最佳实施例描述了本发明,但是只要不背离本发明的范围,在与本发明相关的领域的技术人员可对这些实施例作出各种变化和替换。因此,本发明的技术范围不仅包括上述那些实施例,而且包括属于所附权利要求书的范围之内的所有实施例。
Claims (36)
1.一种图像处理方法,它包括以下步骤:
把构成某个图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属于的局部区域中的像素值的方差;以及
当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值。
2.一种图像处理方法,它包括以下步骤:
把构成多切片图像的每一个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;
进行像素值调整,包括当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及
对受到所述像素值调整的所述多切片图像进行最大强度投影。
3.一种图像处理方法,它包括以下步骤:
把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;
在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值的方差;
进行像素值调整,包括当所述相加后的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及
对受到所述像素值调整的所述多切片图像进行最大强度投影。
4.按照权利要求1-3中任何一个的图像处理方法,其特征在于:所述像素值的抑制是通过乘以一个小于1的系数来完成的。
5.按照权利要求1-3中任何一个的图像处理方法,其特征在于:所述像素值的抑制是通过减去一个预定的数值来完成的。
6.一种图像处理方法,它包括以下步骤:
把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;以及
当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值。
7.一种图像处理方法,它包括以下步骤:
把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;
进行像素值调整,包括当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值;以及
对受到所述像素值调整的所述多切片图像进行最大强度投影。
8.一种图像处理方法,它包括以下步骤:
把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;
在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值的方差;
进行像素值调整,包括当所述相加后的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值;以及
对受到所述像素值调整的所述多切片图像进行最大强度投影。
9.按照权利要求6-8中任何一个的图像处理方法,其特征在于:所述像素值的增大是通过乘以某个大于或等于1的系数来完成的。
10.按照权利要求6-8中任何一个的图像处理方法,其特征在于:所述像素值的增大是通过加上预定的数值来完成的。
11.按照权利要求1-3和6-8中任何一个的图像处理方法,其特征在于还包括以下步骤:对遍及整个图像定义的多个局部区域中的每一个确定像素值的剩余平方和;确定所述剩余平方和的直方图;然后根据给定所述直方图峰值的剩余平方和来确定所述噪声方差。
12.按照权利要求1-3和6-8中任何一项的图像处理方法,其特征在于:所述图像是血流图像。
13.一种图像处理设备,它包括:
方差计算装置,用于把构成图像的每个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中的像素值的方差;以及
像素值调整装置,用于当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值。
14.一种图像处理设备,它包括:
方差计算装置,用于把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;
像素值调整装置,用于当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及
最大强度投影装置,用于对受到所述像素值调整的所述多切片图像进行最大强度投影。
15.一种图像处理设备,它包括:
方差计算装置,用于把构成多切片图像的每个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中的像素值的方差;
加法装置,用于在所述确定的方差上加上相邻切片的图像中相应的感兴趣的像素所属的局部区域中的像素值的方差;
像素值调整装置,用于当所述相加后的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及
最大强度投影装置,用于对受到所述像素值调整的所述多切片图像执行最大强度投影。
16.按照权利要求13-15中任何一个的图像处理设备,其特征在于:所述像素值调整装置通过乘以一个小于1的系数来执行所述像素值的抑制。
17.按照权利要求13-15中任何一个的图像处理方法,其特征在于:所述像素值调整装置通过减去预定的数值来执行所述像素值的抑制。
18.一种图像处理设备,它包括:
方差计算装置,用于把构成图像的每个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中的像素值的方差;以及
像素值调整装置,用于当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值。
19.一种图像处理设备,它包括:
方差计算装置,用于把构成多切片图像的每一个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中的像素值的方差;
像素值调整装置,用于当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值;以及
最大强度投影装置,用于对受到所述像素值调整的多切片图像进行最大强度投影。
20.一种图像处理设备,它包括:
方差计算装置,用于把构成多切片图像的每一个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中的像素值的方差;
加法装置,用于在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值的方差;
像素值调整装置,用于当所述相加后的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值;以及
最大强度投影装置,用来对受到所述像素值调整的所述多切片图像进行最大强度投影。
21.按照权利要求18-20中任何一个的图像处理设备,其特征在于:所述像素值调整装置通过乘以一个等于或大于1的系数来执行所述像素值的增大。
22.按照权利要求18-20中任何一个的图像处理设备,其特征在于:所述像素值调整装置通过加上预定的数值来执行所述像素值的增大。
23.按照权利要求13-15和18-20中任何一项的图像处理设备,其特征在于还包括:
作为用于确定所述噪声方差的装置的噪声方差计算装置,用于对遍及整个图像的多个定义的局部区域中的每一个、确定像素值的剩余平方和;确定所述剩余平方和的直方图;然后根据给出所述直方图峰值的剩余平方和确定所述噪声方差。
24.按照权利要求13-15和18-20中任何一项的图像处理设备,其特征在于:所述图像是血流图像。
25.一种记录介质,它用于以计算机可读的方式记录使计算机实现以下功能的程序:
把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;以及
当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值。
26.一种记录介质,它以计算机可读的方式记录有使计算机实现以下功能的程序:
把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中像素值的方差;
进行像素值调整,包括当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及
对受到所述像素值调整的所述多切片图像进行最大强度投影。
27.一种记录介质,用于以计算机可读的方式记录使计算机实现以下功能的程序:
把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;
在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中的像素值的方差;
进行像素值调整,包括当所述相加后的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及
对受到所述像素值调整的所述多切片图像进行最大强度投影。
28.一种记录介质,它用于以计算机可读的方式记录使计算机实现以下功能的程序:
把构成图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中像素值的方差;以及
当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值。
29.一种记录介质,它用于以计算机可读的方式记录使计算机实现以下功能的程序:
把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属的局部区域中的像素值的方差;
进行像素值调整,包括当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值;以及
对受到所述像素值调整的所述多切片图像进行最大强度投影。
30.一种记录介质,它以计算机可读的方式记录使计算机实现以下功能的程序:
把构成多切片图像的每个像素定义为感兴趣的像素,确定所述感兴趣的像素所属局部区域中的像素值的方差;
在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中像素值的方差;
进行像素值调整,当所述相加后的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值;以及
对受到所述像素值调整的所述多切片图像进行最大强度投影。
31.一种用来根据从对象收集的信号产生图像的成像设备,它包括:
方差计算装置,用来把构成图像的每个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中像素值的方差;以及
像素值调整装置,用于当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值。
32.一种用来根据从对象收集的信号产生图像的成像设备,它包括:
方差计算装置,用来把构成多切片图像的每一个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中像素值的方差;
像素值调整装置,用于当所述确定的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及
最大强度投影装置,用来对受到所述像素值调整的所述多切片图像进行最大强度投影。
33.一种用来根据从对象收集的信号产生图像的成像设备,它包括:
方差计算装置,用来把构成多切片图像的每一个像素定义为感兴趣的像素、确定所述感兴趣的像素所属局部区域中像素值的方差;
加法装置,用来在所述确定的方差上加上相邻切片的图像中相应的感兴趣像素所属的局部区域中像素值的方差;
像素值调整装置,用于当所述相加后的方差明显大于噪声方差时、维持所述感兴趣的像素的像素值,否则抑制所述感兴趣的像素的像素值;以及
最大强度投影装置,用来对受到所述像素值调整的所述多切片图像进行最大强度投影。
34.一种用来根据从对象收集的信号产生图像的成像设备,它包括:
方差计算装置,用来把构成图像的每个像素定义为感兴趣的像素、确定所述感兴趣的像素所属的局部区域中像素值的方差;以及
像素值调整装置,用于当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值。
35.一种用来根据从对象收集的信号产生图像的成像设备,它包括:
方差计算装置,用来把构成多切片图像的每一个像素定义为感兴趣的像素、确定所述感兴趣的像素所属局部区域中像素值的方差;
像素值调整装置,用于当所述确定的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值;以及
最大强度投影装置,用来对受到所述像素值调整的所述多切片图像进行最大强度投影。
36.一种用来根据从对象收集的信号产生图像的成像设备,它包括:
方差计算装置,用来把构成多切片图像的每一个像素定义为感兴趣的像素、确定所述感兴趣的像素所属局部区域中像素值的方差;
加法装置,用来在所述确定的方差上加上相邻切片的图像中相应的感兴趣的像素所属的局部区域中像素值的方差;
像素值调整装置,用于当所述相加后的方差明显大于噪声方差时、增大所述感兴趣的像素的像素值,否则维持所述感兴趣的像素的像素值;以及
最大强度投影装置,用来对受到所述像素值调整的所述多切片图像进行最大强度投影。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2000369444A JP4149126B2 (ja) | 2000-12-05 | 2000-12-05 | 画像処理方法、画像処理装置および画像撮影装置 |
JP369444/00 | 2000-12-05 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1357302A true CN1357302A (zh) | 2002-07-10 |
CN1287731C CN1287731C (zh) | 2006-12-06 |
Family
ID=18839478
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB011435194A Expired - Fee Related CN1287731C (zh) | 2000-12-05 | 2001-12-05 | 图像处理方法和设备、记录介质和成像设备 |
Country Status (5)
Country | Link |
---|---|
US (1) | US6985613B2 (zh) |
EP (1) | EP1225544A3 (zh) |
JP (1) | JP4149126B2 (zh) |
KR (1) | KR100458776B1 (zh) |
CN (1) | CN1287731C (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101472058B (zh) * | 2007-12-29 | 2011-04-20 | 比亚迪股份有限公司 | 图像噪声去除装置和方法 |
CN101770636B (zh) * | 2008-12-30 | 2012-09-26 | 比亚迪股份有限公司 | 图像噪声去除方法及装置 |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1223553A3 (en) * | 2000-10-17 | 2003-09-24 | Fuji Photo Film Co., Ltd. | Apparatus for suppressing noise by adapting filter characteristics to input image signal based on characteristics of input image signal |
JP4045422B2 (ja) * | 2002-07-23 | 2008-02-13 | 株式会社日立メディコ | 画像処理装置 |
US7417636B2 (en) | 2003-05-08 | 2008-08-26 | Siemens Medical Solutions Usa, Inc. | Method and apparatus for automatic setting of rendering parameter for virtual endoscopy |
US7616792B2 (en) * | 2003-11-18 | 2009-11-10 | Siemens Medical Solutions Usa, Inc. | Scale selection for anisotropic scale-space: application to volumetric tumor characterization |
FR2869723A1 (fr) * | 2004-04-29 | 2005-11-04 | Thomson Licensing Sa | Element de transition sans contact entre un guide d'ondes et une ligne mocroruban |
DE102004025417B4 (de) * | 2004-05-24 | 2008-01-24 | Siemens Ag | Verfahren zur Verbesserung der Bildhomogenität von Bilddaten aus phasenzyklierten Steady-State-Sequenzen |
US7840048B2 (en) * | 2004-05-26 | 2010-11-23 | Guardian Technologies International, Inc. | System and method for determining whether there is an anomaly in data |
US7496218B2 (en) * | 2004-05-26 | 2009-02-24 | Ramsay Thomas E | System and method for identifying objects of interest in image data |
US20080159605A1 (en) | 2005-03-15 | 2008-07-03 | Ramsay Thomas E | Method for characterizing an image source utilizing predetermined color spaces |
US20060269135A1 (en) * | 2005-03-15 | 2006-11-30 | Ramsay Thomas E | System and method for identifying objects of interest in image data |
US8538099B2 (en) * | 2005-03-23 | 2013-09-17 | General Electric Company | Method and system for controlling image reconstruction |
DE102005033851A1 (de) * | 2005-07-20 | 2007-02-01 | Sci-Worx Gmbh | Verfahren zur Reduktion von Rauschartefakten in digitalen Bildern und Einrichtung hierfür |
US20080265889A1 (en) * | 2005-10-07 | 2008-10-30 | Koninklijke Philips Electronics N. V. | Multiple-Channel Transmit Magnetic Resonance |
JP4859516B2 (ja) * | 2006-04-27 | 2012-01-25 | 株式会社ソニー・コンピュータエンタテインメント | 画像処理装置および画像復元方法 |
TWI319676B (en) * | 2006-10-18 | 2010-01-11 | Quanta Comp Inc | Image processing apparatus and method |
EP2309285B1 (en) | 2008-07-24 | 2019-08-21 | Toshiba Medical Systems Corporation | Magnetic resonance imaging apparatus for contrast enhancement of flow images |
JP5414335B2 (ja) * | 2008-07-24 | 2014-02-12 | 株式会社東芝 | 磁気共鳴イメージング装置 |
DE102008058740B4 (de) * | 2008-11-24 | 2012-05-31 | Siemens Aktiengesellschaft | Korrektur von Artefakten in der "Time of Flight"- MR-Angiographie |
JP5036785B2 (ja) * | 2009-10-23 | 2012-09-26 | キヤノン株式会社 | 光断層画像生成方法及び光断層画像生成装置 |
JP5653098B2 (ja) * | 2010-07-09 | 2015-01-14 | Hoya株式会社 | ノイズ除去装置 |
IT1403150B1 (it) | 2010-11-24 | 2013-10-04 | St Microelectronics Srl | Procedimento e dispositivo per depurare dal rumore un segnale video digitale, relativo prodotto informatico. |
JP5703068B2 (ja) * | 2011-02-28 | 2015-04-15 | 綜合警備保障株式会社 | 顔検出方法および顔検出装置、ならびに、プログラム |
JP5438066B2 (ja) * | 2011-05-25 | 2014-03-12 | 日立アロカメディカル株式会社 | 超音波画像処理装置およびプログラム |
JP6073661B2 (ja) * | 2012-11-21 | 2017-02-01 | 株式会社日立製作所 | 磁気共鳴イメージング装置および磁気共鳴イメージング方法 |
US11300646B2 (en) | 2014-11-10 | 2022-04-12 | Canon Medical Systems Corporation | Magnetic resonance imaging apparatus, image processing apparatus, and image processing method |
JP6647836B2 (ja) * | 2014-11-10 | 2020-02-14 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置、画像処理装置及び画像処理方法 |
JP7455508B2 (ja) | 2018-12-26 | 2024-03-26 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置および医用複素数画像処理装置 |
CN113506260B (zh) * | 2021-07-05 | 2023-08-29 | 贝壳找房(北京)科技有限公司 | 一种人脸图像质量评估方法、装置、电子设备及存储介质 |
Family Cites Families (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4761819A (en) * | 1987-02-27 | 1988-08-02 | Picker International, Inc. | Adaptive noise reduction filter for reconstructed images |
EP0355331A3 (en) * | 1988-07-01 | 1991-08-14 | G.D. Searle & Co. | Poly(aminoalkyl)aminocarbonyl aminoacyl aminodiol amino acid derivatives as anti-hypertensive agents |
JPH0751134B2 (ja) * | 1988-12-16 | 1995-06-05 | アロカ株式会社 | 超音波診断画像処理装置 |
US5038388A (en) * | 1989-05-15 | 1991-08-06 | Polaroid Corporation | Method for adaptively sharpening electronic images |
JPH0736818B2 (ja) * | 1990-05-07 | 1995-04-26 | アロカ株式会社 | 超音波診断画像処理装置 |
JP2952687B2 (ja) * | 1990-07-11 | 1999-09-27 | ジーイー横河メディカルシステム株式会社 | 画像処理装置 |
JPH0736820B2 (ja) * | 1990-09-14 | 1995-04-26 | アロカ株式会社 | 画像輪郭強調処理装置 |
US5148479A (en) | 1991-03-20 | 1992-09-15 | International Business Machines Corp. | Authentication protocols in communication networks |
US5233299A (en) * | 1991-03-25 | 1993-08-03 | General Electric Company | Projection methods for producing two-dimensional images from three-dimensional data |
JPH04329930A (ja) * | 1991-04-30 | 1992-11-18 | Shimadzu Corp | 医療用診断画像作成装置 |
GB9114512D0 (en) * | 1991-07-04 | 1991-08-21 | Cyzer John A | Improved label printing apparatus |
DE69331719T2 (de) * | 1992-06-19 | 2002-10-24 | Agfa Gevaert Nv | Verfahren und Vorrichtung zur Geräuschunterdrückung |
JPH06217958A (ja) * | 1993-01-22 | 1994-08-09 | Toshiba Medical Eng Co Ltd | 磁気共鳴イメージング装置 |
EP0610604A1 (en) * | 1993-02-11 | 1994-08-17 | Agfa-Gevaert N.V. | Multiple processing of radiographic images based on a pyramidal image decomposition |
EP0610607A1 (en) * | 1993-02-11 | 1994-08-17 | Agfa-Gevaert N.V. | Method of evaluating image processing performed on a radiographic image |
JP3350176B2 (ja) * | 1993-09-24 | 2002-11-25 | 株式会社東芝 | 磁気共鳴イメージング装置 |
US5933540A (en) * | 1995-05-11 | 1999-08-03 | General Electric Company | Filter system and method for efficiently suppressing noise and improving edge definition in a digitized image |
US5832134A (en) * | 1996-11-27 | 1998-11-03 | General Electric Company | Data visualization enhancement through removal of dominating structures |
US6043655A (en) * | 1997-01-09 | 2000-03-28 | Kabushiki Kaisha Toshiba | MR imaging utilizing the time of flight effect |
JP3896188B2 (ja) | 1997-06-13 | 2007-03-22 | 株式会社日立製作所 | 放射線治療計画のための画像処理装置 |
FR2772225A1 (fr) * | 1997-12-09 | 1999-06-04 | Philips Electronics Nv | Procede de traitement d'une image ou le bruit depend du signal |
US6718068B1 (en) * | 2000-03-10 | 2004-04-06 | Eastman Kodak Company | Noise reduction method utilizing statistical weighting, apparatus, and program for digital image processing |
-
2000
- 2000-12-05 JP JP2000369444A patent/JP4149126B2/ja not_active Expired - Fee Related
-
2001
- 2001-10-08 US US09/973,148 patent/US6985613B2/en not_active Expired - Fee Related
- 2001-12-04 KR KR10-2001-0076051A patent/KR100458776B1/ko not_active IP Right Cessation
- 2001-12-05 EP EP01310192A patent/EP1225544A3/en not_active Withdrawn
- 2001-12-05 CN CNB011435194A patent/CN1287731C/zh not_active Expired - Fee Related
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101472058B (zh) * | 2007-12-29 | 2011-04-20 | 比亚迪股份有限公司 | 图像噪声去除装置和方法 |
CN101770636B (zh) * | 2008-12-30 | 2012-09-26 | 比亚迪股份有限公司 | 图像噪声去除方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
EP1225544A3 (en) | 2004-12-15 |
JP2002183709A (ja) | 2002-06-28 |
CN1287731C (zh) | 2006-12-06 |
EP1225544A2 (en) | 2002-07-24 |
US6985613B2 (en) | 2006-01-10 |
US20020094114A1 (en) | 2002-07-18 |
KR20020044074A (ko) | 2002-06-14 |
JP4149126B2 (ja) | 2008-09-10 |
KR100458776B1 (ko) | 2004-12-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1287731C (zh) | 图像处理方法和设备、记录介质和成像设备 | |
EP2501290B1 (en) | Scan plan field of view adjustor, determiner, and/or quality assessor | |
US20200143521A1 (en) | System and method for image correction | |
US9087259B2 (en) | Organ-specific enhancement filter for robust segmentation of medical images | |
RU2527211C2 (ru) | Реконструкция времяпролетной позитронно-эмиссионной томографии с помощью содержания изображения, формируемого поэтапно на основе времяпролетной информации | |
US7512284B2 (en) | Volumetric image enhancement system and method | |
JP4309062B2 (ja) | 画像化方法 | |
US7248749B2 (en) | Method and apparatus for signal-to-noise ratio dependent image processing | |
US7822254B2 (en) | Automatic positioning of matching multi-planar image reformatting (MPR) views of multiple 3D medical images | |
CN111867474A (zh) | 使用深度学习根据低剂量pet成像进行全剂量pet图像估计 | |
US20100021034A1 (en) | System and method for the generation of attenuation correction maps from mr images | |
CN1698540A (zh) | 图像处理方法、图像处理系统以及x-射线ct系统 | |
CN1237487C (zh) | 图像处理方法和装置、成像装置 | |
CN1591480A (zh) | 进行自动聚焦断层综合的系统和方法 | |
JPH10133309A (ja) | 記録部材における放射線画像の欠陥を補正する方法 | |
US9760991B2 (en) | System and method for image intensity bias estimation and tissue segmentation | |
DE102012201412A1 (de) | Verfahren zum Berechnen eines Wertes eines Absorptionsparameters der Positronen-Emissions-Tomographie | |
EP1716537B1 (en) | Apparatus and method for the processing of sectional images | |
CN110073413A (zh) | 用于迭代图像重建的信息自适应正则化 | |
Mouton et al. | A novel intensity limiting approach to metal artefact reduction in 3D CT baggage imagery | |
CN1305008C (zh) | 弥散张量核磁成像脑缺血病灶区的自动分割方法 | |
US20060066911A1 (en) | Edge detection and correcting system and method | |
Leynes et al. | Synthetic CT generation using MRI with deep learning: How does the selection of input images affect the resulting synthetic CT? | |
CN1580808A (zh) | 失真校正的磁共振测量和磁共振设备 | |
US5663566A (en) | Negativity bias reduction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20061206 Termination date: 20141205 |
|
EXPY | Termination of patent right or utility model |