JP2004171301A - 画像処理方法及び装置 - Google Patents

Abstract

【課題】対象画像のコントラストを効率的且つ効果的に改善する。
【解決手段】Ｘ線撮影装置１００における画像処理回路１１２は、対象画像のコントラストを改善するための処理を行なう。まず周波数成分分解回路１１３が画像を複数の周波数帯域の成分値に分解し、周波数帯域毎の係数を得る。係数変換回路１１４は、複数の周波数帯域の少なくとも一つの周波数帯域について、周波数帯域毎の係数値の頻度分布が所定の頻度分布となるように係数値を変換する。そして、復元回路１１５は、係数変換回路１１４ににより変換された係数値を用いて画像を生成する。
【選択図】 図１

Description

【０００１】
【発明の属する技術分野】
本発明は、画像処理装置及び方法に関する。特に、対象画像のコントラストを改善する画像処理装置及び方法に関する。
【０００２】
【従来の技術】
近年のデジタル技術の進歩により放射線画像等をデジタル画像信号に変換し、該デジタル画像信号に対して周波数処理などの画像処理を施し、ＣＲＴ等に表示出力したり、あるいはプリンタによってフィルム上に出力することが行われている。ところで、このような周波数処理においては、対象画像を複数の周波数帯の画像係数に分解し、周波数帯毎の画像係数値を個別的に増減することにより画像処理が行われる。
【０００３】
一方、濃度ヒストグラム一様化の方法（例えば、非特許文献１参照）を用いて、画像の濃度（画素値）のヒストグラム（頻度分布）を一様化することが、一般に知られている。この方法を用いると、与えられた濃度レンジを効率よく使用できるため、画像のコントラストが改善される。
【０００４】
【非特許文献１】
画像理解のためのディジタル画像処理〔１〕、ＰＰ１３３−１３５、鳥脇純一郎著、昭晃堂発行、初版第４刷。
【０００５】
【発明が解決しようとする課題】
上述の方法は、周波数係数の値を変更することで、所望の周波数処理効果を得ようとするものであるが、単に係数の値の強弱を変更しているだけであり、係数変換後の係数値のヒストグラムを制御する思想はない。すなわち、従来の係数値の増減とは、周波数帯毎に個別に係数値の増減を行うのみで、各周波数帯の係数値ヒストグラムを考慮しておらず、所望の効果を得るために試行錯誤的に係数を調整する必要があった。また、単純に周波数係数の値を変更すると、処理後の画像のダイナミックレンジに影響を与えたり、画像中のエッジ（輪郭）部分にオーバーシュート等のアーティファクトが生じたりする場合がある。
【０００６】
また、非特許文献１で示されているように画素値ヒストグラムを一様化することにより、画像のコントラストを改善することができるが、この方法は階調変換処理を行うものであって周波数処理ではなく、よって周波数帯毎の成分を調整することはできない。
【０００７】
本発明者らは周波数帯毎の係数値の頻度分布を目的に応じて制御することにより、所望のコントラストの画像が得られることを見出した。しかしながら、上述した各従来技術においては、係数変換後に得られる係数値の頻度分布を考慮した係数変換は行われていない。
【０００８】
本発明は、対象画像のコントラストを効率的且つ効果的に改善できる画像処理装置及び方法を提供することを目的とする。
【０００９】
【課題を解決するための手段】
上記の目的を達成するための本発明による画像処理装置は以下の構成を備える。すなわち、
画像を複数の周波数帯域の成分に分解する分解手段と、
前記複数の周波数帯域の少なくとも一つの周波数帯域について、周波数帯域毎の成分値の頻度分布が所定の頻度分布となるように成分値を変換する変換手段と、
前記変換手段により得られた成分値を用いて画像を生成する生成手段とを備える。
【００１０】
また、上記の目的を達成するための本発明による画像処理方法は、
画像を複数の周波数帯域の成分に分解する分解工程と、
前記複数の周波数帯域の少なくとも一つの周波数帯域について、周波数帯域毎の成分値の頻度分布が所定の頻度分布となるように成分値を変換する変換工程と、
前記変換工程により得られた成分値を用いて画像を生成する生成工程とを備える。
【００１１】
【発明の実施の形態】
以下、添付の図面を参照して本発明の好適な実施形態について説明する。
【００１２】
本実施形態では、周波数処理を行う画像処理装置であって、周波数帯毎の係数値の頻度分布が所定の態様となるように係数値を変換する画像処理装置が開示される。
【００１３】
図１は、本実施形態のＸ線撮影装置１００の構成を示すブロック図である。図１に示すように、Ｘ線撮影装置１００は、撮影された画像をフイルム上又はモニタ上に出力する際の効果的な画像処理を行う機能を有するＸ線の撮影装置である。Ｘ線撮影装置１００は、データ収集回路１０５、前処理回路１０６、ＣＰＵ１０８、メインメモリ１０９、操作パネル１１０、画像表示器１１１、画像処理回路１１２を備えており、これらはＣＰＵバス１０７を介して互いにデータ授受が可能に接続されている。
【００１４】
また、Ｘ線撮影装置１００において、データ収集回路１０５と前処理回路１０６は相互に接続されており、データ収集回路１０５には２次元Ｘ線センサ１０４及びＸ線発生回路１０１とが接続されている。更に、画像処理回路１１２は周波数成分分解回路１１３、係数変換回路１１４、復元回路１１５、変更回路１１６を含んで構成されており、各回路はＣＰＵバス１０７に接続されている。
【００１５】
なお、上述したＸ線撮影装置１００の各部の機能は、後述の説明から明らかにされる。
【００１６】
さて、上述の様なＸ線撮影装置１００において、まず、メインメモリ１０９は、ＣＰＵ１０８での処理に必要な各種のデータなどを記憶すると共に、ＣＰＵ１０８のワーキング・メモリとして機能する。ＣＰＵ１０８は、メインメモリ１０９を用いて、操作パネル１１０からの操作にしたがった装置全体の動作制御等を行う。これによりＸ線撮影装置１００は、以下のように動作する。
【００１７】
まず、操作パネル１１０を介してユーザから撮影指示が入力されると、この撮影指示はＣＰＵ１０８によりデータ収集回路１０５に伝えられる。データ収集回路１０８は、撮影指示を受けると、Ｘ線発生回路１０１及び２次元Ｘ線センサ１０４を制御してＸ線撮影を実行させる。
【００１８】
Ｘ線撮影では、先ず、Ｘ線発生回路１０１が、被検査体１０３に対してＸ線ビーム１０２を放射する。Ｘ線発生回路１０１から放射されたＸ線ビーム１０２は、被検査体１０３を減衰しながら透過して、２次元Ｘ線センサ１０４に到達する。そして、２次元Ｘ線センサ１０４によりＸ線画像信号が出力される。本実施形態では、被写体１０３を人体とする。すなわち、２次元Ｘ線センサ１０４から出力されるＸ線画像は人体画像となる。
【００１９】
データ収集回路１０５は、２次元Ｘ線センサ１０４から出力されたＸ線画像信号を所定のディジタル信号に変換してＸ線画像データとして前処理回路１０６に供給する。前処理回路１０６は、データ収集回路１０５からの信号（Ｘ線画像データ）に対して、オフセット補正処理やゲイン補正処理等の前処理を行う。この前処理回路１０６で前処理が行われたＸ線画像データは原画像データとして、ＣＰＵ１０８の制御により、ＣＰＵバス１０７を介して、メインメモリ１０９、画像処理回路１１２に転送される。
【００２０】
画像処理回路１１２では、周波数成分分解回路１１３が原画像に対して離散ウェーブレット変換（ＤＷＴ）を施し、各周波数帯の係数（ウェーブレット変換係数）を得る。係数変換回路１１４は、周波数成分分解回路１１３で得られた各周波数帯毎の係数値の頻度分布に基づいて、当該頻度分布を所定形状の頻度分布とするための係数変換テーブルを生成し、この係数変換テーブルを用いて、周波成分分解回路１１３で算出された周波数係数を変換する。最後に、復元回路１１５が、係数変換回路１１４によって得られた係数に対して逆離散ウェーブレット変換（逆ＤＷＴ）を施し、周波数処理後の画像を得る。また、変更回路１１６は操作パネル１１０からの指示等に基づいて係数変換回路１１４の変換特性（変換テーブル）を変更する。
【００２１】
以上のような構成を備えた本実施形態のＸ線撮影装置の動作について、以下により詳細に説明する。図２は本実施形態によるＸ線撮影装置１００の処理の流れを示すフローチャートである。
【００２２】
上述の如く前処理回路１０６によって得られた原画像はＣＰＵバス１０７を介して画像処理装置１１２に転送される。まずステップＳ２０１において、画像処理装置１１２の周波数成分分解回路１１３が原画像ｆ（ｘ、ｙ）に対して２次元の離散ウェーブレット変換処理を行い、周波数帯毎の係数を計算して出力する。
【００２３】
離散ウェーブレット変換について図６の（ａ）及び（ｂ）を参照して説明する。図６（ａ）は周波数成分分解回路１１４の構成を示す図であり、図６（ｂ）は２次元のＤＷＴ処理により得られる２レベルの変換係数群の構成例を示す図である。
【００２４】
入力された画像信号は遅延素子およびダウンサンプラの組み合わせにより、偶数アドレスおよび奇数アドレスの信号に分離され、２つのフィルタｐおよびｕによりフィルタ処理が施される。図６（ａ）のｓおよびｄは、各々１次元の画像信号に対して１レベルの分解を行った際のローパス係数およびハイパス係数を表しており、次式により計算されるものとする。
【００２５】
ｄ_ｎ＝ｘ_２ｎ＋１ − ｆｌｏｏｒ（（ｘ_２ｎ＋ｘ_２ｎ＋２）／２） （３）
ｓ_ｎ＝ｘ_２ｎ ＋ ｆｌｏｏｒ（（ｄ_ｎ−１ ＋ ｄ_ｎ））／４） （４）
ただし、ｘ_ｎは変換対象となる原画像の画素値である。
【００２６】
以上の処理により、画像データに対する１次元の離散ウェーブレット変換処理が行われる。２次元の離散ウェーブレット変換は、上述した１次元の変換を画像の水平・垂直方向に対して順次行うものであり、その詳細は公知であるのでここでは説明を省略する。
【００２７】
図６（ｂ）は２次元の変換処理により得られる２レベルの変換係数群の構成例であり、画像信号は異なる複数の周波数帯域の係数群ＨＨ１，ＨＬ１，ＬＨ１，…，ＬＬに分解される。ここで、ＨＨ１，ＨＬ１，ＬＨ１，…，ＬＬ等は、それぞれサブバンドとも呼ばれ、それぞれ周波数帯毎の係数群を示す。
【００２８】
次に、ステップＳ２０２において、ステップＳ２０１で得られたサブバンドのうち高周波係数を示すＨＬ，ＬＨ，ＨＬの各サブバンドに対して、係数値のヒストグラムを平坦化（一様化）するための解析を行なう。そして、ステップＳ２０３において、その解析結果に基づいて係数変換テーブルを作成する。なお、本実施形態では、ＨＨ１，ＨＬ１，ＬＨ１，ＨＨ２，ＨＬ２，ＬＨ２の６つのサブバンドに対して係数の変換を行なうが、ＨＨ１，ＨＬ１，ＬＨ１に対してのみ係数の変換を行なうなど、目的に合った全部又は一部のサブバンド（周波数帯）の係数に対して変換処理を行えばよい。
【００２９】
例えば、係数値ヒストグラムを一様化する場合には、サブバンド毎に係数値の累積ヒストグラムを作成し、その累積ヒストグラムの形状を正規化して係数変換テーブル（係数変換関数）とすればよい。即ち、横軸に係数の値をとり、縦軸に係数値の累積頻度をとって累積ヒストグラム（累積頻度分布）を作成し、当該累積ヒストグラムの横軸を入力係数値、縦軸を出力係数値に置き換えて係数変換テーブルとすればよい。このようにして作成した係数変換テーブルの一例が図３の係数変換テーブル３０１である。図３は、周波数成分分解手段で算出された周波数帯毎の係数値を変換する係数変換テーブルの一例を示し、横軸が入力係数、縦軸が出力係数を示す。この図において、３０１が係数変換テーブルに相当する。尚、実際の係数値は正負の値をとるため、変換関数（変換テーブル）は本来奇関数であるが、図３ではその第１象限のみを示している。また、累積頻度分布をそのまま用いずに、それをスムージング（平滑化）したものに基づいて係数変換テーブルを作成してもよい。この場合、当該スムージングは累積頻度分布を多項式で近似する方法等により実現できる。
【００３０】
尚、変換後の係数値が所与の最大値を超えない範囲で、図３の変換テーブルを縦軸のスケールを伸縮させることができる。このようにすれば、係数値ヒストグラムの平坦性を維持しつつ、係数値の絶対値を調整することができる。縦軸のスケールを小さくすれば変換後の係数の値が大きくなり、逆に縦軸のスケールを大きくすると変換後の係数の値は小さくなる。これにより周波数処理効果を調整することが可能となる。一般に、係数を増加する場合には鮮鋭化の効果が得られ、係数を低減する場合には平滑化の効果が得られる。
【００３１】
ステップＳ２０５では、以上のようにして生成された変換テーブル３０１を用いて、ステップＳ２０１で生成された係数を変換する。係数変換テーブル３０１を用いて、サブバンドの係数値を変換すると、変換後の係数値のヒストグラムは、理想的には、図４の４０１に示すようなものとなる。図４は、係数変換テーブル３０１を用いて係数変換した場合の、変換後の係数値の頻度分布の１例を示す図である。図において、横軸が係数の値、縦軸が係数値の出現頻度を示す。ここで、前述の縦軸のスケールを変更すると、線４０１の高さ（頻度）及び幅（係数値範囲）が変化する。尚、このスケール変更のような係数変換特性（変換テーブル）の変更は、上述のように、操作パネル１１０からの指示等に基づいて、変更回路１１６により実行可能である。
【００３２】
このように、本実施形態の係数変換処理は、変換後の係数値の頻度分布が所定の形状（例えば平坦）になるように係数変換を行う。特に、係数値の頻度分布を図４に示すように一様にして、後述する復元処理を行うと、コントラストの強調された画像を得られることをが実験的に確認された。つまり、このような係数値ヒストグラム平坦化によれば、画像の周波数処理を効果的かつ効率的に行うことができる。
【００３３】
ステップＳ２０６では、ステップＳ２０２〜Ｓ２０５において係数変換回路１１４によって変更された後の全係数を用いて、復元回路１１５が逆離散ウェーブレット変換処理を行う。
【００３４】
図６（ｃ）は復元回路１１５の構成を示す図である。入力された係数はｕおよびｐの２つのフィルタで処理され、アップサンプリングされた後に重ね合わされて、出力画像信号ｘ’が生成される。これらの処理は次式、
ｘ’_２ｎ＝ｓ’_ｎ − ｆｌｏｏｒ（（ｄ’_ｎ−１＋ ｄ’_ｎ）／４） （５）
ｘ’_２ｎ＋１＝ｄ’_ｎ ＋ ｆｌｏｏｒ（（ｘ’_２ｎ＋ ｘ’_２ｎ＋２）／２） （６）
により表わされる。
【００３５】
以上の処理により、係数に対する１次元の逆離散ウェーブレット変換処理が行われる。２次元の逆離散ウェーブレット変換は、１次元の逆変換を画像の水平・垂直方向に対して順次行うものであり、その詳細は公知であるのでここでは説明を省略する。
【００３６】
なお、上記では各周波数帯の変換後の係数値の頻度分布が平坦化されるように制御したが、変換後の係数値の頻度分布の形状はこれに限られるものではない。例えば、当該頻度分布の形状として、中央部分が凸の山形の形状や、右上がり或いは右下がりの直線形状、理想的な実際の画像の係数値頻度分布に基づく形状等を適宜採用することができる。尚、このような頻度分布形状は、処理の対象や目的等に応じて実験的に決定することができる。
【００３７】
更に、頻度分布制御の変形例として、図５に示すような変換テーブルを用いて、所定の部分的な係数値範囲の係数値ヒストグラムだけを選択的に平坦化することもできる。図５は、図３と同様に係数変換テーブルの一例を示す図である。図５において、５０１はヒストグラムが平坦化される係数値範囲を示し、５０２は係数変換テーブルを示し、係数値範囲５０１に対応する部分はヒストグラムを平坦化する形状をなし、他の係数値範囲に対応する部分は係数を変換しないよう傾きが１の直線をなしている。一般に、高周波のサブバンドにおいて、係数値範囲５０１は連続的な構造物（人体における解剖学的構造物等）における微小構造などの有効情報に対応する範囲を示し、それより大きい係数値範囲は構造物の輪郭等のエッジ部分に対応する範囲を示す。換言すれば、係数値範囲５０１にのみ、ヒストグラム一様化のための係数変換を作用させる。
【００３８】
従って、図５の係数変換テーブルを用いた場合、復元画像において、係数値範囲５０１対応する画像領域ではコントラストが好適に強調され、かつ、他の係数値範囲に対応する画像領域としてのエッジ部では係数が不変に保たれているので、エッジ構造が保存される（擬輪郭が生じない）という効果がある。すなわち、いわゆるオーバーシュート等のアーティファクトを抑えつつ、効果的且つ効率的な周波数処理が行える。また、大きな係数値を不変に保つため、復元画像のダイナミックレンジの変動が抑制される。更に、このような係数変換テーブルを用いる場合、上述のように縦軸のスケールを変更することにより、係数値範囲５０１に対応する画像領域に対するコントラスト強調効果を調整することができる。また、係数値範囲５０１を対象画像、対象サブバンド、処理目的等に応じて変更することも有効である。
【００３９】
以上の様に本実施形態によれば、変換後の係数値の頻度分布が所定形状をなすように係数変換を行うことにより、目的とするコントラスト強調効果を効率的に得ることが出来る。また、ヒストグラム平坦化方法に従った係数変換テーブルを用いて係数変換し、復元処理した場合、一般的に良好なコントラスト強調効果を得ることが出来る。
【００４０】
さらに、絶対値が所定値以下の係数値範囲に対してはヒストグラム平坦化の方法による変換特性とし、絶対値が所定値を超える係数値範囲に対しては係数値を不変に保つようにした係数変換テーブルを用いることにより、処理後の画像中のエッジ部分に擬輪郭等のアーティフアクトを生じさせずに、微細構造などに対応する有効部分を効果的且つ効率的にコントラスト強調することができる。さらに、大きな係数値を不変に保つことで、復元画像のダイナミックレンジの変動を抑制できる効果がある。ダイナミックレンジの変動は大きな係数値の変更により生じる場合が多いからである。また、係数変換テーブルの縦軸スケールの変更により、コントラスト強調効果を調節することが出来る。
【００４１】
なお、画像の周波数成分分解は、離散ウェーブレット変換に限られず、ラプラシアンピラミッド変換等、他の方法を用いても実行可能である。さらに、いわゆる移動平均法を用いて高周波成分及び低周波成分を算出する方法を用いてもよい。この場合には、高周波画像をヒストグラム一様化の方法を用いて変換することになる。高周波成分が周波数成分分解回路１１３で作成され、その高周波成分が係数変換回路１１４によってヒストグラム一様化の方法で作成された変換テーブルを用いて変換され、係数変換回路１１４で変換された高周波成分が復元回路１１５で低周波成分に加算されることにより、処理後の画像が得られる。即ち、本実施形態は離散ウェーブレット変換を用いて説明されたが、その技術的思想は、他の周波数成分分解処理を用いた場合にも適用することができる。
【００４２】
なお、本発明の目的は、前述した実施形態の機能を実現するソフトウェアのプログラムコードを記録した記憶媒体を、システムあるいは装置に供給し、そのシステムあるいは装置のコンピュータ（またはＣＰＵやＭＰＵ）が記憶媒体に格納されたプログラムコードを読出し実行することによっても達成されることは言うまでもない。
【００４３】
この場合、記憶媒体から読出されたプログラムコード自体が前述した実施形態の機能を実現することになり、そのプログラムコードを記憶した記憶媒体は本発明を構成することになる。
【００４４】
プログラムコードを供給するための記憶媒体としては、例えば、フレキシブルディスク，ハードディスク，光ディスク，光磁気ディスク，ＣＤ−ＲＯＭ，ＣＤ−Ｒ，磁気テープ，不揮発性のメモリカード，ＲＯＭなどを用いることができる。
【００４５】
また、コンピュータがプログラムコードを読出し実行することにより、前述した実施形態の機能が実現されるだけでなく、そのプログラムコードの指示に基づき、コンピュータ上で稼働しているＯＳ（オペレーティングシステム）などが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も本発明の実施形態に含まれることは言うまでもない。
【００４６】
さらに、記憶媒体から読出されたプログラムコードが、コンピュータに挿入された機能拡張ボードやコンピュータに接続された機能拡張ユニットに備わるメモリに書込まれた後、そのプログラムコードの指示に基づき、その機能拡張ボードや機能拡張ユニットに備わるＣＰＵなどが実際の処理の一部または全部を行い、その処理によって前述した実施形態の機能が実現される場合も本発明の実施形態に含まれることは言うまでもない。
【００４７】
【発明の効果】
以上説明したように、本発明によれば、対象画像のコントラストを効率的且つ効果的に改善することができる。
【図面の簡単な説明】
【図１】本実施形態によるＸ線撮影装置の構成を示すブロック図である。
【図２】本実施形態のＸ線撮影装置における画像処理の手順を示すフローチャートである。
【図３】ヒストグラム一様化方法を用いた係数変換テーブルの一例を示す図である。
【図４】ヒストグラム一様化方法を用いた係数変換テーブルで係数変換された係数の理想的な例を示す図である。
【図５】ヒストグラム一様化方法を用いた係数変換テーブルと直線状形状の係数変換テーブルの複合形による変換テーブルを示す図である。
【図６】離散ウェーブレット変換を用いた係数分解回路と復元回路を示す図である。

Claims (18)

1. 画像を複数の周波数帯域の成分に分解する分解手段と、
   前記複数の周波数帯域の少なくとも一つの周波数帯域について、周波数帯域毎の成分値の頻度分布が所定の頻度分布となるように成分値を変換する変換手段と、
   前記変換手段により得られた成分値を用いて画像を生成する生成手段とを備えることを特徴とする画像処理装置。
2. 前記変換手段は、周波数帯毎の成分値の頻度分布が実質的に平坦となるように成分値を変換することを特徴とする請求項１に記載の画像処理装置。
3. 前記変換手段は、前記成分値の累積頻度分布に基づいて前記成分値を変換することを特徴とする請求項２に記載の画像処理装置。
4. 前記変換手段による前記成分値の変換特性を変更する変更手段を更に備えることを特徴とする請求項１乃至３のいずれかに記載の画像処理装置。
5. 前記変更手段は、変換後の成分値が一律に所定の割合で変更されるように、前記変換特性を変更することを特徴とする請求項４に記載の画像処理装置。
6. 前記分解手段はウェーブレット変換を用いて画像を複数の周波数帯域の成分に分解することを特徴とする請求項１乃至５のいずれかに記載の画像処理装置。
7. 前記変換手段は、前記成分値の取り得る値範囲のうちの一部の範囲の成分値の頻度分布が所定の頻度分布となるように成分値を変換することを特徴とする請求項１乃至６のいずれかに記載の画像処理装置。
8. 前記変換手段は、前記成分値の取り得る値範囲のうちの一部の範囲の成分値を不変とすることを特徴とする請求項１乃至７のいずれかに記載の画像処理装置。
9. 画像を複数の周波数帯域の成分に分解する分解工程と、
   前記複数の周波数帯域の少なくとも一つの周波数帯域について、周波数帯域毎の成分値の頻度分布が所定の頻度分布となるように成分値を変換する変換工程と、
   前記変換工程により得られた成分値を用いて画像を生成する生成工程とを備えることを特徴とする画像処理方法。
10. 前記変換工程は、周波数帯毎の成分値の頻度分布が実質的に平坦となるように成分値を変換することを特徴とする請求項９に記載の画像処理方法。
11. 前記変換工程は、前記成分値の累積頻度分布に基づいて前記成分値を変換することを特徴とする請求項１０に記載の画像処理方法。
12. 前記変換工程による前記成分値の変換特性を変更する変更工程を更に備えることを特徴とする請求項９乃至１１のいずれかに記載の画像処理方法。
13. 前記変更工程は、変換後の成分値が一律に所定の割合で変更されるように、前記変換特性を変更することを特徴とする請求項１２に記載の画像処理方法。
14. 前記分解工程はウェーブレット変換を用いて画像を複数の周波数帯域の成分に分解することを特徴とする請求項９乃至１３のいずれかに記載の画像処理方法。
15. 前記変換工程は、前記成分値の取り得る値範囲のうちの一部の範囲の成分値の頻度分布が所定の頻度分布となるように成分値を変換することを特徴とする請求項９乃至１４のいずれかに記載の画像処理方法。
16. 前記変換工程は、前記成分値の取り得る値範囲のうちの一部の範囲の成分値を不変とすることを特徴とする請求項９乃至１５のいずれかに記載の画像処理方法。
17. 請求項９乃至１６のいずれかに記載の画像処理方法をコンピュータに実行させるためのプログラム。
18. 請求項９乃至１６のいずれかに記載の画像処理方法をコンピュータに実行させるためのプログラムを格納する記憶媒体。

Images