画像処理方法及びコンピュータプログラム
【目的】入力画像に角度依存情報を付与することによって一次元的な模様を強調して表示させる。
【構成】まず、入力画像のデータをコンピュータに入力する(Sa1)。次に、コントラスト強度として任意の定数mを入力する(Sa2)。次に、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Dとして区画し、演算対象領域D内における注目画素とその周囲の任意の画素との2点の画素の傾きの強さを前記注目画素を中心とする全周に亘って合計し演算対象領域内の画素数で除算することによりその平均値を算出する(Sa3)。さらに、ステップ(Sa3)で求めた勾配の振幅値βにステップ(Sa2)で入力したコントラスト強度m(但し、mは正の定数である。)を乗算し、任意のオフセット値γを加算し(Sa4)、最終的な出力画像を得る(Sa5)。
【構成】まず、入力画像のデータをコンピュータに入力する(Sa1)。次に、コントラスト強度として任意の定数mを入力する(Sa2)。次に、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Dとして区画し、演算対象領域D内における注目画素とその周囲の任意の画素との2点の画素の傾きの強さを前記注目画素を中心とする全周に亘って合計し演算対象領域内の画素数で除算することによりその平均値を算出する(Sa3)。さらに、ステップ(Sa3)で求めた勾配の振幅値βにステップ(Sa2)で入力したコントラスト強度m(但し、mは正の定数である。)を乗算し、任意のオフセット値γを加算し(Sa4)、最終的な出力画像を得る(Sa5)。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、デジタル画像に含まれる特定の部分を抽出し、より強調して表示させたり、逆に、特定の部分を除去したりすることを可能にするデジタル画像の画像処理方法及びコンピュータプログラムに関するものである。
【背景技術】
【0002】
デジタル画像において、注目画素とその周囲の画素とを比較した際の両者のデジタル化された画像信号(例えば輝度値)の変化に関する情報を表すベクトルは一般に「微分画像」又は「勾配画像」と呼ばれている(本明細書では、「微分画像」として統一する。)。微分画像は、注目画素とその周辺の任意の点の2点(2画素)を比較して、座標及びデジタル値のそれぞれの変化から算出された「勾配ベクトル」を表している。勾配ベクトルを利用した画像処理アルゴリズムの一つとして、「アイリスフィルタ」が知られている。アイリスフィルタでは、その勾配ベクトルを求めたあと、そこからデジタル化された画像信号(例えば輝度値)の「集中度」を計算する。この集中度やその他の計算値を閾値処理することにより、マンモグラフィ診断画像の読影を支援するシステムに関するものが開示されている(例えば、特許文献1〜4等参照)
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開平09−167238号公報
【特許文献2】特開平10−289318号公報
【特許文献3】特開2005−52295号公報
【特許文献4】特開2006−130049号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
この種の読影システムは、腫瘤陰影の輝度値が放射状に減少する(すなわち、腫瘤陰影は中心ほど明るく周辺ほど暗いという特徴を有する)というモデルに基づいて設計されている。
【0005】
しかし、実際の腫瘤陰影には様々な形状があり、必ずしも上述のようなモデルで全てが解決できるほど単純なものではない。上述の読影システムは、輝点が腫瘤の中心をピークとする「正規分布」のような、極めて典型的な腫瘤の形状についてはよく機能するが、現時点では、総合的には実用レベルにはほど遠く、あくまで診断の補助的なツールとして位置づけられるにすぎない。
【0006】
すなわち、マンモグラフィによる読影は、十分な経験を積んだ医師であっても判断に迷うものが、少なくない。日本医学放射線学会では、乳がんの発見を目的とするマンモグラフィー診断の読影法などについて、ガイドラインを策定して標準化を図っているが、これも腫瘤の判定がいかに難しいかを示すものに他ならない。
【0007】
特に、乳腺のマンモグラフィー診断における腫瘤の発見は、単なる乳腺の重なりであるか真の腫瘍であるかの区別が非常に難しく、腫瘤のカテゴリーの判断樹に従って、形状、境界及び辺縁、画像上の相対的な輝度などを総合的に判断し、他の診断方法による結果も考慮して最終的な診断を下しているのが実情である。
【0008】
本発明は、上記に鑑みてなされたものであり、デジタル画像における特徴部分を強調して表示させることを可能にする、新規な画像処理方法及び画像処理プログラムを提供することを主たる技術的課題とする。
【課題を解決するための手段】
【0009】
本発明に係る画像処理方法は、大きく分けて3つの新しい処理方法を提供する。
本発明に係る第1及び第2の発明は、コンピュータにより入力画像に含まれる各画素の画像信号の大きさを所定の目的に従って調整する画像処理方法であり、第3の発明は、コンピュータにより入力画像に含まれる各画素の画像信号の特徴部分を抽出し、これを閉曲線で囲うための画像処理方法である。
【0010】
第1の発明は、入力画像のデータをコンピュータに入力するステップ(Sa1)と、コントラスト強度として任意の定数mを入力するステップ(Sa2)と、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Dとして区画し、演算対象領域D内における注目画素とその周囲の任意の画素との2点の画素の傾きの強さを前記注目画素を中心とする全周に亘って合計し演算対象領域内の画素数で除算することによりその平均値を算出するステップ(Sa3)と、ステップ(Sa3)で求めた勾配の振幅値βにステップ(Sa2)で入力したコントラスト強度m(但し、mは正の定数である。)を乗算し、任意のオフセット値γを加算するステップ(Sa4)とを具備することを特徴とする。
【0011】
本発明では、ステップ(Sa3)において、勾配の振幅値βが注目画素と周辺画素との傾きの強さを注目画素を中心とする全周に亘って合計し演算対象領域内の画素数で除算することによりその平均値を算出するため、勾配の角度成分が打ち消され、角度情報を含まない、振幅値のみの差分画像が生成される。しかも、振幅値はコントラスト強度mを乗じることで一層強調されるため、注目画素とその周辺画素との画像信号の差が大きければ大きいほど強く表されることになる。
【0012】
画像信号とは輝度値や色相(例えば、赤・緑・青)の階調を表すスカラー量であるから、画像信号を輝度値として、本発明を適用すると、出力画像は、特に、入力画像中の線状のパターンや輪郭などの一次元的な模様が強調されて表示されることになる。
【0013】
第1の発明では、演算対象領域Dを注目画素の周囲5画素以上とすることにより、より広範囲の周辺画素の画像信号を考慮した振幅値を計算することが可能となる。
【0014】
第2の発明は、入力画像のデータをコンピュータに入力するステップ(Sb1)と、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から求めた勾配ベクトルの角度成分θを算出するステップ(Sb2)と、勾配ベクトルの角度成分θの絶対値に所定の定数を乗算することにより、演算対象領域C内の中心画像信号αを算出するステップ(Sb3)とを具備することを特徴とする。
【0015】
本発明では、ステップ(Sb3)において、中心画像信号αがθの関数すなわち角度依存性を持つものとして表される。このため、画像信号を輝度値とすると、出力信号は輝度が角度依存性を持って表されるため、入力画像中に輝度の勾配がある部分が陰影を伴って立体的に表現され、入力画像の特徴部分が強調されて表示される。
【0016】
第3の発明は、入力画像のデータをコンピュータに入力するステップ(Sc1)と、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から勾配ベクトルの大きさと演算対象領域Cの走査方向に向かう単位ベクトルと勾配ベクトルとのなす角θとを算出するステップ(Sc2)と、前記勾配ベクトルと前記単位ベクトルの内積を求めて規格化することにより勾配ベクトルの集中度Zを算出するステップ(Sc3)と、ラプラシアンフィルタ、ラプラシアンオブガウシアン(LOG)フィルタその他の2階微分に基づくフィルタを乗算するステップ(Sc4)と、前記ステップ(Sc4)に対して2値化処理を行なって出力画像を得るステップ(Sc5)とを具備することを特徴とする。
【0017】
第3の発明は、公知のアイリスフィルタによるアルゴリズム(Sc2〜Sc3)を用いて集中度を算出した後、ラプラシアンフィルタ、ラプラシアンオブガウシアン(LOG)フィルタその他の2階微分に基づくフィルタによってエッジ検出を行うことにより集中度の高い部分の輪郭が閉曲線で囲われる。
【0018】
第3の発明では、アイリスフィルタによる集中度の演算ステップ((ステップSc2〜Sc3)を2回以上繰り返すことにより、画像信号の値が特に集中している特徴点を抽出するようにしてもよい。
【0019】
本発明に係る第1〜第3の画像処理方法は、パーソナルコンピュータ等のコンピュータ上で実行される画像処理用アプリケーションプログラムとして実現することができる。
【発明の効果】
【0020】
本発明に係る画像処理方法によると、入力画像の特徴を抽出したり、際だって表示させることができる。また、抽出結果を公知の閾値処理によって除去すれば特定の特徴を有するパターンを入力画像から除去することも可能である。
【図面の簡単な説明】
【0021】
【図1】本発明の第1の実施形態に係る画像処理方法の手順を示す図
【図2】第2の実施形態に係る画像処理方法の演算対象領域D
【図3】(a)第1の実施形態に係るマンモグラフィの原画像である入力画像の例(b)第1の実施形態に係る出力画像の例(c)第1の実施形態に係る出力画像の他の例
【図4】本発明の第2の実施形態に係る画像処理方法の手順を示す図
【図5】(a)勾配ベクトルを求めるための演算対象領域Cの例(b)(a)において、x方向の増加量=1、y方向の増加量が−1であった場合の勾配ベクトルの例
【図6】(a)第1の実施形態に係る入力画像の例(b)第1の実施形態に係る出力画像の例(c)第1の実施形態に係る出力画像の他の例
【図7】(a)第2の実施形態に係るマンモグラフィの原画像である入力画像の例(b)第2の実施形態に係る出力画像の例(c)第2の実施形態に係る出力画像の他の例
【図8】本発明の第3の実施形態に係る画像処理方法の手順を示す図
【図9】ラプラシアンフィルタの例(8近傍ラプラシアンフィルタ)
【図10】(a)第3の実施形態に係るマンモグラフィの原画像である入力画像の例(b)第3の実施形態に係る出力画像の例(c)第3の実施形態に係る出力画像の他の例
【図11】(a)ラプラシアンフィルタによるエッジ検出を行った出力画像の例(b)(a)の模式図
【図12】本発明の第4の実施形態における画像処理装置の構成を示すブロック図
【発明を実施するための形態】
【0022】
(第1の実施形態)−微分画像の振幅値による特徴点の抽出処理について−
図1は、本発明の第1の実施形態に係る画像処理方法の手順を示す図である。
入力画像のデータを前記コンピュータに入力するステップ(Sb1)に続いて、コントラスト強度として任意の定数mを入力するステップ(Sb2)を実行する。この定数mは、出力画像に付与するコントラストの強度を規定する正の定数である。
【0023】
図2は、入力画像中の複数の画素からなる閉領域(x,y)を区画した演算対象領域Dを表している。この例では、この領域の中心(注目画素)から周囲5画素を、演算対象領域としている。注目画素x(i=0,j=0)としたとき、その振幅値βを式1のように定義する。
【0024】
(式1)
【0025】
一般に、微分画像の振幅値は、x成分とy成分とに分離して求められる。このような振幅値の計算方法によると、注目画素と周囲の画素との勾配ベクトルの角度成分を算出し、勾配の方向を設定しなければならない。しかし、振幅と角度を同時に画像表示することは困難を伴う。
【0026】
そこで、第1の実施形態に係る発明では、演算対象領域D内における注目画素とその周囲の任意の画素との2点の画素における画像信号の勾配(すなわち差分)の振幅値βを勾配の振幅値βが注目画素と周辺画素との傾きの強さを注目画素を中心とする全周に亘って合計し演算対象領域内の画素数で除算することによりその平均値を算出することにより求めている。
【0027】
このような手法により、勾配の方向成分が打ち消されるため、注目画素の振幅値のみが数値化され、角度情報のない振幅値のみの勾配が生成されることになる。また、演算対象領域が広ければ広いほど、中心画素からより遠い画素を含むことになり、より広い範囲を考慮した振幅値を計算することが可能になる。
【0028】
画像信号が輝度値の場合、振幅値βは注目画素の輝度値そのものを表す。このβの値が正の時は注目画素が周辺より明るく表示され、逆にβの値が負の時は注目画素が周辺より暗く表示されることを表している。
【0029】
そして、この演算対象領域Dを前記入力画像内で所定の画素単位で移動させることにより、演算対象領域内における注目画素とその周囲の任意の画素との2点の画素における画像信号の振幅値βを順次算出していく。
【0030】
次のステップ(Sb4)では、ステップ(Sb3)で求めた振幅値βにステップSb2で入力したコントラスト強度mを乗算し、任意のオフセット値を加算する。このステップは演算対象領域Dを入力画像内でx及びy方向に移動させることによりステップ(Sb3)〜ステップ(Sb4)を繰り返したあと(繰り返し(ステップ(Sb5))、最後に入力画像データに各画素に対して計算した振幅値βを入力画像のデータに適応させることにより出力画像を得る(ステップ(Sb6))。
【0031】
図3(a)は、マンモグラフィの原画像である入力画像を表している。図3(b)は、図3(a)の画像に、式1に従って算出した振幅値βが正の値をもつ画素のみを残し、その他は輝度値を0にクリッピングした後、振幅値βの値を3倍した出力画像を表している。
【0032】
図3(c)は、図3(a)の画像に、式1に従って算出した振幅値βを3倍し、オフセットとして全画素に127を加算した出力画像を表している。
【0033】
図3(b)、(c)では、原画像では目立たないスピキュラ、乳腺、血管が浮かび上がって表示されている。このように、振幅値βを計算し、この値を一次関数による線形変換を行うことで入力画像では目立たない、線状のパターンや輪郭などの1次元的な模様が際だって表示される。
【0034】
(第2の実施形態)−勾配ベクトルを用いた立体視化処理について−
図4は、本発明の第2の実施形態に係る画像処理方法の手順を示す図である。「勾配ベクトル」とは、一般に、角度情報を持つ微分画像を表すためのベクトル量と定義される。以下の説明では、画像信号は輝度を表し、画像信号の値は輝度値を表すものとする。また、説明の便宜上、輝度値は8ビットで表し、すなわち0〜255までの値をとるものとする。
【0035】
図5(a)は、入力画像に設けられた演算対象領域(マスク)Cを例示している。演算対象領域Cは一例として中心画素から周囲1画素分の領域としている。このとき、演算対象領域Cの中心画素の輝度値をα、中心画素の周囲の輝度値がそれぞれf1〜f8と表される。なお、中心画素から周囲何画素分の領域を演算対象領域としてもよい。
【0036】
勾配ベクトルの角度成分(傾きθ)は一般にx方向(右方向)の増加量に対するy方向(上方向)への増加量のアークタンジェントで表される。例えば、図5(a)の例では、式2のように記述される。
【0037】
(式2)
【0038】
図5(b)は、x方向の増加量=1、y方向の増加量が−1であった場合の勾配ベクトルを例示している。θ=tan−1(−1)=315度となる。もちろん、θの取り方によって、−45度と記載してもよい。
【0039】
図6(a)は、入力画像の一部を例示したものである。この円形の入力画像は、中心部ほど輝度値が高く、周辺部ほど低い分布を持っている。この表示方法によると、平面上のどの方向から見ても中心部に近づくほど一様に輝度値が高くなるため、平面視図となる。
【0040】
図6(b)は、式3に従って、0度から360度の角度情報を0〜255の輝度値に線形変換し、画像表示した様子を示している。
【0041】
(式3)
【0042】
このように、勾配ベクトルから求めた角度情報を輝度値に線形変換することで角度依存性をもった立体的形状の図形が得られる。ただし、図6(b)に示す画像は、θの定義域を0度から360度までとしたため、輝度変化の断層面が形成された。
【0043】
そこで、図6(c)では、勾配ベクトルが0度(360度)を示すところを最大輝度値255とし、そこから±180度に向かうに従って輝度値を小さくしていき、180度で中心輝度値αの値が最小となるように角度依存情報を付加した。すなわち、中心輝度値αの大きさとθとの関係は、式4のように表される。
【0044】
(式4)
【0045】
図6(c)は図の矢印Aの方向から光を当てたときの立体視をイメージさせる印影が現れている。なお、図6(b)は現実世界のイメージとは一致しないものの、入力画像における勾配をもつ輝度集中点などの特徴点がよく現れるという意味において、図6(c)同様に大きな効果があるといえる。
【0046】
また、角度依存情報を付与した別の効果として、輪郭が鮮明に表示されるという点が挙げられる。図6(a)では、輪郭がぼんやりとしているのに対し、図6(b)、(c)では輪郭が鮮明に現れている。
【0047】
次に、これをマンモグラフィ画像に適応した例を示す。
図7(a)は、マンモグラフィの原画像である入力画像を表している。
【0048】
図7(b)、(c)は、図7(a)の画像に、それぞれ式3、式4に従って角度依存性を付与した出力画像である。図7(b)と(c)とでは、θの定義域が異なるがいずれも角度依存情報が付与されている点で共通する。
【0049】
図7(a)中の輝度集中点はスピキュラを伴う腫瘤である。図7(b)、(c)では、この輝度集中点が立体的に表示され、特徴点がより鮮明に現れた。
【0050】
また、特筆すべき他の点として、乳腺と輪郭が原画像よりも鮮明に現れたことが挙げられる。
【0051】
このように、平面的な入力画像において輝度勾配がある部分が立体的に表現され、入力画像の特徴部分が際だって表示される。
【0052】
(第3の実施形態)−ラプラシアンフィルタによるエッジ検出処理について−
この実施形態では、従来から知られるアイリスフィルタの出力画像に対して、ラプラシアンフィルタを適応することによりエッジ検出を行う方法について説明する。エッジ検出は、コンピュータにより入力画像に含まれる各画素の画像信号の特徴部分を抽出し、これを閉曲線で囲っていわゆる「マーキング」を行うことを目的とするものである。
【0053】
図8は、本発明の第3の実施形態に係る画像処理方法の手順を示す図である。ステップ(Sc2)〜ステップ(Sc3)は、輝度値などの画像信号の集中度を求めるための公知のアルゴリズムであり、「アイリスフィルタ」と呼ばれている。
【0054】
まず、入力画像のデータをコンピュータに入力し、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域Cを前記入力画像内で所定の画素単位で移動させることにより、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から求めた勾配ベクトルと演算対象領域Cの走査方向に向かう単位ベクトルと前記勾配ベクトルと前記単位ベクトルのなす各θとから、両ベクトルの内積を求めて規格化した値が勾配ベクトルの集中度Zとなる。
【0055】
次のステップ(Sc4)は、ラプラシアンフィルタを乗算するアルゴリズムである。
図9は、ラプラシアンフィルタの一例を示している。この図に示すように、ラプラシアンフィルタのオペレータは、3×3のマスクで中心画素の値が−8、その周囲が全て1であるいわゆる「8近傍ラプラシアンフィルタ」を用いることができる。
【0056】
図9は、8近傍ラプラシアンフィルタを表している。このように、上下左右の画素の4方向と、斜め方向の画素の4方向を加えた8方向とで処理を行う。このラプラシアンフィルタを乗算したのち、さらに2値化処理を行ってエッジ部分を強調する。
【0057】
ステップ(Sc5)は前記ステップ(Sc4)に対して2値化処理を行なって出力画像を得るステップである。集中度の計算値が所定の値以下の画素では全ての画像信号を0にクリッピングし、それ以外の画素ではそのまま若しくはより高い値に線形変換する。
【0058】
図10(a)は、マンモグラフィの原画像である入力画像を表している。図10(b)及び(c)は、図10(a)の画像に、上述のアルゴリズムによるアイリスフィルタ処理(ステップSc2〜Sc3)を行った出力画像を閾値処理することにより、所定の輝度値を持つ画素のみを抽出した画像を示している。図10(b)では、集中度が0以下の画素は0にクリッピングし、その後全ての画素における輝度値を320倍した出力画像を示している。また、図10(c)では、図10(b)を入力画像として、輝度値が100以下の画素は0にクリッピングすることにより、集中度の高い領域のみを表示した出力画像を示している。図10(c)で残った部分は腫瘤陰影である。但し、実施には腫瘤陰影のみが残るが腫瘤の位置を判りやすくするため、図10(c)では乳房の輪郭を点線で表示している。
【0059】
この腫瘤陰影をマーキングするために、出力画像2に対して、ラプラシアンフィルタによるエッジ検出を行う。
【0060】
このラプラシアンオペレータを図9の(b)、(c)に対して行ったあと、2値化処理を行ってエッジ部分を強調すると、図11における腫瘤陰影の周囲にマーキングされた図11(a)のような出力画像が得られる。図11(b)は、マーキング部と乳房の外形を実線又は破線で示したものである。
【0061】
図11(b)に示される外側のマーク101は、入力画像にアイリスフィルタを1回かけた出力画像(図9(b))にラプラシアンフィルタのオペレータ処理を行った出力画像を表し、内側のマーク102は、入力画像にアイリスフィルタを2回かけた出力画像(図9(c))にラプラシアンフィルタのオペレータ処理を行った出力画像を表している。
【0062】
外側のマーク101は腫瘤陰影の候補をもれなく拾い上げるために表されたものであり、内側のマーク102は拾い上げた腫瘤陰影の中で、より腫瘤である危険性が高いものを示すものである。
【0063】
(第4の実施形態)
第4の実施形態は、パーソナルコンピュータ等のコンピュータ上で実行される画像処理用アプリケーションプログラムとして、本発明の画像処理方法を実現するコンピュータプログラムを提供する形態を示している。
【0064】
図12は、本発明の第4の実施形態における画像処理装置の構成を示すブロック図である。パーソナルコンピュータに本発明に係る各画像処理方法の各ステップを実行させるプログラムをインストールした画像処理装置20は、装置全体を制御するCPU(Central Processing Unit)21、RAM(Random Access Memory)等の一時記憶手段22、ハードディスクなどの記憶手段23、記録媒体から各データを読み取る光学ディスクドライブ等の補助記憶手段24、マウス及びキーボード等の入力手段25、モニタやプリンタ等の出力手段26を備えている。記録媒体に記録された本発明に係るコンピュータプログラムを記憶手段23に読み取り、制御手段21の制御により、オペレーティングシステム上で実行させることにより、パーソナルコンピュータは本発明に係る画像処理装置20として動作する。
【0065】
画像処理装置20では、記憶手段22に記憶させたコンピュータプログラムを制御手段21の制御に基づいて第1乃至第3の実施形態で説明した各ステップを実行することにより、記憶手段23などに記憶している画像や外部とのインターフェース27を通じて入力を受け付けた画像を入力画像に対し、画像処理を行い、出力手段26から出力画像を出力する。
【産業上の利用可能性】
【0066】
本発明に係る画像処理方法は、マンモグラフィなど医療診断の支援システムを始め、目立たない部分や特徴部分をより見やすく表示するといった用途への応用が期待され、従って、本発明の産業上利用可能性は極めて大きい。
【符号の説明】
【0067】
20 本発明のプログラムをインストールした画像処理装置
21 CPU(Central Processing Unit)
22 一時記憶手段(RAM)
23 記憶手段(ハードディスク)
24 補助記憶手段(光学ディスクドライブ等)
25 入力手段(マウス及びキーボード等)
26 出力手段(モニタやプリンタ等)
27 入出力インターフェース
【技術分野】
【0001】
本発明は、デジタル画像に含まれる特定の部分を抽出し、より強調して表示させたり、逆に、特定の部分を除去したりすることを可能にするデジタル画像の画像処理方法及びコンピュータプログラムに関するものである。
【背景技術】
【0002】
デジタル画像において、注目画素とその周囲の画素とを比較した際の両者のデジタル化された画像信号(例えば輝度値)の変化に関する情報を表すベクトルは一般に「微分画像」又は「勾配画像」と呼ばれている(本明細書では、「微分画像」として統一する。)。微分画像は、注目画素とその周辺の任意の点の2点(2画素)を比較して、座標及びデジタル値のそれぞれの変化から算出された「勾配ベクトル」を表している。勾配ベクトルを利用した画像処理アルゴリズムの一つとして、「アイリスフィルタ」が知られている。アイリスフィルタでは、その勾配ベクトルを求めたあと、そこからデジタル化された画像信号(例えば輝度値)の「集中度」を計算する。この集中度やその他の計算値を閾値処理することにより、マンモグラフィ診断画像の読影を支援するシステムに関するものが開示されている(例えば、特許文献1〜4等参照)
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特開平09−167238号公報
【特許文献2】特開平10−289318号公報
【特許文献3】特開2005−52295号公報
【特許文献4】特開2006−130049号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
この種の読影システムは、腫瘤陰影の輝度値が放射状に減少する(すなわち、腫瘤陰影は中心ほど明るく周辺ほど暗いという特徴を有する)というモデルに基づいて設計されている。
【0005】
しかし、実際の腫瘤陰影には様々な形状があり、必ずしも上述のようなモデルで全てが解決できるほど単純なものではない。上述の読影システムは、輝点が腫瘤の中心をピークとする「正規分布」のような、極めて典型的な腫瘤の形状についてはよく機能するが、現時点では、総合的には実用レベルにはほど遠く、あくまで診断の補助的なツールとして位置づけられるにすぎない。
【0006】
すなわち、マンモグラフィによる読影は、十分な経験を積んだ医師であっても判断に迷うものが、少なくない。日本医学放射線学会では、乳がんの発見を目的とするマンモグラフィー診断の読影法などについて、ガイドラインを策定して標準化を図っているが、これも腫瘤の判定がいかに難しいかを示すものに他ならない。
【0007】
特に、乳腺のマンモグラフィー診断における腫瘤の発見は、単なる乳腺の重なりであるか真の腫瘍であるかの区別が非常に難しく、腫瘤のカテゴリーの判断樹に従って、形状、境界及び辺縁、画像上の相対的な輝度などを総合的に判断し、他の診断方法による結果も考慮して最終的な診断を下しているのが実情である。
【0008】
本発明は、上記に鑑みてなされたものであり、デジタル画像における特徴部分を強調して表示させることを可能にする、新規な画像処理方法及び画像処理プログラムを提供することを主たる技術的課題とする。
【課題を解決するための手段】
【0009】
本発明に係る画像処理方法は、大きく分けて3つの新しい処理方法を提供する。
本発明に係る第1及び第2の発明は、コンピュータにより入力画像に含まれる各画素の画像信号の大きさを所定の目的に従って調整する画像処理方法であり、第3の発明は、コンピュータにより入力画像に含まれる各画素の画像信号の特徴部分を抽出し、これを閉曲線で囲うための画像処理方法である。
【0010】
第1の発明は、入力画像のデータをコンピュータに入力するステップ(Sa1)と、コントラスト強度として任意の定数mを入力するステップ(Sa2)と、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Dとして区画し、演算対象領域D内における注目画素とその周囲の任意の画素との2点の画素の傾きの強さを前記注目画素を中心とする全周に亘って合計し演算対象領域内の画素数で除算することによりその平均値を算出するステップ(Sa3)と、ステップ(Sa3)で求めた勾配の振幅値βにステップ(Sa2)で入力したコントラスト強度m(但し、mは正の定数である。)を乗算し、任意のオフセット値γを加算するステップ(Sa4)とを具備することを特徴とする。
【0011】
本発明では、ステップ(Sa3)において、勾配の振幅値βが注目画素と周辺画素との傾きの強さを注目画素を中心とする全周に亘って合計し演算対象領域内の画素数で除算することによりその平均値を算出するため、勾配の角度成分が打ち消され、角度情報を含まない、振幅値のみの差分画像が生成される。しかも、振幅値はコントラスト強度mを乗じることで一層強調されるため、注目画素とその周辺画素との画像信号の差が大きければ大きいほど強く表されることになる。
【0012】
画像信号とは輝度値や色相(例えば、赤・緑・青)の階調を表すスカラー量であるから、画像信号を輝度値として、本発明を適用すると、出力画像は、特に、入力画像中の線状のパターンや輪郭などの一次元的な模様が強調されて表示されることになる。
【0013】
第1の発明では、演算対象領域Dを注目画素の周囲5画素以上とすることにより、より広範囲の周辺画素の画像信号を考慮した振幅値を計算することが可能となる。
【0014】
第2の発明は、入力画像のデータをコンピュータに入力するステップ(Sb1)と、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から求めた勾配ベクトルの角度成分θを算出するステップ(Sb2)と、勾配ベクトルの角度成分θの絶対値に所定の定数を乗算することにより、演算対象領域C内の中心画像信号αを算出するステップ(Sb3)とを具備することを特徴とする。
【0015】
本発明では、ステップ(Sb3)において、中心画像信号αがθの関数すなわち角度依存性を持つものとして表される。このため、画像信号を輝度値とすると、出力信号は輝度が角度依存性を持って表されるため、入力画像中に輝度の勾配がある部分が陰影を伴って立体的に表現され、入力画像の特徴部分が強調されて表示される。
【0016】
第3の発明は、入力画像のデータをコンピュータに入力するステップ(Sc1)と、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から勾配ベクトルの大きさと演算対象領域Cの走査方向に向かう単位ベクトルと勾配ベクトルとのなす角θとを算出するステップ(Sc2)と、前記勾配ベクトルと前記単位ベクトルの内積を求めて規格化することにより勾配ベクトルの集中度Zを算出するステップ(Sc3)と、ラプラシアンフィルタ、ラプラシアンオブガウシアン(LOG)フィルタその他の2階微分に基づくフィルタを乗算するステップ(Sc4)と、前記ステップ(Sc4)に対して2値化処理を行なって出力画像を得るステップ(Sc5)とを具備することを特徴とする。
【0017】
第3の発明は、公知のアイリスフィルタによるアルゴリズム(Sc2〜Sc3)を用いて集中度を算出した後、ラプラシアンフィルタ、ラプラシアンオブガウシアン(LOG)フィルタその他の2階微分に基づくフィルタによってエッジ検出を行うことにより集中度の高い部分の輪郭が閉曲線で囲われる。
【0018】
第3の発明では、アイリスフィルタによる集中度の演算ステップ((ステップSc2〜Sc3)を2回以上繰り返すことにより、画像信号の値が特に集中している特徴点を抽出するようにしてもよい。
【0019】
本発明に係る第1〜第3の画像処理方法は、パーソナルコンピュータ等のコンピュータ上で実行される画像処理用アプリケーションプログラムとして実現することができる。
【発明の効果】
【0020】
本発明に係る画像処理方法によると、入力画像の特徴を抽出したり、際だって表示させることができる。また、抽出結果を公知の閾値処理によって除去すれば特定の特徴を有するパターンを入力画像から除去することも可能である。
【図面の簡単な説明】
【0021】
【図1】本発明の第1の実施形態に係る画像処理方法の手順を示す図
【図2】第2の実施形態に係る画像処理方法の演算対象領域D
【図3】(a)第1の実施形態に係るマンモグラフィの原画像である入力画像の例(b)第1の実施形態に係る出力画像の例(c)第1の実施形態に係る出力画像の他の例
【図4】本発明の第2の実施形態に係る画像処理方法の手順を示す図
【図5】(a)勾配ベクトルを求めるための演算対象領域Cの例(b)(a)において、x方向の増加量=1、y方向の増加量が−1であった場合の勾配ベクトルの例
【図6】(a)第1の実施形態に係る入力画像の例(b)第1の実施形態に係る出力画像の例(c)第1の実施形態に係る出力画像の他の例
【図7】(a)第2の実施形態に係るマンモグラフィの原画像である入力画像の例(b)第2の実施形態に係る出力画像の例(c)第2の実施形態に係る出力画像の他の例
【図8】本発明の第3の実施形態に係る画像処理方法の手順を示す図
【図9】ラプラシアンフィルタの例(8近傍ラプラシアンフィルタ)
【図10】(a)第3の実施形態に係るマンモグラフィの原画像である入力画像の例(b)第3の実施形態に係る出力画像の例(c)第3の実施形態に係る出力画像の他の例
【図11】(a)ラプラシアンフィルタによるエッジ検出を行った出力画像の例(b)(a)の模式図
【図12】本発明の第4の実施形態における画像処理装置の構成を示すブロック図
【発明を実施するための形態】
【0022】
(第1の実施形態)−微分画像の振幅値による特徴点の抽出処理について−
図1は、本発明の第1の実施形態に係る画像処理方法の手順を示す図である。
入力画像のデータを前記コンピュータに入力するステップ(Sb1)に続いて、コントラスト強度として任意の定数mを入力するステップ(Sb2)を実行する。この定数mは、出力画像に付与するコントラストの強度を規定する正の定数である。
【0023】
図2は、入力画像中の複数の画素からなる閉領域(x,y)を区画した演算対象領域Dを表している。この例では、この領域の中心(注目画素)から周囲5画素を、演算対象領域としている。注目画素x(i=0,j=0)としたとき、その振幅値βを式1のように定義する。
【0024】
(式1)
【0025】
一般に、微分画像の振幅値は、x成分とy成分とに分離して求められる。このような振幅値の計算方法によると、注目画素と周囲の画素との勾配ベクトルの角度成分を算出し、勾配の方向を設定しなければならない。しかし、振幅と角度を同時に画像表示することは困難を伴う。
【0026】
そこで、第1の実施形態に係る発明では、演算対象領域D内における注目画素とその周囲の任意の画素との2点の画素における画像信号の勾配(すなわち差分)の振幅値βを勾配の振幅値βが注目画素と周辺画素との傾きの強さを注目画素を中心とする全周に亘って合計し演算対象領域内の画素数で除算することによりその平均値を算出することにより求めている。
【0027】
このような手法により、勾配の方向成分が打ち消されるため、注目画素の振幅値のみが数値化され、角度情報のない振幅値のみの勾配が生成されることになる。また、演算対象領域が広ければ広いほど、中心画素からより遠い画素を含むことになり、より広い範囲を考慮した振幅値を計算することが可能になる。
【0028】
画像信号が輝度値の場合、振幅値βは注目画素の輝度値そのものを表す。このβの値が正の時は注目画素が周辺より明るく表示され、逆にβの値が負の時は注目画素が周辺より暗く表示されることを表している。
【0029】
そして、この演算対象領域Dを前記入力画像内で所定の画素単位で移動させることにより、演算対象領域内における注目画素とその周囲の任意の画素との2点の画素における画像信号の振幅値βを順次算出していく。
【0030】
次のステップ(Sb4)では、ステップ(Sb3)で求めた振幅値βにステップSb2で入力したコントラスト強度mを乗算し、任意のオフセット値を加算する。このステップは演算対象領域Dを入力画像内でx及びy方向に移動させることによりステップ(Sb3)〜ステップ(Sb4)を繰り返したあと(繰り返し(ステップ(Sb5))、最後に入力画像データに各画素に対して計算した振幅値βを入力画像のデータに適応させることにより出力画像を得る(ステップ(Sb6))。
【0031】
図3(a)は、マンモグラフィの原画像である入力画像を表している。図3(b)は、図3(a)の画像に、式1に従って算出した振幅値βが正の値をもつ画素のみを残し、その他は輝度値を0にクリッピングした後、振幅値βの値を3倍した出力画像を表している。
【0032】
図3(c)は、図3(a)の画像に、式1に従って算出した振幅値βを3倍し、オフセットとして全画素に127を加算した出力画像を表している。
【0033】
図3(b)、(c)では、原画像では目立たないスピキュラ、乳腺、血管が浮かび上がって表示されている。このように、振幅値βを計算し、この値を一次関数による線形変換を行うことで入力画像では目立たない、線状のパターンや輪郭などの1次元的な模様が際だって表示される。
【0034】
(第2の実施形態)−勾配ベクトルを用いた立体視化処理について−
図4は、本発明の第2の実施形態に係る画像処理方法の手順を示す図である。「勾配ベクトル」とは、一般に、角度情報を持つ微分画像を表すためのベクトル量と定義される。以下の説明では、画像信号は輝度を表し、画像信号の値は輝度値を表すものとする。また、説明の便宜上、輝度値は8ビットで表し、すなわち0〜255までの値をとるものとする。
【0035】
図5(a)は、入力画像に設けられた演算対象領域(マスク)Cを例示している。演算対象領域Cは一例として中心画素から周囲1画素分の領域としている。このとき、演算対象領域Cの中心画素の輝度値をα、中心画素の周囲の輝度値がそれぞれf1〜f8と表される。なお、中心画素から周囲何画素分の領域を演算対象領域としてもよい。
【0036】
勾配ベクトルの角度成分(傾きθ)は一般にx方向(右方向)の増加量に対するy方向(上方向)への増加量のアークタンジェントで表される。例えば、図5(a)の例では、式2のように記述される。
【0037】
(式2)
【0038】
図5(b)は、x方向の増加量=1、y方向の増加量が−1であった場合の勾配ベクトルを例示している。θ=tan−1(−1)=315度となる。もちろん、θの取り方によって、−45度と記載してもよい。
【0039】
図6(a)は、入力画像の一部を例示したものである。この円形の入力画像は、中心部ほど輝度値が高く、周辺部ほど低い分布を持っている。この表示方法によると、平面上のどの方向から見ても中心部に近づくほど一様に輝度値が高くなるため、平面視図となる。
【0040】
図6(b)は、式3に従って、0度から360度の角度情報を0〜255の輝度値に線形変換し、画像表示した様子を示している。
【0041】
(式3)
【0042】
このように、勾配ベクトルから求めた角度情報を輝度値に線形変換することで角度依存性をもった立体的形状の図形が得られる。ただし、図6(b)に示す画像は、θの定義域を0度から360度までとしたため、輝度変化の断層面が形成された。
【0043】
そこで、図6(c)では、勾配ベクトルが0度(360度)を示すところを最大輝度値255とし、そこから±180度に向かうに従って輝度値を小さくしていき、180度で中心輝度値αの値が最小となるように角度依存情報を付加した。すなわち、中心輝度値αの大きさとθとの関係は、式4のように表される。
【0044】
(式4)
【0045】
図6(c)は図の矢印Aの方向から光を当てたときの立体視をイメージさせる印影が現れている。なお、図6(b)は現実世界のイメージとは一致しないものの、入力画像における勾配をもつ輝度集中点などの特徴点がよく現れるという意味において、図6(c)同様に大きな効果があるといえる。
【0046】
また、角度依存情報を付与した別の効果として、輪郭が鮮明に表示されるという点が挙げられる。図6(a)では、輪郭がぼんやりとしているのに対し、図6(b)、(c)では輪郭が鮮明に現れている。
【0047】
次に、これをマンモグラフィ画像に適応した例を示す。
図7(a)は、マンモグラフィの原画像である入力画像を表している。
【0048】
図7(b)、(c)は、図7(a)の画像に、それぞれ式3、式4に従って角度依存性を付与した出力画像である。図7(b)と(c)とでは、θの定義域が異なるがいずれも角度依存情報が付与されている点で共通する。
【0049】
図7(a)中の輝度集中点はスピキュラを伴う腫瘤である。図7(b)、(c)では、この輝度集中点が立体的に表示され、特徴点がより鮮明に現れた。
【0050】
また、特筆すべき他の点として、乳腺と輪郭が原画像よりも鮮明に現れたことが挙げられる。
【0051】
このように、平面的な入力画像において輝度勾配がある部分が立体的に表現され、入力画像の特徴部分が際だって表示される。
【0052】
(第3の実施形態)−ラプラシアンフィルタによるエッジ検出処理について−
この実施形態では、従来から知られるアイリスフィルタの出力画像に対して、ラプラシアンフィルタを適応することによりエッジ検出を行う方法について説明する。エッジ検出は、コンピュータにより入力画像に含まれる各画素の画像信号の特徴部分を抽出し、これを閉曲線で囲っていわゆる「マーキング」を行うことを目的とするものである。
【0053】
図8は、本発明の第3の実施形態に係る画像処理方法の手順を示す図である。ステップ(Sc2)〜ステップ(Sc3)は、輝度値などの画像信号の集中度を求めるための公知のアルゴリズムであり、「アイリスフィルタ」と呼ばれている。
【0054】
まず、入力画像のデータをコンピュータに入力し、入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域Cを前記入力画像内で所定の画素単位で移動させることにより、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から求めた勾配ベクトルと演算対象領域Cの走査方向に向かう単位ベクトルと前記勾配ベクトルと前記単位ベクトルのなす各θとから、両ベクトルの内積を求めて規格化した値が勾配ベクトルの集中度Zとなる。
【0055】
次のステップ(Sc4)は、ラプラシアンフィルタを乗算するアルゴリズムである。
図9は、ラプラシアンフィルタの一例を示している。この図に示すように、ラプラシアンフィルタのオペレータは、3×3のマスクで中心画素の値が−8、その周囲が全て1であるいわゆる「8近傍ラプラシアンフィルタ」を用いることができる。
【0056】
図9は、8近傍ラプラシアンフィルタを表している。このように、上下左右の画素の4方向と、斜め方向の画素の4方向を加えた8方向とで処理を行う。このラプラシアンフィルタを乗算したのち、さらに2値化処理を行ってエッジ部分を強調する。
【0057】
ステップ(Sc5)は前記ステップ(Sc4)に対して2値化処理を行なって出力画像を得るステップである。集中度の計算値が所定の値以下の画素では全ての画像信号を0にクリッピングし、それ以外の画素ではそのまま若しくはより高い値に線形変換する。
【0058】
図10(a)は、マンモグラフィの原画像である入力画像を表している。図10(b)及び(c)は、図10(a)の画像に、上述のアルゴリズムによるアイリスフィルタ処理(ステップSc2〜Sc3)を行った出力画像を閾値処理することにより、所定の輝度値を持つ画素のみを抽出した画像を示している。図10(b)では、集中度が0以下の画素は0にクリッピングし、その後全ての画素における輝度値を320倍した出力画像を示している。また、図10(c)では、図10(b)を入力画像として、輝度値が100以下の画素は0にクリッピングすることにより、集中度の高い領域のみを表示した出力画像を示している。図10(c)で残った部分は腫瘤陰影である。但し、実施には腫瘤陰影のみが残るが腫瘤の位置を判りやすくするため、図10(c)では乳房の輪郭を点線で表示している。
【0059】
この腫瘤陰影をマーキングするために、出力画像2に対して、ラプラシアンフィルタによるエッジ検出を行う。
【0060】
このラプラシアンオペレータを図9の(b)、(c)に対して行ったあと、2値化処理を行ってエッジ部分を強調すると、図11における腫瘤陰影の周囲にマーキングされた図11(a)のような出力画像が得られる。図11(b)は、マーキング部と乳房の外形を実線又は破線で示したものである。
【0061】
図11(b)に示される外側のマーク101は、入力画像にアイリスフィルタを1回かけた出力画像(図9(b))にラプラシアンフィルタのオペレータ処理を行った出力画像を表し、内側のマーク102は、入力画像にアイリスフィルタを2回かけた出力画像(図9(c))にラプラシアンフィルタのオペレータ処理を行った出力画像を表している。
【0062】
外側のマーク101は腫瘤陰影の候補をもれなく拾い上げるために表されたものであり、内側のマーク102は拾い上げた腫瘤陰影の中で、より腫瘤である危険性が高いものを示すものである。
【0063】
(第4の実施形態)
第4の実施形態は、パーソナルコンピュータ等のコンピュータ上で実行される画像処理用アプリケーションプログラムとして、本発明の画像処理方法を実現するコンピュータプログラムを提供する形態を示している。
【0064】
図12は、本発明の第4の実施形態における画像処理装置の構成を示すブロック図である。パーソナルコンピュータに本発明に係る各画像処理方法の各ステップを実行させるプログラムをインストールした画像処理装置20は、装置全体を制御するCPU(Central Processing Unit)21、RAM(Random Access Memory)等の一時記憶手段22、ハードディスクなどの記憶手段23、記録媒体から各データを読み取る光学ディスクドライブ等の補助記憶手段24、マウス及びキーボード等の入力手段25、モニタやプリンタ等の出力手段26を備えている。記録媒体に記録された本発明に係るコンピュータプログラムを記憶手段23に読み取り、制御手段21の制御により、オペレーティングシステム上で実行させることにより、パーソナルコンピュータは本発明に係る画像処理装置20として動作する。
【0065】
画像処理装置20では、記憶手段22に記憶させたコンピュータプログラムを制御手段21の制御に基づいて第1乃至第3の実施形態で説明した各ステップを実行することにより、記憶手段23などに記憶している画像や外部とのインターフェース27を通じて入力を受け付けた画像を入力画像に対し、画像処理を行い、出力手段26から出力画像を出力する。
【産業上の利用可能性】
【0066】
本発明に係る画像処理方法は、マンモグラフィなど医療診断の支援システムを始め、目立たない部分や特徴部分をより見やすく表示するといった用途への応用が期待され、従って、本発明の産業上利用可能性は極めて大きい。
【符号の説明】
【0067】
20 本発明のプログラムをインストールした画像処理装置
21 CPU(Central Processing Unit)
22 一時記憶手段(RAM)
23 記憶手段(ハードディスク)
24 補助記憶手段(光学ディスクドライブ等)
25 入力手段(マウス及びキーボード等)
26 出力手段(モニタやプリンタ等)
27 入出力インターフェース
【特許請求の範囲】
【請求項1】
コンピュータにより入力画像に含まれる各画素の画像信号の大きさを所定の目的に従って調整する画像処理方法であって、
入力画像のデータをコンピュータに入力するステップ(Sa1)と、 コントラスト強度として任意の定数mを入力するステップ(Sa2)と、
入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Dとして区画し、演算対象領域D内における注目画素とその周囲の任意の画素との2点の画素の傾きの強さを前記注目画素を中心とする全周に亘って合計し前記演算対象領域内の画素数で除算することによりその平均値を算出するステップ(Sa3)と、
ステップ(Sa3)で求めた勾配の振幅値βにステップ(Sa2)で入力したコントラスト強度m(但し、mは正の定数である。)を乗算し、任意のオフセット値γを加算するステップ(Sa4)とを具備する画像処理方法。
【請求項2】
コンピュータにより入力画像に含まれる各画素の画像信号の大きさを所定の目的に従って調整する画像処理方法であって、
入力画像のデータをコンピュータに入力するステップ(Sb1)と、
入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から求めた勾配ベクトルの角度成分θを算出するステップ(Sb2)と、
勾配ベクトルの角度成分θの絶対値に所定の定数を乗算することにより、演算対象領域C内の中心画像信号αを算出するステップ(Sb3)とを具備することを特徴とする画像処理方法。
【請求項3】
コンピュータにより入力画像に含まれる各画素の画像信号の特徴部分を抽出し、これを閉曲線で囲うための画像処理方法であって、
入力画像のデータをコンピュータに入力するステップ(Sc1)と、
入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から勾配ベクトルの大きさと演算対象領域Cの走査方向に向かう単位ベクトルと勾配ベクトルとのなす角θとを算出するステップ(Sc2)と、
前記勾配ベクトルと前記単位ベクトルとの内積を求めて規格化することにより勾配ベクトルの集中度Zを算出するステップ(Sc3)と、
ラプラシアンフィルタ、ラプラシアンオブガウシアン(LOG)フィルタその他の2階微分に基づくフィルタを乗算するステップ(Sc4)と、
前記ステップ(Sc4)に対して2値化処理を行なって出力画像を得るステップ(Sc5)とを具備する画像処理方法。
【請求項4】
アイリスフィルタによる集中度の演算ステップ((ステップSc2〜Sc3)を2回以上繰り返すことにより、画像信号の値が特に集中している特徴点を抽出することを特徴とする請求項3記載の画像処理方法。
【請求項5】
請求項1乃至4のいずれか1項に記載の各ステップを、コンピュータに実行させるためのコンピュータプログラム。
【請求項1】
コンピュータにより入力画像に含まれる各画素の画像信号の大きさを所定の目的に従って調整する画像処理方法であって、
入力画像のデータをコンピュータに入力するステップ(Sa1)と、 コントラスト強度として任意の定数mを入力するステップ(Sa2)と、
入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Dとして区画し、演算対象領域D内における注目画素とその周囲の任意の画素との2点の画素の傾きの強さを前記注目画素を中心とする全周に亘って合計し前記演算対象領域内の画素数で除算することによりその平均値を算出するステップ(Sa3)と、
ステップ(Sa3)で求めた勾配の振幅値βにステップ(Sa2)で入力したコントラスト強度m(但し、mは正の定数である。)を乗算し、任意のオフセット値γを加算するステップ(Sa4)とを具備する画像処理方法。
【請求項2】
コンピュータにより入力画像に含まれる各画素の画像信号の大きさを所定の目的に従って調整する画像処理方法であって、
入力画像のデータをコンピュータに入力するステップ(Sb1)と、
入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から求めた勾配ベクトルの角度成分θを算出するステップ(Sb2)と、
勾配ベクトルの角度成分θの絶対値に所定の定数を乗算することにより、演算対象領域C内の中心画像信号αを算出するステップ(Sb3)とを具備することを特徴とする画像処理方法。
【請求項3】
コンピュータにより入力画像に含まれる各画素の画像信号の特徴部分を抽出し、これを閉曲線で囲うための画像処理方法であって、
入力画像のデータをコンピュータに入力するステップ(Sc1)と、
入力画像中の複数の画素からなる閉領域(x,y)を演算対象領域Cとして区画し、演算対象領域C内における注目画素とその周囲の任意の画素との座標及びその2点の画素から勾配ベクトルの大きさと演算対象領域Cの走査方向に向かう単位ベクトルと勾配ベクトルとのなす角θとを算出するステップ(Sc2)と、
前記勾配ベクトルと前記単位ベクトルとの内積を求めて規格化することにより勾配ベクトルの集中度Zを算出するステップ(Sc3)と、
ラプラシアンフィルタ、ラプラシアンオブガウシアン(LOG)フィルタその他の2階微分に基づくフィルタを乗算するステップ(Sc4)と、
前記ステップ(Sc4)に対して2値化処理を行なって出力画像を得るステップ(Sc5)とを具備する画像処理方法。
【請求項4】
アイリスフィルタによる集中度の演算ステップ((ステップSc2〜Sc3)を2回以上繰り返すことにより、画像信号の値が特に集中している特徴点を抽出することを特徴とする請求項3記載の画像処理方法。
【請求項5】
請求項1乃至4のいずれか1項に記載の各ステップを、コンピュータに実行させるためのコンピュータプログラム。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【公開番号】特開2010−158317(P2010−158317A)
【公開日】平成22年7月22日(2010.7.22)
【国際特許分類】
【出願番号】特願2009−1187(P2009−1187)
【出願日】平成21年1月6日(2009.1.6)
【出願人】(000116024)ローム株式会社 (3,539)
【出願人】(507022802)Takumi Vision株式会社 (14)
【Fターム(参考)】
【公開日】平成22年7月22日(2010.7.22)
【国際特許分類】
【出願日】平成21年1月6日(2009.1.6)
【出願人】(000116024)ローム株式会社 (3,539)
【出願人】(507022802)Takumi Vision株式会社 (14)
【Fターム(参考)】
[ Back to top ]