JP2002024820A - 画像処理方法および装置、記録媒体並びに画像撮影装置 - Google Patents
画像処理方法および装置、記録媒体並びに画像撮影装置Info
- Publication number
- JP2002024820A JP2002024820A JP2000183049A JP2000183049A JP2002024820A JP 2002024820 A JP2002024820 A JP 2002024820A JP 2000183049 A JP2000183049 A JP 2000183049A JP 2000183049 A JP2000183049 A JP 2000183049A JP 2002024820 A JP2002024820 A JP 2002024820A
- Authority
- JP
- Japan
- Prior art keywords
- image
- function
- squares
- histogram
- residual sum
- 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 abstract description 22
- 230000006870 function Effects 0.000 claims description 95
- 238000005315 distribution function Methods 0.000 claims description 50
- 238000004364 calculation method Methods 0.000 claims description 24
- 238000000034 method Methods 0.000 claims description 17
- 239000006185 dispersion Substances 0.000 abstract 2
- 230000005284 excitation Effects 0.000 description 15
- 238000001914 filtration Methods 0.000 description 13
- 238000001208 nuclear magnetic resonance pulse sequence Methods 0.000 description 13
- 238000013480 data collection Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 9
- 230000003068 static effect Effects 0.000 description 9
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 238000003384 imaging method Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 238000002591 computed tomography Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 101100325793 Arabidopsis thaliana BCA2 gene Proteins 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
-
- 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
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Facsimile Image Signal Circuits (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Picture Signal Circuits (AREA)
Abstract
(57)【要約】
【課題】 画像ノイズの分散を求める画像処理方法およ
び装置、そのような画像処理機能をコンピュータに実現
させるプログラムを記録した媒体、並びに、そのような
画像処理装置を備えた画像撮影装置を実現する。 【解決手段】 画像の全体にわたって設定した複数の局
所領域ごとに画素値の残差平方和を求め（５０２〜５０
８）、そのヒストグラムを求め（５１０）、ヒストグラ
ムのピークを与える残差平方和の値に基づいてノイズの
分散を求める（５１２，５１４）。
び装置、そのような画像処理機能をコンピュータに実現
させるプログラムを記録した媒体、並びに、そのような
画像処理装置を備えた画像撮影装置を実現する。 【解決手段】 画像の全体にわたって設定した複数の局
所領域ごとに画素値の残差平方和を求め（５０２〜５０
８）、そのヒストグラムを求め（５１０）、ヒストグラ
ムのピークを与える残差平方和の値に基づいてノイズの
分散を求める（５１２，５１４）。
Description
【０００１】
【発明の属する技術分野】本発明は、画像処理方法およ
び装置、記録媒体並びに画像撮影装置に関し、特に、画
像のノイズ（ｎｏｉｓｅ）の分散を求める画像処理方法
および装置、そのような画像処理機能をコンピュータに
実現させるプログラムを記録した媒体、並びに、そのよ
うな画像処理装置を備えた画像撮影装置に関する。
び装置、記録媒体並びに画像撮影装置に関し、特に、画
像のノイズ（ｎｏｉｓｅ）の分散を求める画像処理方法
および装置、そのような画像処理機能をコンピュータに
実現させるプログラムを記録した媒体、並びに、そのよ
うな画像処理装置を備えた画像撮影装置に関する。
【０００２】
【従来の技術】磁気共鳴撮影（ＭＲＩ：Ｍａｇｎｅｔｉ
ｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置では、
マグネットシステム（ｍａｇｎｅｔ ｓｙｓｔｅｍ）の
内部空間、すなわち、静磁場を形成した空間に撮影対象
を搬入し、勾配磁場および高周波磁場を印加して対象内
に磁気共鳴信号を発生させ、その受信信号に基づいて断
層像を生成（再構成）する。
ｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置では、
マグネットシステム（ｍａｇｎｅｔ ｓｙｓｔｅｍ）の
内部空間、すなわち、静磁場を形成した空間に撮影対象
を搬入し、勾配磁場および高周波磁場を印加して対象内
に磁気共鳴信号を発生させ、その受信信号に基づいて断
層像を生成（再構成）する。
【０００３】断層像の細部構造をよりよく観察できるよ
うにするために、画像のノイズ（ｎｏｉｓｅ）を除去す
るフィルタリング（ｆｉｌｔｅｒｉｎｇ）が行われる。
フィルタリングはローパスフィルタリング（ｌｏｗ−ｐ
ａｓｓ ｆｉｌｔｅｒｉｎｇ）を基本とするが、それだ
けでは画像の鮮鋭度（シャープネス：ｓｈａｒｐｎｅｓ
ｓ）が低下するので、シャープネスを確保するための処
理を付加したフィルタリングが採用される。
うにするために、画像のノイズ（ｎｏｉｓｅ）を除去す
るフィルタリング（ｆｉｌｔｅｒｉｎｇ）が行われる。
フィルタリングはローパスフィルタリング（ｌｏｗ−ｐ
ａｓｓ ｆｉｌｔｅｒｉｎｇ）を基本とするが、それだ
けでは画像の鮮鋭度（シャープネス：ｓｈａｒｐｎｅｓ
ｓ）が低下するので、シャープネスを確保するための処
理を付加したフィルタリングが採用される。
【０００４】
【発明が解決しようとする課題】シャープネス確保処理
を付加したフィルタリングでは、ノイズが偶然に形成し
たに過ぎないテクスチャ（ｔｅｘｔｕｒｅ）を強調し、
解剖学的に意味のない構造（偽構造）を作り出すという
副作用があるので、ノイズだけの部分をシャープネス確
保処理の対象から外すために、画像の局所ごとに、ノイ
ズだけが存在するのか真の構造も存在するのかを判別す
ることが必要になる。
を付加したフィルタリングでは、ノイズが偶然に形成し
たに過ぎないテクスチャ（ｔｅｘｔｕｒｅ）を強調し、
解剖学的に意味のない構造（偽構造）を作り出すという
副作用があるので、ノイズだけの部分をシャープネス確
保処理の対象から外すために、画像の局所ごとに、ノイ
ズだけが存在するのか真の構造も存在するのかを判別す
ることが必要になる。
【０００５】ノイズと真の構造とを区別するには、画像
に含まれるノイズの分散を予め知り、局所の画素値の分
散がノイズの分散に対して有意差が無ければノイズであ
るとし、そうでなければ構造であると判定することが考
えられるが、ノイズと画像信号が混在している画像から
ノイズだけを正確に分離することはできず、したがって
その分散を知ることができなかった。
に含まれるノイズの分散を予め知り、局所の画素値の分
散がノイズの分散に対して有意差が無ければノイズであ
るとし、そうでなければ構造であると判定することが考
えられるが、ノイズと画像信号が混在している画像から
ノイズだけを正確に分離することはできず、したがって
その分散を知ることができなかった。
【０００６】そこで、本発明の課題は、画像ノイズの分
散を求める画像処理方法および装置、そのような画像処
理機能をコンピュータに実現させるプログラムを記録し
た媒体、並びに、そのような画像処理装置を備えた画像
撮影装置を実現することである。
散を求める画像処理方法および装置、そのような画像処
理機能をコンピュータに実現させるプログラムを記録し
た媒体、並びに、そのような画像処理装置を備えた画像
撮影装置を実現することである。
【０００７】
【課題を解決するための手段】（１）上記の課題を解決
するための１つの観点での発明は、画像の全体にわたっ
て設定した複数の局所領域ごとに画素値の残差平方和を
求め、前記残差平方和のヒストグラムを求め、前記ヒス
トグラムのピークを与える残差平方和に基づいてノイズ
の分散を求める、ことを特徴とする画像処理方法であ
る。
するための１つの観点での発明は、画像の全体にわたっ
て設定した複数の局所領域ごとに画素値の残差平方和を
求め、前記残差平方和のヒストグラムを求め、前記ヒス
トグラムのピークを与える残差平方和に基づいてノイズ
の分散を求める、ことを特徴とする画像処理方法であ
る。
【０００８】この観点での発明では、画像の全体にわた
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムのピークを与える残差平方和
に基づいてノイズの分散を求めるので、画像ノイズの分
散を正しく求めることができる。
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムのピークを与える残差平方和
に基づいてノイズの分散を求めるので、画像ノイズの分
散を正しく求めることができる。
【０００９】（２）上記の課題を解決するための他の観
点での発明は、画像の全体にわたって設定した複数の局
所領域ごとに画素値の残差平方和を求め、前記残差平方
和のヒストグラムを求め、前記ヒストグラムを関数にフ
ィッティングし、前記フィッティングした関数のピーク
値を与える変数に基づいてノイズの分散を求める、こと
を特徴とする画像処理方法である。
点での発明は、画像の全体にわたって設定した複数の局
所領域ごとに画素値の残差平方和を求め、前記残差平方
和のヒストグラムを求め、前記ヒストグラムを関数にフ
ィッティングし、前記フィッティングした関数のピーク
値を与える変数に基づいてノイズの分散を求める、こと
を特徴とする画像処理方法である。
【００１０】この観点での発明では、画像の全体にわた
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムを関数にフィッティングし、
その関数のピーク値を与える変数に基づいてノイズの分
散を求めるので、画像ノイズの分散を一層正しく求める
ことができる。
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムを関数にフィッティングし、
その関数のピーク値を与える変数に基づいてノイズの分
散を求めるので、画像ノイズの分散を一層正しく求める
ことができる。
【００１１】（３）上記の課題を解決するための他の観
点での発明は、前記関数はガウス分布関数である、こと
を特徴とする（２）に記載の画像処理方法である。この
観点での発明では、ガウス分布関数にフィッティングす
るので、画像の実質部における画像ノイズについて、関
数のピークを与える変数を正しく求めることができる。
点での発明は、前記関数はガウス分布関数である、こと
を特徴とする（２）に記載の画像処理方法である。この
観点での発明では、ガウス分布関数にフィッティングす
るので、画像の実質部における画像ノイズについて、関
数のピークを与える変数を正しく求めることができる。
【００１２】（４）上記の課題を解決するための他の観
点での発明は、前記関数はレイリー分布関数である、こ
とを特徴とする（２）に記載の画像処理方法である。こ
の観点での発明では、レイリー分布関数にフィッティン
グするので、画像のバックグラウンド部における画像ノ
イズについて、関数のピークを与える変数を正しく求め
ることができる。
点での発明は、前記関数はレイリー分布関数である、こ
とを特徴とする（２）に記載の画像処理方法である。こ
の観点での発明では、レイリー分布関数にフィッティン
グするので、画像のバックグラウンド部における画像ノ
イズについて、関数のピークを与える変数を正しく求め
ることができる。
【００１３】（５）上記の課題を解決するための他の観
点での発明は、前記関数はガウス分布関数およびレイリ
ー分布関数である、ことを特徴とする（２）に記載の画
像処理方法である。
点での発明は、前記関数はガウス分布関数およびレイリ
ー分布関数である、ことを特徴とする（２）に記載の画
像処理方法である。
【００１４】この観点での発明では、ガウス分布関数お
よびレイリー分布関数にそれぞれフィッティングするの
で、画像の実質部およびバックグラウンド部における画
像ノイズについて、関数のピークを与える変数を正しく
求めることができる。
よびレイリー分布関数にそれぞれフィッティングするの
で、画像の実質部およびバックグラウンド部における画
像ノイズについて、関数のピークを与える変数を正しく
求めることができる。
【００１５】（６）上記の課題を解決するための他の観
点での発明は、前記画像の画素値は複素数の絶対値であ
る、ことを特徴とする（１）ないし（５）のうちのいず
れか１つに記載の画像処理方法である。
点での発明は、前記画像の画素値は複素数の絶対値であ
る、ことを特徴とする（１）ないし（５）のうちのいず
れか１つに記載の画像処理方法である。
【００１６】この観点での発明では、複素数の絶対値を
画素値とする画像についてノイズの分散を求めることが
できる。 （７）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の実数部である、ことを
特徴とする（１）ないし（３）のうちのいずれか１つに
記載の画像処理方法である。
画素値とする画像についてノイズの分散を求めることが
できる。 （７）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の実数部である、ことを
特徴とする（１）ないし（３）のうちのいずれか１つに
記載の画像処理方法である。
【００１７】この観点での発明では、複素数の実数部を
画素値とする画像についてノイズの分散を求めることが
できる。 （８）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の虚数部である、ことを
特徴とする（１）ないし（３）のうちのいずれか１つに
記載の画像処理方法である。
画素値とする画像についてノイズの分散を求めることが
できる。 （８）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の虚数部である、ことを
特徴とする（１）ないし（３）のうちのいずれか１つに
記載の画像処理方法である。
【００１８】この観点での発明では、複素数の虚数部を
画素値とする画像についてノイズの分散を求めることが
できる。 （９）上記の課題を解決するための他の観点での発明
は、前記画像は磁気共鳴を利用して撮影した画像であ
る、ことを特徴とする（１）ないし（８）のうちのいず
れか１つに記載の画像処理方法である。
画素値とする画像についてノイズの分散を求めることが
できる。 （９）上記の課題を解決するための他の観点での発明
は、前記画像は磁気共鳴を利用して撮影した画像であ
る、ことを特徴とする（１）ないし（８）のうちのいず
れか１つに記載の画像処理方法である。
【００１９】この観点での発明では、磁気共鳴を利用し
て撮影した画像についてノイズの分散を求めることがで
きる。 （１０）上記の課題を解決するための他の観点での発明
は、画像の全体にわたって設定した複数の局所領域ごと
に画素値の残差平方和を求める残差平方和計算手段と、
前記残差平方和のヒストグラムを求めるヒストグラム計
算手段と、前記ヒストグラムのピークを与える残差平方
和に基づいてノイズの分散を求めるノイズ分散計算手段
と、を具備することを特徴とする画像処理装置である。
て撮影した画像についてノイズの分散を求めることがで
きる。 （１０）上記の課題を解決するための他の観点での発明
は、画像の全体にわたって設定した複数の局所領域ごと
に画素値の残差平方和を求める残差平方和計算手段と、
前記残差平方和のヒストグラムを求めるヒストグラム計
算手段と、前記ヒストグラムのピークを与える残差平方
和に基づいてノイズの分散を求めるノイズ分散計算手段
と、を具備することを特徴とする画像処理装置である。
【００２０】この観点での発明では、画像の全体にわた
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムのピークを与える残差平方和
に基づいてノイズの分散を求めるので、画像ノイズの分
散を正しく求めることができる。
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムのピークを与える残差平方和
に基づいてノイズの分散を求めるので、画像ノイズの分
散を正しく求めることができる。
【００２１】（１１）上記の課題を解決するための他の
観点での発明は、画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求める残差平方和
計算手段と、前記残差平方和のヒストグラムを求めるヒ
ストグラム計算手段と、前記ヒストグラムを関数にフィ
ッティングするフィッティング手段と、前記フィッティ
ングした関数のピーク値を与える変数に基づいてノイズ
の分散を求めるノイズ分散計算手段と、を具備すること
を特徴とする画像処理装置である。
観点での発明は、画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求める残差平方和
計算手段と、前記残差平方和のヒストグラムを求めるヒ
ストグラム計算手段と、前記ヒストグラムを関数にフィ
ッティングするフィッティング手段と、前記フィッティ
ングした関数のピーク値を与える変数に基づいてノイズ
の分散を求めるノイズ分散計算手段と、を具備すること
を特徴とする画像処理装置である。
【００２２】この観点での発明では、画像の全体にわた
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムを関数にフィッティングし、
その関数のピーク値を与える変数に基づいてノイズの分
散を求めるので、画像ノイズの分散を一層正しく求める
ことができる。
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムを関数にフィッティングし、
その関数のピーク値を与える変数に基づいてノイズの分
散を求めるので、画像ノイズの分散を一層正しく求める
ことができる。
【００２３】（１２）上記の課題を解決するための他の
観点での発明は、前記関数はガウス分布関数である、こ
とを特徴とする（１１）に記載の画像処理装置である。
この観点での発明では、ガウス分布関数にフィッティン
グするので、画像の実質部における画像ノイズについ
て、関数のピークを与える変数を正しく求めることがで
きる。
観点での発明は、前記関数はガウス分布関数である、こ
とを特徴とする（１１）に記載の画像処理装置である。
この観点での発明では、ガウス分布関数にフィッティン
グするので、画像の実質部における画像ノイズについ
て、関数のピークを与える変数を正しく求めることがで
きる。
【００２４】（１３）上記の課題を解決するための他の
観点での発明は、前記関数はレイリー分布関数である、
ことを特徴とする（１１）に記載の画像処理装置であ
る。この観点での発明では、レイリー分布関数にフィッ
ティングするので、画像のバックグラウンド部における
画像ノイズについて、関数のピークを与える変数を正し
く求めることができる。
観点での発明は、前記関数はレイリー分布関数である、
ことを特徴とする（１１）に記載の画像処理装置であ
る。この観点での発明では、レイリー分布関数にフィッ
ティングするので、画像のバックグラウンド部における
画像ノイズについて、関数のピークを与える変数を正し
く求めることができる。
【００２５】（１４）上記の課題を解決するための他の
観点での発明は、前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする（１１）に記載
の画像処理装置である。
観点での発明は、前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする（１１）に記載
の画像処理装置である。
【００２６】この観点での発明では、ガウス分布関数お
よびレイリー分布関数にそれぞれフィッティングするの
で、画像の実質部およびバックグラウンド部における画
像ノイズについて、関数のピークを与える変数を正しく
求めることができる。
よびレイリー分布関数にそれぞれフィッティングするの
で、画像の実質部およびバックグラウンド部における画
像ノイズについて、関数のピークを与える変数を正しく
求めることができる。
【００２７】（１５）上記の課題を解決するための他の
観点での発明は、前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする（１０）ないし（１４）のうち
のいずれか１つに記載の画像処理装置である。
観点での発明は、前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする（１０）ないし（１４）のうち
のいずれか１つに記載の画像処理装置である。
【００２８】この観点での発明では、複素数の絶対値を
画素値とする画像についてノイズの分散を求めることが
できる。 （１６）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の実数部である、ことを
特徴とする（１０）ないし（１２）のうちのいずれか１
つに記載の画像処理装置である。
画素値とする画像についてノイズの分散を求めることが
できる。 （１６）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の実数部である、ことを
特徴とする（１０）ないし（１２）のうちのいずれか１
つに記載の画像処理装置である。
【００２９】この観点での発明では、複素数の実数部を
画素値とする画像についてノイズの分散を求めることが
できる。 （１７）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の虚数部である、ことを
特徴とする（１０）ないし（１２）のうちのいずれか１
つに記載の画像処理装置である。
画素値とする画像についてノイズの分散を求めることが
できる。 （１７）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の虚数部である、ことを
特徴とする（１０）ないし（１２）のうちのいずれか１
つに記載の画像処理装置である。
【００３０】この観点での発明では、複素数の虚数部を
画素値とする画像についてノイズの分散を求めることが
できる。 （１８）上記の課題を解決するための他の観点での発明
は、前記画像は磁気共鳴を利用して撮影した画像であ
る、ことを特徴とする（１０）ないし（１７）のうちの
いずれか１つに記載の画像処理装置である。
画素値とする画像についてノイズの分散を求めることが
できる。 （１８）上記の課題を解決するための他の観点での発明
は、前記画像は磁気共鳴を利用して撮影した画像であ
る、ことを特徴とする（１０）ないし（１７）のうちの
いずれか１つに記載の画像処理装置である。
【００３１】この観点での発明では、磁気共鳴を利用し
て撮影した画像についてノイズの分散を求めることがで
きる。 （１９）上記の課題を解決するための他の観点での発明
は、画像の全体にわたって設定した複数の局所領域ごと
に画素値の残差平方和を求める残差平方和計算機能と、
前記残差平方和のヒストグラムを求めるヒストグラム計
算機能と、前記ヒストグラムのピークを与える残差平方
和に基づいてノイズの分散を求めるノイズ分散計算機能
と、をコンピュータに実現させるプログラムをコンピュ
ータにより読み取り可能なように記録したことを特徴と
する記録媒体である。
て撮影した画像についてノイズの分散を求めることがで
きる。 （１９）上記の課題を解決するための他の観点での発明
は、画像の全体にわたって設定した複数の局所領域ごと
に画素値の残差平方和を求める残差平方和計算機能と、
前記残差平方和のヒストグラムを求めるヒストグラム計
算機能と、前記ヒストグラムのピークを与える残差平方
和に基づいてノイズの分散を求めるノイズ分散計算機能
と、をコンピュータに実現させるプログラムをコンピュ
ータにより読み取り可能なように記録したことを特徴と
する記録媒体である。
【００３２】この観点での発明では、記録媒体に記録し
たプログラムが、画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求め、そのヒスト
グラムのピークを与える残差平方和に基づいてノイズの
分散を求める機能を、コンピュータに実現させるので、
画像ノイズの分散を正しく求めることができる。
たプログラムが、画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求め、そのヒスト
グラムのピークを与える残差平方和に基づいてノイズの
分散を求める機能を、コンピュータに実現させるので、
画像ノイズの分散を正しく求めることができる。
【００３３】（２０）上記の課題を解決するための他の
観点での発明は、画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求める残差平方和
計算機能と、前記残差平方和のヒストグラムを求めるヒ
ストグラム計算機能と、前記ヒストグラムを関数にフィ
ッティングするフィッティング機能と、前記フィッティ
ングした関数のピーク値を与える変数に基づいてノイズ
の分散を求めるノイズ分散計算機能と、をコンピュータ
に実現させるプログラムをコンピュータにより読み取り
可能なように記録したことを特徴とする記録媒体であ
る。
観点での発明は、画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求める残差平方和
計算機能と、前記残差平方和のヒストグラムを求めるヒ
ストグラム計算機能と、前記ヒストグラムを関数にフィ
ッティングするフィッティング機能と、前記フィッティ
ングした関数のピーク値を与える変数に基づいてノイズ
の分散を求めるノイズ分散計算機能と、をコンピュータ
に実現させるプログラムをコンピュータにより読み取り
可能なように記録したことを特徴とする記録媒体であ
る。
【００３４】この観点での発明では、記録媒体に記録し
たプログラムが、画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求め、そのヒスト
グラムを関数にフィッティングし、その関数のピーク値
を与える変数に基づいてノイズの分散を求める機能を、
コンピュータに実現させるので、画像ノイズの分散を一
層正しく求めることができる。
たプログラムが、画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求め、そのヒスト
グラムを関数にフィッティングし、その関数のピーク値
を与える変数に基づいてノイズの分散を求める機能を、
コンピュータに実現させるので、画像ノイズの分散を一
層正しく求めることができる。
【００３５】（２１）上記の課題を解決するための他の
観点での発明は、前記関数はガウス分布関数である、こ
とを特徴とする（２０）に記載の記録媒体である。この
観点での発明では、記録媒体に記録したプログラムが、
ガウス分布関数にフィッティングする機能を、コンピュ
ータに実現させるので、画像の実質部における画像ノイ
ズについて、関数のピークを与える変数を正しく求める
ことができる。
観点での発明は、前記関数はガウス分布関数である、こ
とを特徴とする（２０）に記載の記録媒体である。この
観点での発明では、記録媒体に記録したプログラムが、
ガウス分布関数にフィッティングする機能を、コンピュ
ータに実現させるので、画像の実質部における画像ノイ
ズについて、関数のピークを与える変数を正しく求める
ことができる。
【００３６】（２２）上記の課題を解決するための他の
観点での発明は、前記関数はレイリー分布関数である、
ことを特徴とする（２０）に記載の記録媒体である。こ
の観点での発明では、記録媒体に記録したプログラム
が、レイリー分布関数にフィッティングする機能を、コ
ンピュータに実現させるので、画像のバックグラウンド
部における画像ノイズについて、関数のピークを与える
変数を正しく求めることができる。
観点での発明は、前記関数はレイリー分布関数である、
ことを特徴とする（２０）に記載の記録媒体である。こ
の観点での発明では、記録媒体に記録したプログラム
が、レイリー分布関数にフィッティングする機能を、コ
ンピュータに実現させるので、画像のバックグラウンド
部における画像ノイズについて、関数のピークを与える
変数を正しく求めることができる。
【００３７】（２３）上記の課題を解決するための他の
観点での発明は、前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする（２０）に記載
の記録媒体である。
観点での発明は、前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする（２０）に記載
の記録媒体である。
【００３８】この観点での発明では、記録媒体に記録し
たプログラムが、ガウス分布関数およびレイリー分布関
数にそれぞれフィッティングする機能を、コンピュータ
に実現させるので、画像の実質部およびバックグラウン
ド部における画像ノイズについて、関数のピークを与え
る変数を正しく求めることができる。
たプログラムが、ガウス分布関数およびレイリー分布関
数にそれぞれフィッティングする機能を、コンピュータ
に実現させるので、画像の実質部およびバックグラウン
ド部における画像ノイズについて、関数のピークを与え
る変数を正しく求めることができる。
【００３９】（２４）上記の課題を解決するための他の
観点での発明は、前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする（１９）ないし（２３）のうち
のいずれか１つに記載の記録媒体である。
観点での発明は、前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする（１９）ないし（２３）のうち
のいずれか１つに記載の記録媒体である。
【００４０】この観点での発明では、記録媒体に記録し
たプログラムが、複素数の絶対値を画素値とする画像に
ついてノイズの分散を求める機能を、コンピュータに実
現させることができる。
たプログラムが、複素数の絶対値を画素値とする画像に
ついてノイズの分散を求める機能を、コンピュータに実
現させることができる。
【００４１】（２５）上記の課題を解決するための他の
観点での発明は、前記画像の画素値は複素数の実数部で
ある、ことを特徴とする（１９）ないし（２１）のうち
のいずれか１つに記載の記録媒体である。
観点での発明は、前記画像の画素値は複素数の実数部で
ある、ことを特徴とする（１９）ないし（２１）のうち
のいずれか１つに記載の記録媒体である。
【００４２】この観点での発明では、記録媒体に記録し
たプログラムが、複素数の実数部を画素値とする画像に
ついてノイズの分散を求める機能を、コンピュータに実
現させることができる。
たプログラムが、複素数の実数部を画素値とする画像に
ついてノイズの分散を求める機能を、コンピュータに実
現させることができる。
【００４３】（２６）上記の課題を解決するための他の
観点での発明は、前記画像の画素値は複素数の虚数部で
ある、ことを特徴とする（１９）ないし（２１）のうち
のいずれか１つに記載の記録媒体である。
観点での発明は、前記画像の画素値は複素数の虚数部で
ある、ことを特徴とする（１９）ないし（２１）のうち
のいずれか１つに記載の記録媒体である。
【００４４】この観点での発明では、記録媒体に記録し
たプログラムが、複素数の虚数部を画素値とする画像に
ついてノイズの分散を求める機能を、コンピュータに実
現させることができる。
たプログラムが、複素数の虚数部を画素値とする画像に
ついてノイズの分散を求める機能を、コンピュータに実
現させることができる。
【００４５】（２７）上記の課題を解決するための他の
観点での発明は、前記画像は磁気共鳴を利用して撮影し
た画像である、ことを特徴とする（１９）ないし（２
６）のうちのいずれか１つに記載の記録媒体である。
観点での発明は、前記画像は磁気共鳴を利用して撮影し
た画像である、ことを特徴とする（１９）ないし（２
６）のうちのいずれか１つに記載の記録媒体である。
【００４６】この観点での発明では、記録媒体に記録し
たプログラムが、磁気共鳴を利用して撮影した画像につ
いてノイズの分散を求める機能を、コンピュータに実現
させることができる。
たプログラムが、磁気共鳴を利用して撮影した画像につ
いてノイズの分散を求める機能を、コンピュータに実現
させることができる。
【００４７】（２８）上記の課題を解決するための他の
観点での発明は、対象から収集した信号に基づいて画像
を生成する画像撮影装置であって、画像の全体にわたっ
て設定した複数の局所領域ごとに画素値の残差平方和を
求める残差平方和計算手段と、前記残差平方和のヒスト
グラムを求めるヒストグラム計算手段と、前記ヒストグ
ラムのピークを与える残差平方和に基づいてノイズの分
散を求めるノイズ分散計算手段と、を具備することを特
徴とする画像撮影装置である。
観点での発明は、対象から収集した信号に基づいて画像
を生成する画像撮影装置であって、画像の全体にわたっ
て設定した複数の局所領域ごとに画素値の残差平方和を
求める残差平方和計算手段と、前記残差平方和のヒスト
グラムを求めるヒストグラム計算手段と、前記ヒストグ
ラムのピークを与える残差平方和に基づいてノイズの分
散を求めるノイズ分散計算手段と、を具備することを特
徴とする画像撮影装置である。
【００４８】この観点での発明では、画像の全体にわた
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムのピークを与える残差平方和
に基づいてノイズの分散を求めるので、画像ノイズの分
散を正しく求めることができる。
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムのピークを与える残差平方和
に基づいてノイズの分散を求めるので、画像ノイズの分
散を正しく求めることができる。
【００４９】（２９）上記の課題を解決するための他の
観点での発明は、対象から収集した信号に基づいて画像
を生成する画像撮影装置であって、画像の全体にわたっ
て設定した複数の局所領域ごとに画素値の残差平方和を
求める残差平方和計算手段と、前記残差平方和のヒスト
グラムを求めるヒストグラム計算手段と、前記ヒストグ
ラムを関数にフィッティングするフィッティング手段
と、前記フィッティングした関数のピーク値を与える変
数に基づいてノイズの分散を求めるノイズ分散計算手段
と、を具備することを特徴とする画像撮影装置である。
観点での発明は、対象から収集した信号に基づいて画像
を生成する画像撮影装置であって、画像の全体にわたっ
て設定した複数の局所領域ごとに画素値の残差平方和を
求める残差平方和計算手段と、前記残差平方和のヒスト
グラムを求めるヒストグラム計算手段と、前記ヒストグ
ラムを関数にフィッティングするフィッティング手段
と、前記フィッティングした関数のピーク値を与える変
数に基づいてノイズの分散を求めるノイズ分散計算手段
と、を具備することを特徴とする画像撮影装置である。
【００５０】この観点での発明では、画像の全体にわた
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムを関数にフィッティングし、
その関数のピーク値を与える変数に基づいてノイズの分
散を求めるので、画像ノイズの分散を一層正しく求める
ことができる。
って設定した複数の局所領域ごとに画素値の残差平方和
を求め、そのヒストグラムを関数にフィッティングし、
その関数のピーク値を与える変数に基づいてノイズの分
散を求めるので、画像ノイズの分散を一層正しく求める
ことができる。
【００５１】（３０）上記の課題を解決するための他の
観点での発明は、前記関数はガウス分布関数である、こ
とを特徴とする（２９）に記載の画像撮影装置である。
この観点での発明では、ガウス分布関数にフィッティン
グするので、画像の実質部における画像ノイズについ
て、関数のピークを与える変数を正しく求めることがで
きる。
観点での発明は、前記関数はガウス分布関数である、こ
とを特徴とする（２９）に記載の画像撮影装置である。
この観点での発明では、ガウス分布関数にフィッティン
グするので、画像の実質部における画像ノイズについ
て、関数のピークを与える変数を正しく求めることがで
きる。
【００５２】（３１）上記の課題を解決するための他の
観点での発明は、前記関数はレイリー分布関数である、
ことを特徴とする（２９）に記載の画像撮影装置であ
る。この観点での発明では、レイリー分布関数にフィッ
ティングするので、画像のバックグラウンド部における
画像ノイズについて、関数のピークを与える変数を正し
く求めることができる。
観点での発明は、前記関数はレイリー分布関数である、
ことを特徴とする（２９）に記載の画像撮影装置であ
る。この観点での発明では、レイリー分布関数にフィッ
ティングするので、画像のバックグラウンド部における
画像ノイズについて、関数のピークを与える変数を正し
く求めることができる。
【００５３】（３２）上記の課題を解決するための他の
観点での発明は、前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする（２９）に記載
の画像撮影装置である。
観点での発明は、前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする（２９）に記載
の画像撮影装置である。
【００５４】この観点での発明では、ガウス分布関数お
よびレイリー分布関数にそれぞれフィッティングするの
で、画像の実質部およびバックグラウンド部における画
像ノイズについて、関数のピークを与える変数を正しく
求めることができる。
よびレイリー分布関数にそれぞれフィッティングするの
で、画像の実質部およびバックグラウンド部における画
像ノイズについて、関数のピークを与える変数を正しく
求めることができる。
【００５５】（３３）上記の課題を解決するための他の
観点での発明は、前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする（２８）ないし（３２）のうち
のいずれか１つに記載の画像撮影装置である。
観点での発明は、前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする（２８）ないし（３２）のうち
のいずれか１つに記載の画像撮影装置である。
【００５６】この観点での発明では、複素数の絶対値を
画素値とする画像についてノイズの分散を求めることが
できる。 （３４）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の実数部である、ことを
特徴とする（２８）ないし（３０）のうちのいずれか１
つに記載の画像撮影装置である。
画素値とする画像についてノイズの分散を求めることが
できる。 （３４）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の実数部である、ことを
特徴とする（２８）ないし（３０）のうちのいずれか１
つに記載の画像撮影装置である。
【００５７】この観点での発明では、複素数の実数部を
画素値とする画像についてノイズの分散を求めることが
できる。 （３５）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の虚数部である、ことを
特徴とする（２８）ないし（３０）のうちのいずれか１
つに記載の画像撮影装置である。
画素値とする画像についてノイズの分散を求めることが
できる。 （３５）上記の課題を解決するための他の観点での発明
は、前記画像の画素値は複素数の虚数部である、ことを
特徴とする（２８）ないし（３０）のうちのいずれか１
つに記載の画像撮影装置である。
【００５８】この観点での発明では、複素数の虚数部を
画素値とする画像についてノイズの分散を求めることが
できる。 （３６）上記の課題を解決するための他の観点での発明
は、前記信号は磁気共鳴信号である、ことを特徴とする
（２８）ないし（３５）のうちのいずれか１つに記載の
画像撮影装置である。
画素値とする画像についてノイズの分散を求めることが
できる。 （３６）上記の課題を解決するための他の観点での発明
は、前記信号は磁気共鳴信号である、ことを特徴とする
（２８）ないし（３５）のうちのいずれか１つに記載の
画像撮影装置である。
【００５９】この観点での発明では、磁気共鳴を利用し
て撮影した画像についてノイズの分散を求めることがで
きる。
て撮影した画像についてノイズの分散を求めることがで
きる。
【００６０】
【発明の実施の形態】以下、図面を参照して本発明の実
施の形態を詳細に説明する。なお、本発明は実施の形態
に限定されるものではない。図１に画像撮影装置のブロ
ック（ｂｌｏｃｋ）図を示す。本装置は本発明の実施の
形態の一例である。本装置の構成によって、本発明の装
置に関する実施の形態の一例が示される。本装置の動作
によって、本発明の方法に関する実施の形態の一例が示
される。
施の形態を詳細に説明する。なお、本発明は実施の形態
に限定されるものではない。図１に画像撮影装置のブロ
ック（ｂｌｏｃｋ）図を示す。本装置は本発明の実施の
形態の一例である。本装置の構成によって、本発明の装
置に関する実施の形態の一例が示される。本装置の動作
によって、本発明の方法に関する実施の形態の一例が示
される。
【００６１】図１に示すように、本装置はマグネットシ
ステム１００を有する。マグネットシステム１００は主
磁場コイル（ｃｏｉｌ）部１０２、勾配コイル部１０６
およびＲＦ（ｒａｄｉｏ ｆｒｅｑｕｅｎｃｙ）コイル
部１０８を有する。これら各コイル部は概ね円筒状の形
状を有し、互いに同軸的に配置されている。マグネット
システム１００の概ね円柱状の内部空間（ボア：ｂｏｒ
ｅ）に、撮影する対象３００がクレードル（ｃｒａｄｌ
ｅ）５００に搭載されて図示しない搬送手段により搬入
および搬出される。
ステム１００を有する。マグネットシステム１００は主
磁場コイル（ｃｏｉｌ）部１０２、勾配コイル部１０６
およびＲＦ（ｒａｄｉｏ ｆｒｅｑｕｅｎｃｙ）コイル
部１０８を有する。これら各コイル部は概ね円筒状の形
状を有し、互いに同軸的に配置されている。マグネット
システム１００の概ね円柱状の内部空間（ボア：ｂｏｒ
ｅ）に、撮影する対象３００がクレードル（ｃｒａｄｌ
ｅ）５００に搭載されて図示しない搬送手段により搬入
および搬出される。
【００６２】主磁場コイル部１０２はマグネットシステ
ム１００の内部空間に静磁場を形成する。静磁場の方向
は概ね対象３００の体軸の方向に平行である。すなわち
いわゆる水平磁場を形成する。主磁場コイル部１０２は
例えば超伝導コイルを用いて構成される。なお、超伝導
コイルに限らず常伝導コイル等を用いて構成しても良い
のはもちろんである。
ム１００の内部空間に静磁場を形成する。静磁場の方向
は概ね対象３００の体軸の方向に平行である。すなわち
いわゆる水平磁場を形成する。主磁場コイル部１０２は
例えば超伝導コイルを用いて構成される。なお、超伝導
コイルに限らず常伝導コイル等を用いて構成しても良い
のはもちろんである。
【００６３】勾配コイル部１０６は静磁場強度に勾配を
持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス（ｓｌｉｃｅ）勾配磁場、リードアウト
（ｒｅａｄ ｏｕｔ）勾配磁場およびフェーズエンコー
ド（ｐｈａｓｅ ｅｎｃｏｄｅ）勾配磁場の３種であ
り、これら３種類の勾配磁場に対応して勾配コイル部１
０６は図示しない３系統の勾配コイルを有する。
持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス（ｓｌｉｃｅ）勾配磁場、リードアウト
（ｒｅａｄ ｏｕｔ）勾配磁場およびフェーズエンコー
ド（ｐｈａｓｅ ｅｎｃｏｄｅ）勾配磁場の３種であ
り、これら３種類の勾配磁場に対応して勾配コイル部１
０６は図示しない３系統の勾配コイルを有する。
【００６４】ＲＦコイル部１０８は静磁場空間に対象３
００の体内のスピンを励起するための高周波磁場を形成
する。以下、高周波磁場を形成することをＲＦ励起信号
の送信という。ＲＦコイル部１０８は、また、励起され
たスピンが生じる電磁波すなわち磁気共鳴信号を受信す
る。
００の体内のスピンを励起するための高周波磁場を形成
する。以下、高周波磁場を形成することをＲＦ励起信号
の送信という。ＲＦコイル部１０８は、また、励起され
たスピンが生じる電磁波すなわち磁気共鳴信号を受信す
る。
【００６５】ＲＦコイル部１０８は図示しない送信用の
コイルおよび受信用のコイルを有する。送信用のコイル
および受信用のコイルは、同じコイルを兼用するかある
いはそれぞれ専用のコイルを用いる。
コイルおよび受信用のコイルを有する。送信用のコイル
および受信用のコイルは、同じコイルを兼用するかある
いはそれぞれ専用のコイルを用いる。
【００６６】勾配コイル部１０６には勾配駆動部１３０
が接続されている。勾配駆動部１３０は勾配コイル部１
０６に駆動信号を与えて勾配磁場を発生させる。勾配駆
動部１３０は、勾配コイル部１０６における３系統の勾
配コイルに対応して、図示しない３系統の駆動回路を有
する。
が接続されている。勾配駆動部１３０は勾配コイル部１
０６に駆動信号を与えて勾配磁場を発生させる。勾配駆
動部１３０は、勾配コイル部１０６における３系統の勾
配コイルに対応して、図示しない３系統の駆動回路を有
する。
【００６７】ＲＦコイル部１０８にはＲＦ駆動部１４０
が接続されている。ＲＦ駆動部１４０はＲＦコイル部１
０８に駆動信号を与えてＲＦ励起信号を送信し、対象３
００の体内のスピンを励起する。
が接続されている。ＲＦ駆動部１４０はＲＦコイル部１
０８に駆動信号を与えてＲＦ励起信号を送信し、対象３
００の体内のスピンを励起する。
【００６８】ＲＦコイル部１０８にはデータ収集部１５
０が接続されている。データ収集部１５０はＲＦコイル
部１０８が受信した受信信号を取り込み、それをビュー
データ（ｖｉｅｗ ｄａｔａ）として収集する。
０が接続されている。データ収集部１５０はＲＦコイル
部１０８が受信した受信信号を取り込み、それをビュー
データ（ｖｉｅｗ ｄａｔａ）として収集する。
【００６９】勾配駆動部１３０、ＲＦ駆動部１４０およ
びデータ収集部１５０には制御部１６０が接続されてい
る。制御部１６０は、勾配駆動部１３０ないしデータ収
集部１５０をそれぞれ制御して撮影を遂行する。
びデータ収集部１５０には制御部１６０が接続されてい
る。制御部１６０は、勾配駆動部１３０ないしデータ収
集部１５０をそれぞれ制御して撮影を遂行する。
【００７０】データ収集部１５０の出力側はデータ処理
部１７０に接続されている。データ処理部１７０は、例
えばコンピュータ（ｃｏｍｐｕｔｅｒ）等を用いて構成
される。データ処理部１７０は図示しないメモリ（ｍｅ
ｍｏｒｙ）を有する。メモリはデータ処理部１７０用の
プログラムおよび各種のデータを記憶している。本装置
の機能は、データ処理部１７０がメモリに記憶されたプ
ログラムを実行することによりを実現される。
部１７０に接続されている。データ処理部１７０は、例
えばコンピュータ（ｃｏｍｐｕｔｅｒ）等を用いて構成
される。データ処理部１７０は図示しないメモリ（ｍｅ
ｍｏｒｙ）を有する。メモリはデータ処理部１７０用の
プログラムおよび各種のデータを記憶している。本装置
の機能は、データ処理部１７０がメモリに記憶されたプ
ログラムを実行することによりを実現される。
【００７１】データ処理部１７０は、データ収集部１５
０から取り込んだデータをメモリに記憶する。メモリ内
にはデータ空間が形成される。データ空間は２次元フ−
リエ（Ｆｏｕｒｉｅｒ）空間を構成する。２次元フ−リ
エ空間をｋスペース（ｋ−ｓｐａｃｅ）ともいう。デー
タ処理部１７０は、これら２次元フ−リエ空間のデータ
を２次元逆フ−リエ変換して対象３００の画像を生成
（再構成）する。
０から取り込んだデータをメモリに記憶する。メモリ内
にはデータ空間が形成される。データ空間は２次元フ−
リエ（Ｆｏｕｒｉｅｒ）空間を構成する。２次元フ−リ
エ空間をｋスペース（ｋ−ｓｐａｃｅ）ともいう。デー
タ処理部１７０は、これら２次元フ−リエ空間のデータ
を２次元逆フ−リエ変換して対象３００の画像を生成
（再構成）する。
【００７２】２次元逆フ−リエ変換により再構成した画
像の画素値は複素数となる。複素数の絶対値を用いて絶
対値画像を構成する。複素数の実数部（リアルパート：
ｒｅａｌ ｐａｒｔ）を用いて実数部画像を構成するこ
とができる。複素数の虚数部（イマジナリパート：ｉｍ
ａｇｉｎａｒｙ ｐａｒｔ）を用いて虚数部画像を構成
することができる。実数部および虚数部はいずれも正負
の値をとることができる。これらの画像を正負画像とも
いう。
像の画素値は複素数となる。複素数の絶対値を用いて絶
対値画像を構成する。複素数の実数部（リアルパート：
ｒｅａｌ ｐａｒｔ）を用いて実数部画像を構成するこ
とができる。複素数の虚数部（イマジナリパート：ｉｍ
ａｇｉｎａｒｙ ｐａｒｔ）を用いて虚数部画像を構成
することができる。実数部および虚数部はいずれも正負
の値をとることができる。これらの画像を正負画像とも
いう。
【００７３】データ処理部１７０は再構成した画像につ
いてノイズの分散を求めるための画像処理を行う機能を
有する。データ処理部１７０の画像処理機能については
後にあらためて説明する。
いてノイズの分散を求めるための画像処理を行う機能を
有する。データ処理部１７０の画像処理機能については
後にあらためて説明する。
【００７４】データ処理部１７０は、本発明における画
像処理装置の実施の形態の一例である。本装置の構成に
よって、本発明の装置に関する実施の形態の一例が示さ
れる。本装置の動作によって、本発明の方法に関する実
施の形態の一例が示される。
像処理装置の実施の形態の一例である。本装置の構成に
よって、本発明の装置に関する実施の形態の一例が示さ
れる。本装置の動作によって、本発明の方法に関する実
施の形態の一例が示される。
【００７５】データ処理部１７０は制御部１６０に接続
されている。データ処理部１７０は制御部１６０の上位
にあってそれを統括する。データ処理部１７０には表示
部１８０および操作部１９０が接続されている。表示部
１８０は、グラフィックディスプレー（ｇｒａｐｈｉｃ
ｄｉｓｐｌａｙ）等で構成される。操作部１９０はポ
インティングデバイス（ｐｏｉｎｔｉｎｇ ｄｅｖｉｃ
ｅ）を備えたキーボード（ｋｅｙｂｏａｒｄ）等で構成
される。
されている。データ処理部１７０は制御部１６０の上位
にあってそれを統括する。データ処理部１７０には表示
部１８０および操作部１９０が接続されている。表示部
１８０は、グラフィックディスプレー（ｇｒａｐｈｉｃ
ｄｉｓｐｌａｙ）等で構成される。操作部１９０はポ
インティングデバイス（ｐｏｉｎｔｉｎｇ ｄｅｖｉｃ
ｅ）を備えたキーボード（ｋｅｙｂｏａｒｄ）等で構成
される。
【００７６】表示部１８０は、データ処理部１７０から
出力される再構成画像および各種の情報を表示する。操
作部１９０は、操作者によって操作され、各種の指令や
情報等をデータ処理部１７０に入力する。操作者は表示
部１８０および操作部１９０を通じてインタラクティブ
（ｉｎｔｅｒａｃｔｉｖｅ）に本装置を操作する。
出力される再構成画像および各種の情報を表示する。操
作部１９０は、操作者によって操作され、各種の指令や
情報等をデータ処理部１７０に入力する。操作者は表示
部１８０および操作部１９０を通じてインタラクティブ
（ｉｎｔｅｒａｃｔｉｖｅ）に本装置を操作する。
【００７７】図２に、他の方式の画像撮影装置のブロッ
ク図を示す。本装置は本発明の実施の形態の一例であ
る。本装置の構成によって、本発明の装置に関する実施
の形態の一例が示される。
ク図を示す。本装置は本発明の実施の形態の一例であ
る。本装置の構成によって、本発明の装置に関する実施
の形態の一例が示される。
【００７８】図２に示す装置は、図１に示した装置とは
方式を異にするマグネットシステム１００’を有する。
マグネットシステム１００’以外は図１に示した装置と
同様な構成になっており、同様な部分に同一の符号を付
して説明を省略する。
方式を異にするマグネットシステム１００’を有する。
マグネットシステム１００’以外は図１に示した装置と
同様な構成になっており、同様な部分に同一の符号を付
して説明を省略する。
【００７９】マグネットシステム１００’は主磁場マグ
ネット部１０２’、勾配コイル部１０６’およびＲＦコ
イル部１０８’を有する。これら主磁場マグネット部１
０２’および各コイル部は、いずれも空間を挟んで互い
に対向する１対のものからなる。また、いずれも概ね円
盤状の形状を有し中心軸を共有して配置されている。マ
グネットシステム１００’の内部空間（ボア）に、対象
３００がクレードル５００に搭載されて図示しない搬送
手段により搬入および搬出される。
ネット部１０２’、勾配コイル部１０６’およびＲＦコ
イル部１０８’を有する。これら主磁場マグネット部１
０２’および各コイル部は、いずれも空間を挟んで互い
に対向する１対のものからなる。また、いずれも概ね円
盤状の形状を有し中心軸を共有して配置されている。マ
グネットシステム１００’の内部空間（ボア）に、対象
３００がクレードル５００に搭載されて図示しない搬送
手段により搬入および搬出される。
【００８０】主磁場マグネット部１０２’はマグネット
システム１００’の内部空間に静磁場を形成する。静磁
場の方向は概ね対象３００の体軸方向と直交する。すな
わちいわゆる垂直磁場を形成する。主磁場マグネット部
１０２’は例えば永久磁石等を用いて構成される。な
お、永久磁石に限らず超伝導電磁石あるいは常伝導電磁
石等を用いて構成しても良いのはもちろんである。
システム１００’の内部空間に静磁場を形成する。静磁
場の方向は概ね対象３００の体軸方向と直交する。すな
わちいわゆる垂直磁場を形成する。主磁場マグネット部
１０２’は例えば永久磁石等を用いて構成される。な
お、永久磁石に限らず超伝導電磁石あるいは常伝導電磁
石等を用いて構成しても良いのはもちろんである。
【００８１】勾配コイル部１０６’は静磁場強度に勾配
を持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス勾配磁場、リードアウト勾配磁場およびフ
ェーズエンコード勾配磁場の３種であり、これら３種類
の勾配磁場に対応して勾配コイル部１０６’は図示しな
い３系統の勾配コイルを有する。
を持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス勾配磁場、リードアウト勾配磁場およびフ
ェーズエンコード勾配磁場の３種であり、これら３種類
の勾配磁場に対応して勾配コイル部１０６’は図示しな
い３系統の勾配コイルを有する。
【００８２】ＲＦコイル部１０８’は静磁場空間に対象
３００の体内のスピンを励起するためのＲＦ励起信号を
送信する。ＲＦコイル部１０８’は、また、励起された
スピンが生じる磁気共鳴信号を受信する。ＲＦコイル部
１０８’は図示しない送信用のコイルおよび受信用のコ
イルを有する。送信用のコイルおよび受信用のコイル
は、同じコイルを兼用するかあるいはそれぞれ専用のコ
イルを用いる。
３００の体内のスピンを励起するためのＲＦ励起信号を
送信する。ＲＦコイル部１０８’は、また、励起された
スピンが生じる磁気共鳴信号を受信する。ＲＦコイル部
１０８’は図示しない送信用のコイルおよび受信用のコ
イルを有する。送信用のコイルおよび受信用のコイル
は、同じコイルを兼用するかあるいはそれぞれ専用のコ
イルを用いる。
【００８３】図３に、磁気共鳴撮影に用いるパルスシー
ケンス（ｐｕｌｓｅ ｓｅｑｕｅｎｃｅ）の一例を示
す。このパルスシーケンスは、グラディエントエコー
（ＧＲＥ：Ｇｒａｄｉｅｎｔ Ｅｃｈｏ）法のパルスシ
ーケンスである。
ケンス（ｐｕｌｓｅ ｓｅｑｕｅｎｃｅ）の一例を示
す。このパルスシーケンスは、グラディエントエコー
（ＧＲＥ：Ｇｒａｄｉｅｎｔ Ｅｃｈｏ）法のパルスシ
ーケンスである。
【００８４】すなわち、（１）はＧＲＥ法におけるＲＦ
励起用のα°パルスのシーケンスであり、（２）、
（３）、（４）および（５）は、同じくそれぞれ、スラ
イス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコ
ード勾配ＧｐおよびグラディエントエコーＭＲのシーケ
ンスである。なお、α°パルスは中心信号で代表する。
パルスシーケンスは時間軸ｔに沿って左から右に進行す
る。
励起用のα°パルスのシーケンスであり、（２）、
（３）、（４）および（５）は、同じくそれぞれ、スラ
イス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコ
ード勾配ＧｐおよびグラディエントエコーＭＲのシーケ
ンスである。なお、α°パルスは中心信号で代表する。
パルスシーケンスは時間軸ｔに沿って左から右に進行す
る。
【００８５】同図に示すように、α°パルスによりスピ
ンのα°励起が行われる。フリップアングル（ｆｌｉｐ
ａｎｇｌｅ）α°は９０°以下である。このときスラ
イス勾配Ｇｓが印加され所定のスライスについての選択
励起が行われる。
ンのα°励起が行われる。フリップアングル（ｆｌｉｐ
ａｎｇｌｅ）α°は９０°以下である。このときスラ
イス勾配Ｇｓが印加され所定のスライスについての選択
励起が行われる。
【００８６】α°励起後、フェーズエンコード勾配Ｇｐ
によりスピンのフェーズエンコードが行われる。次に、
リードアウト勾配Ｇｒにより先ずスピンをディフェーズ
（ｄｅｐｈａｓｅ）し、次いでスピンをリフェーズ（ｒ
ｅｐｈａｓｅ）して、グラディエントエコーＭＲを発生
させる。グラディエントエコーＭＲの信号強度は、α°
励起からエコータイム（ｅｃｈｏ ｔｉｍｅ）ＴＥ後の
時点で最大となる。グラディエントエコーＭＲはデータ
収集部１５０によりビューデータとして収集される。
によりスピンのフェーズエンコードが行われる。次に、
リードアウト勾配Ｇｒにより先ずスピンをディフェーズ
（ｄｅｐｈａｓｅ）し、次いでスピンをリフェーズ（ｒ
ｅｐｈａｓｅ）して、グラディエントエコーＭＲを発生
させる。グラディエントエコーＭＲの信号強度は、α°
励起からエコータイム（ｅｃｈｏ ｔｉｍｅ）ＴＥ後の
時点で最大となる。グラディエントエコーＭＲはデータ
収集部１５０によりビューデータとして収集される。
【００８７】このようなパルスシーケンスが周期ＴＲ
（ｒｅｐｅｔｉｔｉｏｎ ｔｉｍｅ）で６４〜５１２回
繰り返される。繰り返しのたびにフェーズエンコード勾
配Ｇｐを変更し、毎回異なるフェーズエンコードを行
う。これによって、ｋスペースを埋める６４〜５１２ビ
ューのビューデータが得られる。
（ｒｅｐｅｔｉｔｉｏｎ ｔｉｍｅ）で６４〜５１２回
繰り返される。繰り返しのたびにフェーズエンコード勾
配Ｇｐを変更し、毎回異なるフェーズエンコードを行
う。これによって、ｋスペースを埋める６４〜５１２ビ
ューのビューデータが得られる。
【００８８】磁気共鳴撮影用パルスシーケンスの他の例
を図４に示す。このパルスシーケンスは、スピンエコー
（ＳＥ：Ｓｐｉｎ Ｅｃｈｏ）法のパルスシーケンスで
ある。
を図４に示す。このパルスシーケンスは、スピンエコー
（ＳＥ：Ｓｐｉｎ Ｅｃｈｏ）法のパルスシーケンスで
ある。
【００８９】すなわち、（１）はＳＥ法におけるＲＦ励
起用の９０°パルスおよび１８０°パルスのシーケンス
であり、（２）、（３）、（４）および（５）は、同じ
くそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇ
ｒ、フェーズエンコード勾配ＧｐおよびスピンエコーＭ
Ｒのシーケンスである。なお、９０°パルスおよび１８
０°パルスはそれぞれ中心信号で代表する。パルスシー
ケンスは時間軸ｔに沿って左から右に進行する。
起用の９０°パルスおよび１８０°パルスのシーケンス
であり、（２）、（３）、（４）および（５）は、同じ
くそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇ
ｒ、フェーズエンコード勾配ＧｐおよびスピンエコーＭ
Ｒのシーケンスである。なお、９０°パルスおよび１８
０°パルスはそれぞれ中心信号で代表する。パルスシー
ケンスは時間軸ｔに沿って左から右に進行する。
【００９０】同図に示すように、９０°パルスによりス
ピンの９０°励起が行われる。このときスライス勾配Ｇ
ｓが印加され所定のスライスについての選択励起が行わ
れる。９０°励起から所定の時間後に、１８０°パルス
による１８０°励起すなわちスピン反転が行われる。こ
のときもスライス勾配Ｇｓが印加され、同じスライスに
ついての選択的反転が行われる。
ピンの９０°励起が行われる。このときスライス勾配Ｇ
ｓが印加され所定のスライスについての選択励起が行わ
れる。９０°励起から所定の時間後に、１８０°パルス
による１８０°励起すなわちスピン反転が行われる。こ
のときもスライス勾配Ｇｓが印加され、同じスライスに
ついての選択的反転が行われる。
【００９１】９０°励起とスピン反転の間の期間に、リ
ードアウト勾配Ｇｒおよびフェーズエンコード勾配Ｇｐ
が印加される。リードアウト勾配Ｇｒによりスピンのデ
ィフェーズが行われる。フェーズエンコード勾配Ｇｐに
よりスピンのフェーズエンコードが行われる。
ードアウト勾配Ｇｒおよびフェーズエンコード勾配Ｇｐ
が印加される。リードアウト勾配Ｇｒによりスピンのデ
ィフェーズが行われる。フェーズエンコード勾配Ｇｐに
よりスピンのフェーズエンコードが行われる。
【００９２】スピン反転後、リードアウト勾配Ｇｒでス
ピンをリフェーズしてスピンエコーＭＲを発生させる。
スピンエコーＭＲの信号強度は、９０°励起からＴＥ後
の時点で最大となる。スピンエコーＭＲはデータ収集部
１５０によりビューデータとして収集される。このよう
なパルスシーケンスが周期ＴＲで６４〜５１２回繰り返
される。繰り返しのたびにフェーズエンコード勾配Ｇｐ
を変更し、毎回異なるフェーズエンコードを行う。これ
によって、ｋスペースを埋める６４〜５１２ビューのビ
ューデータが得られる。
ピンをリフェーズしてスピンエコーＭＲを発生させる。
スピンエコーＭＲの信号強度は、９０°励起からＴＥ後
の時点で最大となる。スピンエコーＭＲはデータ収集部
１５０によりビューデータとして収集される。このよう
なパルスシーケンスが周期ＴＲで６４〜５１２回繰り返
される。繰り返しのたびにフェーズエンコード勾配Ｇｐ
を変更し、毎回異なるフェーズエンコードを行う。これ
によって、ｋスペースを埋める６４〜５１２ビューのビ
ューデータが得られる。
【００９３】なお、撮影に用いるパルスシーケンスはＧ
ＲＥ法またはＳＥ法に限るものではなく、例えば、ＦＳ
Ｅ（Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法、ファーストリ
カバリＦＳＥ（Ｆａｓｔ Ｒｅｃｏｖｅｒｙ Ｆａｓｔ
Ｓｐｉｎ Ｅｃｈｏ）法、エコープラナー・イメージ
ング（ＥＰＩ：Ｅｃｈｏ Ｐｌａｎａｒ Ｉｍａｇｉｎ
ｇ）等、他の適宜の技法のものであって良い。
ＲＥ法またはＳＥ法に限るものではなく、例えば、ＦＳ
Ｅ（Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法、ファーストリ
カバリＦＳＥ（Ｆａｓｔ Ｒｅｃｏｖｅｒｙ Ｆａｓｔ
Ｓｐｉｎ Ｅｃｈｏ）法、エコープラナー・イメージ
ング（ＥＰＩ：Ｅｃｈｏ Ｐｌａｎａｒ Ｉｍａｇｉｎ
ｇ）等、他の適宜の技法のものであって良い。
【００９４】データ処理部１７０は、ｋスペースのビュ
ーデータを２次元逆フ−リエ変換して対象３００の断層
像を再構成する。再構成した画像はメモリに記憶し、ま
た、表示部１８０で表示する。
ーデータを２次元逆フ−リエ変換して対象３００の断層
像を再構成する。再構成した画像はメモリに記憶し、ま
た、表示部１８０で表示する。
【００９５】磁気共鳴撮影した画像は局所的には構造が
一様となることが多いという特質を有する。構造が一様
な局所領域にノイズが存在する場合、画素値の分布は、
領域の画素値の平均値を中心とするガウス（Ｇａｕｓ
ｓ）分布となり、その標準偏差σがノイズの分散を表
す。このような性質に着目し、次のような画像処理によ
って画像のノイズの分散を求める。
一様となることが多いという特質を有する。構造が一様
な局所領域にノイズが存在する場合、画素値の分布は、
領域の画素値の平均値を中心とするガウス（Ｇａｕｓ
ｓ）分布となり、その標準偏差σがノイズの分散を表
す。このような性質に着目し、次のような画像処理によ
って画像のノイズの分散を求める。
【００９６】図５に、データ処理部１７０による画像処
理動作のフローチャート（ｆｌｏｗｃｈａｒｔ）を示
す。同図に示すように、ステップ（ｓｔｅｐ）５０２
で、画像における局所領域を指定する。局所領域とは、
次のステップでの計算に用いる画素値が属する領域であ
る。最初の領域としては、例えば画像の中央における局
所領域が指定される。
理動作のフローチャート（ｆｌｏｗｃｈａｒｔ）を示
す。同図に示すように、ステップ（ｓｔｅｐ）５０２
で、画像における局所領域を指定する。局所領域とは、
次のステップでの計算に用いる画素値が属する領域であ
る。最初の領域としては、例えば画像の中央における局
所領域が指定される。
【００９７】局所領域としてＮｘＮの画素マトリクス
（ｍａｔｒｉｘ）が採用される。Ｎは例えば９である。
なお、マトリクスサイズ（ｍａｔｒｉｘ ｓｉｚｅ）こ
れに限るものではなく、適宜のマトリクスサイズとして
良い。また、ＮｘＮの画素マトリクスに限るものではな
く、画素を中心とする適宜の領域として良い。以下、局
所領域を単に領域ともいう。
（ｍａｔｒｉｘ）が採用される。Ｎは例えば９である。
なお、マトリクスサイズ（ｍａｔｒｉｘ ｓｉｚｅ）こ
れに限るものではなく、適宜のマトリクスサイズとして
良い。また、ＮｘＮの画素マトリクスに限るものではな
く、画素を中心とする適宜の領域として良い。以下、局
所領域を単に領域ともいう。
【００９８】次に、ステップ５０４で、領域に属する画
素値の残差平方和Ｓを求める。すなわち、
素値の残差平方和Ｓを求める。すなわち、
【００９９】
【数１】
【０１００】ここで、Ｐｉは画素値、
【０１０１】
【数２】
【０１０２】は、Ｐｉを中心とするＮｘＮ領域の画素値
の平均値である。また、ｋは例えば８１である。次に、
ステップ５０６で、全ての局所領域について以上の処理
を済ませたか否かを判定し、未済の場合はステップ５０
８で局所領域を変更する。これによって例えば隣のＮｘ
Ｎの領域が新たな局所領域となる。
の平均値である。また、ｋは例えば８１である。次に、
ステップ５０６で、全ての局所領域について以上の処理
を済ませたか否かを判定し、未済の場合はステップ５０
８で局所領域を変更する。これによって例えば隣のＮｘ
Ｎの領域が新たな局所領域となる。
【０１０３】この新たな局所領域について、ステップ５
０４の処理を行い、画素値の残差平方和を求める。以下
同様にして、画像における全ての局所領域について画素
値の残差平方和を求める。
０４の処理を行い、画素値の残差平方和を求める。以下
同様にして、画像における全ての局所領域について画素
値の残差平方和を求める。
【０１０４】ステップ５０２〜５０８の処理を行うデー
タ処理部１７０は、本発明における残差平方和計算手段
の実施の形態の一例である。また、ステップ５０２〜５
０８の処理は、本発明における残差平方和計算機能の実
施の形態の一例である。
タ処理部１７０は、本発明における残差平方和計算手段
の実施の形態の一例である。また、ステップ５０２〜５
０８の処理は、本発明における残差平方和計算機能の実
施の形態の一例である。
【０１０５】上記のようにして求めた残差平方和は、
【０１０６】
【数３】
【０１０７】分布となり、その平均値は、
【０１０８】
【数４】
【０１０９】となる。ｋが大きい場合は、
【０１１０】
【数５】
【０１１１】分布はガウス分布に近づき、そのピークの
位置が、ほぼ
位置が、ほぼ
【０１１２】
【数６】
【０１１３】となる。次に、ステップ５１０で、残差平
方和Ｓのヒストグラム（ｈｉｓｔｏｇｒａｍ）を形成す
る。ステップ５１０の処理を行うデータ処理部１７０
は、本発明におけるヒストグラム計算手段の実施の形態
の一例である。また、ステップ５１０の処理は、本発明
におけるヒストグラム計算機能の実施の形態の一例であ
る。
方和Ｓのヒストグラム（ｈｉｓｔｏｇｒａｍ）を形成す
る。ステップ５１０の処理を行うデータ処理部１７０
は、本発明におけるヒストグラム計算手段の実施の形態
の一例である。また、ステップ５１０の処理は、本発明
におけるヒストグラム計算機能の実施の形態の一例であ
る。
【０１１４】図６に、画像が絶対値画像である場合の、
残差平方和Ｓのヒストグラムを概念的に示す。同図に示
すように、ヒストグラムは３つの分布曲線ａ，ｂ，ｃか
らなる。
残差平方和Ｓのヒストグラムを概念的に示す。同図に示
すように、ヒストグラムは３つの分布曲線ａ，ｂ，ｃか
らなる。
【０１１５】分布曲線ａはガウス分布曲線であり、一様
構造部分におけるノイズによるものである。分布曲線ｂ
はレイリー（Ｒａｙｌｅｉｇｈ）分布曲線であり、ＦＯ
Ｖ（Ｆｉｅｌｄ ｏｆ Ｖｉｅｗ）のうち対象３００が
存在しない部分、すなわち、バックグラウンド（ｂａｃ
ｋｇｒｏｕｎｄ）におけるノイズによるものである。絶
対値画像であることにより、バックグラウンドのノイズ
によるものはガウス分布ではなくレイリー分布となる。
分布曲線ｃは対象の微細構造によるもので、前二者とは
異なり不特定な分布を示す。
構造部分におけるノイズによるものである。分布曲線ｂ
はレイリー（Ｒａｙｌｅｉｇｈ）分布曲線であり、ＦＯ
Ｖ（Ｆｉｅｌｄ ｏｆ Ｖｉｅｗ）のうち対象３００が
存在しない部分、すなわち、バックグラウンド（ｂａｃ
ｋｇｒｏｕｎｄ）におけるノイズによるものである。絶
対値画像であることにより、バックグラウンドのノイズ
によるものはガウス分布ではなくレイリー分布となる。
分布曲線ｃは対象の微細構造によるもので、前二者とは
異なり不特定な分布を示す。
【０１１６】このようなヒストグラムについて、ステッ
プ５１２で、ピーク位置検出を行う。これによって、ガ
ウス分布曲線ａについてはピーク位置ｓ１が検出され、
レイリー分布曲線ｂについてはピーク位置ｓ２が検出さ
れる。
プ５１２で、ピーク位置検出を行う。これによって、ガ
ウス分布曲線ａについてはピーク位置ｓ１が検出され、
レイリー分布曲線ｂについてはピーク位置ｓ２が検出さ
れる。
【０１１７】ヒストグラムは実際は離散的（ディスクリ
ート：ｄｉｓｃｒｅｔｅ）な値をとるのでピーク検出に
先立って、ステップ５１２で、関数フィッティング（ｆ
ｉｔｔｉｎｇ）を行うのが、ピーク位置を精度良く求め
る点で好ましい。フィッティングに用いる関数として
は、例えば、ガウス分布関数およびレイリー分布関数を
それぞれ用いるが、それに限るものではなく、他の適宜
の関数であって良い。
ート：ｄｉｓｃｒｅｔｅ）な値をとるのでピーク検出に
先立って、ステップ５１２で、関数フィッティング（ｆ
ｉｔｔｉｎｇ）を行うのが、ピーク位置を精度良く求め
る点で好ましい。フィッティングに用いる関数として
は、例えば、ガウス分布関数およびレイリー分布関数を
それぞれ用いるが、それに限るものではなく、他の適宜
の関数であって良い。
【０１１８】フィッティングを行うデータ処理部１７０
は、本発明におけるフィッティング手段の実施の形態の
一例である。また、ステップ５１２の処理は、本発明に
おけるフィッティング機能の実施の形態の一例である。
は、本発明におけるフィッティング手段の実施の形態の
一例である。また、ステップ５１２の処理は、本発明に
おけるフィッティング機能の実施の形態の一例である。
【０１１９】次に、ステップ５１４で、ノイズの分散を
算出する。ノイズの分散の算出はピーク位置ｓ１または
ｓ２に基づいて行われる。ｓ１，ｓ２とσの間にはそれ
ぞれ、
算出する。ノイズの分散の算出はピーク位置ｓ１または
ｓ２に基づいて行われる。ｓ１，ｓ２とσの間にはそれ
ぞれ、
【０１２０】
【数７】
【０１２１】
【数８】
【０１２２】の関係があるので、これら関係からσの値
を求める。（２），（３）式のどちらから求めてもσの
値は同じになる。求めたσの値はメモリに記憶して、別
途画像のフィルタリングを行う際に利用する。
を求める。（２），（３）式のどちらから求めてもσの
値は同じになる。求めたσの値はメモリに記憶して、別
途画像のフィルタリングを行う際に利用する。
【０１２３】分布曲線ｃの状態によってはガウス分布曲
線ａのピーク位置ｓ１が精度良く検出できない場合があ
る。その場合はレイリー分布曲線ｂのピーク位置ｓ２に
基づいてσの値を求める。また、バックグラウンド部分
の面積の比率が大きい画像では、レイリー分布曲線ｂの
方がノイズの分散を精度良く求めるのに適する。
線ａのピーク位置ｓ１が精度良く検出できない場合があ
る。その場合はレイリー分布曲線ｂのピーク位置ｓ２に
基づいてσの値を求める。また、バックグラウンド部分
の面積の比率が大きい画像では、レイリー分布曲線ｂの
方がノイズの分散を精度良く求めるのに適する。
【０１２４】ステップ５１２，５１４の処理を行うデー
タ処理部１７０は、本発明におけるノイズ分散計算手段
の実施の形態の一例である。また、ステップ５１２，５
１４の処理は、本発明におけるノイズ分散計算機能の実
施の形態の一例である。
タ処理部１７０は、本発明におけるノイズ分散計算手段
の実施の形態の一例である。また、ステップ５１２，５
１４の処理は、本発明におけるノイズ分散計算機能の実
施の形態の一例である。
【０１２５】以上は、絶対値画像の場合であるが、処理
対象の画像が正負画像、すなわち、実数部画像または虚
数部画像である場合は、バックグラウンド部分における
ノイズは０を中心とする正負の値をとる。
対象の画像が正負画像、すなわち、実数部画像または虚
数部画像である場合は、バックグラウンド部分における
ノイズは０を中心とする正負の値をとる。
【０１２６】このため、ステップ５１０で形成したヒス
トグラムは、例えば図７に示すようになり、レイリー分
布を持たないのとなる。そのような場合は、ステップ５
１４では、ガウス分布曲線ａのピーク位置ｓ１に基づい
てノイズの分散を求める。
トグラムは、例えば図７に示すようになり、レイリー分
布を持たないのとなる。そのような場合は、ステップ５
１４では、ガウス分布曲線ａのピーク位置ｓ１に基づい
てノイズの分散を求める。
【０１２７】以上のようにして求めたノイズの分散は、
対象画像の局所的な構造に応じてフィルタリング方法を
適宜に切り換えるための判定基準として用いることがで
きる。
対象画像の局所的な構造に応じてフィルタリング方法を
適宜に切り換えるための判定基準として用いることがで
きる。
【０１２８】すなわち、注目画素を含む局所領域におけ
る画素値の分散がノイズの分散に対して有意差がない場
合は、その局所領域の画像は格別な構造を持たず、画素
値の分散はノイズによるものである可能性が高い。その
ような場合は、例えば、その領域の画素値をローパスフ
ィルタリングすることにより注目画素の画素値を確定す
る。
る画素値の分散がノイズの分散に対して有意差がない場
合は、その局所領域の画像は格別な構造を持たず、画素
値の分散はノイズによるものである可能性が高い。その
ような場合は、例えば、その領域の画素値をローパスフ
ィルタリングすることにより注目画素の画素値を確定す
る。
【０１２９】これに対して、注目画素を含む局所領域に
おける画素値の分散がノイズの分散に対して有意差があ
る場合は、その局所領域の画像はエッジ等の特定の構造
を持ち、画素値の分散は画像の構造によるものである可
能性が高い。そのような場合は、例えば、その領域の画
素値についてシャープネス確保処理付フィルタリングを
行うことにより注目画素の画素値を確定する。
おける画素値の分散がノイズの分散に対して有意差があ
る場合は、その局所領域の画像はエッジ等の特定の構造
を持ち、画素値の分散は画像の構造によるものである可
能性が高い。そのような場合は、例えば、その領域の画
素値についてシャープネス確保処理付フィルタリングを
行うことにより注目画素の画素値を確定する。
【０１３０】図８に、ノイズの分散を基準値としてフィ
ルタリングの切換を行う場合のフロー図を示す。同図に
示すように、ステップ８０２で、画像中の注目画素を指
定する。最初の注目画素は例えば画像の中央部の画素と
する。
ルタリングの切換を行う場合のフロー図を示す。同図に
示すように、ステップ８０２で、画像中の注目画素を指
定する。最初の注目画素は例えば画像の中央部の画素と
する。
【０１３１】次に、ステップ８０４で注目画素を含む局
所領域の画素値の分散を計算する。注目画素を含む局所
領域は、図９に示すように、注目画素ｉを中心とする例
えば５ｘ５のマトリクスとする。
所領域の画素値の分散を計算する。注目画素を含む局所
領域は、図９に示すように、注目画素ｉを中心とする例
えば５ｘ５のマトリクスとする。
【０１３２】次に、ステップ８０６で、画素値の分散が
ノイズの分散より大きいか否かを判定する。画素値の分
散がノイズの分散より大きいときは、ステップ８０８
で、シャープネス確保処理付フィルタリングを行い、そ
の結果を注目画素の画素値とする。画素値の分散がノイ
ズの分散より大きくないときは、ステップ８１０で、ロ
ーパスフィルタリングを行い、その結果を注目画素の画
素値とする。
ノイズの分散より大きいか否かを判定する。画素値の分
散がノイズの分散より大きいときは、ステップ８０８
で、シャープネス確保処理付フィルタリングを行い、そ
の結果を注目画素の画素値とする。画素値の分散がノイ
ズの分散より大きくないときは、ステップ８１０で、ロ
ーパスフィルタリングを行い、その結果を注目画素の画
素値とする。
【０１３３】次に、ステップ８１２で全ての注目画素に
ついて上記の処理を済ませたか否かを判定し、未済の場
合は、ステップ８１４で、注目画素を例えば隣の画素に
変更してステップ８０４以降の処理を行う。以下、同様
の処理を繰り返し、全ての注目画素の画素値を確定す
る。そして、ステップ７１６で、確定した画素値で画像
を形成する。
ついて上記の処理を済ませたか否かを判定し、未済の場
合は、ステップ８１４で、注目画素を例えば隣の画素に
変更してステップ８０４以降の処理を行う。以下、同様
の処理を繰り返し、全ての注目画素の画素値を確定す
る。そして、ステップ７１６で、確定した画素値で画像
を形成する。
【０１３４】以上のような機能をコンピュータに実現さ
せるプログラムが、記録媒体に、コンピュータで読み取
り可能なように記録される。記録媒体としては、例え
ば、磁気記録媒体、光記録媒体、光磁気記録媒体および
その他の方式の適宜の記録媒体が用いられる。記録媒体
は半導体記憶媒体であっても良い。本書では記憶媒体は
記録媒体と同義である。
せるプログラムが、記録媒体に、コンピュータで読み取
り可能なように記録される。記録媒体としては、例え
ば、磁気記録媒体、光記録媒体、光磁気記録媒体および
その他の方式の適宜の記録媒体が用いられる。記録媒体
は半導体記憶媒体であっても良い。本書では記憶媒体は
記録媒体と同義である。
【０１３５】以上、画像処理を磁気共鳴撮影装置のデー
タ処理部で行う例で説明したが、画像処理は、例えばＥ
ＷＳ（Ｅｎｇｉｎｅｅｒｉｎｇ Ｗｏｒｋ Ｓｔａｔｉ
ｏｎ）やＰＣ（ｐｅｒｓｏｎａｌ ｃｏｍｐｕｔｅｒ）
等、磁気共鳴撮影装置とは別体のデータ処理装置で行う
ようにしても良いのはもちろんである。
タ処理部で行う例で説明したが、画像処理は、例えばＥ
ＷＳ（Ｅｎｇｉｎｅｅｒｉｎｇ Ｗｏｒｋ Ｓｔａｔｉ
ｏｎ）やＰＣ（ｐｅｒｓｏｎａｌ ｃｏｍｐｕｔｅｒ）
等、磁気共鳴撮影装置とは別体のデータ処理装置で行う
ようにしても良いのはもちろんである。
【０１３６】また、画像撮影装置が磁気共鳴撮影装置で
ある例で説明したが、それに限るものではなく、例えば
Ｘ線ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）
装置、Ｘ線撮影装置、ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍ
ｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）、ガンマカメラ
（γ ｃａｍｅｒａ）等、他の方式の画像撮影装置であ
って良い。
ある例で説明したが、それに限るものではなく、例えば
Ｘ線ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）
装置、Ｘ線撮影装置、ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍ
ｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）、ガンマカメラ
（γ ｃａｍｅｒａ）等、他の方式の画像撮影装置であ
って良い。
【０１３７】また、医用画像を処理する例で説明した
が、処理対象は医用画像に限るものではなく、例えば光
学器械で撮影したディジタル画像等、多様な画像のノイ
ズ分散の計算に一般的に適用することができる。
が、処理対象は医用画像に限るものではなく、例えば光
学器械で撮影したディジタル画像等、多様な画像のノイ
ズ分散の計算に一般的に適用することができる。
【０１３８】
【発明の効果】以上詳細に説明したように、本発明によ
れば、画像ノイズの分散を求める画像処理方法および装
置、そのような画像処理機能をコンピュータに実現させ
るプログラムを記録した媒体、並びに、そのような画像
処理装置を備えた画像撮影装置を実現することができ
る。
れば、画像ノイズの分散を求める画像処理方法および装
置、そのような画像処理機能をコンピュータに実現させ
るプログラムを記録した媒体、並びに、そのような画像
処理装置を備えた画像撮影装置を実現することができ
る。
【図１】本発明の実施の形態の一例の装置のブロック図
である。
である。
【図２】本発明の実施の形態の一例の装置のブロック図
である。
である。
【図３】図１または図２に示した装置が実行するパルス
シーケンスの一例を示す図である。
シーケンスの一例を示す図である。
【図４】図１または図２に示した装置が実行するパルス
シーケンスの一例を示す図である。
シーケンスの一例を示す図である。
【図５】図１または図２に示した装置が行う画像処理の
フローチャートである。
フローチャートである。
【図６】ヒストグラムの概念図である。
【図７】ヒストグラムの概念図である。
【図８】図１または図２に示した装置が行う画像処理の
フローチャートである。
フローチャートである。
【図９】注目画素と局所領域の関係を示す図である。
１００，１００’ マグネットシステム １０２ 主磁場コイル部 １０２’ 主磁場マグネット部 １０６，１０６’ 勾配コイル部 １０８，１０８’ ＲＦコイル部 １３０ 勾配駆動部 １４０ ＲＦ駆動部 １５０ データ収集部 １６０ 制御部 １７０ データ処理部 １８０ 表示部 １９０ 操作部 ３００ 対象 ５００ クレードル
フロントページの続き (51)Int.Cl.7 識別記号 ＦＩ テーマコート゛(参考） Ｈ０４Ｎ 5/21 Ｈ０４Ｎ 1/40 １０１Ｃ (72)発明者 荻野 徹男 東京都日野市旭が丘四丁目７番地の127 ジーイー横河メディカルシステム株式会社 内 Ｆターム(参考） 4C096 AB07 AD14 DC05 DC40 5B057 AA07 BA06 CE02 CE06 DA20 DC05 DC30 5C021 PA56 PA66 PA67 PA72 PA74 PA77 YA01 ZA00 ZA01 5C077 LL02 MM30 PQ19 5L096 AA06 BA06 BA13 CA18 FA35 FA70
Claims (36)
- 【請求項１】 画像の全体にわたって設定した複数の局
所領域ごとに画素値の残差平方和を求め、 前記残差平方和のヒストグラムを求め、 前記ヒストグラムのピークを与える残差平方和に基づい
てノイズの分散を求める、ことを特徴とする画像処理方
法。 - 【請求項２】 画像の全体にわたって設定した複数の局
所領域ごとに画素値の残差平方和を求め、 前記残差平方和のヒストグラムを求め、 前記ヒストグラムを関数にフィッティングし、 前記フィッティングした関数のピーク値を与える変数に
基づいてノイズの分散を求める、ことを特徴とする画像
処理方法。 - 【請求項３】 前記関数はガウス分布関数である、こと
を特徴とする請求項２に記載の画像処理方法。 - 【請求項４】 前記関数はレイリー分布関数である、こ
とを特徴とする請求項２に記載の画像処理方法。 - 【請求項５】 前記関数はガウス分布関数およびレイリ
ー分布関数である、ことを特徴とする請求項２に記載の
画像処理方法。 - 【請求項６】 前記画像の画素値は複素数の絶対値であ
る、ことを特徴とする請求項１ないし請求項５のうちの
いずれか１つに記載の画像処理方法。 - 【請求項７】 前記画像の画素値は複素数の実数部であ
る、ことを特徴とする請求項１ないし請求項３のうちの
いずれか１つに記載の画像処理方法。 - 【請求項８】 前記画像の画素値は複素数の虚数部であ
る、ことを特徴とする請求項１ないし請求項３のうちの
いずれか１つに記載の画像処理方法。 - 【請求項９】 前記画像は磁気共鳴を利用して撮影した
画像である、ことを特徴とする請求項１ないし請求項８
のうちのいずれか１つに記載の画像処理方法。 - 【請求項１０】 画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求める残差平方和
計算手段と、 前記残差平方和のヒストグラムを求めるヒストグラム計
算手段と、 前記ヒストグラムのピークを与える残差平方和に基づい
てノイズの分散を求めるノイズ分散計算手段と、を具備
することを特徴とする画像処理装置。 - 【請求項１１】 画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求める残差平方和
計算手段と、 前記残差平方和のヒストグラムを求めるヒストグラム計
算手段と、 前記ヒストグラムを関数にフィッティングするフィッテ
ィング手段と、 前記フィッティングした関数のピーク値を与える変数に
基づいてノイズの分散を求めるノイズ分散計算手段と、
を具備することを特徴とする画像処理装置。 - 【請求項１２】 前記関数はガウス分布関数である、こ
とを特徴とする請求項１１に記載の画像処理装置。 - 【請求項１３】 前記関数はレイリー分布関数である、
ことを特徴とする請求項１１に記載の画像処理装置。 - 【請求項１４】 前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする請求項１１に記
載の画像処理装置。 - 【請求項１５】 前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする請求項１０ないし請求項１４の
うちのいずれか１つに記載の画像処理装置。 - 【請求項１６】 前記画像の画素値は複素数の実数部で
ある、ことを特徴とする請求項１０ないし請求項１２の
うちのいずれか１つに記載の画像処理装置。 - 【請求項１７】 前記画像の画素値は複素数の虚数部で
ある、ことを特徴とする請求項１０ないし請求項１２の
うちのいずれか１つに記載の画像処理装置。 - 【請求項１８】 前記画像は磁気共鳴を利用して撮影し
た画像である、ことを特徴とする請求項１０ないし請求
項１７のうちのいずれか１つに記載の画像処理装置。 - 【請求項１９】 画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求める残差平方和
計算機能と、 前記残差平方和のヒストグラムを求めるヒストグラム計
算機能と、 前記ヒストグラムのピークを与える残差平方和に基づい
てノイズの分散を求めるノイズ分散計算機能と、をコン
ピュータに実現させるプログラムをコンピュータにより
読み取り可能なように記録したことを特徴とする記録媒
体。 - 【請求項２０】 画像の全体にわたって設定した複数の
局所領域ごとに画素値の残差平方和を求める残差平方和
計算機能と、 前記残差平方和のヒストグラムを求めるヒストグラム計
算機能と、 前記ヒストグラムを関数にフィッティングするフィッテ
ィング機能と、 前記フィッティングした関数のピーク値を与える変数に
基づいてノイズの分散を求めるノイズ分散計算機能と、
をコンピュータに実現させるプログラムをコンピュータ
により読み取り可能なように記録したことを特徴とする
記録媒体。 - 【請求項２１】 前記関数はガウス分布関数である、こ
とを特徴とする請求項２０に記載の記録媒体。 - 【請求項２２】 前記関数はレイリー分布関数である、
ことを特徴とする請求項２０に記載の記録媒体。 - 【請求項２３】 前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする請求項２０に記
載の記録媒体。 - 【請求項２４】 前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする請求項１９ないし請求項２３の
うちのいずれか１つに記載の記録媒体。 - 【請求項２５】 前記画像の画素値は複素数の実数部で
ある、ことを特徴とする請求項１９ないし請求項２１の
うちのいずれか１つに記載の記録媒体。 - 【請求項２６】 前記画像の画素値は複素数の虚数部で
ある、ことを特徴とする請求項１９ないし請求項２１の
うちのいずれか１つに記載の記録媒体。 - 【請求項２７】 前記画像は磁気共鳴を利用して撮影し
た画像である、ことを特徴とする請求項１９ないし請求
項２６のうちのいずれか１つに記載の記録媒体。 - 【請求項２８】 対象から収集した信号に基づいて画像
を生成する画像撮影装置であって、 画像の全体にわたって設定した複数の局所領域ごとに画
素値の残差平方和を求める残差平方和計算手段と、 前記残差平方和のヒストグラムを求めるヒストグラム計
算手段と、 前記ヒストグラムのピークを与える残差平方和に基づい
てノイズの分散を求めるノイズ分散計算手段と、を具備
することを特徴とする画像撮影装置。 - 【請求項２９】 対象から収集した信号に基づいて画像
を生成する画像撮影装置であって、 画像の全体にわたって設定した複数の局所領域ごとに画
素値の残差平方和を求める残差平方和計算手段と、 前記残差平方和のヒストグラムを求めるヒストグラム計
算手段と、 前記ヒストグラムを関数にフィッティングするフィッテ
ィング手段と、 前記フィッティングした関数のピーク値を与える変数に
基づいてノイズの分散を求めるノイズ分散計算手段と、
を具備することを特徴とする画像撮影装置。 - 【請求項３０】 前記関数はガウス分布関数である、こ
とを特徴とする請求項２９に記載の画像撮影装置。 - 【請求項３１】 前記関数はレイリー分布関数である、
ことを特徴とする請求項２９に記載の画像撮影装置。 - 【請求項３２】 前記関数はガウス分布関数およびレイ
リー分布関数である、ことを特徴とする請求項２９に記
載の画像撮影装置。 - 【請求項３３】 前記画像の画素値は複素数の絶対値で
ある、ことを特徴とする請求項２８ないし請求項３２の
うちのいずれか１つに記載の画像撮影装置。 - 【請求項３４】 前記画像の画素値は複素数の実数部で
ある、ことを特徴とする請求項２８ないし請求項３０の
うちのいずれか１つに記載の画像撮影装置。 - 【請求項３５】 前記画像の画素値は複素数の虚数部で
ある、ことを特徴とする請求項２８ないし請求項３０の
うちのいずれか１つに記載の画像撮影装置。 - 【請求項３６】 前記信号は磁気共鳴信号である、こと
を特徴とする請求項２８ないし請求項３５のうちのいず
れか１つに記載の画像撮影装置。
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000183049A JP4090671B2 (ja) | 2000-06-19 | 2000-06-19 | 画像処理方法、画像処理装置および画像撮影装置 |
| US09/846,631 US6804383B2 (en) | 2000-06-19 | 2001-05-01 | Image processing method and apparatus, recording medium, and imaging apparatus |
| EP01305292A EP1168244A3 (en) | 2000-06-19 | 2001-06-18 | Image processing method and apparatus, recording medium, and imaging apparatus |
| KR10-2001-0034301A KR100444094B1 (ko) | 2000-06-19 | 2001-06-18 | 이미지 처리 장치 및 그 방법, 촬상 장치 및 기록 매체 |
| CNB011219661A CN1237487C (zh) | 2000-06-19 | 2001-06-19 | 图像处理方法和装置、成像装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000183049A JP4090671B2 (ja) | 2000-06-19 | 2000-06-19 | 画像処理方法、画像処理装置および画像撮影装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002024820A true JP2002024820A (ja) | 2002-01-25 |
| JP4090671B2 JP4090671B2 (ja) | 2008-05-28 |
Family
ID=18683692
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000183049A Expired - Fee Related JP4090671B2 (ja) | 2000-06-19 | 2000-06-19 | 画像処理方法、画像処理装置および画像撮影装置 |
Country Status (5)
| Country | Link |
|---|---|
| US (1) | US6804383B2 (ja) |
| EP (1) | EP1168244A3 (ja) |
| JP (1) | JP4090671B2 (ja) |
| KR (1) | KR100444094B1 (ja) |
| CN (1) | CN1237487C (ja) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005012771A (ja) * | 2003-05-14 | 2005-01-13 | Ge Medical Systems Global Technology Co Llc | 放射線撮像画像のコントラスト及び輝度を設定するための方法及び装置 |
| JP2006340992A (ja) * | 2005-06-10 | 2006-12-21 | Hitachi Medical Corp | 画像処理プログラム、装置及び方法 |
| JP2010074812A (ja) * | 2008-08-18 | 2010-04-02 | Aloka Co Ltd | 画像処理装置 |
| JP2020103365A (ja) * | 2018-12-26 | 2020-07-09 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置および医用複素数画像処理装置 |
| JP2022118536A (ja) * | 2021-02-02 | 2022-08-15 | 富士フイルムヘルスケア株式会社 | 医用画像処理装置および医用画像処理方法 |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4344964B2 (ja) * | 1999-06-01 | 2009-10-14 | ソニー株式会社 | 画像処理装置および画像処理方法 |
| US7015401B2 (en) * | 2001-11-23 | 2006-03-21 | Aiptek International, Inc. | Image processing system with handwriting input function and the method for forming the same |
| JP4407243B2 (ja) * | 2003-11-10 | 2010-02-03 | ソニー株式会社 | 照合処理装置とその方法 |
| KR101097512B1 (ko) * | 2004-11-23 | 2011-12-22 | 엘지디스플레이 주식회사 | 액정표시장치 및 구동방법 |
| CN1328901C (zh) * | 2005-01-26 | 2007-07-25 | 北京中星微电子有限公司 | 一种去除图像噪声的方法 |
| CN101283286A (zh) * | 2005-10-07 | 2008-10-08 | 皇家飞利浦电子股份有限公司 | 多信道发射磁共振 |
| US7660481B2 (en) | 2005-11-17 | 2010-02-09 | Vital Images, Inc. | Image enhancement using anisotropic noise filtering |
| TWI319676B (en) * | 2006-10-18 | 2010-01-11 | Quanta Comp Inc | Image processing apparatus and method |
| DE102008023915A1 (de) * | 2008-05-16 | 2009-12-10 | Siemens Aktiengesellschaft | Verfahren zur Einstellung von wenigstens einer Stellgröße eines Entrauschungsfilters in medizinischen Bildern |
| DE102008058740B4 (de) * | 2008-11-24 | 2012-05-31 | Siemens Aktiengesellschaft | Korrektur von Artefakten in der "Time of Flight"- MR-Angiographie |
| CN111311525A (zh) * | 2019-11-20 | 2020-06-19 | 重庆邮电大学 | 一种基于直方图概率修正的图像梯度场双区间均衡化算法 |
| CN114859352B (zh) * | 2022-07-07 | 2022-10-28 | 自然资源部第二海洋研究所 | Sar卫星海洋观测图像自适应拉伸方法及装置 |
Family Cites Families (14)
| 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 |
| JPH02294786A (ja) * | 1989-05-09 | 1990-12-05 | Hitachi Ltd | 画像輝度レベル値の自動正規化方法 |
| EP0574969B1 (en) * | 1992-06-19 | 2002-03-20 | Agfa-Gevaert | A method and apparatus for noise reduction |
| JP2800633B2 (ja) * | 1993-04-30 | 1998-09-21 | 富士ゼロックス株式会社 | 画像符号化装置 |
| JP3410778B2 (ja) * | 1993-09-14 | 2003-05-26 | 株式会社東芝 | リカーシブフィルタ、ｘ線診断装置 |
| US5803082A (en) * | 1993-11-09 | 1998-09-08 | Staplevision Inc. | Omnispectramammography |
| US5594807A (en) * | 1994-12-22 | 1997-01-14 | Siemens Medical Systems, Inc. | System and method for adaptive filtering of images based on similarity between histograms |
| JPH119595A (ja) * | 1997-06-27 | 1999-01-19 | Aloka Co Ltd | 超音波診断装置 |
| US5876342A (en) * | 1997-06-30 | 1999-03-02 | Siemens Medical Systems, Inc. | System and method for 3-D ultrasound imaging and motion estimation |
| JPH1176232A (ja) * | 1997-09-11 | 1999-03-23 | Hitachi Medical Corp | 超音波診断装置 |
| JP2000126182A (ja) * | 1998-10-27 | 2000-05-09 | Mitani Sangyo Co Ltd | 腫瘍診断方法 |
| JP4313869B2 (ja) * | 1998-11-12 | 2009-08-12 | 株式会社東芝 | 超音波診断装置 |
| JP4344964B2 (ja) * | 1999-06-01 | 2009-10-14 | ソニー株式会社 | 画像処理装置および画像処理方法 |
| KR100317674B1 (ko) * | 1999-06-19 | 2001-12-24 | 윤종용 | 순환 노이즈 제거 시스템 |
-
2000
- 2000-06-19 JP JP2000183049A patent/JP4090671B2/ja not_active Expired - Fee Related
-
2001
- 2001-05-01 US US09/846,631 patent/US6804383B2/en not_active Expired - Fee Related
- 2001-06-18 KR KR10-2001-0034301A patent/KR100444094B1/ko not_active IP Right Cessation
- 2001-06-18 EP EP01305292A patent/EP1168244A3/en not_active Withdrawn
- 2001-06-19 CN CNB011219661A patent/CN1237487C/zh not_active Expired - Fee Related
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2005012771A (ja) * | 2003-05-14 | 2005-01-13 | Ge Medical Systems Global Technology Co Llc | 放射線撮像画像のコントラスト及び輝度を設定するための方法及び装置 |
| JP2006340992A (ja) * | 2005-06-10 | 2006-12-21 | Hitachi Medical Corp | 画像処理プログラム、装置及び方法 |
| JP4707471B2 (ja) * | 2005-06-10 | 2011-06-22 | 株式会社日立メディコ | 画像処理プログラム、装置及び方法 |
| JP2010074812A (ja) * | 2008-08-18 | 2010-04-02 | Aloka Co Ltd | 画像処理装置 |
| JP2020103365A (ja) * | 2018-12-26 | 2020-07-09 | キヤノンメディカルシステムズ株式会社 | 磁気共鳴イメージング装置および医用複素数画像処理装置 |
| JP2022118536A (ja) * | 2021-02-02 | 2022-08-15 | 富士フイルムヘルスケア株式会社 | 医用画像処理装置および医用画像処理方法 |
| JP7395523B2 (ja) | 2021-02-02 | 2023-12-11 | 富士フイルムヘルスケア株式会社 | 医用画像処理装置および医用画像処理方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1329325A (zh) | 2002-01-02 |
| US6804383B2 (en) | 2004-10-12 |
| EP1168244A2 (en) | 2002-01-02 |
| US20020009216A1 (en) | 2002-01-24 |
| KR20010113535A (ko) | 2001-12-28 |
| EP1168244A3 (en) | 2004-12-15 |
| JP4090671B2 (ja) | 2008-05-28 |
| KR100444094B1 (ko) | 2004-08-09 |
| CN1237487C (zh) | 2006-01-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP4149126B2 (ja) | 画像処理方法、画像処理装置および画像撮影装置 | |
| JP2002024820A (ja) | 画像処理方法および装置、記録媒体並びに画像撮影装置 | |
| JP4191884B2 (ja) | 画像処理方法、画像処理装置および画像撮影装置 | |
| JPH074352B2 (ja) | 核磁気共鳴装置の操作方法 | |
| EP2730935A1 (en) | Motional error correction in functional magnetic resonance imaging | |
| JP2003111746A (ja) | Ｔ２／ｔ２＊の重み付けを低減した中位の磁場強度における水分および脂肪のｍｒ画像の区別および識別 | |
| JP3454760B2 (ja) | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 | |
| JP4180936B2 (ja) | 磁気共鳴撮影装置 | |
| JP2003019124A (ja) | スピン励起方法、磁気共鳴撮影方法および磁気共鳴撮影装置 | |
| US6486667B1 (en) | Combination of fluid-attenuated inversion-recovery complex images acquired using magnetic resonance imaging | |
| JP3884227B2 (ja) | 磁気共鳴撮影装置 | |
| JPH0687847B2 (ja) | 磁気共鳴映像装置 | |
| JP2002102201A (ja) | 磁気共鳴信号獲得方法および装置、記録媒体並びに磁気共鳴撮影装置 | |
| JP2002042128A (ja) | 画像処理方法および装置、記録媒体並びに画像撮影装置 | |
| JP2001299722A (ja) | 勾配コイル製造方法および勾配コイル並びに磁気共鳴撮影装置 | |
| KR101502103B1 (ko) | 자기공명영상장치 및 자화강조영상법 | |
| JP2001218749A (ja) | 磁気共鳴撮影装置および記録媒体 | |
| JP4363625B2 (ja) | 画像処理方法および磁気共鳴撮影装置 | |
| JP2004261487A (ja) | 磁気共鳴撮影装置 | |
| JP2002052005A (ja) | 磁気共鳴イメージング方法 | |
| JP2003144410A (ja) | 液量計測装置 | |
| JPH07255697A (ja) | 磁気共鳴イメージングの撮影方法及びその装置 | |
| JP4558866B2 (ja) | 位相分布測定方法および装置、位相補正方法および装置、並びに、磁気共鳴撮像装置 | |
| JPH05123314A (ja) | 磁気共鳴イメージング装置におけるマルチスライス撮像方法 | |
| JP2004254884A (ja) | 磁気共鳴撮影装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A625 | Written request for application examination (by other person) |
Free format text: JAPANESE INTERMEDIATE CODE: A625 Effective date: 20040714 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070517 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20070529 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070828 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20071002 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20071227 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080205 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080227 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110307 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120307 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120307 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120307 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130307 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130307 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130307 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140307 Year of fee payment: 6 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |