JP2002092588A - 画像処理装置及び画像処理方法並びに記憶媒体 - Google Patents

Abstract

(57)【要約】 【課題】 矩形領域毎の画像データの情報に基づいてウ
ェーブレット変換の変換係数を得、この変換係数に基づ
いて画像データを得ることで、ノイズ除去を行うこと。 【解決手段】 前処理された原画像全体を複数の矩形
領域に分割する（Ｓ３０１）。分割されたそれぞれの矩
形領域４０１のウェーブレット係数を出力する（Ｓ３０
２）。また、それぞれの矩形領域４０１に対して、テク
スチャ解析を行い、その解析結果を出力する（Ｓ３０
３）。次に、テクスチャ解析の結果に基づいて、それぞ
れの領域４０１のウェーブレット係数に対して係数変換
を施す（Ｓ３０４）。変換係数の中でも、高周波数成分
に対して係数変換する。次に、変換されたそれぞれの領
域のウェーブレット係数に対して、逆離散ウェーブレッ
ト変換を施すことで、ノイズの除去された画像を出力す
る（Ｓ３０６）。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、画像データを補正
する画像処理装置及び画像処理方法並びに記憶媒体に関
するものである。

【０００２】

【従来の技術】近年のデジタル技術の進歩により放射線
画像をデジタル画像信号に変換し、該デジタル画像信号
に対して画像処理を施し、ＣＲＴ等に表示、あるいはプ
リント出力することが行われている。ところで、放射線
画像の撮影においては、患者に与える影響などから撮影
時のＸ線量は少ないことが望ましい。しかし、Ｘ線量を
少なくして撮影された画像には量子ノイズが多く含まれ
ることが知られており、このような量子ノイズは診断を
行う際の妨害となる可能性がある。そのため、従来から
このようなノイズを除去する処理が検討されており、例
えば、単純なメディアンフィルタによるノイズ除去処理
や、平滑化画像を用いて高周波成分を抽出し、ノイズを
除去する方法（いわゆるフィルタ処理）などが行われて
きた。また近年、入力画像を複数の周波数帯域に分割
し、各周波数帯域ごとに独立した処理を施すことにより
ノイズの除去を行う多重周波数処理が検討されている。

【０００３】

【発明が解決しようとする課題】平滑化画像を用いて高
周波成分を抽出しノイズの除去を行うフィルタ処理で
は、単一の周波数帯域を扱うため、ノイズ成分が広い周
波数帯域に分布している場合には効果的にノイズ除去を
行えないという問題がある。これを避けるために、大き
さの異なる（つまり周波数の異なる）複数のフィルタを
同時に利用することも行われているが、処理に必要な計
算コストが大幅に増加することや、ノイズ除去のために
フィルタの周波数特性を最適に保つためには、被写体に
応じたフィルタサイズの調整が必須となり、汎用性に乏
しいという問題がある。

【０００４】また、ノイズ除去に多重周波数処理を用い
ることにより上記の問題は大きく軽減されるが、同一の
画像内を一律に処理するため、例えば、抜けの良い領域
と悪い領域、細かいテクスチャの多い領域とテクスチャ
のほとんど存在しない領域などをどちらも最適に処理す
ることは難しいという問題がある。

【０００５】本発明は以上の問題点に対して鑑みてなさ
れたものであり、矩形領域毎の画像データの情報に基づ
いてウェーブレット変換の変換係数を得、この変換係数
に基づいて画像データを得ることで、ノイズ除去を行う
ことを目的とする。

【０００６】

【課題を解決するための手段】本発明の目的を達成する
ために、例えば本発明の画像処理装置は以下の構成を備
える。すなわち、画像データを矩形領域に分割する分割
手段と、前記分割手段による矩形領域内の画像データに
基づいて画素情報を得る画像解析手段と、前記分割手段
による矩形領域内の画像データに対してウェーブレット
変換を行い、変換係数を得るウェーブレット変換手段
と、前記画像解析手段による前記画素情報に応じて、矩
形領域単位の変換係数に対して補正を行う補正手段と、
前記補正手段により補正された変換係数を含む全変換係
数に対して、逆ウェーブレット変換を行う逆ウェーブレ
ット変換手段とを備える。

【０００７】本発明の目的を達成するために、例えば本
発明の画像処理装置は以下の構成を備える。すなわち、
画像データを矩形領域に分割する分割手段と、前記分割
手段による矩形領域内の画像データに基づいて画素情報
を得る画像解析手段と、複数の異なった特性を有するウ
ェーブレット変換フィルタを格納する格納手段と、前記
画像解析手段による前記画素情報に応じて、矩形領域毎
のウェーブレット変換フィルタを前記格納手段から選択
する選択手段と、前記分割手段による矩形領域内の画像
データに対して、前記選択手段によるウェーブレット変
換フィルタを用いてウェーブレット変換を行い、変換係
数を得るウェーブレット変換手段と、前記ウェーブレッ
ト変換手段による変換係数に対して、逆ウェーブレット
変換を行う逆ウェーブレット変換手段とを備える。

【０００８】

【発明の実施の形態】以下添付図面に従って、本発明を
好適な実施形態に従って詳細に説明する。

【０００９】［第１の実施形態］図１は、本実施形態に
おける画像処理装置としてのＸ線撮影装置１００の概略
構成を示す図である。

【００１０】すなわち、Ｘ線撮影装置１００は、撮影さ
れた画像のノイズ除去を行う機能を有するＸ線の撮影装
置であり、前処理回路１０６、ＣＰＵ１０８、メインメ
モリ１０９、操作パネル１１０、画像表示器１１１、画
像処理回路１１２を有し、画像処理回路１１２は、テク
スチャ解析回路１１３、領域分割回路１１５、離散ウェ
ーブレット（ＤＷＴ）変換回路１１６、係数変換回路１
１７、特徴抽出回路１１８、階調変換回路１１９を備え
ており、ＣＰＵバス１０７を介して互いにデータの受け
渡しを行うことができる。

【００１１】また、Ｘ線撮影装置１００は、前処理回路
１０６に接続されたデータ収集回路１０５と、データ収
集回路１０５に接続された２次元Ｘ線センサ１０４及び
Ｘ線発生回路１０１とを備えており、これらの各回路も
ＣＰＵバス１０７に接続されている。

【００１２】以上の構成を備えるＸ線撮像装置１００
が、被写体１０３を撮像し、後述する各種の画像処理を
行う一連の処理の大まかなフローチャートを図２に示
す。なお図２に示したフローチャートに従ったプログラ
ムコードはメインメモリ１０９、もしくは図不示のＲＯ
Ｍなどに格納され、ＣＰＵ１０８により読み出され、実
行されるものとする。また、メインメモリ１０９には、
ＣＰＵ１０８での処理に必要な各種のデータなどが記憶
されるものであると共に、ＣＰＵ１０８の作業用として
のワークメモリを含む。また、ＣＰＵ１０８は、メイン
メモリ１０９を用いて、操作パネル１１０からの操作に
従った装置１００全体の動作制御等を行う。これにより
Ｘ線撮影装置１００は、以下のように動作する。

【００１３】先ず、Ｘ線発生回路１０１は、被写体１０
３に対してＸ線ビーム１０２を放射する。Ｘ線発生回路
１０１から放射されたＸ線ビーム１０２は、被写体１０
３を減衰しながら透過して、２次元Ｘ線センサ１０４に
到達し、２次元Ｘ線センサ１０４により放射線画像（Ｘ
線画像）として出力される（ステップＳ２０１）。ここ
では、２次元Ｘ線センサ１０４から出力されるＸ線画像
を、例えば人体画像とする。

【００１４】データ収集回路１０５は、２次元Ｘ線セン
サ１０４から出力されたＸ線画像を電気信号に変換して
前処理回路１０６に供給する。前処理回路１０６は、デ
ータ収集回路１０５からの信号（Ｘ線画像信号）に対し
て、オフセット補正処理やゲイン補正処理等の前処理を
行う（ステップＳ２０２）。この前処理回路１０６で前
処理が行われたＸ線画像信号は原画像として、ＣＰＵ１
０８の制御により、ＣＰＵバス１０７を介して、メイン
メモリ１０９、画像処理回路１１２に転送される。

【００１５】ここで、画像処理回路１１２を構成する各
部について説明する。１１３は放射線画像内の矩形領域
内のテクスチャ解析を行うテクスチャ解析回路、１１５
は放射線画像を複数の矩形領域に分割する領域分割回
路、１１６は矩形領域に対して離散ウェーブレット変換
を施し、ウェーブレット変換係数を得る離散ウェーブレ
ット変換回路、１１７はウェーブレット係数を後述する
規則に従い、変換する係数変換回路である。また画像処
理回路１１２は、階調変換に必要な特徴量を抽出する特
徴量抽出回路１１８と、得られた特徴量に応じて階調変
換を行う階調変換回路１１９を備える。

【００１６】次にステップＳ２０３で、画像処理回路１
１２で行われるノイズ除去処理について、具体的な処理
のフローチャートを図３に示し、以下同図を用いて説明
する。なお図３に示したフローチャートに従ったプログ
ラムコードはメインメモリ１０９、もしくは図不示のＲ
ＯＭに格納され、ＣＰＵ１０８により読み出され、実行
されるものとする。また、図４に矩形領域分割処理の様
子を示し、図４も用いて説明する。

【００１７】前処理回路１０６で前処理された原画像は
ＣＰＵバス１０７を介して画像処理装置１１２に転送さ
れる。画像処理装置１１２では領域分割回路１１５が、
図４（ｂ）に示すとおり、原画像全体を複数の矩形領域
に分割する（ステップＳ３０１）。分割されたそれぞれ
の矩形領域４０１に対して、離散ウェーブレット変換回
路１１６が離散ウェーブレット変換を施し、ウェーブレ
ット係数を出力する（ステップＳ３０２）。またそれと
同時に、分割されたそれぞれの矩形領域４０１に対し
て、テクスチャ解析回路１１３がテクスチャ解析を行
い、その解析結果を出力する（ステップＳ３０３）。こ
こではテクスチャ解析の具体的な方法として、注目矩形
領域内の画素値の分散値ｄを算出する。

【００１８】次に係数変換回路１１７は、ステップＳ３
０３におけるテクスチャ解析の結果に基づいて、ステッ
プＳ３０２で得られたそれぞれの領域４０１のウェーブ
レット係数に対して係数変換を施す（ステップＳ３０
４）。本実施形態では変換係数の中でも、高周波数成分
（例えば最高レベルのＨＨに含まれる変換係数）に対し
て係数変換する。係数変換の具体的なアルゴリズムは特
に限定しないが、本実施形態では、図６に示すフローチ
ャートに従った処理により係数変換を行う。

【００１９】まずステップＳ６０１において、上述の分
散値ｄが閾値（分散値ｄに応じて決まる）θより大きい
か否かを判断し、ｄ＞θの場合、処理をステップＳ６０
２に移行し、高周波数成分に含まれる変換係数に対して
ｎビット（ｎ＞１）下方向にビットシフト（変換係数を
ｈとすると、ｈ＞＞ｎ）を行う。

【００２０】これは、ｄ＞θの場合、注目矩形領域内の
画素値の分布が一定量（θ）以上一定であることを示
す。つまり、ノイズにより画像的には一様な矩形領域と
なる。よって、この注目矩形領域内の画素値の高周波数
成分を下げることで、ノイズの影響を緩和することで、
注目矩形内の画像をより鮮明にすることができる。

【００２１】次にステップＳ６０３では、高周波数成分
に含まれる全ての変換係数に対してビットシフトを行っ
たか否かの判断を行う。

【００２２】ステップＳ３０４における係数変換の処理
は以上の処理の他にも、離散ウェーブレット変換の分解
レベルに応じた係数変換ＬＵＴ（ルックアップテーブ
ル）を作成し、このテーブルに応じて係数変換を行って
も良い。

【００２３】最後に、ステップＳ３０４で係数変換回路
１１７によって変換されたそれぞれの領域のウェーブレ
ット係数を用いて、離散ウェーブレット変換回路１１６
が逆離散ウェーブレット変換を施し、それぞれの領域４
０１の画像データを再構成し（ステップＳ３０５）、そ
れぞれの領域４０１を結合して、ノイズの除去された画
像を出力する（ステップＳ３０６）。

【００２４】離散ウェーブレット変換回路１１６による
離散ウェーブレット変換（ステップＳ３０２）及び逆離
散ウェーブレット変換（ステップＳ３０５）は一般によ
く知られる変換処理であり以下のように動作する。

【００２５】離散ウェーブレット変換回路１１６は入力
した画像信号に対して２次元の離散ウェーブレット変換
処理を行い、変換係数を計算して出力するものである。
メインメモリ１０９に記憶された画像データは、離散ウ
ェーブレット変換回路１１６により順次読み出されて変
換処理が行われ、変換処理による変換係数を再びメイン
メモリ１０９に書きこまれる。

【００２６】なお、この変換係数をメインメモリ１０９
に格納する際、元の画像データが格納された領域には上
書きせず、他の領域に格納する。つまり、この画像デー
タには一切の変更を加えない。それは、ステップＳ３０
３におけるテクスチャ解析処理において、この画像デー
タを参照するからである。このようにすれば、仮にステ
ップＳ３０３における処理がステップＳ３０２における
処理と同時に行われず、例えばステップＳ３０２におけ
る処理が終わってから、ステップＳ３０３における処理
が実行される場合でも、画像データは損なわれることな
く存在するので、ステップＳ３０２における処理とステ
ップＳ３０３における処理を同時に行った場合の結果と
同じ結果が得られるのは自明である。

【００２７】次に、本実施形態における離散ウェーブレ
ット変換回路１１６の処理の構成を図５（ａ）に示す。
同図において、入力された画像信号は遅延素子およびダ
ウンサンプラの組み合わせにより、偶数アドレスおよび
奇数アドレスの信号に分離され、２つのフィルタｐおよ
びｕによりフィルタ処理が施される。同図ｓおよびｄ
は、各々1次元の画像信号に対して１レベルの分解を行
った際のローパス係数およびハイパス係数を表してお
り、次式により計算されるものとする。

【００２８】 ｄ（ｎ）＝ｘ（２ｎ＋１） −ｆｌｏｏｒ（（ｘ（２ｎ）＋ｘ（２ｎ＋２））／２） （式１） ｓ（ｎ）＝ｘ（２ｎ） ＋ｆｌｏｏｒ（（ｄ（ｎ−１）＋ｄ（ｎ））／４） （式２） ただし、ｘ（ｎ）は変換対象となる画像信号である。

【００２９】以上の処理により、画像信号に対する１次
元の離散ウェーブレット変換処理が行われる。２次元の
離散ウェーブレット変換は、１次元の変換を画像の水平
・垂直方向に対して順次行うものであり、その詳細は公
知であるのでここでは説明を省略する。

【００３０】図５（ｂ）は２次元の変換処理により得ら
れる２レベルの変換係数群の構成例であり、画像信号は
異なる周波数帯域の係数列ＨＨ１，ＨＬ１，ＬＨ
１，．．．，ＬＬに分解される。同様に逆離散ウェーブ
レット変換については次のように行われる。メインメモ
リ１０９に記憶された変換係数は離散ウェーブレット変
換回路１１６により順次読み出されて変換処理が行わ
れ、変換係数に基づいた画像データは再びメインメモリ
１０９に書きこまれる。

【００３１】本実施形態における離散ウェーブレット変
換回路１１６による逆離散ウェーブレット変換処理の構
成は図５（ｃ）に示すものとする。入力された変換係数
はｕおよびｐの2つのフィルタ処理を施され、アップサ
ンプリングされた後に重ね合わされて画像信号x'が出力
される。これらの処理は次式により行われる。

【００３２】 ｘ’（２ｎ）＝ｓ’（ｎ） −ｆｌｏｏｒ（（ｄ’（ｎ−１）＋ｄ’（ｎ））／４） （式３） ｘ’（２ｎ＋１）＝ｄ’（ｎ） ＋ｆｌｏｏｒ（（ｘ’（２ｎ）＋ｘ’（２ｎ＋２））／２） （式４） 以上の処理により、変換係数に対する１次元の逆離散ウ
ェーブレット変換処理が行われる。２次元の逆離散ウェ
ーブレット変換は、１次元の逆変換を画像の水平・垂直
方向に対して順次行うものであり、その詳細は公知であ
るのでここでは説明を省略する。

【００３３】なおステップＳ２０４においては前述の通
り、特徴量抽出回路１１８により階調変換に必要な特徴
量を抽出し、ステップＳ２０５においては、ステップＳ
２０４において得られた特徴量に応じて階調変換を行
う。なおステップＳ２０４，Ｓ２０５における各処理は
本実施形態では特に限定しないため、ここでは詳細な説
明は省略する。

【００３４】以上説明したとおり、本実施形態における
画像処理装置及び画像処理方法によれば、離散ウェーブ
レット変換による多重周波数処理を行うため、単一の周
波数帯域を扱う、いわゆるフィルタ処理と比べて、ノイ
ズ成分が広い周波数帯域に分布している場合にもノイズ
除去を行える効果がある。また、大きさの異なる（つま
り周波数の異なる）複数のフィルタを同時に利用する処
理に比べて、処理に必要な計算コストを低減する効果が
有り、フィルタサイズの調整など煩雑な作業を行わずに
汎用性の高い処理が行える効果がある。また画像を複数
の矩形領域に分割し、テクスチャ解析結果に基づいた個
別のパラメータでノイズ除去処理を行うため、同一の画
像内を一律に処理する場合と比べて、各領域に応じた処
理を施すことが可能となる効果があり、結果として従来
手法に比べて比較的高性能でノイズ除去された画像を得
ることが可能となる効果がある。

【００３５】［第２の実施形態］第１の実施形態におい
て、テクスチャ解析部１１３によるテクスチャ解析結果
は、係数変換回路１１７における係数変換処理にのみ使
われていた。本実施形態ではこのテクスチャ解析結果
を、離散ウェーブレット変換において用いるフィルタの
選択に用いる。つまり、矩形領域内の画像に対してテク
スチャ解析を行った際、このテクスチャ解析結果を用い
て、この矩形領域に対する離散ウェーブレット変換フィ
ルタを選択する。

【００３６】使用するフィルタの違いは周波数帯の分解
特性の違いとしてその特徴が表れ、例えば９−７型フィ
ルタ（詳細は公知なので省略する）は高周波成分を細か
く分離する特徴があり、例えばＨａｒｒ型フィルタ（詳
細は公知なので省略する）は周波数帯の分離幅が広く、
同一次元の分解処理で９−７型よりもより低周波までの
分解処理が可能である特徴がある。このような周波数帯
の分解特性に特徴のある離散ウェーブレット変換フィル
タを複数用意しておき、テクスチャ解析の結果に応じて
使い分ける。なお、これらを含む複数のフィルタは、離
散ウェーブレット変換回路１１６内の不図示のメモリ内
に格納しておいても良いし、メインメモリ１０９内に格
納しておいても良い。

【００３７】また本実施形態においてのテクスチャ解析
の具体的な方法として、一例を挙げると、例えば、注目
矩形領域内の画素値の分散値ｄを求める。そして、上述
のフィルタの選択は、閾値θを用いて次にように実現可
能である。

【００３８】 ｉｆ （ｄ＞θ） ｔｈｅｎ 「９−７型フィルタ使用」 ｅｌｓｅ ｔｈｅｎ 「Ｈａａｒ型フィルタ使用」 上記は極めて簡単な例であり、テクスチャ特徴量として
画素値の分散値のみを利用しているが、その他のテクス
チャ特徴量を用いても良いし、複数のテクスチャ特徴量
を組み合わせて利用してももちろん良い。

【００３９】なお、本実施形態における上述の画像処理
方法のフローチャートは、図２と同じであり、ステップ
Ｓ２０３における具体的な処理として、図３において、
ステップＳ３０２，Ｓ３０３の処理を直列に行い、しか
も、処理の順番を、ステップＳ３０３，Ｓ３０２の順に
行い、ステップＳ３０４における処理は行わないとした
処理のフローチャートとなる。

【００４０】［他の実施形態］なお、本発明は、複数の
機器（例えばホストコンピュータ、インタフェイス機
器、リーダ、プリンタなど）から構成されるシステムに
適用しても、一つの機器からなる装置（例えば、複写
機、ファクシミリ装置など）に適用してもよい。

【００４１】また、本発明の目的は、前述した実施形態
の機能を実現するソフトウェアのプログラムコードを記
録した記憶媒体（または記録媒体）を、システムあるい
は装置に供給し、そのシステムあるいは装置のコンピュ
ータ（またはCPUやMPU）が記憶媒体に格納されたプログ
ラムコードを読み出し実行することによっても、達成さ
れることは言うまでもない。この場合、記憶媒体から読
み出されたプログラムコード自体が前述した実施形態の
機能を実現することになり、そのプログラムコードを記
憶した記憶媒体は本発明を構成することになる。また、
コンピュータが読み出したプログラムコードを実行する
ことにより、前述した実施形態の機能が実現されるだけ
でなく、そのプログラムコードの指示に基づき、コンピ
ュータ上で稼働しているオペレーティングシステム(OS)
などが実際の処理の一部または全部を行い、その処理に
よって前述した実施形態の機能が実現される場合も含ま
れることは言うまでもない。

【００４２】さらに、記憶媒体から読み出されたプログ
ラムコードが、コンピュータに挿入された機能拡張カー
ドやコンピュータに接続された機能拡張ユニットに備わ
るメモリに書込まれた後、そのプログラムコードの指示
に基づき、その機能拡張カードや機能拡張ユニットに備
わるCPUなどが実際の処理の一部または全部を行い、そ
の処理によって前述した実施形態の機能が実現される場
合も含まれることは言うまでもない。

【００４３】本発明を上記記憶媒体に適用する場合、そ
の記憶媒体には、先に説明した（図２，３、６に示す）
フローチャート、及び／又は第２の実施形態で説明した
第２の実施形態のフローチャートに対応するプログラム
コードが格納されることになる。

【００４４】

【発明の効果】以上説明したとおり、本発明によれば、
矩形領域毎の画像データの情報に基づいてウェーブレッ
ト変換の変換係数を得、この変換係数に基づいて画像デ
ータを得ることで、ノイズ除去を行うことができた。

【図面の簡単な説明】

【図１】本発明の第１の実施形態における画像処理装置
としてのＸ線撮像装置１００の概略構成を示す図であ
る。

【図２】Ｘ線撮像装置１００が、被写体１０３を撮像
し、各種の画像処理を行う一連の処理の大まかなフロー
チャートである。

【図３】ステップＳ２０３で、画像処理回路１１２で行
われるノイズ除去処理について、具体的な処理のフロー
チャートである。

【図４】矩形領域分割処理の様子を示した図である。

【図５】（ａ）は離散ウェーブレット変換回路１１６の
処理の構成を示す図であり、（ｂ）は２次元の離散ウェ
ーブレット変換処理により得られる２レベルの変換係数
群の構成例を示す図であり、（ｃ）は離散ウェーブレッ
ト変換回路１１６による逆離散ウェーブレット変換処理
の構成を示す図である。

【図６】係数変換を行う処理のフローチャートである。

───────────────────────────────────────────────────── フロントページの続き (72)発明者 佐藤 眞 東京都大田区下丸子３丁目30番２号 キヤ ノン株式会社内 Ｆターム(参考） 4C093 AA01 AA26 CA06 CA32 CA36 EB17 EE01 FC18 FC19 FD03 FD04 FD05 FD09 FD13 FF03 FF08 FF09 FF19 FF21 FF50 FH02 5B057 AA08 BA03 CA02 CA08 CA12 CA16 CB02 CB08 CB12 CB16 CC02 CE02 CE20 CH07 CH18 5C076 AA01 AA36 AA40 BA04 BA06 5C077 LL02 LL19 PP01 PP21 PP49 PQ08 PQ12 PQ22 PQ23

Claims (13)

【特許請求の範囲】
1. 【請求項１】 画像データを矩形領域に分割する分割手
   段と、 前記分割手段による矩形領域内の画像データに基づいて
   画素情報を得る画像解析手段と、 前記分割手段による矩形領域内の画像データに対してウ
   ェーブレット変換を行い、変換係数を得るウェーブレッ
   ト変換手段と、 前記画像解析手段による前記画素情報に応じて、矩形領
   域単位の変換係数に対して補正を行う補正手段と、 前記補正手段により補正された変換係数を含む全変換係
   数に対して、逆ウェーブレット変換を行う逆ウェーブレ
   ット変換手段とを備えることを特徴とする画像処理装
   置。
2. 【請求項２】 画像データを矩形領域に分割する分割手
   段と、 前記分割手段による矩形領域内の画像データに基づいて
   画素情報を得る画像解析手段と、 複数の異なった特性を有するウェーブレット変換フィル
   タを格納する格納手段と、 前記画像解析手段による前記画素情報に応じて、矩形領
   域毎のウェーブレット変換フィルタを前記格納手段から
   選択する選択手段と、 前記分割手段による矩形領域内の画像データに対して、
   前記選択手段によるウェーブレット変換フィルタを用い
   てウェーブレット変換を行い、変換係数を得るウェーブ
   レット変換手段と、 前記ウェーブレット変換手段による変換係数に対して、
   逆ウェーブレット変換を行う逆ウェーブレット変換手段
   とを備えることを特徴とする画像処理装置。
3. 【請求項３】 更に、被写体を撮像することで、当該被
   写体の前記画像データを得る撮像手段を備えることを特
   徴とする請求項１又は２に記載の画像処理装置。
4. 【請求項４】 前記撮像手段は、Ｘ線ビームにより前記
   被写体を撮像し、得られる画像データはＸ線画像データ
   であることを特徴とする請求項３に記載の画像処理装
   置。
5. 【請求項５】 前記画像解析手段は、矩形領域内の画像
   データに対してテクスチャ解析を行うことを特徴とする
   請求項１乃至４のいずれか１項に記載の画像処理装置。
6. 【請求項６】 前記画像解析手段は前記テクスチャ解析
   として、矩形領域内の画素値の分散値と当該分散値に応
   じた閾値を算出することを特徴とする請求項５に記載の
   画像処理装置。
7. 【請求項７】 前記画像解析手段は更に、離散ウェーブ
   レット変換の分解レベルに応じたテーブルを作成するテ
   ーブル作成手段を含み、前記変更手段は当該テーブルを
   用いて、前記矩形領域毎に前記変換係数に対して変更を
   加えることを特徴とする請求項１乃至４のいずれか１項
   に記載の画像処理装置。
8. 【請求項８】 前記テーブルは変換係数ルックアップテ
   ーブルであることを特徴とする請求項７に記載の画像処
   理装置。
9. 【請求項９】 前記補正手段は、前記変換係数の中で
   も、比較的高周波数成分に含まれる変換係数に対して変
   更を加えることを特徴とする請求項１に記載の画像処理
   装置。
10. 【請求項１０】 画像データを矩形領域に分割する分割
    工程と、 前記分割工程による矩形領域内の画像データに基づいて
    画素情報を得る画像解析工程と、 前記分割工程による矩形領域内の画像データに対してウ
    ェーブレット変換を行い、変換係数を得るウェーブレッ
    ト変換工程と、 前記画像解析工程による前記画素情報に応じて、矩形領
    域単位の変換係数に対して補正を行う補正工程と、 前記補正工程により補正された変換係数を含む全変換係
    数に対して、逆ウェーブレット変換を行う逆ウェーブレ
    ット変換工程とを備えることを特徴とする画像処理方
    法。
11. 【請求項１１】 画像データを矩形領域に分割する分割
    工程と、 前記分割工程による矩形領域内の画像データに基づいて
    画素情報を得る画像解析工程と、 複数の異なった特性を有するウェーブレット変換フィル
    タを所定の格納手段に格納する格納工程と、 前記画像解析工程による前記画素情報に応じて、矩形領
    域毎のウェーブレット変換フィルタを前記所定の格納手
    段から選択する選択工程と、 前記分割工程による矩形領域内の画像データに対して、
    前記選択工程によるウェーブレット変換フィルタを用い
    てウェーブレット変換を行い、変換係数を得るウェーブ
    レット変換工程と、 前記ウェーブレット変換工程による変換係数に対して、
    逆ウェーブレット変換を行う逆ウェーブレット変換工程
    とを備えることを特徴とする画像処理方法。
12. 【請求項１２】 コンピュータが読み取り可能なプログ
    ラムコードを格納する記憶媒体であって、 画像データを矩形領域に分割する分割工程のプログラム
    コードと、 前記分割工程による矩形領域内の画像データに基づいて
    画素情報を得る画像解析工程のプログラムコードと、 前記分割工程による矩形領域内の画像データに対してウ
    ェーブレット変換を行い、変換係数を得るウェーブレッ
    ト変換工程のプログラムコードと、 前記画像解析工程による前記画素情報に応じて、矩形領
    域単位の変換係数に対して補正を行う補正工程のプログ
    ラムコードと、 前記補正工程により補正された変換係数を含む全変換係
    数に対して、逆ウェーブレット変換を行う逆ウェーブレ
    ット変換工程のプログラムコードとを備えることを特徴
    とする記憶媒体。
13. 【請求項１３】 コンピュータが読み取り可能なプログ
    ラムコードを格納する記憶媒体であって、 画像データを矩形領域に分割する分割工程のプログラム
    コードと、 前記分割工程による矩形領域内の画像データに基づいて
    画素情報を得る画像解析工程のプログラムコードと、 複数の異なった特性を有するウェーブレット変換フィル
    タを所定の格納手段に格納する格納工程のプログラムコ
    ードと、 前記画像解析工程による前記画素情報に応じて、矩形領
    域毎のウェーブレット変換フィルタを前記所定の格納手
    段から選択する選択工程のプログラムコードと、 前記分割工程による矩形領域内の画像データに対して、
    前記選択工程によるウェーブレット変換フィルタを用い
    てウェーブレット変換を行い、変換係数を得るウェーブ
    レット変換工程のプログラムコードと、 前記ウェーブレット変換工程による変換係数に対して、
    逆ウェーブレット変換を行う逆ウェーブレット変換工程
    のプログラムコードとを備えることを特徴とする記憶媒
    体。