JP2005192154A - 画像処理装置、画像処理方法、プログラム及び記憶媒体 - Google Patents

Abstract

【課題】 階調処理の前後で画像処理効果が変動するという課題があった。
【解決手段】 調整回路１１４が画像の画素値幅の幅を調整し、階調変換することで画素値幅に依存せず画像全体の濃度が安定する。また、ＬＵＴ補正回路１１８が画像の画素値幅の伸縮率に応じて周波数数成分変換変換テーブルを伸縮させることにより、画像の画素値幅伸縮の影響によるコントラスト及び周波数処理効果のばらつきを抑えることができる。
【選択図】 図１

Description

本発明は、画像に画像処理を施すための画像処理装置、画像処理方法、プログラム及び記憶媒体に関し、特に画像処理効果を調整する画像処理装置、画像処理方法、プログラム及び記憶媒体に関する。

近年のデジタル技術の進歩により放射線画像等をデジタル画像信号に変換し、該デジタル画像信号に対して周波数処理などの画像処理を施し、ＣＲＴ等に表示出力したり、あるいはプリンタによってフィルム上に出力することが行われている。ところで、このような周波数処理においては、対象画像を複数の周波数帯の画像成分に分解し、該分解された成分値を変換テーブルを用いて増減することで周波数処理が行われる。

また、該周波数処理されたデジタル画像信号は、ＣＲＴ・フィルム等の出力媒体のダイナミックレンジに最適に収まるように階調処理が行われる場合がある（例えば、非特許文献１参照）。例えば、画像中の被写体の画素値幅が出力媒体のダイナミックレンジに収まるように変換テーブル（階調変換曲線とも呼ぶ）を調整することによって行われる。この方法によれば、出力媒体に与えられたダイナミックレンジ内に被写体画像全体が収まるように処理されるため、出力媒体のダイナミックレンジを最大限に活用でき、診断し易い出力画像が得られる。また、画素値幅の異なる被写体においても出力媒体のダイナミックレンジ内に収まるように調整されるため同等の濃度分布を持つ出力画像が得られる。
コンピュータ画像処理入門、ＰＰ９５−９７、田村秀行監修、総研出版発行

しかしながら、上述の階調処理は被写体の画素値幅に依存して変換テーブルを変更するため階調処理の前後で画像のコントラストが変化するという問題がある。具体的には、被写体の画素値幅が狭い場合においては階調処理後の画像の当該画素値幅が広がるためコントラストが高くなり、逆に被写体の画素値幅が広い場合においては階調処理後の画像の当該画素値幅が狭まるためコントラストが不足するという問題が生じる。

本発明は上記のような課題を解決するためになされたもので、階調処理の前後でも画像処理効果が変動しない、安定した画像処理効果の得られる画像処理装置、画像処理方法、プログラム及び記憶媒体を提供することを目的とする。

本発明の目的を達成するために、例えば本発明のＸ線撮影装置は以下の構成を備える。すなわち、少なくとも２以上の周波数帯域の周波数成分に画像を分解する分解手段と、前記周波数成分の値を変更するための基本変換テーブルを作成する基本ＬＵＴ作成手段と、前記画像に基づいて定まる画素値幅から前記基本変換テーブルの伸縮率を設定する伸縮率設定手段と、前記基本変換テーブルを前記設定された伸縮率で変更するＬＵＴ補正手段とを備える。

本発明の目的を達成するために、例えば本発明のＸ線撮影方法は以下の構成を備える。すなわち、少なくとも２以上の周波数帯域の周波数成分に画像を分解する分解工程と、前記周波数成分の値を変更するための基本変換テーブルを作成する基本変換テーブル作成工程と、前記画像に基づいて定まる画素値幅から前記基本変換テーブルの伸縮率を設定する伸縮率設定工程と、前記基本変換テーブルを前記設定された伸縮率で変更するＬＵＴ補正工程とを備える。

本発明の目的を達成するために、例えば本発明のプログラムは以下の構成を備える。すなわち、画像処理方法をコンピュータに実行させるためのプログラムであって、少なくとも２以上の周波数帯域の周波数成分に画像を分解する分解工程と、前記周波数成分の値を変更するための基本変換テーブルを作成する基本ＬＵＴ作成工程と、前記画像に基づいて定まる画素値幅から前記基本変換テーブルの伸縮率を設定する伸縮率設定工程と、前記基本変換テーブルを前記設定された伸縮率で変更するＬＵＴ補正工程を備える。

本発明の目的を達成するために、例えば本発明のコンピュータ可読記憶媒体は以下の構成を備える。すなわち、画像処理方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記憶媒体であって、少なくとも２以上の周波数帯域の周波数成分に画像を分解する分解工程と、前記周波数成分の値を変更するための基本変換テーブルを作成する基本ＬＵＴ作成作成工程と、前記画像に基づいて定まる画素値幅から前記基本変換テーブルの伸縮率を設定する伸縮率設定工程と、前記基本変換テーブルを前記設定された伸縮率で変更するＬＵＴ補正工程を備える。

階調処理の前後でも画像処理効果が変動しない、安定した画像処理効果の得られる画像処理装置、画像処理方法、プログラム及び記憶媒体を提供する。

以下、添付の図面を参照して本発明の好適な実施形態について説明する。

実施の形態１．
図１は、本実施形態のＸ線撮影装置１００の構成を示すブロック図である。図１に示すように、Ｘ線撮影装置１００は、撮影された画像をフィルム上又はモニタ上に出力する際の画像処理を行う機能を有するＸ線の撮影装置である。Ｘ線撮影装置１００は、データ収集回路１０５、前処理回路１０６、ＣＰＵ１０８、メインメモリ１０９、操作パネル１１０、画像表示器１１１、画像処理回路１１２を備えており、これらはＣＰＵバス１０７を介して互いにデータ授受が可能に接続されている。

また、Ｘ線撮影装置１００において、データ収集回路１０５と前処理回路１０６は相互に接続されており、データ収集回路１０５には２次元Ｘ線センサ１０４及びＸ線発生回路１０１とが接続されている。

上述の様なＸ線撮影装置１００において、
まず、メインメモリ１０９は、ＣＰＵ（制御手段）１０８での処理に必要な各種のデータ、プログラムなどが記憶されるものであると共に、ＣＰＵ（制御手段）１０８の作業用としてのワークメモリを含む。

ＣＰＵ（制御手段）１０８は、メインメモリ１０９を用いて、操作パネル１１０からの操作にしたがった装置全体の動作制御等を行う。これによりＸ線撮影装置１００は、以下のように動作する。

まず、操作パネル１１０を介してユーザから撮影指示が入力されると、この撮影指示はＣＰＵ１０８によりデータ収集回路１０５に伝えられる。ＣＰＵ（制御手段）１０８は、撮影指示を受けると、Ｘ線発生回路１０１及び２次元Ｘ線センサ１０４を制御してＸ線撮影を実行させる。

Ｘ線撮影では、まずＸ線発生回路１０１が、被写体１０３に対してＸ線ビーム１０２を放射する。Ｘ線発生回路１０１から放射されたＸ線ビーム１０２は、被写体１０３を減衰しながら透過して、２次元Ｘ線センサ１０４に到達する。そして、２次元Ｘ線センサ１０４によりＸ線画像信号が出力される。本実施形態では、被写体１０３を人体とする。すなわち、２次元Ｘ線センサ１０４から出力されるＸ線画像は人体画像となる。

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

１１２は画像処理回路の構成を示すブロック図であり、１１２において、１１３は入力画像から被写体の画素値幅を解析する解析回路、１１４は解析回路１１３で解析された画素値幅に基づいて画像の画素値幅を調整する調整回路、１１５は調整回路１１４で調整された画像を周波数分解する分解回路、１１６は分解回路１１５で分解された各周波数成分の値を変更するための基本となる変換テーブルを作成する基本ＬＵＴ作成回路、１１７は解析回路１１３で解析された画素値幅に基づいて変換テーブルの伸縮率を設定する伸縮率設定回路、１１８は伸縮率設定回路で設定された伸縮率に基づいて基本ＬＵＴ作成回路１１６で作成された変換テーブルを伸縮するＬＵＴ補正回路、１１９はＬＵＴ補正回路１１８で補正された変換テーブルによって分解回路１１５で分解された各周波数成分の値を変換する変換回路、１２０は変換回路１１９で変換された各周波数成分の値から画像を復元する再構成回路、１２１は再構成回路１１９で再構成された画像の階調を最適化する階調処理回路を備える。

以上のような構成を備えた本実施形態のＸ線撮影装置の動作について、以下により詳細に説明する。図２は本実施形態によるＸ線装置１００の処理の流れを示すフローチャートである。

上述の如く前処理回路１０６によって得られた原画像はＣＰＵバス１０７を介して画像処理１１２に転送される。まずステップＳ２０１において、画像処理装置１１２の解析回路１１３が原画像の画素値幅を解析する。ここでは、出力媒体に表示したい領域の画素値幅を算出する。本実施形態では入力画像内の被写体の画素値幅を算出する。算出方法は特に限定するものでないが、例えば、入力画像全体の画素値から最大画素値の９０％の値をＴｈ１として算出し、次式により入力画像ｆ（ｘ，ｙ）から素抜け及び素抜け領域と一定間隔内で接する体領域を削除した後の画像ｆ１（ｘ，ｙ）を求める。ここで、ｄ１およびｄ２は入力画像ｆ（ｘ，ｙ）から素抜け領域と一定間隔内で接する体領域を削除する際の一定間隔の幅を決定するための定数を示し、入力画像ｆ（ｘ，ｙ）の大きさ等に基づき設定される。

以上のように素抜けを除去した画像ｆ１（ｘ，ｙ）から画素値０を除く、最大値及び最小値から被写体の画素値幅を算出する。この方法は本出願人によりすでに提案されているものである（特開２００１−０９４８２８参照）。

次に、ステップｓ２０２において調整回路１１４が解析回路１１３で算出された被写体の画素値幅に基づき画像の伸縮処理を行う。本実施形態では、図３のように異なる被写体の最大値及び最小値が同等の値となるように線形変換する。図３において、Ｍｉｎが被写体の最小値、Ｍａｘが被写体の最大値を示す。下部の図は横軸が画素値であり縦軸が出現頻度を示している。上側の図は調整回路１１４で用いる階調変換テーブル（階調変換曲線とも呼ぶ）であり、横軸が入力値で、縦軸が出力値を示している。図３において出力レンジとは出力媒体のダイナミックレンジを意味するものであり、例えばフィルムなら濃度、モニタなら輝度をその単位とする。

また、図４のように被写体内の画素値領域毎に階調変換テーブルの傾きを変えてもよい。

次に、ステップｓ２０３において分解回路１１５が調整回路１１４によって画素値幅調整された画像を複数の周波数帯の成分に分解する。本実施形態ではラプラシアンピラミッドアルゴリズムを用いて周波数分解を行う。

ラプラシアンピラミッドの分解について図５の（ａ）及び（ｂ）を参照して説明する。図５（ａ）のｇおよびｂはおのおの１レベルの分解を行った際の画像の低解像度近似画像および高周波成分を表しており、低解像度近似画像ｇは画像信号ｘをローパスフィルタによりフィルタ処理し、ダウンサンプリングして得られる。また高周波成分ｂは画像信号ｘと低解像度近似画像ｇをアップサンプリングしさらにローパスフィルタによりフィルタ処理した画像との差分をとることにより得られる。ローパスフィルタは例えば図５（ｂ）に示すようなフィルタを用いる。ラプラシアンピラミッドはこの処理を低解像度近似画像ｇに対して繰り返し行なうことで各周波数帯の画像成分を得るものであり、その詳細は公知であるのでここでは説明を省略する。画像成分を変更し、後述の復元処理を行うことで画像処理効果を得ることができる。鮮鋭化処理、ダイナミックレンジ調整処理、ノイズ低減処理などの処理効果を周波数成分を変更することで得ることができる。この様に画像処理効果を変更するために用いる成分を周波数成分と呼ぶこととする。

尚、本実施形態では画像の周波数成分分解をラプラシアンピラミッドアルゴリズムによって行うが、これに限るものではなく離散ウェーブレット変換他の方法を用いても実行可能である。さらに、いわゆる移動平均法を用いて高周波成分及び低周波成分を算出する方法を用いても良い。すなわち、本実施形態はラプラシアンピラミッドアルゴリズムを用いて説明されたが、その技術的思想は、他の周波数成分分解処理を用いた場合にも適用することができる。

また、ウエーブレット変換を用いた周波数分解では、周波数成分を周波数係数と呼ぶ場合もあるが、周波数係数が画像処理効果を変更するために用いられる場合には周波数成分と呼ぶこととする。

また同様に、移動平均を用いた周波数分解では、周波数成分を画像成分又は画素値と呼ぶ場合もあるが、同様に、移動平均を用いた周波数分解では、周波数成分を画像処理効果を変更するために用いられる場合には周波数成分と呼ぶこととする。

以上のように画像を周波数成分に分解し、各周波数帯の成分値を変換テーブルにより変更することで周波数処理を行う。以下に周波数処理方法について詳細に説明する。

まず、ステップｓ２０４において、基本ＬＵＴ作成回路１１６が分解回路１１５で分解された複数の周波数帯のうち所望とする周波数帯の成分値を変更するための基本変換テーブル（変換テーブルを変換曲線と呼ぶ場合もあるが以後「変換テーブル」と呼ぶ）とを作成する。なお、基本変換テーブルとは、周波数成分の値を変換するための変換テーブルのことである。ここで、作成する基本変換テーブルは目的とする画像処理によって決定されるもので、特に限定されるものでない。例えば、図６のように入力値を非線形に増加させるような変換テーブルでも良いし、図７のように入力値を非線形に減少させるような基本変換テーブルを用いても良い。尚、実際の入力値は正負の値をとるため、基本変換テーブルは本来奇関数であるが、図６及び図７はその第１象限のみを示している。ここで、一般に入力値を増加する場合は鮮鋭化、入力値を減少する場合には平滑化の効果が得られる。また、周波数処理を行わない場合においてはリニアの変換テーブルを用いれば良い。尚、ここで作成される基本変換テーブルは１つではなく、対象となる周波数帯毎に異なる基本変換テーブルを作成しても良い。

次に、ステップｓ２０５において、伸縮率設定回路１１７が調整回路１１４で調整された画素値幅に基づいて基本変換テーブルの伸縮率を算出するものである。本実施形態では、伸縮率は調整回路１１４で伸縮された画像の伸縮率と同じ値が設定される。具体的には、調整回路１１４によって画像のダイナミックが０．５倍に圧縮された場合では伸縮率として０．５が設定される。また、画像の画素値幅が１．５倍に伸長された場合は伸縮率として１．５が設定される。尚、図４のように被写体の領域毎に伸縮率が異なる場合においては、領域毎に対応する伸縮率が設定される。

次に、ステップｓ２０６において、ＬＵＴ補正回路１１８が伸縮率設定回路１１７で設定された伸縮率に応じて、基本ＬＵＴ作成回路１１６で作成された基本変換テーブルを伸縮する。本実施形態では、伸縮率に応じて基本変換テーブルを出力軸方向（水平方向）に伸縮する。具体的には図８のように伸縮率が１．５ならば入力軸のスケールを１．５倍にする。また、伸縮率が０．５ならば入力軸のスケールを０．５倍にすれば良い。こうすることで、画像伸縮に応じて基本変換テーブルのスケールも同様に補正されるので、変更後の基本変換テーブルの出力値は常に一定のレンジ幅となる。これにより、画像伸縮に依存しないで画像処理効果が常に一定となるものである。

次に、ステップｓ２０７において、変換回路１１９が周波数帯の成分値をＬＵＴ補正回路１１８で補正された変換テーブルを用いて変換することで周波数処理を行う。ここで、図４の場合のように被写体の領域毎に伸縮率が異なる場合においては、各領域に対応する周波数成分毎に別々の補正変換テーブルを適用すれば良い。

以上のように、周波数処理が行われるが、周波数処理に用いられる基本変換テーブルは画像の伸縮率に応じて補正される。よって、画像伸縮でコントラストが変動しても周波数処理された周波数帯域のコントラストは常に一定となり、安定した周波数効果が得られる。また、画像の伸縮によるコントラスト変動の補正だけを行う場合においては、基本変換テーブルをリニアに設定し、所望とする周波数帯に変換処理を施せば良い。

次にステップｓ２０８において、再構成回路１２０が変換回路１１９で変換された周波数成分を用いてラプラシアンピラミッドの再構成処理を行う。図５（ｃ）が再構成回路１２０の構成を示す図である。入力された高周波成分係数ｂと低解像度近似画像ｇをアップサンプリングしローバスフィルタによりフィルタ処理をしたものを重ね合わせることにより画像信号ｘ’が出力される。この処理を各レベルに対して繰り返し行なうことで画像を復元することができるものであり、その詳細は公知であるのでここでは説明を省略する。

次にステップｓ２０９において、階調処理回路１２１が出力媒体に応じた変換テーブルを用いて最適な階調処理が行われる。変換テーブルについては特に限定しないが、例えばＸ線フィルムに出力する際にはフィルムの特性曲線と同等のＳ字曲線を用いれば良い。尚、ここで階調処理される画像はあらかじめ画素値幅とコントラストを補正したものであり、同一の変換テーブルを用いて処理すれば同等の濃度分布とコントラストを持つ出力画像を得ることができる。

以上の様に本実施形態によれば、画像の画素値幅の幅を調整し、階調変換することで画素値幅に依存せず画像全体の濃度が安定する効果がある。また、画像の画素値幅の伸縮率に応じて周波数数成分変換変換テーブルを伸縮させることにより、画像の画素値幅伸縮の影響によるコントラスト及び周波数処理効果のばらつきを抑える効果がある。

また、本実施の形態では調整回路１１４で画像を伸縮した後に画像処理を行った例を示したが、復元処理後に調整回路１１４で画像を伸縮してもよい。この場合には、ステップｓ２０６において、ＬＵＴ補正回路１１８が伸縮率設定回路１１７で設定された伸縮率に応じて、基本ＬＵＴ作成回路１１６で作成された基本変換テーブルを伸縮する。伸縮率に応じて基本変換テーブルを入力軸方向（垂直方向）に伸縮する。これにより、画像伸縮に依存しないで画像処理効果が常に一定となるものである。

伸縮率を設定し基本変換テーブルを変更することで変換する周波数成分の振幅レベルに応じて柔軟な周波数処理を行える効果がある。

また、画像の画素値幅の幅を調整することで画像全体の濃度が安定する効果がある。さらに、画像の画素値幅の伸縮率に応じて基本変換テーブルを伸縮させることにより、画像の画素値幅伸縮の影響による画像処理効果のばらつきを抑える効果がある。

伸縮率を設定し基本変換テーブルを補正することで変換する周波数成分の振幅レベルに応じて柔軟な周波数処理を行える効果がある。また、基本変換テーブルを入力軸方向あるいは出力軸方向の少なくともどちらか一方向に対して伸縮することで補正するため、非線形変換テーブルの場合においても効率よく補正が行える効果がある。

尚、本発明の目的は、実施形態の装置又はシステムの機能を実現するソフトウェアのプログラムコードを記憶した記憶媒体を、装置又はシステムに供給し、その装置又はシステムのコンピュータ（ＣＰＵ又はＭＰＵ等）が記憶媒体に格納されたプログラムコードを読みだして実行することによっても、達成されることは言うまでもない。

この場合、記憶媒体から読み出されたプログラムコード自体が実施形態の機能を実現することとなり、そのプログラムコードを記憶した記憶媒体及び当該プログラムコードは本発明を構成することとなる。

プログラムコードを供給するための記憶媒体としては、ＲＯＭ、フロッピー（登録商標）ディスク、ハードディスク、光ディスク、光磁気ディスク、ＣＤ−ＲＯＭ、ＣＤ−Ｒ、磁気テープ、不揮発性のメモリカード等を用いることができる。

また、コンピュータが読み出したプログラムコードを実行することにより、実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼動しているＯＳ等が実際の処理の一部又は全部を行い、その処理によって実施形態１〜２の機能が実現される場合も本発明の実施の態様に含まれることは言うまでもない。

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

このようなプログラム又は当該プログラムを格納した記憶媒体に本発明が適用される場合、当該プログラムは、例えば、上述の図２に示されるフローチャートに対応したプログラムコードから構成される。

この発明の実施の形態１による画像処理装置のブロック図である。 この発明の実施の形態１による画像処理装置の処理手順を示すフローチャートである。 画素値幅伸縮処理の一例を示した図である。 画素値幅伸縮処理の一例を示した図である。 ラプラシアンピラミッドの分解及び再構成を説明する図である。 基本変換テーブルの一例を示した図である。 基本変換テーブルの一例を示した図である。 基本変換テーブルの伸縮処理を説明する図である。

Claims (11)

1. 少なくとも２以上の周波数帯域の周波数成分に画像を分解する分解手段と、
   前記周波数成分の値を変更するための基本変換テーブルを作成する基本ＬＵＴ作成手段と、
   前記画像に基づいて定まる画素値幅から前記基本変換テーブルの伸縮率を設定する伸縮率設定手段と、
   前記基本変換テーブルを前記設定された伸縮率で変更するＬＵＴ補正手段とを備えることを特徴とする画像処理装置。
2. 前記画像中から被写体の画素値幅を解析する解析手段を更に有し、前記伸縮率設定手段は、該解析手段で解析された画素値幅に基づいて前記基本変換テーブルの伸縮率を設定することを特徴とする請求項１に記載の画像処理装置。
3. 前記ＬＵＴ補正手段は、前記伸縮率設定手段で設定された伸縮率に基づいて基本ＬＵＴ作成手段によって作成された基本変換テーブルを入力軸方向あるいは出力軸方向の少なくともどちらか一方向に対して伸縮することを特徴とする請求項１又は２に記載の画像処理装置。
4. 前記ＬＵＴ補正手段で変更された変換テーブルを用いて周波数帯域の成分を変換する変換手段を備えることを特徴とする請求項１乃至３に記載の画像処理装置。
5. 前記変換手段によって変換された周波数帯域の周波数成分を用いて画像を再構成する再構成手段を備えることを特徴とする請求項４に記載の画像処理装置。
6. 画像の画素値幅を調整する調整手段を更に有し、前記画像は該調整手段で画素値幅を調整した後の画像であることを特徴とする請求項１乃至５に記載の画像処理装置。
7. 前記再構成手段によって再構成された画像の画素値幅を調整する第２の調整手段を備えることを特徴とする請求項５に記載の画像処理装置。
8. 前記分解手段はラプラシアンピラミッドアルゴリズムを用いていることを特徴とする請求項１乃至７に記載の画像処理装置。
9. 少なくとも２以上の周波数帯域の周波数成分に画像を分解する分解工程と、
   前記周波数成分の値を変更するための基本変換テーブルを作成する基本変換テーブル作成工程と、
   前記画像に基づいて定まる画素値幅から前記基本変換テーブルの伸縮率を設定する伸縮率設定工程と、
   前記基本変換テーブルを前記設定された伸縮率で変更するＬＵＴ補正工程とを備えることを特徴とする画像処理方法。
10. 画像処理方法をコンピュータに実行させるためのプログラムであって、
    少なくとも２以上の周波数帯域の周波数成分に画像を分解する分解工程と、
    前記周波数成分の値を変更するための基本変換テーブルを作成する基本ＬＵＴ作成工程と、
    前記画像に基づいて定まる画素値幅から前記基本変換テーブルの伸縮率を設定する伸縮率設定工程と、
    前記基本変換テーブルを前記設定された伸縮率で変更するＬＵＴ補正工程をコンピュータに実行させるためのプログラム。
11. 画像処理方法をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記憶媒体であって、少なくとも２以上の周波数帯域の周波数成分に画像を分解する分解工程と、
    前記周波数成分の値を変更するための基本変換テーブルを作成する基本ＬＵＴ作成作成工程と、
    前記画像に基づいて定まる画素値幅から前記基本変換テーブルの伸縮率を設定する伸縮率設定工程と、
    前記基本変換テーブルを前記設定された伸縮率で変更するＬＵＴ補正工程
    をコンピュータに実行させるためのプログラムを記録したコンピュータ読み取り可能な記憶媒体。

Images