(57)【特許請求の範囲】
【請求項１】
ヒストグラム領域を検知しディジタル放射線撮影画像のトーンスケールの再現改善を行うトーンスケール再現改善方法において、
適切なサイズを有するサンプルを用いてディジタル放射線撮影画像を疑似ランダムにサンプリングし、目的の領域を描写するステップと、
テクスチャ解析技術を用いて各サンプルを処理し、複数のテクスチャ特性を抽出するステップと、
該抽出テクスチャ特性を用いて、予め調整したニューラル回路網分類器で前記サンプルを分類しそのクラスを判定するステップと、
同一クラスに属するピクセル値を累計し、各クラスに対する個別ヒストグラムを作成し、各ヒストグラムを使用してトーンスケール再現を最適化するステップと、
を有することを特徴とするトーンスケール再現改善方法。
【発明の詳細な説明】
【０００１】
【産業上の利用分野】
本発明は、一般に、ディジタル画像処理に関し、特に、ヒストグラム領域を検知してディジタル放射線撮影画像のトーンスケール再現を改善する技術に関する。
【０００２】
【従来の技術】
医療撮像において、正確に診断できるか否かは、画像内の小さな低コントラストのディテール検出にかかっている。また、このようなディテール検出確率は、画像の可視度によってまったく異なる。撮像される物体（身体部分）、露光形状、露光技術、検出特性などの多くの因子によって、この可視度に影響を与えることができる。
【０００３】
例えば、従来のスクリーン／フィルム放射線撮影では、一般に出力トーンすなわちトーンスケール再現機能は、フィルム製造業者によって既にフィルム設計に盛り込まれており、多種多様なフィルムが試験の種類、露光技術および観察者の好みによって異なる「ルックス」を得るのに利用することができる。一方、ＵＳ−Ａ−ＲＥ．３１，８４７に記載されたように、コンピュータ放射線撮影（ＣＲ）では、画像捕捉段階と表示段階を別々にすることによって、特に必要な場合には画像または画像の一部を任意の出力コントラストで表示することができる。最終画像を表示する前にコンピュータでこの出力トーンスケールを調整するＣＲ装置の能力は、このような装置の効果的な特徴の一つである。ただし、試験タイプ、露光条件、撮像方式、イメージ捕捉装置のノイズ感度、出力装置のダイナミックレンジなどのその他の要因に影響を受けるので、適切なトーンスケールを選択または生成が簡単に行えないこともある。
【０００４】
一般に、市場で入手できるコンピュータ放射線撮影装置では、エッジ強調に続いて、トーンスケールを行う。トーンスケールの目的は、適切な変換関数を定義し、フィルムまたはソフトコピーすなわちＣＲＴディスプレイのいずれかである表示媒体の最大ダイナミックレンジに重要な画像情報をマップすることにある。さまざまな方法がＣＲ画像のトーンスケール曲線を生成するために提案されてきた。特に、ヒストグラムを基にしたトーンスケール変換は、効率が良いためこのような画像を強調するものとして最も広く用いられる方法の一つである。ＵＳ−Ａ−５，０６８，７８８、ＵＳ−Ａ−５，１９８，６９９、ＵＳ−Ａ−４，８３９，８０７、ＵＳ−Ａ−５，０４０，２２５、ＵＳ−Ａ−４，９１４，２９５、ＵＳ−Ａ−４，９５２，８０５は、ディジタル医療画像処理用にヒストグラムに基づく撮像技術の典型である。
【０００５】
一般に、このようなヒストグラムに基づく方法は、不要な領域に対応するヒストグラムのピークが所望の領域に対応するピークとは十分開きがあり、別のものとして識別することができる場合に、最も効果がある。ただし、画像に渡って不要な領域を不均一に照明してピークが重複する場合、または、コリメータ後側の分散によって、その領域の信号レベルが物体（身体部分）の高Ｘ線吸収構造に近付く場合、所望の領域と不要な領域との区別は困難になる。この結果、画像の不要な領域によって過度に影響された場合、トーンスケール変換は最適状態には及ばなくなる。
【０００６】
ディジタル投影放射線撮影技術の分野では、放射線専門医が空間イメージ分割（すなわち身体部分検知）とヒストグラム分析を合わせた画像処理技術を使用して、視覚化に最適の画像を提供することが提案されている（ＵＳ−Ａ−５，１６４，９９３参照）。これらの技術には、身体部分をフォアグラウンドとバックグラウンドから分割し、さらに、ルールベースヒストグラム分析で撮像される解剖学的領域に対応する当該の適切なヒストグラム領域を検知しなければならない。
【０００７】
【発明が解決しようとする課題】
このように、従来の重要な画像情報を表示媒体の最大ダイナミックレンジにマップしてディジタル画像処理方法を提供しイメージヒストグラムの適切な分析を容易に行うディジタル放射線撮影には問題がある。
【０００８】
【課題を解決するための手段および作用】
本発明は、上述の従来の問題を解決する方法を提供する。本発明によれば、ヒストグラム領域を検知しディジタル放射線撮影画像のトーンスケールの再現改善を行う方法にトーンスケール再現改善方法おいて、適切なサイズを有するサンプルを用いてディジタル放射線撮影画像を疑似ランダムにサンプリングし、目的の領域を描写するステップと、テクスチャ（基調）解析技術を用いて各サンプルを処理し、複数のテクスチャ特性を抽出するステップと、該抽出テクスチャ特性を用いて、予め調整したニューラル回路網分類器で前記サンプルを分類しそのクラスを判定するステップと、同一クラスに属するピクセル値を累計し、各クラスに対する個別ヒストグラムを作成し、各ヒストグラムを使用してトーンスケール再現を最適化するステップとを有することを特徴とする。
【０００９】
【実施例】
図１は、本発明を実施することができる一般的な装置１０を示す。装置１０は、入力画像１４をディジタル放射線撮影画像に変換する画像捕捉部（スキャナ）１２を有する。ディジタル画像処理部１６は、画像の外観、診断精度やその他の特性を改善することを目的とした画像処理アルゴリズムでディジタル画像信号を処理する。ディジタル画像処理部１６は、ディジタル処理演算装置、メモリ、入出力装置などの既知の構成要素を有するディジタルコンピュータでも良い。処理された出力画像２０は、画像ディスプレイ１８で確認される。本発明を用いて処理すべきディジタル放射線撮影画像信号は、記憶蛍光体（計算機放射線撮影）装置、フィルムディジタイザ、画像増倍管などから、異なるさまざまな装置で得ることができる。さらに、画像は、計算機断層撮影（ＣＴ）スキャナ、磁気共鳴画像（ＭＲＩ）スキャナ、超音波（ＵＴ）スキャナなどの診断用画像スキャナから得ることができる。ディスプレイ装置で最も広く用いられているものは、レーザプリンタ、陰極線管プリンタなどのフィルム書込装置や、ビデオモニタに装備するソフトディスプレイである。
【００１０】
本発明では、ＵＳ−Ａ−５，１６４，９３３に記載された画像処理アルゴリズムを改善する画像分類技術を提供する。本発明に係る技術は、ディジタル放射線撮影画像がこのＪａｎｇとＳｃｈａｅｔｚｉｎｇの出願に開示されたアルゴリズムに従って分割された後、使用することが好ましい。本発明に係る分類技術を実施した後、ディジタル放射線撮影画像をＵＳ−Ａ−５，１６４，９３３で開示されたトーンスケールアルゴリズムで処理する（図２のブロック２１、２２、２３）。
【００１１】
このＪａｎｇとＳｃｈａｅｔｚｉｎｇの分割アルゴリズムは、図３にさらに詳しく記載されている。この出願は、開示された分割技術をさらに徹底して説明するために参照するものである。一般に、ＣＲ画像では分割方法によって撮像される身体部分周囲の直接Ｘ線撮影領域（バックグラウンド）と、ビーム制限装置を用いて照射領域のサイズを制限した領域（フォアグラウンド）の検出および識別を行う。図３に示すように、分割は、エッジ検出ステップ２４と、ブロック生成ステップ２６と、ブロック分類ステップ２８と、ブロック改善ステップ３０とビットマップ生成ステップ３２の、５つのステップで行われる。
【００１２】
エッジ検出ステップ２４では、さまざまなエッジ検出アルゴリズムを用いる。ただし、ノイズに対して感度が低く、フォアグラウンド／バックグラウンドと物体間のコントラストに釣り合うエッジ強度を測定する能力があり、効率良く実施できることから、形態的エッジ検出器が好ましい。エッジ検出ステップ２４の出力は、高密エッジ領域の振幅が大きいピクセルと、エッジ密度がほとんどまたは全くない領域の振幅が小さいピクセルとを含むエッジ強度マップである。このため、局部のエッジ含有率は、画像のさまざまな部分のディテールレベルの標識であり、従って、これらの部分がフォアグラウンド／バックグラウンド領域寄りなのか物体寄りなのかを示す標識である。
【００１３】
エッジ検出ステップ２４後、ブロック生成ステップ２６では、画像を一連の重複しないピクセル連続ブロックに分解する。画像分析とフォアグラウンド、バックグラウンドおよび物体への分類は、ブロック単位で行われる。
【００１４】
次に、ブロック分類ステップ２８と、ブロック改善ステップ３０を行う。ブロックが基本的に同種であり相対的に低エッジ密度であるか、異種であり相対的に高エッジ密度であるかと、各ブロックが何を含むか（フォアグラウンドピクセル、バックグラウンドピクセル、物体ピクセル、フォアグラウンドピクセルおよびバックグラウンドピクセルまたはその一方と物体ピクセルの混合ピクセル）とに依存して、各ブロックが取り得る状態は１０通りある。隣接ブロックからの類似情報に加えたエッジ含有率に関する情報は、構文的推論部に供給される。この構文的推論部では、一連の臨床的または経験的に判定された決定規則を用いて画像の各ブロックをフォアグラウンド／バックグラウンドまたは物体の範疇に分類する。物体の末梢周辺のブロックについては、ブロック改善ステップ３０でそのブロック部分を後者の２つのブロックに再分類する。ビットマップ生成ステップ３２では、本発明に係る技術の好ましい実施例で用いられるテンプレートとして機能する二値（バイナリ）画像を生成する。
【００１５】
一般に、本発明に係る技術では、どのヒストグラム領域が骨、軟組織、テクスチャ属性に基づくバックグラウンドおよびフォアグラウンドに対応するかを評価することによって画像ヒストグラムを正確に分析することを容易にする。このような知識を用いると、例えば、軟組織領域のコントラストを利用して骨領域のコントラストを上げることができるようになる。
【００１６】
次に、図４について説明する。図４は、本発明に係る好ましい実施例を示す線図／ブロック図である。図示のように、図形３４は、バックグラウンド、フォアグラウンド、当該の領域を含む元の放射線撮影画像である。図形３６は、上述の分割処理から生じるバイナリ画像である。本発明に係る画像処理技術によれば、当該の描写領域（両手）は、図形３８に示したように、任意にサンプリングされる。各サンプルは予め処理されて、テクスチャ特性が抽出される（ステップ４０）。次に、タイプすなわち骨、組織、バックグラウンドまたはフォアグラウンドに関しニューラル分類技術で抽出された特性を用いて分類される（ステップ４２）。同一クラスに属するピクセル値は累計されて、骨、組織、バックグラウンド、フォアグラウンドに対して個別のヒストグラム推定値となる（ステップ４４）。次に、これらのヒストグラム推定値を使用して、骨領域または組織領域に対するトーンスケールアルゴリズムを最適化し、視覚化または診断のためにこれらの領域の表示を強調することができる。
【００１７】
本発明に係る各ステップを次に詳細に説明する。図４の任意サンプリングステップ３８は、各ピクセル毎にディジタル値１２ビットを有することができる分割された身体部分から、３９×３９のピクセル下位領域などのピクセル下位領域を抽出することによって実行される。本発明に係るサンプリング技術は、Ｐｅｔｅｒ  Ｇ．  Ａｎｄｅｒｓｏｎ，  Ｌｉｎｅａｒ  Ｐｉｘｅｌ  Ｓｈｕｆｆｌｉｎｇ  ｆｏｒ  Ｉｍａｇｅ  Ｐｒｏｃｅｓｓｉｎｇ，  ａｎ  Ｉｎｔｒｏｄｕｃｔｉｏｎ，  Ｊｏｕｒｎａｌ  ｏｆ  Ｅｌｅｃｔｒｏｎｉｃ  Ｉｍａｇｉｎｇ，  （２），  Ａｐｒｉｌ，  １９９３，  Ｐａｇｅｓ  １４７−１５４  の論文に述べられた線形ピクセルシャフリング技術であることが好ましい。図５は、好例となるソフトウェアルーチンである。
【００１８】
図４の処理およびテクスチャ抽出ステップすなわち事前処理及び特性算出ステップ４０を、図６から図１２までを参照に説明する。図６に示すように、第１の事前処理ステップ５０では、１６などの因数で各サンプルに対するデータを圧縮する。例えば、あるサンプルの各１２ビットピクセルは、８ビットピクセルに圧縮される。次に、局部適応ヒストグラム等化処理を各圧縮サンプルに行う。好例となる技術は、Ｐｒａｔｔ，  Ｗ．，  （１９３３），  Ｄｉｇｉｔａｌ  Ｉｍａｇｅ  Ｐｒｏｃｅｓｓｉｎｇ，  Ｗｉｌｅｙ，  Ｎｅｗ  Ｙｏｒｋ，  Ｐａｇｅ  ２７５−２８４に記載されている。例えば、各９×９ウィンドウ位置で、ヒストグラムを算出する（ピクセルコード値の度数表）（図７、９（ａ），（ｂ））。次に、等化ヒストグラムが、ヒストグラムに基づく累積度数表として算出される。さらに、中央の３×３ピクセル領域を等化ヒストグラムに基づいて再度マッピングする（図９（ｃ），（ｄ）参照）。３ピクセルだけ移動し、連続するウィンドウでこの処理を繰り返す（図７、９（ａ），（ｂ））。その後、３２×３２ピクセルなどの各サンプルの中間領域を抽出する（図６、ステップ５４）。図１０および図１１は、局部適応ヒストグラム等化前後の骨サンプルヒストグラムをそれぞれ示す。
【００１９】
本発明に係る方法のつぎのステップは、特性抽出を行う。物体または画像のテクスチャは、任意の、周期的な、平滑な、粗雑な、荒い、などの形容詞で説明することができ、統計的技術または構造的技術を用いて表現することができる。応力線に続く補強支柱を成す海綿質骨が柱構造であるため、整形外科的な骨にはユニークなテクスチャ属性がある。放射線撮影画像処理で、これらの属性を記録すると、組織、バックグラウンドまたはフォアグラウンドに比較して、ユニークなテクスチャとなる。本発明に係る好ましい実施例では、灰色レベル同時生起（ＧＬＣ）行列（Ｈａｒａｌｉｃｋ，  Ｒ．，  Ｓｈａｎｍｕｇａｍ，  Ｋ．，  Ｄｉｎｓｔｅｉｎ，  Ｉ．，  （１９７３），  Ｔｅｘｔｕｒａｌ  Ｆｅａｔｕｒｅｓ  ｆｏｒ，  Ｉｍａｇｅ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ，  ＩＥＥＥ  Ｔｒａｎｓ．  Ｓｙｓｔ．  Ｍａｎ．  Ｃｙｂｅｒ．，  （ＳＭＣ−３），  Ｎｏｖ．  Ｐａｇｅｓ  １４１−１５２と、Ｈａｒａｌｉｃｋ，  Ｒ．，  （１９７９），  Ｓｔａｔｉｓｔｉｃａｌ  ａｎｄ  Ｓｔｒｕｃｔｕｒａｌ  Ａｐｐｒｏａｃｅｓ  ｔｏ  Ｔｅｘｔｕｒｅ，  Ｐｒｏｃ．  ｏｆ  ＩＥＥＥ，  （６７），  Ｍａｙ，  Ｐａｇｅｓ  ７８６−８０４を参照）と、隣接灰色レベル従属（ＮＧＬＤ）行列（Ｓｕｎ，  Ｃ．  Ｗｅｅ，  Ｗ．，  （１９８２），  Ｎｅｉｇｈｂｏｒｉｎｇ  Ｇｒａｙ  Ｌｅｖｅｌ  Ｄｅｐｅｎｄｅｎｃｅ  Ｍａｔｒｉｘ  ｆｏｒ  Ｔｅｘｔｕｒｅ  Ｃｌａｓｓｉｆｉｃａｔｉｏｎ，  Ｃｏｍｐｕｔｅｒ  ｖｉｓｏｎ  Ｇｒａｐｈｉｃｓ  ａｎｄ  Ｉｍａｇｅ  Ｐｒｏｃｅｓｓｉｎｇ，  （２３），  Ｐａｇｅｓ  ３４１−３５２  を参照）２つの統計学的なテクスチャ解析方法を使用する。
【００２０】
一般に、ＧＬＣ行列は、イメージの灰色トーンが互いにすなわち相互生起しなければならないという平均的空間関係に基づいている。基本となる前提は、テクスチャ情報がさまなざまなピクセルの差異および配向に基づいて算出される一連の灰色トーン空間依存行列によって適切に指定されることである。このような灰色レベル同時生起行列は、二次結合条件各釣り密度関数の評定値ｐ（ｉ，ｊ，ｄ，ｏ）であり、一般に０度、４５度、９０度、  １３５度の４方向で算出される。各密度関数は、サンプル間の間隔ｄであり、方向が角度ｏである場合、灰色レベルｉから灰色レベルｊに至る確率である。画像がＮｇ灰色レベルを有する場合、灰色レベル同時生起行列はＮｇ×Ｎｇになる。
【００２１】
与えられた例では、すべてのテクスチャ特性は３２×３２ピクセルサンプルから算出されたこれらの行列から抽出される。例えば、ｄは１から３であり、抽出された特性は、エネルギー、コントラスト、相関関係、局部類似性およびエントロピーである。
【００２２】
ＧＬＣ行列方法の制限は、各方向に関し一度に行列および対応する特性を４回算出しなければならないことにある。ＮＧＬＤ行列方法では、一度に一方向算出する代わりに、ピクセルと近傍すべての関係を一度に考慮することでこの制限を解消している。これによって、角度に依存せず、線形灰色レベル変換に基づく不変性を容易にし、計算による負荷を低減する。
【００２３】
ＮＧＬＤ行列は、画像の全要素と特定の距離ｄを開けた近傍のグレートーン関係から形成される。次に、すべてのテクスチャ特性は、このＮＧＬＤ行列から求められる。行列のサイズは、Ｎｇ×Ｎｎであり、Ｎｇは可能な灰色レベル数であり、Ｎｎは、可能な近傍数である。これは、画像の灰色度変化度数と見做すことができ、画像のヒストグラムと意味的に似ている。与えられた例では、スモールナンバー（ｓｍａｌｌ  ｎｕｍｂｅｒ）強調、ラージナンバー（ｌａｒｇｅ  ｎｕｍｂｅｒ）強調、数値不均一性、二次モーメントおよびエントロピーを用いた。
【００２４】
次に図１２について説明する。図１２は、特性抽出ステップで用いるソフトウェアルーチンである。
【００２５】
サンプルを分類する前に、ニューラルネットワーク分類器を放射線撮影画像のライブラリから手作業で選択した既知のデータ（上述の事前処理した骨サンプル、軟組織サンプル、バックグラウンドサンプル、フォアグラウンドサンプルと、算出されたテクスチャ特性）に関して調整する（ＵＳ−Ａ−５，０４８，０９７参照）。
【００２６】
学習過程は以下の通りである。
【００２７】
目標：分類器を調整し、いままでに調査したことがない未知のサンプルを骨、組織、バックグラウンド、フォアグラウンドから成る４クラスの１つとして認識する。
【００２８】
手順：操作条件下で認められる可能性があるサンプルの母集団全体を代表するサンプリングとして役立つと思われるさまざまな試験および画像条件に関する何千ものサンプルを収集する。各サンプルを選択する際、適切なクラスすなわち、骨、組織、バックグラウンド、フォアグラウンドを示すラベルを添付する。
【００２９】
各サンプルを操作条件下で使用すべき方法に一致する方法で事前処理する。
【００３０】
上述のテクスチャ特性を算出し、標準化し、入力ベクトルを生成する。
【００３１】
出力ベクトルをセットアップする。
【００３２】
例
骨は、｛１，０，０，０｝としてマップされる。
【００３３】
組織は、｛０，１，０，０｝としてマップされる。
【００３４】
バックグラウンドは、｛０，０，１，０｝としてマップされる。
【００３５】
フォアグラウンドは、｛０，０，０，１｝としてマップされる。
【００３６】
次に、これらの入力ベクトルと出力ベクトルの組合わせを、指揮された学習アルゴリズムに提供する。好ましい実施例は、三層グローバル接続ニューラルネットワークである。後退誤差伝播法は、最急下降法を用いる学習アルゴリズムとして好ましい。最急下降法は、標的と出力ベクトル間の誤差を反復形式で最小化するのに利用される。図１３は、ニューラルネットワーク分類器を示す図である。
【００３７】
図４に示すように、次のステップでは、予め調整したニューラルネットワーク分類器でサンプルを分類する（ステップ４２）。サンプルが骨タイプからなると認められる場合、中央ピクセルおよび一般には３×３の周辺の近傍は、すべて骨として分類される。対応するピクセルデータは、他の骨ピクセル値とともに累計され、骨ヒストグラムを構成する（図４、ステップ４４）。同様に、その他のクラスについても行う。図１４（ａ）および図１４（ｂ）は、それぞれ、ヒストグラムの構成例と、ヒストグラム構成を実行するソフトウェアルーチンを示す。図１５は、上述の４クラスに対して構成されたヒストグラムを示す三次元グラフ図である。
【００３８】
特定のクラスに対するヒストグラムをさらに正確に識別する場合の本発明に係る分類技術の成果を図１６から図１９に示す。図１６および図１７は、それぞれ実際の骨サンプルヒストグラムと組織サンプルヒストグラムである。図１８および図１９は、それぞれ、本発明に係る技術で分類された骨サンプルヒストグラムと組織サンプルヒストグラムである。実際のヒストグラムと算出したヒストグラムの相関関係は明白である。
【００３９】
図２０は、身体部分全体のヒストグラムと、それに基づくトーンスケール曲線の他、骨のみのヒストグラムと、それに基づくトーンスケール曲線である。骨領域を知っていれば、左のポイントと、特に、右のポイントについてはＵＳ−Ａ−５，１６４，９９３のトーンスケールアルゴリズムを用いて、推測しなければならない。判断を誤ると、コントラストは下がり、画像表示に一貫性が無くなる。骨領域と組織領域の境界が分かれば、左のポイントと右のポイントを正確に選択することができるので、画像表示に矛盾がなくなる。
【００４０】
【発明の効果】
本発明に係る技術を用いることによって、画像内のピクセルすべてを分析しなければならない複雑な規則に基づくシステムを開発せずに、さまざまな分布、特に骨と軟組織の正確な近似値を自動的に得ることができる利点がある。上述したようなシステムを開発するには時間がかかり、また故障しやすいという問題がある。本発明のシステムに係る別の利益は、モニターに対して容易に個別化できることである。例えば、放射線医師が骨の腫瘍を探している場合、骨領域内のトーンスケールのダイナミックレンジをすべて用いて選択し、軟組織またはバックグラウンドに対しては一切用いないと思われる。あるいは、骨折の場合は、ダイナミックレンジの２０％までを軟組織の腫脹検査に、８０％を骨折の位置を決定するために選択することができる。また、分割ステップで、バックグラウンド領域またはフォアグラウンド領域が分割マスクに含まれるようなエラーが生じた場合は、ヒストグラム解析方法は失敗するが、テクスチャ特性を用いてヒストグラムを分解した場合には、この問題を解決することができるという利点もある。
【図面の簡単な説明】
【図１】本発明で用いることができる撮像装置のブロック図である。
【図２】ディジタル放射線撮影画像の画像処理に用いる本発明に係る一実施例のブロック図である。
【図３】好ましい画像分割手続きを示すブロック図である。
【図４】本発明に係る好ましい実施例の部分的に線図を含めた部分的ブロック図である。
【図５】画像サンプル処理ルーチンのフローチャート図である。
【図６】画像サンプル事前処理ステップを示すブロック図である。
【図７】局部適応ヒストグラム等化を説明する線図である。
【図８】本発明に係る事前処理ステップを実行する際に利用するフローチャート図である。
【図９】（ａ），（ｂ），（ｃ），（ｄ）は、本発明に係る技術の説明に用いられるグラフ図である。
【図１０】本発明に係る技術の説明に用いられるグラフ図である。
【図１１】本発明に係る技術の説明に用いられるグラフ図である。
【図１２】テクスチャ特性抽出ステップを実行した場合に使用するフローチャート図である。
【図１３】ニューラルネットワーク分類器を示す図である。
【図１４】（ａ）はヒストグラム構成の一例であり、（ｂ）はヒストグラム構成ステップのフローチャート図である。
【図１５】本発明を説明する際に使用するグラフ図である。
【図１６】本発明を説明する際に使用するグラフ図である。
【図１７】本発明を説明する際に使用するグラフ図である。
【図１８】本発明を説明する際に使用するグラフ図である。
【図１９】本発明を説明する際に使用するグラフ図である。
【図２０】本発明を説明する際に使用するグラフ図である。
【符号の説明】
１０  装置
１２  画像捕捉部（スキャナ）
１４  入力画像
１６  ディジタル画像処理部
１８  画像ディスプレイ
２０  出力画像
２１  分割ステップ
２２  分類ステップ
２３  トーンスケールステップ
２４  エッジ検出ステップ
２６  ブロック生成ステップ
２８  ブロック分類ステップ
３０  ブロック改善ステップ
３２  ビットマップ生成ステップ
３４、３６  図形
３８  任意サンプリングステップ
４０  処理およびテクスチャ抽出ステップ
４２  サンプル分類ステップ
４４  分布構成ステップ
４６  トーンスケールステップ
５０  事前処理ステップ
５２  局部適応ヒストグラム等化処理ステップ
５４  中央領域抽出ステップ
【図１】
0003619270.tif 000002
【図２】
0003619270.tif 000003
【図３】
0003619270.tif 000004
【図４】
0003619270.tif 000005
【図５】
0003619270.tif 000006
【図６】
0003619270.tif 000007
【図７】
0003619270.tif 000008
【図８】
0003619270.tif 000009
【図９】
0003619270.tif 000010
【図１０】
0003619270.tif 000011
【図１１】
0003619270.tif 000012
【図１２】
0003619270.tif 000013
【図１３】
0003619270.tif 000014
【図１４】
0003619270.tif 000015
【図１５】
0003619270.tif 000016
【図１６】
0003619270.tif 000017
【図１７】
0003619270.tif 000018
【図１８】
0003619270.tif 000019
【図１９】
0003619270.tif 000020
【図２０】
0003619270.tif 000021