(57)【要約】  入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、もしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して色変換処理を行ない、順応ファクタを変換対象画素の色によって変化させる。
2008050761.tif 000002
【特許請求の範囲】
【請求項１】
入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法であって、
入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、もしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して前記色変換処理を行ない、
前記順応ファクタを変換対象画素の色によって変化させる画像処理方法。
【請求項２】
入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法であって、
出力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して前記色変換処理を行ない、
前記順応ファクタを変換対象画素の色によって変化させる画像処理方法。
【請求項３】
入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法であって、
入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタもしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタ、および、出力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して前記色変換処理を行ない、
前記それぞれの順応ファクタを変換対象画素の色によって変化させる画像処理方法。
【請求項４】
請求項１から請求項３のいずれかに記載の画像処理方法に於いて、
前記順応ファクタを変換対象画素の色相により変える画像処理方法。
【請求項５】
請求項１から請求項３のいずれかに記載の画像処理方法に於いて、
前記順応ファクタを変換対象画素の彩度に応じて変える画像処理方法。
【請求項６】
入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法であって、
入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、もしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して前記色変換処理を行ない、
前記入力画像に於ける色の分布を解析し、該解析した色の分布に基づいて、前記順応ファクタを変化させる画像処理方法。
【請求項７】
入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法であって、
出力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して前記色変換処理を行ない、
前記入力画像に於ける色の分布を解析し、該解析した色の分布に基づいて、前記順応ファクタを変化させる画像処理方法。
【請求項８】
入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法であって、
入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタもしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタ、および、出力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して前記色変換処理を行ない、
前記入力画像に於ける色の分布を解析し、該解析した色の分布に基づいて、前記それぞれの順応ファクタを変化させる画像処理方法。
【請求項９】
請求項６から請求項８のいずれかに記載の画像処理方法に於いて、
入力画像を複数領域に分割し、該分割領域毎の代表色を決定し、該代表色に基づいて前記順応ファクタを変化させる画像処理方法。
【請求項１０】
請求項６から請求項９のいずれかに記載の画像処理方法に於いて、
前記解析した色の分布に基づき、画像全体に適応する順応ファクタを決定する画像処理方法。
【請求項１１】
請求項６から請求項９のいずれかに記載の画像処理方法に於いて、
前記解析した色の分布に基づき、分割領域毎に適応する順応ファクタを変化させる画像処理方法。
【請求項１２】
請求項６から請求項１１のいずれかに記載の画像処理方法に於いて、
前記解析した色の分布に基づき、色相の分布を求め、該色相の分布に基づいて前記順応ファクタを変化させる画像処理方法。
【請求項１３】
請求項６から請求項１１のいずれかに記載の画像処理方法に於いて、
前記解析した色の分布に基づき、彩度の分布を求め、該彩度の分布に基づいて前記順応ファクタを変化させる画像処理方法。
【請求項１４】
請求項４または請求項１２のいずれかに記載の画像処理方法に於いて、
前記順応ファクタを前記入力画像取得時のシーンの照明の色温度若しくは前記入力画像観察時の視環境における照明の色温度又は前記出力画像観察時の視環境における照明の色温度により変える画像処理方法。
【請求項１５】
請求項４、１２、１４のいずれかに記載の画像処理方法に於いて、
前記順応ファクタを前記変換対象画素の彩度に応じて変える画像処理方法。
【請求項１６】
請求項１、３、６、８のいずれかに記載の画像処理方法に於いて、
入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境に関する情報を取得し、
前記照明条件、若しくは視環境に関する情報も考慮して前記順応ファクタを決定する画像処理方法。
【請求項１７】
請求項２、３、７、８のいずれかに記載の画像処理方法に於いて、
出力画像観察時の視環境に関する情報を取得し、
該視環境に関する情報も考慮して前記順応ファクタを決定する画像処理方法。
【請求項１８】
請求項１６または請求項１７のいずれかに記載の画像処理方法に於いて、
前記順応ファクタを変換対象画素の色相により変える画像処理方法。
【請求項１９】
請求項１８に記載の画像処理方法に於いて、
前記順応ファクタを前記入力画像取得時のシーンの照明の色温度もしくは前記入力画像観察時の視環境における照明の色温度又は前記出力画像観察時の視環境における照明の色温度により変える画像処理方法。
【請求項２０】
請求項１８または請求項１９に記載の画像処理方法に於いて、
前記順応ファクタを前記変換対象画素の彩度に応じて変える画像処理方法。
【請求項２１】
請求項１から３、５、６から１１、１３のいずれかに記載の画像処理方法に於いて、
入力画像取得時のシーンもしくは入力画像観察時の視環境に於ける第１の照明輝度を取得し、
出力画像観察時の視環境に於ける第２の照明輝度を取得し、
前記第１の照明輝度と前記第２の照明輝度の違いに基づいて、無彩色に対する順応ファクタと有彩色に対する順応ファクタを変える画像処理方法。
【請求項２２】
請求項２１に記載の画像処理方法に於いて、
前記第１の照明輝度が、前記第２の照明輝度より大きい場合に、前記無彩色に対する順応ファクタと前記有彩色に対する順応ファクタを変える画像処理方法。
【請求項２３】
請求項１８から２０のいずれかに記載の画像処理方法に於いて、
前記決定した順応ファクタを色相にわたって平均した値は、入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境、及び出力画像観察時の視環境によって変化する画像処理方法。
【請求項２４】
請求項２３記載の画像処理方法に於いて、
前記決定した順応ファクタを色相にわたって平均した値は、CIECAM02で定義される順応ファクタに対応する画像処理方法。
【請求項２５】
請求項１８から２０、２３、２４のいずれかに記載の画像処理方法に於いて、
前記順応ファクタDは、入力画像取得時のシーン、もしくは入力画像観察時もしくは出力画像観察時の視環境に於ける順応輝度La、および周囲の明るさで決まるパラメータFに関する関数α(La,F)、β(La,F)、および、変換対象画素もしくは前記解析した色の分布に基づき求めた色相hに関する関数γ(h)を用いて、
D=α(La,F)+β(La,F)×γ(h)
から算出する画像処理方法。
【請求項２６】
請求項１８から２０、２３、２４のいずれかに記載の画像処理方法に於いて、
前記順応ファクタDは、入力画像取得時のシーン、もしくは入力画像観察時もしくは出力画像観察時の視環境に於ける照明輝度Ld、および周囲の明るさで決まるパラメータFに関する関数α'(Ld,F)、β'(Ld,F)、および、変換対象画素もしくは前記解析した色の分布に基づき求めた色相hに関する関数γ(h)を用いて、
D=α'(Ld,F)+β'(Ld,F)×γ(h)
から算出する画像処理方法。
【請求項２７】
請求項２５に記載の画像処理方法に於いて、
常にα(La,F)=0である画像処理方法。
【請求項２８】
請求項２６に記載の画像処理方法に於いて、
常にα'(Ld,F)=0である画像処理方法。
【請求項２９】
請求項２５に記載の画像処理方法に於いて、
常にβ(La,F)=1である画像処理方法。
【請求項３０】
請求項２６に記載の画像処理方法に於いて、
常にβ'(Ld,F)=1である画像処理方法。
【請求項３１】
請求項２５に記載の画像処理方法に於いて、
Laの値が所定の閾値La_limit以下のとき、β(La,F)= 0とする画像処理方法。
【請求項３２】
請求項２６に記載の画像処理方法に於いて、
Ldの値が所定の閾値Ld_limit以下のとき、β'(Ld,F)=0とする画像処理方法。
【請求項３３】
請求項３１に記載の画像処理方法に於いて、
La_limitは100Cd/m2以下である画像処理方法。
【請求項３４】
請求項３１、３３のいずれかに記載の画像処理方法に於いて、
La_limitは前記出力画像の観察者特性により変更する画像処理方法。
【請求項３５】
請求項３２に記載の画像処理方法に於いて、
Ld_limitは前記出力画像の観察者特性により変更する画像処理方法。
【請求項３６】
請求項２５から３５のいずれかに記載の画像処理方法に於いて、
D=α(La,F)+β(La,F)×γ(h)
もしくは
D=α'(Ld,F)+β'(Ld,F)×γ(h)
より求めたDが１より大きいとき、D=１とし、Dが0より小さいとき、D=0とする画像処理方法。
【請求項３７】
請求項１４、１５、１９、２０、および２３から３６のいずれかに記載の画像処理方法に於いて、
前記入力画像取得時のシーン、もしくは入力画像観察時もしくは出力画像観察時の視環境に於ける照明の色温度が低い場合には、黄色から赤に対する順応ファクタを他の色より高く設定し、青に対する順応ファクタを他の色より低く設定する画像処理方法。
【請求項３８】
請求項１４、１５、１９、２０、および２３から３６のいずれかに記載の画像処理方法に於いて、
前記入力画像取得時のシーン、もしくは入力画像観察時もしくは出力画像観察時の視環境に於ける照明の色温度が高い場合には、黄色から赤に対する順応ファクタを他の色より低く設定し、青に対する順応ファクタを他の色より高く設定する画像処理方法。
【請求項３９】
請求項３１から３５のいずれかに記載の画像処理方法に於いて、
前記La_limit、もしくはLd_limitは、周囲の明るさによって変える画像処理方法。
【請求項４０】
入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法であって、
入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、入力画像観察時の視環境に順応している度合いを示す順応ファクタ、および、出力画像観察時の視環境に順応している度合いを示す順応ファクタの少なくとも１つを考慮して前記色変換処理を行い、
変換対象画素もしくは変換対象画素と周辺画素から決まる代表色によって前記順応ファクタを変化させるか否かを、視環境もしくは入力画像によって切り替える画像処理方法。
【請求項４１】
請求項４０に記載の画像処理方法に於いて、
前記切り替えは、順応ファクタの色相依存の有無を切り替える画像処理方法。
【請求項４２】
請求項４０から４１のいずれかに記載の画像処理方法に於いて、
前記変換対象画素もしくは変換対象画素と周辺画素から決まる代表色の彩度が所定の閾値以下の場合、前記順応ファクタの色相依存を無しにする画像処理方法。
【請求項４３】
請求項４０から４２のいずれかに記載の画像処理方法に於いて、
順応輝度が所定の閾値以下の場合、前記順応ファクタの色相依存を無しにする画像処理方法。
【請求項４４】
請求項４０から４３のいずれかに記載の画像処理方法に於いて、
入力画像取得時のシーンの照明若しくは入出力画像観察時の視環境における照明の種類によって、前記順応ファクタの色相依存を無しにする画像処理方法。
【請求項４５】
請求項４４に記載の画像処理方法に於いて、
入力画像取得時のシーンの照明若しくは入出力画像観察時の視環境における照明の色温度が所定範囲のとき、前記順応ファクタの色相依存を無しにする画像処理方法。
【請求項４６】
請求項４０から４５のいずれかに記載の画像処理方法に於いて、
前記順応ファクタの色相依存の有無の切り替え範囲を前記出力画像の観察者特性により変更する画像処理方法。
【請求項４７】
入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法であって、
入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタもしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタ、および、出力画像観察時の視環境に順応している度合いを示す順応ファクタの少なくとも１つを考慮して前記色変換処理を行ない、
前記順応ファクタを、前記入力画像の色に基づいて変化させる画像処理方法。
【請求項４８】
請求項１から請求項４７のいずれかに記載の画像処理方法に於いて、
色変換は、マトリックスにより変換する画像処理方法。
【請求項４９】
請求項1から請求項４８のいずれかに記載の面像処理方法に於いて、
入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、もしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタ、およぴ出力画像観察時の視環境に順応している度合いを示す順応ファクタは、見えのモードによって変える画像処理方法。
【請求項５０】
請求項４９に記載の面像処理方法に於いて、
見えのモードは少なくとも光源色、物体色のいずれかで切り替える画像処理方法。
【請求項５１】
請求項１から５０の何れかに記載の画像処理方法をコンピュータに実行させる画像処理プログラム。
【請求項５２】
請求項５１に記載の画像処理プログラムを搭載する画像処理装置。
【請求項５３】
カメラであって、
被写体を撮像する撮像部と、
表示部と、
制御部とを備え、
前記制御部は、請求項１から請求項５０のいずれかに記載の画像処理方法により前記撮像部により撮像され生成された画像データを色変換処理し前記表示部に表示するカメラ。
【請求項５４】
請求項５３に記載のカメラに於いて、
測光部をさらに備え、
前記制御部は、前記測光部により取得された測光データ、及び画像表示部の照明条件に基づいて、出力画像観察時の視環境を設定するカメラ。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、視環境の違いを考慮した色変換を行う画像処理方法、画像処理プログラム、画像処理装置、カメラに関する。
【背景技術】
【０００２】
照明の白色点や輝度、周囲の明るさと言ったいわゆる視環境が異なる条件では、色の見えが異なって見える。白色点の違いを考慮した色順応変換としてvon Kriesの色順応式が知られており、更に輝度や周囲の明るさも含めて視環境の違いを考慮した色の見えモデルとしてCIECAM02が知られている。視覚は必ずしも視環境の照明に完全順応している訳ではない為、白色点に順応している割合を示す順応ファクタD（0.0≦D≦1.0）を考慮した色順応変換が使用されている。【０００３】
順応ファクタDは視環境に従って決定される。例えばCIEが定めた色の見えモデルCIECAM97sやCIECAM02によれば、照明の輝度レベルや周囲輝度（周囲が明るいか、暗いか）に応じて決定された値を使用する（特許文献１）。また、特許文献２によれば、照明の相関色温度に基づき順応ファクタを変化させる方法も提案されている。【０００４】
【特許文献１】特許第３６３４６３３号明細書
【特許文献２】特開平１１−０５５６８８号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
上述の様に、従来の色の見えモデルによれば、視環境により順応ファクタDは決まり、実際に注目している色には依存しないとされている。しかし、視環境により順応ファクタDを変えても最適な見えが得られないという問題があった。【０００６】
また人間の視覚特性として、輝度が高いとコントラストが高く見えることが知られている(Stevens Effect)。CIECAM02色の見えモデルでは、この様な見えに関する効果を考慮した変換モデルとなっているが、処理が大変複雑であるという問題があった。一方、von Kriesの色順応変換の場合、CIECAM02より簡単な演算で済むが、白色点変換のみを考慮している為、輝度の違いによるコントラスト感を再現することができないという問題があった。
【課題を解決するための手段】
【０００７】
本発明の第１の態様によると、入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、もしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して色変換処理を行ない、順応ファクタを変換対象画素の色によって変化させる。
本発明の第２の態様によると、入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、出力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して色変換処理を行ない、順応ファクタを変換対象画素の色によって変化させる。
本発明の第３の態様によると、入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタもしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタ、および、出力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して色変換処理を行ない、それぞれの順応ファクタを変換対象画素の色によって変化させる。
本発明の第４の態様によると、第１から第３のいずれかの態様の画像処理方法に於いて、順応ファクタを変換対象画素の色相により変えるのが好ましい。
本発明の第５の態様によると、第１から第３のいずれかの態様の画像処理方法に於いて、順応ファクタを変換対象画素の彩度に応じて変えるのが好ましい。
本発明の第６の態様によると、入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、もしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して色変換処理を行ない、入力画像に於ける色の分布を解析し、該解析した色の分布に基づいて、順応ファクタを変化させる。
本発明の第７の態様によると、入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、出力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して色変換処理を行ない、入力画像に於ける色の分布を解析し、該解析した色の分布に基づいて、順応ファクタを変化させる。
本発明の第８の態様によると、入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタもしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタ、および、出力画像観察時の視環境に順応している度合いを示す順応ファクタを考慮して色変換処理を行ない、入力画像に於ける色の分布を解析し、該解析した色の分布に基づいて、それぞれの順応ファクタを変化させる。
本発明の第９の態様によると、第６から第８のいずれかの態様の画像処理方法に於いて、入力画像を複数領域に分割し、該分割領域毎の代表色を決定し、該代表色に基づいて順応ファクタを変化させるのが好ましい。
本発明の第１０の態様によると、第６から第９のいずれかの態様の画像処理方法に於いて、解析した色の分布に基づき、画像全体に適応する順応ファクタを決定するのが好ましい。
本発明の第１１の態様によると、第６から第９のいずれかの態様の画像処理方法に於いて、解析した色の分布に基づき、分割領域毎に適応する順応ファクタを変化させるのが好ましい。
本発明の第１２の態様によると、第６から第１１のいずれかの態様の画像処理方法に於いて、解析した色の分布に基づき、色相の分布を求め、該色相の分布に基づいて順応ファクタを変化させるのが好ましい。
本発明の第１３の態様によると、第６から第１１のいずれかの態様の画像処理方法に於いて、解析した色の分布に基づき、彩度の分布を求め、該彩度の分布に基づいて順応ファクタを変化させるのが好ましい。
本発明の第１４の態様によると、第４または第１２のいずれかの態様の画像処理方法に於いて、順応ファクタを入力画像取得時のシーンの照明の色温度若しくは入力画像観察時の視環境における照明の色温度又は出力画像観察時の視環境における照明の色温度により変えるのが好ましい。
本発明の第１５の態様によると、第４、１２、１４のいずれかの態様の画像処理方法に於いて、順応ファクタを変換対象画素の彩度に応じて変えるのが好ましい。
本発明の第１６の態様によると、第１、３、６、８のいずれかの態様の画像処理方法に於いて、入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境に関する情報を取得し、照明条件、若しくは視環境に関する情報も考慮して順応ファクタを決定するのが好ましい。
本発明の第１７の態様によると、第２、３、７、８のいずれかの態様の画像処理方法に於いて、出力画像観察時の視環境に関する情報を取得し、該視環境に関する情報も考慮して順応ファクタを決定するのが好ましい。
本発明の第１８の態様によると、第１６または第１７のいずれかの態様の画像処理方法に於いて、順応ファクタを変換対象画素の色相により変えるのが好ましい。
本発明の第１９の態様によると、第１８の態様の画像処理方法に於いて、順応ファクタを入力画像取得時のシーンの照明の色温度もしくは入力画像観察時の視環境における照明の色温度又は出力画像観察時の視環境における照明の色温度により変えるのが好ましい。
本発明の第２０の態様によると、第１８または第１９の態様の画像処理方法に於いて、順応ファクタを変換対象画素の彩度に応じて変えるのが好ましい。
本発明の第２１の態様によると、第１から３、５、６から１１、１３のいずれかの態様の画像処理方法に於いて、入力画像取得時のシーンもしくは入力画像観察時の視環境に於ける第１の照明輝度を取得し、出力画像観察時の視環境に於ける第２の照明輝度を取得し、第１の照明輝度と第２の照明輝度の違いに基づいて、無彩色に対する順応ファクタと有彩色に対する順応ファクタを変えるのが好ましい。
本発明の第２２の態様によると、第２１の態様の画像処理方法に於いて、第１の照明輝度が、第２の照明輝度より大きい場合に、無彩色に対する順応ファクタと有彩色に対する順応ファクタを変えるのが好ましい。
本発明の第２３の態様によると、第１８から２０のいずれかの態様の画像処理方法に於いて、決定した順応ファクタを色相にわたって平均した値は、入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境、及び出力画像観察時の視環境によって変化するのが好ましい。
本発明の第２４の態様によると、第２３記載の画像処理方法に於いて、決定した順応ファクタを色相にわたって平均した値は、CIECAM02で定義される順応ファクタに対応するのが好ましい。
本発明の第２５の態様によると、第１８から２０、２３、２４のいずれかの態様の画像処理方法に於いて、順応ファクタDは、入力画像取得時のシーン、もしくは入力画像観察時もしくは出力画像観察時の視環境に於ける順応輝度La、および周囲の明るさで決まるパラメータFに関する関数α(La,F)、β(La,F)、および、変換対象画素もしくは解析した色の分布に基づき求めた色相hに関する関数γ(h)を用いて、D=α(La,F)+β(La,F)×γ(h)から算出するのが好ましい。
本発明の第２６の態様によると、第１８から２０、２３、２４のいずれかの態様の画像処理方法に於いて、順応ファクタDは、入力画像取得時のシーン、もしくは入力画像観察時もしくは出力画像観察時の視環境に於ける照明輝度Ld、および周囲の明るさで決まるパラメータFに関する関数α'(Ld,F)、β'(Ld,F)、および、変換対象画素もしくは解析した色の分布に基づき求めた色相hに関する関数γ(h)を用いて、D=α'(Ld,F)+β'(Ld,F)×γ(h)から算出するのが好ましい。
本発明の第２７の態様によると、第２５の態様の画像処理方法に於いて、常にα(La,F)=0であるのが好ましい。
本発明の第２８の態様によると、第２６の態様の画像処理方法に於いて、常にα'(Ld,F)=0であるのが好ましい。
本発明の第２９の態様によると、第２５の態様の画像処理方法に於いて、常にβ(La,F)=1であるのが好ましい。
本発明の第３０の態様によると、第２６の態様の画像処理方法に於いて、常にβ'(Ld,F)=1であるのが好ましい。
本発明の第３１の態様によると、第２５の態様の画像処理方法に於いて、Laの値が所定の閾値La_limit以下のとき、β(La,F)= 0とするのが好ましい。
本発明の第３２の態様によると、第２６の態様の画像処理方法に於いて、Ldの値が所定の閾値Ld_limit以下のとき、β'(Ld,F)=0とするのが好ましい。
本発明の第３３の態様によると、第３１の態様の画像処理方法に於いて、La_limitは100Cd/m2以下であるのが好ましい。
本発明の第３４の態様によると、第３１、３３のいずれかの態様の画像処理方法に於いて、La_limitは出力画像の観察者特性により変更するのが好ましい。
本発明の第３５の態様によると、第３２の態様の画像処理方法に於いて、Ld_limitは出力画像の観察者特性により変更するのが好ましい。
本発明の第３６の態様によると、第２５から３５のいずれかの態様の画像処理方法に於いて、D=α(La,F)+β(La,F)×γ(h)、もしくは、D=α'(Ld,F)+β'(Ld,F)×γ(h)より求めたDが１より大きいとき、D=１とし、Dが0より小さいとき、D=0とするのが好ましい。
本発明の第３７の態様によると、第１４、１５、１９、２０、および２３から３６のいずれかの態様の画像処理方法に於いて、入力画像取得時のシーン、もしくは入力画像観察時もしくは出力画像観察時の視環境に於ける照明の色温度が低い場合には、黄色から赤に対する順応ファクタを他の色より高く設定し、青に対する順応ファクタを他の色より低く設定するのが好ましい。
本発明の第３８の態様によると、第１４、１５、１９、２０、および２３から３６のいずれかの態様の画像処理方法に於いて、入力画像取得時のシーン、もしくは入力画像観察時もしくは出力画像観察時の視環境に於ける照明の色温度が高い場合には、黄色から赤に対する順応ファクタを他の色より低く設定し、青に対する順応ファクタを他の色より高く設定するのが好ましい。
本発明の第３９の態様によると、第３１から３５のいずれかの態様の画像処理方法に於いて、La_limit、もしくはLd_limitは、周囲の明るさによって変えるのが好ましい。
本発明の第４０の態様によると、入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、入力画像観察時の視環境に順応している度合いを示す順応ファクタ、および、出力画像観察時の視環境に順応している度合いを示す順応ファクタの少なくとも１つを考慮して色変換処理を行い、変換対象画素もしくは変換対象画素と周辺画素から決まる代表色によって順応ファクタを変化させるか否かを、視環境もしくは入力画像によって切り替える。
本発明の第４１の態様によると、第４０の態様の画像処理方法に於いて、切り替えは、順応ファクタの色相依存の有無を切り替えるのが好ましい。
本発明の第４２の態様によると、第４０から４１のいずれかの態様の画像処理方法に於いて、変換対象画素もしくは変換対象画素と周辺画素から決まる代表色の彩度が所定の閾値以下の場合、順応ファクタの色相依存を無しにするのが好ましい。
本発明の第４３の態様によると、第４０から４２のいずれかの態様の画像処理方法に於いて、順応輝度が所定の閾値以下の場合、順応ファクタの色相依存を無しにするのが好ましい。
本発明の第４４の態様によると、第４０から４３のいずれかの態様の画像処理方法に於いて、入力画像取得時のシーンの照明若しくは入出力画像観察時の視環境における照明の種類によって、順応ファクタの色相依存を無しにするのが好ましい。
本発明の第４５の態様によると、第４４の態様の画像処理方法に於いて、入力画像取得時のシーンの照明若しくは入出力画像観察時の視環境における照明の色温度が所定範囲のとき、順応ファクタの色相依存を無しにするのが好ましい。
本発明の第４６の態様によると、第４０から４５のいずれかの態様の画像処理方法に於いて、順応ファクタの色相依存の有無の切り替え範囲を出力画像の観察者特性により変更するのが好ましい。
本発明の第４７の態様によると、入力画像データに対して色順応を考慮した色変換処理を行う画像処理方法は、入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタもしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタ、および、出力画像観察時の視環境に順応している度合いを示す順応ファクタの少なくとも１つを考慮して色変換処理を行ない、順応ファクタを、入力画像の色に基づいて変化させる。
本発明の第４８の態様によると、第１から第４７のいずれかの態様の画像処理方法に於いて、色変換は、マトリックスにより変換するのが好ましい。
本発明の第４９の態様によると、第1から第４８のいずれかの態様の面像処理方法に於いて、入力画像取得時のシーンの照明条件に順応している度合いを示す順応ファクタ、もしくは入力画像観察時の視環境に順応している度合いを示す順応ファクタ、およぴ出力画像観察時の視環境に順応している度合いを示す順応ファクタは、見えのモードによって変えるのが好ましい。
本発明の第５０の態様によると、第４９の態様の面像処理方法に於いて、見えのモードは少なくとも光源色、物体色のいずれかで切り替えるのが好ましい。
本発明の第５１の態様によると、画像処理プログラムは、第１から５０の何れかの態様の画像処理方法をコンピュータに実行させる。
本発明の第５２の態様によると、画像処理装置は、第５１の態様の画像処理プログラムを搭載する。
本発明の第５３の態様によると、カメラは、被写体を撮像する撮像部と、表示部と、制御部とを備え、制御部は、第１から第５０のいずれかの態様の画像処理方法により撮像部により撮像され生成された画像データを色変換処理し表示部に表示する。
本発明の第５４の態様によると、第５３の態様のカメラに於いて、測光部をさらに備え、制御部は、測光部により取得された測光データ、及び画像表示部の照明条件に基づいて、出力画像観察時の視環境を設定するのが好ましい。
【発明の効果】
【０００８】
視環境や照明条件が異なる場合に、精度良く近い色の見えを再現できる。
【図面の簡単な説明】
【０００９】
【図１】本発明の第１の実施の形態である画像処理装置を示す図である。
【図２】最適な順応ファクタDh*と色相角h*の関係を示す図である。
【図３】第１の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図４】第２の実施の形態の、顕微鏡１００とパーソナルコンピュータ１とデジタルカメラ２とからなる顕微鏡観察システムを示す図である。
【図５】第２の実施の形態におけるパーソナルコンピュータ１が実行する画像処理プログラムのフローチャートを示す図である。
【図６】第３の実施の形態において、色変換マトリックスを求める処理のフローチャートを示す図である。
【図７】第３の実施の形態のデジタルカメラ２０の構成を示す図である。
【図８】デジタルカメラ２０の制御部２１が行う処理のフローチャートを示す図である。
【図９】第４の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図１０】関数α(C*)を示すグラフである。
【図１１】第５の実施の形態におけるパーソナルコンピュータ１が実行する画像処理プログラムのフローチャートを示す図である。
【図１２】関数D(C*)を示すグラフである。
【図１３】関数D'(C*)を示すグラフである。
【図１４】第６の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図１５】第７の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図１６】第８の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図１７】第８の実施の形態の変形例のパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図１８】第９の実施の形態のデジタルカメラ３０の構成を示す図である。
【図１９】輝度比対順応ファクタDachrの非線形関数の例を示す図である。
【図２０】色温度毎における最適な順応ファクタDh*と色相角h*の関係を示す図である。
【図２１】彩度毎における最適な順応ファクタDh*と色相角h*の関係を示す図である。
【図２２】視環境毎における最適な順応ファクタDh*と色相角h*の関係を示す図である。
【図２３】第１３の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図２４】第１４の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図２５】第１４の実施の形態の変形例１におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図２６】第１４の実施の形態の変形例２におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図２７】第１４の実施の形態の変形例３におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図２８】第１４の実施の形態の変形例４におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図２９】第１４の実施の形態のその他の変形例におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図３０】周囲の明るさに応じた複数のLa_limitを示す図である。
【図３１】第１８の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【図３２】第１９の実施の形態における順応輝度毎における最適な順応ファクタと色相角の関係を示す図である。
【図３３】第１９の実施の形態の変形例における周囲の明るさに応じた複数の順応輝度閾値を示す図である。
【図３４】第１８〜第２０の実施の形態の変形例５におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。
【発明を実施するための最良の形態】
【００１０】
−第１の実施の形態−
第１の実施の形態では、不完全順応の度合いを表す順応ファクタを、変換対象色の色相に依存して変えることにより、色の見え再現をより精度よく行うことが特徴である。【００１１】
図１は、本発明の一実施の形態である画像処理装置を示す図である。画像処理装置は、パーソナルコンピュータ１で実現される。パーソナルコンピュータ１は、デジタルカメラ２、ＣＤ−ＲＯＭなどの記録媒体３、他のコンピュータ４などと接続され、各種の画像データの提供を受ける。パーソナルコンピュータ１は、提供された画像データに対して、以下に説明する画像処理を行う。コンピュータ４は、インターネットやその他の電気通信回線５を経由して接続される。【００１２】
パーソナルコンピュータ１が画像処理のために実行するプログラムは、図１の構成と同様に、ＣＤ−ＲＯＭなどの記録媒体や、インターネットやその他の電気通信回線を経由した他のコンピュータから提供され、パーソナルコンピュータ１内にインストールされる。パーソナルコンピュータ１は、ＣＰＵ（不図示）、メモリ（不図示）、周辺回路（不図示）、モニタ６などから構成され、ＣＰＵがインストールされたプログラムを実行する。【００１３】
プログラムがインターネットやその他の電気通信回線を経由して提供される場合は、プログラムは、電気通信回線、すなわち、伝送媒体を搬送する搬送波上の信号に変換して送信される。このように、プログラムは、記録媒体や搬送波などの種々の形態のコンピュータ読み込み可能なコンピュータプログラム製品として供給される。【００１４】
本実施の形態のパーソナルコンピュータ１は、撮影シーンの照明条件で撮影された画像データを、人間の色順応特性を考慮して、撮影シーンの照明条件とは異なる画像観察時の視環境に応じた色の見えを再現した画像データに変換する。【００１５】
（色の見え実験結果）
まず、色の見え実験について説明する。色の見えに熟練した被験者に、色温度の異なる２つの視環境に於いて、色の見えが近く見える色を選択してもらう実験である。色温度の異なる視環境にて、順応ファクタを変えて色変換した複数の色の中から、見えが最も近くなる色を選んでもらう色の見え実験を行うと、色相に依存して最適な順応ファクタDが変化する結果が得られる。【００１６】
色温度T1の視環境下に参照色A(h*)を表示する。参照色A(h*)は、色相h*により表される色である。なお、色相h*は色相角で表される。色温度T2の視環境に、参照色A(h*)から色温度T1,T2の違いを考慮してCIECAM02にて色変換した結果を表示する。このとき、順応ファクタDは照明輝度や周囲の明るさ等の視環境から予測される値ではなく、D=0.0〜1.0の間で変えた複数の色B(h*)(D)を選択肢として表示する。被験者は、色温度T1の視環境下でA(h*)を見たときの色の見えと、色温度T2の視環境下でB(h*)(D)（D=0.0 …,1.0）を見たときの色の見えを比較し、両者の見えが最も近い色B(h*)(Dh*)を選択する。この様な実験を、参照色A(h*)の輝度や彩度は変えずに色相角h*のみ変えた複数の色について行う。【００１７】
更に複数被験者に複数回ずつ実験を行ってもらい、結果を集計すると図２の様な結果が得られた。図２は最適な順応ファクタDh*と色相角h*の関係を示しており、横軸は色相角h*（単位はdeg）、縦軸は夫々の実験で最も見えが近いと選ばれた最適順応ファクタDh*の被験者平均値である。本実施の形態では、h*=30deg（赤）、h*=90deg（黄）、h*=150deg（緑）、h*=210deg（シアン）、h*=270deg（青）、h*=330deg（マゼンタ）の６色についての実験を行った。図２の結果から色相に依存して最適な順応ファクタが変化しており、h*=90deg（黄）でDh*を大きく、h*=270deg（青）でDh*を小さくした方がより色の見えが近くなることが分かる。【００１８】
（変換処理）
次に、色変換処理について説明する。図３は、第１の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。予め上述の様な実験を行って、最適な順応ファクタと色相の関係Dh*(h*)をテーブルや近似式の形で画像処理を行うプログラム内に保有する。あるいは、プログラムが参照するデータとしてパーソナルコンピュータ１の記憶装置やメモリ等に保存しておく。【００１９】
まずステップＳ２１にて、撮影画像を取得する。撮影画像は、複数の画素から構成され、カメラの分光感度に固有の色空間RGBで記述されている。【００２０】
次にステップＳ２２に進む。ステップＳ２２では、撮影シーンの照明に関する情報を取得する。即ち、撮影シーンの照明白色点(XW,YW,ZW)と輝度YW0と周囲輝度YSW0である。例えば分光放射輝度計などを用いて白色点、及び輝度を測定した結果をユーザ（観察者）がキーボードから入力した値を読み込む。もしくはデジタルカメラに備わっているオートホワイトバランス機能を用いて、撮影時の照明白色点(XW,YW,ZW)を算出し、露出を自動調整する為に用いる測光機能の結果を解析して撮影時の照明輝度YW0を算出してもよい。周囲輝度YSW0は、測光範囲を広げて測光した結果を使用して算出してもよいし、撮影者が予め入力し画像に付加された情報やキーボードから入力した値を使用してもよい。【００２１】
ステップＳ２３では、画像観察時の視環境に関する情報を取得する。標準視環境にてモニタ観察することを前提とすれば、sRGB標準視環境（白色点=D65、輝度Ydevice0'=80Cd/m2、周囲輝度YSW0'=4.1Cd/m2）である。sRGB標準視環境と異なる視環境で観察する場合や画像をプリントした印刷物で観察する場合には、パーソナルコンピュータ１は、ユーザ（観察者）がキーボード等で入力した画像観察に使用するデバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'、および周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'を読み込み、デバイス白色点と周囲照明白色点の混合として画像観察視環境に於ける照明白色点(XW',YW',ZW')、輝度YW0'を決定する。【００２２】
具体的には、画像観察に使用するデバイス光源と周囲照明光源の影響の比率Rmix'を用いて、次のように算出する。
【数１】
2008050761.tif 000003
【数２】
2008050761.tif 000004
【数３】
2008050761.tif 000005
【数４】
2008050761.tif 000006
【数５】
2008050761.tif 000007
尚、MCAT02はXYZから錐体応答LMSへの変換マトリックスである。
【数６】
2008050761.tif 000008
【００２３】
次にステップＳ２４に進む。ステップＳ２４〜２８の処理は、全画素について処理が終了するまで繰り返す。カメラRGBはカメラの分光感度で決まる固有の色空間RGBで記述されているので、予め決まっているカメラRGB→XYZマトリックスMRGB→XYZを使用して、CIE 1931 XYZ色空間へ変換する。
【数７】
2008050761.tif 000009
【００２４】
次にステップＳ２５に進み、変換対象画素の色相角h*を算出する。色相角h*の算出は、CIE XYZからCIE LABに変換し、a*,b*の座標でなす角を0〜360°で算出した値を使用する。【００２５】
ステップＳ２６では、予め色の見え実験結果から決定しておいた最適な順応ファクタと色相の関係Dh*(h*)より、ステップＳ２５で算出したh*の時の最適な色の見えを与える順応ファクタDh*(h*)を算出する。順応ファクタは周囲の明るさや順応輝度によっても変化する為、撮影時の視環境と画像観察時の視環境では最適な順応ファクタが異なる。従って画像観察時の順応ファクタDh*'(h*)も算出する。【００２６】
ステップＳ２７では、ステップＳ２２、Ｓ２３で読み込んだ視環境の違いと、ステップＳ２６で算出したDh*(h*) 、Dh*'(h*)を使用して、撮影シーンの照明条件下に於ける色XYZと色の見えが近くなる観察照明条件下の色X'Y'Z'を算出する。ここで使用する色の見えモデルとしては、CIECAM02やvon Kries色順応式など視環境の違いを考慮した色順応変換式であればよい。CIECAM02を使用する場合は、XYZからJChへの変換では、定義式から算出される順応ファクタDの代わりに、ステップＳ２６で算出したDh*(h*)を使用し、JChからX'Y'Z'への変換では、Dh*'(h*)を使用すればよい。Von Kriesの色順応式を使用する場合には、次式にて色変換を行う。
【数８】
2008050761.tif 000010
但し、
【数９】
2008050761.tif 000011
【数１０】
2008050761.tif 000012
【数１１】
2008050761.tif 000013
【数１２】
2008050761.tif 000014
【数１３】
2008050761.tif 000015
【００２７】
ステップＳ２８に進み、X'Y'Z'から出力画像の色空間R'G'B'に変換する。MXYZ→sRGBはXYZから出力画像の色空間RGBに変換するマトリックスであり、出力画像の色空間をsRGBとすれば、規格で決まっているマトリックスを使用すればよい。
【数１４】
2008050761.tif 000016
【００２８】
ステップＳ２９では、全画素処理が終了していなければ、上記ステップＳ２４に戻し、処理を繰り返す。全画素終了したら、ステップＳ３０に進み画像を出力して終了する。【００２９】
以上のようにして、画像撮影時のシーンの照明条件と、撮影画像を観察する照明条件が異なる場合に、精度良く近い色の見えを再現できる。上述の様に、本実施の形態の画像処理装置では、変換対象画素の色すなわち色相により最適な順応ファクタを使用して変換している為、正確な色の見え予測をすることができる。【００３０】
尚、本実施の形態は、CIECAM02の様な色の見えモデルや白色点変換のみを考慮したvon Kriesの色順応変換など、色順応を考慮した多様な色変換方法に適用可能である。【００３１】
尚、ステップＳ２３では、画像観察視環境に於ける照明白色点として、CIE 1931 XYZ表色系の値(XW',YW',ZW')を算出している。しかし式（４）を省略して、錐体応答の空間に於ける値(LW',MW',SW')を、画像観察視環境に於ける照明白色点として取得してもよい。この場合、ステップＳ２７の式（１３）も省略することができ、ステップＳ２３で取得した(LW',MW',SW')をそのまま式（１１）に使用することができる。従って処理を軽くすることが出来る。【００３２】
上記第１の実施の形態では、モニタ６に画像表示する場合を例に説明したが、出力画像をプロジェクタやプリントで表示してもよい。プリントで表示する場合の入力データの例はステップＳ２３の説明に有る通りである。プロジェクタに表示する場合、出力画像観察時の視環境は、ユーザがキーボード等で入力した値を使用する。デバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'は、スペックに書かれているデータを用いるか、ユーザが実測値を入力する。周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'は、暗室すなわち暗い部屋で観察する場合はYSW0'=0、(XSW',YSW',ZSW')= (Xdevice',Ydevice',Zdevice')を入力する。暗室でない場合は、ユーザが周囲視環境を実測し、実測値を入力する。【００３３】
上記第１の実施の形態では、色相方向について行った実験結果を利用する例を説明した。しかし、彩度や明度方向にも色相方向と同様の実験を行い、順応ファクタ算出に用いてもよい。この場合、CIE LABにおいて、彩度はC*=(a*2+b*2)1/2、明度はL*の値を使用する。【００３４】
また、色相、明度、彩度としてCIE LABのh*, L*, C*を使用する例で説明したが、CIE LAB以外の色空間に於けるパラメータを使用してもよい。例えば、CIECAM02の相対的な明度、彩度、色相のパラメータであるJ, C, h等を使用してもよい。【００３５】
−第２の実施の形態−
第２の実施の形態では、CIECAM02の様な複雑な計算を用いずに、輝度差の大きい視環境間の色の見え再現、特に輝度差によってコントラストが強く見える見え特性を簡単に実現できることが特徴である。【００３６】
第２の実施の形態では、顕微鏡像（顕微鏡によって拡大された像のことを言う）をデジタルカメラ２で撮影した撮影画像が、モニタ表示やプリント物として観察した際、実際の肉眼で見た顕微鏡観察像と近い見えになるように色変換する方法について説明する。【００３７】
図４は、本実施の形態の、顕微鏡１００とパーソナルコンピュータ１とデジタルカメラ２とからなる顕微鏡観察システムを示す図である。パーソナルコンピュータ１とデジタルカメラ２は、第１の実施の形態と同様である。パーソナルコンピュータ１とデジタルカメラ２とは、ケーブル７により接続される。【００３８】
顕微鏡１００では、照明用の光源１１により標本（被写体）１２が照明され、対物レンズ１３と接眼レンズ１４により拡大された標本像が観察できる。図４に示すように、顕微鏡１００の上部にはデジタルカメラ２が取り付けられている。パーソナルコンピュータ１は、デジタルカメラ２で撮影された画像を取り込み、取り込んだ画像をモニタ６にリアルタイムで表示する。【００３９】
顕微鏡１００は、視野切り替えスイッチ１５を有し、接眼レンズ１４側とデジタルカメラ２側とに視野を切り替える。これにより、接眼レンズ１４を通して目視で標本１２を直接観察したり、デジタルカメラ２で撮影した顕微鏡像の画像がパーソナルコンピュータ１のモニタ６に表示され、モニタ６上で標本１２を観察したりすることができる。視野切り替えスイッチ１５は、完全に接眼レンズ１４側とデジタルカメラ２側に切り替えるだけでなく、光を半分に分けることも可能である。これにより、接眼レンズ１４側からとモニタ６側からと同時に観察することも可能となる。【００４０】
本実施の形態では、顕微鏡１００の接眼レンズ１４を介して直接観察した時の顕微鏡像と、デジタルカメラ２で撮影しモニタ６で表示した顕微鏡像の画像の色の見えが近くなるようにパーソナルコンピュータ１内で所定の色変換処理が行われる。【００４１】
顕微鏡１００とモニタ６の表示画像の色の見えに影響を及ぼす要因として、主に、光源１１とモニタ６の色温度の違い、光源１１とモニタ６の照明輝度の違い、接眼レンズ１４を通して顕微鏡像を観察する時とモニタ６観察時の周囲照明の違いなどがある。モニタ６の色温度は、標準視環境であればD65と決まっているが、顕微鏡照明の色温度は、光源１１に用いているランプの特性で決まり、3000K〜5000K程度とモニタより低い色温度の場合が多い。従って、照明用の光源１１とモニタ６の色温度の違いを考慮してモニタ６に表示しないと、接眼レンズ１４により肉眼で直接見る顕微鏡観察像と色が異なって見えることになる。【００４２】
また、モニタ６の場合、輝度もsRGB標準視環境では80Cd/m2と決まっており、よく使われるLCDの場合最大でも200Cd/ m2程度までしか設定できない。実際には、最大輝度に設定していても経年劣化で輝度が落ちてくるので200Cd/ m2以下で使用することが多い。一方、顕微鏡を接眼レンズ１４を介して観察している場合、輝度は300Cd/ m2以上で観察していることが多く、場合によっては1000Cd/ m2以上とモニタ観察時より一桁近く高い輝度で観察している条件もある。従って顕微鏡観察像はモニタ観察時より輝度が１．５倍以上高く、コントラストが高く感じるという違いがある。【００４３】
また、周囲照明による影響も異なっている。接眼レンズ１４による顕微鏡観察時は、顕微鏡像しか見えない為、周囲の影響（周囲が明るいのか暗いのか）は余り受けていない。しかし、モニタ６による観察時は、通常明るい実験室で使用している場合が多い。このような理由からも、接眼レンズ１４で直接観察している場合に比べ、モニタ表示画像はフレアの影響等によりコントラストが低く見える。また、モニタ６の色温度と周囲照明の色温度が違う場合には、視覚は周囲照明の色温度にも影響を受けるため、周囲照明の色温度の考慮も必要になる。【００４４】
本実施の形態では、この様な視環境の違いを考慮して、接眼レンズ１４を通した直接観察の場合とモニタ６を介した間接観察の場合とで、標本１２の顕微鏡観察像の色の見えが近くなる様にパーソナルコンピュータ１内部で以下に説明する画像処理を行う。なお、プリント物で観察する場合は、パーソナルコンピュータ１に接続された不図示のプリンターでプリントアウトした画像が観察される。【００４５】
図５は、第２の実施の形態におけるパーソナルコンピュータ１が実行する画像処理プログラムのフローチャートを示す図である。【００４６】
まず、第１の実施の形態（図３）と同様に、ステップＳ２１にて、撮影画像を取得する。撮影画像はカメラの分光感度に固有の色空間RGBで記述されている。【００４７】
ステップＳ４１では、撮影シーンの照明に関する情報を取得する。即ち、撮影シーンの照明白色点と輝度であり、本実施の形態では、光源１１の顕微鏡照明の白色点(XW,YW,ZW)、及び輝度YW0に関する情報である。例えば分光放射輝度計などを用いて顕微鏡照明の白色点、及び輝度を測定する。もしくはデジタルカメラに備わっているオートホワイトバランス機能を用いて、撮影時の照明白色点(XW,YW,ZW)を推定し、露出を自動調整する為に用いる測光機能の結果を解析して撮影時の照明輝度YW0を推定してもよい。顕微鏡観察の場合には周囲照明の影響が小さいので、無視してもよい。【００４８】
ステップＳ４２では、画像観察時の視環境に関する情報を取得する。標準視環境にてモニタ観察することを前提とすれば、sRGB標準視環境（白色点=D65、輝度Ydevice0'=80Cd/m2、周囲輝度YSW0'=4.1Cd/m2）である。sRGB標準視環境と異なる視環境で観察する場合や画像をプリントした印刷物で観察する場合には、パーソナルコンピュータ１は、ユーザ（観察者）がキーボード等で入力した画像観察に使用するデバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'、および周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'を読み込み、デバイス白色点と周囲照明白色点の混合として画像観察視環境に於ける照明白色点(XW',YW',ZW')、輝度YW0'を決定する。【００４９】
具体的には、画像観察に使用するデバイス光源と周囲照明光源の影響の比率Rmix'を用いて、次のように算出する。
【数１５】
2008050761.tif 000017
【数１６】
2008050761.tif 000018
【数１７】
2008050761.tif 000019
【数１８】
2008050761.tif 000020
【数１９】
2008050761.tif 000021
【００５０】
尚、MCAT02は第１の実施の形態で用いたXYZから錐体応答LMSへの変換マトリックスと同じである。【００５１】
ステップＳ４３では、ステップＳ４１で取得した撮影シーンの条件から、無彩色に対する順応ファクタDachr、有彩色に対する順応ファクタDchrを算出する。DchrはCIECAM02で定められた次式等で算出される値を使用する。
【数２０】
2008050761.tif 000022
【００５２】
ここで、順応輝度LAは、通常LA=YW0/5を用いる。またFはデバイス輝度と周囲輝度の比に応じて決まるパラメータである。一方無彩色（白やグレー）に対する順応ファクタDachr=1.0と設定する。輝度が高い場合、照明の色が変わっても白は白く感じる（完全順応に近くなる）人が多いという視覚特性を考慮したもので、結果として生物顕微鏡観察像で言えば標本の有る部分（染色されている部分＝有彩色）と無い部分（染色のない部分＝無彩色）でコントラストが付いて見える。【００５３】
次にステップＳ４４では、ステップＳ４２で取得した画像観察時の視環境から、画像観察時の視環境に於ける無彩色に対する順応ファクタDachr'、及び有彩色に対する順応ファクタDchr'を算出する。有彩色に対する順応ファクタDchr'は、ステップＳ４３同様、画像観察時の順応輝度LA'と周囲の明るさで決まるパラメータF'を用いて次式で算出できる。
【数２１】
2008050761.tif 000023
【数２２】
2008050761.tif 000024
【００５４】
また無彩色に対する順応ファクタDachr'=1.0と設定する。【００５５】
ステップＳ４５にて、ステップＳ４１〜Ｓ４４にて取得・算出したデータを用いて白色点変換マトリックス（無彩色に対するマトリックスMachr、有彩色に対するマトリックスMchr）を算出する。白色点変換は人間の錐体応答に近い色空間であるLMS空間で行う。
【数２３】
2008050761.tif 000025
【数２４】
2008050761.tif 000026
【００５６】
このとき、Mchr、およびMachrを次式で算出する。
【数２５】
2008050761.tif 000027
【数２６】
2008050761.tif 000028
但し、
【数２７】
2008050761.tif 000029
【数２８】
2008050761.tif 000030
【数２９】
2008050761.tif 000031
【数３０】
2008050761.tif 000032
【００５７】
次にステップＳ４６に進む。ステップＳ４６では、ステップＳ２１で取得した画像の各画素のカメラRGBを取得する。ステップＳ４７に進み、カメラRGBはカメラの分光感度で決まる固有の色空間RGBで記述されているので、予め決まっているカメラRGB→XYZマトリックスMRGB→XYZを使用して、機器に依存しない色空間であるCIE 1931 XYZ色空間へ変換する。
【数３１】
2008050761.tif 000033
【００５８】
ステップＳ４８に進み、上述と同じMCAT02を使用してLMS空間へ変換する。
【数３２】
2008050761.tif 000034
【００５９】
ステップＳ４９に進み、対象画素が有彩色か否（無彩色）かの判断を行う。有彩色か否（無彩色）かの判断は、CIE XYZからCIE LABに変換し、彩度C*=(a*2+b*2)1/2と閾値C*0を用いて、C*≧C*0の場合有彩色、C*<C*0の時無彩色と判断する。【００６０】
ステップＳ４９にて肯定された場合は、ステップＳ５０に進み、ステップＳ４５で作成したMchrを使用して白色点変換を行う。
【数３３】
2008050761.tif 000035
【００６１】
ステップＳ４９にて否定された場合は、ステップＳ５１に進み、ステップＳ４５で作成したMachrを使用して白色点変換を行う。
【数３４】
2008050761.tif 000036
【００６２】
ステップＳ５２では、ステップＳ４８の逆変換、すなわちLMS色空間からXYZ色空間に色変換する。
【数３５】
2008050761.tif 000037
【００６３】
ステップＳ５３では、X'Y'Z'から出力画像の色空間R'G'B'に変換する。MXYZ→sRGBはXYZから出力画像の色空間RGBに変換するマトリックスであり、出力画像の色空間をsRGBとすれば、規格で決まっているマトリックスを使用すればよい。
【数３６】
2008050761.tif 000038
【００６４】
ステップＳ５４では、全画素処理が終了していなければ、上記ステップＳ４６に戻し、処理を繰り返す。全画素終了したら、ステップＳ５５に進み画像を出力して終了する。【００６５】
以上の様にして、顕微鏡観察時のシーンの照明輝度と撮影画像観察時の照明輝度が、1.5倍以上と大きく異なる場合でも、彩度に応じて順応ファクタを変えることで、簡単に色の見えを再現することができる。更に顕微鏡観察時とモニタ観察時の輝度比が大きい場合には、本発明がより効果的に作用する。また、白など無彩色に対しては順応ファクタを高くし、有彩色については順応ファクタを低く（より光源色を反映した色再現）するようにした。これにより、輝度の高さによりコントラストが強く見える人間の視覚特性を反映した色変換が再現可能となる。【００６６】
尚、第１の実施の形態同様、ステップＳ４２では、画像観察視環境に於ける照明白色点として、CIE 1931 XYZ表色系の値(XW',YW',ZW')を算出している。しかし式（18）を省略して、錐体応答の空間に於ける値(LW',MW',SW')を、画像観察視環境に於ける照明白色点として取得してもよい。この場合、ステップＳ４５の式（24）も省略することができ、ステップＳ４２で取得した(LW',MW',SW')をそのまま式（28）、（30）に使用することができる。従って処理を軽くすることが出来る。【００６７】
（第２の実施の形態の変形例）
尚、ステッップＳ４７からＳ５３で使用したマトリックスをまとめた以下の２つのマトリックス
【数３７】
2008050761.tif 000039
【数３８】
2008050761.tif 000040
を予めステップＳ４５で算出しておいてもよい。この場合、ステップＳ４７、Ｓ４８、Ｓ５２、Ｓ５３の処理を省略でき、ステップＳ４６の後はステップＳ４９に進む。ステップＳ４９で肯定された場合は、ステップＳ５０に進み、Mchrの代わりにMchr'を使用して、RGBからR'G'B'に変換し、ステップＳ５４に進む。ステップＳ４９で否定された場合は、ステップＳ５１に進み、Machrの代わりにMachr'を使用して、RGBからR'G'B'に変換し、ステップＳ５４に進む。このようにすれば、ステップＳ４７、Ｓ４８、Ｓ５２、Ｓ５３の処理を省略できるので、処理を軽くできる。【００６８】
−第３の実施の形態−
第３の実施の形態では、撮影シーンの照明条件と、画像観察時の照明条件が予め決まっている場合に、決まったマトリックスを予め求めデジタルカメラに組み込んでおく例について説明する。このようにすれば、カメラ内では単一のマトリックスで処理すればよいので、カメラ内での処理を軽くできる。【００６９】
図６は、第３の実施の形態において、色変換マトリックスを求める処理のフローチャートを示す図である。色変換マトリックスを求める処理は、第２の実施の形態の図４と同様な構成で、パーソナルコンピュータ１が実行する。【００７０】
まず、ステップＳ４１からステップＳ４５までは第２の実施の形態と同様であるので説明を省略する。【００７１】
ステップＳ７１からステップＳ７７までは照明の色（白色）と撮影シーンに頻出する色を合わせてN色の予め定めた色パッチについて処理を行う。N色のパッチの選択は、生物顕微鏡用のカメラの様に、頻出する色に特徴が有る場合には、実被写体から頻出する色（有彩色）を抽出する。この場合、無彩色のパッチとして照明の色を取得する。もしくは市販のカラーチャートを使用して、無彩色、有彩色をそれぞれ取得してもよい。【００７２】
まず、ステップＳ７１では、顕微鏡にて上記N色のi番目の色パッチの撮影画像を読み込み、色パッチのカメラRGBを取得する。ステップＳ７２からステップＳ７７の処理は、第２の実施の形態のステップＳ４７からステップＳ５２の処理において、画素をパッチに置き換えた内容と同様であるので、説明を省略する。【００７３】
ステップＳ７１からステップＳ７７までの処理を全パッチに対して繰り返す。ステップＳ７８にて全パッチの処理が終了したか否かの判断が肯定されたら、ステップＳ７９に進む。【００７４】
ステップＳ７９では、ステップＳ７２で算出した(X(i),Y(i),Z(i))(i=0,1,...,N-1)を用いて
【数３９】
2008050761.tif 000041
にて得られる(Xcalc'(i),Ycalc'(i),Zcalc'(i))(i=0,1,...,N-1)と、ステップＳ７７で算出した(X'(i),Y'(i),Z'(i))(i=0,1,...,N-1)の色差平均が最小になるよう、３×３のマトリックスMXYZ→X'Y'Z'を最適化する。すなわち、マトリックスMXYZ→X'Y'Z'のパラメータを逐次変更しながら、色差平均が最小になるようなマトリックスMXYZ→X'Y'Z'を求める。【００７５】
次に、ステップＳ８０に進み、次式よりRGB→R'G'B'の変換マトリックスMRGB→R'G'B'を算出する。
【数４０】
2008050761.tif 000042
【００７６】
尚、MXYZ→sRGBはXYZから出力画像の色空間RGBに変換するマトリックスである。出力画像の色空間をsRGBとすれば、規格で決まっているマトリックスを使用する。【００７７】
次に、上記で求めた変換マトリックスMRGB→R'G'B'をデジタルカメラに搭載し、デジタルカメラ内部で色変換を行う例について説明する。図７は、本実施の形態のデジタルカメラ２０の構成を示す図である。デジタルカメラ２０は、制御部２１、撮像部２２、メモリ２３などから構成される。デジタルカメラ２０にはモニタ２４が接続され、デジタルカメラ２０から出力される顕微鏡像を観察することができる。【００７８】
制御部２１は、マイクロプロセッサおよび周辺回路などから構成され、メモリ２３内に格納されたプログラムを実行することによりデジタルカメラ２０の各種の制御および後述する画像処理を行う。メモリ２３は、プログラムの格納や画像処理時の一時的なバッファーの役割をなすとともに前述した変換マトリックスMRGB→R'G'B'が格納されている。従って、メモリ２３の一部は不揮発性メモリで構成されている。【００７９】
撮像部２２は、ＣＣＤなどの撮像素子（イメージセンサ）から構成され、撮影レンズなどの撮影光学系（不図示）を介して入射する被写体からの光束を受光して電気信号に変換し、変換された電気信号をデジタル変換して画像データとして出力する。撮像部２２は、ＲＧＢカラーフィルターを有し、ＲＧＢ表色系の画像データを出力する。【００８０】
図８は、デジタルカメラ２０の制御部２１が行う処理のフローチャートを示す図である。【００８１】
まず、ステップＳ１０１にて、撮像部２２が撮影（撮像）した撮影画像を取得する。撮影画像はカメラの分光感度に固有の色空間RGBで記述されている。【００８２】
次にステップＳ１０２に進み、各画素について、図６のステップＳ８０で予め算出して予めデジタルカメラ２０のメモリ２３に格納してあるマトリックスMRGB→R'G'B'を用いてRGB→R'G'B'に変換する。
【数４１】
2008050761.tif 000043
【００８３】
ステップＳ１０３にて、全画素についてステップＳ１０２の色変換が終了したことが肯定されると、ステップＳ１０４に進み、変換済み画像をモニタ２４に出力して終了する。【００８４】
以上の様にして、第２の実施の形態と同様に、輝度の高さによりコントラストが強く見える人間の視覚特性を反映した色変換が再現可能となる。しかも、デジタルカメラ２０内で軽い簡単な処理で実現できる。【００８５】
なお、図８の処理は、図６の処理に引き続きパーソナルコンピュータ１内で処理するようにしてもよい。【００８６】
また、本実地の形態では、色変換方法としてマトリックスを用いた例について説明した。しかしマトリックスの代わりにLUTで色変換を行ってもよい。この場合、図６のステップＳ７９でXYZ→X'Y'Z'の変換LUTを作成し、ステップＳ８０でRGB→XYZの変換と、X'Y'Z'→R'G'B'を使用してRGB→R'G'B'の変換LUTを作成する。図８のステップＳ１０２では、MRGB→R'G'B'の代わりに前記RGB→R'G'B'の変換LUTを使用して色変換を行う。【００８７】
処理は、図６のステップＳ７９の「XYZ→X'Y'Z'の変換マトリックスMXYZ→X'Y'Z'作成」を「XYZ→X'Y'Z'変換LUT作成」に、ステップＳ８０の「変換マトリックスMRGB→R'G'B'作成」 を「RGB→R'G'B'変換LUT作成」に書き換え、図８のステップＳ１０２を「各画素について、MRGB→R'G'B'用いてRGB→R'G'B'」を「各画素について、RGB→R'G'B'変換LUT用いてRGB→R'G'B'」に書き換えればよい。【００８８】
−第４の実施の形態−
第２の実施の形態では、まず有彩色と無彩色に夫々異なる順応ファクタを使用して２つの色変換を決定する。次に、変換対象画像の全画素について有彩色か無彩色かに応じて色変換を行うことにより、輝度が高い場合にコントラストが高く見える視覚特性を簡単に再現する例について説明した。【００８９】
第４の実施の形態も、被写体の輝度が高い場合に精度良く、かつ簡単な色変換を行うことが目的であるが、第４の実施の形態では、変換対象画素の状態に応じて２種類の色変換マトリックスを組み合わせて色変換を行う。図９は、第４の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。図９の処理は、第２の実施の形態の図４と同様な構成で、パーソナルコンピュータ１が実行する。【００９０】
ステップＳ２１にて、撮影画像を取得する。撮影画像はカメラの分光感度に固有の色空間RGBで記述されている。【００９１】
次にステップＳ４１に進む。ステップＳ４１〜Ｓ４５は第２の実施の形態の図５と同様の処理であるので、説明を省略する。【００９２】
次にステップＳ６１に進み、ステップＳ２１で取得した画像の各画素について、ステップＳ６１〜Ｓ６７の処理を行う。ステップＳ６１は図５のステップＳ４７と同様の処理であるので、説明を省略する。【００９３】
次にステップＳ６２に進み、変換対象画素の彩度C*を計算する。彩度は、XYZからCIE LABに変換し、C*=(a*2+b*2)1/2を算出する。【００９４】
ステップＳ６３に進み、次式より変換対象画素に対する変換マトリックスMを算出する。α(C*)を用いて次式のように重みづけがなされる。
【数４２】
2008050761.tif 000044
図１０は、関数α(C*)の例を示すグラフである。α(C*)は、図１０の様に、α(0)=0（無彩色）かつC*が大きいとき（有彩色）でα(C*)=1となる関数である。α(0)からα(C*0)の間は0から１の間で単調増加する。C*0は、予め無彩色の範囲として定めた彩度の範囲である。【００９５】
次にステップＳ６４に進む。ステップＳ６４は図５のステップＳ４８と同様の処理であるので、説明を省略する。次にステップＳ６５に進み、ステップＳ６３で算出したマトリックスMを用いて白色点変換を行う。
【数４３】
2008050761.tif 000045
【００９６】
次にステップＳ６６に進む。ステップＳ６６は図５のステップＳ５２と同様の処理なので説明を省略する。【００９７】
次にステップＳ６７に進み、X'Y'Z'から出力画像の色空間R'G'B'に変換する。MXYZ→sRGBはXYZから出力画像の色空間RGBに変換するマトリックスであり、出力画像の色空間をsRGBとすれば、規格で決まっているマトリックスを使用すればよい。
【数４４】
2008050761.tif 000046
【００９８】
（第４の実施の形態の変形例）
尚、ステッップＳ６４からＳ６７で使用したマトリックスをまとめた以下のマトリックス
【数４５】
2008050761.tif 000047
をステップＳ６３の後に算出し、ステップＳ６４からＳ６７の処理の代わりに
【数４６】
2008050761.tif 000048
の１つのステップとして、処理を軽くしてもよい。【００９９】
ステップＳ６８では、全画素処理が終了していなければ、上記ステップＳ６１に戻し、ステップＳ６１〜ステップＳ６７の処理を繰り返す。全画素終了したら、ステップＳ３０に進み画像を出力して終了する。【０１００】
本実施の形態の様に、白など無彩色に対しては順応ファクタを高くし、有彩色については順応ファクタを低く（より光源色を反映した色再現）することで、輝度の高さによりコントラストが強く見える人間の視覚特性を、マトリックス演算の様な簡単な色変換で再現可能である。また第２の実施の形態と違って２つのマトリックスで色変換せず、変換対象画素の彩度に応じて２種類のマトリックスを組み合わせているので、より正確な色再現が可能である。【０１０１】
上記第４の実施の形態においても、撮影シーンの照明条件と画像観察時の照明条件が予め決まっている場合に、第３の実施の形態と同様にして、決まったマトリックスを予め求めデジタルカメラに組み込んでおくことが可能である。この場合、図６の処理において、図６のステップＳ７２からステップＳ７７の処理を、図９のステップＳ６１からステップＳ６６と同様な処理に置き換え、ステップＳ６１の処理にて、「各画素について」を「各パッチについて」に読み替えて実行すればよい。このようにすれば、第３の実施の形態と同様に、カメラ内では単一のマトリックスで処理すればよいので、カメラ内での処理を軽くできる。【０１０２】
−第５の実施の形態−
第２の実施の形態では、まず有彩色と無彩色に夫々異なる順応ファクタを使用して２つの色変換を決定した。次に変換対象画像の全画素について有彩色が無彩色かに応じて色変換を行うことにより、輝度が高い場合にコントラストが高く見える視覚特性を簡単に再現する例について説明した。また第４の実施の形態では、変換対象画素の状態に応じて２種類の色変換マトリックスを組み合わせて色変換を行う例について説明した。【０１０３】
第５の実施の形態も、第２、第３、第４の実施の形態同様、輝度による見えへの影響を考慮しつつ、かつ簡単な色変換を行うことが目的であるが、第５の実施の形態では、変換対象画素の状態に応じて最適な順応ファクタを算出し、算出した順応ファクタに応じて変更した色変換マトリックスを用いて色変換を行う。【０１０４】
図１１は、第５の実施の形態におけるパーソナルコンピュータ１が実行する画像処理プログラムのフローチャートを示す図である。図１１の処理は、第２の実施の形態の図４と同様な構成で、パーソナルコンピュータ１が実行する。【０１０５】
ステップＳ２１、Ｓ４１〜Ｓ４４は、第２の実施の形態の図５と同様の処理であるので、説明を省略する。【０１０６】
ステップＳ９１に進み、変換対象画素の彩度C*の関数として、撮影シーンの照明条件に於ける順応ファクタD(C*)、画像観察時の視環境に於ける順応ファクタD'(C*)を作成する。図１２は関数D(C*)の一例を示すグラフであり、図１３は関数D'(C*)の一例を示すグラフである。D(C*)、D'(C*)は図１２、図１３の様に、D(0)=Dachr、D'(0)=Dachr'、（無彩色）、かつC*が大きいとき（有彩色）でD(C*)=Dchr、D'(C*)=Dchr'となる関数であり、例えば以下の様な関数を用いる。
【数４７】
2008050761.tif 000049
【数４８】
2008050761.tif 000050
C*0は、予め無彩色の範囲として定めた彩度の範囲である。【０１０７】
ステップＳ６１〜Ｓ６２は、第４の実施の形態の図９と同様の処理であるので、説明を省略する。【０１０８】
次にステップＳ９２に進み、ステップＳ６２で算出した変換対象画素の彩度C*とステップＳ９１で算出した順応ファクタD(C*)、D'(C*)を用いて、変換対象画素に於ける撮影シーンの照明条件に対する順応ファクタD、画像観察時の視環境に対する順応ファクタD'を算出する。【０１０９】
ステップＳ９３では、変換対象画素の彩度に応じた色変換マトリックスMを算出する。
【数４９】
2008050761.tif 000051
但し、
【数５０】
2008050761.tif 000052
【数５１】
2008050761.tif 000053
ただし、
【数５２】
2008050761.tif 000054
【数５３】
2008050761.tif 000055
【０１１０】
ステップＳ６４〜Ｓ６８、及びステップＳ３０は、第４の実施の形態の図９と同様の処理であるので、説明を省略する。【０１１１】
尚、D(C*)、D'(C*)は上述の式(47)、(48)や、図１２、図１３に示したものに限らない。【０１１２】
また、上記第５の実施の形態においても、撮影シーンの照明条件と画像観察時の照明条件が予め決まっている場合に、第３の実施の形態と同様にして、決まったマトリックスを予め求めデジタルカメラに組み込んでおくことが可能である。この場合、図６の処理において、図６のステップＳ４５の処理を、図９のステップＳ９１と同様な処理に置き換え、図６のステップＳ７２からＳ７７の処理を、図９のステップＳ６１からＳ６６と同様な処理に置き換え、ステップＳ６１の処理にて、「各画素について」を「各パッチについて」に読み替えて実行すればよい。このようにすれば、第３の実施の形態と同様に、カメラ内では単一のマトリックスで処理すればよいので、カメラ内での処理を軽くできる。【０１１３】
尚、第２から第５の実施の形態では、XYZからLMSへの変換マトリックスとして、MCAT02を使用したが、色順応変換に適した色空間に変換するマトリックスであれば、他の色変換マトリックスを使用してもよい。【０１１４】
尚、第２から第５の実施の形態では、無彩色の順応ファクタは常に1.0としたが、顕微鏡観察時とモニタ観察時の輝度比が大きい程、コントラストの違いが顕著になる。そこで、無彩色の順応ファクタDachrを常に1.0とするのではなく、輝度比が特に大きい場合（例えば一桁以上異なる場合、即ち輝度10倍の場合）に1.0とし、輝度比がその間の場合（1倍〜10倍の場合）には、Dachrを有彩色の順応ファクタDchrと1.0の間の値で調整してもよい。調整は、簡単には線形補間でもよいし、図１９の様に別の形の関数で決定してもよいし、その他の関数を使用してもよい。図１９は、輝度比対順応ファクタDchrの非線形関数の例を示す図である。予め輝度比を変えた実験から決定したDachrの最適値をテーブルとして記憶しておき、実際の輝度比に応じてテーブルから近い条件のDachrを読み出し、補間して求めてもよい。【０１１５】
尚、第２の実施の形態同様、第３から第５の実施の形態に於いても、ステップＳ４２では、画像観察視環境に於ける照明白色点として、CIE 1931 XYZ表色系の値(XW',YW',ZW')を算出している。しかし式（18）を省略して、錐体応答の空間に於ける値(LW',MW',SW')を、画像観察視環境に於ける照明白色点として取得してもよい。第５の実施の形態の場合、図１１のステップＳ９３の式（53）も省略することができ、ステップＳ４２で取得した(LW',MW',SW')をそのまま式（51）に使用することができる。従って処理を軽くすることが出来る。【０１１６】
尚、第２〜第５の実施の形態では、彩度C*に依存して順応ファクタやマトリックスを変化させる例で説明したが、彩度C*の代わりにHSV表色系の次式のSなど、彩度を表すパラメータであれば他のパラメータを用いてもよい。
【数５４】
2008050761.tif 000056
【０１１７】
C*を用いる場合XYZ→CIE LABの変換が必要になるが、Sを用いればRGBから算出でき計算が簡略化できる。またCIE LABのC*は輝度により数値範囲が変化するが、Sの場合0.0≦S≦1.0となる為、第４の実施の形態のα(C*)や第５の実施の形態のD(C*)、D'(C*)で使用する無彩色の範囲（上記実施の形態に於けるC*0）の扱いが容易になる。【０１１８】
−第６の実施の形態−
第１の実施の形態では、撮影した画像をモニタやプリントで観察する際、撮影シーンの照明条件と画像観察時の視環境の違いを考慮して色変換する場合を例に説明した。しかし、上記実施の形態で説明した内容は、視環境の異なる条件下の色の見えを色の見えモデル（色順応式も含む）を用いて予測する場合であれば適用可能である。【０１１９】
例えば上記実施の形態の他に、第１の視環境（例えばモニタ表示）にて観察している画像を第２の視環境（印刷したハードコピーの観察）にて観察する場合に、両者の見えが近くなる様に色の見えを考慮して変換する場合や、撮影時に予め標準視環境で観察されることを前提に色順応変換されていた画像を実際の視環境に合わせて更に色順応変換する場合などにも適用できる。【０１２０】
第６の実施の形態では、第１の視環境（例えばモニタ表示）にて観察する為に作成された画像を第２の視環境（印刷したハードコピーの観察）にて観察する場合について説明する。第１の実施の形態の図１と同様にパーソナルコンピュータ１が処理を実行する。【０１２１】
本実施の形態のパーソナルコンピュータ１は、すでに第１の視環境に関する情報に基づいて作成された画像データを、人間の色順応特性を考慮して、第１の視環境とは異なる第２の視環境に関する情報に応じた色の見えを再現する画像データに変換する。【０１２２】
前述した第１の視環境に関する情報とは、入力画像データ作成時に使用された、入力画像データ観察に使用する機器や周囲観察環境の照明に関する条件である。第２の視環境に関する情報とは、実際に出力画像データ観察に使用する機器や周囲観察環境の照明に関する条件である。これら、第１の視環境に関する情報や第２の視環境に関する情報は、照明光源の輝度や白色点、周囲の明るさなど、見えに影響を与える視野周辺の環境に関する情報である。【０１２３】
図１４は、第６の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。【０１２４】
ステップＳ８１では、第１の視環境で観察されることを前提に色変換されていた第１の画像データを取得する。ステップＳ８２では、取得した画像データから第１の画像データの作成に用いられた第１の視環境に関する条件を読み込み、第１の視環境として設定する。画像データに第１の視環境に関するデータが無い場合は、sRGB標準視環境（白色点=D65、輝度Ydevice0=80Cd/m2、周囲照明白色点=D50、周囲輝度YSW0=4.1Cd/m2）を読み込み、デバイス白色点と周囲照明白色点の混合として第１の視環境に於ける照明白色点(XW,YW,ZW)、輝度YW0を算出する。【０１２５】
sRGB標準視環境と異なる視環境で観察するために作成された画像の場合には、ユーザ（観察者）がキーボードで入力した画像作成に使用したデバイスの白色点(Xdevice,Ydevice,Zdevice)、輝度Ydevice0、および周囲照明白色点(XSW,YSW,ZSW)、周囲照明輝度YSW0を読み込み、デバイス白色点と周囲照明白色点の混合として第１の画像観察視環境に於ける照明白色点(XW,YW,ZW)、輝度YW0を算出する。算出方法は、第１の実施の形態の図３のステップＳ２３と同様であるので、説明を省略する。【０１２６】
ステップＳ２３では、出力画像観察時の照明に関する条件を読み込み、第２の視環境として設定する。ステップＳ２３の処理は、第１の実施の形態の図３のステップＳ２３と同様であるので説明を省略する。ステップＳ２４〜Ｓ３０の処理は、第１の実施の形態の図３のステップＳ２４〜Ｓ３０と同様であるので説明を省略する。【０１２７】
尚、ステップＳ２６で算出するDh*(h*) 、Dh*'(h*)は、それぞれ第１の視環境、ならびに第２の視環境と変換対照画素の色相h*に応じて算出する。【０１２８】
以上のようにして、ある視環境で観察されることを想定して生成された画像データを、その視環境とは異なる視環境においても適切に撮影シーンの色の見えを再現し、観察することができる。この場合、第１の実施の形態と同様に、変換対象画素の色すなわち色相により最適な順応ファクタを使用して変換している為、正確な色の見えを再現することができる。【０１２９】
−第７の実施の形態−
第１の実施の形態では、撮影画像を入力として、撮影シーンを見ているときの見えを出力画像観察時に再現する例で説明した。また第６の実施の形態では、第１の視環境で観察することを前提に作成された画像を、第１の視環境とは異なる第２の視環境にて観察する場合を例に説明した。【０１３０】
第７の実施の形態では、第１の視環境で観察している画像を別の第２の視環境で観察する場合について説明する。例えば、モニタ観察画像をプリント出力する場合やプロジェクタで投影して観察する場合に、モニタ観察画像と見えを合わせる場合である。【０１３１】
図１５は、第７の実施の形態におけるパーソナルコンピュータが実行する処理のフローチャートを示す図である。図１５の処理は、第１の実施の形態の図１と同様にパーソナルコンピュータ１が処理を実行する。【０１３２】
ステップＳ２１、ステップＳ２４〜Ｓ３０は第１の実施の形態（図３）と同様の処理なので説明を省略する。【０１３３】
ステップＳ１３２では、ユーザがキーボード等で入力した第１の画像観察時の視環境に関する情報を取得する。「CRT観察（sRGB標準視環境）」が選択されると、デバイスの白色点(Xdevice,Ydevice,Zdevice)としてD65のXYZ値、輝度Ydevice0=80Cd/m2、周囲照明白色点(XSW,YSW,ZSW)としてD50のXYZ値、周囲照明輝度YSW0=4.1Cd/m2を設定する。もしくは、ユーザが入力した(Xdevice,Ydevice,Zdevice) 、Ydevice0 、(XSW,YSW,ZSW)、YSW0の各値を取得する。【０１３４】
取得した前記４つの値から、入力画像観察時の視環境に於いて順応している照明白色点の錐体応答(LW,MW,SW)と輝度YW0を決定する。具体的には、画像観察に使用するデバイス光源と周囲照明光源の影響の比率Rmixを用いて、次のように算出する。
【数５５】
2008050761.tif 000057
【数５６】
2008050761.tif 000058
【数５７】
2008050761.tif 000059
【数５８】
2008050761.tif 000060
【０１３５】
ステップＳ１３３では、ユーザがキーボード等で入力した出力画像観察時の視環境に関する情報を取得する。即ちユーザ（観察者）がキーボード等で入力した出力画像観察に使用するデバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'、および周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'を読み込み、デバイス白色点と周囲照明白色点の混合として出力画像観察時の視環境に於いて順応している照明白色点の錐体応答(LW',MW',SW')、輝度YW0'を算出する。【０１３６】
例えば、出力画像をプリントした印刷物で観察する場合には、ユーザ（観察者）がキーボード等で入力した出力画像観察に使用する照明の白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'、および周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'を使用する。出力画像をプロジェクタに表示する場合は、ユーザがキーボード等で入力した以下の値を使用する。【０１３７】
デバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'は、スペックに書かれているデータを用いるか、ユーザ実測値を入力する。周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'は、暗室で観察する場合はYSW0'=0、(XSW',YSW',ZSW')= (Xdevice',Ydevice',Zdevice')を入力する。暗室でない場合は、ユーザが周囲視環境を実測し、実測値を入力する。【０１３８】
(LW',MW',SW')、YW0'の算出は、第１の実施の形態同様、出力画像観察に使用するデバイス光源と周囲照明光源の影響の比率Rmix'を用いて、次のように算出する。
【数５９】
2008050761.tif 000061
【数６０】
2008050761.tif 000062
【数６１】
2008050761.tif 000063
【数６２】
2008050761.tif 000064
【０１３９】
ステップＳ２６で算出するDh*(h*) 、Dh*'(h*)は、それぞれ第１の視環境、ならびに第２の視環境と変換対照画素の色相h*に応じて算出する。【０１４０】
以上のようにして、異なる視環境間で観察する場合に、色の再現を一致させることができる。この場合、第１の実施の形態と同様に、変換対象画素の色すなわち色相により最適な順応ファクタを使用して変換している為、正確な色の見えを再現することができる。【０１４１】
−第８の実施の形態−
第１から第７の実施の形態は、各画素の色変換を行う際、各画素の色（色相、彩度等）に応じて順応ファクタDを変化させる例で説明した。しかし、第８の実施の形態では、画素毎に変化させるのではなく、画像データにおいて複数の近傍画素を一つの領域として捉え、その領域毎に代表色を決定し、代表色の色（色相、彩度等）によってその領域内の画素の変換に適応する順応ファクタDを決定する。すなわち、画像データの色（色相、彩度等）の分布を解析し、解析した色の分布に基づいて順応ファクタを決定する。【０１４２】
領域分割の方法としては、エッジ抽出結果により被写体毎に領域分割する場合や、画像をブロックに分割する方法などが考えられる。エッジ抽出結果により被写体毎に領域分割している場合、Dの変化と被写体の輪郭は一致する為、Dの変化による境界が目立つなど不自然さが目立たない。【０１４３】
ブロック分割してブロック領域毎にDを変える場合、領域の境界では使用するDに不連続が生じる。この様に分割領域境界近傍においてDの変化による不自然さが生じる場合、境界領域の複数画素についてはDが連続につながるよう隣接する領域の代表順応ファクタDから補間して求めてもよい。【０１４４】
図１６は、第８の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。図１６の処理は、第１の実施の形態の図１と同様にパーソナルコンピュータ１が処理を実行する。【０１４５】
ステップＳ２１〜Ｓ２３は第１の実施の形態（図３）と同様の処理なので説明を省略する。【０１４６】
ステップＳ１４１では、撮影画像を複数領域に分割する。領域の分割は、例えばブロック分割や、エッジ抽出により抽出された複数の閉空間とその外側という分割方法がある。ステップＳ１４２では、各領域毎に、代表色を決定する。代表色は、領域の中心近傍画素から抽出したり、領域内の画素値を平均した色を使用したり、領域内に頻出する色を求めたりして、代表色を決定する。【０１４７】
ステップＳ１４３では代表色の色相角h*を算出する。ステップＳ１４４では、代表色の色相角h*からその領域の画素に適応する順応ファクタDh*(h*) 、Dh*'(h*)を算出する。すなわち、順応ファクタDh*(h*) 、Dh*'(h*)の２次元分布を算出する。【０１４８】
ステップＳ２４〜Ｓ２８は第１の実施の形態（図３）と同様の処理なので説明を省略する。各画素の色変換で用いる順応ファクタは、ステップＳ１４４で算出した変換対象領域毎に決定したDh*(h*) 、Dh*'(h*)や、近傍の代表色との間で補間して求めたDh*(h*) 、Dh*'(h*)を使用する。【０１４９】
ステップＳ１４５では、領域内の全画素の変換が終了したか否か判断し、否であればステップＳ２４に戻って処理を続ける。肯定されたらステップＳ１４６に進む。ステップＳ１４６では、分割した全領域についてステップＳ１４３〜ステップＳ１４５の処理が終了したか否かの判断を行う。否であればステップＳ１４３に戻って処理を続ける。肯定されたらステップＳ３０に進み、画像を出力して終了する。【０１５０】
以上のようにして、領域毎に領域の代表色に応じた最適順応ファクタを決めるようにした。これにより、画素毎に色相判断をする必要がなくなり、処理を簡略化することができる。しかも、変換対象画素の大体の色に応じた最適な順応ファクタを使用して色変換をするので、正確な色の見えも再現することができる。【０１５１】
（第８の実施の形態の変形例）
ステップＳ１４１の画像分割を行わず、画像全体で代表色を１つ決定してもよい。図１７は、この変形例のパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。図１７の処理は、第１の実施の形態の図１と同様にパーソナルコンピュータ１が処理を実行する。【０１５２】
ステップＳ２１〜Ｓ２３は第１の実施の形態（図３）と同様なので説明を省略する。ステップＳ１５２では、画像全体の代表色を決定する。代表色の決定方法は、第８の実施の形態（図１６）のステップＳ１４２と同様なので、説明を省略する。ステップＳ１４３〜Ｓ１４４、Ｓ２４、Ｓ２７、Ｓ２８は、第８の実施の形態（図１６）と同様なので、説明を省略する。【０１５３】
ステップＳ１４５では、全画素の変換が終了したか否か判断し、否であればステップＳ２４に戻って処理を続ける。肯定されたらステップＳ３０に進み、画像を出力して終了する。この様にすれば、領域毎にも色相判断をする必要がなくなり、より一層処理を簡略化することができる。【０１５４】
尚、第８の実施の形態では、第１の実施の形態の様に、撮影画像を画像観察時の視環境に応じて、撮影時のシーンと出力画像観察時の見えが近くなるように色変換する場合に、領域毎に代表色を決定し代表色に基づいてDを変化させる例を説明した。しかし第６の実施の形態の様に、第１の視環境で観察する前提で作成された画像を、第２の視環境で観察した場合に、撮影時のシーンと見えが近くなる様に色変換する場合にも、第８の実施の形態で説明した様に、領域毎に代表色を決定し代表色に基づいてDを変化させてもよい。また、第８の実施の形態の変形例のように画像全体で代表色を１つ決定してもよい。【０１５５】
また第７の実施の形態の様に、第１の視環境で入力画像を観察した見えと、第２の視環境で出力画像を観察した見えが近くなる様に出力画像を色変換する場合にも、第８の実施の形態で説明した様に、領域毎に代表色を決定し代表色に基づいてDを変化させてもよい。また、第８の実施の形態の変形例のように画像全体で代表色を１つ決定してもよい。【０１５６】
−第９の実施の形態−
第９の実施の形態は、撮影画像をデジタルカメラの背面液晶で表示する場合を例に説明する。図１８は、本実施の形態のデジタルカメラ３０の構成を示す図である。デジタルカメラ３０は、制御部３１、撮像部３２、メモリ３３、モニタ３４、測光部３５などから構成される。モニタ３４は、デジタルカメラ３０の背面に設けられた液晶モニタである。【０１５７】
制御部３１は、マイクロプロセッサおよび周辺回路などから構成され、メモリ３３内に格納されたプログラムを実行することによりデジタルカメラ３０の各種の制御および以下で説明する画像処理を行う。撮像部３２は、ＣＣＤなどの撮像素子（イメージセンサ）から構成され、撮影レンズなどの撮影光学系（不図示）を介して入射する被写体からの光束を受光して電気信号に変換し、変換された電気信号をデジタル変換して画像データとして出力する。撮像部３２は、ＲＧＢカラーフィルターを有し、ＲＧＢ表色系の画像データを出力する。【０１５８】
デジタルカメラ３０の制御部３１が実行する画像処理は図１６と同様なので、図１６を用いて説明する。【０１５９】
ステップＳ２１では、撮像部３２により撮影された撮影画像を取得する。ステップＳ２２では撮影シーンに関する情報を取得する。詳細は第８の実施の形態と同様なので説明を省略する。ステップＳ２３では画像観察時の視環境に関する情報を取得する。【０１６０】
ここで画像観察デバイスはデジタルカメラ３０のモニタ３４である。従って、予めメモリ３３に記憶されたモニタ３４の白色点、並びに輝度を、デバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'として夫々読み込む。また、測光部３５による測光機能を利用して周囲輝度を算出し、周囲照明輝度YSW0'として設定する。周囲の白色点(XSW',YSW',ZSW')は、直前に撮影した画像のWBデータから光源を決定し、光源の白色点を読み込む。【０１６１】
ユーザが「太陽光、曇天、電球、蛍光灯」などのメニューから観察時の光源を選択した場合は、選択した光源に対応する白色点を読み込む。デバイス白色点と周囲照明白色点の混合として画像観察視環境に於ける照明白色点(XW',YW',ZW')、輝度YW0'を算出する方法は、第８の実施の形態と同様なので説明を省略する。【０１６２】
ステップＳ１４１〜ステップＳ３０は第８の実施の形態と同様である。ステップＳ３０で出力した画像をモニタ３４に表示すれば、視環境に関する情報に関してユーザが入力を省略しても、簡単にデジタルカメラ３０の背面に設けられたモニタ３４で実シーンと見えの近い画像が表示される。【０１６３】
第９の実施の形態では、撮影画像をデジタルカメラの背面液晶で表示する場合に、第８の実施の形態で説明した画像処理方法を用いて、撮影シーンと見えが近い画像をカメラのモニタに表示できることを示した。同様に、第１の実施の形態と同様の処理をカメラ内で行ってもよい。第９の実施の形態との違いは、図１６のステップＳ１４１からＳ１４６の代わりに、第１の実施の形態（図３）のステップＳ２４からＳ２９を行えばよい。又、第８の実施の形態の変形例（図１７）と同様の処理をカメラ内で行う場合は、第９の実施の形態の図１６のステップＳ１４１からＳ１４６の代わりに図１７のステップＳ１５２からＳ１４５を行えばよい。【０１６４】
なお、順応ファクタを用いた色順応変換の式は、第１の実施の形態〜第８の実施の形態で説明した方法（例えば第１の実施の形態の式(8)〜(13)）に限らない。即ち、入力画像撮影時、もしくは、入力画像観察時および出力画像観察時に視覚が順応している白の錐体応答（例えば、第１の実施の形態の式(10)、(11)）が下記式で表されればよい。
LW,D=f(LW,D(h*,C*,L*))
MW,D=f(MW,D(h*,C*,L*))
SW,D=f(SW,D(h*,C*,L*))
LW,D'=f(LW',D'(h*,C*,L*))
MW,D'=f(MW',D'(h*,C*,L*))
SW,D'=f(SW',D'(h*,C*,L*))
尚、(LW,MW,SW)、(LW',MW',SW')はそれぞれ画像撮影時もしくは入力画像観察時、及び出力画像観察時に視覚が順応している白に対する錐体応答であり、第１の実施の形態で言えば式（１２）、（１３）に相当する。C*は変換対象画素または代表色の彩度、h*は変換対象画素または代表色の色相、L*は変換対象画素または代表色の明度を示す。【０１６５】
−第１０の実施の形態−
第１、第６、第７、第８、第９の実施の形態は、それぞれ変換対象画素の色相もしくは代表色の色相を基に、図２で示されるデータに基づき、順応ファクタを決定していた。一方、本実施の形態では、異なる照明色温度(Tc)毎に、図２に示すような最適な順応ファクタDh*(h*)と色相角h*の関係を示すデータを予め取得して、撮影画像取得時などの入力画像取得時もしくは入力画像観察時の視環境における照明色温度(Tc)と画像観察時の視環境における表示デバイスの照明の色温度(Tc')と色相h*に基づいて順応ファクタDh*(h*,Tc)及びDh*'(h*,Tc')を変化させた。すなわち、画像データの色相と入力画像取得時における照明の色温度と画像観察時の視環境における照明の色温度に基づいて、順応ファクタを決定する。【０１６６】
図２０は、色温度4000Kと10000Kにおける順応ファクタDh*(h*)と色相角h*の最適な関係を示すグラフDh*(h*,Tc)である。これは先に述べた（色の見え実験結果）と同様な手法により、更に色温度を変えることで、求められた知見である。このように、色温度が変われば、同じ色相であっても順応ファクタDh*(h*)が変わる場合がある。そこで、最適な順応ファクタを算出する際、本実施の形態では第1の実施の形態のステップＳ２６の代わりに、順応ファクタDh*(h*,Tc)、Dh*'(h*,Tc')により最適な順応ファクタを算出することにした。【０１６７】
第１の実施の形態との変更点は、１．ステップＳ２２で撮影シーンの照明に関する情報取得で取得された照明の白色点(Xw,Yw,Zw)から色温度(Tc)を算出する工程を追加することと、２．ステップＳ２３で画像観察時の視環境に関する情報取得で取得された照明に関する白色点(Xw',Yw',Zw')から色温度(Tc')を算出する工程が追加されることと、３．ステップＳ２６で算出される最適な順応ファクタDh*(h*)の代わりに、Dh*(h*,Tc)を用いる点、４．ステップＳ２６で算出される最適な順応ファクタDh*'(h*)の代わりに、Dh*'(h*,Tc')を用いる点である。【０１６８】
なお、このような変形例は、第６、第７、第８の実施の形態にも適用可能である。本実施の形態を第６の実施の形態への適用に際し、今までの説明と異なる点は、取得した画像データから第１の画像データの作成に用いられた第１の視環境に関する情報を取得した後、取得した第一の視環境に関する情報の１つである照明の白色点(Xw,Yw,Zw)から色温度(Tc)を算出する点である。【０１６９】
前述の第1の実施の形態に適用した場合の変更点２．〜４．（ステップＳ２３、Ｓ２６）は、同様に留意して処理する。【０１７０】
第７の実施の形態（図１５）への適用に際し留意する点は、１．ステップＳ１３２で入力画像観察時の視環境に関する情報を取得後、取得した照明白色点から色温度(Tc)を算出する点と、２．ステップＳ１３３で出力画像観察時の視環境に関する情報取得後、取得した照明白色点から色温度(Tc')を算出する点である。前述の第１の実施の形態に適用する際の変更点３、４（ステップＳ２６における変更点）は、同様に留意して処理する。【０１７１】
第８の実施の形態への適用に際し、留意する点は、１．ステップＳ１４４で算出される最適な順応ファクタDh*(h*)の代わりに、Dh*(h*,Tc)を用いる点と、２．ステップＳ１４４で算出される最適な順応ファクタDh*'(h*)の代わりに、Dh*'(h*,Tc')を用いる点が相違する点である。【０１７２】
前述の第１の実施の形態に適用する際の変更点１、２（ステップＳ２２、Ｓ２３における変更点）は、同様に留意して処理する。【０１７３】
なお、第９の実施の形態に適用するに際し留意する点は、第８の実施の形態に適用するに際し留意する点と同じなので、説明を省略する。【０１７４】
ところで、白色点から色温度の求め方は、例えば、黒体輻射されるスペクトルに該当する座標をXYZ色座標に表記し、それぞれの白色点と最も近い位置に座標位置がある当該スペクトルに対応する色温度を選び出す。これを色温度とする。他にも、第９の実施の形態ではカメラで補正されるホワイトバランス量から色温度を算出しても良い。【０１７５】
−第１１の実施の形態−
第１、第６、第７、第８、第９の実施の形態は、それぞれ変換対象画素の色（色相等）もしくは、代表色の色（色相等）を基に、図２で示されるデータに基づき、順応ファクタを決定していた。一方、本実施の形態では、異なる彩度(C*)毎に、図２に示すような最適な順応ファクタDh*(h*)と色相角h*の関係を示すデータを予め取得して、変換対象画素もしくは代表色の色相角h*と彩度C*に基づいて順応ファクタDh*(h*,C*)を変化させた。すなわち、画像データの色相と彩度に基づいて、順応ファクタを決定する。【０１７６】
図２１は、各彩度(C*=25,C*=50)における順応ファクタDh*(h*)と色相角h*の最適な関係を示すグラフである。これは先に述べた（色の見え実験結果）と同様な手法により、更に彩度を変えることで、求められた知見である。このように、彩度が変われば、同じ色相であっても順応ファクタDh*(h*)が変わる場合がある。そこで、最適な順応ファクタを算出する際、本実施の形態では第1の実施の形態のステップＳ２６の代わりに、順応ファクタDh*(h*,C*)により最適な順応ファクタを算出することにした。【０１７７】
第1の実施の形態との変更点は、１．ステップＳ２５で変換対象画素の色相角(h*)のほかに彩度(C*)を算出する工程を追加することと、２．ステップＳ２６で算出される最適な順応ファクタDh*(h*)の代わりに、Dh*(h*,C*)を用いる点、３．ステップＳ２６で算出される最適な順応ファクタDh*'(h*)の代わりに、Dh*'(h*,C*')を用いる点である。【０１７８】
なお、このような変形例は、第６、第７、第８、第９の実施の形態にも適用可能である。【０１７９】
ところで、第８及び第９の実施の形態、及び第８の実施の形態の変形例への適用に際し、留意する点は、１．ステップＳ１４３で代表色における色相角h*を算出するのと合わせて、彩度C*も算出することと、２．ステップＳ１４３で求められた代表色の色相角h*と彩度C*に基づいて、ステップＳ１４４では最適順応ファクタをDh*(h*)、Dh*'(h*)の代わりにDh*(h*,C*)とDh*'(h*,C*')を用いて算出する点である。【０１８０】
第６と第７の実施の形態に適用する際留意する点は、第１の実施の形態に適用するに際し留意する点と同じなので、説明を省略する。【０１８１】
−第１２の実施の形態−
第１、第６、第７、第８、第９の実施の形態は、それぞれ変換対象画素の色（色相等）もしくは、代表色の色（色相等）を基に、図２で示されるデータに基づき、順応ファクタを決定していた。第１０の実施の形態では、照明の色温度毎に、図２に示すような最適な順応ファクタDh*(h*)と色相角h*の関係を示すデータを予め取得して、変換対象画素もしくは代表色の色相角h*と入力画像取得時、もしくは入出力画像観察時の照明の色温度Tc,Tc'に基づいて順応ファクタを決定した。第１１の実施の形態では、異なる彩度(C*)毎に、図２に示すような最適な順応ファクタDh*(h*)と色相角h*の関係を示すデータを予め取得して、変換対象画素もしくは代表色の色相h*と彩度C*に基づいて順応ファクタを決定した。【０１８２】
一方、最適順応ファクタDh*(h*,C*)を色温度ごとに色相角(h*)と彩度(C*)に対して取得しておき、撮影画像取得時もしくは入力画像観察時の視環境における照明色温度(Tc)と出力画像観察時の視環境における表示デバイスの照明の色温度(Tc')と変換対象画素もしくは代表色の色相角(h*)と彩度(C*)に基づいて、最適順応ファクタDh*(h*,C*,Tc)、Dh*'(h*,C*,Tc')を決定しても良い。【０１８３】
第１の実施の形態との変更点は、１．ステップＳ２２で撮影シーンの照明に関する情報取得で取得された照明の白色点(Xw,Yw,Zw)から色温度(Tc)を算出する工程を追加することと、２．ステップＳ２３で画像観察時の視環境に関する情報取得で取得された照明に関する白色点(Xw',Yw',Zw')から色温度(Tc')を算出する工程が追加されることと、３．ステップＳ２５で変換対象画素の色相角(h*)のほかに彩度(C*)を算出する工程を追加することと、４．ステップＳ２６で算出される最適な順応ファクタDh*(h*)の代わりに、Dh*(h*,C*,Tc)を用いることと、５．ステップＳ２６で算出される最適な順応ファクタDh*'(h*)の代わりに、Dh*'(h*,C*,Tc')を用いることである。【０１８４】
なお、色温度Tc,Tc'の求め方など、その他の留意する点は、第１０、第１１の実施の形態と同様である。【０１８５】
なお、このような変形例は、第６、第７、第８、第９の実施の形態にも適用可能である。【０１８６】
本実施の形態を第６の実施の形態への適用に際し留意する点は、１．ステップＳ８２で取得した画像データから第１の画像データの作成に用いられた第１の視環境に関する情報を取得した後、取得した第一の視環境に関する情報の１つである照明の白色点(Xw,Yw,Zw)から色温度(Tc)を算出すること、２．ステップＳ２３で画像観察時の視環境に関する情報取得で取得された照明に関する白色点(Xw',Yw',Zw')から色温度(Tc')を算出する工程が追加されることと、３．ステップＳ２５で変換対象画素の色相角(h*)のほかに彩度(C*)を算出する工程を追加することと、４．ステップＳ２６で算出される最適な順応ファクタDh*(h*)の代わりに、Dh*(h*,C*,Tc)を用いることと、５．ステップＳ２６で算出される最適な順応ファクタDh*'(h*)の代わりに、Dh*'(h*,C*,Tc')を用いることである。【０１８７】
第７の実施の形態への適用に際し留意する点は、１．ステップＳ１３２で入力画像観察時の視環境に関する情報を取得後、取得した照明白色点から色温度(Tc)を算出することと、２．ステップＳ１３３で出力画像観察時の視環境に関する情報取得後、取得した照明白色点から色温度(Tc')を算出することと、３．ステップＳ２５で変換対象画素の色相角(h*)のほかに彩度(C*)を算出する工程を追加することと、４．ステップＳ２６で算出される最適な順応ファクタDh*(h*)の代わりに、Dh*(h*,C*,Tc)を用いることと、５．ステップＳ２６で算出される最適な順応ファクタDh*'(h*)の代わりに、Dh*'(h*,C*,Tc')を用いることである。【０１８８】
第８及び第９の実施の形態、及び第８の実施の形態の変形例への適用に際し留意する点は、第１の実施の形態に適用した場合に留意する点１、２に加え、３．ステップＳ１４３で代表色における色相角h*を算出するのと合わせて、彩度C*も算出することと、４．ステップＳ１４３で求められた代表色の色相角h*と彩度C*に基づいて、ステップＳ１４４では最適順応ファクタをDh*(h*)、Dh*'(h*)の代わりにDh*(h*,C*,Tc)とDh*'(h*,C*,Tc')を用いて算出することである。【０１８９】
−第１３の実施の形態−
第１、第６、第７、第８、第９の実施の形態は、それぞれ変換対象画素の色（色相等）もしくは、代表色の色（色相等）を基に、図２で示されるデータに基づき、順応ファクタを決定していた。第１０の実施の形態では、照明の色温度毎に、図２に示すような最適な順応ファクタDh*(h*)と色相角h*の関係を示すデータを予め取得して、変換対象画素もしくは代表色の色相角h*と入力画像取得時、もしくは入出力画像観察時の照明の色温度Tc,Tc'に基づいて順応ファクタを決定した。第１１の実施の形態では、異なる彩度(C*)毎に、図２に示すような最適な順応ファクタDh*(h*)と色相角h*の関係を示すデータを予め取得して、変換対象画素もしくは代表色の色相h*と彩度C*に基づいて順応ファクタを決定した。更に第１２の実施の形態では、照明の色温度毎に、最適な順応ファクタDh*(h*,C*)のデータを取得しておき、変換対象画素もしくは代表色の色相h*と彩度C*と、入力画像取得時もしくは入出力画像観察時の照明の色温度Tc,Tc'に基づいて順応ファクタを決定した。【０１９０】
第１３の実施の形態では、第１の実施の形態の図２に示すような最適な順応ファクタDh*(h*)と色相角h*の関係を示すデータを、視環境（例えば順応輝度LA、周囲の明るさ等）を変えて図２２の様に予め取得しておき、変換対象画素もしくは代表色の色（色相等）に加え、視環境も考慮する。【０１９１】
処理の流れを図２３を用いて説明する。【０１９２】
ステップＳ２１、２２は、図３と同様の処理であるので説明を省略する。【０１９３】
ステップＳ２１１では、ステップＳ２２で求めた撮影シーンの照明に関する情報を用いて、順応輝度LAと周囲の明るさで決まるパラメータF等を求める。ここでFは、YW0とYSW0の比を用いて、CIECAM02で定義されているテーブルから補間して求める。またLAは、LA=YW0/5から求める。尚、LAやFの算出は、別の方法で求めてもよい。例えばLAは、画像全体の平均輝度から推定してもよく、画像全体の平均輝度が白の輝度値に対してm[%]であれば、LA=m・YW0/100として算出してもよい。【０１９４】
ステップＳ２３は、図３と同様の処理であるので説明を省略する。【０１９５】
ステップＳ２１２では、ステップＳ２３で求めた画像観察時の視環境に関する情報を用いて、順応輝度LA'、周囲の明るさで決まるパラメータF'等を求める。ここでF'は、YW0'とYSW0'の比を用いて、CIECAM02で定義されているテーブルから補間して求める。またLA'は、LA'=YW0'/5から求める。尚、LA'やF'の算出は、別の方法で求めてもよい。例えばLA'は、画像全体の平均輝度から推定してもよく、画像全体の平均輝度が白の輝度値に対してm[%]であれば、LA'=m・YW0'/100として算出してもよい。【０１９６】
ステップＳ２4、Ｓ２５は、図３と同様の処理であるので、説明を省略する。【０１９７】
ステップＳ２１３では、予め色の見え実験結果から決定しておいた最適な順応ファクタと視環境、および色相の関係Dh*(LA,F,h*)から、ステップＳ２２からＳ２１２で算出した視環境のとき、ステップＳ２５で算出したh*に対して最適な色の見えを与える順応ファクタDh*(LA,F,h*)、Dh*'(LA',F',h*)を夫々算出する。【０１９８】
ステップＳ２７〜Ｓ３０は、図３と同様の処理であるので説明を省略する。【０１９９】
尚、第６から第１２の実施の形態に適用することも可能である。【０２００】
第６の実施の形態に適用する場合の留意点は、１．図１４（第６の実施の形態）のステップＳ８２の後に第１３の実施の形態の図２３のステップＳ２１１を追加する点、２．図１４（第６の実施の形態）のステップＳ２３の後に第１３の実施の形態の図２３のステップＳ２１２を追加する点、３．図１４（第６の実施の形態）のステップＳ２６を第１３の実施の形態の図２３のステップＳ２１３に変更する点である。【０２０１】
第７の実施の形態に適用する場合の留意点は、１．図１５（第７の実施の形態）のステップＳ１３２の後に第１３の実施の形態の図２３のステップＳ２１１を追加する点、２．図１５（第７の実施の形態）のステップＳ１３３の後に第１３の実施の形態の図２３のステップＳ２１２を追加する点、３．図１５（第７の実施の形態）のステップＳ２６を第１３の実施の形態の図２３のステップＳ２１３に変更する点である。【０２０２】
第８、９の実施の形態、及び第８の実施の形態の変形例に適用する場合の留意点は、１．図１６、図１７（第８の実施の形態、第８の実施の形態の変形例）のステップＳ２２の後に第１３の実施の形態の図２３のステップＳ２１１を追加する点、２．図１６、図１７（第８の実施の形態、第８の実施の形態の変形例）のステップＳ２３の後に第１３の実施の形態の図２３のステップＳ２１２を追加する点、３．図１６、図１７（第８の実施の形態、第８の実施の形態の変形例）のステップＳ１４４を第１３の実施の形態の図２３のステップＳ２１３に変更する点である。【０２０３】
第１０から１２の実施の形態に適用する場合の留意点は、第１の実施の形態に適用する場合（本実施の形態）と同様である。【０２０４】
−第１４の実施の形態−
第１３の実施の形態では、変換対象画素もしくは代表色の色（色相h*等）と視環境（例えば順応輝度LA、周囲の明るさ等）に依存して順応ファクタDh*を変化させた。【０２０５】
第１４の実施の形態では、h*に依存してDh*を変える場合に、全色相に渡って平均した順応ファクタの値は、入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境、及び出力画像観察時の視環境によって決まる値とすることで、視環境を変えた実験を多数行わなくても、ある一つの視環境で順応ファクタの色相依存関係を確かめておけばよくなり、見えの実験の負荷削減やパーソナルコンピュータやカメラのメモリ等の容量負荷削減が可能である。なお、全色相に渡って平均した順応ファクタの値は、入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境、及び出力画像観察時の視環境によって変化することになる。【０２０６】
CIECAM02の定義によれば、順応ファクタDは、順応輝度LAと周囲の明るさで決まるパラメータFから
【数６３】
2008050761.tif 000065
で求められる。【０２０７】
一方、LAを変えて第１の実施の形態で説明した様な見えの実験を行い、LAの色相依存を求めると、図２２の様に、LAによらず色相依存が見られるが、全色相で平均した値はLAによって変化し、上述のCIECAM02の定義から求められる値と近い値となる。【０２０８】
そこで本実施の形態では、まず第１の実施の形態で説明した様な見えの実験を、ある視環境X（順応輝度LA_X、周囲の明るさで決まるパラメータFX）下に於いて行い、視環境X下の順応ファクタの色相依存DX(h*)を図２の様に求める。次に、DX(h*)から順応輝度や周囲の明るさに依らない色相依存成分D0(h*)のみを抽出する。
【数６４】
2008050761.tif 000066
【数６５】
2008050761.tif 000067
【０２０９】
算出したD0(h*)は、予め、第１４の実施の形態で処理を行うパーソナルコンピュータ１の記憶媒体に保存しておく。【０２１０】
第１４の実施の形態におけるパーソナルコンピュータ１が実行する処理の流れを、図２４を用いて説明する。【０２１１】
ステップＳ２１、Ｓ２２は第１の実施の形態（図３）と同様の処理なので、説明を省略する。【０２１２】
ステップＳ２２１では、ステップＳ２２で求めた撮影シーンの照明に関する情報を用いて、順応輝度LAと周囲の明るさで決まるパラメータFを求める。ここでFは、YW0とYSW0の比を用いて、CIECAM02で定義されているテーブルから補間して求める。またLAは、LA=YW0/5から求める。尚、LAやFの算出は、別の方法で求めてもよい。例えばLAは、画像全体の平均輝度から推定してもよく、画像全体の平均輝度が白の輝度値に対してm[%]であれば、LA=m・YW0/100として算出してもよい。次に、LAとFから、D(LA,F)を算出する。
【数６６】
2008050761.tif 000068
【０２１３】
ステップＳ２３は、第１の実施の形態（図３）と同様の処理なので、説明を省略する。【０２１４】
ステップＳ２２２では、ステップＳ２３で求めた画像観察時の視環境に関する情報を用いて、順応輝度LA'、周囲の明るさで決まるパラメータF'を求める。ここでF'は、YW0'とYSW0'の比を用いて、CIECAM02で定義されているテーブルから補間して求める。またLA'は、LA'=YW0'/5から求める。尚、LA'やF'の算出は、別の方法で求めてもよい。例えばLA'は、画像全体の平均輝度から推定してもよく、画像全体の平均輝度が白の輝度値に対してm[%]であれば、LA'=m・YW0'/100として算出してもよい。次に、LA'とF'から、D(LA',F')を算出する。
【数６７】
2008050761.tif 000069
【０２１５】
ステップＳ２4〜Ｓ２５は第１の実施の形態（図３）と同様の処理なので、説明を省略する。【０２１６】
ステップＳ２２３では、
【数６８】
2008050761.tif 000070
【数６９】
2008050761.tif 000071
を算出する。【０２１７】
ステップＳ２７〜Ｓ３０は第１の実施の形態（図３）と同様なので説明を省略する。【０２１８】
（第１４の実施の形態の変形例１）
第１４の実施の形態では、第１の実施の形態に於けるDが、ただ色相に依存するだけではなく、全色相の平均Dが視環境で決まる値となる様にした。同様の処理を第８の実施の形態にも適用することができる。パーソナルコンピュータ１が実行する処理の流れを、図２５を用いて説明する。【０２１９】
ステップＳ２１、Ｓ２２、Ｓ２２１，Ｓ２３、Ｓ２２２は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。ステップＳ１４１〜ステップＳ１４３は第８の実施の形態の図１６と同様の処理なので、説明を省略する。【０２２０】
ステップＳ２２３は第１４の実施の形態の図２４のステップＳ２２３同様、代表色の色相角h*と、ステップＳ２２１、Ｓ２２２で算出したD(LA,F)、D(LA',F')を用いてDh*0(h*)、Dh*0'(h*)を算出する。【０２２１】
ステップＳ２４、Ｓ２７、Ｓ２８、Ｓ１４５、Ｓ１４６、Ｓ３０は第８の実施の形態の図１６と同様の処理なので、説明を省略する。即ち、第８の実施の形態に対し、ステップＳ２２の次にステップＳ２２１、ステップＳ２３の次にステップＳ２２２、ステップＳ１４４の代わりにステップＳ２２３に変更する点に留意すればよい。【０２２２】
尚、第９の実施の形態に適用する場合の留意点は同様である。【０２２３】
（第１４の実施の形態の変形例２）
第１４の実施の形態の変形例１では、第８の実施の形態に於けるDが、ただ色相に依存するだけではなく、全色相の平均Dが視環境で決まる値となる様にした。同様の処理を第８の実施の形態の変形例にも適用することができる。パーソナルコンピュータ１が実行する処理の流れを、図２６を用いて説明する。【０２２４】
ステップＳ２１、Ｓ２２、Ｓ２２１，Ｓ２３、Ｓ２２２は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。ステップＳ１５２、Ｓ１４３は第８の実施の形態の変形例の図１７と同様の処理なので、説明を省略する。【０２２５】
ステップＳ２２３は第１４の実施の形態の図２４のステップＳ２２３同様、代表色の色相角h*と、ステップＳ２２１、Ｓ２２２で算出したD(LA,F)、D(LA',F')を用いてDh*0(h*)、Dh*0'(h*)を算出する。【０２２６】
ステップＳ２４、Ｓ２７、Ｓ２８、Ｓ１４５、Ｓ３０は第８の実施の形態の変形例の図１７と同様の処理なので、説明を省略する。即ち、第８の実施の形態の変形例に対し、ステップＳ２２の次にステップＳ２２１、ステップＳ２３の次にステップＳ２２２、ステップＳ１４４の代わりにステップＳ２２３に変更する点に留意すればよい。【０２２７】
（第１４の実施の形態の変形例３）
第１４の実施の形態の変形例１では、第８の実施の形態に於けるDが、ただ色相に依存するだけではなく、全色相の平均Dが視環境で決まる値となる様にした。同様の処理を第６の実施の形態にも適用することができる。パーソナルコンピュータ１が実行する処理の流れを、図２７を用いて説明する。【０２２８】
ステップＳ８１、Ｓ８２は、第６の実施の形態の図１４と同様なので説明を省略する。ステップＳ２２１は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。ステップＳ２３は、第６の実施の形態の図１４と同様なので説明を省略する。ステップＳ２２２は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。ステップＳ２４、Ｓ２５は、第６の実施の形態の図１４と同様なので説明を省略する。【０２２９】
ステップＳ２２３は第１４の実施の形態の図２４のステップＳ２２３同様、変換対象画素の色相角h*と、ステップＳ２２１、Ｓ２２２で算出したD(LA,F)、D(LA',F')を用いてDh*0(h*)、Dh*0'(h*)を算出する。【０２３０】
ステップＳ２７〜Ｓ３０は第６の実施の形態の図１４と同様の処理なので、説明を省略する。即ち、第６の実施の形態に対し、ステップＳ８２の次にステップＳ２２１、ステップＳ２３の次にステップＳ２２２、ステップＳ２６の代わりにステップＳ２２３に変更する点に留意すればよい。【０２３１】
（第１４の実施の形態の変形例４）
第１４の実施の形態の変形例１では、第８の実施の形態に於けるDが、ただ色相に依存するだけではなく、全色相の平均Dが視環境で決まる値となる様にした。同様の処理を第７の実施の形態にも適用することができる。パーソナルコンピュータ１が実行する処理の流れを、図２８を用いて説明する。【０２３２】
ステップＳ２１、Ｓ１３２は、第７の実施の形態の図１５と同様なので説明を省略する。ステップＳ２２１は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。ステップＳ１３３は、第７の実施の形態の図１５と同様なので説明を省略する。ステップＳ２２２は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。ステップＳ２４、Ｓ２５は、第７の実施の形態の図１５と同様なので説明を省略する。【０２３３】
ステップＳ２２３は第１４の実施の形態の図２４のステップＳ２２３同様、変換対象画素の色相角h*と、ステップＳ２２１、Ｓ２２２で算出したD(LA,F)、D(LA',F')を用いてDh*0(h*)、Dh*0'(h*)を算出する。【０２３４】
ステップＳ２７〜Ｓ３０は第７の実施の形態の図１５と同様の処理なので、説明を省略する。即ち、第７の実施の形態に対し、ステップ１３２の次にステップＳ２２１、ステップＳ１３３の次にステップＳ２２２、ステップＳ２６の代わりにステップＳ２２３に変更する点に留意すればよい。【０２３５】
（第１４の実施の形態の変形例５）
第１１の実施の形態に適用する場合には、見えの実験では予めD0(h*,C*)を求め、順応輝度LA、周囲の明るさで決まるパラメータF、順応輝度LA'、周囲の明るさで決まるパラメータF'からD(LA,F)、D(LA',F')を算出し、変換対象画素や代表色の色相に応じて順応ファクタを算出する際、Dh*(h*,C*)、Dh*'(h*,C*)の代わりに、Dh*0(h*,C*)=D0(h*,C*)・D(LA,F)、Dh*0'(h*,C*)=D0(h*,C*)・D(LA',F')を算出することに留意すれば、上述の変形例と同様の手順で適用可能である。【０２３６】
（第１４の実施の形態の変形例６）
第１０の実施の形態に適用する場合には、見えの実験では予めD0(h*,Tc)を求め、順応輝度LA、周囲の明るさで決まるパラメータF、順応輝度LA'、周囲の明るさで決まるパラメータF'からD(LA,F)、D(LA',F')を算出し、変換対象画素や代表色の色相に応じて順応ファクタを算出する際、Dh*(h*,Tc)、Dh*'(h*,Tc')の代わりに、Dh*0(h*,Tc)=D0(h*,Tc)・D(LA,F)、Dh*0'(h*,Tc')=D0(h*,Tc')・D(LA',F')を算出することに留意すれば、上述の変形例と同様の手順で適用可能である。【０２３７】
（第１４の実施の形態の変形例７）
第１２の実施の形態に適用する場合には、見えの実験では予めD0(h*,C*,Tc)を求め、順応輝度LA、周囲の明るさで決まるパラメータF、順応輝度LA'、周囲の明るさで決まるパラメータF'からD(LA,F)、D(LA',F')を算出し、変換対象画素や代表色の色相に応じて順応ファクタを算出する際、Dh*(h*,C*,Tc)、Dh*'(h*,C*,Tc')の代わりに、Dh*0(h*,C*,Tc)=D0(h*,C*,Tc)・D(LA,F)、Dh*0'(h*,C*,Tc')=D0(h*,C*,Tc')・D(LA',F')を算出することに留意すれば、上述の変形例と同様の手順で適用可能である。【０２３８】
（第１４の実施の形態のその他の変形例）
尚、順応輝度LA、順応輝度LA'の算出では、全画像全体で均一の値としたが、画素毎に、若しくは領域毎に、算出してもよい。例えば変換対象画素周囲の所定数画素平均の相対輝度値や代表色の相対輝度値yを用いて、白色の輝度値Ywに対してLA=y・YW0/YW、LA'=y・YW0'/YW'として算出してもよい。第１４の実施の形態にこの変更を適用すると、パーソナルコンピュータ１が実行する処理の流れは、図２９の様になる。【０２３９】
ステップＳ２１、Ｓ２２は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。【０２４０】
ステップＳ２７１では、ステップＳ２２で求めた撮影シーンの照明に関する情報を用いて、周囲の明るさで決まるパラメータFを求める。Fの求め方は第１４の実施の形態のステップＳ２２１と同様であるので説明を省略する。ステップＳ２３は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。【０２４１】
ステップＳ２７２では、ステップＳ２３で求めた画像観察時の視環境に関する情報を用いて、周囲の明るさで決まるパラメータF'を求める。F'の求め方は第１４の実施の形態のステップＳ２２２と同様であるので説明を省略する。ステップＳ２４、Ｓ２５は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。【０２４２】
ステップＳ２７３では、変換対象画素周囲の所定数画素平均の相対輝度値yを用いて、白色の輝度値Ywに対してLA=y・YW0/YW、LA'=y・YW0'/YW'を算出し、ステップＳ２７１、Ｓ２７２で求めたF、F'を用いて、D(LA,F)、D(LA',F')を算出する。
【数７０】
2008050761.tif 000072
【数７１】
2008050761.tif 000073
【０２４３】
ステップＳ２２３、Ｓ２７〜Ｓ３０は第１４の実施の形態の図２４と同様の処理なので、説明を省略する。即ち、第１４の実施の形態に対し、ステップＳ２２１の代わりにステップＳ２７１、ステップＳ２２２の代わりにステップＳ２７２、ステップＳ２２３の前にステップＳ２７３を挿入する点に留意すればよい。【０２４４】
尚、順応輝度LA,LA'の求め方は、上述の手法以外でもよく、変換対象画素に応じた順応輝度を求められれば他の手法でもよい。【０２４５】
同様に、第１４の実施の形態の変形例１〜７の場合にも、変換対象画素周囲の所定数画素平均の相対輝度値や代表色の相対輝度値yを用いて順応輝度を算出してもよい。処理の流れは、ステップＳ２２１の代わりにステップＳ２７１、ステップＳ２２２の代わりにステップＳ２７２、ステップＳ２２３の前にステップＳ２７３を挿入する点に留意すればよいので、説明は省略する。【０２４６】
同様に、第１３の実施の形態の場合にも、変換対象画素周囲の所定数画素平均の相対輝度値yを用いて順応輝度を算出してもよい。処理の流れは、ステップＳ２１１の代わりに図２９のステップＳ２７１、Ｓ２１２の代わりに図２９のステップＳ２７２、ステップＳ２１３の前に図２９のステップＳ２７３のLA、LA'算出行程のみを行う点に留意すればよい。【０２４７】
−第１５の実施の形態−
第１４の実施の形態では、h*に依存して順応ファクタを変える場合に、全色相の平均順応ファクタ値は、入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境、及び出力画像観察時の視環境によって決まる値とした。この時、全色相で平均した順応ファクタの値D(LA,F)に対し、色相依存の順応ファクタをDh*0(h*)=D0(h*)・D(LA,F)、Dh*0'(h*)=D0(h*)・D(LA',F')で求めた。【０２４８】
第１５の実施の形態でも、全色相の平均順応ファクタ値は、入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境、及び出力画像観察時の視環境によって決まる値D(LA,F)とする。但し、色相依存の順応ファクタの求め方が異なるだけである。【０２４９】
本実施の形態でも、第１の実施の形態で説明した様な見えの実験を、ある視環境X（順応輝度LA_X、周囲の明るさで決まるパラメータFX）下に於いて行い、視環境X下の順応ファクタの色相依存DX(h*)を図２の様に求める。次に、DX(h*)から順応輝度や周囲の明るさに依らない色相依存成分D1(h*)のみを抽出する。
【数７２】
2008050761.tif 000074
【数７３】
2008050761.tif 000075
【０２５０】
算出したD1(h*)は、予め、第１５の実施の形態で処理を行うパーソナルコンピュータ１の記憶媒体に保存しておく。【０２５１】
第１５の実施の形態におけるパーソナルコンピュータ１が実行する処理の流れは、第１４の実施の形態と似ているので、第１４の実施の形態の処理の流れを説明した図２４を用いて説明する。第１５の実施の形態におけるパーソナルコンピュータ１が実行する処理の流れと、第１４の実施の形態におけるパーソナルコンピュータ１が実行する処理の流れの違いは、ステップＳ２２３のみである。【０２５２】
ステップＳ２２３では、Dh*0(h*)=D0(h*)・D(LA,F)、Dh*0'(h*)=D0(h*)・D(LA',F')の代わりに
【数７４】
2008050761.tif 000076
【数７５】
2008050761.tif 000077
を算出する。【０２５３】
その他の処理は第１４の実施の形態における図２４と同様なので説明を省略する。即ち、第１４の実施の形態に対し、予め見えの実験で算出する順応ファクタの色相依存成分を、割算より処理の軽い引き算に変更し、それに伴い、パーソナルコンピュータ１のステップＳ２２３で行う順応ファクタの算出方法を変更する点に留意すればよい。【０２５４】
また、第１４の実施の形態の変形例１〜７においても、見えの実験では予めD1(h*)を
【数７６】
2008050761.tif 000078
【数７７】
2008050761.tif 000079
から求めておく。そして、パーソナルコンピュータ１が実行する処理では、ステップＳ２２３において、Dh*0(h*)=D0(h*)・D(LA,F)、Dh*0'(h*)=D0(h*)・D(LA',F')の代わりに
【数７８】
2008050761.tif 000080
【数７９】
2008050761.tif 000081
としてもよい。D1(h*)の求め方や、ステップＳ２２３の処理の変更方法は、第１５の実施の形態と同様である。また、パーソナルコンピュータ１が実行するその他の処理は、第１４の実施の形態の変形例１〜７の夫々と同様なので説明を省略する。【０２５５】
また第１４の実施の形態のその他の変形例で説明した様に、画像全体で均一の順応輝度を算出する代わりに、変換対象画素周囲の所定数画素平均の相対輝度値や代表色の相対輝度値yを用いて順応輝度を算出してもよい。留意する点は第１４の実施の形態のその他の変形例で説明したことと同様であり、ステップＳ２２３においてDh*0(h*)=D0(h*)・D(LA,F)、Dh*0'(h*)=D0(h*)・D(LA',F')の代わりに
【数８０】
2008050761.tif 000082
【数８１】
2008050761.tif 000083
を算出する点に留意すればよいので、処理の詳しい説明を省略する。【０２５６】
−第１６の実施の形態−
第１４の実施の形態では、順応ファクタを色相に依存する成分D0(h*)と視環境に依存成分D(LA,F)を用いて、Dh*0(h*)=D0(h*)・D(LA,F)の形で算出した。第１５の実施の形態では、順応ファクタを色相に依存する成分D1(h*)と視環境に依存成分D(LA,F)を用いて、Dh*1(h*)=D1(h*)+D(LA,F)の形で算出した。第１４、１５の実施の形態では、ある一つの視環境に於ける順応ファクタの色相依存特性を調べることで、簡単に他の視環境に於ける順応ファクタの推定が可能である。【０２５７】
一方、第１の実施の形態で説明した様な見えの実験を、ある視環境X（順応輝度LA_X、周囲の明るさで決まるパラメータFX）下に於いて行い、視環境X下の順応ファクタの色相依存DX(h*)を図２の様に求める。次に、第１５の実施の形態の様に、DX(h*)から色相依存成分D2(h*)を抽出する。
【数８２】
2008050761.tif 000084
【数８３】
2008050761.tif 000085
【０２５８】
同様に、視環境Xとは異なる視環境Y（順応輝度LA_Y、周囲の明るさで決まるパラメータFY）下に於いて行い、視環境X下の順応ファクタの色相依存DY(h*)を図２の様に求める。次に、DY(h*)から色相依存成分D2'(h*)を抽出する。
【数８４】
2008050761.tif 000086
【数８５】
2008050761.tif 000087
【０２５９】
ここで、視環境X下で求めたD2(h*)と視環境Y下で求めたD2'(h*)が完全に一致していれば、第１５の実施の形態でも十分精度よく順応ファクタを算出できる。しかし図２２より、平均値を差し引いた色相依存成分（例えば振幅など）が、視環境により変化していることが判る。従って、正確にはD2(h*)とD2'(h*)は一致しない。【０２６０】
そこで、第１６の実施の形態では、複数の視環境に於ける順応ファクタの色相依存性を調べることで、順応ファクタをD(h*)=α(La,F)+β(La,F)×γ(h*)の形で求めることとし、より正確に他の視環境に於ける順応ファクタの色相依存特性を推定できる様にした。【０２６１】
まず上述の様に複数の視環境下に於いて見えの実験を行い、
【数８６】
2008050761.tif 000088
を求める。色相依存成分D3(LA,F,h*)のうち、振幅が視環境によって変化する分をD4(LA,F)として決定し、
【数８７】
2008050761.tif 000089
算出したD4(LA,F) 、D5(h*)は、予め、第１５の実施の形態で処理を行うパーソナルコンピュータ１の記憶媒体に保存しておく。【０２６２】
第１６の実施の形態におけるパーソナルコンピュータ１が実行する処理の流れは、第１４の実施の形態と似ているので、第１４の実施の形態の処理の流れを説明した図２４を用いて説明する。第１６の実施の形態におけるパーソナルコンピュータ１が実行する処理の流れと、第１４の実施の形態におけるパーソナルコンピュータ１が実行する処理の流れの違いは、ステップＳ２２３のみである。【０２６３】
ステップＳ２２３では、Dh*0(h*)=D0(h*)・D(LA,F)、Dh*0'(h*)=D0(h*)・D(LA',F')の代わりに
【数８８】
2008050761.tif 000090
【数８９】
2008050761.tif 000091
を算出する。【０２６４】
その他の処理は第１４の実施の形態における図２４と同様なので説明を省略する。即ち、第１４の実施の形態に対し、予め見えの実験で算出する順応ファクタの色相依存成分と視環境依存成分をより正確に反映できるように変更し、それに伴い、パーソナルコンピュータ１のステップＳ２２３で行う順応ファクタの算出方法を変更する点に留意すればよい。【０２６５】
また、第１４の実施の形態の変形例１〜７と同様の変形が可能である。即ち、パーソナルコンピュータ１が実行する処理では、ステップＳ２２３において、Dh*0(h*)=D0(h*)・D(LA,F)、Dh*0'(h*)=D0(h*)・D(LA',F')の代わりに
【数９０】
2008050761.tif 000092
【数９１】
2008050761.tif 000093
としてもよい。【０２６６】
D5(h*)の求め方や、ステップＳ２２３の処理の変更方法は、第１６の実施の形態と同様である。また、パーソナルコンピュータ１が実行するその他の処理は、第１４の実施の形態の変形例１〜７の夫々と同様なので説明を省略する。【０２６７】
また第１４の実施の形態のその他の変形例で説明した様に、画像全体で均一の順応輝度を算出する代わりに、変換対象画素周囲の所定数画素平均の相対輝度値や代表色の相対輝度値yを用いて順応輝度を算出してもよい。留意する点は第１４の実施の形態のその他の変形例で説明したことと同様であり、ステップＳ２２３においてDh*0(h*)=D0(h*)・D(LA,F)、Dh*0'(h*)=D0(h*)・D(LA',F')の代わりに
【数９２】
2008050761.tif 000094
【数９３】
2008050761.tif 000095
を算出する点に留意すればよいので、処理の詳しい説明を省略する。【０２６８】
−第１７の実施の形態−
第１３から１６の実施の形態では、視環境に応じて順応ファクタの色相依存を変化させた。しかし、視環境が暗い場合、色を認識する錐体の他に桿体という細胞が明暗の認識に関わっている。即ち、桿体の働くと言われている100Cd/m2以下では、桿体の応答が影響し、順応ファクタの色相依存が顕著でなくなる傾向にある。これは図２２にて、順応輝度が低い時に振幅が小さくなっていることから確認できる。【０２６９】
そこで、色相依存の影響が無視できる順応輝度の閾値La_limitを予め定め、それ以下の順応輝度の場合は、色相依存を考慮しなくてもよい。例えば、第１６の実施の形態に於いて、予め定めた閾値第La_limitに対して、
La＜La_limitの時、D4(LA,F)=0
La≧La_limitの時、D4(LA,F)≠0
としておけばよい。【０２７０】
尚、桿体の働きは、順応輝度Laだけではなく、周囲の明るさによっても変化する為、閾値La_limitは周囲の明るさにも応じて適切に設定する方が望ましい。例えば図３０の様に、周囲の明るさに応じて複数条件でLa_limitを決定しておき、観察視環境やシーンの照明条件に応じて、適切な条件に於けるLa_limitを選択したり、条件が近い２つのLa_limitから補間して適切なLa_limitを算出したりしてもよい。【０２７１】
尚、桿体の飽和する輝度が100Cd/m2程度であることから、La_limit≦100Cd/m2程度で設定するのが好ましい。【０２７２】
（第１７の実施の形態の変形例）
閾値La_limitは、視細胞に依存する為、観察者の特性によっても最適な閾値は変化する。従って、観察者毎に、色相依存の顕著な色（例えば黄色と青など）にて最適な順応ファクタ評価実験を、順応輝度を変えて複数行い、その結果に応じて、観察者毎に最適な閾値La_limitを設定してもよい。【０２７３】
また、閾値La_limitだけでなく、順応の程度も被験者によりバラツキは生じるので、より観察者に適した見えを提供する為には、第１４の実施の形態に於けるD0(h*)、第１５の実施の形態に於けるD1(h*)、第１６の実施の形態に於けるD4(LA,F)、D5(h*)、第１４から１６の実施の形態に於けるD(LA,F)等を観察者に応じて、調整してもよい。【０２７４】
尚、第１４から１７の実施の形態に於いて、CIECAM02で定義されるD(LA,F)を基準に色相依存成分を追加する形とした。定義では、視環境を変えると、0.6≦D(LA,F)≦1.0と変化する為、ある特定の視環境から他の視環境に於ける最適順応ファクタ値を推定する第１４から第１７の実施の形態等の場合、算出した順応ファクタが1より大きく、もしくは0より小さくなる可能性がある。その場合には、上述の方法で算出された順応ファクタが0から1の間に収まるよう処理を加える方がよい。【０２７５】
尚、第１４から１７の実施の形態に於いて、ある視環境X（順応輝度LA_X、周囲の明るさで決まるパラメータFX）下の順応ファクタの色相依存DX(h*)を求めたとき、最適順応ファクタの色相方向の平均；
【数９４】
2008050761.tif 000096
もしくは
【数９５】
2008050761.tif 000097
（Nは、色の見え実験に最適順応ファクタを評価した色の色相方向の分割数）
が、CIECAM02の定義から求められる順応ファクタ値
【数９６】
2008050761.tif 000098
と一致していれば、第１４から１７の実施の形態を適用して求めた他の視環境下に於ける順応ファクタDh*0(h*)等の色相方向の平均もCIECAM02の定義から求められる順応ファクタ値と一致する。従って、見えの実験で得られたDX(h*)の色相方向の平均がCIECAM02の定義から求められる順応ファクタ値から外れる場合は、実験誤差として、DX(h*)の平均を補正する処理を行ってもよい。【０２７６】
尚、実施例１４から１７の実施の形態、およびその変形例では、入力画像取得時のシーンの照明条件、もしくは入力画像観察時の視環境、および出力画像観察時の視環境によって決まる、全色相で平均した順応ファクタとして、CIECAM02で定義される
【数９７】
2008050761.tif 000099
【数９８】
2008050761.tif 000100
を用いた。しかし、入力画像取得時のシーンの照明条件もしくは入力画像観察時の視環境における順応輝度LA、および出力画像観察時の視環境における順応輝度LA'や、周囲の明るさのパラメータF、F'など視環境に関するパラメータで決まればよく、その他の関係式で算出してもよい。【０２７７】
また前述した様に、入力画像取得時のシーンの照明条件もしくは入力画像観察時の視環境における順応輝度LA、および出力画像観察時の視環境における順応輝度LA'は、それぞれ、シーンの照明輝度もしくは入力画像観察時の視環境における照明輝度YW0、および出力画像観察時の視環境における照明輝度YW0'に依存する。従って、D(LA,F)、D(LA',F')の代わりにYW0、YW0'の関数としてD(YW0,F)、D(YW0',F')を用いてもよい。【０２７８】
上記第１から第５の実施の形態、および、第８から第１７の実施の形態では、撮影シーンの照明に関する情報、および、画像観察時の視環境に関する情報を取得して、それぞれにおいて色によって順応ファクタを異ならせる場合を説明した。また、第６の実施の形態および第7の実施の形態では、異なる２つの視環境において、それぞれにおいて色によって順応ファクタを異ならせる場合を説明した。しかし、撮影シーンおよび画像観察時の視環境のいずれか一方に基づいてのみ、また、異なる２つの視環境のいずれか一方においてのみ、色によって順応ファクタを異ならせる処理を行ってもよい。【０２７９】
尚、上述の説明では、明度、彩度、色相としてCIELABのL*, C*, h*を用いた場合で説明したが、明度、彩度、色相に関するパラメータが得られれば、用いる色空間はCIELABに限らず、CIECAM02の相対的な色空間であるJ, C, hや絶対的な色空間であるQ, M, h等、その他の空間でもよい。【０２８０】
尚、第１の実施の形態、及び第６から第１７の実施の形態は、CIECAM02の様な色の見えモデル、白色点変換のみを考慮したvon Kriesの色順応変換など、色順応を考慮した多様な色変換方法に適用可能である。【０２８１】
−第１８の実施の形態−
第１の実施の形態では、不完全順応の度合いを表す順応ファクタを、変換対象色の色相に依存して変えることにより、色の見え再現をより精度よく行う例を示した。第１８の実施の形態では、色相依存の考慮の有無を変換対象画素の彩度により切り替えるようにしたことが特徴である。【０２８２】
（色の見え実験結果）
第１の実施の形態で説明した色の見えの実験を、ある視環境X（順応輝度LA_X、周囲の明るさで決まるパラメータFX）下に於いて行い、視環境X下の順応ファクタの色相依存DX(h*)を第１の実施の形態の図２の様に求める。次に、DX(h*)から順応輝度や周囲の明るさに依らない色相依存成分D1(h*)のみを抽出する。
【数９９】
2008050761.tif 000101
【数１００】
2008050761.tif 000102
【０２８３】
D1(h*)を用いて、任意の視環境（順応輝度LA、周囲の明るさで決まるパラメータF）下に於ける順応ファクタの色相依存性を
【数１０１】
2008050761.tif 000103
として推定することができる。尚、D2(LA,F)は、色相依存の振幅変化を複数の視環境に於ける見えの実験から決定しておく。【０２８４】
一方、彩度が低くなりグレー軸に近づくと、色相の区別がなくなる為、上述の様な色相依存が顕著でなくなる。そこで、予め、彩度を変えて図２のデータを取得しておき、色相依存が無視できる最大彩度C0*を決定しておく。色相依存が無視できるとは、図２の色相依存特性のグラフがフラットになるような状態、もしくは、グラフの振幅ΔDh*に対して、ある順応ファクタDで変換した時の色と、順応ファクタ(D+ΔDh*)で変換した時の色の差（例えばCIE LABのΔE）がある閾値以下となる状態を言う。【０２８５】
（変換処理）
次に、色変換処理について説明する。第１８の実施の形態に於ける処理は、第１の実施の形態の図１と同様にパーソナルコンピュータ１が処理を実行する。図３１は、第１８の実施の形態におけるパーソナルコンピュータ１が実行する処理のフローチャートを示す図である。予め上述の様な実験を行って、最適な順応ファクタと色相の関係Dh*(h*)をテーブルや近似式の形で画像処理を行うプログラム内に保有する。あるいは、プログラムが参照するデータとしてパーソナルコンピュータ１の記憶装置やメモリ等に保存しておく。【０２８６】
第１の実施の形態の図３の処理と同様なステップには同様なステップ番号を付す。まず、第１の実施の形態と同様に、ステップＳ２１にて、撮影画像を取得する。撮影画像は、複数の画素から構成され、カメラの分光感度に固有の色空間RGBで記述されている。【０２８７】
次にステップＳ２２に進む。ステップＳ２２では、撮影シーンの照明に関する情報を取得する。即ち、撮影シーンの照明白色点(XW,YW,ZW)と輝度YW0と周囲輝度YSW0である。例えば分光放射輝度計などを用いて白色点、及び輝度を測定した結果をユーザ（観察者）がキーボードから入力した値を読み込む。もしくはデジタルカメラに備わっているオートホワイトバランス機能を用いて、撮影時の照明白色点(XW,YW,ZW)を算出し、露出を自動調整する為に用いる測光機能の結果を解析して撮影時の照明輝度YW0を算出してもよい。周囲輝度YSW0は、測光範囲を広げて測光した結果を使用して算出してもよいし、撮影者が予め入力し画像に付加された情報やキーボードから入力した値を使用してもよい。【０２８８】
ステップＳ２２１では、ステップＳ２２で求めた撮影シーンの照明に関する情報を用いて、順応輝度LAと周囲の明るさで決まるパラメータFを求める。ここでFは、YW0とYSW0の比を用いて、CIECAM02で定義されているテーブルから補間して求める。またLAは、LA=YW0/5から求める。尚、LAやFの算出は、別の方法で求めてもよい。例えばLAは、画像全体の平均輝度から推定してもよく、画像全体の平均輝度が白の輝度値に対してm[%]であれば、LA=m・YW0/100として算出してもよい。次に、LAとFを式(2)に代入し、D(LA,F)を算出しておく。
【数１０２】
2008050761.tif 000104
【０２８９】
ステップＳ２３では、画像観察時の視環境に関する情報を取得する。標準視環境にてモニタ観察することを前提とすれば、sRGB標準視環境（白色点=D65、輝度Ydevice0'=80 Cd/m2、周囲照明白色点＝D50、周囲輝度YSW0'=4.1 Cd/m2）である。sRGB標準視環境と異なる視環境で観察する場合や画像をプリントした印刷物で観察する場合には、パーソナルコンピュータ１は、ユーザ（観察者）がキーボード等で入力した画像観察に使用するデバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'、および周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'を読み込み、デバイス白色点と周囲照明白色点の混合として画像観察視環境に於ける照明白色点(XW',YW',ZW')、輝度YW0'を決定する。【０２９０】
具体的には、画像観察に使用するデバイス光源と周囲照明光源の影響の比率Rmix'を用いて、次のように算出する。
【数１０３】
2008050761.tif 000105
【数１０４】
2008050761.tif 000106
【数１０５】
2008050761.tif 000107
【数１０６】
2008050761.tif 000108
【数１０７】
2008050761.tif 000109
【数１０８】
2008050761.tif 000110
【０２９１】
ステップＳ２２２では、ステップＳ２３で求めた画像観察時の視環境に関する情報を用いて、順応輝度LA'、周囲の明るさで決まるパラメータF'を求める。ここでF'は、YW0'とYSW0'の比を用いて、CIECAM02で定義されているテーブルから補間して求める。またLA'は、LA'=YW0'/5から求める。尚、LA'やF'の算出は、別の方法で求めてもよい。例えばLA'は、画像全体の平均輝度から推定してもよく、画像全体の平均輝度が白の輝度値に対してm[%]であれば、LA'=m・YW0'/100として算出してもよい。次に、LA'とF'を式(2)に代入し、D(LA',F')を算出しておく。
【数１０９】
2008050761.tif 000111
【０２９２】
次にステップＳ２４に進む。ステップＳ２４〜２８の処理は、全画素について処理が終了するまで繰り返す。カメラRGBはカメラの分光感度で決まる固有の色空間RGBで記述されているので、予め決まっているカメラRGB→XYZマトリックスMRGB→XYZを使用して、CIE 1931 XYZ色空間へ変換する。
【数１１０】
2008050761.tif 000112
【０２９３】
次にステップＳ２５に進み、変換対象画素の色相角h*とC*を算出する。色相角h*の算出は、CIE XYZからCIE LABに変換し、a*,b*の座標でなす角を0〜360°で算出した値を使用する。彩度はC*=(a*2+b*2)1/2から算出する。【０２９４】
ステップＳ３０１に進み、彩度C*が予め定めた色相依存の小さい範囲C*<C0*に入るか否かの判定を行う。【０２９５】
ステップＳ３０１で肯定された場合は、ステップＳ３０２に進み、撮影シーンの照明条件に於ける順応ファクタDと出力画像観察視環境に於ける順応ファクタD'を算出する。具体的には、Ｓ２２１で求めたLAとFからD(LA,F)を、Ｓ２２２で求めたLA'とF'からD(LA',F')を使用する。
【数１１１】
2008050761.tif 000113
【数１１２】
2008050761.tif 000114
【０２９６】
一方、ステップＳ３０１で否定された場合は、ステップＳ３０３に進み、予め見えの実験で求めておいた色相依存を考慮した順応ファクタの関係式を用いて、撮影シーンの照明条件に於ける順応ファクタDと出力画像観察視環境に於ける順応ファクタD'を算出する。
【数１１３】
2008050761.tif 000115
【数１１４】
2008050761.tif 000116
【０２９７】
ステップＳ２７では、ステップＳ２２、Ｓ２３で読み込んだ視環境の違いと、ステップＳ３０２、またはＳ３０３で算出したD、D'を使用して、撮影シーンの照明条件下に於ける色XYZと色の見えが近くなる観察照明条件下の色X'Y'Z'を算出する。ここで使用する色の見えモデルとしては、CIECAM02やvon Kries色順応式など視環境の違いを考慮した色順応変換式であればよい。CIECAM02を使用する場合は、XYZからJChへの変換では、定義式から算出される順応ファクタDの代わりに、ステップＳ３０２、またはＳ３０３で算出したDを使用し、JChからX'Y'Z'への変換ではD'を使用すればよい。Von Kriesの色順応式を使用する場合には、次式にて色変換を行う。
【数１１５】
2008050761.tif 000117
【０２９８】
但し、
【数１１６】
2008050761.tif 000118
【数１１７】
2008050761.tif 000119
【数１１８】
2008050761.tif 000120
【数１１９】
2008050761.tif 000121
【数１２０】
2008050761.tif 000122
【０２９９】
ステップＳ２８に進み、X'Y'Z'から出力画像の色空間R'G'B'に変換する。MXYZ→sRGBはXYZから出力画像の色空間RGBに変換するマトリックスであり、出力画像の色空間をsRGBとすれば、規格で決まっているマトリックスを使用すればよい。
【数１２１】
2008050761.tif 000123
【０３００】
ステップＳ２９では、全画素処理が終了していなければ、上記ステップＳ２４に戻し、処理を繰り返す。全画素終了したら、ステップＳ３０に進み画像を出力して終了する。【０３０１】
以上のようにして、画像撮影時のシーンの照明条件と、撮影画像を観察する照明条件が異なる場合に、精度良く近い色の見えを再現できる。上述の様に、本実施の形態の画像処理装置では、変換対象画素の色すなわち色相により最適な順応ファクタを使用して変換している為、正確な色の見え予測をすることができる。【０３０２】
−第１９の実施の形態−
第１８の実施の形態では、不完全順応の度合いを表す順応ファクタを、変換対象色の色相に依存して変えることにより、色の見え再現をより精度よく行うが、このとき、色相依存の考慮の有無を変換対象画素の彩度により切り替えた。第１９の実施の形態では、順応輝度によって色相依存の有無を切り替える。【０３０３】
視環境が暗い場合、色を認識する錐体の他に桿体という細胞が明暗の認識に関わっている。即ち、桿体の働くと言われている100Cd/m2以下では、桿体の応答が影響し、順応ファクタの色相依存が顕著でなくなる傾向にある。そこで、本実施の形態では、色相依存の影響が無視できる順応輝度の閾値LA0を予め定め、それ以下の順応輝度の場合は、色相依存を考慮しないようにする。【０３０４】
（色の見え実験）
第１８の実施の形態と同様の見えの実験を行い、任意の視環境（順応輝度LA、周囲の明るさで決まるパラメータF）下に於ける順応ファクタの色相依存性を
【数１２２】
2008050761.tif 000124
として決定しておく。【０３０５】
一方、順応輝度を下げると、図３２の様に色相依存による順応ファクタの変化が小さくなる。従って予め色相依存の変化が、ある程度以下（例えば、一般に色の差が判別できると言われる色差ΔE=３以下と言われるので、色相依存で取り得る順応ファクタの最大値と最小値に対する色の差ΔE<3など）では順応ファクタの色相依存を無視しても見えには影響を与えない。そこで、この閾値LA0を予め上述の見えの実験から決定しておく。若しくは、明暗の認識のみを行う桿体が働いて認識している様な順応輝度（約100 Cd/m2以下）では色の認識が弱くなるので、色相依存の影響も小さくなる。色の認識を行う錐体は約1 Cd/m2以上から有効であることから、1 Cd/m2≦LA0≦100 Cd/m2の範囲で適当な値をLA0として設定すればよい。【０３０６】
（変換処理）
第１９の実施の形態に於ける処理は、第１の実施の形態の図１と同様にパーソナルコンピュータ１が処理を実行する。【０３０７】
次に、色変換処理について説明する。パーソナルコンピュータ１が行う変換処理は第１８の実施の形態と同様なので、図３１を用いて説明する。【０３０８】
ステップＳ２１からＳ２５は、第１８の実施の形態の図３１と同様である。ステップＳ３０１に進み、順応輝度LAが予め定めた色相依存の小さい範囲LA<LA0に入るか否かの判定を行う。ステップＳ３０１で肯定された場合は、ステップＳ３０２に進む。ステップＳ３０２の処理は第１８の実施の形態の図３１と同様である。一方、ステップＳ３０１で否定された場合は、ステップＳ３０３に進む。ステップＳ３０３の処理は第１の実施の形態の図３と同様である。ステップＳ２７からＳ３０の処理は第１８の実施の形態の図３１と同様である。【０３０９】
なお、上記Dh*(LA,F,h*)の式(24)は、D(h*)=α(La,F)+β(La,F)×γ(h*)の形で記載することができる。β(La,F)は、式(24)のD2(LA,F)に相当し、色相依存の振幅変化を複数の視環境に於ける見えの実験から決定しておくものであるが、ステップＳ３０２の処理では、β(La,F)=0とする。【０３１０】
尚、桿体の働きは、順応輝度Laだけではなく、周囲の明るさによっても変化する為、閾値LA0は周囲の明るさにも応じて適切に設定する方が望ましい。例えば図３３の様に、周囲の明るさに応じて複数条件でLA0を決定しておき、観察視環境やシーンの照明条件に応じて、適切な条件に於けるLA0を選択したり、条件が近い２つのLA0から補間して適切なLA0を算出したりしてもよい。【０３１１】
−第２０の実施の形態−
第１８の実施の形態では、不完全順応の度合いを表す順応ファクタを、変換対象色の色相に依存して変えることにより、色の見え再現をより精度よく行うが、このとき、色相依存の考慮の有無を変換対象画素の彩度により切り替えた。第１９の実施の形態では、順応輝度によって色相依存の有無を切り替えた。第２０の実施の形態では、撮影シーンの照明、もしくは入出力画像観察時の視環境の照明の色温度により、色相依存の有無を切り替える。【０３１２】
（色の見え実験）
第１８の実施の形態と同様の見えの実験を行い、任意の視環境（順応輝度LA、周囲の明るさで決まるパラメータF）下に於ける順応ファクタの色相依存性を
【数１２３】
2008050761.tif 000125
として決定しておく。【０３１３】
一方、順応ファクタを考慮した白
【数１２４】
2008050761.tif 000126
に於いて、照明の白色点が等エネルギー白色に近い場合、
【数１２５】
2008050761.tif 000127
となる為、Dの値によらず順応を考慮した白(LW,D,MW,D,SW,D)の値は変化しにくくなる。【０３１４】
そこで、Dを変えても色の変化に与える影響が小さい照明条件の範囲を予め決定しておく。決定方法は、色の見え実験から実際に色相依存が顕著でなくなる範囲を決定したり、順応ファクタDの取り得る範囲（0.0から1.0）に於いて色の変化（色差ΔE）が、無視しても見えに影響を与えない色温度の範囲とするなどである。また照明条件の範囲とは、例えば色温度Tに対し、T1<T<T2と範囲と決定したり、照明白色点の色度座標(xw,yw)の範囲を２次元で決定しておけばよい。なお、色相依存の小さい色温度の下限値T1は4500Kで、色相依存の小さい色温度の上限値T2は7000Kであることが好ましく、T1は5000KでT2は6000Kに設定することがより好ましい。【０３１５】
（変換処理）
第２０の実施の形態に於ける処理は、第１の実施の形態の図１と同様にパーソナルコンピュータ１が処理を実行する。【０３１６】
次に、色変換処理について説明する。パーソナルコンピュータ１が行う変換処理は第１８の実施の形態と同様なので、図３１を用いて説明する。【０３１７】
ステップＳ２１からＳ２５は、第１８の実施の形態の図３１と同様である。第１８の実施の形態に加えて留意する点は、ステップＳ２２で撮影シーンの照明に関する情報取得で取得された照明の白色点(Xw,Yw,Zw)から色温度(T)を算出する工程を追加することと、ステップＳ２３で画像観察時の視環境に関する情報取得で取得された照明に関する白色点(Xw',Yw',Zw')から色温度(T')を算出する工程が追加されることである。【０３１８】
次に、ステップＳ３０１に進み、ステップＳ２２、Ｓ２３で算出した色温度T,T'が、予め定めた色相依存の小さい範囲T1<(T,T')<T2に入るか否かの判定を行う。ステップＳ３０１で肯定された場合は、ステップＳ３０２に進む。ステップＳ３０２の処理は第１の実施の形態の図３と同様である。一方、ステップＳ３０１で否定された場合は、ステップＳ３０３に進む。ステップＳ３０３の処理は第１８の実施の形態の図３１と同様である。ステップＳ２７からＳ３０の処理は第１８の実施の形態の図３１と同様である。【０３１９】
ところで、白色点から色温度の求め方は、例えば、黒体輻射されるスペクトルに該当する座標をxy色度座標に表記し、それぞれの白色点と最も近い位置に座標位置がある当該スペクトルに対応する色温度を選び出す。これを色温度とする。【０３２０】
−変形例−
尚、第１８から第２０の実施の形態は、第１の実施の形態などと同様に、CIECAM02の様な色の見えモデルや白色点変換のみを考慮したvon Kriesの色順応変換など、色順応を考慮した多様な色変換方法に適用可能である。【０３２１】
尚、第１８の実施の形態のステップＳ２３では、画像観察視環境に於ける照明白色点として、CIE 1931 XYZ表色系の値(XW',YW',ZW')を算出している。しかし式(8)を省略して、錐体応答の空間に於ける値(LW',MW',SW')を、画像観察視環境に於ける照明白色点として取得してもよい。この場合、ステップＳ２７の式(22)も省略することができ、ステップＳ２３で取得した(LW',MW',SW')をそのまま式(20)に使用することができる。従って処理を軽くすることが出来る。【０３２２】
上記第１８の実施の形態では、モニタ６に画像表示する場合を例に説明したが、出力画像をプロジェクタやプリントで表示してもよい。プリントで表示する場合の入力データの例はステップＳ２３の説明に有る通りである。プロジェクタに表示する場合、出力画像観察時の視環境は、ユーザがキーボード等で入力した値を使用する。デバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'は、スペックに書かれているデータを用いるか、ユーザが実測値を入力する。周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'は、暗室すなわち暗い部屋で観察する場合はYSW0'=0、(XSW',YSW',ZSW')=(Xdevice',Ydevice',Zdevice')を入力する。暗室でない場合は、ユーザが周囲視環境を実測し、実測値を入力する。【０３２３】
また、色相、明度、彩度としてCIE LABのh*,L*,C*を使用する例で説明したが、CIE LAB以外の色空間に於けるパラメータを使用してもよい。例えば、CIECAM02の相対的な明度、彩度、色相のパラメータであるJ,C,h等を使用してもよい。【０３２４】
尚、第１８から第２０の実施の形態を組み合わせて、彩度、順応輝度、照明の色温度の３つの条件の組み合わせで、順応ファクタの色相依存の考慮有無を決定してもよい。【０３２５】
尚、第１８から第２０の実施の形態では、閾値を複数被験者の実験結果から決定しているが、実際には閾値は被験者によってばらつく。そこで、より画像観察者の見えに忠実に再現する為には、上述の閾値C0*、LA0、T1、T2などを観察者毎に最適化してもよい。【０３２６】
例えば、第１９の実施の形態の順応輝度の閾値を観察者特性に応じて調整する場合は、予め平均的な値として定めた順応輝度の閾値LA0の前後で複数変えた条件下で、
【数１２６】
2008050761.tif 000128
【数１２７】
2008050761.tif 000129
のDを用いて変換した場合に、どちらの色の見えが近いかの評価をしてもらう。このとき、全ての色に付いて行う必要はなく、色相依存の影響が顕著な色（黄色と青など）だけでもよい。【０３２７】
色相依存の顕著な色に付いて、Dh*(LA,F,h*)=D(LA,F)による変換でも、見えが近いと選ばれる順応輝度の閾値LA0'を決定し、観察者固有の閾値としてパーソナルコンピュータのメモリ等に記憶しておく。実際の処理の流れは、第１９の実施の形態と同様であり、ステップS３０１での判定ではLA0の代わりにLA0'を使用する点に留意すればよい。【０３２８】
（第１８〜第２０の実施の形態の変形例１）
上記実施の形態では、撮影した画像をモニタやプリントで観察する際、撮影シーンの照明条件と画像観察時の視環境の違いを考慮して色変換する場合を例に説明した。しかし、上記実施の形態で説明した内容は、視環境の異なる条件下の色の見えを色の見えモデル（色順応式も含む）を用いて予測する場合であれば適用可能である。【０３２９】
例えば上記実施の形態の他に、第１の視環境（例えばモニタ表示）にて観察する為に作成された画像を第２の視環境（印刷したハードコピーの観察）にて観察する場合に、両者の見えが近くなる様に色の見えを考慮して変換する場合や、撮影時に予め標準視環境で観察されることを前提に色順応変換されていた画像を実際の視環境に合わせて更に色順応変換する場合などにも適用できる。【０３３０】
この場合、Ｓ２２では、撮影シーンの照明に関する情報の代わりに、取得した画像データから第１の画像データの作成に用いられた第１の視環境に関する情報を取得する。第１の視環境に関する情報は、入力画像データ作成時に使用された、入力画像データ観察に使用する機器や周囲観察環境の照明に関する条件である。画像データに第１の視環境に関するデータが無い場合は、sRGB標準視環境（白色点=D65、輝度Ydevice0=80 Cd/m2、周囲照明白色点=D50、周囲輝度YSW0=4.1 Cd/m2）を読み込み、デバイス白色点と周囲照明白色点の混合として第１の視環境に於ける照明白色点(XW,YW,ZW)、輝度YW0を算出する。【０３３１】
sRGB標準視環境と異なる視環境で観察するために作成された画像の場合には、ユーザ（観察者）がキーボードで入力した画像作成に使用したデバイスの白色点(Xdevice,Ydevice,Zdevice)、輝度Ydevice0、および周囲照明白色点(XSW,YSW,ZSW)、周囲照明輝度YSW0を読み込み、デバイス白色点と周囲照明白色点の混合として第１の画像観察視環境に於ける照明白色点(XW,YW,ZW)、輝度YW0を算出してもよい。【０３３２】
Ｓ２３では、第２の視環境に関する情報を取得する。第２の視環境に関する情報とは、実際に出力画像データ観察に使用する機器や周囲観察環境の照明に関する条件である。これら、第１の視環境に関する情報や第２の視環境に関する情報は、照明光源の輝度や白色点、周囲の明るさなど、見えに影響を与える視野周辺の環境に関する情報である。【０３３３】
他の処理は、第１８から第２０の実施の形態の図３１と同様なので説明を省略する。【０３３４】
（第１８〜第２０の実施の形態の変形例２）
第１の視環境で観察している画像を別の第２の視環境で観察する場合について説明する。例えば、モニタ観察画像をプリント出力する場合やプロジェクタで投影して観察する場合に、モニタ観察画像と見えを合わせる場合である。【０３３５】
この場合、Ｓ２２では、撮影シーンの照明に関する情報の代わりに、ユーザがキーボード等で入力した第１の画像観察時の視環境に関する情報を取得する。「CRT観察（sRGB標準視環境）」が選択されると、デバイスの白色点(Xdevice,Ydevice,Zdevice)としてD65のXYZ値、輝度Ydevice0=80 Cd/m2、周囲照明白色点(XSW,YSW,ZSW)としてD50のXYZ値、周囲照明輝度YSW0=4.1 Cd/m2を設定する。もしくは、ユーザが入力した(Xdevice,Ydevice,Zdevice)、Ydevice0、(XSW,YSW,ZSW)、YSW0の各値を取得する。【０３３６】
取得した前記４つの値から、入力画像観察時の視環境に於いて順応している照明白色点の錐体応答(LW,MW,SW)と輝度YW0を決定する。【０３３７】
Ｓ２３では、ユーザがキーボード等で入力した出力画像観察時の視環境に関する情報を取得する。即ちユーザ（観察者）がキーボード等で入力した出力画像観察に使用するデバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'、および周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'を読み込み、デバイス白色点と周囲照明白色点の混合として出力画像観察時の視環境に於いて順応している照明白色点の錐体応答(LW',MW',SW')、輝度YW0'を算出する。【０３３８】
例えば、出力画像をプリントした印刷物で観察する場合には、ユーザ（観察者）がキーボード等で入力した出力画像観察に使用する照明の白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'、および周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'を使用する。出力画像をプロジェクタに表示する場合は、ユーザがキーボード等で入力した以下の値を使用する。【０３３９】
デバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'は、スペックに書かれているデータを用いるか、ユーザ実測値を入力する。周囲照明白色点(XSW',YSW',ZSW')、周囲照明輝度YSW0'は、暗室で観察する場合はYSW0'=0、(XSW',YSW',ZSW')=(Xdevice',Ydevice',Zdevice')を入力する。暗室でない場合は、ユーザが周囲視環境を実測し、実測値を入力する。【０３４０】
(LW',MW',SW')、YW0'の算出は、第１８の実施の形態同様、出力画像観察に使用するデバイス光源と周囲照明光源の影響の比率Rmix'を用いて算出する。【０３４１】
（第１８〜第２０の実施の形態の変形例３）
第１８から第２０の実施の形態は、各画素の色変換を行う際、各画素の色（色相、彩度等）に応じて順応ファクタDを変化させる例で説明した。しかし、画素毎に変化させるのではなく、画像データにおいて複数の近傍画素を一つの領域として捉え、その領域毎に代表色を決定し、代表色の色（色相、彩度等）によってその領域内の画素の変換に適応する順応ファクタDを決定してもよい。すなわち、画像データの色（色相、彩度等）の分布を解析し、解析した色の分布に基づいて順応ファクタを決定してもよい。【０３４２】
領域分割の方法としては、エッジ抽出結果により被写体毎に領域分割する場合や、画像をブロックに分割する方法などが考えられる。エッジ抽出結果により被写体毎に領域分割している場合、Dの変化と被写体の輪郭は一致する為、Dの変化による境界が目立つなど不自然さが目立たない。【０３４３】
ブロック分割してブロック領域毎にDを変える場合、領域の境界では使用する順応ファクタに不連続が生じる。この様に分割領域境界近傍においてDの変化による不自然さが生じる場合、境界領域の複数画素についてはDが連続につながるよう隣接する領域の代表順応ファクタDから補間して求めてもよい。【０３４４】
（第１８〜第２０の実施の形態の変形例４）
撮影画像をデジタルカメラの背面液晶で表示する場合を例に説明する。本変形例のデジタルカメラは、第９の実施の形態のデジタルカメラ３０と同様な構成であるので、第９の実施の形態の図１８を参照して説明する。デジタルカメラ３０は、制御部３１、撮像部３２、メモリ３３、モニタ３４、測光部３５などから構成される。モニタ３４は、デジタルカメラ３０の背面に設けられた液晶モニタである。【０３４５】
制御部３１は、マイクロプロセッサおよび周辺回路などから構成され、メモリ３３内に格納されたプログラムを実行することによりデジタルカメラ３０の各種の制御および以下で説明する画像処理を行う。撮像部３２は、ＣＣＤなどの撮像素子（イメージセンサ）から構成され、撮影レンズなどの撮影光学系（不図示）を介して入射する被写体からの光束を受光して電気信号に変換し、変換された電気信号をデジタル変換して画像データとして出力する。撮像部３２は、ＲＧＢカラーフィルターを有し、ＲＧＢ表色系の画像データを出力する。【０３４６】
デジタルカメラ３０の制御部３１が実行する画像処理は図３１と同様なので、図３１を用いて説明する。【０３４７】
ステップＳ２１では、撮像部３２により撮影された撮影画像を取得する。Ｓ２２１は第１の実施の形態の図３１と同様である。ステップＳ２２では撮影シーンに関する情報を取得する。詳細は第１８の実施の形態と同様なので説明を省略する。ステップＳ２３では画像観察時の視環境に関する情報を取得する。Ｓ２２２は第１８の実施の形態の図３１と同様である。【０３４８】
ここで画像観察デバイスはデジタルカメラ３０のモニタ３４である。従って、予めメモリ３３に記憶されたモニタ３４の白色点、並びに輝度を、デバイスの白色点(Xdevice',Ydevice',Zdevice')、輝度Ydevice0'として夫々読み込む。また、測光部３５による測光機能を利用して周囲輝度を算出し、周囲照明輝度YSW0'として設定する。周囲の白色点(XSW',YSW',ZSW')は、直前に撮影した画像のWBデータから光源を決定し、光源の白色点を読み込む。【０３４９】
ユーザが「太陽光、曇天、電球、蛍光灯」などのメニューから観察時の光源を選択した場合は、選択した光源に対応する白色点を周囲照明の白色点として読み込む。デバイス白色点と周囲照明白色点の混合として画像観察視環境に於ける照明白色点(XW',YW',ZW')、輝度YW0'を算出する方法は、第１８の実施の形態と同様なので説明を省略する。【０３５０】
ステップＳ２４〜ステップＳ３０は第１８の実施の形態と同様である。ステップＳ３０で出力した画像をモニタ３４に表示すれば、視環境に関する情報に関してユーザが入力を省略しても、簡単にデジタルカメラ３０の背面に設けられたモニタ３４で実シーンと見えの近い画像が表示される。【０３５１】
（第１８〜第２０の実施の形態の変形例５）
順応輝度LA、順応輝度LA'の算出では、全画像全体で均一の値としたが、画素毎に、若しくは領域毎に、算出してもよい。例えば変換対象画素周囲の所定数画素平均の相対輝度値や代表色の相対輝度値yを用いて、白色の輝度値Ywに対してLA=y・YW0/ YW、LA'=y・YW0'/ YW'として算出してもよい。第１８の実施の形態にこの変更を適用すると、パーソナルコンピュータ１が実行する処理の流れは、図３４の様になる。【０３５２】
ステップＳ２１、Ｓ２２は第１８の実施の形態の図３１と同様の処理なので、説明を省略する。ステップＳ３１１では、ステップＳ２２で求めた撮影シーンの照明に関する情報を用いて、周囲の明るさで決まるパラメータFを求める。Fの求め方は第１８の実施の形態のステップＳ２２１と同様であるので説明を省略する。ステップＳ２３は第１８の実施の形態の図３１と同様の処理なので、説明を省略する。【０３５３】
ステップＳ３１２では、ステップＳ２３で求めた画像観察時の視環境に関する情報を用いて、周囲の明るさで決まるパラメータF'を求める。F'の求め方は第１８の実施の形態のステップＳ２２２と同様であるので説明を省略する。ステップＳ２４、Ｓ２５は第１８の実施の形態の図３１と同様の処理なので、説明を省略する。ステップＳ３１３では、変換対象画素周囲の所定数画素平均の相対輝度値yを用いて、白色の輝度値Ywに対してLA=y・YW0/ YW、LA'=y・YW0'/ YW'を算出し、ステップＳ３１１、Ｓ３１２で求めたF、F'を用いて、D(LA,F)、D(LA',F')を算出する。
【数１２８】
2008050761.tif 000130
【数１２９】
2008050761.tif 000131
【０３５４】
ステップＳ３０１〜Ｓ３０は第１８の実施の形態の図３１と同様の処理なので、説明を省略する。即ち、第１８の実施の形態に対し、ステップＳ２２１の代わりにステップＳ３１１、ステップＳ２２２の代わりにステップＳ３１２、ステップＳ２５の前にステップＳ３１３を挿入する点に留意すればよい。同様に、第１９から第２０の実施の形態にも適用できる。【０３５５】
尚、順応輝度LA、LA'の求め方は、上述の手法以外でもよく、変換対象画素に応じた順応輝度を求められれば他の手法でもよい。【０３５６】
上記第１８の実施の形態では、撮影シーンの照明に関する情報、および、画像観察時の視環境に関する情報を取得して、それぞれにおいて色によって順応ファクタを異ならせる場合を説明した。しかし、撮影シーンおよび画像観察時の視環境のいずれか一方に基づいてのみ、また、異なる２つの視環境のいずれか一方においてのみ、色によって順応ファクタを異ならせる処理を行ってもよい。【０３５７】
尚、第１８〜第２０の実施の形態の説明では、明度、彩度、色相としてCIELABのL*,C*,h*を用いた場合で説明したが、明度、彩度、色相に関するパラメータが得られれば、第１〜第１７の実施の形態と同様に、用いる色空間はCIELABに限らず、CIECAM02の相対的な色空間であるJ,C,hや絶対的な色空間であるQ,M,h等、その他の空間でもよい。【０３５８】
尚、第１８から第２０の実施の形態では、順応ファクタは変換対象画素、若しくは変換対象画素と周辺画素から算出した代表色の色相、及び視環境（順応輝度、周囲の明るさ）のみに依存するとし、視環境や画像の色（色相）により、順応ファクタの色相依存の有無を切り替えた。しかし、その他の影響（例えば彩度依存性、明度依存性やその他の視環境など）を考慮した場合にも、視環境により影響の大小が異なる場合には、閾値を設定して、その閾値以下では画像の色による依存性を無視することで処理を軽くすることが出来る。【０３５９】
観察対象がディスプレイの様な自発光の色である光源色モードか、反射光を見ている物体色モードかといった見えのモードの違いで、見えが変わる場合がある。第１から第２０の実施の形態では、順応ファクタは変換対象画素や代表色の色相や彩度、観察環境の順応輝度、色温度、明るさ等で変化させたが、前記の様な見えのモードの違いで使用する順応ファクタを変化させてもよい。【０３６０】
第１の実施の形態に適応する場合には、色の見え実験を、ディスプレイの様な自発光の色票に対して（光源色モード）と、印刷された色票に対して（物体色モード）の最低２種類の見えのモードに対して行い、図２の様な最適順応ファクタDh*と色相角h*の関係を示すデータを、光源色モード、物体色モードの最低2種類のモードに対して用意しておく。そして、図３のＳ２３の画像表示条件がディスプレイに表示する例の様に光源色に相当する場合は、光源色モードに対するDh*を使用し、図３のＳ２３の画像表示条件が印刷物観察の場合は物体色モードに対するDh*を使用するように、入出力環境の見えのモードに応じて使用する順応ファクタを切り替えればよい。【０３６１】
また、図３のＳ２２では、通常のシーンは反射光を見ている場合が多いので物体色モードを使用し、光源色に相当する被写体を撮影した場合はユーザー入力により光源色モードに対するDh*を使用する様に切り替えればよい。他の実施の形態でも同様に、入出力視環境に於ける対象被写体の見えのモードに応じて、使用する順応ファクタを切り替えて使用すると尚良い。顕微鏡の試料観察の様に透過光観察の場合には、見えのモードとして光源色モードを使用しても良いし、透過光を対象とした透過光色モードを作成しておき、こちらを使用しても良い。【０３６２】
上記では、種々の実施の形態および変形例を説明したが、本発明はこれらの内容に限定されるものではない。本発明の技術的思想の範囲内で考えられるその他の態様も本発明の範囲内に含まれる。【０３６３】
次の優先権基礎出願の開示内容は引用文としてここに組み込まれる。
日本国特許出願２００６年第２８７９１８号（２００６年１０月２３日出願）
日本国特許出願２００６年第３１７３９１号（２００６年１１月２４日出願）
日本国特許出願２００６年第３１７７４８号（２００６年１１月２４日出願）
【図１】
2008050761.tif 000132
【図２】
2008050761.tif 000133
【図３】
2008050761.tif 000134
【図４】
2008050761.tif 000135
【図５】
2008050761.tif 000136
【図６】
2008050761.tif 000137
【図７】
2008050761.tif 000138
【図８】
2008050761.tif 000139
【図９】
2008050761.tif 000140
【図１０】
2008050761.tif 000141
【図１１】
2008050761.tif 000142
【図１２】
2008050761.tif 000143
【図１３】
2008050761.tif 000144
【図１４】
2008050761.tif 000145
【図１５】
2008050761.tif 000146
【図１６】
2008050761.tif 000147
【図１７】
2008050761.tif 000148
【図１８】
2008050761.tif 000149
【図１９】
2008050761.tif 000150
【図２０】
2008050761.tif 000151
【図２１】
2008050761.tif 000152
【図２２】
2008050761.tif 000153
【図２３】
2008050761.tif 000154
【図２４】
2008050761.tif 000155
【図２５】
2008050761.tif 000156
【図２６】
2008050761.tif 000157
【図２７】
2008050761.tif 000158
【図２８】
2008050761.tif 000159
【図２９】
2008050761.tif 000160
【図３０】
2008050761.tif 000161
【図３１】
2008050761.tif 000162
【図３２】
2008050761.tif 000163
【図３３】
2008050761.tif 000164
【図３４】
2008050761.tif 000165
【国際調査報告】
2008050761.tif 000166
2008050761.tif 000167
2008050761.tif 000168
2008050761.tif 000169