(57)【要約】【課題】反射率成分を強調することで画像の視認性を高める際に、照明成分を変化させない場合に比較して、反射率成分の安定性を向上する。【解決手段】画像処理装置１０では、色変換部１１が、原画像を輝度画像及び色度画像に変換する色変換を行い、照明成分推定部１２が、輝度画像から、原画像の照明成分を画素値とする照明成分画像を生成し、安定化反射率算出部１３が、輝度画像と照明成分画像と照明特性シフト量とに基づいて、輝度画像の安定化された反射率成分を画素値とする反射率画像を生成し、輝度再現部１４が、輝度画像と反射率画像とに基づいて輝度再現画像を生成し、色逆変換部１５が、輝度再現画像及び色度画像に対して、色変換部１１が行う色変換とは逆の変換を行う。【選択図】図１
2015226227.tif 000002
【特許請求の範囲】
【請求項１】
原画像を、当該原画像の輝度成分を画素値とする輝度画像及び当該原画像の色度成分を画素値とする色度画像に変換する色変換を行う色変換手段と、
前記輝度画像から、当該輝度画像の照明成分を画素値とする照明画像を生成する照明画像生成手段と、
前記輝度画像と、前記照明画像と、当該照明画像の画素値を変化させる度合いを示す変化度合い情報とに基づいて、当該輝度画像の安定化された反射率成分を画素値とする反射率画像を生成する反射率画像生成手段と、
前記輝度画像と、前記反射率画像とに基づいて、当該輝度画像の反射率成分を強調することにより当該輝度画像の視認性が高まるように再現された輝度再現画像を生成する輝度再現画像生成手段と、
前記輝度再現画像及び前記色度画像に対して、前記色変換手段が行う色変換とは逆の変換を行う色逆変換手段と
を備えたことを特徴とする画像処理装置。
【請求項２】
前記反射率画像生成手段は、前記輝度画像の画素値と、前記変化度合い情報により示される度合いに応じて前記照明画像の画素値を増加させることにより得られた画素値の逆数とに基づいて、当該輝度画像の安定化された反射率成分を算出することを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記原画像に起因する暗部領域の特性に基づいて、前記変化度合い情報を生成する変化度合い情報生成手段を更に備えたことを特徴とする請求項１又は請求項２に記載の画像処理装置。
【請求項４】
前記変化度合い情報生成手段は、前記原画像に起因する暗部領域の特性として、当該暗部領域における画素値のばらつき度合いを用いることを特徴とする請求項３に記載の画像処理装置。
【請求項５】
原画像を当該原画像の再現範囲を縮小するように変更することにより、変更画像を生成する変更画像生成手段と、
前記変更画像を、当該変更画像の輝度成分を画素値とする輝度画像及び当該変更画像の色度成分を画素値とする色度画像に変換する色変換を行う色変換手段と、
前記輝度画像の安定化された反射率成分を強調することにより当該輝度画像の視認性が高まるように再現された輝度再現画像を生成する輝度再現画像生成手段と、
前記輝度再現画像及び前記色度画像に対して、前記色変換手段が行う色変換とは逆の変換を行う色逆変換手段と、
前記色逆変換手段が前記逆の変換を行うことで得られた逆変換画像を当該逆変換画像の再現範囲を拡大するように変更することにより、再現画像を生成する再現画像生成手段と
を備えたことを特徴とする画像処理装置。
【請求項６】
前記変更画像生成手段は、前記原画像の再現範囲を高明度側に縮小し、
前記輝度再現画像生成手段は、前記輝度画像の画素値と、当該輝度画像の照明成分を表す画素値の逆数とに基づいて、当該輝度画像の安定化された反射率成分を算出することを特徴とする請求項５に記載の画像処理装置。
【請求項７】
前記再現画像生成手段は、前記逆変換画像の再現範囲を低明度側に拡大することを特徴とする請求項５又は請求項６に記載の画像処理装置。
【請求項８】
前記変更画像生成手段が前記原画像の再現範囲を縮小する度合いを、当該原画像の暗部領域の特性に基づいて決定する縮小度合い決定手段を更に備えたことを特徴とする請求項５乃至請求項７の何れかに記載の画像処理装置。
【請求項９】
前記縮小度合い決定手段は、前記原画像の暗部領域の特性として、当該暗部領域における画素値のばらつき度合いを用いることを特徴とする請求項８に記載の画像処理装置。
【請求項１０】
前記再現画像生成手段が前記逆変換画像の再現範囲を拡大する度合いを、当該逆変換画像における画素値の出現頻度に基づいて決定する拡大度合い決定手段を更に備えたことを特徴とする請求項５乃至請求項７の何れかに記載の画像処理装置。
【請求項１１】
コンピュータに、
原画像を、当該原画像の輝度成分を画素値とする輝度画像及び当該原画像の色度成分を画素値とする色度画像に変換する色変換を行う機能と、
前記輝度画像から、当該輝度画像の照明成分を画素値とする照明画像を生成する機能と、
前記輝度画像と、前記照明画像と、当該照明画像の画素値を変化させる度合いを示す変化度合い情報とに基づいて、当該輝度画像の安定化された反射率成分を画素値とする反射率画像を生成する機能と、
前記輝度画像と、前記反射率画像とに基づいて、当該輝度画像の反射率成分を強調することにより当該輝度画像の視認性が高まるように再現された輝度再現画像を生成する機能と、
前記輝度再現画像及び前記色度画像に対して、前記色変換とは逆の変換を行う機能と
を実現させるためのプログラム。
【請求項１２】
コンピュータに、
原画像を当該原画像の再現範囲を縮小するように変更することにより、変更画像を生成する機能と、
前記変更画像を、当該変更画像の輝度成分を画素値とする輝度画像及び当該変更画像の色度成分を画素値とする色度画像に変換する色変換を行う機能と、
前記輝度画像の安定化された反射率成分を強調することにより当該輝度画像の視認性が高まるように再現された輝度再現画像を生成する機能と、
前記輝度再現画像及び前記色度画像に対して、前記色変換とは逆の変換を行う機能と、
前記逆の変換を行うことで得られた逆変換画像を当該逆変換画像の再現範囲を拡大するように変更することにより、再現画像を生成する機能と
を実現させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置及びプログラムに関する。
【背景技術】
【０００２】
入力画像の各領域のそれぞれの輝度を第１平滑フィルタによってそれぞれ平滑化した第１平滑化処理信号を取得し、入力画像の各領域の中で輝度が低い領域については、第１平滑フィルタと同じ平滑化を行って調整信号を得、調整信号から第１平滑化処理信号を除算し、入力画像のコントラスト成分信号を算出し、第１平滑化処理信号とコントラスト成分信号とを乗算して合成画像を取得する画像処理装置は知られている（例えば、特許文献１参照）。
【０００３】
逆光画像部分であるか否かを、補正領域判定処理によって生成されるレティネックス画像データと原画像データとの差によって判定し、かかる逆光画像部分にデノイズ処理を施すことにより、ノイズを低減（除去）し、デノイズ処理を行った画像データに対しレティネックス処理を実行する画像処理プログラムも知られている（例えば、特許文献２参照）。
【０００４】
入力画像をフィルタにより照明光成分と反射光成分に分離し、分離された照明光成分を、その明るさを調整する明るさ調整パラメータによって調整することにより、補正後照明光成分を算出し、分離された反射光成分を、その反射率を調整する反射率調整パラメータによって調整することにより、補正後反射率を算出し、補正後照明光成分および補正後反射率に基づいて、出力画像データを生成する画像表示装置も知られている（例えば、特許文献３参照）。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００８−２３６０７２号公報
【特許文献２】特開２００７−６７７４６号公報
【特許文献３】特開２０１３−２１０７０９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００６】
反射率成分を強調することで画像の視認性を高める際に、照明成分を変化させないとすると、反射率成分の安定性が低下してしまう。
【０００７】
そこで、本発明は、反射率成分を強調することで画像の視認性を高める際に、照明成分を変化させない場合に比較して、反射率成分の安定性を向上することを目的とする。
【課題を解決するための手段】
【０００８】
請求項１に記載の発明は、原画像を、当該原画像の輝度成分を画素値とする輝度画像及び当該原画像の色度成分を画素値とする色度画像に変換する色変換を行う色変換手段と、前記輝度画像から、当該輝度画像の照明成分を画素値とする照明画像を生成する照明画像生成手段と、前記輝度画像と、前記照明画像と、当該照明画像の画素値を変化させる度合いを示す変化度合い情報とに基づいて、当該輝度画像の安定化された反射率成分を画素値とする反射率画像を生成する反射率画像生成手段と、前記輝度画像と、前記反射率画像とに基づいて、当該輝度画像の反射率成分を強調することにより当該輝度画像の視認性が高まるように再現された輝度再現画像を生成する輝度再現画像生成手段と、前記輝度再現画像及び前記色度画像に対して、前記色変換手段が行う色変換とは逆の変換を行う色逆変換手段とを備えたことを特徴とする画像処理装置である。
請求項２に記載の発明は、前記反射率画像生成手段は、前記輝度画像の画素値と、前記変化度合い情報により示される度合いに応じて前記照明画像の画素値を増加させることにより得られた画素値の逆数とに基づいて、当該輝度画像の安定化された反射率成分を算出することを特徴とする請求項１に記載の画像処理装置である。
請求項３に記載の発明は、前記原画像に起因する暗部領域の特性に基づいて、前記変化度合い情報を生成する変化度合い情報生成手段を更に備えたことを特徴とする請求項１又は請求項２に記載の画像処理装置である。
請求項４に記載の発明は、前記変化度合い情報生成手段は、前記原画像に起因する暗部領域の特性として、当該暗部領域における画素値のばらつき度合いを用いることを特徴とする請求項３に記載の画像処理装置である。
請求項５に記載の発明は、原画像を当該原画像の再現範囲を縮小するように変更することにより、変更画像を生成する変更画像生成手段と、前記変更画像を、当該変更画像の輝度成分を画素値とする輝度画像及び当該変更画像の色度成分を画素値とする色度画像に変換する色変換を行う色変換手段と、前記輝度画像の安定化された反射率成分を強調することにより当該輝度画像の視認性が高まるように再現された輝度再現画像を生成する輝度再現画像生成手段と、前記輝度再現画像及び前記色度画像に対して、前記色変換手段が行う色変換とは逆の変換を行う色逆変換手段と、前記色逆変換手段が前記逆の変換を行うことで得られた逆変換画像を当該逆変換画像の再現範囲を拡大するように変更することにより、再現画像を生成する再現画像生成手段とを備えたことを特徴とする画像処理装置である。
請求項６に記載の発明は、前記変更画像生成手段は、前記原画像の再現範囲を高明度側に縮小し、前記輝度再現画像生成手段は、前記輝度画像の画素値と、当該輝度画像の照明成分を表す画素値の逆数とに基づいて、当該輝度画像の安定化された反射率成分を算出することを特徴とする請求項５に記載の画像処理装置である。
請求項７に記載の発明は、前記再現画像生成手段は、前記逆変換画像の再現範囲を低明度側に拡大することを特徴とする請求項５又は請求項６に記載の画像処理装置である。
請求項８に記載の発明は、前記変更画像生成手段が前記原画像の再現範囲を縮小する度合いを、当該原画像の暗部領域の特性に基づいて決定する縮小度合い決定手段を更に備えたことを特徴とする請求項５乃至請求項７の何れかに記載の画像処理装置である。
請求項９に記載の発明は、前記縮小度合い決定手段は、前記原画像の暗部領域の特性として、当該暗部領域における画素値のばらつき度合いを用いることを特徴とする請求項８に記載の画像処理装置である。
請求項１０に記載の発明は、前記再現画像生成手段が前記逆変換画像の再現範囲を拡大する度合いを、当該逆変換画像における画素値の出現頻度に基づいて決定する拡大度合い決定手段を更に備えたことを特徴とする請求項５乃至請求項７の何れかに記載の画像処理装置である。
請求項１１に記載の発明は、コンピュータに、原画像を、当該原画像の輝度成分を画素値とする輝度画像及び当該原画像の色度成分を画素値とする色度画像に変換する色変換を行う機能と、前記輝度画像から、当該輝度画像の照明成分を画素値とする照明画像を生成する機能と、前記輝度画像と、前記照明画像と、当該照明画像の画素値を変化させる度合いを示す変化度合い情報とに基づいて、当該輝度画像の安定化された反射率成分を画素値とする反射率画像を生成する機能と、前記輝度画像と、前記反射率画像とに基づいて、当該輝度画像の反射率成分を強調することにより当該輝度画像の視認性が高まるように再現された輝度再現画像を生成する機能と、前記輝度再現画像及び前記色度画像に対して、前記色変換とは逆の変換を行う機能とを実現させるためのプログラムである。
請求項１２に記載の発明は、コンピュータに、原画像を当該原画像の再現範囲を縮小するように変更することにより、変更画像を生成する機能と、前記変更画像を、当該変更画像の輝度成分を画素値とする輝度画像及び当該変更画像の色度成分を画素値とする色度画像に変換する色変換を行う機能と、前記輝度画像の安定化された反射率成分を強調することにより当該輝度画像の視認性が高まるように再現された輝度再現画像を生成する機能と、前記輝度再現画像及び前記色度画像に対して、前記色変換とは逆の変換を行う機能と、前記逆の変換を行うことで得られた逆変換画像を当該逆変換画像の再現範囲を拡大するように変更することにより、再現画像を生成する機能とを実現させるためのプログラムである。
【発明の効果】
【０００９】
請求項１の発明によれば、反射率成分を強調することで画像の視認性を高める際に、照明成分を変化させない場合に比較して、反射率成分の安定性を向上することができる。
請求項２の発明によれば、反射率成分を照明成分の逆数に基づいて算出する際に、照明成分を変化させない場合に比較して、反射率成分の安定性を向上することができる。
請求項３の発明によれば、反射率成分の安定性を向上するための照明成分の変化の度合いを、原画像に起因する暗部領域の特性に基づいて決定することができる。
請求項４の発明によれば、反射率成分の安定性を向上するための照明成分の変化の度合いを、原画像に起因する暗部領域における画素値のばらつき度合いに基づいて決定することができる。
請求項５の発明によれば、反射率成分を強調することで画像の視認性を高める際に、画像の再現範囲を変更しない場合に比較して、反射率成分の安定性を向上することができる。
請求項６の発明によれば、反射率成分を照明成分の逆数に基づいて算出する際に、画像の再現範囲を変更しない場合に比較して、反射率成分の安定性を向上することができる。
請求項７の発明によれば、原画像の再現範囲で反射率成分の安定性を向上することができる。
請求項８の発明によれば、反射率成分の安定性を向上するための再現範囲の縮小の度合いを、原画像の暗部領域の特性に基づいて決定することができる。
請求項９の発明によれば、反射率成分の安定性を向上するための再現範囲の縮小の度合いを、原画像の暗部領域における画素値のばらつき度合いに基づいて決定することができる。
請求項１０の発明によれば、原画像の再現範囲で反射率成分の安定性を向上するための再現範囲の拡大の度合いを、画素値の出現頻度に基づいて決定することができる。
請求項１１の発明によれば、反射率成分を強調することで画像の視認性を高める際に、照明成分を変化させない場合に比較して、反射率成分の安定性を向上することができる。
請求項１２の発明によれば、反射率成分を強調することで画像の視認性を高める際に、画像の再現範囲を変更しない場合に比較して、反射率成分の安定性を向上することができる。
【図面の簡単な説明】
【００１０】
【図１】本発明の第１の実施の形態における画像処理装置の機能構成例を示したブロック図である。
【図２】Ｌ（ｘ，ｙ）が変化したときのＬ（ｘ，ｙ）の逆数の変化を示したグラフである。
【図３】Ｌ（ｘ，ｙ）が変化したときの（Ｌ（ｘ，ｙ）＋Ｌｏ）の逆数の変化をＬ（ｘ，ｙ）の逆数の変化と比較して示したグラフである。
【図４】反射率が不安定な場合に原画像に基づいて再現される視認性向上画像及びその暗部領域を示した図である。
【図５】反射率が安定している場合に原画像に基づいて再現される視認性向上画像を示した図である。
【図６】本発明の第１の実施の形態における画像処理装置の動作例を示したフローチャートである。
【図７】本発明の第２の実施の形態における画像処理装置の機能構成例を示したブロック図である。
【図８】（ａ）は視認性の悪い風景画像の例を示した図であり、（ｂ）は視認性の悪い展示画像の例を示した図である。
【図９】（ａ）は風景画像に適した照明特性シフト量により視認性が向上した風景画像の例を示した図であり、（ｂ）は風景画像に適した照明特性シフト量により暗部ノイズが発生した展示画像の例を示した図である。
【図１０】（ａ）は展示画像に適した照明特性シフト量により視認性が向上しなかった風景画像の例を示した図であり、（ｂ）は展示画像に適した照明特性シフト量により視認性が向上した展示画像の例を示した図である。
【図１１】（ａ）は色のばらつきが大きい画像の例を示した図であり、（ｂ）は色のばらつきが小さい画像の例を示した図である。
【図１２】ＲＧＢの分散の平均値σ２ＲＧＢと照明特性シフト量Ｌｏとを対応付けた関数の例を示した図である。
【図１３】（ａ）は視認性が向上された風景画像の例を示した図であり、（ｂ）は視認性が向上された展示画像の例を示した図である。
【図１４】本発明の第２の実施の形態の変形例における画像処理装置の機能構成例を示したブロック図である。
【図１５】本発明の第３の実施の形態における画像処理装置の機能構成例を示したブロック図である。
【図１６】原画像の色成分ＲＧＢをシフト画像の色成分Ｒ’Ｇ’Ｂ’に変換する関数の例を示した図である。
【図１７】原画像の色域の範囲とシフト画像の色域の範囲とをＨＳＶ色空間内に示した図である。
【図１８】本発明の第３の実施の形態における画像処理装置の視認性再現部の機能構成例を示したブロック図である。
【図１９】本発明の第３の実施の形態で用いるＬ（ｘ，ｙ）の範囲を示した図である。
【図２０】シフト再現画像の色成分Ｒ’Ｇ’Ｂ’をカラー再現画像の色成分Ｒ”Ｇ”Ｂ”に変換する関数の例を示した図である。
【図２１】本発明の第３の実施の形態における画像処理装置の動作例を示したフローチャートである。
【図２２】本発明の第４の実施の形態における画像処理装置の機能構成例を示したブロック図である。
【図２３】ＲＧＢの分散の平均値σ２ＲＧＢとカラー特性シフト量ｄ１とを対応付けた関数の例を示した図である。
【図２４】シフト再現画像におけるＲＧＢの画素頻度の例を示した図である。
【図２５】シフト再現画像におけるＲＧＢの画素頻度上にカラー特性逆シフト量ｄ２の計算例を示した図である。
【図２６】シフト再現画像の再現範囲を拡大したときのＲＧＢの画素頻度の例を示した図である。
【図２７】本発明の実施の形態における画像処理装置のハードウェア構成例を示したブロック図である。
【発明を実施するための形態】
【００１１】
以下、添付図面を参照して、本発明の実施の形態について詳細に説明する。
【００１２】
［発明の背景］
近年、ＩＣＴ（Information and Communication Technology）デバイスを持ち歩くユーザが増え、様々な場所で写真の撮影や画像の閲覧を行うことが多くなった。ＩＣＴデバイスには、場所を問わず手軽に画像の閲覧等を行えるという利点があるが、環境光の変化等、持ち運んだ先に大きく左右されるという欠点もある。こうしたことから、近年のモニタ環境では、従来のモニタ環境との違いとして「手軽に使え」「使用場所が多様」であることが挙げられるものの、画像が見易いかどうかの「視認性」の方が重視されてきている。
【００１３】
「視認性」とは、視対象がはっきり見えるか否かの特性である。画像の視認性を高める方法には、ガンマ補正、ヒストグラムイコライゼーション等がある。ガンマ補正では、暗部や対象となる領域を盛り上げるカーブを生成し、画素値に適用することで、暗部を明るくする。ヒストグラムイコライゼーションでは、画像のヒストグラムの偏りをなくすカーブを生成し、画素値に適用することで、ヒストグラムが平滑化される再現を行う。
【００１４】
また、人間の視覚認知特性を利用した視認性向上の方法には、レティネックス原理を利用したものもある。レティネックス原理に基づく再現は、人間のシーンの知覚は、物体の反射率によるところが大きいという心理物理学的な知見を応用したものである。具体的には、レティネックスモデルによる画像再現処理では、画素値を照明成分と反射率成分とに分解し、反射率成分を強調することで視認性の向上を行う。
【００１５】
このように、レティネックス原理に基づく再現は、画像の視認性を高める。しかしながら、その一方で、暗部ノイズを強調するという問題がある。
【００１６】
かかる問題に対し、従来は、発生するノイズに対して対策を講じていたが、本実施の形態では、レティネックスモデルの計算原理に対して対策を講ずる。レティネックス原理では、計算が不安定になる部分があり、これがノイズ発生の原因にもなる。本実施の形態は、視認性を向上させる再現を行う際に、ノイズを抑え、かつ、画質を安定化するものである。
【００１７】
［第１の実施の形態］
図１は、本発明の第１の実施の形態における画像処理装置１０の機能構成例を表すブロック図である。図示するように、第１の実施の形態における画像処理装置１０は、色変換部１１と、照明成分推定部１２と、安定化反射率算出部１３と、輝度再現部１４と、色逆変換部１５とを備える。
【００１８】
色変換部１１は、原画像であるカラー画像を、輝度成分を画素値とする輝度画像と、色度成分を画素値とする色度画像とに分離する。カラー画像は一般にＲＧＢやＣＹＭＫ等の色成分で表された画像である。一方、輝度画像及び色度画像としては、ＨＳＶ、Ｌ＊ａ＊ｂ＊、ＹＣｂＣｒ等で表された画像が挙げられる。本実施の形態では、輝度画像及び色度画像をＨＳＶ画像として説明を行う。このように、輝度画像及び色度画像をＨＳＶ画像とすると、色変換部１１は、輝度画像としてＶ画像を生成し、色度画像としてＨ（色相）Ｓ（彩度）画像を生成する。
【００１９】
本実施の形態では、原画像を輝度画像及び色度画像に変換する色変換を行う色変換手段の一例として、色変換部１１を設けている。
【００２０】
ここで、Ｖ画像の位置（ｘ，ｙ）に対する画素値をＩ（ｘ，ｙ）とすると、レティネックスモデルでは、以下のように、Ｉ（ｘ，ｙ）を照明成分Ｌ（ｘ，ｙ）と反射率成分Ｒ（ｘ，ｙ）とに分ける。
【００２１】
【数１】
2015226227.tif 000003
【００２２】
数式１のように１つの値を２つの値に分離する問題は、一般に不良設定問題と呼ばれる。何らかの方法を用いて、照明成分及び反射率成分の何れか一方を推定しなければ、もう一方は求めることができない。
【００２３】
そこで、照明成分推定部１２は、色変換部１１で生成されたＩ（ｘ，ｙ）に基づいて、Ｌ（ｘ，ｙ）を推定する。視覚認知の１つの特性として、網膜に入った光の１点（画像では１つの画素）の知覚量は、その１点の周辺輝度の平均に影響することが知られている。周辺輝度は推定された照明成分に相当し、照明推定モデルは以下の関数のコンボリューションとなる。
【００２４】
【数２】
2015226227.tif 000004
【００２５】
ここで、ｋは画像処理のフィルタサイズの画素分で積分した場合に結果が１になるように正規化する係数を表す。また、σは平滑化度合い（スケール）を表し、σが大きいほどぼかしの強いコンボリューション画像となる。尚、上記の関数は一例であり、結果として画像が平滑化されるフィルタであれば如何なるものを用いてもよい。例えば、数式１を変形した関数によるフィルタで、エッジ保存を行う平滑化フィルタとして知られるバイラテラルフィルタがあるが、これを用いてもよいものとする。
【００２６】
レティネックス原理によれば、人の視覚特性は、注目領域の周辺から照明光を推定していることが知られている。従って、平滑化した画像は、推定された照明光を表すことになる。しかしながら、シーンによって適するスケールが異なるため、例えば、照明光の推定は、以下のようにスケール１からスケールＮまでのＮ層の画像の加重合計をとるのが望ましい。
【００２７】
【数３】
2015226227.tif 000005
【００２８】
ここで、Ｌ（ｘ，ｙ）は推定される照明成分を表し、Ｇｎ（ｘ，ｙ）はスケールｎに対する数式２を表し、Ｉ（ｘ，ｙ）は輝度画像の画素値を表し、Ｗｎはスケールｎに対する重みを表し、「×」を「○」で囲んだ記号は畳み込みを表す。尚、Ｗｎは、簡易的に１／Ｎとしてもよいし、層に応じて可変としてもよい。そして、このように算出されたＬ（ｘ，ｙ）を画素値とする画像を照明成分画像とする。
【００２９】
本実施の形態では、輝度画像の照明成分を画素値とする照明画像の一例として、照明成分画像を用いており、照明画像を生成する照明画像生成手段の一例として、照明成分推定部１２を設けている。
【００３０】
ところで、数式１によれば、推定された照明成分に基づいて、以下のように反射率が算出される。
【００３１】
【数４】
2015226227.tif 000006
【００３２】
数式４から分かるように、Ｒ（ｘ，ｙ）はＬ（ｘ，ｙ）の逆数を用いた計算により求められる。そこで、Ｌ（ｘ，ｙ）をＬとし、Ｌ（ｘ，ｙ）の逆数をＬ’として、Ｌが変化したときのＬ’の変化について考える。図２に、「Ｌ’＝１／Ｌ」のグラフを示す。図２に破線で囲んで示すように、Ｌが小さい領域、つまり、照明成分の値が小さい領域（周辺輝度が低い領域）では、Ｌの少しの変化に対して、Ｌ’が大きく変化する。画素値Ｉ（ｘ，ｙ）が０又は０に近い値であれば、数式４が０になるように上手く計算できることもある。しかしながら、Ｉ（ｘ，ｙ）及びＬ（ｘ，ｙ）が共に小さく、１／Ｌが発散するような場合は、反射率の計算が破綻してしまう。
【００３３】
そこで、安定化反射率算出部１３は、安定した反射率を算出するための照明特性シフト量により、反射率を算出する。具体的には、以下のように、色変換部１１で生成された輝度画像Ｉ（ｘ，ｙ）と、照明成分推定部１２で生成された照明成分画像Ｌ（ｘ，ｙ）とに加え、照明特性シフト量Ｌｏを用いて、安定化された反射率（以下、「安定化反射率」という）Ｒｓを算出する。
【００３４】
【数５】
2015226227.tif 000007
【００３５】
ここで、Ｌｏは画像の特性に応じて変えることのできるパラメータである。また、この安定化反射率Ｒｓを画素値とする画像が反射率画像である。
【００３６】
さて、数式５についても、Ｌ（ｘ，ｙ）をＬとし、（Ｌ（ｘ，ｙ）＋Ｌｏ）の逆数をＬ’として、Ｌが変化したときのＬ’の変化について考える。図３に、「Ｌ’＝１／（Ｌ＋Ｌｏ）」のグラフを、「Ｌ’＝１／Ｌ」のグラフと比較して示している。図３から分かるように、照明特性シフト量Ｌｏを用いることで、「Ｌ’＝１／（Ｌ＋Ｌｏ）」のグラフは、「Ｌ’＝１／Ｌ」のグラフを左側にシフトしたものとなっている。これにより、Ｌが小さい領域でも、Ｌの少しの変化に対して、Ｌ’が大きく変化することがない。
【００３７】
本実施の形態では、照明画像の画素値を変化させる度合いを示す変化度合い情報の一例として、照明特性シフト量を用いており、輝度画像と照明画像と変化度合い情報とに基づいて反射率画像を生成する反射率画像生成手段の一例として、安定化反射率算出部１３を設けている。
【００３８】
輝度再現部１４は、安定化反射率算出部１３が算出した安定化反射率Ｒｓを強調することにより、色変換部１１で生成された輝度画像の視認性を高めた再現画像である輝度再現画像を生成する。
一般に、輝度再現画像の画素値Ｉ＾（ｘ，ｙ）は、例えば、以下の式で求められる。尚、本明細書では、ハット記号を、数式中では文字の真上に付すが、文中では文字の後ろに付すものとする。
【００３９】
【数６】
2015226227.tif 000008
【００４０】
ここで、αは反射率の強調度合いを示すパラメータであり、０≦α≦１とする。Ｉ＾（ｘ，ｙ）は、α＝０の場合、輝度画像の画素値を表し、α＝１の場合、反射率画像の画素値を表す。
【００４１】
これに対し、本実施の形態の安定化反射率Ｒｓ（ｘ，ｙ）による再現式は、以下のようなものとすればよい。
【００４２】
【数７】
2015226227.tif 000009
【００４３】
例えば、数式６により視認性を高める再現を行うと、不安定な反射率の影響により、その暗部領域が図４の破綻部拡大図に示すような状態になることがある。即ち、不安定な計算が暗部のノイズの原因となっている。これに対し、本実施の形態によれば、逆数の計算が安定することにより、反射率も安定化されるので、図５に示すように、暗部ノイズの影響が抑えられる。
【００４４】
本実施の形態では、輝度再現画像を生成する輝度再現画像生成手段の一例として、輝度再現部１４を設けている。
【００４５】
色逆変換部１５は、このようにして視認性を高めた輝度再現画像が生成されると、この輝度再現画像と、色変換部１１で生成された色度画像とに対して、色変換部１１とは逆の色変換を行う。即ち、この第１の実施の形態の一連の処理で得られたＨＳＶ色空間の再現画像を、ＲＧＢ色空間の再現画像に変換することにより、最終的な再現画像であるカラー再現画像を得る。
【００４６】
本実施の形態では、色変換手段が行う色変換とは逆の変換を行う色逆変換手段の一例として、色逆変換部１５を設けている。
【００４７】
図６は、本発明の第１の実施の形態における画像処理装置１０の動作例を表すフローチャートである。
【００４８】
カラー画像が入力されると、まず、色変換部１１が、カラー画像に対してカラー画像の色空間から輝度及び色度の色空間への色変換を行うことにより、輝度画像及び色度画像を生成する（ステップ１０１）。
【００４９】
次に、照明成分推定部１２が、ステップ１０１で生成された輝度画像に基づいて、数式３により照明成分画像を生成する（ステップ１０２）。次いで、安定化反射率算出部１３が、ステップ１０１で生成された輝度画像と、ステップ１０２で生成された照明成分画像と、予め定められた照明特性シフト量Ｌｏとに基づいて、数式５により反射率画像を生成する（ステップ１０３）。その後、輝度再現部１４が、ステップ１０１で生成された輝度画像と、ステップ１０３で生成された反射率画像とに基づいて、数式７により輝度再現画像を生成する（ステップ１０４）。
【００５０】
最後に、色逆変換部１５が、ステップ１０４で生成された輝度再現画像と、ステップ１０１で生成された色度画像とに対して、色変換部１１が行ったのとは逆の色変換、つまり、輝度及び色度の色空間からカラー画像の色空間への色変換を行うことにより、カラー再現画像を生成する（ステップ１０５）。
【００５１】
［第２の実施の形態］
図７は、本発明の第２の実施の形態における画像処理装置１０の機能構成例を表すブロック図である。図示するように、第２の実施の形態における画像処理装置１０は、色変換部１１と、照明成分推定部１２と、安定化反射率算出部１３と、輝度再現部１４と、色逆変換部１５と、暗部特性解析部１６とを備える。ここで、色変換部１１、照明成分推定部１２、安定化反射率算出部１３、輝度再現部１４、及び色逆変換部１５については、第１の実施の形態と同じなので説明を省略し、以下では、暗部特性解析部１６についてのみ説明する。
【００５２】
暗部特性解析部１６は、カラー画像の暗部特性から、画像に応じた照明特性シフト量Ｌｏを算出する。
【００５３】
画像によっては、照明特性シフト量Ｌｏを変更した方が望ましい場合がある。図８（ａ）は視認性の悪い風景画像の例であり、図８（ｂ）は視認性の悪い展示画像の例である。図８（ａ）の風景画像では、特に木の葉の部分が暗部であり、図８（ｂ）の展示画像では、特に壁の部分が暗部であるとする。このような視認性の悪い画像では、壁等の均一な領域が暗部の場合、及び、木の葉等のテクスチャが存在する領域が暗部の場合の両方において、暗部は殆ど何も見えない状態になっている。
【００５４】
このような場合に、図８（ａ）の風景画像及び図８（ｂ）の展示画像に対して、同一の照明特性シフト量Ｌｏを用いて視認性を高める再現を行ったとする。すると、例えば、図９（ａ）に示すように、照明特性シフト量Ｌｏが０に近くても視認性が良くなるような風景画像では、木の葉がくっきりと見えるようになる。一方で、図９（ｂ）に示すように、展示画像の壁においては、暗部ノイズが引き立っている。また、図１０（ｂ）に示すように、照明特性シフト量Ｌｏを、展示画像で暗部ノイズが出ないような値に合わせたとする。この場合は、図１０（ａ）に示すように、風景画像で木の葉が見え難い等の不都合が生じる。即ち、同一の照明特性シフト量Ｌｏで両方の画像の暗部ノイズを除去することができない場合がある。
【００５５】
これは、同じ暗部領域でも、図１１（ａ）に示す木の葉のようなテクスチャ領域と、図１１（ｂ）に示す壁のような領域とでは、色のばらつき方の特性に違いがあるからだと考えられる。即ち、どちらも暗部であるが、色のばらつきの範囲、ばらつき度合い等が明らかに異なっている。尚、図１１（ａ），（ｂ）は何れも暗部領域の図であるが、特性を把握し易くするため、明度を上げて示している。
【００５６】
そこで、暗部特性解析部１６は、図１１（ａ），（ｂ）に示すような特性に合わせた照明特性シフト量Ｌｏを算出する。例えば、図１１（ａ）は、暗部にノイズではないテクスチャがある領域を示しており、このような領域では、ＲＧＢ値は０から３０ほどの範囲でばらついている。また、図１１（ｂ）は、レティネックス強調を行うと明らかにノイズが増幅する領域、つまり、本来は視認性を向上する再現を行っても画素のばらつきがあってほしくない領域を示しており、このような領域では、ＲＧＢ値は０から１５ほどの範囲でばらついている。
【００５７】
以下、暗部特性解析部１６によるカラー画像の解析例を示す。
【００５８】
まず、暗部特性解析部１６は、以下の式を満たす画素を抽出する。
【００５９】
【数８】
2015226227.tif 000010
【００６０】
ここで、θＲＧＢは暗部とみなすＲＧＢ値の閾値である。
【００６１】
次に、暗部特性解析部１６は、以下のように、Ｒ値、Ｇ値、Ｂ値のそれぞれの分散を算出し、これらの分散の平均値σ２ＲＧＢを求める。
【００６２】
【数９】
2015226227.tif 000011
【００６３】
ここで、Ｎは数式８を満たす暗部に含まれる画素の数とする。数式９で求めたσ２ＲＧＢを目安として、暗部でどれだけ画素値がばらついているかが分かる。図１１（ａ），（ｂ）の例で示した通り、この値が小さければ、壁のような領域と考えられるので、視認性を高める再現を行った後にノイズが目立つ可能性が高くなる。また、逆に、この値が大きければ、木の葉のようなテクスチャと考えられる。
【００６４】
σ２ＲＧＢが小さい場合は、計算安定性を確保するために、照明特性シフト量Ｌｏを大きくすることが望ましい。また、σ２ＲＧＢが大きくなるにつれて、照明特性シフト量Ｌｏを小さくすればよい。従って、例えば、図１２のように、σ２ＲＧＢとＬｏとを対応付ければよい。ここで、図１２中のｔＲＧＢをどの程度の値とすればよいか等は、壁とみなされる領域やテクスチャである領域の統計値から予め決めておけばよい。このようにシーンに応じて決定された照明特性シフト量Ｌｏを用いることで、図１３（ａ），（ｂ）に示すように、どちらの画像に対しても良好な視認性が確保される。
【００６５】
ところで、これまでは、暗部特性解析部１６がカラー画像の暗部を解析した結果に基づいて照明特性シフト量Ｌｏを決めるものとして説明したが、この限りではない。より一般化して、カラー画像に起因する暗部を解析した結果に基づいて照明特性シフト量Ｌｏを決めるものとしてもよい。一例として、図１４に示すように、暗部特性解析部１６は、色変換部１１で生成された輝度画像及び色度画像の暗部を解析した結果に基づいて照明特性シフト量Ｌｏを決めるものであってもよい。この場合は、輝度Ｖ、色度ＨＳの値を用いて、以下のようにσ２Ｖ、σ２Ｈ、σ２Ｓを求める。
【００６６】
【数１０】
2015226227.tif 000012
【００６７】
そして、これらの値を数式９のＲＧＢの例のように組み合わせる等して、Ｌｏとの対応付けをとればよい。具体的には、σ２Ｖ、σ２Ｈ、σ２Ｓの平均値σ２ＨＳＶを求め、σ２ＨＳＶを数式９の平均値σ２ＲＧＢのようにＬｏと対応付ければよい。或いは、σ２Ｖ、σ２Ｈ、σ２Ｓの重み付け平均値をσ２ＨＳＶとしてもよいし、σ２Ｖ、σ２Ｈ、σ２Ｓの何れか１つをσ２ＨＳＶとしてもよい。本実施の形態において、暗部特性解析部１６は、暗部の画素のばらつき度合いとシフト量とを対応付けるものであるので、この意義を損なわないものであれば如何なるものでもよいものとする。
【００６８】
尚、この第２の実施の形態では、暗部特性解析部１６が、暗部における画素値のばらつき度合いを暗部の特性として求め、これに応じて照明特性シフト量Ｌｏを決定するようにしたが、これには限らない。例えば、暗部における彩度ばらつき度合い等の他の特性に応じて照明特性シフト量Ｌｏを決定するようにしてもよい。
【００６９】
本実施の形態では、変化度合い情報を生成する変化度合い情報生成手段の一例として、暗部特性解析部１６を設けている。
【００７０】
［第３の実施の形態］
第３の実施の形態も、第１の実施の形態及び第２の実施の形態と解決する課題は同じで、暗部ノイズを引き立たせないことである。第１の実施の形態及び第２の実施の形態では、照明成分の特性に対して、計算安定性を確保する処理を行っていたが、第３の実施の形態では、原画像であるカラー画像に対して、後の反射率成分の計算が安定して行えるように特性を変える処理を行う。
【００７１】
図１５は、本発明の第３の実施の形態における画像処理装置１０の機能構成例を表すブロック図である。図示するように、第３の実施の形態における画像処理装置１０は、カラー特性シフト部２１と、色変換部２２と、視認性再現部２３と、色逆変換部２４と、カラー特性逆シフト部２５とを備える。ここで、色変換部２２は、第１の実施の形態及び第２の実施の形態の色変換部１１でカラー画像を後述するシフト画像に置き換えたに過ぎず、色逆変換部２４は、第１の実施の形態及び第２の実施の形態の色逆変換部１５でカラー再現画像を後述するシフト再現画像に置き換えたにすぎない。従って、これらの詳細な説明は省略し、以下では、カラー特性シフト部２１、視認性再現部２３及びカラー特性逆シフト部２５についてのみ説明する。
【００７２】
カラー特性シフト部２１は、カラー特性シフト量ｄ１を用いて、カラー画像の色成分に対して、色空間の色域を縮小するようなシフトを行う。図１６に、原画像の色成分をＲＧＢとし、シフトを行った後のシフト画像の色成分をＲ’Ｇ’Ｂ’とした場合のＲＧＢとＲ’Ｇ’Ｂ’との対応付けを示す。ＲＧＢ及びＲ’Ｇ’Ｂ’の範囲は０〜２５５なので、ｍａｘは２５５とすればよい。図１６は式で表すと以下のようになる。
【００７３】
【数１１】
2015226227.tif 000013
【００７４】
このようにすることで、元のカラー画像の階調性を保ったまま、後の視認性再現の計算を暗部領域がない状態で行う。尚、シフト画像の色域の範囲は、図１６の変換を行うことにより、例えば、図１７に太線で示される範囲となり、原画像の色域の範囲よりも高明度側に縮小されていることが分かる。ここで、図１７はＨＳＶ画像における色域を表す図である。
【００７５】
本実施の形態では、変更画像の一例として、シフト画像を用いており、変更画像を生成する変更画像生成手段の一例として、カラー特性シフト部２１を設けている。
【００７６】
視認性再現部２３は、色変換部２２で生成された輝度画像に対して、視認性を向上するための画像処理を行うことにより、輝度再現画像を生成する。
【００７７】
本実施の形態では、輝度再現画像を生成する輝度再現画像生成手段の一例として、視認性再現部２３を設けている。
【００７８】
図１８は、視認性再現部２３の機能構成例を示したブロック図である。照明成分推定部２３１は、第１の実施の形態における照明成分推定部１２と同じでよい。また、反射率算出部２３２は、数式４の計算を実行するものであればよく、輝度再現部２３３は、数式６の計算を実行するものであればよい。
【００７９】
このように、反射率算出部２３２は、数式４の計算を実行するが、本実施の形態では、図１６の変換を行っているので、反射率算出部２３２で用いるＬ（ｘ，ｙ）の範囲は、図１９に「照明成分の範囲」として示した範囲となる。これにより、図１９の「Ｌ’＝１／Ｌ」がとり得るＬ’の値のうち、「照明成分の範囲」内のＬに対応するＬ’の値のみが反射率に寄与することになる。従って、安定性確保という意味では、第１の実施の形態のように照明成分をシフトする場合と類似した効果がある。
【００８０】
色逆変換部２４は、このようにして得られた輝度再現画像（Ｖ画像）と、色変換部２２で得られた色度画像（ＨＳ画像）とに対して、ＨＳＶ色空間からＲＧＢ色空間への色変換を行うことにより、シフト再現画像を生成する。このとき、シフト再現画像の色域の範囲は、図１７に示すシフト画像の色域の範囲のままである。
【００８１】
そこで、カラー特性逆シフト部２５は、カラー特性逆シフト量ｄ２を用いて、シフト再現画像の色成分に対して、図２０に示すような逆シフトを行う。
【００８２】
図２０が表す直線は、０≦ＲＧＢ＜ｄ２の範囲では、以下のようにすればよい。
【００８３】
【数１２】
2015226227.tif 000014
【００８４】
また、ｄ２≦ＲＧＢ≦２５５の範囲では、以下のようにすればよい。
【００８５】
【数１３】
2015226227.tif 000015
【００８６】
このような式で、色域の範囲が低明度側に、つまり、図１７における原画像の色域の範囲の方向に広げられる。これにより、視認性が高められ、かつ、計算安定性を確保した反射率によりノイズが抑制されたカラー再現画像が得られる。
【００８７】
本実施の形態では、逆変換画像の一例として、シフト再現画像を用いており、再現画像の一例として、カラー再現画像を用いており、逆変換画像をその再現範囲を拡大するように変更することにより再現画像を生成する再現画像生成手段の一例として、カラー特性逆シフト部２５を設けている。
【００８８】
図２１は、本発明の第３の実施の形態における画像処理装置１０の動作例を表すフローチャートである。
【００８９】
カラー画像が入力されると、まず、カラー特性シフト部２１が、カラー画像に対して図１６の変換を行うことにより、シフト画像を生成する（ステップ２０１）。
【００９０】
すると、色変換部２２が、シフト画像に対してシフト画像の色空間から輝度及び色度の色空間への色変換を行うことにより、輝度画像及び色度画像を生成する（ステップ２０２）。
【００９１】
次に、視認性再現部２３が、ステップ２０２で生成された輝度画像に基づいて、数式３、数式４及び数式６により輝度再現画像を生成する（ステップ２０３）。
【００９２】
その後、色逆変換部２４が、ステップ２０３で生成された輝度再現画像と、ステップ２０２で生成された色度画像とに対して、色変換部２２が行ったのとは逆の色変換、つまり、輝度及び色度の色空間からシフト画像の色空間への色変換を行うことにより、シフト再現画像を生成する（ステップ２０４）。
【００９３】
最後に、カラー特性逆シフト部２５が、ステップ２０４で生成されたシフト再現画像に対して図２０の変換を行うことにより、カラー再現画像を生成する（ステップ２０５）。
【００９４】
［第４の実施の形態］
図２２は、本発明の第４の実施の形態における画像処理装置１０の機能構成例を表すブロック図である。図示するように、第４の実施の形態における画像処理装置１０は、カラー特性シフト部２１と、色変換部２２と、視認性再現部２３と、色逆変換部２４と、カラー特性逆シフト部２５と、暗部特性解析部２６と、頻度解析部２７とを備える。ここで、カラー特性シフト部２１、色変換部２２、視認性再現部２３、色逆変換部２４、及びカラー特性逆シフト部２５については、第３の実施の形態と同じなので説明を省略し、以下では、暗部特性解析部２６及び頻度解析部２７についてのみ説明する。
【００９５】
暗部特性解析部１６は、カラー特性シフト部２１で用いるカラー特性シフト量ｄ１をカラー画像の暗部特性から計算する。この暗部特性の解析の基本的な考え方は、第２の実施の形態における考え方と同じである。例えば、第２の実施の形態で述べたように、数式８及び数式９によりσ２ＲＧＢを算出する。σ２ＲＧＢが大きければ、視認性を高めてほしいテクスチャ領域である可能性が高い。また、σ２ＲＧＢが小さければ、計算安定性を確保してほしい領域である可能性が高い。従って、カラー特性シフト量ｄ１は、例えば図２３のような対応関係により求められる。ここで、図２３中のｔＲＧＢをどの程度の値とすればよいか等は、第２の実施の形態と同じく、壁とみなされる領域やテクスチャである領域の統計値から予め決めておけばよい。
【００９６】
本実施の形態では、原画像の再現範囲を縮小する度合いを決定する縮小度合い決定手段の一例として、暗部特性解析部２６を設けている。
【００９７】
頻度解析部２７は、再現範囲が図１７に太線で示した範囲であるシフト再現画像の状態を参照し、範囲を広げる量であるカラー特性逆シフト量ｄ２を算出する。例えば、色逆変換部２４でＨＳＶ色空間からＲＧＢ色空間へ変換された場合のＲＧＢの画素頻度（ＲＧＢヒストグラム）が図２４のようになっていたとする。この場合は、例えば、以下のように、Ｒの最小値、Ｇの最小値、Ｂの最小値のうちの最小値ｍｉｎＲＧＢを求める。
【００９８】
【数１４】
2015226227.tif 000016
【００９９】
そして、ｍｉｎＲＧＢを用いて、以下のように、ｄ２を求める。
【０１００】
【数１５】
2015226227.tif 000017
【０１０１】
ここで、ｄ２は、画質設計にもよるが、数式２１のように、ｍｉｎＲＧＢに基づいて決定するとよい。図２５は、ｍｉｎＲＧＢからΔｄを減算することによりｄ２を求める場合の例である。
【０１０２】
このようにして求めたカラー特性逆シフト量ｄ２を用いて図２０のように変換すると、ＲＧＢの画素頻度は図２６に示すようなものとなる。これにより、原画像の色域範囲で、視認性を高め、かつ、ノイズの発生を抑えた視認性再現が可能となる。
【０１０３】
本実施の形態では、逆変換画像の再現範囲を拡大する度合いを決定する拡大度合い決定手段の一例として、頻度解析部２７を設けている。
【０１０４】
［第５の実施の形態］
本発明の第５の実施の形態は、第３の実施の形態及び第４の実施の形態に示す視認性再現部２３を、第１の実施の形態又は第２の実施の形態の照明成分推定部１２、安定化反射率算出部１３及び輝度再現部１４としたものである。
【０１０５】
第１の実施の形態及び第２の実施の形態と、第３の実施の形態及び第４の実施の形態とは、同じ計算安定性の効果をもたらすものである。しかしながら、第１の実施の形態及び第２の実施の形態が、輝度成分と照明成分とから反射率を求める際の計算安定性を確保するものであるのに対し、第３の実施の形態及び第４の実施の形態は、カラー画像の色成分に対して後の計算を安定させるための変換を行うものである。効果は類似するが、画像が表すシーンによっては、何れか一方の構成を利用した方がよい場合がある。或いは、両方の構成を利用した方がよい場合もある。
【０１０６】
例えば、輝度成分に対しての暗部ノイズを主に考慮すればよい場合は、第１の実施の形態及び第２の実施の形態に重きをおけばよく、暗部ノイズで輝度以外のカラーノイズが気になる場合は、第３の実施の形態及び第４の実施の形態に重きをおけばよい。両方の構成を利用することで、どちらの対策にもなる。
【０１０７】
［画像処理装置のハードウェア構成］
本実施の形態における画像処理装置１０は、例えばＰＣにインストールされた画像処理ソフトウェアとしても実現され得るが、典型的には、画像読取り及び画像形成を行う画像処理装置１０として実現される。
【０１０８】
図２７は、このような画像処理装置１０のハードウェア構成例を示した図である。図示するように、画像処理装置１０は、ＣＰＵ（Central Processing Unit）３１と、ＲＡＭ（Random Access Memory）３２と、ＲＯＭ（Read Only Memory）３３と、ＨＤＤ（Hard Disk Drive）３４と、操作パネル３５と、画像読取部３６と、画像形成部３７と、通信インターフェース（以下、「通信Ｉ／Ｆ」と表記する）３８とを備える。
【０１０９】
ＣＰＵ３１は、ＲＯＭ３３等に記憶された各種プログラムをＲＡＭ３２にロードして実行することにより、後述する各機能を実現する。
【０１１０】
ＲＡＭ３２は、ＣＰＵ３１の作業用メモリ等として用いられるメモリである。
【０１１１】
ＲＯＭ３３は、ＣＰＵ３１が実行する各種プログラム等を記憶するメモリである。
【０１１２】
ＨＤＤ３４は、画像読取部３６が読み取った画像データや画像形成部３７における画像形成にて用いる画像データ等を記憶する例えば磁気ディスク装置である。
【０１１３】
操作パネル３５は、各種情報の表示やユーザからの操作入力の受付を行うタッチパネルである。ここで、操作パネル３５は、各種情報が表示されるディスプレイと、指やスタイラスペン等で指示された位置を検出する位置検出シートとからなる。
【０１１４】
画像読取部３６は、紙等の記録媒体に記録された画像を読み取る。ここで、画像読取部３６は、例えばスキャナであり、光源から原稿に照射した光に対する反射光をレンズで縮小してＣＣＤ（Charge Coupled Devices）で受光するＣＣＤ方式や、ＬＥＤ光源から原稿に順に照射した光に対する反射光をＣＩＳ（Contact Image Sensor）で受光するＣＩＳ方式のものを用いるとよい。
【０１１５】
画像形成部３７は、記録媒体に画像を形成する。ここで、画像形成部３７は、例えばプリンタであり、感光体に付着させたトナーを記録媒体に転写して像を形成する電子写真方式や、インクを記録媒体上に吐出して像を形成するインクジェット方式のものを用いるとよい。
【０１１６】
通信Ｉ／Ｆ３８は、ネットワークを介して他の装置との間で各種情報の送受信を行う。
【０１１７】
尚、本実施の形態を実現するプログラムは、通信手段により提供することはもちろん、ＣＤ−ＲＯＭ等の記録媒体に格納して提供することも可能である。
【符号の説明】
【０１１８】
１０…画像処理装置、１１，２２…色変換部、１２…照明成分推定部、１３…安定化反射率算出部、１４…輝度再現部、１５，２４…色逆変換部、１６，２６…暗部特性解析部、２１…カラー特性シフト部、２３…視認性再現部、２５…カラー特性逆シフト部、２７…頻度解析部
【図１】
2015226227.tif 000018
【図２】
2015226227.tif 000019
【図３】
2015226227.tif 000020
【図４】
2015226227.tif 000021
【図５】
2015226227.tif 000022
【図６】
2015226227.tif 000023
【図７】
2015226227.tif 000024
【図８】
2015226227.tif 000025
【図９】
2015226227.tif 000026
【図１０】
2015226227.tif 000027
【図１２】
2015226227.tif 000028
【図１３】
2015226227.tif 000029
【図１４】
2015226227.tif 000030
【図１５】
2015226227.tif 000031
【図１６】
2015226227.tif 000032
【図２０】
2015226227.tif 000033
【図２１】
2015226227.tif 000034
【図２２】
2015226227.tif 000035
【図２３】
2015226227.tif 000036
【図２４】
2015226227.tif 000037
【図２５】
2015226227.tif 000038
【図２６】
2015226227.tif 000039
【図２７】
2015226227.tif 000040
【図１１】
2015226227.tif 000041
【図１７】
2015226227.tif 000042
【図１８】
2015226227.tif 000043
【図１９】
2015226227.tif 000044