JPH1049670A - 画像処理方法 - Google Patents
画像処理方法Info
- Publication number
- JPH1049670A JPH1049670A JP8204490A JP20449096A JPH1049670A JP H1049670 A JPH1049670 A JP H1049670A JP 8204490 A JP8204490 A JP 8204490A JP 20449096 A JP20449096 A JP 20449096A JP H1049670 A JPH1049670 A JP H1049670A
- Authority
- JP
- Japan
- Prior art keywords
- value
- image
- input image
- predetermined
- minimum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000012545 processing Methods 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims description 69
- 238000009499 grossing Methods 0.000 claims abstract description 77
- 238000003672 processing method Methods 0.000 claims description 22
- 239000013598 vector Substances 0.000 claims description 21
- 238000012935 Averaging Methods 0.000 claims description 5
- 230000003044 adaptive effect Effects 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 9
- 238000005070 sampling Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 5
- 210000001015 abdomen Anatomy 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000001186 cumulative effect Effects 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 230000033772 system development Effects 0.000 description 1
Landscapes
- Facsimile Image Signal Circuits (AREA)
- Image Processing (AREA)
Abstract
さない画像処理を行うことである。 【解決手段】入力画像の各点において、ステップ101
でエッジ部に有効な非線形平滑化処理を行いエッジ部に
有効な画像を作成し、ステップ102で平坦部に有効な
画像を作成する。ステップ103で画像のエッジ部であ
るか平坦部であるか認識し、ステップ104で、画像の
エッジ部はエッジ部に有効な画像を、平坦部は平坦部に
有効な画像が出力画像となるようにデータの切り張り処
理を行う。エッジ部に有効な処理は、参照点数の少ない
1次元の非線形平滑化処理を行い、平坦部に有効な画像
は参照点数の大きい1次元の非線形平滑化処理を行う。
Description
して所定の処理を施す画像処理に関する発明である。
理方式」には、画像の各点で、複数の所定方向について
それぞれ画像の輝度変化の大きさを評価し、輝度変化が
最小となる方向(最小変化方向)を画像の各点毎に求
め、最小変化方向に1次元の非線形平滑化を行うものが
記載されている。
傍の点がないために処理できない画像の端の部分は除
く)で、皆同じ非線形処理が行われることが特徴であ
る。
像上の全ての点で、皆同じ非線形処理を行うようになっ
ており、一枚の画像でもそれぞれかけ離れた特徴を持つ
画像の平坦部とエッジ部について、どちらも同じ非線形
処理を行うことになっている。従って、上記従来技術は
画像の中で大きく特徴の異なるエッジ部と平坦部に着目
すると、それぞれの特徴に会った最適な処理からかけ離
れた処理をするという問題点があった。本発明の課題
は、上記従来技術の問題点を考慮し、上記従来技術の処
理よりもさらに画質を改善することにある。
ため、入力画像のエッジ部と平坦部を分離して認識する
処理を行い、エッジ部はエッジ部、平坦部は平坦部に適
した処理を行った値が出力画像値として得られるように
した。
あり、具体的には下記の方法を実施するプログラムまた
はハードを作成することにより実現できる。以下、本発
明の実施例を(1)〜(5)節に分けて説明する。
(1)節では本発明の実施例の手順の概略を説明し、
(2)〜(4)節で概略で示したそれぞれの部分の手順
を詳しく説明し、(5)節でその他の実施例を簡単に説
明する。
をステップ番号として以下に説明する。
非線形平滑化処理を入力画像全体に行い、結果画像A
(i,j)を得る。
線形平滑化処理を入力画像全体に行い、結果画像B(i,
j)を得る。
坦部を区別し、分離して認識する処理を行う。認識結果
に基づいて、エッジ部と平坦部を区別するマップM(i,
j)を作成する。例えば、エッジ部と認識した部分のマッ
プ値は、M(i,j)=1とし、平坦部と認識した部分のマ
ップ値は、M(i,j)=0とする。ただし、i,jは画像のi
行j列を意味する。
たエッジ部と平坦部を区別するマップM(i,j)の値に応
じてエッジ部に有効な処理画像A(i,j)か平坦部に有効
な処理画像B(i,j)かの切り張り処理を行う。すなわち
以下の数1の処理を行い出力画像G(i,j)を得る。 M(i,j)=1の場合(エッジ部の場合): G(i,j) = A(i,j)、 M(i,j)=0の場合(平坦部の場合): G(i,j) = B(i,j)…(数1)。
02は順番が交換でき、ステップ102の方を先に行っ
てもよい。また、ステップ103も、ステップ101や
ステップ102の結果を参照しない場合は、ステップ1
01やステップ102と順番を交換できる。
を得る他の方法もある。例えば、以上の手順では、それ
ぞれの処理画像A(i,j)、 B(i,j)を作成していたが、
両画像を作成せずに、初めに入力画像のエッジ部と平坦
部を区別し、分離して認識する処理を行い、着目した画
素がエッジ部と判断した点はエッジ部に有効な非線形平
滑化処理を行い着目点の出力画像値を得、平坦部と判断
した場合は、平坦部に有効な非線形平滑化処理を行い着
目点の出力画像値を得て、次の着目点に進む手順を行う
こともできる。
ばどのような処理も用いることができる。一例として上
記従来技術の方法を用いた例を図2に示し、図中の番号
をステップ番号として以下説明する。
中心に複数の所定の方向について輝度変化を評価し、各
方向の輝度変化値を求める。図6は、所定の輝度変化を
評価する方向として8方向取り、各方向に5サンプリン
グ点分の輝度変化の評価をする場合について、各方向と
サンプリング点の例を図示したものである。評価する方
向の数、サンプリング点の数はこの他にも各種考えられ
る。なお、図6で示した黒丸の点は、元々のデータ値が
ある点であるが、白丸で示した点は元々のデータ値がな
いので周りの点の値から補間を要する。補間法は各種あ
り、例えば、最も近い2点の値の平均値とする方法や、
最も近い2点のうち、中心から見て上下か左右方向に近
い方の点の値とする方法などがあり、どのような方法を
用いることができる。
方法は各種あるが、一例として以下の数2を挙げる。
けるd方向の輝度変化を捉えた値であり、dは図6で示
された8方向について各方向に順番を付けた数値であ
り、例えば右上と左下方向をd=0と定義し、それから
時計周りに順次1〜7番と定義する。 absは絶対値を意
味し、 dirval(I(i,j),d,m)は、 入力画像Iのi行j
列の点からd方向にmサンプリング点行った点の値を意
味する。なお、 S(d,i,j)はi行j列の点における各輝
度変化を成分値とする輝度変化を示すベクトルと見なす
ことができるため、 i行j列の点における輝度変化ベ
クトルと呼ぶことにする。
大きさを評価する計算方法は各種ある。例えば数2では
dirval(I(i,j),d,m)-I(i,j))の絶対値を取っていたが
dirval(I(i,j),d,m)-I(i,j))の2乗とすることも考え
られる。どちらの評価式を用いても以後の結果に大差は
生じなかった。
各方向の輝度変化値の中から最小値を求め、その方向
(最小変化方向)を求める。図6の8方向の輝度変化の
中から最小値となる方向(最小変化方向)を求める場合
を式で示すと以下の数3になる。 dmin(i,j) = mindir(S(d,i,j)|0<=d<=7) …(数3) ただし、mindir(S(d,i,,j)|0<=d<=7)は、i行j列の輝
度変化ベクトルS(d,i,,j)の中から、0<=d<=7で示され
る8方向捉えた輝度変化値( S(0,i,,j)〜 S(7,i,,
j))の中で最小値となるものの方向(dの値)を求める
ことを意味する。最小値を持つ方向が2つ以上ある場合
は、方向を数字として表した時の数値の少ない方向とし
てもよいし、逆に大きい数値の方向としてもどちらでも
構わない。
元の非線形平滑化を行う。具体的な方法は各種あるがこ
こでは上記従来技術で示した方法を数4〜数7に示し、
別の方法を数8〜数11に示す。
元の非線形平滑化処理を示す。平滑化の方向は、図6で
示された8方向の中から最小変化方向dmin(i,j)が選ば
れており、平滑化の点数は、着目した画素 (i,j)点を
中心に最小変化方向dmin(i,j)に5点取る場合について
以下示す。
は、入力画像Iのi行j列の点から数3で求めた最小変
化方向であるdmin(i,j)方向にmサンプリング点行った
点における輝度値を意味する。Hmin(i,j)は、i行j列
の点における平滑化後の出力値を意味し、σ0は入力画
像のノイズの標準偏差に比例する量を計算した値を意味
し、αは平滑化の程度を決める所定の値である。なお、
σ0の計算の仕方は各種あり、例えばノイズ部とみなせ
る画像の端の小さな領域の輝度の2乗平均の平方根を用
いてもよいし、画像各点における近傍の点の中で最も変
位の小さい値の平均値を用いてもよい。
滑化の他にもいろいろな非線形平滑化を行うことができ
る。たとえば、以下の数8〜数11の1次元非線形平滑
化を行うこともできる。
の出力値を意味し、σ0は入力画像のノイズの標準偏差
に比例する量を計算した値を意味し、αは平滑化の程度
を決める所定の値であり、βは平滑化の寄与を決める所
定のパラメータ値である。
形平滑化を行う方法は各種ある。例えば所定の点数の1
次元データの中から中間値を選ぶこともできる。また、
入力データが0と1の2値画像で出力も0が1を取るよ
うな場合は、1次元データの中から数の多い方の値を選
ぶことなどもできる。
2までの5点であるため、上記数4から数7までは、最
小変化方向に5点の1次元非線形平滑化をしたものとな
っている。同様に数8から数11で示した平滑化も5点
の1次元の平滑化をしたものとなっている。平滑化を行
う点数は5点に限らず、どのように取ることもできる
が、エッジ部では、平坦部よりも短い距離の平滑化をし
た方が好ましく、ここでは一例として5点とした。通常
の臨床撮影に使われる頭部や腹部のMRI画像で評価し
た所、エッジ部に有効な平滑化は、3点または5点であ
り、3点の場合はエッジ部の微妙な変化をよく残す反
面、平滑化の効果が少なかった。5点では、エッジ部の
微妙な変化を3点の時ほど残さないが、ほぼ満足できる
程度に画像の特徴をつかむために必要なエッジ部の微妙
な変化を残しており、平滑化の効果も高く、満足いく結
果であった。骨部のMRI画像で評価した所、平滑化は
あまりしない方が好ましく、5点より3点の方が好まし
かった。その他にも、撮影条件の違いによる画像のS/
Nや撮影部位などで最適な点数も異なることが予想され
た。
101の具体的方法について述べたが、上記方法では、
S/Nの低く、画像のエッジがノイズであいまいになっ
た画像を処理すると、画像のエッジ部で希にぶつぶつし
た人工色の強い感触の得られる欠点がある。S/Nの低
く画像のエッジがノイズであいまいになるような画像に
たいしては次に述べる図3の処理がエッジ部のなめらか
さを出す上で有効であった。以下、図3の処理について
簡単に述べる。
たステップ201と同一の処理を行う。
変化ベクトルを平均してボケた輝度変化ベクトル S'
(d,i,j)を作成する。近傍の取り方や平均の仕方などに
ついては各種考えられる。一例として簡単な平均を行う
場合を以下の数12に示す。
2の説明と同様に、ボケた輝度変化ベクトル S'(d,i,
j)について、各画素毎に、各方向のボケた輝度変化ベク
トルの値の中から最小値を求め、その方向(ボケ最小変
化方向)を求める。数3と同様に式で表すと、ボケ最小
変化方向dmin'(i,j) は以下の数13になる。
様にボケ最小変化方向dmin'(i,j)に1次元の非線形平滑
化を行う。上式の数5、数7、数8、数11中にあるdm
in(i,j)の代わりにdmin'(i,j)とするだけで、本ステッ
プ303の式になる。
は、 S/Nの低く画像のエッジがノイズであいまいに
なるような画像にたいしても、画像のエッジ部がなめら
かに繋がった処理画像になるという利点があった。
像の端の領域で近傍の点のデータがないことから参照デ
ータがないことから生じる特例処理について触れていな
いのでここで記す。画像の端の領域で必要な特例処理の
具体的な仕方は各種ある。例えば、5点の平滑化を行う
場合、画像の端から2点分の領域はデータがない部分が
生じるが、その領域は、原画の値を出力画像の値とする
方法がある。その外、画像の両端がサイクリックに接続
されているものと仮定して、サイクリックに接続した時
のデータ値を参照して処理する方法などもある。
どのような処理も用いることができる。一例として上記
(2)節で説明した1次元非線形平滑化を用いた場合に
ついて説明する。
号をステップ番号として以下に示す。
様に輝度変化ベクトルを求める。ただし、ステップ20
1では各方向5点の輝度変化を評価したが、ここでは、
各方向9点の輝度変化を評価する。従って、数2のmの
範囲は−2<=m<=2(1次元の評価5点)であった
が、−4<=m<=4(1次元の評価9点)として輝度
変化ベクトルを求めることができる。
様に、輝度変化ベクトルから最小変化方向を求める。
様に、最小変化方向に1次元の非線形平滑化を行う。た
だし、ステップ203では1次元の非線形の平滑化を5
点で行っていたが、ここでは9点で行う。
合では、数5と数7について和を取る範囲を−2<=m
<=2(5点)から−4<=m<=4(9点)にし、数
5の分母を5から9にすれば平滑化の点数が9点の1次
元非線形平滑化を行う式になる。
を用いる場合は、数8と数10のmの値を−4<=m<
=4(9点)にすることで本ステップの処理の式とな
る。
点数をエッジ部の同様の処理よりも大きく設定するほど
平滑化の効果が増大するが、あまり大きな点数を取りす
ぎても本来あるべき濃淡が消え、逆効果の方が大きい場
合も生じる。通常臨床で用いられる頭部や腹部のMRI
画像に対して処理を行い、評価した結果では、平坦部の
平滑化は9点か11の場合が良好であった。9点の場合
は無難な画像が得られ、11点の場合は多少平滑化しす
ぎている感じもしたが、最良な点数については、一般に
は処理を行う対象とする画像で異なることが予想され
る。
5で示す処理では、図4のステップ401で輝度変化を
評価する点数を各方向9点としておたが、ここでは、ス
テップ201と同様にして輝度変化を評価する点数を各
方向5点とするものである。他の処理は図4の場合と同
様である。輝度変化を各方向5点で評価したので、次の
ステップ402で求められる最小変化方向は、図5の場
合は、輝度変化の評価を5点で行った場合の最小変化方
向である。従って、次のステップ403でも図5の場合
は、輝度変化の評価を5点で行った場合の最小変化方向
に9点の1次元非線形平滑化を行うため、図4の場合と
は平滑化を行う方向が異なる画素が生じる。図5の処理
の結果は、図4の処理結果に比べ、画質で若干劣ってい
たが、大差がない。図5の処理は、輝度変化をエッジ部
と同じ5点で評価するため、エッジ部で計算した輝度変
化の評価結果を再度利用することができるため、演算を
その分高速化することができる利点がある。
点または9点で行い、その最小変化方向にエッジ部より
大きな9点の1次元非線形平滑化を行う場合を示した
が、その他にも、平坦部に有効な非線形平滑化処理とし
ては最小変化方向の1次元の非線形平滑化だけでなく、
最小方向と45度づれた方向の1次元平滑化を行い、両者
の平均値としたり、平面内の適当な重み付き平均値とす
るなど、1次元方向だけでない平滑化を行うこともでき
る。
て認識する方法であればどのようなものでも用いること
ができる。例えばエッジ画像を利用する方法がある。以
下、処理方法の具体例を説明する。
像E(i,j)を作成する処理を行う。
j+q列の値である。 f(p,q)の具体的な値は、E(i,j)
がエッジ画像となるものであれば、どのような値を用い
てかまわない。たとえば、以下の数15としたり、以下
の数16としたりすることができる。
E(i,j)の値が所定のしきい値よりも大きい場合にはエ
ッジ部分と見做し、マップの値M(i,j)を1とする。 E
(i,j)の値が所定のしきい値よりも低い場合には平坦部
と見做し、マップの値M(i,j)を0とする。ただし、上
記しきい値の設定の仕方は各種ある。
ップ値を表示して適当であると判定するまでインターラ
クティブに変えて適切なしきい値を探す方法を用いるこ
ともできるし、画像のノイズの大きさを評価してノイズ
の標準偏差値の所定倍(例えば3倍)に設定する方法を
用いることもできるし、エッジ画像のヒストグラムを求
めてその累積ヒストグラムが所定の割合(例えば80
%)になる値にすることもできるし、累積ヒストグラム
の80%になる値がノイズの標準偏差の3倍よりも小さ
い場合はしきい値をノイズの3倍とし、それ以外は、累
積ヒストグラムの80%になる値とノイズの標準偏差の
3倍の値との平均値をしきい値にすることなど、しきい
値の設定の仕方は各種方法が考えられる。
てM(i,j)の値が1ならばその近傍の新たなマップ値M'
(i+p,j+q)を1とし、残りの点の新たなマップ値M'(i,
j)は0とする処理を行い、エッジ部位の領域を増やした
処理を行うこともできる。この場合エッジが不安定であ
ってもエッジ部をより確実にエッジ部とみなすことがで
きるという利点がある。
理方法は各種ある。例えば、以下のような方法を用い、
よりエッジ部を明確に認識することもできる。
り、エッジ部に有効な非線形処理画像A(i,j)ができて
いるものとする。また、ステップ101の非線形平滑化
が1次元の非線形平滑化で、非線形の平滑化をする方向
である最小変化方向dmin(i,j) (またはボケ最小変化方
向dmin'(i,j) )がすでに計算されており、ここで参照
できるとする。以下、ステップ101の1次元非線形平
滑化を行う方向がボケ最小変化方向dmin'(i,j)であると
して記す。
行方向qd(i,j)を計算する。方向を0〜7までの8方向
で表す場合、以下の数17で直行方向qd(i,j)を求める
ことができる。
ある。
j)に所定の点数、例えば5点のデータを取り、その中で
最大値と最小値の差をエッジ画像E' (i,j)とする。す
なわち、式で書くと以下の数18となる。
ジ画像E(i,j)と考えて同様の処理を行うとマップの値
M(i,j)が求まる。このようにE' (i,j)を作成して求め
ると、 上記E(i,j)を作成してマップを求めるよりも、
エッジ部を安定して認識できるという効果がある。
かを0と1で分けて認識したが、このほかにも、エッジ
部らしら、平坦部らしさを表した度合いマップM" (i,
j)を作成する方法がある。以下、ステップ101とステ
ップ102は上記と同一で、ステップ103とステップ
104を上記から変えた一例を示す。
画像E(i,j)を基に、以下の数19ようにエッジ部と平
坦部の度合いマップM" (i,j)を作る。
準偏差に比例する量を計算した値を意味し、γはエッジ
と平坦部の度合いを決める所定のパラメータ値である。
なお、σ0の計算の仕方は各種あり、例えばノイズ部と
みなせる画像の端の小さな領域の輝度の2乗平均の平方
根を用いてもよいし、画像各点における近傍の点の中で
最も変位の小さい値の平均値を用いてもよい。
を行う。
部を分離して各々独立な処理ができることから、それぞ
れの部位の特徴に適した処理を行えるようになり、従来
よりも画質を改善した出力画像を得られるという効果が
ある。
した図である。
を詳しく説明した図である。
した図である。
を詳しく説明した図である。
を示す図である。
の例を示す図である。
らに別の例を示す図である。
プ、102…平坦部に有効な非線形平滑化をするステッ
プ、103…エッジ部と平坦部を分離するステップ、1
04…エッジ部と平坦部のデータ切り張り処理をするス
テップ、201…各方向5点の輝度変化を評価するステ
ップ、202…最小変化方向を求めるステップ、203
…最小変化方向に5点の1次元非線形平滑化をするステ
ップ、311…ボケた輝度変化ベクトルを作成するステ
ップ、302…ボケ最小変化方向を求めるステップ、3
03…ボケ最小変化方向の1次元非線形平滑化をするス
テップ、401…各方向9点の輝度変化を評価するステ
ップ、402…最小変化方向を求めるステップ、403
…最小変化方向に9点の1次元非線形平滑化をするステ
ップである。
Claims (10)
- 【請求項1】入力画像のエッジ部と平坦部を分離して認
識し、 認識された結果に応じて、前記入力画像のエッジ部と認
識した部位にはエッジ部用の所定方向の所定の1次元非
線形平滑化を行い、 前記入力画像の平坦部と認識した部位には平坦部用の所
定の平滑化処理を行うことを特徴とする画像処理方法。 - 【請求項2】請求項1に記載の画像処理方法において、 前記エッジ部用の所定方向の1次元非線形平滑化の処理
として、前記入力画像の着目した画素を中心に複数の所
定方向について輝度変化を評価し、各輝度変化の値を基
に特定の1方向を選択し、前記選択された1方向上の該
入力画像値のデータを基に1次元の非線形平滑化処理を
行うことを特徴とした画像処理方法。 - 【請求項3】請求項1に記載の画像処理方法において、 前記エッジ部用の所定方向の1次元非線形平滑化の処理
として、前記入力画像の着目した画素を中心に複数の所
定方向について輝度変化を評価し、各輝度変化値の中か
ら最小値を取る方向である最小変化方向を選択し、前記
最小変化方向上の前記入力画像値のデータを基に1次元
の非線形平滑化処理を行うことを特徴とした画像処理方
法。 - 【請求項4】請求項1に記載の画像処理方法において、 上記エッジ部用の所定方向の所定の1次元非線形平滑化
の処理として、該入力画像について、着目した画素を中
心に複数の所定方向について輝度変化を評価して各方向
の輝度変化値を求め、該各方向の輝度変化値を成分値と
する輝度変化ベクトを作成し、該入力画像の各画素に対
応した該輝度変化ベクトルを各々の画素について作成
し、着目した画素とその周辺の画素の該輝度変化ベクト
ルについて、所定の重み付けベクトル加算を行って着目
画素のボケ輝度変化ベクトルを作成し、該着目画素のボ
ケ輝度変化ベクトルの成分値の中から最小値となる成分
を与えた方向、すなわちボケ最小変化方向を求め、該ボ
ケ最小変化方向上の該入力画像値のデータを基に1次元
の非線形平滑化処理を行うことを特徴とした画像処理方
法。 - 【請求項5】請求項1に記載の画像処理方法において、 上記平坦部用の非線形平滑化の処理として、 該入力画像の着目した画素を中心に複数の所定方向につ
いて輝度変化を評価し、各輝度変化の値を基に特定の1
方向を選択し、該選択された1方向について、エッジ部
の同様の処理よりも長い距離にある該入力画像値のデー
タを基に、エッジ部の同様の処理よりも長い距離の1次
元の非線形平滑化処理を行うことを特徴とした画像処理
方法。 - 【請求項6】請求項1に記載の画像処理方法において、 上記平坦部用の所定の非線形平滑化の処理として、 該入力画像の着目した画素を中心に複数の所定方向につ
いて輝度変化を評価し、各輝度変化の値を基に特定の1
方向を選択し、該選択された1方向について、エッジ部
の同様の処理よりも長い距離にある該入力画像値のデー
タを基に、エッジ部の同様の処理よりも長い距離の1次
元の非線形平滑化処理を行った値と、上記選択された1
方向と所定方向ずれた方向にある該入力画像値のデータ
を基に1次元の非線形平滑化を行った値を平均する処理
を行うことを特徴とした画像処理方法。 - 【請求項7】請求項1に記載の画像処理方法において、 上記エッジ部用の所定の非線形平滑化の処理として、 該入力画像の着目した画素を中心に複数の所定方向につ
いて輝度変化を評価し、各輝度変化の値を基に特定の1
方向を選択し、該選択された1方向上の該入力画像値の
データを基に1次元の非線形平滑化処理を行い、上記平
坦部に適した所定の非線形平滑化の処理として、該入力
画像の着目した画素を中心にした周囲の該入力画像値を
基に、各入力画像値の寄与する重みを変えた適合型重み
付き平均を行うことを特徴とした画像処理方法。 - 【請求項8】1)入力画像の各点において、着目した画
素を中心に複数の所定方向について輝度変化を評価して
各方向の輝度変化値を求め、該各方向の輝度変化値を成
分値とする輝度変化ベクトを作成し、 2)該入力画像の各点について、着目した画素とその周
辺の画素の該輝度変化ベクトルを基に所定の重み付けベ
クトル加算を行って着目画素のボケ輝度変化ベクトルを
作成し、 3)該入力画像の各点について、着目画素のボケ輝度変
化ベクトルの成分値の中から最小値となる成分を与えた
方向、すなわちボケ最小変化方向を求め、 4)該入力画像の各点について、該ボケ最小変化方向上
の該入力画像値のデータを基に1次元の非線形平滑化処
理を行い、ボケ最小変化方向平滑化画像を作成し、 5)該ボケ最小変化方向平滑化画像を基に画像のエッジ
部と平坦部を区別するマップを作成し、 6)該エッジ部と平坦部を区別するマップを参照して、
画像のエッジ部と判断した部分では、該ボケ最小変化方
向平滑化画像を出力画像値とし、 7)該エッジ部と平坦部を区別するマップを参照して、
画像の平坦部分と判断した部分では、平坦部用の所定の
非線形平滑化の処理を行った結果を出力画像値とするこ
とを特徴とした画像処理方法。 - 【請求項9】請求項8に記載の画像処理方法おいて、 前記マップを作成する処理として該ボケ最小変化方向平
滑化画像の各点について、着目した画素に対応する該入
力画像の該ボケ最小変化方向と直行した方向を作り、該
ボケ最小変化方向平滑化画像の着目点を中心として該直
行方向上にある所定の点数分の該ボケ最小変化方向平滑
化画像値のデータ値の中から最大値と最小値の差を求
め、該最大値と最小値の差が所定の値よりも大きい場合
はエッジ部と判定し、該最大値と最小値の差が所定の値
よりも小さい場合は平坦部と判定する処理を行うことを
特徴とした画像処理方法。 - 【請求項10】請求項8に記載の画像処理方法におい
て、 前記画像の平坦部用の所定の非線形平滑化の処理とし
て、 該輝度変化ベクトの成分値の中から最小値となる成分を
与えた方向である最小変化方向を求め、該最小変化方向
上の該入力画像値のデータを基に、エッジ部に行う処理
よりも長い距離の1次元の非線形平滑化処理を行うこと
を特徴とした画像処理方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP20449096A JP3462960B2 (ja) | 1996-08-02 | 1996-08-02 | 画像処理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP20449096A JP3462960B2 (ja) | 1996-08-02 | 1996-08-02 | 画像処理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH1049670A true JPH1049670A (ja) | 1998-02-20 |
JP3462960B2 JP3462960B2 (ja) | 2003-11-05 |
Family
ID=16491398
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP20449096A Expired - Fee Related JP3462960B2 (ja) | 1996-08-02 | 1996-08-02 | 画像処理方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3462960B2 (ja) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001143068A (ja) * | 1999-10-08 | 2001-05-25 | Eastman Kodak Co | 領域成長に基づくデジタル画像ノイズ低減法 |
JP2001155148A (ja) * | 1999-10-08 | 2001-06-08 | Eastman Kodak Co | 形状可変の画素近隣領域を用いたデジタル画像からのクロマノイズの除去 |
US6826291B2 (en) | 1997-09-03 | 2004-11-30 | Hitachi, Ltd. | Method and system for embedding information into contents |
-
1996
- 1996-08-02 JP JP20449096A patent/JP3462960B2/ja not_active Expired - Fee Related
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6826291B2 (en) | 1997-09-03 | 2004-11-30 | Hitachi, Ltd. | Method and system for embedding information into contents |
JP2001143068A (ja) * | 1999-10-08 | 2001-05-25 | Eastman Kodak Co | 領域成長に基づくデジタル画像ノイズ低減法 |
JP2001155148A (ja) * | 1999-10-08 | 2001-06-08 | Eastman Kodak Co | 形状可変の画素近隣領域を用いたデジタル画像からのクロマノイズの除去 |
JP4688269B2 (ja) * | 1999-10-08 | 2011-05-25 | イーストマン コダック カンパニー | 領域成長に基づくデジタル画像ノイズ低減法 |
Also Published As
Publication number | Publication date |
---|---|
JP3462960B2 (ja) | 2003-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2018223994A1 (zh) | 中文打印字符图像合成方法及装置 | |
JP4738469B2 (ja) | 画像処理装置、画像処理プログラムおよび画像処理方法 | |
EP1310912A2 (en) | Image processing method, apparatus and system | |
JPH0877334A (ja) | 顔画像の特徴点自動抽出方法 | |
JP2007213377A (ja) | 顔特徴点検出方法および装置並びにプログラム | |
CN111488836B (zh) | 人脸轮廓修正方法、装置、设备及存储介质 | |
WO2017061106A1 (ja) | 情報処理装置、画像処理システム、画像処理方法及びプログラム記録媒体 | |
JP4639754B2 (ja) | 画像処理装置 | |
JP3708042B2 (ja) | 画像処理方法及びプログラム | |
CN104268845A (zh) | 极值温差短波红外图像的自适应双局部增强方法 | |
JP3636936B2 (ja) | 濃淡画像の2値化方法および濃淡画像の2値化プログラムを記録した記録媒体 | |
JP3477603B2 (ja) | 画像処理方法、画像処理装置および画像処理プログラムを記録した記録媒体 | |
JPH1049670A (ja) | 画像処理方法 | |
JPH1049669A (ja) | 画像処理方法 | |
JP2775122B2 (ja) | イラストデータの自動輪郭抽出ベクトル化処理方法,及びそれに使用する処理装置 | |
JP3533050B2 (ja) | 画像領域分割装置 | |
JPH07334648A (ja) | 画像処理方法及びその装置 | |
JP2011198006A (ja) | オブジェクト検出装置、オブジェクト検出方法、およびオブジェクト検出プログラム | |
CN112419208A (zh) | 一种基于施工图审查矢量图纸编译方法及系统 | |
JP2017138743A (ja) | 画像処理装置、画像処理方法及びプログラム | |
CN112085683A (zh) | 一种显著性检测中的深度图可信度检测方法 | |
JP2003173436A (ja) | 画像処理装置、画像処理方法及びそのプログラムをコンピュータ読み取り可能に記録した記録媒体 | |
EP4246428A1 (en) | Perspective method for physical whiteboard and generation method for virtual whiteboard | |
JP2006072829A (ja) | 画像認識システム及び画像認識方法 | |
JP2840485B2 (ja) | 画像処理方法及び装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080815 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090815 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100815 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100815 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110815 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120815 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120815 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130815 Year of fee payment: 10 |
|
LAPS | Cancellation because of no payment of annual fees |