JP3462957B2

JP3462957B2 - 画像処理方法

Info

Publication number: JP3462957B2
Application number: JP14892696A
Authority: JP
Inventors: 順一田口; 邦彦木戸; 耕一佐野
Original assignee: Hitachi Ltd; Hitachi Medical Corp
Current assignee: Hitachi Ltd; Hitachi Healthcare Manufacturing Ltd
Priority date: 1996-06-11
Filing date: 1996-06-11
Publication date: 2003-11-05
Anticipated expiration: 2016-06-11
Also published as: JPH09330399A

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】画像処理を行える機能を有す
処理装置において、画像のノイズ低減を目的とした平滑
化処理方法を提供する。処理対象となる画像は、例え
ば、磁気共鳴診断装置によって撮影された画像、Ｘ線Ｃ
Ｔ装置によって撮影された画像、人工衛星から伝送され
る画像、写真をスキャナで読んだ画像、コンピュータグ
ラフィックを使って作成した画像など多く存在する。各
々の画像を取得処理する装置で本発明を実施することも
できるし、画像処理を行う専用ハードとして本発明を実
施することもできる。

【０００２】

【従来の技術】

（１）特願平５−１４０３３４木戸邦彦、他：ノイズ
低減フィルター（２）特願平６−３０５８８８木戸邦彦、他：画像処
理方式上記（１）は、画像の各点毎に、最小変化方向を見つ
け、その最小変化方向の所定のサンプリング点について
加算平均をした値を出力値とする。上記（２）は、画像
の各点毎に、最小変化方向を見つけ、その最小変化方向
上の所定のサンプリング点について非線形平滑化をした
値を出力値とする。非線形平滑化の具体的方法として、
以下の数１〜数２に基づいた非線形平滑化処理を行う。

【０００３】

【数１】Ｈ２＝ａ・Ｈ１ + (1-ａ)・Ｇ(i,j)

【０００４】

【数２】ａ＝ ρ0**2／(ρ1**2＋(α・ρ0)**2) ただし、Ｇ(i,j)は画像のｉ行ｊ列の点の輝度、Ｈ１は
上記最小変化方向の所定のサンプリング点について加算
平均した値で、上記従来技術（１）の出力値であり、Ｈ
２は上記従来技術（２）の出力値であり、ρ1は上記最
小変化方向の所定のサンプリング点について偏差を求め
た値であり、ρ0は画像からノイズの程度を決める値
で、例えば画像の最大値より５％以下の点の平均値であ
り、αは平滑化の程度を決める所定の定数である。

【０００５】上記従来技術（１）では画像の構造の入り
組んだエッジ部分で画像構造がボケる欠点があり、上記
従来技術（２）は偏差の大きい部位の平滑化を抑えるこ
とでその欠点を克服している。

【０００６】

【発明が解決しようとする課題】上記従来技術（２）で
は上記最小方向の平滑化を行う所定のサンプリング点数
をあまり多くできないという課題がある。上記従来技術
（２）で、最小方向の平滑化を行うサンプリング点数
は、５点から７点までが適切であり、それ以上多くなる
と平滑化の程度を決める定数αの設定が難しくなり、画
像のボケが目立つ場合と平滑化の効果が得られない場合
に分かれやすくなって、処理効果が安定して得られなく
なる。

【０００７】以上の理由から、上記従来技術（２）で
は、最小方向の平滑化を行うサンプリング点数は、５点
か７点程度を選ぶことになるが、磁気共鳴診断装置で撮
影した腹部のＳ／Ｎが悪い画像に本処理を適用すると、
画像に不自然な線構造が現れ人工感の強い画像になると
いう問題点があった。最小方向の平滑化を行うサンプリ
ング点数を９点にして定数αを適切な値に設定すると画
像の不自然な構造が消えたが、前述のように定数αの設
定が難しくなるという問題点があった。

【０００８】本発明の目的は、上記従来技術（２）の問
題点を考慮し、最小方向の平滑化を行うサンプリング点
数を９点以上にしても、平滑化の程度を決めるパラメー
タ値にある程度幅を持って平滑化の効果が現れる非線形
平滑化方法を提供することにある。

【０００９】また、本発明は、上記平滑化の程度を決め
るパラメータ値にある程度幅を持たせるだけでなく、そ
れと同時に、なるべくエッジ部のボケを抑えた平滑化を
行うことのできる方法を提供することを目的とする。

【００１０】

【課題を解決するための手段】上記問題点を解決するた
め、上記最小方向の１次元非線形平滑化を以下の数３〜
数５に基づいて行う。

【００１１】

【数３】Ｗ(i,j,q,m) = ｗ(i,j,q,m)／ｗｔ

【００１２】

【数４】ｗ(i,j,q,m) = １／（１ + ((Ｉ(i,j,q,m)-Ｇ
(i,j))／(α・ρ))**k ）

【００１３】

【数５】ｗｔ = Σ ｗ(i,j,q,m) ただし、ｉ，ｊは画像の着目点の位置を示し、ｉ行ｊ列
の意味、ｑは画像のｉ行ｊ列点における最小変化方向を
所定の数値で表した値、ｍは画像のｉ行ｊ列点を中心と
してｑ方向の１次元直線上の所定のサンプリング点を順
序づけた順番、Ｇ(i,j)は画像のｉ行ｊ列点における輝
度を表し、Ｉ(i,j,q,m)は画像のｉ行ｊ列点における最
小変化方向であるｑ方向にｍサンプリング点行った点の
輝度を表し、Ｗ(i,j,q,m)は平滑化でＩ(i,j,q,m)を加算
する際のウエイトを示し、ｗ(i,j,q,m)とｗｔはウエイ
トＷ(i,j,q,m)を計算する際の途中項で、ｋは所定の定
数で、**kは、ｋ乗を意味し、ρは画像からノイズの程
度を決める値で、例えば画像の最大値より５％以下の点
の平均値であり、αは平滑化の程度を決める所定の定数
であり、Σはｍについて所定の点数全ての和を取るもの
とする。

【００１４】また、さらに画像のエッジ部でのボケを少
なくした自然な平滑化を行うため、以下の数６〜数１０
に基づく平滑化を行う。

【００１５】

【数６】Ｗ(i,j,q,m) = ｗ(i,j,q,m)／ｗｔ

【００１６】

【数７】ｍ≧１の時、ｗ(i,j,q,m) = ｗ(i,j,q,m-1)／（１+((Ｉ(i,j,q,m)-Ｇ
(i,j))/(α・ρ))**k）

【００１７】

【数８】ｍ＝０の時、ｗ(i,j,q,0) = １

【００１８】

【数９】ｍ≦−１の時、ｗ(i,j,q,m) = ｗ(i,j,q,m+1)／（１+((Ｉ(i,j,q,m)-Ｇ
(i,j))/(α・ρ))**k）

【００１９】

【数１０】ｗｔ = Σ ｗ(i,j,q,m) 上記従来技術（２）の出力結果は、平滑化する１次元方
向の全体の平均値Ｈ１と原画像の輝度との中間の値であ
り、どちらに近い値とするかを平滑化の程度を決めるパ
ラメータαを含む式に依存して決めるが、１次元平滑化
をするサンプリング点数を多く取ると平均値Ｈ１は画像
本来の値からかけ離れて来るので、その分、平滑化の程
度を決めるパラメータαの設定が難しくなる。

【００２０】上記数３〜数５に基づいた平滑化では、着
目点上の画像輝度と平滑化をする１次元データ上の点の
輝度の差に応じて平滑化のウエイトを決めるので、平滑
化の程度を決める定数αをある程度変化させても、変化
の大きい部位の平滑化は抑えられるので、画像本来の値
を保つように平滑化され、極端なパラメータαについて
の依存性が現れない。

【００２１】上記数６〜数１０も、同様に変化の大きい
部位の平滑化を抑えるので、極端なパラメータαについ
ての依存性が現れない。また、平滑化する１次元データ
上の点の位置が方向を検出した着目点の位置から離れる
につれてウエイトの値が小さくなる式であるので、一度
段差が現れた場合、別の構造に移ったと判断して平滑化
を抑えるような働きがあり、画像エッジ部位のボケが少
なくなる。

【００２２】

【発明の実施の形態】本発明は、画像処理機能を有す装
置において画像処理を行う方法の発明であり、具体的に
は下記の方法を実施するプログラムまたはハードを作成
することにより実現できる。以下、手順概要を述べた後
に具体的な手順と方法について述べる。

【００２３】（１）手順概要図１を用い手順の概要を示す。画像上の各点においてそ
れぞれ以下の処理を行う。今、仮にある画像上の１点に
ついて着目したとする。

【００２４】画像上の着目点における８方向１０１上に
所定のサンプリング点を取り、各方向の変化量を計算す
る。その変化量の中から最小値となる方向を求め、最小
変化方向１０２を求める。最小変化方向１０２における
サンプリング点１０３の輝度データを１次元の被平滑化
データとする。この１次元データの輝度の差に応じてウ
エイト付けした平滑化１０４を行い処理値とする。

【００２５】以上の操作を各点について行うが、画像の
端の点は、８方向１０１について中心部と同様の所定の
サンプリング点を取れないため、特殊処理を行う必要が
ある。特殊処理として、画像の両端がつながっていると
して中心部と同様の所定のサンプリング点を取る方法
や、端の点は処理をせずに原画や０を出力する方法、端
のサンプリング点を少なくして処理する方法、原画のま
ま出力する方法など、端の部分の特殊処理方法は各種あ
る。

【００２６】端の特殊処理を除いた画像中心部に対する
処理手順の詳しい流れ、画像上の着目点における８方向
１０１のサンプリング点の取り方、各方向の変化量の計
算の仕方、輝度の差に応じたウエイト付け平滑化１０４
の処理方法については以下各々節を設けて具体的に述べ
る。

【００２７】（２）手順図２と図３に処理の流れの例を示す。各々、端の特殊処
理の必要のない画像中心部分の処理手順について示し
た。図２は各画素毎に各々最終処理結果を求めていく流
れ、図３は計算の途中経過のマップを作り、マップを参
照して処理をする流れの例である。各々図中の番号をス
テップ番号として、以下説明する。

【００２８】はじめに図２の流れについて説明する。

【００２９】ステップ２０１：画像の各点をｉ行ｊ列
で記述し、点特殊処理の必要のない画像領域を以下の数
１１と数１２の領域とする。

【００３０】

【数１１】ｉｓｔａｒｔ≦ｉ≦ｉｅｎｄ

【００３１】

【数１２】ｊｓｔａｒｔ≦ｊ≦ｊｅｎｄただし、ｉｓｔａｒｔ、ｉ、ｉｅｎｄ、ｊｓｔａｒｔ、
ｊ、ｊｅｎｄは整数。

【００３２】この時、初めに、ｉ＝ｉｓｔａｒｔ、ｊ＝
ｊｓｔａｒｔ、と初期値の設定をする。

【００３３】ステップ２０２：着目点（ｉ，ｊ）を中
心として、８方向１０１の輝度変化量を求める。８方向
１０１の取り方や輝度変化量の計算の仕方は後に述べ
る。

【００３４】ステップ２０３：８方向１０１の輝度変
化量のうち、最も小さい値を持つ方向（最小変化方向１
０２）を求める。最小値を持つ方向が２つ以上の場合
は、どちらを選ぶか各種考えられる。あらかじめ所定の
順番で最初の方向を選ぶこともできるし、後の方向を選
ぶこともできるし、ランダムにどちらかを選ぶこともで
きるし、各々選んで、後の処理はその両者について各処
理を行い、両者の処理結果の平均値を処理値とすること
もできる。

【００３５】ステップ２０４：最小変化方向１０２に
ついて、その方向のサンプリング点１０３上の１次元デ
ータを参照して輝度の差に応じたウエイトづけ平滑化１
０４を行う。最小変化方向１０２上のサンプリング点１
０３の取り方は、ステップ１０２で計算した８方向１０
１のサンプリング点と同一のサンプリング点から選ぶこ
ともできるし、それと異なる点を選び、サンプリング点
数も増やすことや減らすこともできる。

【００３６】ステップ２０５：ｉがｉｅｎｄより大きい
かどうか判定し、ｉｅｎｄより大きければ次のステップ
２０６に進み、そうでなければｉを１増やしてステップ
２０２に戻る。

【００３７】ステップ２０６：ｊがｊｅｎｄより大きい
かどうか判定し、ｊｅｎｄより大きければ処理を終わ
り、そうでなければｊを１増やしてステップ２０２に戻
る。

【００３８】次に図３の流れについて説明する。

【００３９】ステップ３０１：数１１と数１２を満た
す各点においてそれぞれ８方向１０１の輝度変化量を求
め、それぞれの点における最小方向１０２を求め、最小
方向マップを作る。方向は、８方向１０１についてそれ
ぞれ順番を決め、例えば０から７の整数値で各方向を表
すことにする。最小方向マップは、各点における最小方
向の数値を記録したマップである。

【００４０】ステップ３０２：最小方向マップを参照
して、各点毎に最小変化方向について所定のサンプリン
グ点１０３の１次元データについて輝度の差に応じたウ
エイト付け平滑化１０４を行う。

【００４１】図２の手順であっても、図３の手順であっ
ても、どちらの手順でも実施できる。その他にも、例え
ば、ステップ３０１の手順を分離して、８方向１０１の
変化量を求めた変化量ベクトルマップを作り、変化量ベ
クトルマップを参照して最小方向マップを作ることもで
きる。なお、変化量ベクトルマップは、画像上の各点で
８方向の変化量を持ったべくトル値を持つので、ベクト
ルマップと呼んだ。

【００４２】（３）８各方向１０１上のサンプリング点
の取り方図４に８方向１０１上のサンプリング点の取り方の一例
を示す。図４は、各方向に９点のサンプリング点を取っ
た場合であり、黒丸で示した点が画像上の画素の上のサ
ンプリング点で、白丸で示した点が画素の上になく補間
を必要とするサンプリング点である。補間の仕方は各種
あり、左右上下のうちの近い２点の平均値とする方法や
近傍の４点からキュービックコンボリューションをする
方法や、最も近い画素点のうちの特定の画素点の値とす
る方法などがある。図５は、図４の白丸で示した点の値
としてその点から最も近い２点のうち、中心に近い方の
点の値を参照する場合について図示したものである。図
４の補間点を最近傍の点から選んだと考えてもよいし、
もともと８方向を図５のように非直線的に選んだと考え
ても良い。

【００４３】８方向１０１のサンプリング点の選び方
は、図４や図５のように各方向に９点選ぶ方法だけでな
く、１１点や１３点などもっと中心から遠い多くの点を
選ぶこともできるし、７点や５点など、中心から近い点
だけの少ないデータを選ぶこともできる。

【００４４】（４）各方向の変化量の計算の仕方各方向の変化量の計算の仕方の具体例を以下の数１３に
示す。

【００４５】

【数１３】Ｓ(i,j,p) = Σ abs(Ｉ(i,j,p,m)-Ｇ(i,j)) ただし、ｉ，ｊは画像の着目点の位置を示し、ｉ行ｊ列
の意味、ｐは画像のｉ行ｊ列点の８方向のうち、着目方
向を所定の数値で表した値、ｍは画像のｉ行ｊ列点を中
心としてｐ方向の１次元直線上の所定のサンプリング点
を順序づけた順番、Ｇ(i,j)は画像のｉ行ｊ列点におけ
る輝度を表し、Ｉ(i,j,p,m)は画像のｉ行ｊ列点におけ
る着目方向であるｐ方向にｍサンプリング点行った点の
輝度を表し、Σはｍについて所定の点数全ての和を示
し、abs(）は、括弧内の値の絶対値を取る関数であり、
Ｓ(i,j,p)は画像のｉ行ｊ列点における着目方向である
ｐ方向の変化量である。

【００４６】以上の数１３の他にも変化量を算出する定
義は各種ある。例えば、以下の数１４を用いることもで
きる。

【００４７】

【数１４】Ｓ(i,j,p) = Σ (Ｉ(i,j,p,m)-Ｇ(i,j))**２ただし、各記号は、数１３と同様であっり、**２は２乗
を示す。

【００４８】（５）輝度の差に応じたウエイト付け平滑
化１０４の処理方法最小変化方向の１次元サンプリング点について行う平滑
化の処理方法は、上記手段の項記載の数３〜数５に記載
の処理を行うこともできるし、数６〜数１０記載の処理
を行うこともできる。その他にも各種方法があり、例え
ば、数６〜数１０の処理のうち、数７を下記の数１５に
変え、数９を下記の数１６に変えて実行することもでき
る。

【００４９】

【数１５】ｍ≧１の時、ｗ(i,j,q,m) = ｗ(i,j,q,m-1)／（１+exp((Ｉ(i,j,q,m)
-Ｇ(i,j))/(α・ρ))）

【００５０】

【数１６】ｍ≦−１の時、ｗ(i,j,q,m) = ｗ(i,j,q,m+1)／（１+exp((Ｉ(i,j,q,m)
-Ｇ(i,j))/(α・ρ))）なお、数４、数７、数９、数１５、数１６には、ρがあ
るが、ρは画像輝度の変化程度を表す統計量であり、画
像からノイズの程度を求めたものであればどのようなも
のも採用できる。例えばρとして画像の最大値より５％
以下の点の平均値を採用することもできるし、画像の４
隅の偏差のうち小さい方を採用することもできるし、画
像の隣合う点の引き算をしたデータのヒストグラムを作
り、ヒストグラムを小さい方から足していって全体の６
０％となる値までの平均値とすることなどもできる。

【００５１】以上本発明の実施例を示したが、上記実施
例の方向は８方向定義して最小方向を検出したが、４方
向や１６方向など方向の定義の仕方は各種ある。その他
にも、例えば、上記実施例では画像の最小変化方向のサ
ンプリング点についてウエイト付け平滑化１０４を行っ
ていたが、各方向の変化量を演算してボケ操作をして大
域的な方向をつかんだ画像の最小変化方向について１次
元のサンプリング点を選び上記ウエイト付け平滑化１０
４の操作を行うこともできるし、最小変化方向ではか
く、最大変化方向の直行方向を選ぶこともできるし、最
大変化方向についてさえも、上記ウエイト付け平滑化１
０４の操作を行うこともできる。以上の様に本発明のウ
エイト付け平滑化１０４は、１次元データであればどの
ようなデータであっても適用することは可能である。例
えば、音声データとか、プラント制御のためのプラント
モニタから送られてくる時系列データなどについても本
発明の上記ウエイト付け平滑化１０４の操作を行うこと
が可能である。１次元データに本発明の上記ウエイト付
け平滑化１０４を行うと、見た目に自然な平滑化ができ
るという効果がある。

【００５２】

【発明の効果】本発明によれば、画像のエッジ部のボケ
を抑え、見た目に自然な平滑化をした画像が得られる効
果がある。

【図面の簡単な説明】

【図１】本発明の概略を示した図である。

【図２】本発明の１手順を示した図である。

【図３】本発明の別の手順を示した図である。

【図４】８方向とサンプリング点を示した図である。

【図５】最近傍法による８方向とサンプリング点を示し
た図である。

【符号の説明】

１０１…画像上の着目点における８方向、１０２
…最小変化方向、１０３…最小変化方向におけるサンプ
リング点、１０４…輝度の差に応じたウエイト付け平滑
化。

───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.⁷ 識別記号ＦＩＨ０４Ｎ 1/409 Ｈ０４Ｎ 1/40 １０１Ｃ (72)発明者佐野耕一神奈川県川崎市麻生区王禅寺1099番地株式会社日立製作所システム開発研究所内 (56)参考文献特開平７−262368（ＪＰ，Ａ) 特開平７−288768（ＪＰ，Ａ) 特開平６−282641（ＪＰ，Ａ) (58)調査した分野(Int.Cl.⁷，ＤＢ名) G06T 1/00 - 7/60 H04N 1/40 - 1/409 A61B 5/055 A61B 6/00 - 6/14

Claims

(57)【特許請求の範囲】

【請求項１】（１）画像の各点においてそれぞれ最小変
化方向を検出し、（２）画像各点でそれぞれ上記最小変化方向にある１次
元データを選択し、（３）上記１次元データについて、着目点上の画像輝度
と該１次元データ上の点の輝度の差、および画像輝度の
変化程度を表す統計量に応じてウエイトを決める適応型
ウエイト付け加算平滑化を行い出力値とする画像処理方
法であって、ｉ，ｊは画像の着目点の位置を示し、ｑは画像のｉ行ｊ
列点における最小変化方向を所定の数値で表した値、ｍ
は画像のｉ行ｊ列点を中心としてｑ方向の１次元直線上
の所定のサンプリング点を順序づけた順番、Ｇ（ｉ，
ｊ）は画像のｉ行ｊ列点における輝度を表し、Ｉ（ｉ，
ｊ，ｑ，ｍ）は画像のｉ行ｊ列点における最小変化方向
であるｑ方向にｍサンプリング点行った点の輝度を表
し、Ｗ（ｉ，ｊ，ｑ，ｍ）は平滑化でＩ（ｉ，ｊ，ｑ，
ｍ）を加算する際のウエイトを示し、ｗ（ｉ，ｊ，ｑ，
ｍ）とｗｔはウエイトＷ（ｉ，ｊ，ｑ，ｍ）を計算する
際の途中項で、ｋは所定の定数で、＊＊ｋは、ｋ乗を意
味し、αは平滑化の程度を決める所定の定数であり、Σ
はｍについて所定の点数全ての和を取るものとした場合
に、前記ウエイトＷは、Ｗ（ｉ，ｊ，ｑ，ｍ）＝ｗ（ｉ，ｊ，ｑ，ｍ）／ｗ
ｔｗ（ｉ，ｊ，ｑ，ｍ）＝１／（１＋（（Ｉ
（ｉ，ｊ，ｑ，ｍ）−Ｇ（ｉ，ｊ））／（α・ρ））＊
＊ｋ）ｗｔ＝ Σ ｗ（ｉ，ｊ，ｑ，ｍ）によって与えられることを特徴とする画像処理方法。
【請求項２】画像の各点においてそれぞれ最小変化方向
を検出する検出処理部と、画像各点でそれぞれ上記最小変化方向にある１次元デー
タを選択する選択処理部と、上記１次元データについて、着目点上の画像輝度と該１
次元データ上の点の輝度の差、および画像輝度の変化程
度を表す統計量に応じてウエイトを決める適応型ウエイ
ト付け加算平滑化を行い出力値とする平滑化処理部とを
備えた画像処理装置であって、ｉ，ｊは画像の着目点の位置を示し、ｑは画像のｉ行ｊ
列点における最小変化方向を所定の数値で表した値、ｍ
は画像のｉ行ｊ列点を中心としてｑ方向の１次元直線上
の所定のサンプリング点を順序づけた順番、Ｇ（ｉ，
ｊ）は画像のｉ行ｊ列点における輝度を表し、Ｉ（ｉ，
ｊ，ｑ，ｍ）は画像のｉ行ｊ列点における最小変化方向
であるｑ方向にｍサンプリング点行った点の輝度を表
し、Ｗ（ｉ，ｊ，ｑ，ｍ）は平滑化でＩ（ｉ，ｊ，ｑ，
ｍ）を加算する際のウエイトを示し、ｗ（ｉ，ｊ，ｑ，
ｍ）とｗｔはウエイトＷ（ｉ，ｊ，ｑ，ｍ）を計算する
際の途中項で、ｋは所定の定数で、＊＊ｋは、ｋ乗を意
味し、αは平滑化の程度を決める所定の定数であり、Σ
はｍについて所定の点数全ての和を取るものとした場合
に、前記ウエイトＷは、Ｗ（ｉ，ｊ，ｑ，ｍ）＝ｗ（ｉ，ｊ，ｑ，ｍ）／ｗ
ｔｗ（ｉ，ｊ，ｑ，ｍ）＝１／（１＋（（Ｉ
（ｉ，ｊ，ｑ，ｍ）−Ｇ（ｉ，ｊ））／（α・ρ））＊
＊ｋ）ｗｔ＝ Σ ｗ（ｉ，ｊ，ｑ，ｍ）によって与えられることを特徴とする画像処理装置。
【請求項３】（１）画像の各点においてそれぞれ最小変
化方向を検出し、（２）画像各点でそれぞれ上記最小変化方向にある１次
元データを選択し、（３）上記１次元データについて、着目点上の画像輝度
と該１次元データ上の点の輝度の差、および画像輝度の
変化程度を表す統計量に応じてウエイトを決める適応型
ウエイト付け加算平滑化を行い出力値とする画像処理方
法であって、前記ウエイトは、平滑化する該１次元データ上の点の位
置が方向を検出した該着目点の位置から離れるにつれて
ウエイトの値が小さくなり、さらに、ｉ，ｊは画像の着目点の位置を示し、ｑは画像
のｉ行ｊ列点における最小変化方向を所定の数値で表し
た値、ｍは画像のｉ行ｊ列点を中心としてｑ方向の１次
元直線上の所定のサンプリング点を順序づけた順番、Ｇ
(ｉ，ｊ)は画像のｉ行ｊ列点における輝度を表し、Ｉ
（ｉ，ｊ，ｑ，ｍ）は画像のｉ行ｊ列点における最小変
化方向であるｑ方向にｍサンプリング点行った点の輝度
を表し、Ｗ（ｉ，ｊ，ｑ，ｍ）は平滑化でＩ（ｉ，ｊ，
ｑ，ｍ）を加算する際のウエイトを示し、ｗ（ｉ，ｊ，
ｑ，ｍ）とｗｔはウエイトＷ（ｉ，ｊ，ｑ，ｍ）を計算
する際の途中項で、ｋは所定の定数で、＊＊ｋは、ｋ乗
を意味し、αは平滑化の程度を決める所定の定数であ
り、Σはｍについて所定の点数全ての和を取るものとし
た場合に、前記ウエイトＷは、Ｗ（ｉ，ｊ，ｑ，ｍ）＝ｗ（ｉ，ｊ，ｑ，ｍ）／ｗ
ｔｍ≧１の時、ｗ（ｉ，ｊ，ｑ，ｍ）＝ｗ（ｉ，ｊ，ｑ，ｍ−１）
／（１＋（（Ｉ（ｉ，ｊ，ｑ，ｍ）−Ｇ（ｉ，ｊ））／
（α・ρ））＊＊ｋ）ｍ＝０の時、ｗ（ｉ，ｊ，ｑ，０）＝１ｍ≦−１の時、ｗ（ｉ，ｊ，ｑ，ｍ）＝ｗ（ｉ，ｊ，ｑ，ｍ＋１）
／（１＋（（Ｉ（ｉ，ｊ，ｑ，ｍ）−Ｇ（ｉ，ｊ））／
（α・ρ））＊＊ｋ）ｗｔ＝ Σ ｗ（ｉ，ｊ，ｑ，ｍ）によって与えられることを特徴とする画像処理方法。
【請求項４】上記請求項３において、所定の定数ｋの値
を２とすることを特徴とする画像処理方法。
【請求項５】画像の各点においてそれぞれ最小変化方向
を検出する検出処理部と、画像各点でそれぞれ上記最小変化方向にある１次元デー
タを選択する選択処理部と、上記１次元データについて、着目点上の画像輝度と該１
次元データ上の点の輝度の差、および画像輝度の変化程
度を表す統計量に応じてウエイトを決める適応型ウエイ
ト付け加算平滑化を行い出力値とする平滑化処理部とを
備えた画像処理装置であって、前記ウエイトは、平滑化する該１次元データ上の点の位
置が方向を検出した該着目点の位置から離れるにつれて
ウエイトの値が小さくなり、さらに、ｉ，ｊは画像の着目点の位置を示し、ｑは画像
のｉ行ｊ列点における最小変化方向を所定の数値で表し
た値、ｍは画像のｉ行ｊ列点を中心としてｑ方向の１次
元直線上の所定のサンプリング点を順序づけた順番、Ｇ
(ｉ，ｊ)は画像のｉ行ｊ列点における輝度を表し、Ｉ
（ｉ，ｊ，ｑ，ｍ）は画像のｉ行ｊ列点における最小変
化方向であるｑ方向にｍサンプリング点行った点の輝度
を表し、Ｗ（ｉ，ｊ，ｑ，ｍ）は平滑化でＩ（ｉ，ｊ，
ｑ，ｍ）を加算する際のウエイトを示し、ｗ（ｉ，ｊ，
ｑ，ｍ）とｗｔはウエイトＷ（ｉ，ｊ，ｑ，ｍ）を計算
する際の途中項で、ｋは所定の定数で、＊＊ｋは、ｋ乗
を意味し、αは平滑化の程度を決める所定の定数であ
り、Σはｍについて所定の点数全ての和を取るものとし
た場合に、前記ウエイトＷは、Ｗ（ｉ，ｊ，ｑ，ｍ）＝ｗ（ｉ，ｊ，ｑ，ｍ）／ｗ
ｔｍ≧１の時、ｗ（ｉ，ｊ，ｑ，ｍ）＝ｗ（ｉ，ｊ，ｑ，ｍ−１）
／（１＋（（Ｉ（ｉ，ｊ，ｑ，ｍ）−Ｇ（ｉ，ｊ））／
（α・ρ））＊＊ｋ）ｍ＝０の時、ｗ（ｉ，ｊ，ｑ，０）＝１ｍ≦−１の時、ｗ（ｉ，ｊ，ｑ，ｍ）＝ｗ（ｉ，ｊ，ｑ，ｍ＋１）
／（１＋（（Ｉ（ｉ，ｊ，ｑ，ｍ）−Ｇ（ｉ，ｊ））／
（α・ρ））＊＊ｋ）ｗｔ＝ Σ ｗ（ｉ，ｊ，ｑ，ｍ）によって与えられることを特徴とする画像処理装置。