JP2002042128A - 画像処理方法および装置、記録媒体並びに画像撮影装置 - Google Patents

画像処理方法および装置、記録媒体並びに画像撮影装置

Info

Publication number
JP2002042128A
JP2002042128A JP2000220416A JP2000220416A JP2002042128A JP 2002042128 A JP2002042128 A JP 2002042128A JP 2000220416 A JP2000220416 A JP 2000220416A JP 2000220416 A JP2000220416 A JP 2000220416A JP 2002042128 A JP2002042128 A JP 2002042128A
Authority
JP
Japan
Prior art keywords
pixel
pixel group
value
group
pixel value
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
Application number
JP2000220416A
Other languages
English (en)
Other versions
JP2002042128A5 (ja
JP4012669B2 (ja
Inventor
Tetsuo Ogino
徹男 荻野
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
GE Medical Systems Global Technology Co LLC
Original Assignee
GE Medical Systems Global Technology Co LLC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by GE Medical Systems Global Technology Co LLC filed Critical GE Medical Systems Global Technology Co LLC
Priority to JP2000220416A priority Critical patent/JP4012669B2/ja
Priority to US09/863,998 priority patent/US7061651B2/en
Priority to TW090116777A priority patent/TWI221263B/zh
Priority to EP01306159A priority patent/EP1176555B1/en
Priority to DE60133694T priority patent/DE60133694T2/de
Priority to KR10-2001-0043809A priority patent/KR100452644B1/ko
Priority to CNB011253886A priority patent/CN1262237C/zh
Publication of JP2002042128A publication Critical patent/JP2002042128A/ja
Publication of JP2002042128A5 publication Critical patent/JP2002042128A5/ja
Application granted granted Critical
Publication of JP4012669B2 publication Critical patent/JP4012669B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • 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
    • 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

Landscapes

  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Biophysics (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Veterinary Medicine (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Image Processing (AREA)
  • Facsimile Image Signal Circuits (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Image Analysis (AREA)
  • Editing Of Facsimile Originals (AREA)

Abstract

(57)【要約】 【課題】 フィルタリングを適正に行う画像処理方法お
よび装置、フィルタリングプログラムを記録した媒体、
並びに、フィルタリングを適正に行う画像撮影装置を実
現する。 【解決手段】 元画像において注目画素を含む局所領域
を設定し（５０４）、その領域において画素グループを
複数の態様で設定し（５０６）、局所領域における元画
像の構造に最も良く適合する画素グループ態様を選択し
（５０８，５１０）、注目画素を含む画素グループの平
均画素値を注目画素の新たな画素値とする。

Description

【発明の詳細な説明】
【０００１】
【発明の属する技術分野】本発明は、画像処理方法およ
び装置、記録媒体並びに画像撮影装置に関し、特に、画
像のノイズ（ｎｏｉｓｅ）を除去する画像処理方法およ
び装置、そのような画像処理機能をコンピュータに実現
させるプログラムを記録した媒体、並びに、そのような
画像処理装置を備えた画像撮影装置に関する。
【０００２】
【従来の技術】磁気共鳴撮影（ＭＲＩ：Ｍａｇｎｅｔｉ
ｃ Ｒｅｓｏｎａｎｃｅ Ｉｍａｇｉｎｇ）装置では、
マグネットシステム（ｍａｇｎｅｔ ｓｙｓｔｅｍ）の
内部空間、すなわち、静磁場を形成した空間に撮影の対
象を搬入し、勾配磁場および高周波磁場を印加して対象
内に磁気共鳴信号を発生させ、その受信信号に基づいて
断層像を生成（再構成）する。
【０００３】断層像の細部構造をよりよく観察できるよ
うにするために、画像のノイズ（ｎｏｉｓｅ）を除去す
るフィルタリング（ｆｉｌｔｅｒｉｎｇ）が行われる。
フィルタリングはローパスフィルタリング（ｌｏｗ−ｐ
ａｓｓ ｆｉｌｔｅｒｉｎｇ）を基本とするが、それだ
けでは画像の鮮鋭度（シャープネス：ｓｈａｒｐｎｅｓ
ｓ）が低下するので、シャープネスを確保するための処
理を付加したフィルタリングが採用される。
【０００４】
【発明が解決しようとする課題】シャープネスを確保す
るための処理を付加したフィルタリングでは、ノイズが
偶然に形成したに過ぎないテクスチャ（ｔｅｘｔｕｒ
ｅ）を強調し、解剖学的に意味のない構造（偽構造）を
作り出すという副作用がある。
【０００５】そこで、本発明の課題は、フィルタリング
を適正に行う画像処理方法および装置、そのような画像
処理機能をコンピュータに実現させるプログラムを記録
した媒体、並びに、そのような画像処理装置を備えた画
像撮影装置を実現することである。
【０００６】
【課題を解決するための手段】（１）上記の課題を解決
するための１つの観点での発明は、元画像において注目
画素を含む局所的な領域を設定し、前記領域において前
記注目画素を含む複数の画素からなる画素グループおよ
びこの画素グループの画素と重複しない複数の画素から
なる画素グループを複数の態様で設定し、前記設定した
複数の画素グループ態様のうち前記領域における前記元
画像の構造に最も良く適合する画素グループ態様を選択
し、前記選択した画素グループ態様における前記注目画
素を含む画素グループの平均画素値を前記注目画素の新
たな画素値として画像を形成する、ことを特徴とする画
像処理方法である。
【０００７】この観点での発明では、元画像の構造に最
も良く適合する画素グループ態様を選択し、その画素グ
ループ態様において、注目画素を含む画素グループの平
均画素値を注目画素の新たな画素値とするので、偽構造
が発生しにくいフィルタリングを行うことができる。
【０００８】（２）上記の課題を解決するための他の観
点での発明は、元画像において注目画素を含む局所的な
領域を設定し、前記領域において前記注目画素を含む複
数の画素からなる画素グループおよびこの画素グループ
の画素と重複しない複数の画素からなる画素グループを
複数の態様で設定し、前記設定した複数の画素グループ
態様のうち前記領域における前記元画像の構造に最も良
く適合する画素グループ態様を選択し、前記選択した画
素グループ態様における前記注目画素を含む画素グルー
プの平均画素値を求め、前記領域の平均画素値を求め、
前記画素グループの平均画素値と前記領域の平均画素値
とを重み付け加算した画素値を前記注目画素の新たな画
素値として画像を形成する、ことを特徴とする画像処理
方法である。
【０００９】この観点での発明では、（１）と同様にし
て求めた画素グループの平均画素値に、領域の平均画素
値を重み付け加算して注目画素の新たな画素値とするの
で、さらに偽構造が発生しにくいフィルタリングを行う
ことができる。
【００１０】（３）上記の課題を解決するための他の観
点での発明は、元画像において注目画素を含む局所的な
領域を設定し、前記領域において前記注目画素を含む複
数の画素からなる画素グループおよびこの画素グループ
の画素と重複しない複数の画素からなる画素グループを
複数の態様で設定し、前記設定した複数の画素グループ
態様のうち前記領域における前記元画像の構造に最も良
く適合する画素グループ態様を選択し、前記選択した画
素グループ態様における前記注目画素を含む画素グルー
プの平均画素値を求め、前記領域の平均画素値を求め、
前記画素グループの平均画素値と前記領域の平均画素値
とを重み付け加算し、前記重み付け加算して得られた画
素値と前記注目画素の画素値とを重み付け加算した画素
値を前記注目画素の新たな画素値として画像を形成す
る、ことを特徴とする画像処理方法である。
【００１１】この観点での発明では、（２）と同様にし
て求めた画素値に元画像の注目画素の画素値を重み付け
加算して注目画素の新たな画素値とするので、元画像の
特徴的な構造を損なわないフィルタリングを行うことが
できる。
【００１２】（４）上記の課題を解決するための他の観
点での発明は、元画像において注目画素を含む局所的な
領域を設定し、前記領域において前記注目画素を含む複
数の画素からなる画素グループおよびこの画素グループ
の画素と重複しない複数の画素からなる画素グループを
複数の態様で設定し、前記設定した複数の画素グループ
態様のうち前記領域における前記元画像の構造に最も良
く適合する画素グループ態様を選択し、前記選択した画
素グループ態様における前記注目画素を含む画素グルー
プの平均画素値を求め、前記領域の平均画素値を求め、
前記画素グループの平均画素値と前記領域の平均画素値
とを重み付け加算し、前記重み付け加算して得られた画
素値と前記注目画素の画素値とを重み付け加算した画素
値を前記注目画素の新たな画素値として画像を形成し、
前記形成した画像と前記元画像を重み付け加算する、こ
とを特徴とする画像処理方法である。
【００１３】この観点での発明では、（３）と同様にし
て求めた画素値によって形成した画像に元画像を重み付
け加算するので、フィルタリングの効果を弱めることが
できる。
【００１４】（５）上記の課題を解決するための他の観
点での発明は、前記形成した画像と前記元画像を重み付
け加算する重み係数は調節可能である、ことを特徴とす
る（４）に記載の画像処理方法である。
【００１５】この観点での発明では、重み係数の調整に
よりフィルタリングの効果を調節することができる。
（６）上記の課題を解決するための他の観点での発明
は、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する重み係数は、前記画素グルー
プ態様ごとに求めた前記画素グループの画素値の残差平
方和の総和の最小値および最大値の関数である、ことを
特徴とする（２）ないし（５）のうちのいずれか１つに
記載の画像処理方法である。
【００１６】この観点での発明では、（２）ないし
（５）のうちのいずれか１つに記載の画像処理におい
て、画素グループの平均画素値と領域の平均画素値とを
重み付け加算する重み係数を、画素グループ態様ごとに
求めた画素グループの画素値の残差平方和の総和の最小
値および最大値の関数としたので、偽構造の発生しにく
いフィルタリングを行うことができる。
【００１７】（７）上記の課題を解決するための他の観
点での発明は、前記関数は、前記残差平方和の総和の最
小値と最大値が等しいとき前記領域の平均画素値に対す
る重みを最大とし、前記残差平方和の総和の最小値が最
大値に対して小さくなるほど前記領域の平均画素値に対
する重みを小さくする関数である、ことを特徴とする
（６）に記載の画像処理方法である。
【００１８】この観点での発明では、（６）に記載の画
像処理において、関数は、残差平方和の総和の最小値と
最大値が等しいとき領域の平均画素値に対する重みを最
大とし、残差平方和の総和の最小値が最大値に対して小
さくなるほど領域の平均画素値に対する重みを小さくす
る関数であるので、偽構造の発生しにくいフィルタリン
グを行うことができる。
【００１９】（８）上記の課題を解決するための他の観
点での発明は、前記重み付け加算して得られた画素値と
前記注目画素の画素値とを重み付け加算する重み係数
は、前記画素グループ態様ごとに求めた前記画素グルー
プの画素値の残差平方和の総和の最小値および前記元画
像のノイズの分散の関数である、ことを特徴とする
（３）ないし（７）のうちのいずれか１つに記載の画像
処理方法である。
【００２０】この観点での発明では、（３）ないし
（７）のうちのいずれか１つに記載の画像処理におい
て、重み付け加算して得られた画素値と元画像の注目画
素の画素値とを重み付け加算する重み係数を、画素グル
ープ態様ごとに求めた画素グループの画素値の残差平方
和の総和の最小値および元画像のノイズの分散の関数と
したので、元画像の特徴的な構造を損なわないフィルタ
リングを行うことができる。
【００２１】（９）上記の課題を解決するための他の観
点での発明は、前記関数は、前記残差平方和の総和の最
小値を前記画素グループ態様における画素数で割ったも
のと前記ノイズの分散が等しいとき前記重み付け加算し
て得られた画素値に対する重みを最大とし、前記残差平
方和の総和の最小値を前記画素グループ態様における画
素数で割ったものと前記ノイズの分散の違いが大きくな
るほど前記重み付け加算して得られた画素値に対する重
みを小さくする関数である、ことを特徴とする（８）に
記載の画像処理方法である。
【００２２】この観点での発明では、（８）に記載の画
像処理において、関数は、残差平方和の総和の最小値を
画素グループ態様における画素数で割ったものとノイズ
の分散が等しいとき重み付け加算して得られた画素値に
対する重みを最大とし、残差平方和の総和の最小値を画
素グループ態様における画素数で割ったものとノイズの
分散の違いが大きくなるほど重み付け加算して得られた
画素値に対する重みを小さくする関数であるので、元画
像の特徴的な構造を損なわないフィルタリングを行うこ
とができる。
【００２３】（１０）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が最小となる
画素グループ態様を選択することによって行う、ことを
特徴とする（１）ないし（９）のうちのいずれか１つに
記載の画像処理方法である。
【００２４】この観点での発明では、（１）ないし
（９）のうちのいずれか１つに記載の画像処理における
画素グループ態様の選択は、画素グループの画素値の残
差平方和の総和が最小となる画素グループ態様を選択す
ることによって行うので、元画像の構造に最も良く適合
する画素グループ態様を選択することができる。
【００２５】（１１）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記元画像
のノイズの分散に画素グループ態様における画素数を掛
けたものに最も近くなる画素グループ態様を選択するこ
とによって行う、ことを特徴とする（１）ないし（９）
のうちのいずれか１つに記載の画像処理方法である。
【００２６】この観点での発明では、（１）ないし
（９）のうちのいずれか１つに記載の画像処理における
画素グループ態様の選択は、画素グループの画素値の残
差平方和の総和が元画像のノイズの分散に画素グループ
態様における画素数を掛けたものに最も近くなる画素グ
ループ態様を選択することによって行うので、元画像の
構造に最も良く適合する画素グループ態様を選択するこ
とができる。
【００２７】（１２）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定し、前記設定した複数態様の
領域において、それぞれ、前記注目画素を含む複数の画
素からなる画素グループおよびこの画素グループの画素
と重複しない複数の画素からなる画素グループを複数の
態様で設定し、前記設定した複数の画素グループ態様の
うち前記設定した複数態様の領域を通じてその領域にお
ける前記元画像の構造に最も良く適合する画素グループ
態様を選択し、前記選択した画素グループ態様における
前記注目画素を含む画素グループの平均画素値を前記注
目画素の新たな画素値として画像を形成する、ことを特
徴とする画像処理方法である。
【００２８】この観点での発明では、複数態様の領域を
通じてその領域における元画像の構造に最も良く適合す
る画素グループ態様を選択し、その画素グループ態様に
おいて、注目画素を含む画素グループの平均画素値を注
目画素の新たな画素値とするので、偽構造が発生しにく
いフィルタリングを行うことができる。
【００２９】（１３）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定し、前記設定した複数態様の
領域において、それぞれ、前記注目画素を含む複数の画
素からなる画素グループおよびこの画素グループの画素
と重複しない複数の画素からなる画素グループを複数の
態様で設定し、前記設定した複数の画素グループ態様の
うち前記設定した複数態様の領域を通じてその領域にお
ける前記元画像の構造に最も良く適合する画素グループ
態様を選択し、前記選択した画素グループ態様における
前記注目画素を含む画素グループの平均画素値を求め、
前記複数態様の領域のうち前記選択した画素グループ態
様が設定されている領域の平均画素値を求め、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算した画素値を前記注目画素の新たな画素値とし
て画像を形成する、ことを特徴とする画像処理方法であ
る。
【００３０】この観点での発明では、（１２）と同様に
して求めた画素グループの平均画素値に、領域の平均画
素値を重み付け加算して注目画素の新たな画素値とする
ので、さらに偽構造が発生しにくいフィルタリングを行
うことができる。
【００３１】（１４）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定し、前記設定した複数態様の
領域において、それぞれ、前記注目画素を含む複数の画
素からなる画素グループおよびこの画素グループの画素
と重複しない複数の画素からなる画素グループを複数の
態様で設定し、前記設定した複数の画素グループ態様の
うち前記設定した複数態様の領域を通じてその領域にお
ける前記元画像の構造に最も良く適合する画素グループ
態様を選択し、前記選択した画素グループ態様における
前記注目画素を含む画素グループの平均画素値を求め、
前記複数態様の領域のうち前記選択した画素グループ態
様が設定されている領域の平均画素値を求め、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算し、前記重み付け加算して得られた画素値と前
記注目画素の画素値とを重み付け加算した画素値を前記
注目画素の新たな画素値として画像を形成する、ことを
特徴とする画像処理方法である。
【００３２】この観点での発明では、（１３）と同様に
して求めた画素値に元画像の注目画素の画素値を重み付
け加算して注目画素の新たな画素値とするので、元画像
の特徴的な構造を損なわないフィルタリングを行うこと
ができる。
【００３３】（１５）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定し、前記設定した複数態様の
領域において、それぞれ、前記注目画素を含む複数の画
素からなる画素グループおよびこの画素グループの画素
と重複しない複数の画素からなる画素グループを複数の
態様で設定し、前記設定した複数の画素グループ態様の
うち前記設定した複数態様の領域を通じてその領域にお
ける前記元画像の構造に最も良く適合する画素グループ
態様を選択し、前記選択した画素グループ態様における
前記注目画素を含む画素グループの平均画素値を求め、
前記複数態様の領域のうち前記選択した画素グループ態
様が設定されている領域の平均画素値を求め、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算し、前記重み付け加算して得られた画素値と前
記注目画素の画素値とを重み付け加算した画素値を前記
注目画素の新たな画素値として画像を形成し、前記形成
した画像と前記元画像を重み付け加算する、ことを特徴
とする画像処理方法である。
【００３４】この観点での発明では、（１４）と同様に
して求めた画素値によって形成した画像に元画像を重み
付け加算するので、フィルタリングの効果を弱めること
ができる。
【００３５】（１６）上記の課題を解決するための他の
観点での発明は、前記形成した画像と前記元画像を重み
付け加算する重み係数は調節可能である、ことを特徴と
する（１５）に記載の画像処理方法である。
【００３６】この観点での発明では、重み係数の調整に
よりフィルタリングの効果を調節することができる。
（１７）上記の課題を解決するための他の観点での発明
は、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する重み係数は、前記設定した複
数態様の領域を通じて前記画像グループ態様ごとに求め
た前記画素グループの画素値の残差平方和の総和の最小
値および最大値の関数である、ことを特徴とする（１
３）ないし（１６）のうちのいずれか１つに記載の画像
処理方法である。
【００３７】この観点での発明では、（１３）ないし
（１６）のうちのいずれか１つに記載の画像処理におい
て、画素グループの平均画素値と領域の平均画素値とを
重み付け加算する重み係数を、複数態様の領域を通じて
画素グループ態様ごとに求めた画素グループの画素値の
残差平方和の総和の最小値および最大値の関数としたの
で、偽構造の発生しにくいフィルタリングを行うことが
できる。
【００３８】（１８）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値と最大値が等しいとき前記領域の平均画素値に対
する重みを最大とし、前記残差平方和の総和の最小値が
最大値に対して小さくなるほど前記領域の平均画素値に
対する重みを小さくする関数である、ことを特徴とする
（１７）に記載の画像処理方法である。
【００３９】この観点での発明では、（１７）に記載の
画像処理において、関数は、残差平方和の総和の最小値
と最大値が等しいとき領域の平均画素値に対する重みを
最大とし、残差平方和の総和の最小値が最大値に対して
小さくなるほど領域の平均画素値に対する重みを小さく
する関数であるので、偽構造の発生しにくいフィルタリ
ングを行うことができる。
【００４０】（１９）上記の課題を解決するための他の
観点での発明は、前記重み付け加算して得られた画素値
と前記注目画素の画素値とを重み付け加算する重み係数
は、前記設定した複数態様の領域を通じて前記画像グル
ープ態様ごとに求めた前記画素グループの画素値の残差
平方和の総和の最小値および前記元画像のノイズの分散
の関数である、ことを特徴とする（１４）ないし（１
８）のうちのいずれか１つに記載の画像処理方法であ
る。
【００４１】この観点での発明では、（１４）ないし
（１８）のうちのいずれか１つに記載の画像処理におい
て、重み付け加算して得られた画素値と元画像の注目画
素の画素値とを重み付け加算する重み係数を、複数態様
の領域を通じて画素グループ態様ごとに求めた画素グル
ープの画素値の残差平方和の総和の最小値および元画像
のノイズの分散の関数としたので、元画像の特徴的な構
造を損なわないフィルタリングを行うことができる。
【００４２】（２０）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値を前記画素グループ態様における画素数で割った
ものと前記ノイズの分散が等しいとき前記重み付け加算
して得られた画素値に対する重みを最大とし、前記残差
平方和の総和の最小値を前記画素グループ態様における
画素数で割ったものと前記ノイズの分散の違いが大きく
なるほど前記重み付け加算して得られた画素値に対する
重みを小さくする関数である、ことを特徴とする（１
９）に記載の画像処理方法である。
【００４３】この観点での発明では、（１９）に記載の
画像処理において、関数は、残差平方和の総和の最小値
を画素グループ態様における画素数で割ったものとノイ
ズの分散が等しいとき重み付け加算して得られた画素値
に対する重みを最大とし、残差平方和の総和の最小値を
画素グループ態様における画素数で割ったものとノイズ
の分散の違いが大きくなるほど重み付け加算して得られ
た画素値に対する重みを小さくする関数であるので、元
画像の特徴的な構造を損なわないフィルタリングを行う
ことができる。
【００４４】（２１）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記設定し
た複数態様の領域を通じて最小となる画像グループ態様
を選択することによって行う、ことを特徴とする（１
２）ないし（２０）のうちのいずれか１つに記載の画像
処理方法である。
【００４５】この観点での発明では、（１２）ないし
（２０）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が複数態様の領域を通じて最小となる
画素グループ態様を選択することによって行うので、元
画像の構造に最も良く適合する画素グループ態様を選択
することができる。
【００４６】（２２）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記元画像
のノイズの分散に画素グループ態様における画素数を掛
けたものに最も近くなる画素グループ態様を選択するこ
とによって行う、ことを特徴とする（１２）ないし（２
０）のうちのいずれか１つに記載の画像処理方法であ
る。
【００４７】この観点での発明では、（１２）ないし
（２０）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が元画像のノイズの分散に画素グルー
プ態様における画素数を掛けたものに最も近くなる画素
グループ態様を選択することによって行うので、元画像
の構造に最も良く適合する画素グループ態様を選択する
ことができる。
【００４８】（２３）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を設定する領域設定手段と、前記領域において前
記注目画素を含む複数の画素からなる画素グループおよ
びこの画素グループの画素と重複しない複数の画素から
なる画素グループを複数の態様で設定する画素グループ
設定手段と、前記設定した複数の画素グループ態様のう
ち前記領域における前記元画像の構造に最も良く適合す
る画素グループ態様を選択する選択手段と、前記選択し
た画素グループ態様における前記注目画素を含む画素グ
ループの平均画素値を前記注目画素の新たな画素値とし
て画像を形成する画像形成手段と、を具備することを特
徴とする画像処理装置である。
【００４９】この観点での発明では、元画像の構造に最
も良く適合する画素グループ態様を選択し、その画素グ
ループ態様において、注目画素を含む画素グループの平
均画素値を注目画素の新たな画素値とするので、偽構造
が発生しにくいフィルタリングを行うことができる。
【００５０】（２４）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を設定する領域設定手段と、前記領域において前
記注目画素を含む複数の画素からなる画素グループおよ
びこの画素グループの画素と重複しない複数の画素から
なる画素グループを複数の態様で設定する画素グループ
設定手段と、前記設定した複数の画素グループ態様のう
ち前記領域における前記元画像の構造に最も良く適合す
る画素グループ態様を選択する選択手段と、前記選択し
た画素グループ態様における前記注目画素を含む画素グ
ループの平均画素値を求める第１の画素値計算手段と、
前記領域の平均画素値を求める第２の画素値計算手段
と、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算した画素値を前記注目画素の新た
な画素値として画像を形成する画像形成手段と、を具備
することを特徴とする画像処理装置である。
【００５１】この観点での発明では、（２３）と同様に
して求めた画素グループの平均画素値に、領域の平均画
素値を重み付け加算して注目画素の新たな画素値とする
ので、さらに偽構造が発生しにくいフィルタリングを行
うことができる。
【００５２】（２５）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を設定する領域設定手段と、前記領域において前
記注目画素を含む複数の画素からなる画素グループおよ
びこの画素グループの画素と重複しない複数の画素から
なる画素グループを複数の態様で設定する画素グループ
設定手段と、前記設定した複数の画素グループ態様のう
ち前記領域における前記元画像の構造に最も良く適合す
る画素グループ態様を選択する選択手段と、前記選択し
た画素グループ態様における前記注目画素を含む画素グ
ループの平均画素値を求める第１の画素値計算手段と、
前記領域の平均画素値を求める第２の画素値計算手段
と、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する加算手段と、前記重み付け加
算して得られた画素値と前記注目画素の画素値とを重み
付け加算した画素値を前記注目画素の新たな画素値とし
て画像を形成する画像形成手段と、を具備することを特
徴とする画像処理装置である。
【００５３】この観点での発明では、（２４）と同様に
して求めた画素値に元画像の注目画素の画素値を重み付
け加算して注目画素の新たな画素値とするので、元画像
の特徴的な構造を損なわないフィルタリングを行うこと
ができる。
【００５４】（２６）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を設定する領域設定手段と、前記領域において前
記注目画素を含む複数の画素からなる画素グループおよ
びこの画素グループの画素と重複しない複数の画素から
なる画素グループを複数の態様で設定する画素グループ
設定手段と、前記設定した複数の画素グループ態様のう
ち前記領域における前記元画像の構造に最も良く適合す
る画素グループ態様を選択する選択手段と、前記選択し
た画素グループ態様における前記注目画素を含む画素グ
ループの平均画素値を求める第１の画素値計算手段と、
前記領域の平均画素値を求める第２の画素値計算手段
と、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する第１の加算手段と、前記重み
付け加算して得られた画素値と前記注目画素の画素値と
を重み付け加算した画素値を前記注目画素の新たな画素
値として画像を形成する画像形成手段と、前記形成した
画像と前記元画像を重み付け加算する第２の加算手段
と、を具備することを特徴とする画像処理装置である。
【００５５】この観点での発明では、（２５）と同様に
して求めた画素値によって形成した画像に元画像を重み
付け加算するので、フィルタリングの効果を弱めること
ができる。
【００５６】（２７）上記の課題を解決するための他の
観点での発明は、前記形成した画像と前記元画像を重み
付け加算する重み係数は調節可能である、ことを特徴と
する（２６）に記載の画像処理装置である。
【００５７】この観点での発明では、重み係数の調整に
よりフィルタリングの効果を調節することができる。 （２８）上記の課題を解決するための他の観点での発明
は、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する重み係数は、前記画素グルー
プ態様ごとに求めた前記画素グループの画素値の残差平
方和の総和の最小値および最大値の関数である、ことを
特徴とする（２４）ないし（２７）のうちのいずれか１
つに記載の画像処理装置である。
【００５８】この観点での発明では、（２４）ないし
（２７）のうちのいずれか１つに記載の画像処理におい
て、画素グループの平均画素値と領域の平均画素値とを
重み付け加算する重み係数を、画素グループ態様ごとに
求めた画素グループの画素値の残差平方和の総和の最小
値および最大値の関数としたので、偽構造の発生しにく
いフィルタリングを行うことができる。
【００５９】（２９）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値と最大値が等しいとき前記領域の平均画素値に対
する重みを最大とし、前記残差平方和の総和の最小値が
最大値に対して小さくなるほど前記領域の平均画素値に
対する重みを小さくする関数である、ことを特徴とする
（２８）に記載の画像処理装置である。
【００６０】この観点での発明では、（２８）に記載の
画像処理において、関数は、残差平方和の総和の最小値
と最大値が等しいとき領域の平均画素値に対する重みを
最大とし、残差平方和の総和の最小値が最大値に対して
小さくなるほど領域の平均画素値に対する重みを小さく
する関数であるので、偽構造の発生しにくいフィルタリ
ングを行うことができる。
【００６１】（３０）上記の課題を解決するための他の
観点での発明は、前記重み付け加算して得られた画素値
と前記注目画素の画素値とを重み付け加算する重み係数
は、前記画素グループ態様ごとに求めた前記画素グルー
プの画素値の残差平方和の総和の最小値および前記元画
像のノイズの分散の関数である、ことを特徴とする（２
５）ないし（２９）のうちのいずれか１つに記載の画像
処理装置である。
【００６２】この観点での発明では、（２５）ないし
（２９）のうちのいずれか１つに記載の画像処理におい
て、重み付け加算して得られた画素値と元画像の注目画
素の画素値とを重み付け加算する重み係数を、画素グル
ープ態様ごとに求めた画素グループの画素値の残差平方
和の総和の最小値および元画像のノイズの分散の関数と
したので、元画像の特徴的な構造を損なわないフィルタ
リングを行うことができる。
【００６３】（３１）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値を前記画素グループ態様における画素数で割った
ものと前記ノイズの分散が等しいとき前記重み付け加算
して得られた画素値に対する重みを最大とし、前記残差
平方和の総和の最小値を前記画素グループ態様における
画素数で割ったものと前記ノイズの分散の違いが大きく
なるほど前記重み付け加算して得られた画素値に対する
重みを小さくする関数である、ことを特徴とする（３
０）に記載の画像処理装置である。
【００６４】この観点での発明では、（３０）に記載の
画像処理において、関数は、残差平方和の総和の最小値
を画素グループ態様における画素数で割ったものとノイ
ズの分散が等しいとき重み付け加算して得られた画素値
に対する重みを最大とし、残差平方和の総和の最小値を
画素グループ態様における画素数で割ったものとノイズ
の分散の違いが大きくなるほど重み付け加算して得られ
た画素値に対する重みを小さくする関数であるので、元
画像の特徴的な構造を損なわないフィルタリングを行う
ことができる。
【００６５】（３２）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が最小となる
画素グループ態様を選択することによって行う、ことを
特徴とする（２３）ないし（３１）のうちのいずれか１
つに記載の画像処理装置である。
【００６６】この観点での発明では、（２３）ないし
（３１）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が最小となる画素グループ態様を選択
することによって行うので、元画像の構造に最も良く適
合する画素グループ態様を選択することができる。
【００６７】（３３）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記元画像
のノイズの分散に画素グループ態様における画素数を掛
けたものに最も近くなる画素グループ態様を選択するこ
とによって行う、ことを特徴とする（２３）ないし（３
１）のうちのいずれか１つに記載の画像処理装置であ
る。
【００６８】この観点での発明では、（２３）ないし
（３１）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が元画像のノイズの分散に画素グルー
プ態様における画素数を掛けたものに最も近くなる画素
グループ態様を選択することによって行うので、元画像
の構造に最も良く適合する画素グループ態様を選択する
ことができる。
【００６９】（３４）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定する領域設定手段と、前記設
定した複数態様の領域において、それぞれ、前記注目画
素を含む複数の画素からなる画素グループおよびこの画
素グループの画素と重複しない複数の画素からなる画素
グループを複数の態様で設定する画素グループ設定手段
と、前記設定した複数の画素グループ態様のうち前記設
定した複数態様の領域を通じてその領域における前記元
画像の構造に最も良く適合する画素グループ態様を選択
する選択手段と、前記選択した画素グループ態様におけ
る前記注目画素を含む画素グループの平均画素値を前記
注目画素の新たな画素値として画像を形成する画像形成
手段と、を具備することを特徴とする画像処理装置であ
る。
【００７０】この観点での発明では、複数態様の領域を
通じてその領域における元画像の構造に最も良く適合す
る画素グループ態様を選択し、その画素グループ態様に
おいて、注目画素を含む画素グループの平均画素値を注
目画素の新たな画素値とするので、偽構造が発生しにく
いフィルタリングを行うことができる。
【００７１】（３５）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定する領域設定手段と、前記設
定した複数態様の領域において、それぞれ、前記注目画
素を含む複数の画素からなる画素グループおよびこの画
素グループの画素と重複しない複数の画素からなる画素
グループを複数の態様で設定する画素グループ設定手段
と、前記設定した複数の画素グループ態様のうち前記設
定した複数態様の領域を通じてその領域における前記元
画像の構造に最も良く適合する画素グループ態様を選択
する選択手段と、前記選択した画素グループ態様におけ
る前記注目画素を含む画素グループの平均画素値を求め
る第１の画素値計算手段と、前記複数態様の領域のうち
前記選択した画素グループ態様が設定されている領域の
平均画素値を求める第２の画素値計算手段と、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算した画素値を前記注目画素の新たな画素値とし
て画像を形成する画像形成手段と、を具備することを特
徴とする画像処理装置である。
【００７２】この観点での発明では、（３４）と同様に
して求めた画素グループの平均画素値に、領域の平均画
素値を重み付け加算して注目画素の新たな画素値とする
ので、さらに偽構造が発生しにくいフィルタリングを行
うことができる。
【００７３】（３６）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定する領域設定手段と、前記設
定した複数態様の領域において、それぞれ、前記注目画
素を含む複数の画素からなる画素グループおよびこの画
素グループの画素と重複しない複数の画素からなる画素
グループを複数の態様で設定する画素グループ設定手段
と、前記設定した複数の画素グループ態様のうち前記設
定した複数態様の領域を通じてその領域における前記元
画像の構造に最も良く適合する画素グループ態様を選択
する選択手段と、前記選択した画素グループ態様におけ
る前記注目画素を含む画素グループの平均画素値を求め
る第１の画素値計算手段と、前記複数態様の領域のうち
前記選択した画素グループ態様が設定されている領域の
平均画素値を求める第２の画素値計算手段と、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算する加算手段と、前記重み付け加算して得られ
た画素値と前記注目画素の画素値とを重み付け加算した
画素値を前記注目画素の新たな画素値として画像を形成
する画像形成手段と、を具備することを特徴とする画像
処理装置である。
【００７４】この観点での発明では、（３５）と同様に
して求めた画素値に元画像の注目画素の画素値を重み付
け加算して注目画素の新たな画素値とするので、元画像
の特徴的な構造を損なわないフィルタリングを行うこと
ができる。
【００７５】（３７）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定する領域設定手段と、前記設
定した複数態様の領域において、それぞれ、前記注目画
素を含む複数の画素からなる画素グループおよびこの画
素グループの画素と重複しない複数の画素からなる画素
グループを複数の態様で設定する画素グループ設定手段
と、前記設定した複数の画素グループ態様のうち前記設
定した複数態様の領域を通じてその領域における前記元
画像の構造に最も良く適合する画素グループ態様を選択
する選択手段と、前記選択した画素グループ態様におけ
る前記注目画素を含む画素グループの平均画素値を求め
る第１の画素値計算手段と、前記複数態様の領域のうち
前記選択した画素グループ態様が設定されている領域の
平均画素値を求める第２の画素値計算手段と、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算する第１の加算手段と、前記重み付け加算して
得られた画素値と前記注目画素の画素値とを重み付け加
算した画素値を前記注目画素の新たな画素値として画像
を形成する画像形成手段と、前記形成した画像と前記元
画像を重み付け加算する第２の加算手段と、を具備する
ことを特徴とする画像処理装置である。
【００７６】この観点での発明では、（３６）と同様に
して求めた画素値によって形成した画像に元画像を重み
付け加算するので、フィルタリングの効果を弱めること
ができる。
【００７７】（３８）上記の課題を解決するための他の
観点での発明は、前記形成した画像と前記元画像を重み
付け加算する重み係数は調節可能である、ことを特徴と
する（３７）に記載の画像処理装置である。
【００７８】この観点での発明では、重み係数の調整に
よりフィルタリングの効果を調節することができる。
（３９）上記の課題を解決するための他の観点での発明
は、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する重み係数は、前記設定した複
数態様の領域を通じて前記画像グループ態様ごとに求め
た前記画素グループの画素値の残差平方和の総和の最小
値および最大値の関数である、ことを特徴とする（３
５）ないし（３８）のうちのいずれか１つに記載の画像
処理装置である。
【００７９】この観点での発明では、（３５）ないし
（３８）のうちのいずれか１つに記載の画像処理におい
て、画素グループの平均画素値と領域の平均画素値とを
重み付け加算する重み係数を、複数態様の領域を通じて
画素グループ態様ごとに求めた画素グループの画素値の
残差平方和の総和の最小値および最大値の関数としたの
で、偽構造の発生しにくいフィルタリングを行うことが
できる。
【００８０】（４０）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値と最大値が等しいとき前記領域の平均画素値に対
する重みを最大とし、前記残差平方和の総和の最小値が
最大値に対して小さくなるほど前記領域の平均画素値に
対する重みを小さくする関数である、ことを特徴とする
（３９）に記載の画像処理装置である。
【００８１】この観点での発明では、（３９）に記載の
画像処理において、関数は、残差平方和の総和の最小値
と最大値が等しいとき領域の平均画素値に対する重みを
最大とし、残差平方和の総和の最小値が最大値に対して
小さくなるほど領域の平均画素値に対する重みを小さく
する関数であるので、偽構造の発生しにくいフィルタリ
ングを行うことができる。
【００８２】（４１）上記の課題を解決するための他の
観点での発明は、前記重み付け加算して得られた画素値
と前記注目画素の画素値とを重み付け加算する重み係数
は、前記設定した複数態様の領域を通じて前記画像グル
ープ態様ごとに求めた前記画素グループの画素値の残差
平方和の総和の最小値および前記元画像のノイズの分散
の関数である、ことを特徴とする（３６）ないし（４
０）のうちのいずれか１つに記載の画像処理装置であ
る。
【００８３】この観点での発明では、（３６）ないし
（４０）のうちのいずれか１つに記載の画像処理におい
て、重み付け加算して得られた画素値と元画像の注目画
素の画素値とを重み付け加算する重み係数を、複数態様
の領域を通じて画素グループ態様ごとに求めた画素グル
ープの画素値の残差平方和の総和の最小値および元画像
のノイズの分散の関数としたので、元画像の特徴的な構
造を損なわないフィルタリングを行うことができる。
【００８４】（４２）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値を前記画素グループ態様における画素数で割った
ものと前記ノイズの分散が等しいとき前記重み付け加算
して得られた画素値に対する重みを最大とし、前記残差
平方和の総和の最小値を前記画素グループ態様における
画素数で割ったものと前記ノイズの分散の違いが大きく
なるほど前記重み付け加算して得られた画素値に対する
重みを小さくする関数である、ことを特徴とする（４
１）に記載の画像処理装置である。
【００８５】この観点での発明では、（４１）に記載の
画像処理において、関数は、残差平方和の総和の最小値
を画素グループ態様における画素数で割ったものとノイ
ズの分散が等しいとき重み付け加算して得られた画素値
に対する重みを最大とし、残差平方和の総和の最小値を
画素グループ態様における画素数で割ったものとノイズ
の分散の違いが大きくなるほど重み付け加算して得られ
た画素値に対する重みを小さくする関数であるので、元
画像の特徴的な構造を損なわないフィルタリングを行う
ことができる。
【００８６】（４３）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記設定し
た複数態様の領域を通じて最小となる画像グループ態様
を選択することによって行う、ことを特徴とする（３
４）ないし（４２）のうちのいずれか１つに記載の画像
処理装置である。
【００８７】この観点での発明では、（３４）ないし
（４２）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が複数態様の領域を通じて最小となる
画素グループ態様を選択することによって行うので、元
画像の構造に最も良く適合する画素グループ態様を選択
することができる。
【００８８】（４４）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記元画像
のノイズの分散に画素グループ態様における画素数を掛
けたものに最も近くなる画素グループ態様を選択するこ
とによって行う、ことを特徴とする（３４）ないし（４
２）のうちのいずれか１つに記載の画像処理装置であ
る。
【００８９】この観点での発明では、（３４）ないし
（４２）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が元画像のノイズの分散に画素グルー
プ態様における画素数を掛けたものに最も近くなる画素
グループ態様を選択することによって行うので、元画像
の構造に最も良く適合する画素グループ態様を選択する
ことができる。
【００９０】（４５）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を設定し、前記領域において前記注目画素を含む
複数の画素からなる画素グループおよびこの画素グルー
プの画素と重複しない複数の画素からなる画素グループ
を複数の態様で設定し、前記設定した複数の画素グルー
プ態様のうち前記領域における前記元画像の構造に最も
良く適合する画素グループ態様を選択し、前記選択した
画素グループ態様における前記注目画素を含む画素グル
ープの平均画素値を前記注目画素の新たな画素値として
画像を形成する、機能をコンピュータに実現させるプロ
グラムをコンピュータにより読み取り可能なように記録
したことを特徴とする記録媒体である。
【００９１】この観点での発明では、記録媒体に記録さ
れたプログラムが、元画像の構造に最も良く適合する画
素グループ態様を選択し、その画素グループ態様におい
て、注目画素を含む画素グループの平均画素値を注目画
素の新たな画素値とする、機能をコンピュータに実現さ
せるので、偽構造が発生しにくいフィルタリングを行う
ことができる。
【００９２】（４６）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を設定し、前記領域において前記注目画素を含む
複数の画素からなる画素グループおよびこの画素グルー
プの画素と重複しない複数の画素からなる画素グループ
を複数の態様で設定し、前記設定した複数の画素グルー
プ態様のうち前記領域における前記元画像の構造に最も
良く適合する画素グループ態様を選択し、前記選択した
画素グループ態様における前記注目画素を含む画素グル
ープの平均画素値を求め、前記領域の平均画素値を求
め、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算した画素値を前記注目画素の新た
な画素値として画像を形成する、機能をコンピュータに
実現させるプログラムをコンピュータにより読み取り可
能なように記録したことを特徴とする記録媒体である。
【００９３】この観点での発明では、記録媒体に記録さ
れたプログラムが、（４５）と同様にして求めた画素グ
ループの平均画素値に、領域の平均画素値を重み付け加
算して注目画素の新たな画素値とする、機能をコンピュ
ータに実現させるので、さらに偽構造が発生しにくいフ
ィルタリングを行うことができる。
【００９４】（４７）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を設定し、前記領域において前記注目画素を含む
複数の画素からなる画素グループおよびこの画素グルー
プの画素と重複しない複数の画素からなる画素グループ
を複数の態様で設定し、前記設定した複数の画素グルー
プ態様のうち前記領域における前記元画像の構造に最も
良く適合する画素グループ態様を選択し、前記選択した
画素グループ態様における前記注目画素を含む画素グル
ープの平均画素値を求め、前記領域の平均画素値を求
め、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算し、前記重み付け加算して得られ
た画素値と前記注目画素の画素値とを重み付け加算した
画素値を前記注目画素の新たな画素値として画像を形成
する、機能をコンピュータに実現させるプログラムをコ
ンピュータにより読み取り可能なように記録したことを
特徴とする記録媒体である。
【００９５】この観点での発明では、記録媒体に記録さ
れたプログラムが、（４６）と同様にして求めた画素値
に元画像の注目画素の画素値を重み付け加算して注目画
素の新たな画素値とする、機能をコンピュータに実現さ
せるので、元画像の特徴的な構造を損なわないフィルタ
リングを行うことができる。
【００９６】（４８）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を設定し、前記領域において前記注目画素を含む
複数の画素からなる画素グループおよびこの画素グルー
プの画素と重複しない複数の画素からなる画素グループ
を複数の態様で設定し、前記設定した複数の画素グルー
プ態様のうち前記領域における前記元画像の構造に最も
良く適合する画素グループ態様を選択し、前記選択した
画素グループ態様における前記注目画素を含む画素グル
ープの平均画素値を求め、前記領域の平均画素値を求
め、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算し、前記重み付け加算して得られ
た画素値と前記注目画素の画素値とを重み付け加算した
画素値を前記注目画素の新たな画素値として画像を形成
し、前記形成した画像と前記元画像を重み付け加算す
る、機能をコンピュータに実現させるプログラムをコン
ピュータにより読み取り可能なように記録したことを特
徴とする記録媒体である。
【００９７】この観点での発明では、記録媒体に記録さ
れたプログラムが、（４７）と同様にして求めた画素値
によって形成した画像に元画像を重み付け加算する、機
能をコンピュータに実現させるので、フィルタリングの
効果を弱めることができる。
【００９８】（４９）上記の課題を解決するための他の
観点での発明は、前記形成した画像と前記元画像を重み
付け加算する重み係数は調節可能である、ことを特徴と
する（４８）に記載の記録媒体である。
【００９９】この観点での発明では、重み係数の調整に
よりフィルタリングの効果を調節することができる。
（５０）上記の課題を解決するための他の観点での発明
は、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する重み係数は、前記画素グルー
プ態様ごとに求めた前記画素グループの画素値の残差平
方和の総和の最小値および最大値の関数である、ことを
特徴とする（４６）ないし（４９）のうちのいずれか１
つに記載の記録媒体である。
【０１００】この観点での発明では、（４６）ないし
（４９）のうちのいずれか１つに記載の画像処理におい
て、画素グループの平均画素値と領域の平均画素値とを
重み付け加算する重み係数を、画素グループ態様ごとに
求めた画素グループの画素値の残差平方和の総和の最小
値および最大値の関数としたので、偽構造の発生しにく
いフィルタリングを行うことができる。
【０１０１】（５１）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値と最大値が等しいとき前記領域の平均画素値に対
する重みを最大とし、前記残差平方和の総和の最小値が
最大値に対して小さくなるほど前記領域の平均画素値に
対する重みを小さくする関数である、ことを特徴とする
（５０）に記載の記録媒体である。
【０１０２】この観点での発明では、（５０）に記載の
画像処理において、関数は、残差平方和の総和の最小値
と最大値が等しいとき領域の平均画素値に対する重みを
最大とし、残差平方和の総和の最小値が最大値に対して
小さくなるほど領域の平均画素値に対する重みを小さく
する関数であるので、偽構造の発生しにくいフィルタリ
ングを行うことができる。
【０１０３】（５２）上記の課題を解決するための他の
観点での発明は、前記重み付け加算して得られた画素値
と前記注目画素の画素値とを重み付け加算する重み係数
は、前記画素グループ態様ごとに求めた前記画素グルー
プの画素値の残差平方和の総和の最小値および前記元画
像のノイズの分散の関数である、ことを特徴とする（４
７）ないし（５１）のうちのいずれか１つに記載の記録
媒体である。
【０１０４】この観点での発明では、（４７）ないし
（５１）のうちのいずれか１つに記載の画像処理におい
て、重み付け加算して得られた画素値と元画像の注目画
素の画素値とを重み付け加算する重み係数を、画素グル
ープ態様ごとに求めた画素グループの画素値の残差平方
和の総和の最小値および元画像のノイズの分散の関数と
したので、元画像の特徴的な構造を損なわないフィルタ
リングを行うことができる。
【０１０５】（５３）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値を前記画素グループ態様における画素数で割った
ものと前記ノイズの分散が等しいとき前記重み付け加算
して得られた画素値に対する重みを最大とし、前記残差
平方和の総和の最小値を前記画素グループ態様における
画素数で割ったものと前記ノイズの分散の違いが大きく
なるほど前記重み付け加算して得られた画素値に対する
重みを小さくする関数である、ことを特徴とする（５
２）に記載の記録媒体である。
【０１０６】この観点での発明では、（５２）に記載の
画像処理において、関数は、残差平方和の総和の最小値
を画素グループ態様における画素数で割ったものとノイ
ズの分散が等しいとき重み付け加算して得られた画素値
に対する重みを最大とし、残差平方和の総和の最小値を
画素グループ態様における画素数で割ったものとノイズ
の分散の違いが大きくなるほど重み付け加算して得られ
た画素値に対する重みを小さくする関数であるので、元
画像の特徴的な構造を損なわないフィルタリングを行う
ことができる。
【０１０７】（５４）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が最小となる
画素グループ態様を選択することによって行う、ことを
特徴とする（４５）ないし（５３）のうちのいずれか１
つに記載の記録媒体である。
【０１０８】この観点での発明では、（４５）ないし
（５３）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が最小となる画素グループ態様を選択
することによって行うので、元画像の構造に最も良く適
合する画素グループ態様を選択することができる。
【０１０９】（５５）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記元画像
のノイズの分散に画素グループ態様における画素数を掛
けたものに最も近くなる画素グループ態様を選択するこ
とによって行う、ことを特徴とする（４５）ないし（５
３）のうちのいずれか１つに記載の記録媒体である。
【０１１０】この観点での発明では、（４５）ないし
（５３）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が元画像のノイズの分散に画素グルー
プ態様における画素数を掛けたものに最も近くなる画素
グループ態様を選択することによって行うので、元画像
の構造に最も良く適合する画素グループ態様を選択する
ことができる。
【０１１１】（５６）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定し、前記設定した複数態様の
領域において、それぞれ、前記注目画素を含む複数の画
素からなる画素グループおよびこの画素グループの画素
と重複しない複数の画素からなる画素グループを複数の
態様で設定し、前記設定した複数の画素グループ態様の
うち前記設定した複数態様の領域を通じてその領域にお
ける前記元画像の構造に最も良く適合する画素グループ
態様を選択し、前記選択した画素グループ態様における
前記注目画素を含む画素グループの平均画素値を前記注
目画素の新たな画素値として画像を形成する、機能をコ
ンピュータに実現させるプログラムをコンピュータによ
り読み取り可能なように記録したことを特徴とする記録
媒体である。
【０１１２】この観点での発明では、記録媒体に記録さ
れたプログラムが、複数態様の領域を通じてその領域に
おける元画像の構造に最も良く適合する画素グループ態
様を選択し、その画素グループ態様において、注目画素
を含む画素グループの平均画素値を注目画素の新たな画
素値とする、機能をコンピュータに実現させるので、偽
構造が発生しにくいフィルタリングを行うことができ
る。
【０１１３】（５７）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定し、前記設定した複数態様の
領域において、それぞれ、前記注目画素を含む複数の画
素からなる画素グループおよびこの画素グループの画素
と重複しない複数の画素からなる画素グループを複数の
態様で設定し、前記設定した複数の画素グループ態様の
うち前記設定した複数態様の領域を通じてその領域にお
ける前記元画像の構造に最も良く適合する画素グループ
態様を選択し、前記選択した画素グループ態様における
前記注目画素を含む画素グループの平均画素値を求め、
前記複数態様の領域のうち前記選択した画素グループ態
様が設定されている領域の平均画素値を求め、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算した画素値を前記注目画素の新たな画素値とし
て画像を形成する、機能をコンピュータに実現させるプ
ログラムをコンピュータにより読み取り可能なように記
録したことを特徴とする記録媒体である。
【０１１４】この観点での発明では、記録媒体に記録さ
れたプログラムが、（５６）と同様にして求めた画素グ
ループの平均画素値に、領域の平均画素値を重み付け加
算して注目画素の新たな画素値とする、機能をコンピュ
ータに実現させるので、さらに偽構造が発生しにくいフ
ィルタリングを行うことができる。
【０１１５】（５８）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定し、前記設定した複数態様の
領域において、それぞれ、前記注目画素を含む複数の画
素からなる画素グループおよびこの画素グループの画素
と重複しない複数の画素からなる画素グループを複数の
態様で設定し、前記設定した複数の画素グループ態様の
うち前記設定した複数態様の領域を通じてその領域にお
ける前記元画像の構造に最も良く適合する画素グループ
態様を選択し、前記選択した画素グループ態様における
前記注目画素を含む画素グループの平均画素値を求め、
前記複数態様の領域のうち前記選択した画素グループ態
様が設定されている領域の平均画素値を求め、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算し、前記重み付け加算して得られた画素値と前
記注目画素の画素値とを重み付け加算した画素値を前記
注目画素の新たな画素値として画像を形成する、機能を
コンピュータに実現させるプログラムをコンピュータに
より読み取り可能なように記録したことを特徴とする記
録媒体である。
【０１１６】この観点での発明では、記録媒体に記録さ
れたプログラムが、（５７）と同様にして求めた画素値
に元画像の注目画素の画素値を重み付け加算して注目画
素の新たな画素値とする、機能をコンピュータに実現さ
せるので、元画像の特徴的な構造を損なわないフィルタ
リングを行うことができる。
【０１１７】（５９）上記の課題を解決するための他の
観点での発明は、元画像において注目画素を含む局所的
な領域を複数の態様で設定し、前記設定した複数態様の
領域において、それぞれ、前記注目画素を含む複数の画
素からなる画素グループおよびこの画素グループの画素
と重複しない複数の画素からなる画素グループを複数の
態様で設定し、前記設定した複数の画素グループ態様の
うち前記設定した複数態様の領域を通じてその領域にお
ける前記元画像の構造に最も良く適合する画素グループ
態様を選択し、前記選択した画素グループ態様における
前記注目画素を含む画素グループの平均画素値を求め、
前記複数態様の領域のうち前記選択した画素グループ態
様が設定されている領域の平均画素値を求め、前記画素
グループの平均画素値と前記領域の平均画素値とを重み
付け加算し、前記重み付け加算して得られた画素値と前
記注目画素の画素値とを重み付け加算した画素値を前記
注目画素の新たな画素値として画像を形成し、前記形成
した画像と前記元画像を重み付け加算する、機能をコン
ピュータに実現させるプログラムをコンピュータにより
読み取り可能なように記録したことを特徴とする記録媒
体である。
【０１１８】この観点での発明では、記録媒体に記録さ
れたプログラムが、（５８）と同様にして求めた画素値
によって形成した画像に元画像を重み付け加算する、機
能をコンピュータに実現させるので、フィルタリングの
効果を弱めることができる。
【０１１９】（６０）上記の課題を解決するための他の
観点での発明は、前記形成した画像と前記元画像を重み
付け加算する重み係数は調節可能である、ことを特徴と
する（５９）に記載の記録媒体である。
【０１２０】この観点での発明では、重み係数の調整に
よりフィルタリングの効果を調節することができる。
（６１）上記の課題を解決するための他の観点での発明
は、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する重み係数は、前記設定した複
数態様の領域を通じて前記画像グループ態様ごとに求め
た前記画素グループの画素値の残差平方和の総和の最小
値および最大値の関数である、ことを特徴とする（５
７）ないし（６０）のうちのいずれか１つに記載の記録
媒体である。
【０１２１】この観点での発明では、（５７）ないし
（６０）のうちのいずれか１つに記載の画像処理におい
て、画素グループの平均画素値と領域の平均画素値とを
重み付け加算する重み係数を、複数態様の領域を通じて
画素グループ態様ごとに求めた画素グループの画素値の
残差平方和の総和の最小値および最大値の関数としたの
で、偽構造の発生しにくいフィルタリングを行うことが
できる。
【０１２２】（６２）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値と最大値が等しいとき前記領域の平均画素値に対
する重みを最大とし、前記残差平方和の総和の最小値が
最大値に対して小さくなるほど前記領域の平均画素値に
対する重みを小さくする関数である、ことを特徴とする
（６１）に記載の記録媒体である。
【０１２３】この観点での発明では、（６１）に記載の
画像処理において、関数は、残差平方和の総和の最小値
と最大値が等しいとき領域の平均画素値に対する重みを
最大とし、残差平方和の総和の最小値が最大値に対して
小さくなるほど領域の平均画素値に対する重みを小さく
する関数であるので、偽構造の発生しにくいフィルタリ
ングを行うことができる。
【０１２４】（６３）上記の課題を解決するための他の
観点での発明は、前記重み付け加算して得られた画素値
と前記注目画素の画素値とを重み付け加算する重み係数
は、前記設定した複数態様の領域を通じて前記画像グル
ープ態様ごとに求めた前記画素グループの画素値の残差
平方和の総和の最小値および前記元画像のノイズの分散
の関数である、ことを特徴とする（５８）ないし（６
２）のうちのいずれか１つに記載の記録媒体である。
【０１２５】この観点での発明では、（５８）ないし
（６２）のうちのいずれか１つに記載の画像処理におい
て、重み付け加算して得られた画素値と元画像の注目画
素の画素値とを重み付け加算する重み係数を、複数態様
の領域を通じて画素グループ態様ごとに求めた画素グル
ープの画素値の残差平方和の総和の最小値および元画像
のノイズの分散の関数としたので、元画像の特徴的な構
造を損なわないフィルタリングを行うことができる。
【０１２６】（６４）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値を前記画素グループ態様における画素数で割った
ものと前記ノイズの分散が等しいとき前記重み付け加算
して得られた画素値に対する重みを最大とし、前記残差
平方和の総和の最小値を前記画素グループ態様における
画素数で割ったものと前記ノイズの分散の違いが大きく
なるほど前記重み付け加算して得られた画素値に対する
重みを小さくする関数である、ことを特徴とする（６
３）に記載の記録媒体である。
【０１２７】この観点での発明では、（６３）に記載の
画像処理において、関数は、残差平方和の総和の最小値
を画素グループ態様における画素数で割ったものとノイ
ズの分散が等しいとき重み付け加算して得られた画素値
に対する重みを最大とし、残差平方和の総和の最小値を
画素グループ態様における画素数で割ったものとノイズ
の分散の違いが大きくなるほど重み付け加算して得られ
た画素値に対する重みを小さくする関数であるので、元
画像の特徴的な構造を損なわないフィルタリングを行う
ことができる。
【０１２８】（６５）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記設定し
た複数態様の領域を通じて最小となる画像グループ態様
を選択することによって行う、ことを特徴とする（５
６）ないし（６４）のうちのいずれか１つに記載の記録
媒体である。
【０１２９】この観点での発明では、（５６）ないし
（６４）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が複数態様の領域を通じて最小となる
画素グループ態様を選択することによって行うので、元
画像の構造に最も良く適合する画素グループ態様を選択
することができる。
【０１３０】（６６）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記元画像
のノイズの分散に画素グループ態様における画素数を掛
けたものに最も近くなる画素グループ態様を選択するこ
とによって行う、ことを特徴とする（５６）ないし（６
４）のうちのいずれか１つに記載の記録媒体である。
【０１３１】この観点での発明では、（５６）ないし
（６４）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が元画像のノイズの分散に画素グルー
プ態様における画素数を掛けたものに最も近くなる画素
グループ態様を選択することによって行うので、元画像
の構造に最も良く適合する画素グループ態様を選択する
ことができる。
【０１３２】（６７）上記の課題を解決するための他の
観点での発明は、対象から信号を収集する信号収集手段
と、前記収集した信号に基づいて元画像を生成する元画
像生成手段と、前記元画像において注目画素を含む局所
的な領域を設定する領域設定手段と、前記領域において
前記注目画素を含む複数の画素からなる画素グループお
よびこの画素グループの画素と重複しない複数の画素か
らなる画素グループを複数の態様で設定する画素グルー
プ設定手段と、前記設定した複数の画素グループ態様の
うち前記領域における前記元画像の構造に最も良く適合
する画素グループ態様を選択する選択手段と、前記選択
した画素グループ態様における前記注目画素を含む画素
グループの平均画素値を前記注目画素の新たな画素値と
して画像を形成する画像形成手段と、を具備することを
特徴とする画像撮影装置である。
【０１３３】この観点での発明では、元画像の構造に最
も良く適合する画素グループ態様を選択し、その画素グ
ループ態様において、注目画素を含む画素グループの平
均画素値を注目画素の新たな画素値とするので、偽構造
が発生しにくいフィルタリングを施した画像を得ること
ができる。
【０１３４】（６８）上記の課題を解決するための他の
観点での発明は、対象から信号を収集する信号収集手段
と、前記収集した信号に基づいて元画像を生成する元画
像生成手段と、前記元画像において注目画素を含む局所
的な領域を設定する領域設定手段と、前記領域において
前記注目画素を含む複数の画素からなる画素グループお
よびこの画素グループの画素と重複しない複数の画素か
らなる画素グループを複数の態様で設定する画素グルー
プ設定手段と、前記設定した複数の画素グループ態様の
うち前記領域における前記元画像の構造に最も良く適合
する画素グループ態様を選択する選択手段と、前記選択
した画素グループ態様における前記注目画素を含む画素
グループの平均画素値を求める第１の画素値計算手段
と、前記領域の平均画素値を求める第２の画素値計算手
段と、前記画素グループの平均画素値と前記領域の平均
画素値とを重み付け加算した画素値を前記注目画素の新
たな画素値として画像を形成する画像形成手段と、を具
備することを特徴とする画像撮影装置である。
【０１３５】この観点での発明では、（６７）と同様に
して求めた画素グループの平均画素値に、領域の平均画
素値を重み付け加算して注目画素の新たな画素値とする
ので、さらに偽構造が発生しにくいフィルタリングを施
した画像を得ることができる。
【０１３６】（６９）上記の課題を解決するための他の
観点での発明は、対象から信号を収集する信号収集手段
と、前記収集した信号に基づいて元画像を生成する元画
像生成手段と、前記元画像において注目画素を含む局所
的な領域を設定する領域設定手段と、前記領域において
前記注目画素を含む複数の画素からなる画素グループお
よびこの画素グループの画素と重複しない複数の画素か
らなる画素グループを複数の態様で設定する画素グルー
プ設定手段と、前記設定した複数の画素グループ態様の
うち前記領域における前記元画像の構造に最も良く適合
する画素グループ態様を選択する選択手段と、前記選択
した画素グループ態様における前記注目画素を含む画素
グループの平均画素値を求める第１の画素値計算手段
と、前記領域の平均画素値を求める第２の画素値計算手
段と、前記画素グループの平均画素値と前記領域の平均
画素値とを重み付け加算する加算手段と、前記重み付け
加算して得られた画素値と前記注目画素の画素値とを重
み付け加算した画素値を前記注目画素の新たな画素値と
して画像を形成する画像形成手段と、を具備することを
特徴とする画像撮影装置である。
【０１３７】この観点での発明では、（６８）と同様に
して求めた画素値に元画像の注目画素の画素値を重み付
け加算して注目画素の新たな画素値とするので、元画像
の特徴的な構造を損なわないフィルタリングを施した画
像を得ることができる。
【０１３８】（７０）上記の課題を解決するための他の
観点での発明は、対象から信号を収集する信号収集手段
と、前記収集した信号に基づいて元画像を生成する元画
像生成手段と、前記元画像において注目画素を含む局所
的な領域を設定する領域設定手段と、前記領域において
前記注目画素を含む複数の画素からなる画素グループお
よびこの画素グループの画素と重複しない複数の画素か
らなる画素グループを複数の態様で設定する画素グルー
プ設定手段と、前記設定した複数の画素グループ態様の
うち前記領域における前記元画像の構造に最も良く適合
する画素グループ態様を選択する選択手段と、前記選択
した画素グループ態様における前記注目画素を含む画素
グループの平均画素値を求める第１の画素値計算手段
と、前記領域の平均画素値を求める第２の画素値計算手
段と、前記画素グループの平均画素値と前記領域の平均
画素値とを重み付け加算する第１の加算手段と、前記重
み付け加算して得られた画素値と前記注目画素の画素値
とを重み付け加算した画素値を前記注目画素の新たな画
素値として画像を形成する画像形成手段と、前記形成し
た画像と前記元画像を重み付け加算する第２の加算手段
と、を具備することを特徴とする画像撮影装置である。
【０１３９】この観点での発明では、（６９）と同様に
して求めた画素値によって形成した画像に元画像を重み
付け加算するので、フィルタリングの効果を弱めること
ができる。
【０１４０】（７１）上記の課題を解決するための他の
観点での発明は、前記形成した画像と前記元画像を重み
付け加算する重み係数は調節可能である、ことを特徴と
する（７０）に記載の画像撮影装置である。
【０１４１】この観点での発明では、重み係数の調整に
よりフィルタリングの効果を調節することができる。
（７２）上記の課題を解決するための他の観点での発明
は、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する重み係数は、前記画素グルー
プ態様ごとに求めた前記画素グループの画素値の残差平
方和の総和の最小値および最大値の関数である、ことを
特徴とする（６８）ないし（７１）のうちのいずれか１
つに記載の画像撮影装置である。
【０１４２】この観点での発明では、（６８）ないし
（７１）のうちのいずれか１つに記載の画像処理におい
て、画素グループの平均画素値と領域の平均画素値とを
重み付け加算する重み係数を、画素グループ態様ごとに
求めた画素グループの画素値の残差平方和の総和の最小
値および最大値の関数としたので、偽構造の発生しにく
いフィルタリングを施した画像を得ることができる。
【０１４３】（７３）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値と最大値が等しいとき前記領域の平均画素値に対
する重みを最大とし、前記残差平方和の総和の最小値が
最大値に対して小さくなるほど前記領域の平均画素値に
対する重みを小さくする関数である、ことを特徴とする
（７２）に記載の画像撮影装置である。
【０１４４】この観点での発明では、（７２）に記載の
画像処理において、関数は、残差平方和の総和の最小値
と最大値が等しいとき領域の平均画素値に対する重みを
最大とし、残差平方和の総和の最小値が最大値に対して
小さくなるほど領域の平均画素値に対する重みを小さく
する関数であるので、偽構造の発生しにくいフィルタリ
ングを行うことができる。
【０１４５】（７４）上記の課題を解決するための他の
観点での発明は、前記重み付け加算して得られた画素値
と前記注目画素の画素値とを重み付け加算する重み係数
は、前記画素グループ態様ごとに求めた前記画素グルー
プの画素値の残差平方和の総和の最小値および前記元画
像のノイズの分散の関数である、ことを特徴とする（６
９）ないし（７３）のうちのいずれか１つに記載の画像
撮影装置である。
【０１４６】この観点での発明では、（６９）ないし
（７３）のうちのいずれか１つに記載の画像処理におい
て、重み付け加算して得られた画素値と元画像の注目画
素の画素値とを重み付け加算する重み係数を、画素グル
ープ態様ごとに求めた画素グループの画素値の残差平方
和の総和の最小値および元画像のノイズの分散の関数と
したので、元画像の特徴的な構造を損なわないフィルタ
リングを施した画像を得ることができる。
【０１４７】（７５）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値を前記画素グループ態様における画素数で割った
ものと前記ノイズの分散が等しいとき前記重み付け加算
して得られた画素値に対する重みを最大とし、前記残差
平方和の総和の最小値を前記画素グループ態様における
画素数で割ったものと前記ノイズの分散の違いが大きく
なるほど前記重み付け加算して得られた画素値に対する
重みを小さくする関数である、ことを特徴とする（７
４）に記載の画像撮影装置である。
【０１４８】この観点での発明では、（７４）に記載の
画像処理において、関数は、残差平方和の総和の最小値
を画素グループ態様における画素数で割ったものとノイ
ズの分散が等しいとき重み付け加算して得られた画素値
に対する重みを最大とし、残差平方和の総和の最小値を
画素グループ態様における画素数で割ったものとノイズ
の分散の違いが大きくなるほど重み付け加算して得られ
た画素値に対する重みを小さくする関数であるので、元
画像の特徴的な構造を損なわないフィルタリングを行う
ことができる。
【０１４９】（７６）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が最小となる
画素グループ態様を選択することによって行う、ことを
特徴とする（６７）ないし（７５）のうちのいずれか１
つに記載の画像撮影装置である。
【０１５０】この観点での発明では、（６７）ないし
（７５）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が最小となる画素グループ態様を選択
することによって行うので、元画像の構造に最も良く適
合する画素グループ態様を選択することができる。
【０１５１】（７７）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記元画像
のノイズの分散に画素グループ態様における画素数を掛
けたものに最も近くなる画素グループ態様を選択するこ
とによって行う、ことを特徴とする（６７）ないし（７
５）のうちのいずれか１つに記載の画像撮影装置であ
る。
【０１５２】この観点での発明では、（６７）ないし
（７５）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が元画像のノイズの分散に画素グルー
プ態様における画素数を掛けたものに最も近くなる画素
グループ態様を選択することによって行うので、元画像
の構造に最も良く適合する画素グループ態様を選択する
ことができる。
【０１５３】（７８）上記の課題を解決するための他の
観点での発明は、対象から信号を収集する信号収集手段
と、前記収集した信号に基づいて元画像を生成する元画
像生成手段と、前記元画像において注目画素を含む局所
的な領域を複数の態様で設定する領域設定手段と、前記
設定した複数態様の領域において、それぞれ、前記注目
画素を含む複数の画素からなる画素グループおよびこの
画素グループの画素と重複しない複数の画素からなる画
素グループを複数の態様で設定する画素グループ設定手
段と、前記設定した複数の画素グループ態様のうち前記
設定した複数態様の領域を通じてその領域における前記
元画像の構造に最も良く適合する画素グループ態様を選
択する選択手段と、前記選択した画素グループ態様にお
ける前記注目画素を含む画素グループの平均画素値を前
記注目画素の新たな画素値として画像を形成する画像形
成手段と、を具備することを特徴とする画像撮影装置で
ある。
【０１５４】この観点での発明では、複数態様の領域を
通じてその領域における元画像の構造に最も良く適合す
る画素グループ態様を選択し、その画素グループ態様に
おいて、注目画素を含む画素グループの平均画素値を注
目画素の新たな画素値とするので、偽構造が発生しにく
いフィルタリングを施した画像を得ることができる。
【０１５５】（７９）上記の課題を解決するための他の
観点での発明は、対象から信号を収集する信号収集手段
と、前記収集した信号に基づいて元画像を生成する元画
像生成手段と、前記元画像において注目画素を含む局所
的な領域を複数の態様で設定する領域設定手段と、前記
設定した複数態様の領域において、それぞれ、前記注目
画素を含む複数の画素からなる画素グループおよびこの
画素グループの画素と重複しない複数の画素からなる画
素グループを複数の態様で設定する画素グループ設定手
段と、前記設定した複数の画素グループ態様のうち前記
設定した複数態様の領域を通じてその領域における前記
元画像の構造に最も良く適合する画素グループ態様を選
択する選択手段と、前記選択した画素グループ態様にお
ける前記注目画素を含む画素グループの平均画素値を求
める第１の画素値計算手段と、前記複数態様の領域のう
ち前記選択した画素グループ態様が設定されている領域
の平均画素値を求める第２の画素値計算手段と、前記画
素グループの平均画素値と前記領域の平均画素値とを重
み付け加算した画素値を前記注目画素の新たな画素値と
して画像を形成する画像形成手段と、を具備することを
特徴とする画像撮影装置である。
【０１５６】この観点での発明では、（７８）と同様に
して求めた画素グループの平均画素値に、領域の平均画
素値を重み付け加算して注目画素の新たな画素値とする
ので、さらに偽構造が発生しにくいフィルタリングを施
した画像を得ることができる。
【０１５７】（８０）上記の課題を解決するための他の
観点での発明は、対象から信号を収集する信号収集手段
と、前記収集した信号に基づいて元画像を生成する元画
像生成手段と、前記元画像において注目画素を含む局所
的な領域を複数の態様で設定する領域設定手段と、前記
設定した複数態様の領域において、それぞれ、前記注目
画素を含む複数の画素からなる画素グループおよびこの
画素グループの画素と重複しない複数の画素からなる画
素グループを複数の態様で設定する画素グループ設定手
段と、前記設定した複数の画素グループ態様のうち前記
設定した複数態様の領域を通じてその領域における前記
元画像の構造に最も良く適合する画素グループ態様を選
択する選択手段と、前記選択した画素グループ態様にお
ける前記注目画素を含む画素グループの平均画素値を求
める第１の画素値計算手段と、前記複数態様の領域のう
ち前記選択した画素グループ態様が設定されている領域
の平均画素値を求める第２の画素値計算手段と、前記画
素グループの平均画素値と前記領域の平均画素値とを重
み付け加算する加算手段と、前記重み付け加算して得ら
れた画素値と前記注目画素の画素値とを重み付け加算し
た画素値を前記注目画素の新たな画素値として画像を形
成する画像形成手段と、を具備することを特徴とする画
像撮影装置である。
【０１５８】この観点での発明では、（７９）と同様に
して求めた画素値に元画像の注目画素の画素値を重み付
け加算して注目画素の新たな画素値とするので、元画像
の特徴的な構造を損なわないフィルタリングを施した画
像を得ることができる。
【０１５９】（８１）上記の課題を解決するための他の
観点での発明は、対象から信号を収集する信号収集手段
と、前記収集した信号に基づいて元画像を生成する元画
像生成手段と、前記元画像において注目画素を含む局所
的な領域を複数の態様で設定する領域設定手段と、前記
設定した複数態様の領域において、それぞれ、前記注目
画素を含む複数の画素からなる画素グループおよびこの
画素グループの画素と重複しない複数の画素からなる画
素グループを複数の態様で設定する画素グループ設定手
段と、前記設定した複数の画素グループ態様のうち前記
設定した複数態様の領域を通じてその領域における前記
元画像の構造に最も良く適合する画素グループ態様を選
択する選択手段と、前記選択した画素グループ態様にお
ける前記注目画素を含む画素グループの平均画素値を求
める第１の画素値計算手段と、前記複数態様の領域のう
ち前記選択した画素グループ態様が設定されている領域
の平均画素値を求める第２の画素値計算手段と、前記画
素グループの平均画素値と前記領域の平均画素値とを重
み付け加算する第１の加算手段と、前記重み付け加算し
て得られた画素値と前記注目画素の画素値とを重み付け
加算した画素値を前記注目画素の新たな画素値として画
像を形成する画像形成手段と、前記形成した画像と前記
元画像を重み付け加算する第２の加算手段と、を具備す
ることを特徴とする画像撮影装置である。
【０１６０】この観点での発明では、（８０）と同様に
して求めた画素値によって形成した画像に元画像を重み
付け加算するので、フィルタリングの効果を弱めること
ができる。
【０１６１】（８２）上記の課題を解決するための他の
観点での発明は、前記形成した画像と前記元画像を重み
付け加算する重み係数は調節可能である、ことを特徴と
する（８１）に記載の画像撮影装置である。
【０１６２】この観点での発明では、重み係数の調整に
よりフィルタリングの効果を調節することができる。 （８３）上記の課題を解決するための他の観点での発明
は、前記画素グループの平均画素値と前記領域の平均画
素値とを重み付け加算する重み係数は、前記設定した複
数態様の領域を通じて前記画像グループ態様ごとに求め
た前記画素グループの画素値の残差平方和の総和の最小
値および最大値の関数である、ことを特徴とする（７
９）ないし（８２）のうちのいずれか１つに記載の画像
撮影装置である。
【０１６３】この観点での発明では、（７９）ないし
（８２）のうちのいずれか１つに記載の画像処理におい
て、画素グループの平均画素値と領域の平均画素値とを
重み付け加算する重み係数を、複数態様の領域を通じて
画素グループ態様ごとに求めた画素グループの画素値の
残差平方和の総和の最小値および最大値の関数としたの
で、偽構造の発生しにくいフィルタリングを施した画像
を得ることができる。
【０１６４】（８４）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値と最大値が等しいとき前記領域の平均画素値に対
する重みを最大とし、前記残差平方和の総和の最小値が
最大値に対して小さくなるほど前記領域の平均画素値に
対する重みを小さくする関数である、ことを特徴とする
（８３）に記載の画像撮影装置である。
【０１６５】この観点での発明では、（８３）に記載の
画像処理において、関数は、残差平方和の総和の最小値
と最大値が等しいとき領域の平均画素値に対する重みを
最大とし、残差平方和の総和の最小値が最大値に対して
小さくなるほど領域の平均画素値に対する重みを小さく
する関数であるので、偽構造の発生しにくいフィルタリ
ングを行うことができる。
【０１６６】（８５）上記の課題を解決するための他の
観点での発明は、前記重み付け加算して得られた画素値
と前記注目画素の画素値とを重み付け加算する重み係数
は、前記設定した複数態様の領域を通じて前記画像グル
ープ態様ごとに求めた前記画素グループの画素値の残差
平方和の総和の最小値および前記元画像のノイズの分散
の関数である、ことを特徴とする（８０）ないし（８
４）のうちのいずれか１つに記載の画像撮影装置であ
る。
【０１６７】この観点での発明では、（８０）ないし
（８４）のうちのいずれか１つに記載の画像処理におい
て、重み付け加算して得られた画素値と元画像の注目画
素の画素値とを重み付け加算する重み係数を、複数態様
の領域を通じて画素グループ態様ごとに求めた画素グル
ープの画素値の残差平方和の総和の最小値および元画像
のノイズの分散の関数としたので、元画像の特徴的な構
造を損なわないフィルタリングを施した画像を得ること
ができる。
【０１６８】（８６）上記の課題を解決するための他の
観点での発明は、前記関数は、前記残差平方和の総和の
最小値を前記画素グループ態様における画素数で割った
ものと前記ノイズの分散が等しいとき前記重み付け加算
して得られた画素値に対する重みを最大とし、前記残差
平方和の総和の最小値を前記画素グループ態様における
画素数で割ったものと前記ノイズの分散の違いが大きく
なるほど前記重み付け加算して得られた画素値に対する
重みを小さくする関数である、ことを特徴とする（８
５）に記載の画像撮影装置である。
【０１６９】この観点での発明では、（８５）に記載の
画像処理において、関数は、残差平方和の総和の最小値
を画素グループ態様における画素数で割ったものとノイ
ズの分散が等しいとき重み付け加算して得られた画素値
に対する重みを最大とし、残差平方和の総和の最小値を
画素グループ態様における画素数で割ったものとノイズ
の分散の違いが大きくなるほど重み付け加算して得られ
た画素値に対する重みを小さくする関数であるので、元
画像の特徴的な構造を損なわないフィルタリングを行う
ことができる。
【０１７０】（８７）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記設定し
た複数態様の領域を通じて最小となる画像グループ態様
を選択することによって行う、ことを特徴とする（７
８）ないし（８６）のうちのいずれか１つに記載の画像
撮影装置である。
【０１７１】この観点での発明では、（７８）ないし
（８６）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が複数態様の領域を通じて最小となる
画素グループ態様を選択することによって行うので、元
画像の構造に最も良く適合する画素グループ態様を選択
することができる。
【０１７２】（８８）上記の課題を解決するための他の
観点での発明は、前記画素グループ態様の選択は、前記
画素グループの画素値の残差平方和の総和が前記元画像
のノイズの分散に画素グループ態様における画素数を掛
けたものに最も近くなる画素グループ態様を選択するこ
とによって行う、ことを特徴とする（７８）ないし（８
６）のうちのいずれか１つに記載の画像撮影装置であ
る。
【０１７３】この観点での発明では、（７８）ないし
（８６）のうちのいずれか１つに記載の画像処理におけ
る画素グループ態様の選択は、画素グループの画素値の
残差平方和の総和が元画像のノイズの分散に画素グルー
プ態様における画素数を掛けたものに最も近くなる画素
グループ態様を選択することによって行うので、元画像
の構造に最も良く適合する画素グループ態様を選択する
ことができる。
【０１７４】（８９）上記の課題を解決するための他の
観点での発明は、前記信号は磁気共鳴信号である、こと
を特徴とする（６７）ないし（８８）のうちのいずれか
１つに記載の画像撮影装置である。
【０１７５】この観点での発明では、フィルタリングを
適正に行う画像処理装置を備えた磁気共鳴撮影装置を実
現することができる。
【０１７６】
【発明の実施の形態】以下、図面を参照して本発明の実
施の形態を詳細に説明する。なお、本発明は実施の形態
に限定されるものではない。図１に画像撮影装置のブロ
ック（ｂｌｏｃｋ）図を示す。本装置は本発明の実施の
形態の一例である。本装置の構成によって、本発明の装
置に関する実施の形態の一例が示される。本装置の動作
によって、本発明の方法に関する実施の形態の一例が示
される。
【０１７７】図１に示すように、本装置はマグネットシ
ステム１００を有する。マグネットシステム１００は主
磁場コイル（ｃｏｉｌ）部１０２、勾配コイル部１０６
およびＲＦ（ｒａｄｉｏ ｆｒｅｑｕｅｎｃｙ）コイル
部１０８を有する。これら各コイル部は概ね円筒状の形
状を有し、互いに同軸的に配置されている。マグネット
システム１００の概ね円柱状の内部空間（ボア：ｂｏｒ
ｅ）に、撮影の対象３００がクレードル（ｃｒａｄｌ
ｅ）５００に搭載されて図示しない搬送手段により搬入
および搬出される。
【０１７８】主磁場コイル部１０２はマグネットシステ
ム１００の内部空間に静磁場を形成する。静磁場の方向
は概ね対象３００の体軸の方向に平行である。すなわち
いわゆる水平磁場を形成する。主磁場コイル部１０２は
例えば超伝導コイルを用いて構成される。なお、超伝導
コイルに限らず常伝導コイル等を用いて構成しても良い
のはもちろんである。
【０１７９】勾配コイル部１０６は静磁場強度に勾配を
持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス（ｓｌｉｃｅ）勾配磁場、リードアウト
（ｒｅａｄ ｏｕｔ）勾配磁場およびフェーズエンコー
ド（ｐｈａｓｅ ｅｎｃｏｄｅ）勾配磁場の３種であ
り、これら３種類の勾配磁場に対応して勾配コイル部１
０６は図示しない３系統の勾配コイルを有する。
【０１８０】ＲＦコイル部１０８は静磁場空間に対象３
００の体内のスピンを励起するための高周波磁場を形成
する。以下、高周波磁場を形成することをＲＦ励起信号
の送信ともいう。ＲＦコイル部１０８は、また、励起さ
れたスピンが生じる電磁波すなわち磁気共鳴信号を受信
する。
【０１８１】ＲＦコイル部１０８は図示しない送信用の
コイルおよび受信用のコイルを有する。送信用のコイル
および受信用のコイルは、同じコイルを兼用するかある
いはそれぞれ専用のコイルを用いる。
【０１８２】勾配コイル部１０６には勾配駆動部１３０
が接続されている。勾配駆動部１３０は勾配コイル部１
０６に駆動信号を与えて勾配磁場を発生させる。勾配駆
動部１３０は、勾配コイル部１０６における３系統の勾
配コイルに対応して、図示しない３系統の駆動回路を有
する。
【０１８３】ＲＦコイル部１０８にはＲＦ駆動部１４０
が接続されている。ＲＦ駆動部１４０はＲＦコイル部１
０８に駆動信号を与えてＲＦ励起信号を送信し、対象３
００の体内のスピンを励起する。
【０１８４】ＲＦコイル部１０８にはデータ収集部１５
０が接続されている。データ収集部１５０はＲＦコイル
部１０８が受信した受信信号を取り込み、それをビュー
データ（ｖｉｅｗ ｄａｔａ）として収集する。
【０１８５】勾配駆動部１３０、ＲＦ駆動部１４０およ
びデータ収集部１５０には制御部１６０が接続されてい
る。制御部１６０は、勾配駆動部１３０ないしデータ収
集部１５０をそれぞれ制御して撮影を遂行する。
【０１８６】マグネットシステム１００、勾配駆動部１
３０、ＲＦ駆動部１４０、データ収集部１５０および制
御部１６０からなる部分は、本発明における信号収集手
段の実施の形態の一例である。
【０１８７】データ収集部１５０の出力側はデータ処理
部１７０に接続されている。データ処理部１７０は、例
えばコンピュータ（ｃｏｍｐｕｔｅｒ）等を用いて構成
される。データ処理部１７０は図示しないメモリ（ｍｅ
ｍｏｒｙ）を有する。メモリはデータ処理部１７０用の
プログラムおよび各種のデータを記憶している。本装置
の機能は、データ処理部１７０がメモリに記憶されたプ
ログラムを実行することによりを実現される。
【０１８８】データ処理部１７０は、データ収集部１５
０から取り込んだデータをメモリに記憶する。メモリ内
にはデータ空間が形成される。データ空間は２次元フ−
リエ（Ｆｏｕｒｉｅｒ）空間を構成する。データ処理部
１７０は、これら２次元フ−リエ空間のデータを２次元
逆フ−リエ変換して対象３００の画像を生成（再構成）
する。以下、２次元フ−リエ空間をｋスペース（ｋ−ｓ
ｐａｃｅ）ともいう。データ処理部１７０は、本発明に
おける元画像生成手段の実施の形態の一例である。
【０１８９】データ処理部１７０は、また、再構成した
画像をフィルタリングする機能を有する。データ処理部
１７０は、本発明の画像処理装置の実施の形態の一例で
ある。データ処理部１７０のフィルタリング機能につい
ては後にあらためて説明する。
【０１９０】データ処理部１７０は制御部１６０に接続
されている。データ処理部１７０は制御部１６０の上位
にあってそれを統括する。データ処理部１７０には表示
部１８０および操作部１９０が接続されている。表示部
１８０は、グラフィックディスプレー（ｇｒａｐｈｉｃ
ｄｉｓｐｌａｙ）等で構成される。操作部１９０はポ
インティングデバイス（ｐｏｉｎｔｉｎｇ ｄｅｖｉｃ
ｅ）を備えたキーボード（ｋｅｙｂｏａｒｄ）等で構成
される。
【０１９１】表示部１８０は、データ処理部１７０から
出力される再構成画像および各種の情報を表示する。操
作部１９０は、操作者によって操作され、各種の指令や
情報等をデータ処理部１７０に入力する。操作者は表示
部１８０および操作部１９０を通じてインタラクティブ
（ｉｎｔｅｒａｃｔｉｖｅ）に本装置を操作する。
【０１９２】図２に、他の方式の画像撮影装置のブロッ
ク図を示す。本装置は本発明の実施の形態の一例であ
る。本装置の構成によって、本発明の装置に関する実施
の形態の一例が示される。本装置の動作によって、本発
明の方法に関する実施の形態の一例が示される。
【０１９３】図２に示す装置は、図１に示した装置とは
方式を異にするマグネットシステム１００’を有する。
マグネットシステム１００’以外は図１に示した装置と
同様な構成になっており、同様な部分に同一の符号を付
して説明を省略する。
【０１９４】マグネットシステム１００’は主磁場マグ
ネット部１０２’、勾配コイル部１０６’およびＲＦコ
イル部１０８’を有する。これら主磁場マグネット部１
０２’および各コイル部は、いずれも空間を挟んで互い
に対向する１対のものからなる。また、いずれも概ね円
盤状の形状を有し中心軸を共有して配置されている。マ
グネットシステム１００’の内部空間（ボア）に、対象
３００がクレードル５００に搭載されて図示しない搬送
手段により搬入および搬出される。
【０１９５】主磁場マグネット部１０２’はマグネット
システム１００’の内部空間に静磁場を形成する。静磁
場の方向は概ね対象３００の体軸方向と直交する。すな
わちいわゆる垂直磁場を形成する。主磁場マグネット部
１０２’は例えば永久磁石等を用いて構成される。な
お、永久磁石に限らず超伝導電磁石あるいは常伝導電磁
石等を用いて構成しても良いのはもちろんである。
【０１９６】勾配コイル部１０６’は静磁場強度に勾配
を持たせるための勾配磁場を生じる。発生する勾配磁場
は、スライス勾配磁場、リードアウト勾配磁場およびフ
ェーズエンコード勾配磁場の３種であり、これら３種類
の勾配磁場に対応して勾配コイル部１０６’は図示しな
い３系統の勾配コイルを有する。
【０１９７】ＲＦコイル部１０８’は静磁場空間に対象
３００の体内のスピンを励起するためのＲＦ励起信号を
送信する。ＲＦコイル部１０８’は、また、励起された
スピンが生じる磁気共鳴信号を受信する。ＲＦコイル部
１０８’は図示しない送信用のコイルおよび受信用のコ
イルを有する。送信用のコイルおよび受信用のコイル
は、同じコイルを兼用するかあるいはそれぞれ専用のコ
イルを用いる。
【０１９８】マグネットシステム１００’、勾配駆動部
１３０、ＲＦ駆動部１４０、データ収集部１５０および
制御部１６０からなる部分は、本発明における信号収集
手段の実施の形態の一例である。
【０１９９】図３に、磁気共鳴撮影に用いるパルスシー
ケンス（ｐｕｌｓｅ ｓｅｑｕｅｎｃｅ）の一例を示
す。このパルスシーケンスは、グラディエントエコー
（ＧＲＥ：Ｇｒａｄｉｅｎｔ Ｅｃｈｏ）法のパルスシ
ーケンスである。
【０２００】すなわち、（１）はＧＲＥ法におけるＲＦ
励起用のα°パルスのシーケンスであり、（２）、
（３）、（４）および（５）は、同じくそれぞれ、スラ
イス勾配Ｇｓ、リードアウト勾配Ｇｒ、フェーズエンコ
ード勾配ＧｐおよびグラディエントエコーＭＲのシーケ
ンスである。なお、α°パルスは中心信号で代表する。
パルスシーケンスは時間軸ｔに沿って左から右に進行す
る。
【０２０１】同図に示すように、α°パルスによりスピ
ンのα°励起が行われる。フリップアングル（ｆｌｉｐ
ａｎｇｌｅ）α°は９０°以下である。このときスラ
イス勾配Ｇｓが印加され所定のスライスについての選択
励起が行われる。
【０２０２】α°励起後、フェーズエンコード勾配Ｇｐ
によりスピンのフェーズエンコードが行われる。次に、
リードアウト勾配Ｇｒにより先ずスピンをディフェーズ
（ｄｅｐｈａｓｅ）し、次いでスピンをリフェーズ（ｒ
ｅｐｈａｓｅ）して、グラディエントエコーＭＲを発生
させる。グラディエントエコーＭＲの信号強度は、α°
励起からエコータイム（ｅｃｈｏ ｔｉｍｅ）ＴＥ後の
時点で最大となる。グラディエントエコーＭＲはデータ
収集部１５０によりビューデータとして収集される。
【０２０３】このようなパルスシーケンスが周期ＴＲ
（ｒｅｐｅｔｉｔｉｏｎ ｔｉｍｅ）で６４〜５１２回
繰り返される。繰り返しのたびにフェーズエンコード勾
配Ｇｐを変更し、毎回異なるフェーズエンコードを行
う。これによって、ｋスペースを埋める６４〜５１２ビ
ューのビューデータが得られる。
【０２０４】磁気共鳴撮影用パルスシーケンスの他の例
を図４に示す。このパルスシーケンスは、スピンエコー
（ＳＥ：Ｓｐｉｎ Ｅｃｈｏ）法のパルスシーケンスで
ある。
【０２０５】すなわち、（１）はＳＥ法におけるＲＦ励
起用の９０°パルスおよび１８０°パルスのシーケンス
であり、（２）、（３）、（４）および（５）は、同じ
くそれぞれ、スライス勾配Ｇｓ、リードアウト勾配Ｇ
ｒ、フェーズエンコード勾配ＧｐおよびスピンエコーＭ
Ｒのシーケンスである。なお、９０°パルスおよび１８
０°パルスはそれぞれ中心信号で代表する。パルスシー
ケンスは時間軸ｔに沿って左から右に進行する。
【０２０６】同図に示すように、９０°パルスによりス
ピンの９０°励起が行われる。このときスライス勾配Ｇ
ｓが印加され所定のスライスについての選択励起が行わ
れる。９０°励起から所定の時間後に、１８０°パルス
による１８０°励起すなわちスピン反転が行われる。こ
のときもスライス勾配Ｇｓが印加され、同じスライスに
ついての選択的反転が行われる。
【０２０７】９０°励起とスピン反転の間の期間に、リ
ードアウト勾配Ｇｒおよびフェーズエンコード勾配Ｇｐ
が印加される。リードアウト勾配Ｇｒによりスピンのデ
ィフェーズが行われる。フェーズエンコード勾配Ｇｐに
よりスピンのフェーズエンコードが行われる。
【０２０８】スピン反転後、リードアウト勾配Ｇｒでス
ピンをリフェーズしてスピンエコーＭＲを発生させる。
スピンエコーＭＲの信号強度は、９０°励起からＴＥ後
の時点で最大となる。スピンエコーＭＲはデータ収集部
１５０によりビューデータとして収集される。このよう
なパルスシーケンスが周期ＴＲで６４〜５１２回繰り返
される。繰り返しのたびにフェーズエンコード勾配Ｇｐ
を変更し、毎回異なるフェーズエンコードを行う。これ
によって、ｋスペースを埋める６４〜５１２ビューのビ
ューデータが得られる。
【０２０９】なお、撮影に用いるパルスシーケンスはＧ
ＲＥ法またはＳＥ法に限るものではなく、例えば、ＦＳ
Ｅ（Ｆａｓｔ Ｓｐｉｎ Ｅｃｈｏ）法、ファーストリ
カバリＦＳＥ（Ｆａｓｔ Ｒｅｃｏｖｅｒｙ Ｆａｓｔ
Ｓｐｉｎ Ｅｃｈｏ）法、エコープラナー・イメージ
ング（ＥＰＩ：Ｅｃｈｏ Ｐｌａｎａｒ Ｉｍａｇｉｎ
ｇ）等、他の適宜の技法のものであって良い。
【０２１０】データ処理部１７０は、ｋスペースのビュ
ーデータを２次元逆フ−リエ変換して対象３００の断層
像を再構成する。再構成した画像はメモリに記憶し、ま
た、表示部１８０で表示する。
【０２１１】画像のノイズを除去するために、データ処
理部１７０において画像のフィルタリングが行われる。
フィルタリングは画像再構成の一環として行うようにし
ても良く、あるいは、再構成画像の観察結果に基づい
て、操作者の選択により行うようにしても良い。
【０２１２】図５に、データ処理部１７０による画像フ
ィルタリング動作のフローチャート（ｆｌｏｗ ｃｈａ
ｒｔ）を示す。同図に示すように、ステップ（ｓｔｅ
ｐ）５０２で、元画像における注目画素を指定する。注
目画素とは、これからフィルタリングによって画素値を
確定する画素であり、元画像中の１つの画素が指定され
る。最初の画素としては、例えば元画像の中央画素等が
指定される。
【０２１３】次に、ステップ５０４で、元画像中に局所
領域を設定する。局所領域は注目画素を含む局所的な領
域であり、図６に示すように、注目画素ｋを中心とする
例えば３ｘ３の画素マトリクス（ｍａｔｒｉｘ）が設定
される。なお、局所領域は３ｘ３の画素マトリクスに限
るものではなく、例えば５ｘ５や７ｘ７の画素マトリク
ス等、適宜の領域で良い。ステップ５０４の処理を行う
データ処理部１７０は、本発明における領域設定手段の
実施の形態の一例である。
【０２１４】次に、ステップ５０６で、画素グループ
（ｇｒｏｕｐ）群を設定する。画素グループ群とは複数
の画素グループの集まり（群）であり、画素グループと
は複数の画素の集まり（グループ）である。画素グルー
プは局所領域内の複数の画素によって構成される。グル
ープ間で画素の重複はない。ステップ５０６の処理を行
うデータ処理部１７０は、本発明における画素グループ
設定手段の実施の形態の一例である。
【０２１５】グループ構成は画素の取り合わせに応じて
複数の態様がある。すなわち、図６において、縦に連な
る３つの画素列をＡ１，Ａ２，Ａ３としたとき、画素列
Ａ１，Ａ２，Ａ３をおのおの１つの画素グループとして
３グループ構成とする態様がある。また、画素列Ａ１を
１つのグループとし、画素列Ａ２，Ａ３を一緒にしたも
のを他のグループとして２グループ構成とする態様があ
り、さらに、画素列Ａ１，Ａ２を一緒にしたものを１つ
のグループとし、画素列Ａ３を他のグループとする２グ
ループ構成の態様がある。
【０２１６】また、図７に示すように、横に連なる３つ
の画素列をＢ１，Ｂ２，Ｂ３としたとき、画素列Ｂ１，
Ｂ２，Ｂ３をおのおの１つの画素グループとして３グル
ープ構成とする態様、画素列Ｂ１を１つのグループと
し、画素列Ｂ２，Ｂ３を一緒にしたものを他のグループ
とする２グループ構成の態様、および、画素列Ｂ１，Ｂ
２を一緒にしたものを１つのグループとし、画素列Ｂ３
を他のグループとする２グループ構成の態様がある。
【０２１７】さらに、図８に示すように、右上がりの４
５°方向に連なる１つの画素列をＣ２とし、その両側の
くの字および逆くの字状に連なる２つの画素列をそれぞ
れＣ１，Ｃ３としたとき、画素列Ｃ１，Ｃ２，Ｃ３をお
のおの１つの画素グループとした３グループ構成の態
様、画素列Ｃ１を１つのグループとし、画素列Ｃ２，Ｃ
３を一緒にしたものを他のグループとした２グループ構
成の態様、および、画素列Ｃ１，Ｃ２を一緒にしたもの
を１つのグループとし、画素列Ｃ３を他のグループとし
た２グループ構成の態様がある。
【０２１８】さらにまた、図９に示すように、左上がり
の４５°方向に連なる１つの画素列をＤ２とし、その両
側のくの字および逆くの字状に連なる２つの画素列をそ
れぞれＤ１，Ｄ３としたとき、画素列Ｄ１，Ｄ２，Ｄ３
をおのおの１つの画素グループとした３グループ構成の
態様、画素列Ｄ１を１つのグループとし、画素列Ｄ２，
Ｄ３を一緒にしたものを他のグループとした２グループ
構成の態様、および、画素列Ｄ１，Ｄ２を一緒にしたも
のを１つのグループとし、画素列Ｄ３を他のグループと
した２グループ構成の態様がある。
【０２１９】ステップ５０６では、このような各態様の
グループ構成をそれぞれ画素グループ群として設定す
る。これにより、全体として１２の画素グループ群が設
定される。
【０２２０】なお、ステップ５０４で局所領域を設定す
る際に、４５°方向の画素列のために、例えば図１０に
示すように菱形の画素マトリクスを局所領域として設定
するようにしても良い。なお、この場合は、図６に示し
た方形の画素マトリクスの局所領域に加えて、図１０に
示す菱形の画素マトリクスの局所領域を設定することに
なる。すなわち複数の態様で局所領域を設定することに
なる。この例では態様の数が２となるが、局所領域の大
きさによってはそれ以上の態様があり得る。
【０２２１】このような菱形マトリックスにおいて、右
上がりの４５°方向の３つの画素列をＥ１，Ｅ２，Ｅ３
としたとき、画素列Ｅ１，Ｅ２，Ｅ３をおのおの１つの
画素グループとする３グループ構成の画素グループ群、
画素列Ｅ１を１つのグループとし、画素列Ｅ２，Ｅ３を
一緒にしたものを他のグループとする２グループ構成の
画素グループ群、および、画素列Ｅ１，Ｅ２を一緒にし
たものを１つのグループとし、画素列Ｅ３を他のグルー
プとする２グループ構成の画素グループ群を設定するこ
とができる。
【０２２２】また、図１１に示すように、左上がりの４
５°方向の３つの画素列をＦ１，Ｆ２，Ｆ３としたと
き、画素列Ｆ１，Ｆ２，Ｆ３をおのおの１つの画素グル
ープとした３グループ構成の画素グループ群、画素列Ｆ
１を１つのグループとし、画素列Ｆ２，Ｆ３を一緒にし
たものを他のグループとした２グループ構成の画素グル
ープ群、および、画素列Ｆ１，Ｆ２を一緒にしたものを
１つのグループとし、画素列Ｆ３を他のグループとした
２グループ構成の画素グループ群を設定することができ
る。
【０２２３】次に、ステップ５０８で、画素グループ群
ごとに、画素グループの画素値の残差平方和の総和を計
算する。画素値の残差平方和の総和は、画素グループご
とに画素値の残差平方和をそれぞれ計算し、それらを画
素グループ群内で合算することによって求める。
【０２２４】画素値の残差平方和の総和の計算には次式
が用いられる。
【０２２５】
【数１】
【０２２６】または、
【０２２７】
【数２】
【０２２８】ここで、Ｐｉは画素値、
【０２２９】
【数３】
【０２３０】は、それぞれ、グループＭ１，Ｍ２，Ｍ３
の画素値の平均値である。（１）式は３グループ構成の
画素グループ群に対して用いられ、（２）式は２グルー
プ構成の画素グループ群に対して用いられる。この例で
は、画素グループ群が１２あることにより１２個のＰＶ
が求まる。
【０２３１】次に、ステップ５１０で、ＰＶが最小とな
る画素グループを選択する。すなわち、前のステップで
求めたＰＶの中から値が最小のものを抽出し、この最小
値を与える画素グループ群を特定する。
【０２３２】なお、局所領域が複数の態様で設定されて
いるときは、全ての態様を通じて最小のＰＶを抽出し、
その値を与える画素グループ群を特定する。ステップ５
０８および５１０の処理を行うデータ処理部１７０は、
本発明における選択手段の実施の形態の一例である。
【０２３３】このようにして選択された画素グループ群
は、局所領域における元画像の構造に最も適合した画素
グループ群となる。その理由は次の通りである。例え
ば、図６に示した画素マトリクスにおいて、元画像の輪
郭すなわちエッジ（ｅｄｇｅ）が画素列Ａ２の上にある
とすると、図１２の（ａ）に示すように、画素列Ａ２の
画素値はエッジを表す値を持ち、画素列Ａ１の画素値は
エッジの左側の構造を表す値を持ち、画素列Ａ３の画素
値はエッジの右側の構造を表す値を持つ。ただし、実際
はそれらにノイズが重畳している。
【０２３４】なお、図１２では、各画素列の画素値をそ
れぞれ黒丸で表し、それらの値の相違を上下方向の位置
の相違によって表し、画素列の位置の相違を左右方向の
位置の相違で表す。また、黒丸を結ぶ線によってエッジ
の横断面のプロファイル（ｐｒｏｆｉｌｅ）を表す。
【０２３５】このような状況においては、画素列Ａ１，
Ａ２，Ａ３をそれぞれ１つの画素グループとする画素グ
ループ群が元画像の構造に適合しており、そのＰＶが他
のいずれの画素グループ群のＰＶよりも小さくなる。
【０２３６】なぜなら、このときのＰＶは、画素列Ａ
１，Ａ２，Ａ３の画素値の残差平方和をそれぞれ求めて
合算したものとなるが、画素列Ａ１の画素値は元画像に
おいて構造が同一の部分（エッジの左側の構造）を表す
画素値にノイズが重畳したものであり、画素列Ａ２の画
素値は元画像の構造の他の同一部分（エッジ）を表す画
素値にノイズが重畳したものであり、画素列Ａ３の画素
値は元画像の構造のさらに他の同一部分（エッジの右側
の構造）を表す画素値にノイズが重畳したものであるか
ら、残差平方和はいずれも実質的にノイズの残差平方和
だけとなり、ＰＶも実質的にノイズの残差平方和の総和
となる。
【０２３７】これに対して、他の画素グループ群では、
画素グループ構成が元画像の構造に適合しないため、構
造を異にする部分にまたがって形成された画素グループ
が含まれる。そのような画素グループではノイズに加え
て構造の相違に基づく残差平方和の増加がある。したが
って、画素グループ群のＰＶが増加する。
【０２３８】図６に示した画素マトリクスにおいて、元
画像のエッジが画素列Ａ１に一致しているときは、図１
２の（ｂ）に示すように、画素列Ａ１の画素値はエッジ
を表す値を持ち、画素列Ａ２の画素値および画素列Ａ３
の画素値はエッジの右側の構造を表す値を持つ。ただ
し、実際はそれらにノイズが重畳している。
【０２３９】このような構造に対しては、画素列Ａ１を
１つの画素グループとし、画素列Ａ２，Ａ３を一緒にし
たものを他の画素グループとする２グループ構成の画素
グループ群が元画像の構造に適合し、そのＰＶが最小に
なる。
【０２４０】また、図６に示した画素マトリクスにおい
て、元画像のエッジが画素列Ａ３に一致しているとき
は、図１２の（ｃ）に示すように、画素列Ａ３の画素値
はエッジを表す値を持ち、画素列Ａ１の画素値および画
素列Ａ２の画素値はエッジの左側の構造を表す値を持
つ。ただし、実際はそれらにノイズが重畳している。
【０２４１】このような構造に対しては、画素列Ａ１，
Ａ２を一緒にしたものを１つの画素グループとし、画素
列Ａ３を他の画素グループとする２グループ構成の画素
グループ群が元画像の構造に適合し、そのＰＶが最小に
なる。
【０２４２】エッジの方向が横方向であるときは、図７
に示した画素列Ｂ１，Ｂ２，Ｂ３の組み合わせからなる
３つの画素グループ群のいずれか１つのＰＶが最小にな
る。エッジの方向が右上がりの４５°方向であるとき
は、図８に示した画素列Ｃ１，Ｃ２，Ｃ３の組み合わせ
からなる３つの画素グループ群のいずれか１つ、また
は、図１０に示した画素列Ｅ１，Ｅ２，Ｅ３の組み合わ
せからなる３つの画素グループ群のいずれか１つのＰＶ
が最小になる。
【０２４３】エッジの方向が左上がりの４５°方向であ
るときは、図９に示した画素列Ｄ１，Ｄ２，Ｄ３の組み
合わせからなる３つの画素グループ群のいずれか１つ、
または、図１１に示した画素列Ｆ１，Ｆ２，Ｆ３の組み
合わせからなる３つの画素グループ群のいずれか１つの
ＰＶが最小になる。
【０２４４】画素グループ群が元画像の構造に適合する
ときのＰＶの最小値は、実質的にノイズの残差平方和に
等しくなる。ノイズの残差平方和は、ノイズの分散に画
素グループ群の画素数を掛けたものに相当する。そこ
で、ＰＶについては、最小値の代わりに、ノイズの分散
に画素グループ群の画素数を掛けて得られる値に最も近
い値を抽出し、この値を持つ画素グループ群を選択する
ようにしても良い。
【０２４５】次に、ステップ５１２で、画素グループの
画素の平均値を計算する。平均値の計算に用いる画素
は、ステップ５１０で選択した画素グループ群におけ
る、注目画素ｋを含む画素グループの画素である。
【０２４６】選択された画素グループ群が、例えば、画
素列Ａ１，Ａ２，Ａ３をそれぞれ１つの画素グループと
する３グループ構成の画素グループ群であるときは、注
目画素ｋを含む画素グループは画素列Ａ２であることに
より、画素列Ａ２の平均画素値が計算される。
【０２４７】選択された画素グループ群が、例えば、画
素列Ａ１を１つの画素グループとし、画素列Ａ２，Ａ３
を一緒にしたものを他の画素グループとする２グループ
構成の画素グループ群であるときは、注目画素ｋを含む
画素グループは画素列Ａ２，Ａ３を一緒にしたものであ
ることにより、画素列Ａ２，Ａ３の平均画素値が計算さ
れる。
【０２４８】選択された画素グループ群が、例えば、画
素列Ａ１，Ａ２を一緒にしたものを１つの画素グループ
とし、画素列Ａ３を他の画素グループとする２グループ
構成の画素グループ群であるときは、注目画素ｋを含む
画素グループは画素列Ａ１，Ａ２を一緒にしたものであ
ることにより、画素列Ａ１，Ａ２の平均画素値が計算さ
れる。
【０２４９】選択された画素グループ群が、画素列Ｂ１
〜Ｂ３、Ｃ１〜Ｃ３、Ｄ１〜Ｄ３、Ｅ１〜Ｅ３またはＦ
１〜Ｆ３のいずれかで構成される画素グループからなる
ものである場合も、上記に準じて注目画素ｋを含む画素
グループの平均画素値が求められる。
【０２５０】平均画素値の計算に用いる画素値は局所領
域における元画像の構造に適合した画素グループの画素
値であるから、局所領域における元画像の構造を反映し
た画素値を得ることができる。また複数の画素値の平均
であるからノイズが除去される。すなわち、元画像の構
造を強調しつつノイズを除去した画素値を得ることがで
きる。このようにして求めた平均画素値を、注目画素ｋ
の新たな画素値とする。
【０２５１】このようにして１つの注目画素の画素値を
確定した後に、ステップ５１４で、全ての注目画素につ
いて以上の処理を済ませたか否かを判定し、未済の場合
はステップ５１６で注目画素を変更する。これによって
例えば隣の画素が新たな注目画素となる。
【０２５２】この新たな注目画素について、ステップ５
０４〜５１２の処理が行われ、その画素値が確定する。
以下同様にして、元画像における全ての注目画素が逐一
処理される。全ての注目画素の画素値を確定した後に、
ステップ５１８で、確定済みの画素値Ｐｆ１によって画
像を形成する。
【０２５３】注目画素は元画像を構成する全画素であ
る。なお、必ずしもそれに限るものではなく、必要に応
じて、例えば元画像における関心領域（ＲＯＩ：Ｒｅｇ
ｉｏｎｏｆ Ｉｎｔｅｒｅｓｔ）等、予め設定した領域
の画素であって良い。ステップ５１２および５１８の処
理を行うデータ処理部１７０は、本発明における画像形
成手段の実施の形態の一例である。
【０２５４】このようにして形成した画像は、ノイズが
少なくかつ局所領域における元画像の構造を適正に強調
したものとなる。すなわち、元画像をフィルタリングし
て品質を高めた画像を得ることができる。フィルタリン
グした画像はメモリに記憶し、また、表示部１８０に表
示する。
【０２５５】上記のようなフィルタリングを行った画像
は、従来より可能性は小さいものの、細部構造を強調す
る副作用として、ノイズの影響を受けた画素の組み合わ
せによっては実際には存在しない構造すなわち偽構造を
現出することがあり得る。
【０２５６】そこで、そのような副作用を除去する処理
を付加したフィルタリングを行う。図１３に、そのよう
なフィルタリングの動作のフローチャートを示す。同図
において図５に示した動作と同様な動作は同一の符号を
付して説明を省略する。なお、ステップ５１２の処理を
行うデータ処理部１７０は、本発明における第１の画素
値計算手段の実施の形態の一例である。
【０２５７】ステップ５１２の処理の次に、ステップ５
２０で、ＰＶの最小値と最大値の比ｗ１を求める。すな
わち、
【０２５８】
【数４】
【０２５９】次に、ステップ５２２で、局所領域の画素
値の平均値
【０２６０】
【数５】
【０２６１】を計算する。平均値を求める局所領域は、
ステップ５１０で選択した画素グループ群が存在する局
所領域である。したがって、選択した画素グループ群
が、例えば、図１０または図１１に示した菱形の局所領
域に存在するときは、この領域の平均画素値が求められ
る。ステップ５２２の処理を行うデータ処理部１７０
は、本発明における第２の画素値計算手段の実施の形態
の一例である。
【０２６２】次に、ステップ５２４で、
【０２６３】
【数６】
【０２６４】すなわち、画素グループの平均画素値と、
【０２６５】
【数７】
【０２６６】すなわち領域の画素値の平均値を、次式に
よって重み付け加算する。以下、画素グループの平均画
素値をグループ平均画素値ともいい、局所領域の平均画
素値を領域平均画素値ともいう。また、局所領域におけ
る元画像の構造を、単に元画像の構造ともいう。
【０２６７】
【数８】
【０２６８】ただし、重み係数
【０２６９】
【数９】
【０２７０】はｗ１の単調増加関数であり、例えば図１
４に実線で示すような特性を持つ。なお、特性曲線は図
示のような上に凸の特性に限るものではなく、破線で示
すような下に凸の特性であって良く、または、一点鎖線
で示すような直線であっても良い。
【０２７１】グループ平均画素値と領域平均画素値を重
み付け加算することにより、加算によって得られる画素
値Ｐｆ２においては、グループ平均画素値の寄与度が重
みに応じて変化する。重みはｗ１すなわち画素値のＰＶ
の最小値と最大値の比の関数であり、比が１に近づくほ
どグループ平均画素値の重みが減じて領域平均画素値の
重みが増す。
【０２７２】画素値のＰＶの最小値と最大値の比が１に
近いということは、画素グループ群のどれをとってもＰ
Ｖの値が同じようなものになることを意味し、これは、
領域にエッジ等の目立った構造がなく画素値の相違はノ
イズによるものである可能性が高いことを意味してい
る。
【０２７３】そこで、そのような場合には領域平均画素
値の重みを大きくしてその寄与度を上げるとともに、グ
ループ平均画素値の重みを小さくしてその寄与度を下げ
る。領域に目立った構造がない可能性が高いことによ
り、領域平均画素値の寄与度を上げても画像の鮮鋭さを
損なうおそれはなく、むしろ、グループ平均画素値の寄
与度を下げることによりその副作用を減殺する効果の方
が大きい。
【０２７４】これに対して、画素値のＰＶの最小値と最
大値の比が０に近いということは、領域にはエッジ等の
目立った構造がある可能性が高いことを意味しているの
で、グループ平均画素値の重みを大きくしてその寄与度
を上げるとともに、領域平均画素値の重みを小さくして
その寄与度を下げる。
【０２７５】以上の処理を全ての注目画素について逐一
行い、得られた画素値Ｐｆ２により、ステップ５２４で
画像を形成する。注目画素は元画像を構成する全画素で
ある。なお、必ずしもそれに限るものではなく、必要に
応じて、例えば元画像における関心領域等、予め設定し
た領域の画素であって良い。ステップ５２４および５２
６の処理を行うデータ処理部１７０は、本発明における
画像形成手段の実施の形態の一例である。
【０２７６】このように処理した画像では、グループ平
均画素値の持つ性質を生かしつつその副作用を除去する
ことができる。したがって、元画像をフィルタリングし
て品質をさらに高めた画像を得ることができる。フィル
タリングした画像はメモリに記憶し、また、表示部１８
０に表示する。
【０２７７】以上のような画像処理では、局所領域にお
ける元画像の特徴的な構造が局所領域を通過するエッジ
であることを前提にしている。このため、元画像の構造
がこの前提に合わないときは、グループ平均画素値は元
画像の構造を適正に反映したものにならない。
【０２７８】すなわち、例えば元画像の構造が、図１５
に示すように、縦、横、斜めのいずれの方向において
も、局所領域の一端から他端まで同一の構造が繋がらな
いときは、注目画素ｋを含む画素グループには構造を異
にする部分の画素が必ず含まれる。したがって、そのよ
うな画素グループの平均値は構造を異にする複数の部分
の平均値となるので、この部分に関しては鮮鋭度が低下
する。
【０２７９】このような現象を回避する処理を付加した
フィルタリングのフローチャートを図１６に示す。同図
において図１３に示した処理と同様な処理は同一の符号
を付し説明を省略する。なお、ステップ５２４処理を行
うデータ処理部１７０は、本発明における加算手段の実
施の形態の一例である。
【０２８０】ステップ５２４で画素値Ｐｆ２を求めた後
に、ステップ５２８で、ＰＶの最小値をＮで割ったもの
とノイズの分散との比ｗ２を求める。すなわち、
【０２８１】
【数１０】
【０２８２】なお、ＮはＰＶの最小値を与える画素グル
ープ群における画素の総数である。これにより、ＰＶの
最小値をＮで割ったものは、ＰＶの最小値を与える画素
グループ群の画素値の分散を表す。以下、これを画素値
の分散の最小値という。
【０２８３】次に、ステップ５３０で、画素値Ｐｆ２と
元画像の注目画素ｋの画素値Ｐｏを次式によって重み付
け加算する。以下、元画像の注目画素ｋの画素値を、単
に元画像の画素値という。
【０２８４】
【数１１】
【０２８５】ただし、重み係数はｗ２の関数であり、例
えば、
【０２８６】
【数１２】
【０２８７】で与えられる。このような関数は、図１７
に示すような特性を持つ。すなわち、ｗ２が１のとき関
数の値も１となり、ｗ２が１より大きい範囲ではｗ２が
大きくなるほど関数値は０に向かって低下する。また、
ｗ２が１より小さい範囲ではｗ２が小さくなるほど関数
値が０に向かって低下する。
【０２８８】このような関数を用いてＰｆ２とＰｏを重
み付け加算することにより、画素値の分散の最小値とノ
イズの分散の比が１より大きくなるにつれて、Ｐｆ２の
重みが減じてＰｏの重みが増す。
【０２８９】画素値の分散の最小値とノイズの分散の比
が１に近いということは、ＰＶが最小となる画素グルー
プ群が元画像の構造の適合することを意味する。そのよ
うな場合にはＰｆ２の重みを大きくしてその寄与度を上
げる。
【０２９０】これに対して、画素値の分散の最小値とノ
イズの分散の比が１より大きいということは、ＰＶが最
小となる画素グループ群でさえも元画像の構造の適合し
ないことを意味する。そこで、そのような場合には元画
像の画素値Ｐｏの重みを大きくしてその寄与度を上げる
とともに、Ｐｆ２の重みを小さくしてその寄与度を下げ
る。これによって、元画像の構造を反映した画素値を得
ることができる。
【０２９１】また、画素値の分散の最小値とノイズの分
散の比が１より小さい、すなわち、画素値の分散の最小
値がノイズの分散より小さいということは、通常の画像
ではありえないことであり、ノイズ等に過剰に適合して
いる等の何らかの異常である可能性が高い。そのような
場合にも、Ｐｆ２の重みを小さくしてその寄与度を下げ
るとともに、元画像の画素値Ｐｏの重みを大きくしてそ
の寄与度を上げる。
【０２９２】以上の処理を全ての注目画素について逐一
行い、得られた画素値Ｐｆ３により、ステップ５３２で
画像を形成する。注目画素は元画像を構成する全画素で
ある。なお、必ずしもそれに限るものではなく、必要に
応じて、例えば元画像における関心領域等、予め設定し
た領域の画素であって良い。ステップ５３０および５３
２の処理を行うデータ処理部１７０は、本発明における
画像形成手段の実施の形態の一例である。
【０２９３】この画像では、元画像の特殊な構造を損な
うことなく表現することができる。このようにして、元
画像をより適切にフィルタリングした画像を得ることが
できる。フィルタリングした画像はメモリに記憶し、ま
た、表示部１８０に表示する。
【０２９４】以上の処理によって得られる画像は、一般
的に、フィルタリングが極めて良く利いた画像となる。
画像を観察する便宜上からは、フィルタリングの程度を
適宜に調節可能にすることが要求される。
【０２９５】図１８に、そのような要求に応える処理を
付加したフィルタリングのフローチャートを示す。同図
において図１６に示した処理と同様な処理は同一の符号
を付し説明を省略する。
【０２９６】なお、ステップ５２４の処理を行うデータ
処理部１７０は、本発明における第１の加算手段の実施
の形態の一例である。また、ステップ５３０およびステ
ップ５３２の処理を行うデータ処理部１７０は、本発明
における画像形成手段の実施の形態の一例である。
【０２９７】ステップ５３２の処理の次に、ステップ５
３４で操作者が重み係数αを設定する。次に、ステップ
５３６で、画素値Ｐｆ３と元画像の画素値Ｐｏを、同一
画素同士で次式によって重み付け加算する。
【０２９８】
【数１３】
【０２９９】ステップ５３６の処理を行うデータ処理部
１７０は、本発明における第２の加算手段の実施の形態
の一例である。次に、ステップ５３８で、画素値Ｐｆ４
により画像を形成する。重み係数αを用いてＰｆ３とＰ
ｏを重み付け加算することにより、αの値に応じて最終
画像のフィルタリングの程度を所望の通りに調節するこ
とができる。この画像をメモリに記憶し、また、表示部
１８０に表示する。
【０３００】以上のような機能をコンピュータに実現さ
せるプログラムが、記録媒体に、コンピュータで読み取
り可能なように記録される。記録媒体としては、例え
ば、磁気記録媒体、光記録媒体、光磁気記録媒体および
その他の方式の適宜の記録媒体が用いられる。記録媒体
は半導体記憶媒体であっても良い。本書では記憶媒体は
記録媒体と同義である。
【０３０１】以上、画像のフィルリングを磁気共鳴撮影
装置のデータ処理部で行う例で説明したが、フィルタリ
ングは、例えばＥＷＳ（Ｅｎｇｉｎｅｅｒｉｎｇ Ｗｏ
ｒｋＳｔａｔｉｏｎ）やＰＣ（ｐｅｒｓｏｎａｌ ｃｏ
ｍｐｕｔｅｒ）等、磁気共鳴撮影装置とは別体のデータ
処理装置で行うようにしても良いのはもちろんである。
【０３０２】また、画像撮影装置が磁気共鳴撮影装置で
ある例で説明したが、それに限るものではなく、例えば
Ｘ線ＣＴ（Ｃｏｍｐｕｔｅｄ Ｔｏｍｏｇｒａｐｈｙ）
装置、Ｘ線撮影装置、ＰＥＴ（Ｐｏｓｉｔｒｏｎ Ｅｍ
ｉｓｓｉｏｎ Ｔｏｍｏｇｒａｐｈｙ）、ガンマカメラ
（γ ｃａｍｅｒａ）等、他の方式の画像撮影装置であ
って良い。
【０３０３】また、医用画像を処理する例で説明した
が、処理対象は医用画像に限るものではなく、例えば光
学器械で撮影したディジタル画像等のノイズ除去にも一
般的に適用することができる。
【０３０４】
【発明の効果】以上詳細に説明したように、本発明によ
れば、フィルタリングを適正に行う画像処理方法および
装置、そのような画像処理機能をコンピュータに実現さ
せるプログラムを記録した媒体、並びに、そのような画
像処理装置を備えた画像撮影装置を実現することができ
る。
【図面の簡単な説明】
【図１】本発明の実施の形態の一例の装置のブロック図
である。
【図２】本発明の実施の形態の一例の装置のブロック図
である。
【図３】図１または図２に示した装置が実行するパルス
シーケンスの一例を示す図である。
【図４】図１または図２に示した装置が実行するパルス
シーケンスの一例を示す図である。
【図５】図１または図２に示した装置が行う画像処理の
フローチャートである。
【図６】局所領域の概念図である。
【図７】局所領域の概念図である。
【図８】局所領域の概念図である。
【図９】局所領域の概念図である。
【図１０】局所領域の概念図である。
【図１１】局所領域の概念図である。
【図１２】画素グループ群の概念図である。
【図１３】図１または図２に示した装置が行う画像処理
のフローチャートである。
【図１４】重み係数を示すグラフである。
【図１５】局所領域における元画像の構造の概念図であ
る。
【図１６】図１または図２に示した装置が行う画像処理
のフローチャートである。
【図１７】重み係数を示すグラフである。
【図１８】図１または図２に示した装置が行う画像処理
のフローチャートである。
【符号の説明】
１００，１００’ マグネットシステム １０２ 主磁場コイル部 １０２’ 主磁場マグネット部 １０６，１０６’ 勾配コイル部 １０８，１０８’ ＲＦコイル部 １３０ 勾配駆動部 １４０ ＲＦ駆動部 １５０ データ収集部 １６０ 制御部 １７０ データ処理部 １８０ 表示部 １９０ 操作部 ３００ 対象 ５００ クレードル
───────────────────────────────────────────────────── フロントページの続き (72)発明者 荻野 徹男 東京都日野市旭が丘四丁目７番地の127 ジーイー横河メディカルシステム株式会社 内 Ｆターム(参考） 4C096 AB07 AB08 DC05 DC11 DC27 5B057 BA07 CC03 CE02 CE06 CH08 CH11 CH18

Claims (89)

    【特許請求の範囲】
  1. 【請求項１】 元画像において注目画素を含む局所的な
    領域を設定し、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定し、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を前記注目画素の新たな
    画素値として画像を形成する、ことを特徴とする画像処
    理方法。
  2. 【請求項２】 元画像において注目画素を含む局所的な
    領域を設定し、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定し、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算した画素値を前記注目画素の新たな画
    素値として画像を形成する、ことを特徴とする画像処理
    方法。
  3. 【請求項３】 元画像において注目画素を含む局所的な
    領域を設定し、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定し、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算し、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する、ことを特徴とする画
    像処理方法。
  4. 【請求項４】 元画像において注目画素を含む局所的な
    領域を設定し、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定し、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算し、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成し、 前記形成した画像と前記元画像を重み付け加算する、こ
    とを特徴とする画像処理方法。
  5. 【請求項５】 前記形成した画像と前記元画像を重み付
    け加算する重み係数は調節可能である、ことを特徴とす
    る請求項４に記載の画像処理方法。
  6. 【請求項６】 前記画素グループの平均画素値と前記領
    域の平均画素値とを重み付け加算する重み係数は、 前記画素グループ態様ごとに求めた前記画素グループの
    画素値の残差平方和の総和の最小値および最大値の関数
    である、ことを特徴とする請求項２ないし請求項５のう
    ちのいずれか１つに記載の画像処理方法。
  7. 【請求項７】 前記関数は、 前記残差平方和の総和の最小値と最大値が等しいとき前
    記領域の平均画素値に対する重みを最大とし、前記残差
    平方和の総和の最小値が最大値に対して小さくなるほど
    前記領域の平均画素値に対する重みを小さくする関数で
    ある、ことを特徴とする請求項６に記載の画像処理方
    法。
  8. 【請求項８】 前記重み付け加算して得られた画素値と
    前記注目画素の画素値とを重み付け加算する重み係数
    は、 前記画素グループ態様ごとに求めた前記画素グループの
    画素値の残差平方和の総和の最小値および前記元画像の
    ノイズの分散の関数である、ことを特徴とする請求項３
    ないし請求項７のうちのいずれか１つに記載の画像処理
    方法。
  9. 【請求項９】 前記関数は、 前記残差平方和の総和の最小値を前記画素グループ態様
    における画素数で割ったものと前記ノイズの分散が等し
    いとき前記重み付け加算して得られた画素値に対する重
    みを最大とし、前記残差平方和の総和の最小値を前記画
    素グループ態様における画素数で割ったものと前記ノイ
    ズの分散の違いが大きくなるほど前記重み付け加算して
    得られた画素値に対する重みを小さくする関数である、
    ことを特徴とする請求項８に記載の画像処理方法。
  10. 【請求項１０】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が最小と
    なる画素グループ態様を選択することによって行う、こ
    とを特徴とする請求項１ないし請求項９のうちのいずれ
    か１つに記載の画像処理方法。
  11. 【請求項１１】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記元
    画像のノイズの分散に画素グループ態様における画素数
    を掛けたものに最も近くなる画素グループ態様を選択す
    ることによって行う、 ことを特徴とする請求項１ないし請求項９のうちのいず
    れか１つに記載の画像処理方法。
  12. 【請求項１２】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定し、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定し、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を前記注目画素の新たな
    画素値として画像を形成する、ことを特徴とする画像処
    理方法。
  13. 【請求項１３】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定し、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定し、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算した画素値を前記注目画素の新たな画
    素値として画像を形成する、ことを特徴とする画像処理
    方法。
  14. 【請求項１４】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定し、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定し、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算し、前記重み付け加算して得られた画
    素値と前記注目画素の画素値とを重み付け加算した画素
    値を前記注目画素の新たな画素値として画像を形成す
    る、ことを特徴とする画像処理方法。
  15. 【請求項１５】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定し、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定し、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算し、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成し、 前記形成した画像と前記元画像を重み付け加算する、こ
    とを特徴とする画像処理方法。
  16. 【請求項１６】 前記形成した画像と前記元画像を重み
    付け加算する重み係数は調節可能である、ことを特徴と
    する請求項１５に記載の画像処理方法。
  17. 【請求項１７】 前記画素グループの平均画素値と前記
    領域の平均画素値とを重み付け加算する重み係数は、 前記設定した複数態様の領域を通じて前記画像グループ
    態様ごとに求めた前記画素グループの画素値の残差平方
    和の総和の最小値および最大値の関数である、ことを特
    徴とする請求項１３ないし請求項１６のうちのいずれか
    １つに記載の画像処理方法。
  18. 【請求項１８】 前記関数は、 前記残差平方和の総和の最小値と最大値が等しいとき前
    記領域の平均画素値に対する重みを最大とし、前記残差
    平方和の総和の最小値が最大値に対して小さくなるほど
    前記領域の平均画素値に対する重みを小さくする関数で
    ある、ことを特徴とする請求項１７に記載の画像処理方
    法。
  19. 【請求項１９】 前記重み付け加算して得られた画素値
    と前記注目画素の画素値とを重み付け加算する重み係数
    は、 前記設定した複数態様の領域を通じて前記画像グループ
    態様ごとに求めた前記画素グループの画素値の残差平方
    和の総和の最小値および前記元画像のノイズの分散の関
    数である、ことを特徴とする請求項１４ないし請求項１
    ８のうちのいずれか１つに記載の画像処理方法。
  20. 【請求項２０】 前記関数は、 前記残差平方和の総和の最小値を前記画素グループ態様
    における画素数で割ったものと前記ノイズの分散が等し
    いとき前記重み付け加算して得られた画素値に対する重
    みを最大とし、前記残差平方和の総和の最小値を前記画
    素グループ態様における画素数で割ったものと前記ノイ
    ズの分散の違いが大きくなるほど前記重み付け加算して
    得られた画素値に対する重みを小さくする関数である、
    ことを特徴とする請求項１９に記載の画像処理方法。
  21. 【請求項２１】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記設
    定した複数態様の領域を通じて最小となる画像グループ
    態様を選択することによって行う、 ことを特徴とする請求項１２ないし請求項２０のうちの
    いずれか１つに記載の画像処理方法。
  22. 【請求項２２】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記元
    画像のノイズの分散に画素グループ態様における画素数
    を掛けたものに最も近くなる画素グループ態様を選択す
    ることによって行う、ことを特徴とする請求項１２ない
    し請求項２０のうちのいずれか１つに記載の画像処理方
    法。
  23. 【請求項２３】 元画像において注目画素を含む局所的
    な領域を設定する領域設定手段と、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定する画素グループ設定手段と、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択する選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を前記注目画素の新たな
    画素値として画像を形成する画像形成手段と、を具備す
    ることを特徴とする画像処理装置。
  24. 【請求項２４】 元画像において注目画素を含む局所的
    な領域を設定する領域設定手段と、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定する画素グループ設定手段と、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択する選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記領域の平均画素値を求める第２の画素値計算手段
    と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算した画素値を前記注目画素の新たな画
    素値として画像を形成する画像形成手段と、を具備する
    ことを特徴とする画像処理装置。
  25. 【請求項２５】 元画像において注目画素を含む局所的
    な領域を設定する領域設定手段と、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定する画素グループ設定手段と、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択する選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記領域の平均画素値を求める第２の画素値計算手段
    と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算する加算手段と、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する画像形成手段と、を具
    備することを特徴とする画像処理装置。
  26. 【請求項２６】 元画像において注目画素を含む局所的
    な領域を設定する領域設定手段と、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定する画素グループ設定手段と、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択する選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記領域の平均画素値を求める第２の画素値計算手段
    と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算する第１の加算手段と、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する画像形成手段と、 前記形成した画像と前記元画像を重み付け加算する第２
    の加算手段と、を具備することを特徴とする画像処理装
    置。
  27. 【請求項２７】 前記形成した画像と前記元画像を重み
    付け加算する重み係数は調節可能である、ことを特徴と
    する請求項２６に記載の画像処理装置。
  28. 【請求項２８】 前記画素グループの平均画素値と前記
    領域の平均画素値とを重み付け加算する重み係数は、 前記画素グループ態様ごとに求めた前記画素グループの
    画素値の残差平方和の総和の最小値および最大値の関数
    である、ことを特徴とする請求項２４ないし請求項２７
    のうちのいずれか１つに記載の画像処理装置。
  29. 【請求項２９】 前記関数は、 前記残差平方和の総和の最小値と最大値が等しいとき前
    記領域の平均画素値に対する重みを最大とし、前記残差
    平方和の総和の最小値が最大値に対して小さくなるほど
    前記領域の平均画素値に対する重みを小さくする関数で
    ある、ことを特徴とする請求項２８に記載の画像処理装
    置。
  30. 【請求項３０】 前記重み付け加算して得られた画素値
    と前記注目画素の画素値とを重み付け加算する重み係数
    は、 前記画素グループ態様ごとに求めた前記画素グループの
    画素値の残差平方和の総和の最小値および前記元画像の
    ノイズの分散の関数である、ことを特徴とする請求項２
    ５ないし請求項２９のうちのいずれか１つに記載の画像
    処理装置。
  31. 【請求項３１】 前記関数は、 前記残差平方和の総和の最小値を前記画素グループ態様
    における画素数で割ったものと前記ノイズの分散が等し
    いとき前記重み付け加算して得られた画素値に対する重
    みを最大とし、前記残差平方和の総和の最小値を前記画
    素グループ態様における画素数で割ったものと前記ノイ
    ズの分散の違いが大きくなるほど前記重み付け加算して
    得られた画素値に対する重みを小さくする関数である、
    ことを特徴とする請求項３０に記載の画像処理装置。
  32. 【請求項３２】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が最小と
    なる画素グループ態様を選択することによって行う、こ
    とを特徴とする請求項２３ないし請求項３１のうちのい
    ずれか１つに記載の画像処理装置。
  33. 【請求項３３】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記元
    画像のノイズの分散に画素グループ態様における画素数
    を掛けたものに最も近くなる画素グループ態様を選択す
    ることによって行う、ことを特徴とする請求項２３ない
    し請求項３１のうちのいずれか１つに記載の画像処理装
    置。
  34. 【請求項３４】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定する領域設定手段と、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定する画素グループ設
    定手段と、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択する
    選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を前記注目画素の新たな
    画素値として画像を形成する画像形成手段と、を具備す
    ることを特徴とする画像処理装置。
  35. 【請求項３５】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定する領域設定手段と、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定する画素グループ設
    定手段と、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択する
    選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求める第２の画
    素値計算手段と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算した画素値を前記注目画素の新たな画
    素値として画像を形成する画像形成手段と、を具備する
    ことを特徴とする画像処理装置。
  36. 【請求項３６】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定する領域設定手段と、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定する画素グループ設
    定手段と、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択する
    選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求める第２の画
    素値計算手段と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算する加算手段と、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する画像形成手段と、を具
    備することを特徴とする画像処理装置。
  37. 【請求項３７】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定する領域設定手段と、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定する画素グループ設
    定手段と、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択する
    選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求める第２の画
    素値計算手段と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算する第１の加算手段と、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する画像形成手段と、 前記形成した画像と前記元画像を重み付け加算する第２
    の加算手段と、 を具備することを特徴とする画像処理装置。
  38. 【請求項３８】 前記形成した画像と前記元画像を重み
    付け加算する重み係数は調節可能である、ことを特徴と
    する請求項３７に記載の画像処理装置。
  39. 【請求項３９】 前記画素グループの平均画素値と前記
    領域の平均画素値とを重み付け加算する重み係数は、 前記設定した複数態様の領域を通じて前記画像グループ
    態様ごとに求めた前記画素グループの画素値の残差平方
    和の総和の最小値および最大値の関数である、ことを特
    徴とする請求項３５ないし請求項３８のうちのいずれか
    １つに記載の画像処理装置。
  40. 【請求項４０】 前記関数は、 前記残差平方和の総和の最小値と最大値が等しいとき前
    記領域の平均画素値に対する重みを最大とし、前記残差
    平方和の総和の最小値が最大値に対して小さくなるほど
    前記領域の平均画素値に対する重みを小さくする関数で
    ある、ことを特徴とする請求項３９に記載の画像処理装
    置。
  41. 【請求項４１】 前記重み付け加算して得られた画素値
    と前記注目画素の画素値とを重み付け加算する重み係数
    は、 前記設定した複数態様の領域を通じて前記画像グループ
    態様ごとに求めた前記画素グループの画素値の残差平方
    和の総和の最小値および前記元画像のノイズの分散の関
    数である、ことを特徴とする請求項３６ないし請求項４
    ０のうちのいずれか１つに記載の画像処理装置。
  42. 【請求項４２】 前記関数は、 前記残差平方和の総和の最小値を前記画素グループ態様
    における画素数で割ったものと前記ノイズの分散が等し
    いとき前記重み付け加算して得られた画素値に対する重
    みを最大とし、前記残差平方和の総和の最小値を前記画
    素グループ態様における画素数で割ったものと前記ノイ
    ズの分散の違いが大きくなるほど前記重み付け加算して
    得られた画素値に対する重みを小さくする関数である、
    ことを特徴とする請求項４１に記載の画像処理装置。
  43. 【請求項４３】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記設
    定した複数態様の領域を通じて最小となる画像グループ
    態様を選択することによって行う、ことを特徴とする請
    求項３４ないし請求項４２のうちのいずれか１つに記載
    の画像処理装置。
  44. 【請求項４４】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記元
    画像のノイズの分散に画素グループ態様における画素数
    を掛けたものに最も近くなる画素グループ態様を選択す
    ることによって行う、 ことを特徴とする請求項３４ないし請求項４２のうちの
    いずれか１つに記載の画像処理装置。
  45. 【請求項４５】 元画像において注目画素を含む局所的
    な領域を設定し、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定し、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を前記注目画素の新たな
    画素値として画像を形成する、機能をコンピュータに実
    現させるプログラムをコンピュータにより読み取り可能
    なように記録したことを特徴とする記録媒体。
  46. 【請求項４６】 元画像において注目画素を含む局所的
    な領域を設定し、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定し、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算した画素値を前記注目画素の新たな画
    素値として画像を形成する、 機能をコンピュータに実現させるプログラムをコンピュ
    ータにより読み取り可能なように記録したことを特徴と
    する記録媒体。
  47. 【請求項４７】 元画像において注目画素を含む局所的
    な領域を設定し、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定し、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算し、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する、 機能をコンピュータに実現させるプログラムをコンピュ
    ータにより読み取り可能なように記録したことを特徴と
    する記録媒体。
  48. 【請求項４８】 元画像において注目画素を含む局所的
    な領域を設定し、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定し、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算し、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成し、 前記形成した画像と前記元画像を重み付け加算する、 機能をコンピュータに実現させるプログラムをコンピュ
    ータにより読み取り可能なように記録したことを特徴と
    する記録媒体。
  49. 【請求項４９】 前記形成した画像と前記元画像を重み
    付け加算する重み係数は調節可能である、ことを特徴と
    する請求項４８に記載の記録媒体。
  50. 【請求項５０】 前記画素グループの平均画素値と前記
    領域の平均画素値とを重み付け加算する重み係数は、 前記画素グループ態様ごとに求めた前記画素グループの
    画素値の残差平方和の総和の最小値および最大値の関数
    である、ことを特徴とする請求項４６ないし請求項４９
    のうちのいずれか１つに記載の記録媒体。
  51. 【請求項５１】 前記関数は、 前記残差平方和の総和の最小値と最大値が等しいとき前
    記領域の平均画素値に対する重みを最大とし、前記残差
    平方和の総和の最小値が最大値に対して小さくなるほど
    前記領域の平均画素値に対する重みを小さくする関数で
    ある、ことを特徴とする請求項５０に記載の記録媒体。
  52. 【請求項５２】 前記重み付け加算して得られた画素値
    と前記注目画素の画素値とを重み付け加算する重み係数
    は、 前記画素グループ態様ごとに求めた前記画素グループの
    画素値の残差平方和の総和の最小値および前記元画像の
    ノイズの分散の関数である、ことを特徴とする請求項４
    ７ないし請求項５１のうちのいずれか１つに記載の記録
    媒体。
  53. 【請求項５３】 前記関数は、 前記残差平方和の総和の最小値を前記画素グループ態様
    における画素数で割ったものと前記ノイズの分散が等し
    いとき前記重み付け加算して得られた画素値に対する重
    みを最大とし、前記残差平方和の総和の最小値を前記画
    素グループ態様における画素数で割ったものと前記ノイ
    ズの分散の違いが大きくなるほど前記重み付け加算して
    得られた画素値に対する重みを小さくする関数である、
    ことを特徴とする請求項５２に記載の記録媒体。
  54. 【請求項５４】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が最小と
    なる画素グループ態様を選択することによって行う、こ
    とを特徴とする請求項４５ないし請求項５３のうちのい
    ずれか１つに記載の記録媒体。
  55. 【請求項５５】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記元
    画像のノイズの分散に画素グループ態様における画素数
    を掛けたものに最も近くなる画素グループ態様を選択す
    ることによって行う、ことを特徴とする請求項４５ない
    し請求項５３のうちのいずれか１つに記載の記録媒体。
  56. 【請求項５６】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定し、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定し、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を前記注目画素の新たな
    画素値として画像を形成する、機能をコンピュータに実
    現させるプログラムをコンピュータにより読み取り可能
    なように記録したことを特徴とする記録媒体。
  57. 【請求項５７】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定し、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定し、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算した画素値を前記注目画素の新たな画
    素値として画像を形成する、機能をコンピュータに実現
    させるプログラムをコンピュータにより読み取り可能な
    ように記録したことを特徴とする記録媒体。
  58. 【請求項５８】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定し、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定し、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算し、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する、機能をコンピュータ
    に実現させるプログラムをコンピュータにより読み取り
    可能なように記録したことを特徴とする記録媒体。
  59. 【請求項５９】 元画像において注目画素を含む局所的
    な領域を複数の態様で設定し、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定し、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択し、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求め、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求め、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算し、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成し、 前記形成した画像と前記元画像を重み付け加算する、機
    能をコンピュータに実現させるプログラムをコンピュー
    タにより読み取り可能なように記録したことを特徴とす
    る記録媒体。
  60. 【請求項６０】 前記形成した画像と前記元画像を重み
    付け加算する重み係数は調節可能である、ことを特徴と
    する請求項５９に記載の記録媒体。
  61. 【請求項６１】 前記画素グループの平均画素値と前記
    領域の平均画素値とを重み付け加算する重み係数は、 前記設定した複数態様の領域を通じて前記画像グループ
    態様ごとに求めた前記画素グループの画素値の残差平方
    和の総和の最小値および最大値の関数である、ことを特
    徴とする請求項５７ないし請求項６０のうちのいずれか
    １つに記載の記録媒体。
  62. 【請求項６２】 前記関数は、 前記残差平方和の総和の最小値と最大値が等しいとき前
    記領域の平均画素値に対する重みを最大とし、前記残差
    平方和の総和の最小値が最大値に対して小さくなるほど
    前記領域の平均画素値に対する重みを小さくする関数で
    ある、ことを特徴とする請求項６１に記載の記録媒体。
  63. 【請求項６３】 前記重み付け加算して得られた画素値
    と前記注目画素の画素値とを重み付け加算する重み係数
    は、 前記設定した複数態様の領域を通じて前記画像グループ
    態様ごとに求めた前記画素グループの画素値の残差平方
    和の総和の最小値および前記元画像のノイズの分散の関
    数である、ことを特徴とする請求項５８ないし請求項６
    ２のうちのいずれか１つに記載の記録媒体。
  64. 【請求項６４】 前記関数は、 前記残差平方和の総和の最小値を前記画素グループ態様
    における画素数で割ったものと前記ノイズの分散が等し
    いとき前記重み付け加算して得られた画素値に対する重
    みを最大とし、前記残差平方和の総和の最小値を前記画
    素グループ態様における画素数で割ったものと前記ノイ
    ズの分散の違いが大きくなるほど前記重み付け加算して
    得られた画素値に対する重みを小さくする関数である、
    ことを特徴とする請求項６３に記載の記録媒体。
  65. 【請求項６５】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記設
    定した複数態様の領域を通じて最小となる画像グループ
    態様を選択することによって行う、ことを特徴とする請
    求項５６ないし請求項６４のうちのいずれか１つに記載
    の記録媒体。
  66. 【請求項６６】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記元
    画像のノイズの分散に画素グループ態様における画素数
    を掛けたものに最も近くなる画素グループ態様を選択す
    ることによって行う、ことを特徴とする請求項５６ない
    し請求項６４のうちのいずれか１つに記載の記録媒体。
  67. 【請求項６７】 対象から信号を収集する信号収集手段
    と、 前記収集した信号に基づいて元画像を生成する元画像生
    成手段と、 前記元画像において注目画素を含む局所的な領域を設定
    する領域設定手段と、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定する画素グループ設定手段と、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択する選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を前記注目画素の新たな
    画素値として画像を形成する画像形成手段と、を具備す
    ることを特徴とする画像撮影装置。
  68. 【請求項６８】 対象から信号を収集する信号収集手段
    と、 前記収集した信号に基づいて元画像を生成する元画像生
    成手段と、 前記元画像において注目画素を含む局所的な領域を設定
    する領域設定手段と、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定する画素グループ設定手段と、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択する選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記領域の平均画素値を求める第２の画素値計算手段
    と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算した画素値を前記注目画素の新たな画
    素値として画像を形成する画像形成手段と、を具備する
    ことを特徴とする画像撮影装置。
  69. 【請求項６９】 対象から信号を収集する信号収集手段
    と、 前記収集した信号に基づいて元画像を生成する元画像生
    成手段と、 前記元画像において注目画素を含む局所的な領域を設定
    する領域設定手段と、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定する画素グループ設定手段と、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択する選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記領域の平均画素値を求める第２の画素値計算手段
    と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算する加算手段と、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する画像形成手段と、を具
    備することを特徴とする画像撮影装置。
  70. 【請求項７０】 対象から信号を収集する信号収集手段
    と、 前記収集した信号に基づいて元画像を生成する元画像生
    成手段と、 前記元画像において注目画素を含む局所的な領域を設定
    する領域設定手段と、 前記領域において前記注目画素を含む複数の画素からな
    る画素グループおよびこの画素グループの画素と重複し
    ない複数の画素からなる画素グループを複数の態様で設
    定する画素グループ設定手段と、 前記設定した複数の画素グループ態様のうち前記領域に
    おける前記元画像の構造に最も良く適合する画素グルー
    プ態様を選択する選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記領域の平均画素値を求める第２の画素値計算手段
    と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算する第１の加算手段と、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する画像形成手段と、 前記形成した画像と前記元画像を重み付け加算する第２
    の加算手段と、 を具備することを特徴とする画像撮影装置。
  71. 【請求項７１】 前記形成した画像と前記元画像を重み
    付け加算する重み係数は調節可能である、 ことを特徴とする請求項７０に記載の画像撮影装置。
  72. 【請求項７２】 前記画素グループの平均画素値と前記
    領域の平均画素値とを重み付け加算する重み係数は、 前記画素グループ態様ごとに求めた前記画素グループの
    画素値の残差平方和の総和の最小値および最大値の関数
    である、ことを特徴とする請求項６８ないし請求項７１
    のうちのいずれか１つに記載の画像撮影装置。
  73. 【請求項７３】 前記関数は、 前記残差平方和の総和の最小値と最大値が等しいとき前
    記領域の平均画素値に対する重みを最大とし、前記残差
    平方和の総和の最小値が最大値に対して小さくなるほど
    前記領域の平均画素値に対する重みを小さくする関数で
    ある、ことを特徴とする請求項７２に記載の画像撮影装
    置。
  74. 【請求項７４】 前記重み付け加算して得られた画素値
    と前記注目画素の画素値とを重み付け加算する重み係数
    は、 前記画素グループ態様ごとに求めた前記画素グループの
    画素値の残差平方和の総和の最小値および前記元画像の
    ノイズの分散の関数である、ことを特徴とする請求項６
    ９ないし請求項７３のうちのいずれか１つに記載の画像
    撮影装置。
  75. 【請求項７５】 前記関数は、 前記残差平方和の総和の最小値を前記画素グループ態様
    における画素数で割ったものと前記ノイズの分散が等し
    いとき前記重み付け加算して得られた画素値に対する重
    みを最大とし、前記残差平方和の総和の最小値を前記画
    素グループ態様における画素数で割ったものと前記ノイ
    ズの分散の違いが大きくなるほど前記重み付け加算して
    得られた画素値に対する重みを小さくする関数である、
    ことを特徴とする請求項７４に記載の画像撮影装置。
  76. 【請求項７６】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が最小と
    なる画素グループ態様を選択することによって行う、こ
    とを特徴とする請求項６７ないし請求項７５のうちのい
    ずれか１つに記載の画像撮影装置。
  77. 【請求項７７】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記元
    画像のノイズの分散に画素グループ態様における画素数
    を掛けたものに最も近くなる画素グループ態様を選択す
    ることによって行う、ことを特徴とする請求項６７ない
    し請求項７５のうちのいずれか１つに記載の画像撮影装
    置。
  78. 【請求項７８】 対象から信号を収集する信号収集手段
    と、 前記収集した信号に基づいて元画像を生成する元画像生
    成手段と、 前記元画像において注目画素を含む局所的な領域を複数
    の態様で設定する領域設定手段と、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定する画素グループ設
    定手段と、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択する
    選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を前記注目画素の新たな
    画素値として画像を形成する画像形成手段と、を具備す
    ることを特徴とする画像撮影装置。
  79. 【請求項７９】 対象から信号を収集する信号収集手段
    と、 前記収集した信号に基づいて元画像を生成する元画像生
    成手段と、 前記元画像において注目画素を含む局所的な領域を複数
    の態様で設定する領域設定手段と、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定する画素グループ設
    定手段と、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択する
    選択手段と、前記選択した画素グループ態様における前
    記注目画素を含む画素グループの平均画素値を求める第
    １の画素値計算手段と、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求める第２の画
    素値計算手段と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算した画素値を前記注目画素の新たな画
    素値として画像を形成する画像形成手段と、を具備する
    ことを特徴とする画像撮影装置。
  80. 【請求項８０】 対象から信号を収集する信号収集手段
    と、 前記収集した信号に基づいて元画像を生成する元画像生
    成手段と、 前記元画像において注目画素を含む局所的な領域を複数
    の態様で設定する領域設定手段と、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定する画素グループ設
    定手段と、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択する
    選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求める第２の画
    素値計算手段と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算する加算手段と、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する画像形成手段と、を具
    備することを特徴とする画像撮影装置。
  81. 【請求項８１】 対象から信号を収集する信号収集手段
    と、 前記収集した信号に基づいて元画像を生成する元画像生
    成手段と、 前記元画像において注目画素を含む局所的な領域を複数
    の態様で設定する領域設定手段と、 前記設定した複数態様の領域において、それぞれ、前記
    注目画素を含む複数の画素からなる画素グループおよび
    この画素グループの画素と重複しない複数の画素からな
    る画素グループを複数の態様で設定する画素グループ設
    定手段と、 前記設定した複数の画素グループ態様のうち前記設定し
    た複数態様の領域を通じてその領域における前記元画像
    の構造に最も良く適合する画素グループ態様を選択する
    選択手段と、 前記選択した画素グループ態様における前記注目画素を
    含む画素グループの平均画素値を求める第１の画素値計
    算手段と、 前記複数態様の領域のうち前記選択した画素グループ態
    様が設定されている領域の平均画素値を求める第２の画
    素値計算手段と、 前記画素グループの平均画素値と前記領域の平均画素値
    とを重み付け加算する第１の加算手段と、 前記重み付け加算して得られた画素値と前記注目画素の
    画素値とを重み付け加算した画素値を前記注目画素の新
    たな画素値として画像を形成する画像形成手段と、 前記形成した画像と前記元画像を重み付け加算する第２
    の加算手段と、を具備することを特徴とする画像撮影装
    置。
  82. 【請求項８２】 前記形成した画像と前記元画像を重み
    付け加算する重み係数は調節可能である、ことを特徴と
    する請求項８１に記載の画像撮影装置。
  83. 【請求項８３】 前記画素グループの平均画素値と前記
    領域の平均画素値とを重み付け加算する重み係数は、 前記設定した複数態様の領域を通じて前記画像グループ
    態様ごとに求めた前記画素グループの画素値の残差平方
    和の総和の最小値および最大値の関数である、ことを特
    徴とする請求項７９ないし請求項８２のうちのいずれか
    １つに記載の画像撮影装置。
  84. 【請求項８４】 前記関数は、前記残差平方和の総和の
    最小値と最大値が等しいとき前記領域の平均画素値に対
    する重みを最大とし、前記残差平方和の総和の最小値が
    最大値に対して小さくなるほど前記領域の平均画素値に
    対する重みを小さくする関数である、ことを特徴とする
    請求項８３に記載の画像撮影装置。
  85. 【請求項８５】 前記重み付け加算して得られた画素値
    と前記注目画素の画素値とを重み付け加算する重み係数
    は、 前記設定した複数態様の領域を通じて前記画像グループ
    態様ごとに求めた前記画素グループの画素値の残差平方
    和の総和の最小値および前記元画像のノイズの分散の関
    数である、ことを特徴とする請求項８０ないし請求項８
    ４のうちのいずれか１つに記載の画像撮影装置。
  86. 【請求項８６】 前記関数は、 前記残差平方和の総和の最小値を前記画素グループ態様
    における画素数で割ったものと前記ノイズの分散が等し
    いとき前記重み付け加算して得られた画素値に対する重
    みを最大とし、前記残差平方和の総和の最小値を前記画
    素グループ態様における画素数で割ったものと前記ノイ
    ズの分散の違いが大きくなるほど前記重み付け加算して
    得られた画素値に対する重みを小さくする関数である、
    ことを特徴とする請求項８５に記載の画像撮影装置。
  87. 【請求項８７】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記設
    定した複数態様の領域を通じて最小となる画像グループ
    態様を選択することによって行う、 ことを特徴とする請求項７８ないし請求項８６のうちの
    いずれか１つに記載の画像撮影装置。
  88. 【請求項８８】 前記画素グループ態様の選択は、 前記画素グループの画素値の残差平方和の総和が前記元
    画像のノイズの分散に画素グループ態様における画素数
    を掛けたものに最も近くなる画素グループ態様を選択す
    ることによって行う、ことを特徴とする請求項７８ない
    し請求項８６のうちのいずれか１つに記載の画像撮影装
    置。
  89. 【請求項８９】 前記信号は磁気共鳴信号である、こと
    を特徴とする請求項６７ないし請求項８８のうちのいず
    れか１つに記載の画像撮影装置。
JP2000220416A 2000-07-21 2000-07-21 画像処理方法および装置、記録媒体並びに画像撮影装置 Expired - Fee Related JP4012669B2 (ja)

Priority Applications (7)

Application Number Priority Date Filing Date Title
JP2000220416A JP4012669B2 (ja) 2000-07-21 2000-07-21 画像処理方法および装置、記録媒体並びに画像撮影装置
US09/863,998 US7061651B2 (en) 2000-07-21 2001-05-23 Image processing method and apparatus, recording medium, and imaging apparatus
TW090116777A TWI221263B (en) 2000-07-21 2001-07-09 Image processing method and apparatus, recording medium, and imaging apparatus
DE60133694T DE60133694T2 (de) 2000-07-21 2001-07-18 Bildverarbeitungsvorrichtung und -verfahren, Aufzeichnungsmedium und Bildgebungsgerät
EP01306159A EP1176555B1 (en) 2000-07-21 2001-07-18 Image processing method and apparatus, recording medium, and imaging apparatus
KR10-2001-0043809A KR100452644B1 (ko) 2000-07-21 2001-07-20 영상 처리 방법과 장치, 기록 매체 및 촬상 장치
CNB011253886A CN1262237C (zh) 2000-07-21 2001-07-21 图像处理方法和设备、记录介质以及成像设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2000220416A JP4012669B2 (ja) 2000-07-21 2000-07-21 画像処理方法および装置、記録媒体並びに画像撮影装置

Publications (3)

Publication Number Publication Date
JP2002042128A true JP2002042128A (ja) 2002-02-08
JP2002042128A5 JP2002042128A5 (ja) 2005-05-26
JP4012669B2 JP4012669B2 (ja) 2007-11-21

Family

ID=18715028

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000220416A Expired - Fee Related JP4012669B2 (ja) 2000-07-21 2000-07-21 画像処理方法および装置、記録媒体並びに画像撮影装置

Country Status (7)

Country Link
US (1) US7061651B2 (ja)
EP (1) EP1176555B1 (ja)
JP (1) JP4012669B2 (ja)
KR (1) KR100452644B1 (ja)
CN (1) CN1262237C (ja)
DE (1) DE60133694T2 (ja)
TW (1) TWI221263B (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8537893B2 (en) 2006-01-23 2013-09-17 Max-Planck-Gesellschaft Zur Forderung Der Wissenschaften E.V. High dynamic range codecs
US7429862B2 (en) * 2006-04-18 2008-09-30 Siemens Aktiengesellschaft PPA MRT method and apparatus with scaled contribution of k-space data to the reconstructed image
US9639915B2 (en) * 2012-08-08 2017-05-02 Samsung Electronics Co., Ltd. Image processing method and apparatus
US11300646B2 (en) 2014-11-10 2022-04-12 Canon Medical Systems Corporation Magnetic resonance imaging apparatus, image processing apparatus, and image processing method
GB2567881B (en) 2017-10-30 2021-02-10 Imagination Tech Ltd Systems and methods for processing a stream of data values
GB2568038B (en) 2017-10-30 2020-12-02 Imagination Tech Ltd Systems and methods for processing a stream of data values
CN111436905B (zh) * 2020-04-03 2024-03-19 图湃(北京)医疗科技有限公司 一种光学相干层析血管造影方法及装置

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2737034A1 (fr) * 1995-07-21 1997-01-24 Philips Electronique Lab Procede de traitement spatial d'une image numerique pour la reduction du bruit, et dispositif mettant en oeuvre ce procede
US5771318A (en) * 1996-06-27 1998-06-23 Siemens Corporate Research, Inc. Adaptive edge-preserving smoothing filter

Also Published As

Publication number Publication date
EP1176555A3 (en) 2004-12-08
US20020026947A1 (en) 2002-03-07
DE60133694D1 (de) 2008-06-05
CN1262237C (zh) 2006-07-05
TWI221263B (en) 2004-09-21
KR20020009443A (ko) 2002-02-01
DE60133694T2 (de) 2009-05-07
CN1342441A (zh) 2002-04-03
EP1176555B1 (en) 2008-04-23
JP4012669B2 (ja) 2007-11-21
US7061651B2 (en) 2006-06-13
KR100452644B1 (ko) 2004-10-12
EP1176555A2 (en) 2002-01-30

Similar Documents

Publication Publication Date Title
JP4191884B2 (ja) 画像処理方法、画像処理装置および画像撮影装置
CN103477238B (zh) 采用来自先验采集的约束的压缩感测mr图像重建
JP2750866B2 (ja) 磁気共鳴画像化装置
JP3976684B2 (ja) 画像における動きの影響を低減する方法および装置
EP1225544A2 (en) Image processing method and apparatus, recording medium and imaging apparatus
JPH07148133A (ja) アンジオグラフィ方法及び装置
JPH1040365A (ja) 磁気共鳴映像法および装置
JP4070268B2 (ja) 画像の分解能を向上させる方法、投影画像の一部を拡大させる方法及び画像の分解能を向上させる装置
KR100444094B1 (ko) 이미지 처리 장치 및 그 방법, 촬상 장치 및 기록 매체
CN105143906A (zh) 金属抗mr成像
US6486667B1 (en) Combination of fluid-attenuated inversion-recovery complex images acquired using magnetic resonance imaging
US20240138700A1 (en) Medical image processing apparatus, method of medical image processing, and nonvolatile computer readable storage medium storing therein medical image processing program
JP2002042128A (ja) 画像処理方法および装置、記録媒体並びに画像撮影装置
JPH05123312A (ja) 磁気共鳴イメージング装置
JP3884227B2 (ja) 磁気共鳴撮影装置
JP5289011B2 (ja) 磁気共鳴イメージング装置および画像再構成方法
JP4363625B2 (ja) 画像処理方法および磁気共鳴撮影装置
JP2019037575A (ja) 磁気共鳴撮像装置及び画像処理方法
JP3023613B2 (ja) Ｍｒｉ装置
JP5566587B2 (ja) 磁気共鳴イメージング装置およびリンギング低減方法
Garach Robust phase sensitive inversion recovery imaging
JP2005137782A (ja) 核磁気共鳴撮像装置およびそれを用いたパルスシーケンス設定方法

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20040726

A625 Written request for application examination (by other person)

Free format text: JAPANESE INTERMEDIATE CODE: A625

Effective date: 20040811

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20070410

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070629

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: 20070828

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070910

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: 20100914

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100914

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20100914

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110914

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110914

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120914

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120914

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120914

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130914

Year of fee payment: 6

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees