画像処理装置、画像処理方法
【課題】 同一フレーム内の画像信号と時間軸方向の変化を利用しつつ、より高い効果のある画質改善が可能な画像処理技術を提供すること。
【解決手段】 第(t−1)フレーム画像について分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を用いて、第tフレーム画像について分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を変更する(S804)。そして、変更した周波数帯域を含む全ての周波数帯域の係数に基づいて画像を再構成し(S807)、再構成された画像を、第tフレーム画像として出力する(S808)。
【解決手段】 第(t−1)フレーム画像について分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を用いて、第tフレーム画像について分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を変更する(S804)。そして、変更した周波数帯域を含む全ての周波数帯域の係数に基づいて画像を再構成し(S807)、再構成された画像を、第tフレーム画像として出力する(S808)。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は複数のフレームから構成される動画像シーケンスを対象とし、対象となる画像を高画質化処理する為の画像処理装置に関する。
【背景技術】
【0002】
従来から、複数フレームで構成される動画像シーケンスのS/N比を改善する画像処理装置として、様々な考案がなされている。S/N比を改善する例としてノイズ除去処理を対象とした場合は、例えば特許文献1に開示されるように、リカーシブフィルタを備えたものがある。この文献においては、リカーシブフィルタのフィルタ係数を画像信号の大きさに応じて変化させることにより、X線画像の特性である画像の濃度とノイズ量との関係に基づいたノイズ除去処理を行うことができる。
【0003】
また、画像を鮮鋭化してより見やすくするための処理としては、例えば特許文献2に開示されるように、ボケマスク処理における強調係数を画像に描出された被写体に応じて変化させる。これにより、不必要なノイズの強調を抑制しつつ画像信号の強調を行うことができる。
【特許文献1】特開平05−087394号公報
【特許文献2】特開平7−21364号公報
【発明の開示】
【発明が解決しようとする課題】
【0004】
前述した従来技術のリカーシブフィルタにおいては、複数のフレームにまたがる画像信号の時間方向加算によりノイズを減少させることができるが、同一フレーム内の情報を有効に活用しているものではない。また、ボケマスク処理においては、画像のエッジ近傍においてオーバーシュートなどのアーチファクトが生ずることが知られている。
【0005】
本発明は以上の問題に鑑みて成されたものであり、同一フレーム内の画像信号と時間軸方向の変化を利用しつつ、より高い効果のある画質改善が可能な画像処理技術を提供することを目的とする。
【課題を解決するための手段】
【0006】
本発明の目的を達成するために、例えば、本発明の画像処理装置は以下の構成を備える。
【0007】
即ち、動画像を構成する各フレームの画像を順次入力する入力手段と、
第1フレーム画像を複数の周波数帯域に分解する分解手段と、
前記第1フレーム画像について前記分解手段が分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を、前記第1フレーム画像よりも前のフレームの画像である第2フレーム画像について前記分解手段が分解した複数の周波数帯域のうち当該所望の周波数帯域に対応する周波数帯域に属する係数を用いて変更する変更手段と、
前記変更手段により変更した周波数帯域を含む全ての周波数帯域の係数に基づいて前記分解手段が行った逆の処理を行うことにより、画像を再構成する再構成手段と、
前記再構成手段により再構成された画像を、前記第1フレーム画像として出力する出力手段と
を備えることを特徴とする。
【0008】
本発明の目的を達成するために、例えば、本発明の画像処理方法は以下の構成を備える。
【0009】
即ち、動画像を構成する各フレームの画像を順次入力する入力工程と、
第1フレーム画像を複数の周波数帯域に分解する分解工程と、
前記第1フレーム画像について前記分解工程で分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を、前記第1フレーム画像よりも前のフレームの画像である第2フレーム画像について前記分解工程で分解した複数の周波数帯域のうち当該所望の周波数帯域に対応する周波数帯域に属する係数を用いて変更する変更工程と、
前記変更工程で変更した周波数帯域を含む全ての周波数帯域の係数に基づいて前記分解工程で行った逆の処理を行うことにより、画像を再構成する再構成工程と、
前記再構成工程で再構成された画像を、前記第1フレーム画像として出力する出力工程と
を備えることを特徴とする。
【発明の効果】
【0010】
本発明の構成により、同一フレーム内の画像信号と時間軸方向の変化を利用しつつ、より高い効果のある画質改善を計ることができる。
【発明を実施するための最良の形態】
【0011】
以下添付図面を参照して、本発明を好適な実施形態に従って詳細に説明する。
【0012】
[第1の実施形態]
図1は、本実施形態に係る画像処理装置100の機能構成を示すブロック図である。同図に示す如く、画像処理装置100は、画像入力部1、周波数分解部2,係数変更部3、画像再構成部4、係数記憶部5、画像出力部6により構成されている。
【0013】
画像入力部1は、動画像を構成する各フレームの画像(画像信号)を入力する。画像の供給源としては例えばディジタルカメラがある。例えば、このディジタルカメラがCCDカメラであり、撮像した画像をアナログ画像信号として出力するものである場合には、この画像入力部1は、このアナログ画像信号に対してA/D変換を施し、ディジタル画像信号として後段の周波数分解部2に送出する。また、画像入力部1は入力された画像信号の形式を変換する以外に、例えば、この画像信号に対して、ディジタルカメラ固有の補正処理を行うようにしても良い。このような補正処理には例えば、ディジタルカメラがCCDカメラである場合にはダーク補正処理などがある。
【0014】
また、画像入力部1には、画像信号は、フレーム単位で、且つスキャンライン順に入力される。即ち、プログレッシブスキャンであることを前提とする。しかし、画像信号の供給源によっては、このような順序でもって画像信号が入力されない場合がある。そこでこのような場合には、画像入力部1は、入力した画像信号の出力順序を入れ替えるようにしても良い。例えば、ディジタルカメラから出力された画像信号がインターレーススキャンである場合、画像入力部1はこの画像信号を受け、プログレッシブスキャンに変換して出力する。
【0015】
本実施形態では、画像入力部1に入力される動画像は、所定のフレームレート(例えば1秒間に30フレーム)で、且つ各フレームの画像の解像度が所定の解像度(例えば縦横1024画素)であるとする。
【0016】
また、本実施形態では、各フレームの画像は色成分が1つの画像、即ち、モノクロ画像であるとしているが、これに限定するものではなく、各フレームの画像が複数の色成分を有しても良い。その場合、後述する各処理は各色成分に対して行われることになる。
【0017】
更に、各フレームの画像が複数の色成分を有する場合、画像入力部1は、色空間変換を行うようにしても良い。例えば、画像入力部1に入力される各フレームの画像がRGBの3つの色成分を有する場合、画像入力部1は、RGBからYUVに色空間を変換し、後続の処理については変換後の各色成分に対して行うようにすればよい。
【0018】
周波数分解部2には画像入力部1からフレーム単位で画像が入力されるので、周波数分解部2は、フレーム毎に画像を複数の周波数帯域に分解する。そして、それぞれのフレームについて分解した周波数帯域の係数をフレーム毎に係数変更部3に送出する。
【0019】
画像を複数の周波数帯域に分解する手法については図2に示す何れの方法を用いても良い。図2(a)は、あるフレームの画像Fnに対して2次元離散ウェーブレット変換処理を施すことで、10個の周波数帯域(LL〜HH1)を得る様子を説明する図である。同図では、3レベルの分解が行われている。なお離散ウェーブレット変換については、例えばS.Mallat, “A Theory of Multiresolution Signal Decomposition: the Wavelet Representation”, IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol 11, No.7, July 1989, 674-693等の文献に詳述されているため、説明は省略する。
【0020】
また、図2(b)は、あるフレームの画像Fnに対してラプラシアンピラミッド分解を行うことで、3つの周波数帯域H1、H2、Lを得る様子を説明する図である。ここでH1およびH2は異なる周波数帯域に係る高周波成分であり、Lは残余分の低周波帯域に係る成分を表す係数群である。H2およびLは、画像Fnに対して解像度が水平および垂直方向の各々について、1/2、1/4となっている。ラプラシアンピラミッド分解については、例えばP.Burt and H.Adelson, “The Laplacian Pyramid as a Compact Image Code”, IEEE Transactions on Communications, Vol.Com-31, No.4, April 1983, 532-540等の文献に詳述されているため説明は省略する。
【0021】
係数変更部3は、周波数分解部2から出力されたフレームtの画像(以下、第tフレーム画像と呼称する)の周波数帯域kに係る係数Cijt(k)を、係数C’ijt(k)に変更する処理を行う。ここで、kは周波数帯域を示す変数で、周波数分解部2が図2(a)に示す如く、画像を10個の周波数帯域に分解した場合には、k=LL〜HH1となる。また、i,jは、それぞれの周波数帯域における係数固有のインデックスである。従って、Cijt(k)は、フレームtにおける画像の周波数帯域kに係る係数群を2次元的に配置した場合に、(i、j)の位置に位置する係数を示すことになる。
【0022】
ここで、係数変更部3が行う係数変更処理について説明する。
【0023】
上述の通り、係数変更部3は、周波数分解部2から出力されたフレームの画像の各周波数帯域に係る係数群を変更するのであるが、変更後の係数群は画像再構成部4に出力すると共に、係数記憶部5に格納する。従って、係数変更部3が、第tフレーム画像の周波数帯域kに係る係数Cijt(k)について変更処理を行う場合、係数記憶部5には既に、第(t−1)フレーム画像の周波数帯域kに係る係数(係数変更部3により変更済み)C’ij(t−1)(k)が格納されている。
【0024】
従って、係数変更部3は、Cijt(k)、C’ij(t−1)(k)を用いて、Cijt(k)を以下のように変更する。
【0025】
C’ijt(k)=Ek{αCijt(k)+(1−α)C’ij(t−1)(k)} d<T
C’ijt(k)=EkCijt(k) d≧T
(式1)
但し、α、Tは、ノイズ除去に係るパラメータであり、画像信号に含まれるノイズ量に応じて最適な値が予め決定され、係数変更部3内に予め保持されている。
【0026】
ここで、αは、変更後の係数を生成するにあたり、現在と過去のフレームに係る係数の混合比率を表している。本実施形態の場合、フレームtにおける係数とフレーム(t−1)における係数とをαで決まる合成比率でもって合成することになる。
【0027】
Tは、係数の差の大きいものを動きのある部分として除外することで、過去のフレームの係数の影響をなくし、不要な残像効果を除くためのパラメータである。従ってdは、以下の式でもって表される。
【0028】
d=|Cijt(k)−C’ij(t−1)(k)|
このdは、係数変更部3により計算される。また、Ekは、周波数帯域kに対する強調係数で、必要な強調の強さに応じて予め設定され、同じく係数変更部3内に予め保持されている。
【0029】
係数変更部3は、以上説明した係数変更処理を、全ての周波数帯域について行う。例えば、周波数分解部2が図2(a)に示す如く、画像を10個の周波数帯域に分解した場合には、k=LL〜HH1について以上説明した係数変更処理を行う。
【0030】
以上のようにして、係数変更部3は、現在のフレームにおけるそれぞれの周波数帯域に属する係数群を、1つ前のフレームにおける変更後の係数群でもって変更する処理を行う。なお、上述の通り、本実施形態では現在のフレームにおける係数群を変更する際には1つ前のフレームにおける変更後の係数群を用いているが、2つ以上前のフレームにおける変更後の係数群を用いるようにしても良い。
【0031】
そして係数変更部3は、変更した係数群を画像再構成部4に出力すると共に、次のフレームにおける係数群を変更する際に用いるために係数記憶部5に格納する。
【0032】
なお、上記(式1)による計算処理を行う前に、以下の式に従って、係数Cijt(k)の縮退処理を行うようにしても良い。
【0033】
Cijt(k)=Cijt(k) Cijt(k)≧Ts
Cijt(k)=0 Cijt(k)<Ts
(式3)
また、Cijt(k)の値を、閾値Tsに満たないからといっていきなり0にするのではなく、Cijt(k)の値が閾値Tsよりも小さい場合に、その値に応じた減衰率でもって減少させるようにしても良い。
【0034】
画像再構成部4は、1フレーム分の画像のそれぞれの周波数帯域に属する変更後の係数群を係数変更部3から受けると、これに対して周波数分解部2が行った処理とは逆の処理を行い、画像を再構成する。例えば、周波数分解部2が2次元離散ウェーブレット変換処理を行ったのであれば、画像再構成部4は、2次元離散逆ウェーブレット変換処理を行う。
【0035】
画像出力部6は、画像再構成部4により得られた再構成画像を出力する。
【0036】
以上説明した各部により、画像入力部1に入力された各フレームの画像において各周波数帯域における係数群を、1つ前のフレームの画像において対応する周波数帯域の係数群を用いて修正することができる。そして修正した係数群に基づいて再構成された画像を順次、画像出力部6により出力することができる。
【0037】
なお、画像出力部6による出力先については特に限定するものではないが、出力先の装置によっては、出力画像のフォーマットを変更する必要がある。そこで、そのような場合には、画像出力部6は、画像のフォーマットを、出力先の装置に応じたフォーマットに変更してから出力する。
【0038】
ここで、図1に示した各部は専用のハードウェアでもって構成し、その結果、画像処理装置100を1つのチップとして構成するようにしても良い。本実施形態では係数記憶部5以外をソフトウェアでもって構成し、係数記憶部5はRAMなどのメモリでもって構成する。この場合、図1に示した画像処理装置は、一般のPC(パーソナルコンピュータ)やWS(ワークステーション)等のコンピュータでもって構成することができる。
【0039】
図3は、画像処理装置として適用するコンピュータのハードウェア構成、及びその周辺機器を示すブロック図である。
【0040】
同図に示す如く、コンピュータ300は、ネットワーク350を介して撮像装置360やファイルサーバ370に接続されており、互いにデータ通信が可能な構成となっている。
【0041】
先ず、コンピュータ300について説明する。
【0042】
301はCPUで、RAM302やROM303に格納されているプログラムやデータを用いてコンピュータ300全体の制御を行うと共に、コンピュータ300を適用した上記画像処理装置100が行う各処理を実行する。
【0043】
302はRAMで、光磁気ディスク305やハードディスク304からロードされたプログラムやデータを一時的に記憶するためのエリアを備える。更にRAM302は、ネットワーク350を介して撮像装置360やファイルサーバ370からダウンロードした画像データ等を一時的に記憶するためのエリアを備える。また、RAM302は、CPU301が各種の処理を実行する際に用いるワークエリアも備える。
【0044】
303はROMで、コンピュータ300の設定データや、ブートプログラムなどを格納する。
【0045】
304はハードディスクで、OS(オペレーティングシステム)や、図1に示した画像入力部1、周波数分解部2,係数変更部3、画像再構成部4、画像出力部6の各部が行う上記各処理をCPU301に実行させるためのプログラムやデータを保持している。そしてこれらはCPU301による制御に従って適宜RAM302にロードされ、CPU301による処理対象となる。また、上記動画像のデータをこのハードディスク304に保存させておくようにしても良い。
【0046】
305は光磁気ディスク305で、情報記憶媒体としての一例であり、ハードディスク304に保存されているプログラムやデータの一部若しくは全部をこの光磁気ディスク305に格納するようにしても良い。
【0047】
306、307はそれぞれマウス、キーボードで、コンピュータ300の操作者が操作することで、各種の指示をCPU301に対して入力することができる。
【0048】
308はプリンタで、画像出力部6による出力先の一例である。
【0049】
309は表示装置で、CRTや液晶画面などにより構成されており、CPU301による処理結果を画像や文字などでもって表示する。例えば、図1に示した各部により処理され、最終的に画像出力部6から出力された各フレームの画像を表示する。
【0050】
310は上述の各部を繋ぐバスである。
【0051】
次に、撮像装置360について説明する。撮像装置360は、例えばディジタルカメラなど、動画像を撮像するものであり、撮像した動画像のデータはネットワーク350を介してコンピュータ300に供給される。なお、動画像のデータは複数フレーム分をまとめてコンピュータ300に供給するようにしても良いし、撮像したフレーム毎に供給するようにしても良い。
【0052】
次に、ファイルサーバ370について説明する。ファイルサーバ370には画像データを含む様々な情報が格納可能であり、本実施形態の場合には、図1に示した各部によって行った上記各処理により得られる各フレームの画像を格納することができる。従って、コンピュータ300は、各フレームの画像について図1に示した各部による処理を行うと、処理済みの各フレームの画像をネットワーク350を介してファイルサーバ370に転送する。これにより、ファイルサーバ370には、コンピュータ300による処理済みの動画像データを登録することができる。
【0053】
次にネットワーク350について説明する。ネットワーク350は、LANやインターネット等のネットワークにより構成されている。なお、このネットワーク350は、無線、有線の何れで構成しても良く、また双方を部分的に組み合わせるようにしても良い。
【0054】
図8は、図1に示した各部が行う処理のフローチャートである。なお、同図のフローチャートに従った処理をコンピュータ300のCPU301に実行させるためのプログラムやデータはハードディスク304や光磁気ディスク305に格納されており、適宜CPU301によりRAM302にロードされる。そして、CPU301がこれを用いて処理を実行することで、コンピュータ300は以下説明する各処理を実行する。
【0055】
なお、同図のフローチャートは、1つのフレーム(第tフレーム)の画像についてのものであるので、複数フレームのそれぞれについて同様の処理を行う場合には、同図のフローチャートに従った処理を各フレームについて行えばよい。
【0056】
先ず、CPU301は周波数分解部2として機能し、RAM302に保持されている第tフレーム画像を複数の周波数帯域に分解する処理を行う(ステップS802)。次にステップS803からステップS806までの処理を分解されたそれぞれの周波数帯域について行う。
【0057】
先ず、ハードディスク304内、若しくはRAM302内に保持している第(t−1)フレーム画像の周波数帯域kに係る係数C’ij(t−1)(k)を読み出す(ステップS803)。そして、CPU301は係数変更部3として機能し、読み出したこの係数を用いて、第tフレーム画像の周波数帯域kに係る係数Cijt(k)を上記(式1)に従って変更する処理を行う(ステップS804)。なお、上述の通り、本ステップの前段で上記(式3)に従った処理等でもって、係数の縮退処理を行うようにしても良い。
【0058】
次に、CPU301は、第tフレーム画像の周波数帯域kに属する全ての係数について変更処理を行ったかを判断し(ステップS805)、行っていない場合には、処理をステップS803に戻す。そして、未だ変更を行っていない第tフレーム画像の周波数帯域kに係る係数Cijt(k)についてステップS803以降の処理を行う。
【0059】
一方、第tフレーム画像の周波数帯域kに属する全ての係数について変更処理を行った場合には処理をステップS806に進める。そして、処理の対象となる全ての周波数帯域(処理対象となる全ての周波数帯域についての変数k)について、ステップS804における変更処理を行ったか否かを判断する(ステップS806)。そして、行っていない場合には処理をステップS803に戻し、第tフレーム画像の複数周波数帯域のうち、未だ処理対象となっていない周波数帯域についてステップS803以降の処理を行う。
【0060】
ここで、処理の対象となる周波数帯域については予め決定されており、処理対象となる周波数帯域を示す情報は予めRAM302に保持されている。処理対象となる周波数帯域は、ノイズ除去を効果的に行うために、高周波成分を多く含む帯域とすることが望ましい。そのために、例えば図2(a)に示す離散ウェーブレット変換の場合には、LLを除く全ての周波数帯域を処理対象となる周波数帯域とすれば良い。また、図2(b)に示すラプラシアンピラミッドの場合には、Lを除く全ての周波数帯域を処理対象となる周波数帯域とすれば良い。
【0061】
一方、全ての対象となる周波数帯域について、ステップS804における変更処理を行った場合には処理をステップS807に進める。そして、CPU301は画像再構成部4として機能し、第tフレーム画像の複数の変更済みの周波数帯域に対して、ステップS802で行った処理の逆処理を行うことで、画像を再構成する(ステップS807)。そしてCPU301は画像出力部6として機能し、再構成された画像を、第tフレーム画像としてハードディスク304やRAM302、ファイルサーバ370等に出力する(ステップS808)。また、表示装置309に出力するようにしても良い。
【0062】
以上の説明により、本実施形態によれば、複数フレームから構成される動画像に対して周波数空間上の処理を効果的に行うことでき、より高いノイズ除去効果および鮮鋭化処理を実現することが可能となる。
【0063】
[第2の実施形態]
本実施形態では、第1の実施形態に係る処理において、画像入力部1に入力された画像に関する情報を更に入力する。これにより、入力される画像の種別や撮影条件などに応じて動画像を処理することができる。
【0064】
図4は、本実施形態に係る画像処理装置400の機能構成を示すブロック図である。同図において図1と同じ部分については同じ番号を付けており、その説明は省略する。同図に示す如く、本実施形態に係る画像処理装置400は、図1に示した第1の実施形態に係る画像処理装置100に撮影条件入力部7を加えた構成を有する。
【0065】
撮影条件入力部7は、画像入力部1に入力される動画像の撮影条件等、動画像に関する情報を入力するものである。「動画像に関する情報」については特に限定するものではないが、本実施形態では、入力画像が医用X線画像である場合について説明する。
【0066】
処理対象を医用X線動画像とした場合、撮影条件としては撮影部位、X線照射量、あるいは撮影方法などがある。これらの情報は、撮影条件入力部7に対して入力される。
【0067】
そして、例えば、撮影条件入力部7は、以下に列挙する撮影方法に応じて、周波数分解部2、画像再構成部4における周波数分解レベルを決定する。
【0068】
(1) 透視
(2) DSA
(3) 静止画撮影
例えば、撮影条件入力部7に「撮影方法が透視である旨を示す情報」が入力された場合には、周波数分解部2、画像再構成部4における周波数分解レベルを3と決定する。また、撮影条件入力部7に「撮影方法がDSAである旨を示す情報」が入力された場合には、周波数分解部2、画像再構成部4における周波数分解レベルを5と決定する。また、撮影条件入力部7に「撮影方法が静止画撮影である旨を示す情報」が入力された場合には、周波数分解部2、画像再構成部4における周波数分解レベルを7と決定する。
【0069】
周波数分解レベルを多く取ることで、より細かい周波数処理が可能となるが、透視撮影においてはその用途が診断ではなく、カテーテルなどの手術デバイスの位置を確認できれば十分である。しかし、表示の遅れを回避するためには、フレームレートは高く取る必要がある。そこで、透視に対しては分解レベルを低く設定することにより計算量を抑え、透視に最適な処理を実現することができる。
【0070】
一方、DSAにおいてフレームレートは例えば毎秒7フレーム程度であるが、微細な血管構造までを描出する必要性から、より高い分解レベルで処理を行うことが望ましい。さらに静止画撮影においては、診断に耐えうる画質とするために計算量が増大したとしても可能な限り高画質を達成する為、高い分解レベルとする必要がある。
【0071】
このように、撮影の方法により周波数分解のレベルを変化させることで、目的に応じた処理を実現することができる。
【0072】
また分解レベルのみではなく、係数変更部3における処理のパラメータを変更するようにしても良い。例えば、上記(式1)におけるαの値を透視の場合は小さく、DSAの場合は大きくなるように変更する。また、透視の場合は、被爆量を抑えるためにX線照射量を小さくするが、これにより画像におけるノイズ量は多くなる。このため前フレームの係数に対する重みを大きく取ることによって、よりノイズ除去効果を高めることができるようになる。
【0073】
また、処理対象周波数帯域の選択を、撮影条件入力部7への入力に基づいて、最も高い周波数帯域から所定の数の帯域についてのみ処理を行うようにしてもよい。特に、透視撮影など画像中に映出された被写体に動きがあるような場合においては、低い周波数帯域に対して上述した処理を行うと残像が発生することがある。したがってこのような撮影においては、処理対象となる周波数帯域として所定数の高周波帯域が選択すればよい。
【0074】
例えば、離散ウェーブレット変換の場合であれば、LLおよびHL3、LH3、HH3を除く帯域を処理対象から除くことによって、残像の影響を抑えながらノイズ除去を行うことができる。
【0075】
[第3の実施形態]
第1の実施形態では、係数変更部3により変更した係数C’ijt(k)を係数記憶部5に格納するようにしたが、変更前の係数Cijt(k)を格納するようにしても良い。その場合、係数変更部3は、上記(式1)の代わりに、以下の(式4)を用いて、係数を変更する処理を行う。
【0076】
C’ijt(k)=Ek{αCijt(k)+(1−α)Cij(t−1)(k)} d<T
C’ijt(k)=EkCijt(k) d≧T
(式4)
即ち、第1の実施形態において係数C’ij(t−1)(k)の代わりに、係数Cij(t−1)(k)を用いるようにしても良い。
【0077】
更に、係数変更部3において上記(式1)に基づいて係数を変更するのか、それとも上記(式4)を用いて係数を変更するのかを、図4に示した撮影条件入力部7に入力される情報に基づいて決定するようにしても良い。
【0078】
例えば、画像入力部1に入力される画像中の被写体が微細な構造を持ち、更にある程度の動きがある場合は、(式1)に基づくリカーシブなフィルタリングを用いると残像効果が大きくなるため観察上好ましくない。しかし、(式4)を用いればその影響を抑制することができる。
【0079】
従って、画像入力部1に入力される画像が医用画像であって、この画像の撮影方式が血管造影などの場合には以下のようにすればよい。即ち、撮影条件入力部7が係数変更部3における処理のパラメータを切り替えると同時に、(式1)を用いる場合には係数記憶部5に対する入力を係数変更部3による出力とし、(式4)を用いる場合には周波数分解部2による出力とする。これにより、被写体および撮影方法に最適な処理を選択することができるようになる。
【0080】
更に、係数記憶部5には、現在係数変更部3によって処理対象となっているフレームの過去複数フレームの係数を保持させるようにしても良く、その場合、(式1)や(式4)において重み付け係数αは、過去のフレーム毎に変化させるようにしてもよい。
【0081】
[第4の実施形態]
第1の実施形態では、図8に示したフローチャートに従った処理は、1つのフレーム画像全体について行うものとして説明した。しかし、図5に示す如く、1つのフレーム画像を複数の矩形(矩形画像)に分割し、それぞれの矩形毎に図8に示したフローチャートに従った処理を行うようにしても良い。図5は、1つのフレーム画像を16個の矩形画像(T1〜T16)に分割した様子を示す図である。
【0082】
この場合、第1の実施形態で説明に用いた第tフレーム画像は、第tフレームにおけるフレーム画像を分割した1つの矩形画像(第t矩形画像)に相当する。また、第(t−1)フレーム画像は、第(t−1)フレームにおけるフレーム画像において、この第t矩形画像に位置的に対応する矩形画像(第(t−1)矩形画像)に相当する。
【0083】
このように構成することで、周波数分解部2において周波数分解を行う際の所要メモリ量を低減することができるようになる。さらに、矩形画像毎に処理条件を変更するようにしてもよい。
【0084】
図6は、本実施形態に係る画像処理装置600の機能構成を示すブロック図である。なお、図1,4と同じ部分については同じ番号を付けており、その説明を省略する。図6に示す如く、本実施形態に係る画像処理装置600は、図4に示した画像処理装置400に更に、画像解析部8を加えた構成となっている。
【0085】
このような構成において、画像入力部1は、外部から入力されたフレーム画像を複数の矩形画像に分割し、矩形画像毎に出力する。画像解析部8は、矩形画像を受けると、この矩形画像を構成する全ての画素の画素値を用いて平均画素値、及び画素値の分散値を求める。そして求めたこれらの統計値に基づいて、この矩形画像中に含まれているノイズの量を判断し(例えば、ノイズの量を多い、普通、少ないの3段階でもって判断し)、判断した結果を係数変更部3に通知する。
【0086】
フレーム画像が医用X線画像の場合、X線吸収の多い臓器が映出された矩形画像については、平均画素値は低くなり、またノイズ量も多くなる為、より効果的にノイズ除去を行う必要がある。従って、ノイズ量が多いと画像解析部8により判断された矩形画像については、係数変更部3における処理パラメータを第2の実施形態で説明したように変更することで、より効果的な処理を行うことが可能となる。
【0087】
[第5の実施形態]
図2に示すような撮像装置360が複数のユニットでもって構成されている場合には、上述の実施形態で説明したような画像処理装置をそれぞれのユニットに接続するようにしても良い。この場合、画像処理装置は1つのチップでもって構成することが好ましい。そして、それぞれのユニットに接続された画像処理装置の出力を最後に合成し、最終的な出力を得る。
【0088】
図7(a)は、上述の実施形態で説明したような画像処理装置を、撮像装置を構成している各ユニットに接続した場合の、この撮像装置の概略構成を示す図である。同図において701,702は平面センサを構成しているユニット(それぞれU1,2)である。即ち、この平面センサは、2つのユニット701,702を張り合わせた構成を有する。
【0089】
711,712はそれぞれ画像処理装置(それぞれIP1,2)で、例えば図6に示すような構成を有する。それぞれの画像処理装置711,712は、ユニット701、702からの画像信号を受け、上述の如くフレーム画像に対する処理を行う。この構成において各画像処理装置711,712は並列的に処理を行い、処理結果を画像合成部9に出力する。
【0090】
画像合成部9は、それぞれの画像処理装置711,712からの出力を合成し、最終的に1つのフレーム画像を出力する。
【0091】
このようにセンサが複数のユニットから構成されている場合は、処理を並列的に行えるように画像処理装置をそれぞれのユニットに接続することで、より高い処理効率を得ることが可能となる。
【0092】
また、図7(b)に示す如く、画像処理装置711に対して複数個のユニット(U1〜U4)、画像処理装置712に対して複数個のユニット(U5〜U8)を接続するような構成であっても良い。この場合、画像処理装置711,712それぞれに対して1つのユニットが順次選択され、選択されたユニットからの画像信号を受けることになる。図7(b)は、上述の実施形態で説明したような画像処理装置に対して複数個のユニットを接続した場合の、撮像装置の概略構成を示す図である。
【0093】
また、上記各実施形態は適宜組み合わせても良い。
【0094】
以上の各実施形態によれば、動画像を処理する場合に、同一フレーム内の画像信号と時間軸方向の変化を利用しつつ、より高い効果のある画質改善を行うことができる。
【0095】
[その他の実施形態]
また、本発明の目的は、以下のようにすることによって達成されることはいうまでもない。即ち、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記録媒体(または記憶媒体)を、システムあるいは装置に供給する。そして、そのシステムあるいは装置のコンピュータ(またはCPUやMPU)が記録媒体に格納されたプログラムコードを読み出し実行する。この場合、記録媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記録した記録媒体は本発明を構成することになる。
【0096】
また、コンピュータが読み出したプログラムコードを実行することにより、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているオペレーティングシステム(OS)などが実際の処理の一部または全部を行う。その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【0097】
さらに、記録媒体から読み出されたプログラムコードが、コンピュータに挿入された機能拡張カードやコンピュータに接続された機能拡張ユニットに備わるメモリに書込まれたとする。その後、そのプログラムコードの指示に基づき、その機能拡張カードや機能拡張ユニットに備わるCPUなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【0098】
本発明を上記記録媒体に適用する場合、その記録媒体には、先に説明したフローチャートに対応するプログラムコードが格納されることになる。
【図面の簡単な説明】
【0099】
【図1】本発明の第1の実施形態に係る画像処理装置100の機能構成を示すブロック図である。
【図2】(a)は、画像Fnに対して2次元離散ウェーブレット変換処理を施すことで10個の周波数帯域を得る様子を説明する図、(b)は、画像Fnに対してラプラシアンピラミッド分解を行うことで3つの周波数帯域を得る様子を説明する図である。
【図3】画像処理装置として適用するコンピュータのハードウェア構成、及びその周辺機器を示すブロック図である。
【図4】本発明の第2の実施形態に係る画像処理装置400の機能構成を示すブロック図である。
【図5】1つのフレーム画像を16個の矩形画像(T1〜T16)に分割した様子を示す図である。
【図6】本発明の第1の実施形態に係る画像処理装置600の機能構成を示すブロック図である。
【図7】(a)は、画像処理装置を、撮像装置を構成している各ユニットに接続した場合の、この撮像装置の概略構成を示す図で、(b)は、画像処理装置に対して複数個のユニットを接続した場合の、撮像装置の概略構成を示す図である。
【図8】図1に示した各部が行う処理のフローチャートである。
【技術分野】
【0001】
本発明は複数のフレームから構成される動画像シーケンスを対象とし、対象となる画像を高画質化処理する為の画像処理装置に関する。
【背景技術】
【0002】
従来から、複数フレームで構成される動画像シーケンスのS/N比を改善する画像処理装置として、様々な考案がなされている。S/N比を改善する例としてノイズ除去処理を対象とした場合は、例えば特許文献1に開示されるように、リカーシブフィルタを備えたものがある。この文献においては、リカーシブフィルタのフィルタ係数を画像信号の大きさに応じて変化させることにより、X線画像の特性である画像の濃度とノイズ量との関係に基づいたノイズ除去処理を行うことができる。
【0003】
また、画像を鮮鋭化してより見やすくするための処理としては、例えば特許文献2に開示されるように、ボケマスク処理における強調係数を画像に描出された被写体に応じて変化させる。これにより、不必要なノイズの強調を抑制しつつ画像信号の強調を行うことができる。
【特許文献1】特開平05−087394号公報
【特許文献2】特開平7−21364号公報
【発明の開示】
【発明が解決しようとする課題】
【0004】
前述した従来技術のリカーシブフィルタにおいては、複数のフレームにまたがる画像信号の時間方向加算によりノイズを減少させることができるが、同一フレーム内の情報を有効に活用しているものではない。また、ボケマスク処理においては、画像のエッジ近傍においてオーバーシュートなどのアーチファクトが生ずることが知られている。
【0005】
本発明は以上の問題に鑑みて成されたものであり、同一フレーム内の画像信号と時間軸方向の変化を利用しつつ、より高い効果のある画質改善が可能な画像処理技術を提供することを目的とする。
【課題を解決するための手段】
【0006】
本発明の目的を達成するために、例えば、本発明の画像処理装置は以下の構成を備える。
【0007】
即ち、動画像を構成する各フレームの画像を順次入力する入力手段と、
第1フレーム画像を複数の周波数帯域に分解する分解手段と、
前記第1フレーム画像について前記分解手段が分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を、前記第1フレーム画像よりも前のフレームの画像である第2フレーム画像について前記分解手段が分解した複数の周波数帯域のうち当該所望の周波数帯域に対応する周波数帯域に属する係数を用いて変更する変更手段と、
前記変更手段により変更した周波数帯域を含む全ての周波数帯域の係数に基づいて前記分解手段が行った逆の処理を行うことにより、画像を再構成する再構成手段と、
前記再構成手段により再構成された画像を、前記第1フレーム画像として出力する出力手段と
を備えることを特徴とする。
【0008】
本発明の目的を達成するために、例えば、本発明の画像処理方法は以下の構成を備える。
【0009】
即ち、動画像を構成する各フレームの画像を順次入力する入力工程と、
第1フレーム画像を複数の周波数帯域に分解する分解工程と、
前記第1フレーム画像について前記分解工程で分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を、前記第1フレーム画像よりも前のフレームの画像である第2フレーム画像について前記分解工程で分解した複数の周波数帯域のうち当該所望の周波数帯域に対応する周波数帯域に属する係数を用いて変更する変更工程と、
前記変更工程で変更した周波数帯域を含む全ての周波数帯域の係数に基づいて前記分解工程で行った逆の処理を行うことにより、画像を再構成する再構成工程と、
前記再構成工程で再構成された画像を、前記第1フレーム画像として出力する出力工程と
を備えることを特徴とする。
【発明の効果】
【0010】
本発明の構成により、同一フレーム内の画像信号と時間軸方向の変化を利用しつつ、より高い効果のある画質改善を計ることができる。
【発明を実施するための最良の形態】
【0011】
以下添付図面を参照して、本発明を好適な実施形態に従って詳細に説明する。
【0012】
[第1の実施形態]
図1は、本実施形態に係る画像処理装置100の機能構成を示すブロック図である。同図に示す如く、画像処理装置100は、画像入力部1、周波数分解部2,係数変更部3、画像再構成部4、係数記憶部5、画像出力部6により構成されている。
【0013】
画像入力部1は、動画像を構成する各フレームの画像(画像信号)を入力する。画像の供給源としては例えばディジタルカメラがある。例えば、このディジタルカメラがCCDカメラであり、撮像した画像をアナログ画像信号として出力するものである場合には、この画像入力部1は、このアナログ画像信号に対してA/D変換を施し、ディジタル画像信号として後段の周波数分解部2に送出する。また、画像入力部1は入力された画像信号の形式を変換する以外に、例えば、この画像信号に対して、ディジタルカメラ固有の補正処理を行うようにしても良い。このような補正処理には例えば、ディジタルカメラがCCDカメラである場合にはダーク補正処理などがある。
【0014】
また、画像入力部1には、画像信号は、フレーム単位で、且つスキャンライン順に入力される。即ち、プログレッシブスキャンであることを前提とする。しかし、画像信号の供給源によっては、このような順序でもって画像信号が入力されない場合がある。そこでこのような場合には、画像入力部1は、入力した画像信号の出力順序を入れ替えるようにしても良い。例えば、ディジタルカメラから出力された画像信号がインターレーススキャンである場合、画像入力部1はこの画像信号を受け、プログレッシブスキャンに変換して出力する。
【0015】
本実施形態では、画像入力部1に入力される動画像は、所定のフレームレート(例えば1秒間に30フレーム)で、且つ各フレームの画像の解像度が所定の解像度(例えば縦横1024画素)であるとする。
【0016】
また、本実施形態では、各フレームの画像は色成分が1つの画像、即ち、モノクロ画像であるとしているが、これに限定するものではなく、各フレームの画像が複数の色成分を有しても良い。その場合、後述する各処理は各色成分に対して行われることになる。
【0017】
更に、各フレームの画像が複数の色成分を有する場合、画像入力部1は、色空間変換を行うようにしても良い。例えば、画像入力部1に入力される各フレームの画像がRGBの3つの色成分を有する場合、画像入力部1は、RGBからYUVに色空間を変換し、後続の処理については変換後の各色成分に対して行うようにすればよい。
【0018】
周波数分解部2には画像入力部1からフレーム単位で画像が入力されるので、周波数分解部2は、フレーム毎に画像を複数の周波数帯域に分解する。そして、それぞれのフレームについて分解した周波数帯域の係数をフレーム毎に係数変更部3に送出する。
【0019】
画像を複数の周波数帯域に分解する手法については図2に示す何れの方法を用いても良い。図2(a)は、あるフレームの画像Fnに対して2次元離散ウェーブレット変換処理を施すことで、10個の周波数帯域(LL〜HH1)を得る様子を説明する図である。同図では、3レベルの分解が行われている。なお離散ウェーブレット変換については、例えばS.Mallat, “A Theory of Multiresolution Signal Decomposition: the Wavelet Representation”, IEEE Transactions on Pattern Analysis and Machine Intelligence, Vol 11, No.7, July 1989, 674-693等の文献に詳述されているため、説明は省略する。
【0020】
また、図2(b)は、あるフレームの画像Fnに対してラプラシアンピラミッド分解を行うことで、3つの周波数帯域H1、H2、Lを得る様子を説明する図である。ここでH1およびH2は異なる周波数帯域に係る高周波成分であり、Lは残余分の低周波帯域に係る成分を表す係数群である。H2およびLは、画像Fnに対して解像度が水平および垂直方向の各々について、1/2、1/4となっている。ラプラシアンピラミッド分解については、例えばP.Burt and H.Adelson, “The Laplacian Pyramid as a Compact Image Code”, IEEE Transactions on Communications, Vol.Com-31, No.4, April 1983, 532-540等の文献に詳述されているため説明は省略する。
【0021】
係数変更部3は、周波数分解部2から出力されたフレームtの画像(以下、第tフレーム画像と呼称する)の周波数帯域kに係る係数Cijt(k)を、係数C’ijt(k)に変更する処理を行う。ここで、kは周波数帯域を示す変数で、周波数分解部2が図2(a)に示す如く、画像を10個の周波数帯域に分解した場合には、k=LL〜HH1となる。また、i,jは、それぞれの周波数帯域における係数固有のインデックスである。従って、Cijt(k)は、フレームtにおける画像の周波数帯域kに係る係数群を2次元的に配置した場合に、(i、j)の位置に位置する係数を示すことになる。
【0022】
ここで、係数変更部3が行う係数変更処理について説明する。
【0023】
上述の通り、係数変更部3は、周波数分解部2から出力されたフレームの画像の各周波数帯域に係る係数群を変更するのであるが、変更後の係数群は画像再構成部4に出力すると共に、係数記憶部5に格納する。従って、係数変更部3が、第tフレーム画像の周波数帯域kに係る係数Cijt(k)について変更処理を行う場合、係数記憶部5には既に、第(t−1)フレーム画像の周波数帯域kに係る係数(係数変更部3により変更済み)C’ij(t−1)(k)が格納されている。
【0024】
従って、係数変更部3は、Cijt(k)、C’ij(t−1)(k)を用いて、Cijt(k)を以下のように変更する。
【0025】
C’ijt(k)=Ek{αCijt(k)+(1−α)C’ij(t−1)(k)} d<T
C’ijt(k)=EkCijt(k) d≧T
(式1)
但し、α、Tは、ノイズ除去に係るパラメータであり、画像信号に含まれるノイズ量に応じて最適な値が予め決定され、係数変更部3内に予め保持されている。
【0026】
ここで、αは、変更後の係数を生成するにあたり、現在と過去のフレームに係る係数の混合比率を表している。本実施形態の場合、フレームtにおける係数とフレーム(t−1)における係数とをαで決まる合成比率でもって合成することになる。
【0027】
Tは、係数の差の大きいものを動きのある部分として除外することで、過去のフレームの係数の影響をなくし、不要な残像効果を除くためのパラメータである。従ってdは、以下の式でもって表される。
【0028】
d=|Cijt(k)−C’ij(t−1)(k)|
このdは、係数変更部3により計算される。また、Ekは、周波数帯域kに対する強調係数で、必要な強調の強さに応じて予め設定され、同じく係数変更部3内に予め保持されている。
【0029】
係数変更部3は、以上説明した係数変更処理を、全ての周波数帯域について行う。例えば、周波数分解部2が図2(a)に示す如く、画像を10個の周波数帯域に分解した場合には、k=LL〜HH1について以上説明した係数変更処理を行う。
【0030】
以上のようにして、係数変更部3は、現在のフレームにおけるそれぞれの周波数帯域に属する係数群を、1つ前のフレームにおける変更後の係数群でもって変更する処理を行う。なお、上述の通り、本実施形態では現在のフレームにおける係数群を変更する際には1つ前のフレームにおける変更後の係数群を用いているが、2つ以上前のフレームにおける変更後の係数群を用いるようにしても良い。
【0031】
そして係数変更部3は、変更した係数群を画像再構成部4に出力すると共に、次のフレームにおける係数群を変更する際に用いるために係数記憶部5に格納する。
【0032】
なお、上記(式1)による計算処理を行う前に、以下の式に従って、係数Cijt(k)の縮退処理を行うようにしても良い。
【0033】
Cijt(k)=Cijt(k) Cijt(k)≧Ts
Cijt(k)=0 Cijt(k)<Ts
(式3)
また、Cijt(k)の値を、閾値Tsに満たないからといっていきなり0にするのではなく、Cijt(k)の値が閾値Tsよりも小さい場合に、その値に応じた減衰率でもって減少させるようにしても良い。
【0034】
画像再構成部4は、1フレーム分の画像のそれぞれの周波数帯域に属する変更後の係数群を係数変更部3から受けると、これに対して周波数分解部2が行った処理とは逆の処理を行い、画像を再構成する。例えば、周波数分解部2が2次元離散ウェーブレット変換処理を行ったのであれば、画像再構成部4は、2次元離散逆ウェーブレット変換処理を行う。
【0035】
画像出力部6は、画像再構成部4により得られた再構成画像を出力する。
【0036】
以上説明した各部により、画像入力部1に入力された各フレームの画像において各周波数帯域における係数群を、1つ前のフレームの画像において対応する周波数帯域の係数群を用いて修正することができる。そして修正した係数群に基づいて再構成された画像を順次、画像出力部6により出力することができる。
【0037】
なお、画像出力部6による出力先については特に限定するものではないが、出力先の装置によっては、出力画像のフォーマットを変更する必要がある。そこで、そのような場合には、画像出力部6は、画像のフォーマットを、出力先の装置に応じたフォーマットに変更してから出力する。
【0038】
ここで、図1に示した各部は専用のハードウェアでもって構成し、その結果、画像処理装置100を1つのチップとして構成するようにしても良い。本実施形態では係数記憶部5以外をソフトウェアでもって構成し、係数記憶部5はRAMなどのメモリでもって構成する。この場合、図1に示した画像処理装置は、一般のPC(パーソナルコンピュータ)やWS(ワークステーション)等のコンピュータでもって構成することができる。
【0039】
図3は、画像処理装置として適用するコンピュータのハードウェア構成、及びその周辺機器を示すブロック図である。
【0040】
同図に示す如く、コンピュータ300は、ネットワーク350を介して撮像装置360やファイルサーバ370に接続されており、互いにデータ通信が可能な構成となっている。
【0041】
先ず、コンピュータ300について説明する。
【0042】
301はCPUで、RAM302やROM303に格納されているプログラムやデータを用いてコンピュータ300全体の制御を行うと共に、コンピュータ300を適用した上記画像処理装置100が行う各処理を実行する。
【0043】
302はRAMで、光磁気ディスク305やハードディスク304からロードされたプログラムやデータを一時的に記憶するためのエリアを備える。更にRAM302は、ネットワーク350を介して撮像装置360やファイルサーバ370からダウンロードした画像データ等を一時的に記憶するためのエリアを備える。また、RAM302は、CPU301が各種の処理を実行する際に用いるワークエリアも備える。
【0044】
303はROMで、コンピュータ300の設定データや、ブートプログラムなどを格納する。
【0045】
304はハードディスクで、OS(オペレーティングシステム)や、図1に示した画像入力部1、周波数分解部2,係数変更部3、画像再構成部4、画像出力部6の各部が行う上記各処理をCPU301に実行させるためのプログラムやデータを保持している。そしてこれらはCPU301による制御に従って適宜RAM302にロードされ、CPU301による処理対象となる。また、上記動画像のデータをこのハードディスク304に保存させておくようにしても良い。
【0046】
305は光磁気ディスク305で、情報記憶媒体としての一例であり、ハードディスク304に保存されているプログラムやデータの一部若しくは全部をこの光磁気ディスク305に格納するようにしても良い。
【0047】
306、307はそれぞれマウス、キーボードで、コンピュータ300の操作者が操作することで、各種の指示をCPU301に対して入力することができる。
【0048】
308はプリンタで、画像出力部6による出力先の一例である。
【0049】
309は表示装置で、CRTや液晶画面などにより構成されており、CPU301による処理結果を画像や文字などでもって表示する。例えば、図1に示した各部により処理され、最終的に画像出力部6から出力された各フレームの画像を表示する。
【0050】
310は上述の各部を繋ぐバスである。
【0051】
次に、撮像装置360について説明する。撮像装置360は、例えばディジタルカメラなど、動画像を撮像するものであり、撮像した動画像のデータはネットワーク350を介してコンピュータ300に供給される。なお、動画像のデータは複数フレーム分をまとめてコンピュータ300に供給するようにしても良いし、撮像したフレーム毎に供給するようにしても良い。
【0052】
次に、ファイルサーバ370について説明する。ファイルサーバ370には画像データを含む様々な情報が格納可能であり、本実施形態の場合には、図1に示した各部によって行った上記各処理により得られる各フレームの画像を格納することができる。従って、コンピュータ300は、各フレームの画像について図1に示した各部による処理を行うと、処理済みの各フレームの画像をネットワーク350を介してファイルサーバ370に転送する。これにより、ファイルサーバ370には、コンピュータ300による処理済みの動画像データを登録することができる。
【0053】
次にネットワーク350について説明する。ネットワーク350は、LANやインターネット等のネットワークにより構成されている。なお、このネットワーク350は、無線、有線の何れで構成しても良く、また双方を部分的に組み合わせるようにしても良い。
【0054】
図8は、図1に示した各部が行う処理のフローチャートである。なお、同図のフローチャートに従った処理をコンピュータ300のCPU301に実行させるためのプログラムやデータはハードディスク304や光磁気ディスク305に格納されており、適宜CPU301によりRAM302にロードされる。そして、CPU301がこれを用いて処理を実行することで、コンピュータ300は以下説明する各処理を実行する。
【0055】
なお、同図のフローチャートは、1つのフレーム(第tフレーム)の画像についてのものであるので、複数フレームのそれぞれについて同様の処理を行う場合には、同図のフローチャートに従った処理を各フレームについて行えばよい。
【0056】
先ず、CPU301は周波数分解部2として機能し、RAM302に保持されている第tフレーム画像を複数の周波数帯域に分解する処理を行う(ステップS802)。次にステップS803からステップS806までの処理を分解されたそれぞれの周波数帯域について行う。
【0057】
先ず、ハードディスク304内、若しくはRAM302内に保持している第(t−1)フレーム画像の周波数帯域kに係る係数C’ij(t−1)(k)を読み出す(ステップS803)。そして、CPU301は係数変更部3として機能し、読み出したこの係数を用いて、第tフレーム画像の周波数帯域kに係る係数Cijt(k)を上記(式1)に従って変更する処理を行う(ステップS804)。なお、上述の通り、本ステップの前段で上記(式3)に従った処理等でもって、係数の縮退処理を行うようにしても良い。
【0058】
次に、CPU301は、第tフレーム画像の周波数帯域kに属する全ての係数について変更処理を行ったかを判断し(ステップS805)、行っていない場合には、処理をステップS803に戻す。そして、未だ変更を行っていない第tフレーム画像の周波数帯域kに係る係数Cijt(k)についてステップS803以降の処理を行う。
【0059】
一方、第tフレーム画像の周波数帯域kに属する全ての係数について変更処理を行った場合には処理をステップS806に進める。そして、処理の対象となる全ての周波数帯域(処理対象となる全ての周波数帯域についての変数k)について、ステップS804における変更処理を行ったか否かを判断する(ステップS806)。そして、行っていない場合には処理をステップS803に戻し、第tフレーム画像の複数周波数帯域のうち、未だ処理対象となっていない周波数帯域についてステップS803以降の処理を行う。
【0060】
ここで、処理の対象となる周波数帯域については予め決定されており、処理対象となる周波数帯域を示す情報は予めRAM302に保持されている。処理対象となる周波数帯域は、ノイズ除去を効果的に行うために、高周波成分を多く含む帯域とすることが望ましい。そのために、例えば図2(a)に示す離散ウェーブレット変換の場合には、LLを除く全ての周波数帯域を処理対象となる周波数帯域とすれば良い。また、図2(b)に示すラプラシアンピラミッドの場合には、Lを除く全ての周波数帯域を処理対象となる周波数帯域とすれば良い。
【0061】
一方、全ての対象となる周波数帯域について、ステップS804における変更処理を行った場合には処理をステップS807に進める。そして、CPU301は画像再構成部4として機能し、第tフレーム画像の複数の変更済みの周波数帯域に対して、ステップS802で行った処理の逆処理を行うことで、画像を再構成する(ステップS807)。そしてCPU301は画像出力部6として機能し、再構成された画像を、第tフレーム画像としてハードディスク304やRAM302、ファイルサーバ370等に出力する(ステップS808)。また、表示装置309に出力するようにしても良い。
【0062】
以上の説明により、本実施形態によれば、複数フレームから構成される動画像に対して周波数空間上の処理を効果的に行うことでき、より高いノイズ除去効果および鮮鋭化処理を実現することが可能となる。
【0063】
[第2の実施形態]
本実施形態では、第1の実施形態に係る処理において、画像入力部1に入力された画像に関する情報を更に入力する。これにより、入力される画像の種別や撮影条件などに応じて動画像を処理することができる。
【0064】
図4は、本実施形態に係る画像処理装置400の機能構成を示すブロック図である。同図において図1と同じ部分については同じ番号を付けており、その説明は省略する。同図に示す如く、本実施形態に係る画像処理装置400は、図1に示した第1の実施形態に係る画像処理装置100に撮影条件入力部7を加えた構成を有する。
【0065】
撮影条件入力部7は、画像入力部1に入力される動画像の撮影条件等、動画像に関する情報を入力するものである。「動画像に関する情報」については特に限定するものではないが、本実施形態では、入力画像が医用X線画像である場合について説明する。
【0066】
処理対象を医用X線動画像とした場合、撮影条件としては撮影部位、X線照射量、あるいは撮影方法などがある。これらの情報は、撮影条件入力部7に対して入力される。
【0067】
そして、例えば、撮影条件入力部7は、以下に列挙する撮影方法に応じて、周波数分解部2、画像再構成部4における周波数分解レベルを決定する。
【0068】
(1) 透視
(2) DSA
(3) 静止画撮影
例えば、撮影条件入力部7に「撮影方法が透視である旨を示す情報」が入力された場合には、周波数分解部2、画像再構成部4における周波数分解レベルを3と決定する。また、撮影条件入力部7に「撮影方法がDSAである旨を示す情報」が入力された場合には、周波数分解部2、画像再構成部4における周波数分解レベルを5と決定する。また、撮影条件入力部7に「撮影方法が静止画撮影である旨を示す情報」が入力された場合には、周波数分解部2、画像再構成部4における周波数分解レベルを7と決定する。
【0069】
周波数分解レベルを多く取ることで、より細かい周波数処理が可能となるが、透視撮影においてはその用途が診断ではなく、カテーテルなどの手術デバイスの位置を確認できれば十分である。しかし、表示の遅れを回避するためには、フレームレートは高く取る必要がある。そこで、透視に対しては分解レベルを低く設定することにより計算量を抑え、透視に最適な処理を実現することができる。
【0070】
一方、DSAにおいてフレームレートは例えば毎秒7フレーム程度であるが、微細な血管構造までを描出する必要性から、より高い分解レベルで処理を行うことが望ましい。さらに静止画撮影においては、診断に耐えうる画質とするために計算量が増大したとしても可能な限り高画質を達成する為、高い分解レベルとする必要がある。
【0071】
このように、撮影の方法により周波数分解のレベルを変化させることで、目的に応じた処理を実現することができる。
【0072】
また分解レベルのみではなく、係数変更部3における処理のパラメータを変更するようにしても良い。例えば、上記(式1)におけるαの値を透視の場合は小さく、DSAの場合は大きくなるように変更する。また、透視の場合は、被爆量を抑えるためにX線照射量を小さくするが、これにより画像におけるノイズ量は多くなる。このため前フレームの係数に対する重みを大きく取ることによって、よりノイズ除去効果を高めることができるようになる。
【0073】
また、処理対象周波数帯域の選択を、撮影条件入力部7への入力に基づいて、最も高い周波数帯域から所定の数の帯域についてのみ処理を行うようにしてもよい。特に、透視撮影など画像中に映出された被写体に動きがあるような場合においては、低い周波数帯域に対して上述した処理を行うと残像が発生することがある。したがってこのような撮影においては、処理対象となる周波数帯域として所定数の高周波帯域が選択すればよい。
【0074】
例えば、離散ウェーブレット変換の場合であれば、LLおよびHL3、LH3、HH3を除く帯域を処理対象から除くことによって、残像の影響を抑えながらノイズ除去を行うことができる。
【0075】
[第3の実施形態]
第1の実施形態では、係数変更部3により変更した係数C’ijt(k)を係数記憶部5に格納するようにしたが、変更前の係数Cijt(k)を格納するようにしても良い。その場合、係数変更部3は、上記(式1)の代わりに、以下の(式4)を用いて、係数を変更する処理を行う。
【0076】
C’ijt(k)=Ek{αCijt(k)+(1−α)Cij(t−1)(k)} d<T
C’ijt(k)=EkCijt(k) d≧T
(式4)
即ち、第1の実施形態において係数C’ij(t−1)(k)の代わりに、係数Cij(t−1)(k)を用いるようにしても良い。
【0077】
更に、係数変更部3において上記(式1)に基づいて係数を変更するのか、それとも上記(式4)を用いて係数を変更するのかを、図4に示した撮影条件入力部7に入力される情報に基づいて決定するようにしても良い。
【0078】
例えば、画像入力部1に入力される画像中の被写体が微細な構造を持ち、更にある程度の動きがある場合は、(式1)に基づくリカーシブなフィルタリングを用いると残像効果が大きくなるため観察上好ましくない。しかし、(式4)を用いればその影響を抑制することができる。
【0079】
従って、画像入力部1に入力される画像が医用画像であって、この画像の撮影方式が血管造影などの場合には以下のようにすればよい。即ち、撮影条件入力部7が係数変更部3における処理のパラメータを切り替えると同時に、(式1)を用いる場合には係数記憶部5に対する入力を係数変更部3による出力とし、(式4)を用いる場合には周波数分解部2による出力とする。これにより、被写体および撮影方法に最適な処理を選択することができるようになる。
【0080】
更に、係数記憶部5には、現在係数変更部3によって処理対象となっているフレームの過去複数フレームの係数を保持させるようにしても良く、その場合、(式1)や(式4)において重み付け係数αは、過去のフレーム毎に変化させるようにしてもよい。
【0081】
[第4の実施形態]
第1の実施形態では、図8に示したフローチャートに従った処理は、1つのフレーム画像全体について行うものとして説明した。しかし、図5に示す如く、1つのフレーム画像を複数の矩形(矩形画像)に分割し、それぞれの矩形毎に図8に示したフローチャートに従った処理を行うようにしても良い。図5は、1つのフレーム画像を16個の矩形画像(T1〜T16)に分割した様子を示す図である。
【0082】
この場合、第1の実施形態で説明に用いた第tフレーム画像は、第tフレームにおけるフレーム画像を分割した1つの矩形画像(第t矩形画像)に相当する。また、第(t−1)フレーム画像は、第(t−1)フレームにおけるフレーム画像において、この第t矩形画像に位置的に対応する矩形画像(第(t−1)矩形画像)に相当する。
【0083】
このように構成することで、周波数分解部2において周波数分解を行う際の所要メモリ量を低減することができるようになる。さらに、矩形画像毎に処理条件を変更するようにしてもよい。
【0084】
図6は、本実施形態に係る画像処理装置600の機能構成を示すブロック図である。なお、図1,4と同じ部分については同じ番号を付けており、その説明を省略する。図6に示す如く、本実施形態に係る画像処理装置600は、図4に示した画像処理装置400に更に、画像解析部8を加えた構成となっている。
【0085】
このような構成において、画像入力部1は、外部から入力されたフレーム画像を複数の矩形画像に分割し、矩形画像毎に出力する。画像解析部8は、矩形画像を受けると、この矩形画像を構成する全ての画素の画素値を用いて平均画素値、及び画素値の分散値を求める。そして求めたこれらの統計値に基づいて、この矩形画像中に含まれているノイズの量を判断し(例えば、ノイズの量を多い、普通、少ないの3段階でもって判断し)、判断した結果を係数変更部3に通知する。
【0086】
フレーム画像が医用X線画像の場合、X線吸収の多い臓器が映出された矩形画像については、平均画素値は低くなり、またノイズ量も多くなる為、より効果的にノイズ除去を行う必要がある。従って、ノイズ量が多いと画像解析部8により判断された矩形画像については、係数変更部3における処理パラメータを第2の実施形態で説明したように変更することで、より効果的な処理を行うことが可能となる。
【0087】
[第5の実施形態]
図2に示すような撮像装置360が複数のユニットでもって構成されている場合には、上述の実施形態で説明したような画像処理装置をそれぞれのユニットに接続するようにしても良い。この場合、画像処理装置は1つのチップでもって構成することが好ましい。そして、それぞれのユニットに接続された画像処理装置の出力を最後に合成し、最終的な出力を得る。
【0088】
図7(a)は、上述の実施形態で説明したような画像処理装置を、撮像装置を構成している各ユニットに接続した場合の、この撮像装置の概略構成を示す図である。同図において701,702は平面センサを構成しているユニット(それぞれU1,2)である。即ち、この平面センサは、2つのユニット701,702を張り合わせた構成を有する。
【0089】
711,712はそれぞれ画像処理装置(それぞれIP1,2)で、例えば図6に示すような構成を有する。それぞれの画像処理装置711,712は、ユニット701、702からの画像信号を受け、上述の如くフレーム画像に対する処理を行う。この構成において各画像処理装置711,712は並列的に処理を行い、処理結果を画像合成部9に出力する。
【0090】
画像合成部9は、それぞれの画像処理装置711,712からの出力を合成し、最終的に1つのフレーム画像を出力する。
【0091】
このようにセンサが複数のユニットから構成されている場合は、処理を並列的に行えるように画像処理装置をそれぞれのユニットに接続することで、より高い処理効率を得ることが可能となる。
【0092】
また、図7(b)に示す如く、画像処理装置711に対して複数個のユニット(U1〜U4)、画像処理装置712に対して複数個のユニット(U5〜U8)を接続するような構成であっても良い。この場合、画像処理装置711,712それぞれに対して1つのユニットが順次選択され、選択されたユニットからの画像信号を受けることになる。図7(b)は、上述の実施形態で説明したような画像処理装置に対して複数個のユニットを接続した場合の、撮像装置の概略構成を示す図である。
【0093】
また、上記各実施形態は適宜組み合わせても良い。
【0094】
以上の各実施形態によれば、動画像を処理する場合に、同一フレーム内の画像信号と時間軸方向の変化を利用しつつ、より高い効果のある画質改善を行うことができる。
【0095】
[その他の実施形態]
また、本発明の目的は、以下のようにすることによって達成されることはいうまでもない。即ち、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記録媒体(または記憶媒体)を、システムあるいは装置に供給する。そして、そのシステムあるいは装置のコンピュータ(またはCPUやMPU)が記録媒体に格納されたプログラムコードを読み出し実行する。この場合、記録媒体から読み出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記録した記録媒体は本発明を構成することになる。
【0096】
また、コンピュータが読み出したプログラムコードを実行することにより、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているオペレーティングシステム(OS)などが実際の処理の一部または全部を行う。その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【0097】
さらに、記録媒体から読み出されたプログラムコードが、コンピュータに挿入された機能拡張カードやコンピュータに接続された機能拡張ユニットに備わるメモリに書込まれたとする。その後、そのプログラムコードの指示に基づき、その機能拡張カードや機能拡張ユニットに備わるCPUなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も含まれることは言うまでもない。
【0098】
本発明を上記記録媒体に適用する場合、その記録媒体には、先に説明したフローチャートに対応するプログラムコードが格納されることになる。
【図面の簡単な説明】
【0099】
【図1】本発明の第1の実施形態に係る画像処理装置100の機能構成を示すブロック図である。
【図2】(a)は、画像Fnに対して2次元離散ウェーブレット変換処理を施すことで10個の周波数帯域を得る様子を説明する図、(b)は、画像Fnに対してラプラシアンピラミッド分解を行うことで3つの周波数帯域を得る様子を説明する図である。
【図3】画像処理装置として適用するコンピュータのハードウェア構成、及びその周辺機器を示すブロック図である。
【図4】本発明の第2の実施形態に係る画像処理装置400の機能構成を示すブロック図である。
【図5】1つのフレーム画像を16個の矩形画像(T1〜T16)に分割した様子を示す図である。
【図6】本発明の第1の実施形態に係る画像処理装置600の機能構成を示すブロック図である。
【図7】(a)は、画像処理装置を、撮像装置を構成している各ユニットに接続した場合の、この撮像装置の概略構成を示す図で、(b)は、画像処理装置に対して複数個のユニットを接続した場合の、撮像装置の概略構成を示す図である。
【図8】図1に示した各部が行う処理のフローチャートである。
【特許請求の範囲】
【請求項1】
動画像を構成する各フレームの画像を順次入力する入力手段と、
第1フレーム画像を複数の周波数帯域に分解する分解手段と、
前記第1フレーム画像について前記分解手段が分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を、前記第1フレーム画像よりも前のフレームの画像である第2フレーム画像について前記分解手段が分解した複数の周波数帯域のうち当該所望の周波数帯域に対応する周波数帯域に属する係数を用いて変更する変更手段と、
前記変更手段により変更した周波数帯域を含む全ての周波数帯域の係数に基づいて前記分解手段が行った逆の処理を行うことにより、画像を再構成する再構成手段と、
前記再構成手段により再構成された画像を、前記第1フレーム画像として出力する出力手段と
を備えることを特徴とする画像処理装置。
【請求項2】
前記分解手段は、ウェーブレット変換を用いることで、画像を複数の周波数帯域に分解することを特徴とする請求項1に記載の画像処理装置。
【請求項3】
前記所望の周波数帯域は、LLを除く全ての周波数帯域であることを特徴とする請求項2に記載の画像処理装置。
【請求項4】
前記分解手段は、ラプラシアンピラミッド分解処理を行うことで、画像を複数の周波数帯域に分解することを特徴とする請求項1に記載の画像処理装置。
【請求項5】
前記所望の周波数帯域は、Lを除く全ての周波数帯域であることを特徴とする請求項4に記載の画像処理装置。
【請求項6】
前記変更手段は、
前記第1フレーム画像における注目周波数帯域の係数を、前記第2フレーム画像において当該注目周波数帯域に対応する周波数帯域の係数を用いて変更する処理を、前記第1フレーム画像におけるそれぞれの周波数帯域について行うことを特徴とする請求項1に記載の画像処理装置。
【請求項7】
前記第2フレーム画像について前記分解手段が分解したそれぞれの周波数帯域に属する係数は、前記第1フレーム画像における注目周波数帯域の係数を変更する為に用いる際には既に前記変更手段でもって変更されたものであることを特徴とする請求項1に記載の画像処理装置。
【請求項8】
更に、
前記入力手段に入力される画像に係る情報に基づいて、
前記第2フレーム画像について前記分解手段が分解したそれぞれの周波数帯域に属する係数として、前記第1フレーム画像における注目周波数帯域の係数を変更する為に用いる際には既に前記変更手段でもって変更されたものを用いるか、当該変更のないものを用いるのかを切り替える手段
を備えることを特徴とする請求項1に記載の画像処理装置。
【請求項9】
更に、前記入力手段に入力される画像に係る情報に基づいて、前記分解手段、前記再構成手段による分解レベルを制御する手段を備えることを特徴とする請求項1に記載の画像処理装置。
【請求項10】
更に、前記入力手段に入力される画像に係る情報に基づいて、前記変更手段による変更処理に係るパラメータを制御する手段を備えることを特徴とする請求項1に記載の画像処理装置。
【請求項11】
前記第2フレーム画像は、前記第1フレーム画像よりも前の複数フレーム画像であることを特徴とする請求項1に記載の画像処理装置。
【請求項12】
前記第1フレーム画像は第1フレームにおける画像を複数の矩形に分割することで得られる矩形画像であり、前記第2フレーム画像は、当該第1フレームよりも前のフレームにおける画像を複数の矩形に分割した場合に、当該矩形画像に対応する矩形画像であることを特徴とする請求項1に記載の画像処理装置。
【請求項13】
更に、
前記第1フレーム画像におけるノイズの量を判断する判断手段と、
前記判断手段による判断結果に応じて、前記変更手段による処理に係るパラメータを制御する手段と
を備えることを特徴とする請求項12に記載の画像処理装置。
【請求項14】
動画像を構成する各フレームの画像を順次入力する入力工程と、
第1フレーム画像を複数の周波数帯域に分解する分解工程と、
前記第1フレーム画像について前記分解工程で分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を、前記第1フレーム画像よりも前のフレームの画像である第2フレーム画像について前記分解工程で分解した複数の周波数帯域のうち当該所望の周波数帯域に対応する周波数帯域に属する係数を用いて変更する変更工程と、
前記変更工程で変更した周波数帯域を含む全ての周波数帯域の係数に基づいて前記分解工程で行った逆の処理を行うことにより、画像を再構成する再構成工程と、
前記再構成工程で再構成された画像を、前記第1フレーム画像として出力する出力工程と
を備えることを特徴とする画像処理方法。
【請求項15】
コンピュータに請求項14に記載の画像処理方法を実行させることを特徴とするプログラム。
【請求項16】
請求項15に記載のプログラムを格納したことを特徴とする、コンピュータ読み取り可能な記憶媒体。
【請求項1】
動画像を構成する各フレームの画像を順次入力する入力手段と、
第1フレーム画像を複数の周波数帯域に分解する分解手段と、
前記第1フレーム画像について前記分解手段が分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を、前記第1フレーム画像よりも前のフレームの画像である第2フレーム画像について前記分解手段が分解した複数の周波数帯域のうち当該所望の周波数帯域に対応する周波数帯域に属する係数を用いて変更する変更手段と、
前記変更手段により変更した周波数帯域を含む全ての周波数帯域の係数に基づいて前記分解手段が行った逆の処理を行うことにより、画像を再構成する再構成手段と、
前記再構成手段により再構成された画像を、前記第1フレーム画像として出力する出力手段と
を備えることを特徴とする画像処理装置。
【請求項2】
前記分解手段は、ウェーブレット変換を用いることで、画像を複数の周波数帯域に分解することを特徴とする請求項1に記載の画像処理装置。
【請求項3】
前記所望の周波数帯域は、LLを除く全ての周波数帯域であることを特徴とする請求項2に記載の画像処理装置。
【請求項4】
前記分解手段は、ラプラシアンピラミッド分解処理を行うことで、画像を複数の周波数帯域に分解することを特徴とする請求項1に記載の画像処理装置。
【請求項5】
前記所望の周波数帯域は、Lを除く全ての周波数帯域であることを特徴とする請求項4に記載の画像処理装置。
【請求項6】
前記変更手段は、
前記第1フレーム画像における注目周波数帯域の係数を、前記第2フレーム画像において当該注目周波数帯域に対応する周波数帯域の係数を用いて変更する処理を、前記第1フレーム画像におけるそれぞれの周波数帯域について行うことを特徴とする請求項1に記載の画像処理装置。
【請求項7】
前記第2フレーム画像について前記分解手段が分解したそれぞれの周波数帯域に属する係数は、前記第1フレーム画像における注目周波数帯域の係数を変更する為に用いる際には既に前記変更手段でもって変更されたものであることを特徴とする請求項1に記載の画像処理装置。
【請求項8】
更に、
前記入力手段に入力される画像に係る情報に基づいて、
前記第2フレーム画像について前記分解手段が分解したそれぞれの周波数帯域に属する係数として、前記第1フレーム画像における注目周波数帯域の係数を変更する為に用いる際には既に前記変更手段でもって変更されたものを用いるか、当該変更のないものを用いるのかを切り替える手段
を備えることを特徴とする請求項1に記載の画像処理装置。
【請求項9】
更に、前記入力手段に入力される画像に係る情報に基づいて、前記分解手段、前記再構成手段による分解レベルを制御する手段を備えることを特徴とする請求項1に記載の画像処理装置。
【請求項10】
更に、前記入力手段に入力される画像に係る情報に基づいて、前記変更手段による変更処理に係るパラメータを制御する手段を備えることを特徴とする請求項1に記載の画像処理装置。
【請求項11】
前記第2フレーム画像は、前記第1フレーム画像よりも前の複数フレーム画像であることを特徴とする請求項1に記載の画像処理装置。
【請求項12】
前記第1フレーム画像は第1フレームにおける画像を複数の矩形に分割することで得られる矩形画像であり、前記第2フレーム画像は、当該第1フレームよりも前のフレームにおける画像を複数の矩形に分割した場合に、当該矩形画像に対応する矩形画像であることを特徴とする請求項1に記載の画像処理装置。
【請求項13】
更に、
前記第1フレーム画像におけるノイズの量を判断する判断手段と、
前記判断手段による判断結果に応じて、前記変更手段による処理に係るパラメータを制御する手段と
を備えることを特徴とする請求項12に記載の画像処理装置。
【請求項14】
動画像を構成する各フレームの画像を順次入力する入力工程と、
第1フレーム画像を複数の周波数帯域に分解する分解工程と、
前記第1フレーム画像について前記分解工程で分解した複数の周波数帯域のうち所望の周波数帯域に属する係数を、前記第1フレーム画像よりも前のフレームの画像である第2フレーム画像について前記分解工程で分解した複数の周波数帯域のうち当該所望の周波数帯域に対応する周波数帯域に属する係数を用いて変更する変更工程と、
前記変更工程で変更した周波数帯域を含む全ての周波数帯域の係数に基づいて前記分解工程で行った逆の処理を行うことにより、画像を再構成する再構成工程と、
前記再構成工程で再構成された画像を、前記第1フレーム画像として出力する出力工程と
を備えることを特徴とする画像処理方法。
【請求項15】
コンピュータに請求項14に記載の画像処理方法を実行させることを特徴とするプログラム。
【請求項16】
請求項15に記載のプログラムを格納したことを特徴とする、コンピュータ読み取り可能な記憶媒体。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【公開番号】特開2007−94555(P2007−94555A)
【公開日】平成19年4月12日(2007.4.12)
【国際特許分類】
【出願番号】特願2005−280372(P2005−280372)
【出願日】平成17年9月27日(2005.9.27)
【出願人】(000001007)キヤノン株式会社 (59,756)
【Fターム(参考)】
【公開日】平成19年4月12日(2007.4.12)
【国際特許分類】
【出願日】平成17年9月27日(2005.9.27)
【出願人】(000001007)キヤノン株式会社 (59,756)
【Fターム(参考)】
[ Back to top ]