JPS62267880A - 像を補正する装置と方法 - Google Patents

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 発明の背景 この発明は作像装置、更に具体的に云えば、像にわたっ
て起る系統的な輝度変動に対して像を）ｉ（ｉ正するこ
とに関する。

この発明は任意の方法によって作られた像を補正する為
に用いることが出来、この発明の範囲もその様に解釈さ
れるべきであるが、説明を具体的にする為に、核磁気共
鳴装置によって発生される像を補正する場合についてこ
の発明を説明する。

核磁気共鳴装置は、１つ又は更に多くの一磁気勾配を重
畳した一定磁界を物体にかけ、物体内の原子核を、正し
く選ばれた無線周波数帯で励振することにより、物体の
像を作る。例えば水素の原子核の磁気モーメントが、予
定の長さの無線周波パルスによって横平面へ回転する。

無線周波パルスが終った後、磁気モーメントが磁界と整
合する状態に戻る。その時、励振パルスと同じ無線周波
数の小さな信号を放出する。放出された無線周波数を検
出して、作像する物体の薄いスライス内にある水素原子
核の密度を感知する。核磁気共鳴作（象の理論と実際の
完全な説明が、米国特許第４，４３１．９６８号及び同
第４．４４４，７６０号に記載されている。

水素の原子核は、関心のある大抵の軟らかい組織に大量
にあるので便利であり、それから発生した像は重要な価
値がある。磁気勾配と無線周波数を正しく遮択すれば、
例え°ばナトリウムの様な他の材料の磁気共鳴に基づく
作像も出来る。

原子核から放出される無線信号は非常に弱い。

受信する無線信号の信号対雑音比を改善する為に、作像
する物体の表面にごく接近して表面コイルを配置するの
が普通である。この表面コイルは、物体から隔てた送信
コイルを共に、受信する為にだけ用いてもよいし、或い
は無線周波信号の送信と受信の両方に用いてもよい。

表面コイルを用いて発生される像は、作像技術による人
為効果（ａｒｔｉｆａｃｔ）を含むことがあり、像の解
釈が複雑になる。１つの問題は、作像する物体を通る時
の無線周波数の系統的な減衰である。

作像する薄いスライスが、例えば人間の胴体に対して横
方向であると、胴体の一層深い部分では、無線周波信号
の減衰が強くなり、表面近くの部分から受信した信号の
強さが、胴体内部の更に深い所の信号よりも大幅に強く
なる。その結果、表面近くの部分は、受信した像では、
物体の中の一層深い所にある部分よりもかなり明るくな
って見える。この様に距離によって決定される像の変動
は、作像装置の産物であって、作像する解剖学的な部分
の結果ではない。

別の問題は、例えば脂肪性組織の様な成る種の組織によ
って生ずる一層明るい像と減衰の増加によって起る。脂
肪性組織の塊が作像するスライスの表面近くにあると、
その像の輝度は、例えば陰極線管の様な作像装置のダイ
ナミックレンジを越えることがある。１つの解決策は、
像の輝度を調節して、最も明るい区域を作像装置のダイ
ナミックレンジの範囲内に抑えることである。多くの場
合、スライスの中にある関心の持たれる物体は、像で明
るくなる脂＋１Ｊｉ区域の外側にある。全体的な像の輝
度を下げて、脂肪区域の輝度を表示装置のダイナミック
レンジ内に抑えようとする−と、それ程明るくない関心
のある区域は更に暗くなり、所望の情報が識別出来なく
なることがある。

輝度が高過ぎる区域は別の欠点がある。像のそれ程明る
くないが、情報を含む部分を解析しようとする者は、像
の中の明るすぎる区域がその近くに存在することによっ
て、目の疲れ及び不快感を経験することがある。

表面近くの脂肪区域の減衰が一層大きいことにより、無
線信号が脂肪区域を通過しなければならない様な、身体
の中の一層深い位置にある部分から受信した無線信号が
減少する。この為、像の輝度には、別の距離依存性を持
つ系統的な低下が起る。

従って、この発明の目的は、系統的な輝度変動に対して
象を補正する装置と方法を提供することである。

この発明の別の目的は、空間的な周波数が低い輝度変動
に対して像を補正する装置と方法を提供することである
。

この発明の別の目的は、像の平滑化した総体輝度マツプ
を作り、元の像には画素毎に、平滑化した総体輝度マツ
プの逆数を乗算して、作像装置の人為効果による輝度変
動を除外した補正像を作る装置と方法を提供することで
ある。

簡Ｉｌｌに云うと、この発明は像の空間周波数の低い輝
度マツプを作る２工程の平滑化方法を提供する。像に画
素毎に輝度マツプの逆数を乗算して、作像する物体によ
るものではなく、作像装置による像の輝度変動を補正す
る。平滑化の１番目の工程が、像にわたって四角の多重
画素の窓を走査し、像の総体輝度マツプを作る。窓の中
の中心の画素の輝度値は、窓の中の全ての画素の値の和
を付す。

別の実施例では、中心画素の輝度値として、窓の中の最
も明るい画素の値を付す。平滑化の２番目の工程が、一
層小さい四角の窓で総体輝度マツプを走査して、総体輝
度マツプの濃淡むらを平滑し、平滑輝度マツプを作る。

平滑により、中心画素には、窓の中の残りの画素の加重
和に等しい値か付せられる。窓の中の各々の画素の重み
の寄与を、その値を計算している画素からの距離に応じ
て変える。線形及び非線形加重を用いることが出来る。

元の像に画素毎に平滑輝度マツプの逆数を乗算して、補
正像を作る。補正像を正規化して、その値を表示装置の
ダイナミックレンジ内に戻す。

この発明の実施例では、夫々ディジタル値によって表わ
される輝度を持つ画素の行列の配列である様な像を補正
する装置を提供する。この装置は、その画素を取巻く第
１組の画素の内の少なくとも若干の画素の輝度の値を加
算することにより、像の中の各々の画素の総体輝度を推
定して、像の総体輝度マツプを発生ずる第１の手段と、
その値を推定している画素からの、それを取囲む画素の
距離に関係する市みを用いて、総体輝度マツプ中の画素
を取囲む第２組の画素の少なくとも若干の画素の輝度の
加重値を加算することにより、総体輝度マツプ中の各々
の画素の平滑輝度を推定して、平滑輝度マツプを発生す
る第２の手段と、像の中の各々の画素に平滑輝度マツプ
の対応する画素の逆数を乗算して、作像装置の人為効果
を除外した補正像を発生する手段と、作像装置のダイナ
ミックレンジに見合う値に、補正像の各々の画素のレベ
ルを復元する手段とを有する。

この発明の特徴として、夫々ディジタル値によって表わ
される輝度を持つ画素の行列の配列である様な像を補正
する方法を提供する。この方法は、その画素を取囲む第
１組の画素の内の少なくとも若干の画素の輝度の値を加
算することにより、像の中の各々の画素の総体輝度を推
定して像の総体輝度マツプを作り、その値を推定してい
る画素からの、それを取囲む画素の距離に関係する重み
を用いて、総体輝度マツプ中の画素を取囲む第２組の画
素の内の少なくとも若干の画素の輝度の加重値を加算す
ることにより、総体輝度マツプ中の各々の画素の平滑輝
度を推定して平滑輝度マツプを作り、像の中の各々の画
素に平滑輝度マツプ中の対応する画素の逆数を乗算して
、作像装置の人為効果をそれから除外した補正像を作り
、作像装置のダイナミックレンジに見合う値に、補正像
の各々の画素のレベルを復元する工程を含む。

この発明の上記並びにその他の目的、特徴及び利点は、
以下図面について説明する所から明らかになろう。図面
全体にわたり、同様な部分には同じ参照数字を用いてい
る。

好ましい実施例の詳しい説明 第１図について説明すると、核磁気共鳴作像装置によっ
て作像しようとする物体１０が、その周縁１４の近（に
脂肪の塊１２を含む。核磁気共鳴技術によって作られる
像のスライス１６が脂肪の塊１２を通っている。前に述
べた様に、像のスライス１６からの無線周波数の放出が
、周縁１４の外部から送信された適当な無線周波帯のパ
ルスによって励振される。像のスライス１６の内、脂肪
の塊１２の中にある部分は、脂肪の塊１２の他の部分よ
りも、一層大きな無線周波出力を発生することが出来る
。像のスライス１６の中心を励振する無線周波数は、周
縁１４の近くにある像のスライス１６の縁を励振する無
線周波数よりも、物体１０の中を一層長い距離通らなけ
ればならないし、同様に放出された無線周波数も物体１
ｏの中を一層長い距離通らなければならない。この為、
像のスライス１６の中心から受信した無線周波数は、周
縁１４の近くから受信したものよりも弱いのが普通であ
る。この結果、周縁１４の近辺に較べてξΦ心は系統的
に像が一層暗くなる。

脂肪の塊１２は、不釣合いに一層明るい像を作る他に、
物体１０の他の部分に較べて、物体１゜に入り且つそれ
から出て行く無線周波数の減衰を強める。成る意味で、
脂肪の塊１２は、像のスライス１６の近くにある部分を
、励振無線周波数の受信及び外部受信機への放出された
無線周波数の伝達の両方に対して遮蔽するものとみなす
ことが出来る。

脂肪の塊１２が不釣合いに明るい部分として作像され、
像のスライス１６の残りの部分は比較的暗い。像のスラ
イス１６の中心に一層近い部分は、周縁１４に一層近い
部分よりも暗い。更に、像のスライス１６の内、脂肪の
塊１２を通過する無線周波数を受取り且つやはり脂肪の
塊１２を通過する無線周波数を放出する部分は、作像し
た時、像のスライス１６の内、脂肪の塊１２から離れた
端にある部分よりも一層暗い。

核磁気共鳴作像で普通行なわれる様に、像がディジタル
・ファイルとして貯蔵され、像は、各々の画素の輝度が
２進数で表わされる画素の配列として表わされている。

この発明は、上に述べた様な原因で起る像の輝度の系統
的な変動に対処する為、見る者にとって関心が一層大き
いかも知れない様な、像のスライス１６の他の部分を洗
い流したり或いは過度に暗（せずに、脂肪の塊１２の様
な身体の最大輝度を制限する形で、画素を表わす２進数
を操作する。

第２図には、像収集装置２０で発生された像を補正する
ディジタル形像補正装置１８が示されている。像収集装
置２０は例えば核磁気共鳴作像装置の様な任意の普通の
形式であってよい。発生されたま＼の像の要素は一般的
にアナログである。

アナログの像要素又は画素をディジタル化して、例えば
４つの２進デイジツト（ビット）から成る対応するディ
ジタル数により、像の中の各々の画素を表わす。

この発明を補正する主な種類の人為効果は、像にわたっ
て滑らかに変化する人為効果である。こういう人為効果
は前に述べた距離依存性及び遮蔽依存性を持つ効果によ
って起る。第１近似として、空間周波数が低い、像の輝
度の変動は、作像する身体の特性ではなく、この様な人
為効果によるものであると仮定することが出来る。低い
空間周波数とは、像が変化する速さを指す。例えば、輝
度を持つ成る点が、高い空間周波数を表わすが、像の一
方の縁から他方の縁まで像の輝度が一定の勾配であるこ
とは、空間周波数が低いことを表わす。

ディジタル形像補正装置１８が、一方の縁から他方の縁
までの像の輝度の緩やかな変化の表示を取出す為に、像
収集装置２０によって発生された像の分解能を低下させ
る。

像の分解能の低下は２工程に分けて行なわれる。

工程１では、総体輝度マツプ装置２２が、その輝度を推
定しようとする画素を中心とする窓の中の予定数の画素
の輝度の和から、各々の画素の輝度の推定値を発生する
。更に好ましい実施例では、各々の画素の推定値は、窓
の中にある一番明るい画素の値である。工程１により、
元の像の総体輝度マツプが発生される。総体輝度マツプ
は、明るいスポットが窓に入ったり、出たりする時の、
輝度の急激な上向き及び下向きの階段を持つむらのある
表示である。

工程２で、平滑装置２４が、その輝度を推定しようとす
る画素を中心とする窓の中の総体輝度マツプの画素の加
重和を求めることにより、各々の画素の値を推定するこ
とによって、総体輝度マツプの輝度の階段を平滑する。

ディジタル形像補正装置１８の説明を続ける前に、第３
図及び第４図について説明する。元の像２６がＸ−Ｙ配
列の画素の行及び列で構成される。

第１の位置にある窓２８を実線で示しであるが、これは
その輝度を推定しようとする中心画素３０を中心として
いる。図示の窓２８がＸ方向にそれを取巻く１３個の画
素３２を持ち、Ｙ方向にも１３個の画素を持っており、
こうして中心画素３０の上下左右に６個の画素３２が取
巻く様にしている。

総体輝度マツプ装置２２（第２図）が、窓２８の中にあ
る、それを取巻く全ての画素３２の輝度値を加算し、そ
の結果を中心画素３０に指定することにより、中心画素
３０の輝度値を推定する。

推定された輝度は、この推定に入る画素の総数で除すこ
とによって、倍率を定めることが出来るが、後で説明す
るこれより後の正規化動作によって、この工程が省略さ
れる。

中心画素３０の輝度を推定した後、窓２８を行に沿って
画素１個だけ歩進させ、第３図に破線で示す位置に持っ
て来る。この動作を繰返して、新しい中心画素３０の輝
度を推定する。

この発明では、密度低下方法（ｒｅｄｕｃｅｄ−ｄｅｎ
ｓ　ｉｔｙ　ｔｅｃｈｎｉｑｕｅ　）を用いて輝度の推
定値を計算することにより、性能を犠牲にせずに計算時
間を大幅に短縮することが出来ることが判った。即ち、
窓２８の中にある、周囲の全ての画素３２の値を加算す
る代りに、好ましい実施例では、取巻く全部の画素３２
より少ない画素を推定に入れる。例えば、この発明では
、中心画素３０の輝度を推定するのに、Ｘ及びＹ方向に
取巻く４番１］毎の画素３２を用いても、性能が目立っ
て低下しない。これによって、計算速度は１６倍になる
。

これまで述べた推定方法は、中心画素３０の推定値を計
算するのに、取巻く画素３２の重みをつけない値を用い
ている。これはむらのある総体輝度マツプになる。窓２
８の右側の緑にある１つの取巻く画素３２が極めて明る
いと仮定し、元の像２６内の他の全てが一様な暗さであ
ると仮定すれば、このむらが理解出来よう。明るい取巻
く画素３２を包み込む様に窓２８を動かす１工程で、全
ての取巻く画素３２の和は階段形の増加をする。

窓２８の１３回の歩進の後、明るい取巻く画素３２が窓
２８の外へ出て、その結果Ｈｋ定値の［１４度に下向き
の階段が生ずる。

窓２８にある画素の数は、例えば元の像２６にある画素
の総数、補正すべき像の劣化の激しさに応じて、中心画
素３０の推定値を計算するのに利用し１）る時間との釣
合いをとって変えることが出来る。更に、計算時間を短
縮する為に選ばれる取巻く画素３２の輝度も同様な因子
に関係する。好ましい実施例では、窓２８は６５Ｘ６５
個の画素であって、４番目毎の画素が取巻く画素３２に
含まれることが好ましい。

更に好ましい実施例では、中心画素３２は、窓２８の中
の最も明るい取巻く画素３２の輝度を付す。これによっ
てこの結果得られる総体輝度マツプのむらが増大するが
、この後の平滑動作でこの様に増加したむらが除かれる
。

平滑装置２４が総体輝度マツプの倍率を加えた加算を行
う。第４図について説明を続けると、総体輝度マツプの
上の窓２８が、中心画素３０と、縮めた配列に配置され
た取巻く画素３２とを含む。

総体輝度マツプ装置２２によって発生される総体輝度マ
ツプから説明を続けると、平滑装置２４が元の像２６の
総体輝度マツプにわたって一度に１画素ずつ、窓２８を
動かし、各々の位置に於ける水平及び垂直の値を加算し
て、中心画素３０の輝度を推定する。各々の取巻く画素
３２の輝度値には、和を求める前に、中心画素３０から
のその距離に応じた蚤みを加える。倍率作用は、取巻く
画素３２の水平及び中心線からの各々の画素の距離によ
って、取巻く画素３２に与えられる輝度値に１１１１位
の減少が生ずる様にする線形傾斜関数にすることが出来
る。便利な１つの倍率方法は、窓２８の周縁に沿った各
々の取巻く画素３０の輝度を６（窓２８の垂直及び／又
は水平中心線からの距離）で除し、残りの取巻く画素３
２の中にある各々の取巻く画素３２の輝度を２（窓２８
の垂直及び／又は水平中心線からのその距離）で除す。

更に曳雑な処理が必要であるが、各々の取巻く画素３２
の輝度を、中心画素３０からの対角線方向の距離を表わ
す数で除してから、和を求めてもよい。線形傾斜関数以
外の関数を用いて、加算過程に入る前の取巻く画素３２
の輝度の倍率を定めることが出来る。例えば、余弦関数
を用いることが出来る。この場合、各々の取巻く画素３
２の輝度には、中心画素３０からのその距離を表わす余
弦関数を乗算する。

平滑装置２４で使われる窓の寸法は、特定の用途に応じ
て選ぶことも出来る。この寸法は、総体輝度マツプ装置
２２の窓２８の寸法より小さいことが好ましい。総体輝
度マツプ２２で１辺６５個の画素の窓２８を用いる時、
平滑装置２４では、１辺１３個の画素の窓２８を用いる
のが好ましい。

総体輝度マツプ装置２２及び平滑装置２４に於ける中心
画素３０の輝度の推定では、特定の和で選択されなかっ
た取巻く画素３２が省略されているが、最終的な結果に
は全ての画素が表示されることに注意されたい。即ち、
窓２８が一度に画素１個だけ水手に歩進する時、和に入
る画素を持つ行のことごとくの画素が、水平走査の端に
表示される。この後の走査の間、各々の水平走査の前に
、窓２８が画素の１行だけ下向きに歩進する。その結果
、１回の水平走査の間に省略された画素が、この後の水
平走査の間に計算に入る。

平滑装置２４の出力は、それから像の全ての細部か脱落
した様な、元の像の空間周波数の低いマツプである。こ
のマツプは、距離、減衰及び現実の解ｚす学的な特徴を
表わさないその他の効果によって生じた像の人為効果の
良好な表示である。

第２図に戻って説明すると、平滑装置２４からの平滑さ
れた像が像補正装置３４に印加され、これが像収集装置
２０からめ元の像にある各々の画素に、平滑装置２４か
らの平滑輝度マツプ内の対応する画素の逆数を乗算する
。こうして、作像装置の人為効果が像から除かれる。

成る作像装置では、表示の為に、元の像の各々の画素に
定数が加算される。この定数は、平滑輝度マツプの対応
する画素で除算する前に、像補正装置３４で各々の画素
から減算することが好ましい。

その値は、１より大きいことも小さいこともあるが、そ
ういう値でその画素を除算した補正像は、表示装置のダ
イナミックレンジを最大限に活用する様な輝度値をもは
や持たないことがある。正規化装置３６が表示装置のダ
イナミックレンジを活用する様な範囲へ、像を復元する
。

正規化装置３６が像補正装置３４からの補正像の最も明
るい画素の輝度を決定する。その後、それが、最も明る
い画素の輝度を表示装置のダイナミックレンジの最大値
に等しくなる様に調節する補正係数を計算する。最後に
、それが補正像内のことごとくの画素にこの補正係数を
乗算する。この結果得られた像は、初めは距離効果及び
減衰効果の影響を受けていた区域で輝度が増加し、初め
は放出の強い材料が存在したこと又は表面近くの場所の
為に明るすぎた区域では輝度が低下する。

像補正装置３４に於ける補正の前に、各々の画素から表
示定数を減算する場合、正規化装置３６を出る前に、各
々の画素に対応する定数を加算する。

成る用途では、補正像の成る画素の暗さが、作像装置の
雑音閾値にあるか或いはそれより低いことがある。この
為、最も明るい画素に対して予定の比より低い輝度を持
つ各々の画素に対口、低レベル閾値を適用することが出
来る。例えば、最も明るい画素の輝度の１／６未満の輝
度（直では雑音しか表示されないことが判った場合、低
レベル閾値は、一層低い輝度を持つ全ての画素に対し、
最も明るい画素の１／６に等しい値を与えることが出来
る。輝度が低い画素を表示しようとすることは、雑音を
表示させることに他ならないから、この様な閾値を使う
方法によっても、同等有益な情報が失われることはない
。実際の低レベル閾値は作像装置によって左右され、最
も明るい画素の値の１７６の値が限界であると考えては
ならない。

この発明の重要な１つの価値は、１１１７正を必要とす
る項目が事前に判っていなくても、表面コイルの欠陥、
距離、減衰及び明るい物体による像の人為効果を補正す
ることが出来ることである。

像収集装置２０によって収集された元の像に対し、ディ
ジタル形像補正装置１８を作用させ、永久的に変更する
ことが出来る。この代りに、破線３８で示す様に、元の
像は変更なしに貯蔵又は表示の為に伝達し、元のデータ
を乱さずに、正規化装置３６によって補正像を発生する
ことが出来る。

元の像を利用出来ることは、この発明にとっては関係の
ない他の種類の像の解析に価値があることがある。

この発明の好ましい実施例を図面について説明したが、
この発明がこの実施例そのものに制限されず、当業者で
あれば、特許請求の範囲によって定められたこの発明の
範囲内で種々の変更を加えることが出来ることは云うま
でもない。

【図面の簡単な説明】

第１図はこの発明を説明する為の核磁気共鳴装置によっ
て作像された胴体のスライスを示す断面図、 第２図はこの発明の１実施例の像補正装置及び方法を示
す機械的な簡略ブロック図、 第３図は表示装置の正面図で、平滑窓の２つの位置を示
す。 第４図は第３図の平滑窓の拡大図で、関連する画素を示
している。 主な符号の説明 ２２：総体輝度マツプ装置 ２４：平滑装置 ３４：像補正装置 ３６：正規化装置

Claims (1)

1. 【特許請求の範囲】 １）夫々がディジタル値によって表わされる輝度を持つ
   画素の行列の配列である様な像を補正する装置に於て、 該像内の各々の画素の総体輝度を、前記画素を取囲む第
   １組の画素の内の少なくとも１つの画素から推定して、
   前記像の総体輝度マップを作る第１の手段と、 該総体輝度マップ内の各々の画素の平滑輝度を推定して
   平滑輝度マップを作る第２の手段と、前記像の各々の画
   素に、前記平滑輝度マップ内の対応する画素の逆数を乗
   算して、作像装置の人為効果を除いた補正像を発生する
   手段と、 作像装置のダイナミックレンジに見合う値に前記補正像
   の各々の画素の輝度レベルを復元する手段とを有する像
   を補正する装置。 ２）特許請求の範囲１）に記載した装置に於て、重みが
   、その平滑値を推定している画素からの、第２組の画素
   の内の各々の画素の距離に線形関係を持つ装置。 ３）特許請求の範囲１）に記載した装置に於て、重みが
   、その平滑値を評価している画素からの、第２組の画素
   内の各々の画素の距離に非線形関係を持つ装置。 ４）特許請求の範囲１）に記載した装置に於て、第１組
   の画素の内の少なくとも若干の画素は、その総体輝度を
   推定する画素を取囲む画素の全部より少ない装置。 ５）特許請求の範囲４）に記載した装置に於て、前記画
   素の全部より少ない画素が、画素の第２の規則的な配列
   をそれから除外した画素の第１の規則的な配列を含んで
   いる装置。 ６）特許請求の範囲４）に記載した装置に於て、第１組
   の画素の内の少なくとも若干の画素が、垂直方向及び水
   平方向の窓の中の４番目毎の画素を含んでいる装置。 ７）特許請求の範囲１）に記載した装置に於て、前記復
   元する手段が、画素の輝度の最低値に閾値を設ける手段
   を含み、こうして雑音を抑圧した装置。 ８）特許請求の範囲１）に記載した装置に於て、前記復
   元する手段が、前記補正像内の最大輝度を持つ画素に補
   正像を正規化することを含む装置。 ９）特許請求の範囲１）に記載した装置に於て、各々の
   画素の総体輝度を推定する手段が、該画素から予定数の
   画素以内にある、それを取囲む最も明るい画素の輝度に
   等しい輝度値を各々の画素に付すことを含む装置。 １０）特許請求の範囲１）に記載した装置に於て、前記
   第２の手段が、その値を推定している画素から取囲む画
   素までの距離に関係する重みを用いて、前記総体輝度マ
   ップ内で前記画素を取囲む第２組の画素の内の少なくと
   も若干の画素の輝度の加重値を加算する手段を含む装置
   。 １１）夫々ディジタル値によって表わされる輝度を持つ
   画素の行列の配列である様な画素を補正する装置に於て
   、 前記像内の各々の画素の総体輝度を、該画素を取囲む第
   １組の画素の内の少なくとも１つの輝度を用いて推定し
   て、前記像の総体輝度マップを作り、 該総体輝度マップ内の各々の画素の平滑輝度を推定して
   平滑輝度マップを作り、 前記像内の各々の画素に、前記平滑輝度マップ内の対応
   する画素の逆数を乗算して、作像装置の人為効果をそれ
   から除外した補正像を作り、作像装置のダイナミックレ
   ンジに見合う輝度値に、補正像の各々の画素のレベルを
   復元する工程を含む方法。 １２）特許請求の範囲１１）に記載した方法に於て、各
   々の画素の総体輝度を推定する工程が、該画素から予定
   数の画素以内にある、それを取囲む最も明るい画素の輝
   度に等しい輝度値を各々の画素に付すことを含む方法。 １３）特許請求の範囲１１）に記載した方法に於て、平
   滑輝度を推定する工程が、その値を推定する画素からそ
   れを取囲む画素までの距離に関係する重みを用いて、前
   記総体輝度マップ内の、前記画素を取囲む第２組の画素
   の内の少なくとも若干の画素の輝度の加重値を加算する
   ことを含む方法。