(57)【特許請求の範囲】
【請求項１】
注目領域を含む画像に対して画像処理を行う画像処理装置であって、
画像を予め定められたサイズのタイルに分割する分割手段と、
注目タイルに含まれる前記注目領域の割合を計算し、該割合が閾値よりも大きい場合には、前記注目タイル内の高周波成分を細かく分離する機能を有する第１の周波数変換部を用いて前記注目タイルに対する周波数変換を行い、前記割合が前記閾値以下である場合には、周波数帯の分離幅が広く、同一次元の分解処理で前記第１の周波数変換部よりも低周波までの分解処理が可能な機能を有する第２の周波数変換部を用いて前記注目タイルに対する周波数変換を行う周波数変換手段と、
前記周波数変換手段で周波数変換された前記注目タイルに含まれる変換係数に対して成分変換を行う成分変換手段と、
前記成分変換手段により成分変換された前記注目タイルに対して前記周波数変換手段とは逆の変換を行う逆周波数変換手段と
を備えることを特徴とする画像処理装置。
【請求項２】
更に画像を撮像する撮像手段を備え、
前記注目領域を含む画像は前記撮像手段により撮像されたものであることを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記撮像手段は、Ｘ線により被写体の画像を撮像することを特徴とする請求項２に記載の画像処理装置。
【請求項４】
注目領域を含む画像に対して画像処理を行う画像処理装置が行う画像処理方法であって、
画像を予め定められたサイズのタイルに分割する分割工程と、
注目タイルに含まれる前記注目領域の割合を計算し、該割合が閾値よりも大きい場合には、前記注目タイル内の高周波成分を細かく分離する機能を有する第１の周波数変換部を用いて前記注目タイルに対する周波数変換を行い、前記割合が前記閾値以下である場合には、周波数帯の分離幅が広く、同一次元の分解処理で前記第１の周波数変換部よりも低周波までの分解処理が可能な機能を有する第２の周波数変換部を用いて前記注目タイルに対する周波数変換を行う周波数変換工程と、
前記周波数変換工程で周波数変換された前記注目タイルに含まれる変換係数に対して成分変換を行う成分変換工程と、
前記成分変換工程で成分変換された前記注目タイルに対して前記周波数変換工程とは逆の変換を行う逆周波数変換工程と
を備えることを特徴とする画像処理方法。
【請求項５】
コンピュータに請求項１乃至３の何れか１項に記載の画像処理装置が有する各手段として機能させるためのコンピュータプログラムを格納した、コンピュータ読み取り可能な記憶媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、注目領域を含む画像に対して画像処理を行う画像処理装置及びその方法並びに記憶媒体に関するものである。
【０００２】
【従来の技術】
近年のデジタル技術の進歩により放射線画像をデジタル画像信号に変換し、該デジタル画像信号に対して周波数処理などの画像処理を施し、ＣＲＴ等に表示、あるいはプリント出力することが行われている。ところで、このような周波数処理は複数の周波数帯の画像成分に分離し、周波数帯毎の画像成分を増加又は減弱することで行われている。
【０００３】
【発明が解決しようとする課題】
しかしながら、上述の方法では、画像全体に同一の周波数分離処理を行っているため、領域ごとの周波数特性が異なる場合でも一律な周波数分解処理をおこなっていた。そのため、領域毎に適する周波数処理を効率よく行えない問題があった。
【０００４】
また、画像の領域に毎に含有される周波数成分はことなるが、全領域に同一のフィルタをかけると圧縮効率がよくない問題があった。
【０００５】
本発明は上記のような課題を解決するためになされたもので、領域毎に効果的に周波数帯毎の画像成分を変換することを目的とする。
【０００６】
【課題を解決するための手段】
本発明の目的を達成するために、例えば本発明の画像処理装置は以下の構成を備える。
【０００７】
即ち、注目領域を含む画像に対して画像処理を行う画像処理装置であって、
画像を予め定められたサイズのタイルに分割する分割手段と、
注目タイルに含まれる前記注目領域の割合を計算し、該割合が閾値よりも大きい場合には、前記注目タイル内の高周波成分を細かく分離する機能を有する第１の周波数変換部を用いて前記注目タイルに対する周波数変換を行い、前記割合が前記閾値以下である場合には、周波数帯の分離幅が広く、同一次元の分解処理で前記第１の周波数変換部よりも低周波までの分解処理が可能な機能を有する第２の周波数変換部を用いて前記注目タイルに対する周波数変換を行う周波数変換手段と、
前記周波数変換手段で周波数変換された前記注目タイルに含まれる変換係数に対して成分変換を行う成分変換手段と、
前記成分変換手段により成分変換された前記注目タイルに対して前記周波数変換手段とは逆の変換を行う逆周波数変換手段と
を備えることを特徴とする。
【００１０】
本発明の目的を達成するために、例えば、本発明の画像処理方法は以下の構成を備える。
【００１１】
即ち、注目領域を含む画像に対して画像処理を行う画像処理装置が行う画像処理方法であって、
画像を予め定められたサイズのタイルに分割する分割工程と、
注目タイルに含まれる前記注目領域の割合を計算し、該割合が閾値よりも大きい場合には、前記注目タイル内の高周波成分を細かく分離する機能を有する第１の周波数変換部を用いて前記注目タイルに対する周波数変換を行い、前記割合が前記閾値以下である場合には、周波数帯の分離幅が広く、同一次元の分解処理で前記第１の周波数変換部よりも低周波までの分解処理が可能な機能を有する第２の周波数変換部を用いて前記注目タイルに対する周波数変換を行う周波数変換工程と、
前記周波数変換工程で周波数変換された前記注目タイルに含まれる変換係数に対して成分変換を行う成分変換工程と、
前記成分変換工程で成分変換された前記注目タイルに対して前記周波数変換工程とは逆の変換を行う逆周波数変換工程と
を備えることを特徴とする。
【００１４】
【発明の実施の形態】
以下添付図面を参照して、本発明を好適な実施形態に従って詳細に説明する。
［第１の実施形態］
図１は、本実施形態の画像処理装置としてのＸ線撮影装置１００の構成を示す図である。すなわちＸ線撮影装置１００は、撮影された画像に対して周波数帯毎に処理を行う機能を有するＸ線の撮影装置であり、前処理回路１０６、ＣＰＵ１０８、メインメモリ１０９、操作パネル１１０、画像表示器１１１、画像処理回路１１２を備えており、ＣＰＵバス１０７を介して互いにデータ授受されるようになされている。
【００１５】
またＸ線撮影装置１００は、前処理回路１０６に接続されたデータ収集回路１０５と、データ収集回路１０５に接続された２次元Ｘ線センサ１０４及びＸ線発生回路１０１とを備えており、これらの各回路はＣＰＵバス１０７にも接続されている。
【００１６】
上述の構成を備えるＸ線撮影装置１００の各部について説明する。まずメインメモリ１０９は、ＣＰＵ１０８での処理に必要な各種のデータなどが記憶されるものであると共に、ＣＰＵ１０８の作業用としてのワークメモリを含む。
【００１７】
ＣＰＵ１０８は、メインメモリ１０９を用いて、操作パネル１１０からの操作に従った装置全体の動作制御等を行う。これによりＸ線撮影装置１００は、以下のように動作する。
【００１８】
先ず、Ｘ線発生回路１０１は、被写体（被検査体）１０３に対してＸ線ビーム１０２を放射する。Ｘ線発生回路１０１から放射されたＸ線ビーム１０２は、被検査体１０３を減衰しながら透過して、２次元Ｘ線センサ１０４に到達し、２次元Ｘ線センサ１０４によりＸ線画像として出力される。ここでは、２次元Ｘ線センサ１０４から出力されるＸ線画像を、例えば人体画像等とする。
【００１９】
データ収集回路１０５は、２次元Ｘ線センサ１０４から出力されたＸ線画像を電気信号に変換して前処理回路１０６に供給する。前処理回路１０６は、データ収集回路１０５からの信号（Ｘ線画像信号）に対して、オフセット補正処理やゲイン補正処理等の前処理を行う。この前処理回路１０６で前処理が行われたＸ線画像信号は原画像として、ＣＰＵ１０８の制御により、ＣＰＵバス１０７を介して、メインメモリ１０９、画像処理回路１１２に転送される。
【００２０】
同回路１１２において、１１３は画像を小領域に分割し、小領域（以後タイル）毎に後述の離散ウェーブレット変換回路１１４によるウェーブレット変換方式を決定する領域分割回路、１１４は複数のウェーブレット変換方式を有する離散ウェーブレット変換回路であり、原画像に対して離散ウェーブレット変換を施し、各周波数帯の画像成分（ウェーブレット変換係数）を得る。
【００２１】
１１８は離散ウェーブレット変換回路１１４で得られた各周波数帯の画像成分を変換する画像成分変換回路であり、１１５は画像成分変換回路１１８で変換した画像成分を逆変換して複数の逆ＤＷＴ変換方式で画像を合成する逆ＤＷＴ変換回路である。
【００２２】
以上の構成を備える本実施形態におけるＸ線撮影装置１００が行う処理を、図２，３，４に示す同処理のフローチャートを用いて以下説明する。
【００２３】
前処理回路１０６で前処理された原画像はＣＰＵバス１０７を介して画像処理回路１１２に転送される（ステップＳ２００）。ここで原画像を例えば図５（ａ）に示した胸部正面画像とする。画像処理回路１１２では、はじめに領域分割回路１１３が原画像をタイルに分割し、分割した領域毎にウェーブレット変換方式（フィルタ）を決定する（ステップＳ２０１）。この処理の詳細を図３に示すフローチャートに従って説明する。
【００２４】
まず、領域分割回路１１３は図５（ｂ）に例示するように、予め決まるサイズの小四角に原画像を分割する（ステップＳ３０１）。次に目的領域として、例えば図５（ｂ）に例示するように、肺部などの解剖学的に特徴のある領域を目的領域５０１として抽出する（ステップＳ３０２）。このような目的領域５０１を抽出する方法としては公知の方法を用いるが、この方法は特に限定されたものではない。
【００２５】
次に、領域分割回路１１３は目的領域とタイルの重なりの割合をタイル毎に計算する（ステップＳ３０３）。そして、タイル中に含まれる目的領域の面積の割合を重なり度として計算する（ステップＳ３０３）。具体的には、タイル中の画素数に対する、目的領域を構成する画素数の割合を計算する。この重なり度の値が予め設定された一定閾値Ｔｈより大きければ（ステップＳ３０４）、そのタイルに対して第一のＤＷＴ変換回路を用いて周波数分解を行なう（この処理を第一のＤＷＴ変換と呼ぶ）こととし（ステップＳ３０６）、Ｔｈ以下ならば（ステップＳ３０４）、第二のＤＷＴ変換回路を用いて周波数分解（この処理を第二のＤＷＴ変換と呼ぶ）を行なうこととする（ステップＳ３０５）。
【００２６】
ここで、各ＤＷＴ変換回路の特性の違いは用いるフィルタ形の相違であり、それぞれ周波数帯の分解特性に特徴がある。例えば第一のＤＷＴ変換回路では高周波成分を細かく分離する特徴がある（例えば９−７型フィルタ：詳細は公知なので省略する。）を使用し、第二のＤＷＴ変換回路では周波数帯の分離幅が広く、同一次元の分解処理で第一のＤＷＴ変換回路よりもより低周波までの分解処理が可能なフィルタ（例えばＨａｒｒ型フィルタ：詳細は公知なので省略する。）を用いる。このような周波数帯の分解特性に特徴のあるＤＷＴ変換回路を複数用意しておき、領域や部位に応じて使い分ける。
【００２７】
そして、ＤＷＴ変換回路１１４では領域分割回路１１３で決められたタイルとウェーブレット変換方式に従い、タイルごとに周波数分解処理を行なう（ステップＳ２０２）。
【００２８】
次にステップＳ２０２の処理の詳細を、図４に示すフローチャートに従って説明する。ＤＷＴ変換回路１１４はタイル毎に領域分割回路１１３で決められたＤＷＴ変換回路を選択し、タイル内の画像信号に対して２次元の離散ウェーブレット変換処理を行い、画像成分（変換係数）を計算して出力するものである。なお、タイル毎の変換係数が求まった際には、どのＤＷＴ変換回路を用いたかを示すフラグ（ＤＷＴ変換回路選択フラグ）を、各タイルに添付しておく。
【００２９】
本実施形態におけるＤＷＴ変換回路１１４の構成を図６（ａ）に示す。同図の構成を有する回路において、入力された画像信号は遅延素子およびダウンサンプラの組み合わせにより、偶数アドレスおよび奇数アドレスの信号に分離され、２つのフィルタｐおよびｕによりフィルタ処理が施される。ｓおよびｄは、各々１次元の画像信号に対して１レベルの分解を行った際のローパス係数およびハイパス係数を表しており、次式により計算されるものとする。
【００３０】
ｄ（ｎ）＝ｘ（２ｎ＋１）−
ｆｌｏｏｒ（（ｘ（２ｎ）＋ｘ（２ｎ＋２））／２）          （式１）
ｓ（ｎ）＝ｘ（２ｎ）＋
ｆｌｏｏｒ（（ｄ（ｎ−１）＋ｄ（ｎ））／４）              （式２）
ただし、ｘ（ｎ）は変換対象となる画像信号である。以上の処理により、画像信号に対する１次元の離散ウェーブレット変換処理が行われる。尚、上述の様に各ＤＷＴ変換回路では用いるフィルタが異なるものである。
【００３１】
２次元の離散ウェーブレット変換は、１次元の変換を画像の水平・垂直方向に対して順次行うものであり、その詳細は公知であるのでここでは説明を省略する。図６（ｂ）は２次元の変換処理により得られる２レベルの変換係数群の構成例であり、画像信号は異なる周波数帯域の画像成分ＨＨ１，ＨＬ１，ＬＨ１，．．．，ＬＬに分解される（ステップＳ４０１）。図６（ｂ）においてＨＨ１，ＨＬ１，ＬＨ１，．．．，ＬＬ等（以下サブバンドと呼ぶ）が周波数帯毎の画像成分を示す。
【００３２】
次に、画像成分変換回路１１８は例えば図７に示す変換曲線を周波数帯毎に作成し（ステップＳ４０２）、画像成分を変換する（ステップＳ４０３）。つまり値がＴｈ１以上の変換係数（同図では入力系数）に対しては、曲線７０２に従ったスケーリングを行っている。
【００３３】
図７は画像成分変換回路１１８で作成された変換曲線の一例であり、横軸が入力係数（入力画像成分）、縦軸が出力係数（出力画像成分）を示す。
【００３４】
ここで図７に示す変換曲線において、変換曲線７０２の傾き、閾値７０１は、小領域毎、サブバンド毎にステップＳ４０２において作成されるものを用いる。例えば同一のＤＷＴ変換回路を用いた場合であってもタイル毎に異なる周波数毎の変換曲線形を用いるものである。この変換曲線の曲線形としては、例えば解剖学的領域と原画像の画素値に応じて変更するものである。これにより、タイル毎にきめ細かく周波数処理を行なうことができるものである。
【００３５】
逆ＤＷＴ変換回路１１５は、タイルに添付された前述のＤＷＴ変換回路選択フラグを参照して、２次元の離散ウェーブレット変換処理の際に用いたフィルタを特定し、特定したフィルタに応じた逆ＤＷＴ変換に用いるフィルタを特定し、以下のようにタイル毎に逆ＤＷＴ変換を行う（ステップＳ４０４）。逆ＤＷＴ変換回路１１５の構成は図６（ｃ）に示すものとする。入力された画像成分はｕおよびｐの２つのフィルタ処理を施され、アップサンプリングされた後に重ね合わされて画像信号ｘ’が出力される。これらの処理は次式により行われる。
【００３６】
ｘ’（２ｎ）＝ｓ’（ｎ）−
ｆｌｏｏｒ（（ｄ’（ｎ−１）＋ｄ’（ｎ））／４）      （式３）
ｘ’（２ｎ＋１）＝ｄ’（ｎ）＋
ｆｌｏｏｒ（（ｘ’（２ｎ）＋ｘ’（２ｎ＋２））／２）  （式４）
以上の処理により、変換係数に対する１次元の逆離散ウェーブレット変換処理が行われる。２次元の逆離散ウェーブレット変換は、１次元の逆変換を画像の水平・垂直方向に対して順次行うものであり、その詳細は公知であるのでここでは説明を省略する。また、各ＤＷＴ変換回路毎にフィルタ形が異なるものであり、例えば第一の逆ＤＷＴ変換回路では９−７型フィルタに対応するフィルタを用い、例えば第二の逆ＤＷＴ変換回路ではＨａｒｒ型フィルタに対応するフィルタを用いるものである。
【００３７】
また、画像を保存する場合には、逆変換する前の画像データを保存するものであるが、この場合には画像成分毎の変換処理をする前のものを保存してもよいし、成分変換処理後の画像データを保存してもよい。成分変換前の画像を保存する場合には、逆変換処理を行なう前に成分変換処理を行なう。これにより、成分変換前の画像を保存する場合には原画像に復元が可能であるが、復元処理において成分変換処理を行なう必要がある。
【００３８】
一方、成分変換後に画像データを保存する場合には、原画像への復元はできないが、高周波成分を遮断することで圧縮効率が上がる効果がある。また、復元処理中に画像変換処理をおこなわなくてもよい。さらに、高周波成分を多く含む領域（例えば肺部）に高周波の分離性のよいフィルタを用いることで、画像データとフィルタの相関度が上がり、画像を表現するためのデータ数は減じ、結果として圧縮率が上がる効果がある。同様に、低周波成分を多く含む領域（例えば腹部）には低周波高周波の分離性のよいフィルタを用いることで、画像データとフィルタの相関度が上がり、画像を表現するためのデータ数は減じ、結果として圧縮率が上がる効果がある。したがって、全体として一律なフィルタをかけるよりも領域毎にことなるフィルタを用いることで全体のデータの圧縮率が上がる効果がある。特に人体の場合は胸部や腹部などの解剖学的領域が異なるとその周波数成分もことなり、解剖学的領域毎に異なるフィルタを用いることでデータ圧縮率などの効率がよくなる効果がる。
【００３９】
ところで、肺部では、血管や細かい骨梁は高周波成分で構成されており、これら部分を強調するためには高周波成分を変換する必要が生じる。しかし、本実施形態では、肺部には高周波成分を細かく分離できるＤＷＴ変換回路を用いることで、的確に目的とする高周波成分を強調することができる。
【００４０】
一方腹部などの臓器部は中低周波数で構成されており、臓器部を強調するためには低中周波の広範囲にわたり強調する必要がある。しかし、本実施形態では腹部には分離周波数帯の幅が広く、低次数の分離回数で低周波まで分離できるフィルタ形を用いるため、的確に低中周波数の画像成分を強調することができる。このように領域ごとの特性に応じてＤＷＴ回路を使い分けることで、的確に目的とする領域の強調処理を効率よく行なうことができる効果がある。
【００４１】
以上の様に本実施形態では、領域毎に異なる分解特性を有するＤＷＴ変換回路を用いることで、効率よく目的とする周波数分解を行なうことができる効果がある。また、目的とする周波数帯幅で領域毎に周波数数分解を行なうことで、目的とする周波数成分をきめ細かく強調することができる効果がある。さらに、解剖学的領域毎に用いるＤＷＴ変換回路を使い分けることで、解剖学的領域毎にきめ細かく周波数処理を行なうことが可能であり、画像全体の周波数処理を的確に行なえる効果あがる。また、小領域毎に周波数強調する方法、強さを変更するので、画像細部にわたってきめ細かく周波数処理を行なうことが可能であり、さらに画像全体で最適な周波数処理を行なうことが可能である。
【００４２】
［第２の実施形態］
本実施形態の画像処理装置としてのＸ線撮影装置は、第１の実施形態におけるＸ線撮影装置に、ＤＷＴ変換回路、逆ＤＷＴ変換回路をもう一つ加えた構成を備え、第１の実施形態と同様に、各タイルに占める目的領域の割合に応じてＤＷＴの変換計画を選択し、選択した変換計画に従って第１の実施形態と同様に各タイルに２次元の離散ウェーブレット変換を施す機能を備える。また、離散ウェーブレット変換を施した後は第１の実施形態と同様に成分変換を施し、各タイルに施した離散ウェーブレット変換の逆変換を施す。
【００４３】
以下、本実施形態におけるＸ線撮影装置について説明する。図８に本実施形態におけるＸ線撮影装置の構成を示す。なお図１示した装置と同じ部分は同じ番号で示している。図１に示した装置と異なる点は、ＤＷＴ変換回路、逆ＤＷＴ変換回路をもう一つ加えた点と、記憶回路１１７付け加えた点である。
【００４４】
以下本装置における処理を、図９，１０に示す同処理のフローチャートを参照して以下説明する。
【００４５】
前処理回路１０６で前処理された原画像はＣＰＵバス１０７を介して画像処理回路１１２に転送される（ステップＳ９００）。ここで原画像を例えば図５（ａ）に示した胸部正面画像とする。画像処理回路１１２では、はじめに領域分割回路１１３が原画像をタイルに分割し、分割した領域（タイル）毎にＤＷＴ変換計画を決定する（ステップＳ９０１）。この処理の詳細を図１０に示すフローチャートに従って説明する。
【００４６】
まず、領域分割回路１１３は図５（ｂ）に例示するように、予め決まるサイズのタイルに原画像を分割する（ステップＳ１００１）。次に目的領域として、例えば図５（ｂ）に例示するように、肺部などの解剖学的に特徴のある領域を目的領域５０１として抽出する（ステップＳ１００２）。このような目的領域５０１を抽出する方法としては公知の方法を用いるが、この方法は特に限定されたものではない。
【００４７】
次に、領域分割回路１１３は目的領域とタイルの重なりの割合をタイル毎に計算する（ステップＳ１００３）。そして、タイルに含まれる目的領域の面積の割合を重なり度として計算する（ステップＳ１００３）。具体的には、タイル中の画素数に対する、目的領域を構成する画素数の割合を計算する。この重なり度の値が予め設定された一定閾値Ｔｈより大きければ（ステップＳ１００４）、そのタイルを第一のＤＷＴ変換計画を選択して周波数分解を行うこととし（ステップＳ１００６）、Ｔｈ以下ならば（ステップＳ１００４）、第二のＤＷＴ変換計画を選択して周波数分解を行うこととする（ステップＳ１００５）。
【００４８】
ここで第一のＤＷＴ変換計画としては例えばＤＷＴ変換の第１，２の分解レベルでは第一のＤＷＴ変換回路（９−７型）を用い、第３，４の分解レベルでは第二のＤＷＴ変換回路（５−３型）を用い、第５，６の分解レベルでは第三のＤＷＴ変換回路（Ｈａｒｒ型）を用いるものである。ここで第三のＤＷＴ変換回路は、第１の実施形態におけるＸ線撮影装置に新たに設けたＤＷＴ変換回路であって、周波数帯の分離幅が広く、同一次元の分解処理で第二のＤＷＴ変換回路よりもより低周波までの分解処理が可能なフィルタである。本実施形態ではこの第三のＤＷＴ変換回路としてＨａｒｒ型を用いているがこれに限定されるものではない。
【００４９】
一方、第二のＤＷＴ変換計画としては例えばＤＷＴ変換の第１，２の分解レベルでは第二のＤＷＴ変換回路（５−３型）を用い、第３，４の分解レベルでは第三のＤＷＴ変換回路（Ｈａｒｒ型）を用いるものである。以上の第一、第二のＤＷＴ変換計画の内容はメインメモリ１０９に格納されているものとする。
【００５０】
ここで分解レベルとはＤＷＴ変換処理を行う回数であって、例えば２回目のＤＷＴ変換処理を第二の分解レベルと呼ぶ。
【００５１】
これにより、胸部には第１のＤＷＴ変換計画が選択され、腹部には第二のＤＷＴ変換計画が選択されることになる。胸部には細かい血管、細かい器官、中位の血管、中位の器官、肋骨などが存在する。胸部において、細かい血管、細かい器官などは高周波数を主の成分とし、中位の血管、中位の器官、肋骨などは高中周波数を主の成分とする。また、肋骨は中低周波数を主の成分とする。従って、胸部画像の画像全体の鮮鋭度を調節するために、周波数帯毎に適するフィルタを用いることが効率的である。例えば、高周波を強調するには９−７型のような高周波成分の分離特性に優れるフィルタを用いるのがよく、高中周波数帯をまたぐ画像成分を強調するには、９−７型より周波数帯域が広く、かつ高中周波数の分離特性に優れる５−３型を用いるのがよい。そして肋骨みたいに中低周波数を主とする画像を強調するにはより周波数帯域が広く、低周波成分の分離特性に優れるＨａｒｒ型を用いるのが効率がよい。
【００５２】
例えば、９−７型のフィルタ一つで周波数分解を低周波数まで分離するには、高次の分解レベルを要し、効率が悪くなるためである。更に、例えば周波数帯域の狭い９−７型で低周波数帯域の画像成分を強調するためには、複数の分解レベルに渡る画像成分を調整する必要もあり、効率が悪くなるためでもある。
【００５３】
一方、腹部は主に中低周波数成分で構成されているため、５−３型、Ｈａｒｒ型のフィルタ順で分解する方が効率がよい。高周波成分を含まない領域に９−７型のフィルタをかけても上述の理由でよくないからである。さらに、第二のＤＷＴ変換計画では周波数帯域の広いフィルタを用いるため、分解レベルを低回数に押さえられる効果もあるためである。さらに、各分解レベル毎に適するフィルタを用いることで画像の圧縮効率もよくなる効果があるためである。
【００５４】
そして、領域分割回路１１３で決定されたタイル毎のＤＷＴ変換計画とこの計画に用いたフィルタを記憶回路１１７に保存する（ステップＳ９０２）。そしてＤＷＴ変換回路１１４では記憶回路１１７に保存されたＤＷＴ変換計画とフィルタに従って、タイル毎に周波数分解処理を行う（ステップＳ９０３）。この際、ＤＷＴ変換計画の内容もメインメモリ１０９から読み込む。周波数分解処理は第１の実施形態と同じ２次元の離散ウェーブレット変換処理を用いる。なお本実施形態では周波数分解処理の際に、例えば第一のＤＷＴ変換計画を例に取ると、例えば第１の分解レベルでは９−７型のフィルタを周波数分解処理に用い、ＨＨ１，ＨＬ１，ＬＨ１，ＬＬ１の４つのサブバンドに分解する。そして第二の分解レベルではＬＬ１に対して９−７型のフィルタを用いて更に周波数分解処理を行い、ＨＨ２，ＨＬ２，ＬＨ２，ＬＬ２の４つのサブバンドに分解する。そして第三の分解レベルではＬＬ２に対して５−３型のフィルタを用いて更に周波数分解処理を行い、ＨＨ３，ＨＬ３，ＬＨ３，ＬＬ３の４つのサブバンドに分解する。以下同様の処理を繰り返す。
【００５５】
そしてタイル毎に画像成分を変換する（ステップＳ９０４）が、この処理は第１の実施形態におけるステップＳ４０３の処理と同じである。
【００５６】
逆ＤＷＴ変換回路１１５は、ステップＳ９０２で記憶回路１１７に保存したタイル毎のＤＷＴ変換計画に用いたフィルタを読み込んで、２次元の離散ウェーブレット変換処理の際に用いたフィルタを特定し、特定したフィルタに応じた逆ＤＷＴ変換に用いるフィルタを特定し、タイル毎に逆ＤＷＴ変換を行う（ステップＳ９０５）。逆ＤＷＴ変換の方法は第１の実施形態で説明した方法と同じである。
【００５７】
［他の実施形態］
なお、本発明は、複数の機器（例えばホストコンピュータ、インタフェイス機器、リーダ、プリンタなど）から構成されるシステムに適用しても、一つの機器からなる装置（例えば、複写機、ファクシミリ装置など）に適用してもよい。
【００５８】
また、本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記憶媒体（または記録媒体）を、システムあるいは装置に供給し、そのシステムあるいは装置のコンピュータ（またはCPUやMPU）が記憶媒体に格納されたプログラムコードを読み出し実行することによっても、達成されることは言うまでもない。この場合、記憶媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記憶した記憶媒体は本発明を構成することになる。また、コンピュータが読み出したプログラムコードを実行することにより、前述した実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているオペレーティングシステム(OS)などが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【００５９】
さらに、記憶媒体から読み出されたプログラムコードが、コンピュータに挿入された機能拡張カードやコンピュータに接続された機能拡張ユニットに備わるメモリに書込まれた後、そのプログラムコードの指示に基づき、その機能拡張カードや機能拡張ユニットに備わるCPUなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【００６０】
本発明を上記記憶媒体に適用する場合、その記憶媒体には、先に説明した（図２，及び／又は図３、及び／又は図４、及び／又は図９、及び／又は図１０に示す）フローチャートに対応するプログラムコードが格納されることになる。
【００６１】
【発明の効果】
以上の説明により、本発明によれば、領域毎に効果的に周波数帯毎の画像成分を変換することができる。
【図面の簡単な説明】
【図１】本発明の第１の実施形態におけるＸ線撮影装置の構成を示す図である。
【図２】本発明の第１の実施形態におけるＸ線撮影装置が行う処理のフローチャートである。
【図３】本発明の第１の実施形態による領域分割の処理手順を示すフローチャートである。
【図４】本発明の第１の実施形態による領域分割の処理手順を示すフローチャートである。
【図５】領域分割の説明図である。
【図６】ＤＷＴ変換回路を示す図である。
【図７】変換曲線と閾値を示す図である。
【図８】本発明の第２の実施形態におけるＸ線撮影装置の構成を示す図である。
【図９】本発明の第２の実施形態におけるＸ線撮影装置が行う処理のフローチャートである。
【図１０】本発明の第２の実施形態におけるＸ線撮影装置が行う処理において、ステップＳ９０１の処理を詳細に説明するフローチャートである。
【図１】
0004428868.tif 000002
【図２】
0004428868.tif 000003
【図３】
0004428868.tif 000004
【図４】
0004428868.tif 000005
【図５】
0004428868.tif 000006
【図６】
0004428868.tif 000007
【図７】
0004428868.tif 000008
【図８】
0004428868.tif 000009
【図９】
0004428868.tif 000010
【図１０】
0004428868.tif 000011