(57)【特許請求の範囲】
【請求項１】
一連のコマ画像に対してコマ画像単位にその構成画素に対するカラー補正量を求める補正量演算手段と、三原色のカラー成分を各座標軸とする色空間を所定間隔で画定して複数の空間領域を構成し、任意のコマ画像の構成画素の属する空間領域に当該コマ画像に対する前記カラー補正量を割り付ける処理を全コマ画像に対して実行する補正量割付手段と、空間領域毎に割り付けられた前記カラー補正量に基づいて領域に対する補正量を求める領域補正量演算手段と、特定のコマ画像の構成画素の属する空間領域から抽出した前記領域補正量を当該特定のコマ画像の構成要素に対する補正量、或いは補正量算出の要因として求める揃い補正量演算手段とを備え、前記揃い補正量演算手段は、各空間領域に属する画素数が当該特定のコマ画像の所定比率以上となる空間領域から前記領域補正量を抽出するように構成される写真画像処理装置。
【請求項２】
前記揃い補正量演算手段により求められた揃い補正量と、対応するコマ画像のカラー補正量とを所定の重み係数で融合した適正カラー補正量を求める適正補正量演算手段を備えて構成される請求項１記載の写真画像処理装置。
【請求項３】
前記補正量割付手段は、一連のコマ画像の空間領域への分布状態に基づいて領域サイズを調整する領域サイズ調整手段を備えている請求項１または２記載の写真画像処理装置。
【請求項４】
一連のコマ画像に対してコマ画像単位にカラー補正量を求める補正量演算ステップと、三原色のカラー成分を各座標軸とする色空間を所定間隔で画定して複数の空間領域を構成し、任意のコマ画像の構成画素の属する空間領域に当該コマ画像に対するカラー補正量を割り付ける処理を全コマ画像に対して実行する補正量割付ステップと、空間領域毎に割り付けられたカラー補正量に基づいて領域に対する補正量を求める領域補正量演算ステップと、特定のコマ画像の構成画素の属する空間領域から抽出した前記領域補正量を当該特定のコマ画像の構成要素に対する補正量、或いは補正量算出の要因として求める揃い補正量演算ステップとを備え、前記揃い補正量演算ステップは、各空間領域に属する画素数が当該特定のコマ画像の所定比率以上となる空間領域から前記領域補正量を抽出するように構成される写真画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ネガフィルム等の写真フィルムを読み取り得られたカラー画像データやデジタルスチルカメラ等で撮影されたカラー画像データに対して、モニタや写真プリントを行なう際に、自然なカラーを再現できるようにＲ（赤）、Ｇ（緑）、Ｂ（青）（以下、「ＲＧＢ」と記す。）のカラーバランスを調整する写真画像処理装置及び写真画像処理方法に関する。
【背景技術】
【０００２】
従来の写真プリンタでは、ネガフィルムに記録された画像を色合いよく感光材料である印画紙に焼き付けるための写真画像処理方法として、エバンスの定理に基づくＬＡＴＤ（Ｌａｒｇｅ  Ａｒｅａ  Ｔｒａｎｓｍｉｔｔａｎｃｅ  Ｄｅｎｓｉｔｙ）露光方式が知られている。この露光方式は、平均的な戸外の被写体は、ネガ全体の色を混ぜ合わせると灰色に近くなるというエバンスの説に基づいて、色に偏りが見られる場合は、ネガフィルムを透過したＲＧＢの積算光が印画紙上でグレーに再現されるようにＲＧＢの各露光量を調節して露光を行う方式である。
【０００３】
具体的には、ネガフィルムに光を照射して透過光を撮像素子で読取ってＲＧＢのカラー画像データを生成し、全構成画素に対するＲＧＢ成分毎の平均画素値を演算導出し、ＲＧＢ各平均画素値がそれぞれグレーに対応する所定の値となるように、アナログ方式の写真プリンタでは調光フィルタを調節して印画紙を露光し、デジタル方式の写真プリンタではＲＧＢ夫々の光源からの露光量を調節していた。
【０００４】
また、従来の写真プリンタでは、前記写真プリンタに搭載されたＹＭＣ等のカラー補正キーを用いたオペレータによるマニュアル補正や、予め記憶されている関数に従って入力画素値を変換して出力するルックアップテーブルに基づく補正を行うことで色の偏りを是正していた。
【０００５】
上述した従来の写真画像処理方法によれば、被写体（人物、背景）の色の偏りにより過補正され、却って見辛い写真プリントが出力されてしまうという問題があった。例えば、芝生を背景に人物を撮影したシーンの場合には、芝生の領域がグレーに仕上がる一方、人物の領域に芝生の補色であるマゼンダが強く現われる。このような状況をカラーフェリアといい、その対策としてＬＡＴＤ露光方式において高彩度画素を除去したり、彩度によって重み付けをした条件付平均値を求める方法などが提案されている。
【特許文献１】特開２０００−３３０２２１号公報
【発明の開示】
【発明が解決しようとする課題】
【０００６】
上述の方法によりカラーバランスが調整された画像データに基づいて得られる写真プリントの品質として重要視される項目の一つに複数コマ間の揃い、つまり、仕上がりの一連性がある。同一シーン間における同一の被写体の色が同一若しくは略同一または同一シーン間における背景の色が同一若しくは略同一の仕上がりにした方が鑑賞し易く、色濃度が異なると違和感が生じるからである。
【０００７】
このため、従来の写真フィルムでは、撮影シーンの類似性の高いフィルム一本分の全コマ画像に対してＬＡＴＤ処理がなされていたが、短い撮影時間帯で多数枚の画像が撮影される頻度の高いデジタルスチルカメラによる写真画像も同様に複数のコマ画像を加味したカラー補正処理が望まれる。
【０００８】
しかし、そのようなフィルム一本分に撮影された多数のコマ画像に基づきカラーバランスが調整された写真プリントと、そのフィルムが数コマ単位でカットされたピースフィルムに基づきリプリントされた写真プリントでは色の差が生じるなどの問題があり、また、デジタルスチルカメラによる写真画像であっても、カメラの時計が停止しているために撮影時刻が不明な画像や、ファイル名やＥｘｉｆ情報が書き換えられたような画像ファイルでは、複数の画像が混在し、撮影機種や撮影日時などの判別ができないため、コマ画像毎にカラー補正処理をせざるを得ず、仕上がりの一連性を確保できないという問題があった。
【０００９】
本発明の目的は、上述の従来欠点に鑑み、コマ画像毎にカラー補正処理して得られた補正量と補正前のコマ画像情報とに基づいて、複数のコマ画像間の仕上がりの一連性を確保するための適正な補正量を獲得する揃い処理が行なえる写真画像処理装置及び写真画像処理方法を提供する点にある。
【課題を解決するための手段】
【００１０】
上述の目的を達成するため、本発明による写真画像処理装置の第一の特徴構成は、特許請求の範囲の書類の請求項１に記載した通り、一連のコマ画像に対してコマ画像単位にその構成画素に対するカラー補正量を求める補正量演算手段と、三原色のカラー成分を各座標軸とする色空間を所定間隔で画定して複数の空間領域を構成し、任意のコマ画像の構成画素の属する空間領域に当該コマ画像に対する前記カラー補正量を割り付ける処理を全コマ画像に対して実行する補正量割付手段と、空間領域毎に割り付けられた前記カラー補正量に基づいて領域に対する補正量を求める領域補正量演算手段と、特定のコマ画像の構成画素の属する空間領域から抽出した前記領域補正量を当該特定のコマ画像の構成要素に対する補正量、或いは補正量算出の要因として求める揃い補正量演算手段とを備え、前記揃い補正量演算手段は、各空間領域に属する画素数が当該特定のコマ画像の所定比率以上となる空間領域から前記領域補正量を抽出するように構成される点にある。
【００１１】
同一シーンや類似シーンとは、「コマ内に同じ被写体や背景が含まれたコマ同士である」と定義することができ、コマ単位に求められたカラー補正量は仕上がりの連続性が考慮されずに個々に補正されているので、元々同じ物体であったものが違う色に変化する可能性がある。これが仕上がりの一連性が確保されていないという状態である。逆に原画像は補正していないのであるから同一物体は同じ色である可能性が高い。
【００１２】
よって、上述の構成によれば、構成された複数の空間領域において、同一シーンや類似シーンの各コマ画像における各画素が同一の空間領域を占める場合は、当該各画素は各コマ画像の同一物体を構成する画素であるとみなすことができる。そして、同一シーンや類似シーンにおける一連のコマ画像における同一領域に着目した場合に、これら一連のコマ画像のカラー補正量が同一領域に存在していると、領域毎のカラー補正量を算出することで、一連のコマ画像に最適な一つの補正値を得ることができる。
【００１３】
さらに、特定のコマ画像における特定の空間領域を構成する画素の数が多いことは、前記特定のコマ画像における前記特定の空間領域に対応する物体の占める面積が大きいことを示しており、各空間領域に属する画素数が前記特定のコマ画像の所定比率以上となる空間領域から前記領域補正量を抽出することで、各コマ画像において、占める面積が大きい物体を重視した揃い補正量演算処理を行うことができる。
【００１４】
同第二の特徴構成は、同請求項２に記載した通り、上述の第一特徴構成に加えて、前記揃い補正量演算手段により求められた揃い補正量と、対応するコマ画像のカラー補正量とを所定の重み係数で融合した適正カラー補正量を求める適正補正量演算手段を備えて構成される点にある。
【００１５】
特定のコマ画像について算出したカラー補正量のみを当該コマ画像に適用した場合は、同一シーンや類似シーンにおける複数のコマ画像で夫々異なるカラー補正量を適用することになるので、複数のコマ画像間における同一物体の色が異なる色となり、同一物体の色の揃いに一連性がなくなる。一方、補正を実行する全てのコマ画像に対して揃い補正量のみを適用した場合は、複数のコマ画像間の色の揃いの一連性は維持できるが、個々のコマ画像に対するカラー補正の精度が落ちる。そこで、上述の構成によれば、所定の重み係数で前記カラー補正量と前記揃い補正量を融合して適性カラー補正量を求めることによって、各コマ画像に対するカラー補正の精度とコマ画像間の色の揃い一連性を維持することができる。
【００１６】
同第三の特徴構成は、同請求項３に記載した通り、上述の第一または第二の特徴構成に加えて、前記補正量割付手段は、一連のコマ画像の空間領域への分布状態に基づいて領域サイズを調整する領域サイズ調整手段を備えている点にある。
【００１７】
上述の構成によれば、領域サイズを調整することによって、コマ画像間の色の揃いの一連性の効果を制御することができる。例えば、前記領域サイズを大きく設定すると、各空間領域に複数のコマ画像の画素が属する傾向が強くなるので、領域補正量は複数のコマ画像を反映させたものとなり、前記色の揃いの一連性の効果が強くなる。一方、前記領域サイズを小さく設定すると、各空間領域には一つのコマ画像の画素のみが属する傾向が強くなるので、領域補正量は特定のコマ画像を反映させたものとなり、前記色の揃いの一連性の効果が弱くなる。
【００１８】
上述の目的を達成するため、本発明による写真画像処理方法の特徴構成は、特許請求の範囲の書類の請求項４に記載した通り、一連のコマ画像に対してコマ画像単位にカラー補正量を求める補正量演算ステップと、三原色のカラー成分を各座標軸とする色空間を所定間隔で画定して複数の空間領域を構成し、任意のコマ画像の構成画素の属する空間領域に当該コマ画像に対するカラー補正量を割り付ける処理を全コマ画像に対して実行する補正量割付ステップと、空間領域毎に割り付けられたカラー補正量に基づいて領域に対する補正量を求める領域補正量演算ステップと、特定のコマ画像の構成画素の属する空間領域から抽出した前記領域補正量を当該特定のコマ画像の構成要素に対する補正量、或いは補正量算出の要因として求める揃い補正量演算ステップとを備え、前記揃い補正量演算ステップは、各空間領域に属する画素数が当該特定のコマ画像の所定比率以上となる空間領域から前記領域補正量を抽出するように構成される点にある。
【発明の効果】
【００１９】
以上説明した通り、本発明によればコマ画像毎にカラー補正処理して得られた補正量と補正前のコマ画像情報とに基づいて、複数のコマ画像間の仕上がりの一連性を確保するための適正な補正量を獲得する揃い処理が行なえる写真画像処理装置及び写真画像処理方法を提供することができるようになった。
【発明を実施するための最良の形態】
【００２０】
以下に本発明による画像処理装置が組み込まれた写真画像処理装置の実施の形態を説明する。図２に示すように、写真画像処理装置１は、印画紙Ｐに対して出力画像データに基づいた露光処理を行ない露光された印画紙を現像処理する写真プリンタ２と、現像済みの写真フィルムＦから画像を読み込むフィルムスキャナ３１やデジタルスチルカメラ等で撮影された画像データが格納されたメモリカード等の画像データ記憶メディアＭから画像データを読み取るメディアドライバ３２や、コントローラ３３としての汎用コンピュータ等を備え、入力された原画像としての写真画像に対するプリントオーダ情報を設定入力するとともに、各種の画像補正処理を行なう操作ステーション３を備えて構成され、前記操作ステーション３で原画像から編集処理されたプリントデータが前記写真プリンタ２に出力されて所望の写真プリントが生成される。
【００２１】
前記写真プリンタ２は、図２及び図３に示すように、ロール状の印画紙Ｐを収容した二系統の印画紙マガジン２１と、前記印画紙マガジン２１から引き出された印画紙Ｐを所定のプリントサイズに切断するシートカッター２２と、切断後の印画紙Ｐの背面にコマ番号等のプリント情報を印字するバックプリント部２３と、前記プリントデータに基づいて印画紙Ｐを露光する露光部２４と、露光後の印画紙Ｐを現像、漂白、定着するための各処理液が充填された複数の処理槽２５ａ、２５ｂ、２５ｃを備えた現像処理部２５が印画紙Ｐの搬送経路に沿って配置され、現像処理後に乾燥処理された印画紙Ｐが排出される横送りコンベア２６と、横送りコンベア２６に集積された複数枚の印画紙（写真プリント）Ｐがオーダー単位で仕分けられるソータ２７を備えて構成される。
【００２２】
前記露光部２４には、搬送機構２８によって副走査方向に搬送される印画紙Ｐに対して、搬送方向に直交する主走査方向に前記プリントデータに基づき変調されたＲＧＢ三色のレーザ光線束を出力して露光する露光ヘッド２４ａが収容されている。
【００２３】
前記搬送経路に沿って配置された前記露光部２４や現像処理部２５に応じたプロセス速度で印画紙Ｐを搬送する複数のローラ対でなる搬送機構２８が配置され、前記露光部２４の前後には印画紙Ｐを複列に搬送可能なチャッカー式搬送機構２８ａが設けられている。
【００２４】
前記操作ステーション３に設けられたコントローラ３３には、汎用のオペレーティングシステムの管理下で動作し、前記写真処理装置１の各種制御が実行されるアプリケーションプログラムがインストールされ、オペレータとの操作インターフェースとしてモニタ３４、キーボード３５、マウス３６等が接続されている。
【００２５】
前記コントローラ３３のハードウェア及びソフトウェアの協働により実行される写真処理プロセスを機能ブロックで説明すると、図４に示すように、前記フィルムスキャナ３１やメディアドライバ３２によって読み取られた写真画像データを受け取り、所定の前処理を行なって後述のメモリ４１に転送する画像入力部４０と、前記モニタ３４の画面にプリントオーダ情報や画像編集情報を表示するとともに、それらに対して必要なデータ入力のための操作用アイコンを表示するグラフィック操作画面を生成し、或いは表示されたグラフィック操作画面に対する前記キーボード３５やマウス３６からの入力操作に基づいて各種の制御コマンドを生成するグラフィックユーザーインターフェース部４２と、前記画像入力部４０から転送される写真画像データ及び後述の画像処理部４７による補正処理後の写真画像データやそのときの補正パラメータ、更には設定されたプリントオーダ情報等が所定領域に区画されて格納されるメモリ４１と、プリントオーダ情報を生成するオーダー処理部４３と、前記メモリ４１に格納された各写真画像データに対してコマ画像毎または所定枚数のコマ画像に濃度補正処理やコントラスト補正処理等を行なう画像処理部４７と、前記グラフィックユーザーインターフェース部４２からの表示コマンドに基づいて前記メモリ４１に展開された画像データや各種の入出力用グラフィックデータ等を前記モニタ３４に表示処理するビデオＲＡＭ等を備えた表示制御部４６と、各種の補正処理が終了した最終の補正画像を前記写真プリンタ２に出力するためのプリントデータを生成するプリントデータ生成部４４と、顧客のオーダーに応じて最終の補正画像をＣＤ−Ｒ等の記憶媒体に書き込むためのファイル形式に変換するフォーマッタ部４５等で構成される。
【００２６】
前記フィルムスキャナ３１は、フィルムＦに記録された画像を低解像度ではあるものの高速で読み取るプレスキャンモードと、低速ではあるものの高解像度で読み取る本スキャンモードの二モードで作動するように構成され、プレスキャンモードで読み込まれた低解像度の画像に対して各種の補正処理が行なわれ、その際に前記メモリ４１に記憶された補正パラメータに基づいて本スキャンモードで読み込まれた高解像度の画像に対する最終の補正処理が実行されて前記プリンタ２に出力される。
【００２７】
同様に、前記メディアドライバ３２から読み込まれた画像ファイルには高解像度の撮影画像とそのサムネイル画像が含まれ、サムネイル画像に対して後述の各種の補正処理が行なわれ、その際に前記メモリ４１に記憶された補正パラメータに基づいて高解像度の撮影画像に対する最終の補正処理が実行される。尚、画像ファイルにサムネイル画像が含まれないときには、前記画像入力部４０で高解像度の撮影画像からサムネイル画像が生成されて前記メモリ４１に転送される。このように、低解像度の画像に対して頻繁に試行錯誤される各種の編集処理が実行されることによりコントローラ３３の演算負荷が低減されるように構成されている。
【００２８】
前記画像処理部４７は、前記メモリ４１に格納された各写真画像データに対してコマ画像毎または所定枚数のコマ画像に濃度補正処理やコントラスト補正処理等を行なうように構成されており、図１に示すように、前記原画像データに対してカラー補正処理を実施するカラー補正処理部４７１と、前記原画像データに対して濃度補正処理を実施する濃度補正処理部４７２と、前記濃度補正処理部４７２で濃度補正のために生成された濃度補正曲線に対してコントラスト補正を施すコントラスト補正部４７３と、撮影レンズに起因する歪を補正する歪補正部４７４と、画像のエッジを強調しノイズを低減する鮮鋭化処理部４７５と、写真プリントのサイズに適した画像サイズに変換する拡縮処理部４７６等を備えて構成されている。
【００２９】
前記カラー補正処理部４７１は、一連のコマ画像に対してコマ画像単位にカラー補正量を求める補正量演算手段１０と、三原色のカラー成分を各座標軸とする色空間を所定間隔で画定して複数の空間領域を構成し、任意のコマ画像の構成画素の属する空間領域に当該コマ画像に対するカラー補正量を割り付ける処理を全コマ画像に対して実行する補正量割付手段１１と、空間領域毎に割り付けられたカラー補正量に基づいて領域に対する補正量を求める領域補正量演算手段１２と、特定のコマ画像の構成画素の属する空間領域から抽出した前記領域補正量を当該特定のコマ画像の構成要素に対する補正量、或いは補正量算出の要因として求める揃い補正量演算手段１３と、前記揃い補正量演算手段１３により求められた揃い補正量と、対応するコマ画像のカラー補正量とを所定の重み係数で融合した適正カラー補正量を求める適正補正量演算手段１４と、前記揃い補正量演算手段１３により求められた適正カラー補正量に基づいて対応するコマ画像をカラー補正するカラー補正手段１５とを備えて構成されている。
【００３０】
前記補正量演算手段１０は、一連のコマ画像に対してコマ画像単位にその構成画素に対するカラー補正量を求めるように構成されている。
【００３１】
以下に挙げる構成は、前記補正量演算手段１０の一例である。つまり、コマ画像単位に全構成画素に対するカラー成分毎の平均画素値を求める平均画素値演算手段１００と、前記コマ画像の各画素を、中心からの距離で彩度が規定され中心周りの角度で色相が規定される色相彩度座標系の画素に変換する座標変換手段１０１と、前記色相彩度座標系に前記平均画素値を中心とする同心円状の複数の円領域を画定し、各円領域に含まれる画素の分布特性を示す円領域特性値を領域毎に求める円領域特性値演算手段１０２と、前記円領域特性値に基づいてカラー補正量を求めるカラー補正量演算手段１０３とを備えた構成である。
【００３２】
前記平均画素値演算手段１００は、入力されたコマ画像単位に全構成画素に対するカラー成分毎の平均画素値を求めるように構成されている。具体的には、前記コマ画像の画素数がｉ個であり、前記コマ画像のｎ番目の画素のＲＧＢ成分を夫々ｒｎ、ｇｎ、ｂｎとした場合に、前記コマ画像の各画素のＲＧＢ成分の平均値ａｖｇ（ｒ）、ａｖｇ（ｇ）、ａｖｇ（ｂ）は、〔数１〕のように表される。
【００３３】
【数１】
0005168809.tif 000002
0005168809.tif 000003
0005168809.tif 000004
【００３４】
前記座標変換手段１０１は、前記コマ画像の各画素を、中心からの距離で彩度が規定され中心周りの角度で色相が規定される色相彩度座標系の画素に変換するように構成されており、詳述すると、前記コマ画像の各画素ｒｎ、ｇｎ、ｂｎを〔数２〕に示す演算式に基づいてｘｙ座標系に変換する。前記ｘｙ座標系は、中心からの距離で彩度が規定され中心周りの角度で色相が規定される色相彩度座標系として示されている。
【００３５】
【数２】
0005168809.tif 000005
【００３６】
尚、〔数２〕で示す演算式は、図５に示すように、ｘｙ座標系におけるｘの正方向を０度の方向とした場合に、Ｒ成分を０度の方向、Ｇ成分を１２０度の方向、及びＢ成分を２４０度の方向として、並びに各成分の大きさを各成分方向への原点からの距離として、各成分の合成ベクトルＣＲＧＢのｘ成分とｙ成分の大きさへの変換を表している。
【００３７】
前記円領域特性値演算手段１０２は、前記色相彩度座標系に前記平均画素値を中心とする同心円状の複数の円領域を画定し、各円領域に含まれる画素の分布特性を示す円領域特性値を領域毎に求めるように構成されている。
【００３８】
詳述すると、図６に示すように、前記平均画素値をｘｙ座標系に変換した値、つまり〔数１〕で導出したａｖｇ（ｒ）、ａｖｇ（ｇ）、ａｖｇ（ｂ）を夫々〔数２〕における画素ｒｎ、ｇｎ、ｂｎに代入して算出したｘｎとｙｎを夫々中心座標ａｖｇ（ｘ）とａｖｇ（ｙ）として前記色相彩度座標系にプロットし、前記中心座標を中心とした所定の間隔の同心円状の１２個の円領域Ｓ１からＳ１２を画定する。
【００３９】
尚、前記円領域において、外側に位置する円領域は、当該円領域よりも内側に位置する円領域を包含している。例えば、円領域Ｓ２は円領域Ｓ１を包含しており、円領域Ｓ１２は円領域Ｓ１１を包含（つまり、円領域Ｓ１からＳ１１の全てを包含）している。また、本例では、円領域の数は１２個であるが、円領域の数は適宜設定されてもよい。また、本例では、円領域間の所定の間隔は等間隔であるが、前記所定の間隔も適宜設定されてもよく、等間隔で異なる大きさであっても各領域の間隔の大きさが互いに異なっていてもよい。
【００４０】
前記円領域特性値は、各円領域に含まれる画素の画素数、及び各円領域に含まれる画素のカラー成分毎の平均値である領域平均画素値と前記平均画素値との差分値で構成される。前記円領域特性値の算出は４ステップで構成されており、以下に詳述する。
【００４１】
第一ステップでは、コマ画像の全構成画素を前記色相彩度座標系にプロットして、同心円上の複数の領域Ｓｊ（ｊは１から１２の整数）毎にプロットされた画素数ｍＳｊをカウントする。
【００４２】
第二ステップでは、前記領域Ｓｊにおいてプロットされた画素のうちｎ番目の画素のＲＧＢ成分を夫々ｒＳｊｎ、ｇＳｊｎ、ｂＳｊｎとした場合に、前記領域Ｓｊを構成しているＲＧＢ各成分の平均画素値ａｖｇ（ｒＳｊｎ）、ａｖｇ（ｇＳｊｎ）、ａｖｇ（ｂＳｊｎ）を〔数３〕に示す演算式に基づいて算出する。
【００４３】
【数３】
0005168809.tif 000006
0005168809.tif 000007
0005168809.tif 000008
【００４４】
第三ステップでは、前記平均画素値演算手段１００で算出した前記コマ画像の各画素のＲＧＢ成分の平均値である前記平均画素値と前記領域Ｓｊを構成しているＲＧＢ各成分の平均画素値との差分値ΔｒＳｊ、ΔｇＳｊ、ΔｂＳｊを〔数４〕に示す演算式に基づいて算出する。
【００４５】
【数４】
0005168809.tif 000009
【００４６】
以上の３ステップより算出された特性値ｍＳｊ、ΔｒＳｊ、ΔｇＳｊ、ΔｂＳｊを０から１の範囲の値へと正規化することによって、円領域特性値ｍｊ、Δｒｊ、Δｇｊ、Δｂｊを算出する。具体的には、前記円領域特性値ｍｊは、前記特性値ｍＳｊを前記コマ画像に含まれる画素の総数で割ることによって算出され、前記円領域特性値Δｒｊ、Δｇｊ、Δｂｊは、前記特性値ΔｒＳｊ、ΔｇＳｊ、ΔｂＳｊを前記コマ画像における濃度設定値の最大値で夫々割ることによって算出される。
【００４７】
尚、前記特性値ｍＳｊ、ΔｒＳｊ、ΔｇＳｊ、ΔｂＳｊをそのまま前記円領域特性値とする構成であってもよい。
【００４８】
尚、図７に示すように、前記領域Ｓｊを横軸、各円領域に含まれる画素の画素数ｍＳｊを縦軸として画素数の累積ヒストグラムを示すと、カラーフェリアが発生しないコマ画像は、前記色相彩度座標系におけるデータの分布が小さくなるために、前記累積ヒストグラムは、内側の領域（図７ではＳ５）で飽和する傾向にあるが、カラーフェリアが発生するコマ画像や異種光源のコマ画像は、前記色相彩度座標系におけるデータの分布が大きくなるために、前記累積ヒストグラムは、より外側の領域（図７では、異種光源はＳ７、フェリアはＳ１０）で飽和する、つまりなかなか飽和しない傾向にある。
【００４９】
前記カラー補正量演算手段１０３は、前記円領域特性値ｍｊ、Δｒｊ、Δｇｊ、Δｂｊ（または、ｍＳｊ、ΔｒＳｊ、ΔｇＳｊ、ΔｂＳｊ）に基づいてカラー補正量Ｒ１、Ｇ１、Ｂ１を求める。前記カラー補正量Ｒ１、Ｇ１、Ｂ１は、元の画素値と新たな画素値との差分であり、前記カラー補正量演算手段１０３は、例えば、誤差逆伝搬学習法により学習されるニューラルネットワーク、または、重回帰分析手段で構成される。
【００５０】
誤差逆伝搬学習法により学習されるニューラルネットワークで構成される場合は、例えば、図８に示すように、入力層、中間層、出力層の三層で構成されたパーセプトロン形で構成されており、前記コマ画像の各領域Ｓｊの円領域特性値ｍｊ、Δｒｊ、Δｇｊ、Δｂｊが入力される入力層と、入力データが所定の結合加重で重み演算された結果が入力される中間層と、中間層の出力が所定の結合加重で重み演算されてカラー補正量Ｒ１、Ｇ１、Ｂ１として出力される出力層の三層から構成され、各結合加重が誤差逆伝搬学習法により予め設定される。
【００５１】
ニューラルネットワークは、具体的には撮影条件の異なる９０００枚のサンプル画像に対してなされたマニュアル補正によるカラー補正量及びサンプル画像に対する円領域特性値を教師データとして、サンプル画像毎に各々一億回の繰り返し学習を行なって最もばらつきが少なかった結果を得たものを設定している。
【００５２】
重回帰分析手段で構成される場合は、例えば、〔数５〕に示すような重回帰式に、従属変数として前記円領域特性値ｍｊ、Δｒｊ、Δｇｊ、Δｂｊを夫々適用することによって、独立変数としてのＲＧＢ成分毎のカラー補正量Ｒ１、Ｇ１、Ｂ１を導出する構成が挙げられる。
【００５３】
【数５】
0005168809.tif 000010
【００５４】
上述の〔数５〕において、Ａ１１からＡ３４は係数、Ｃ１からＣ３は定数項で、多数のサンプル画像に対するマニュアル補正で適正と評価されたときの前記サンプル画像と補正量との関係から算出される値である。
【００５５】
以上説明した前記補正量演算手段１０の一例を採用する効果としては以下が挙げられる。コマ画像単位に行なわれる補正量演算に上述の構成を採用することで、太陽光やタングステン光等の撮影光源を含む撮影条件の差異や、フィルム種の差異等があると、基本的にコマ画像全体にそれらの影響が及び、同一シーンを撮影しても、フィルム上に形成された画像は見た目にも異なる色合いを呈している。しかし、コマ画像に写っている物体間の相対的な色に注目すれば著しい違いがあるわけでない。そこで、ＬＡＴＤ補正を行なうと撮影条件やフィルム種等に基づくコマ画像全体の特徴が無彩色として全体として補正され、平均化が過剰に作用したカラーフェリア現象だけが不都合として残ることになる。
【００５６】
色相彩度座標系に前記平均画素値演算手段１００により求められた平均画素値を中心とする同心円上の複数の円領域を画定し、各領域におけるコマ画像の画素の分布特性を評価すると、カラーフェリアの起こらないコマ画像においては各円領域の分布特性がほぼ等しくなるが、カラーフェリアの傾向が強くなるコマ画像においては各円領域間の分布特性の差が大きくなり、カラーフェリアの状況が特徴的に反映される。つまり、中心側の円領域から外側の円領域に沿って構成画素数の累積ヒストグラムをとると、カラーフェリアの傾向が弱いコマ画像は比較的低彩度領域、つまり中心側の円領域で飽和する傾向があるのに対して、カラーフェリアの傾向が強いコマ画像は比較的高彩度領域、つまり外側の円領域まで飽和することが無いという顕著な相違が見られるのである。この傾向は、画素数の分布特性のみならず画素値も同様である。
【００５７】
そこで、各円領域に含まれる画素の分布特性を示す円領域特性値に基づいて所定の分布特性を示すようにカラー補正量を演算するカラー補正量演算手段１０３を設けて、平均画素値演算手段１００により求められた平均画素値を補正することにより、カラーフェリアの発生を回避した適切なカラー補正量を求めることができるようになるのである。
【００５８】
前記補正量割付手段１１は、三原色のカラー成分を各座標軸とする色空間を所定間隔で画定して複数の空間領域を構成し、任意のコマ画像の構成画素の属する空間領域に当該コマ画像に対するカラー補正量を割り付ける処理を全コマ画像に対して実行するように構成されている。
【００５９】
詳述すると、前記コマ画像が、前記三原色のカラー成分、つまり、赤色成分（Ｒ成分）と緑色成分（Ｇ成分）と青色成分（Ｂ成分）が夫々１２ビットで表現されている場合は、例えば、図９に示すように、前記三原色の１２ビットを夫々８４個の領域、つまり、赤色成分をＲ１からＲ８４、緑色成分をＧ１からＧ８４、青色成分をＢ１からＢ８４に分割することで、領域（１，１，１）から領域（８４，８４，８４）の合計８４×８４×８４個の空間領域が構成される。
【００６０】
構成された前記空間領域に対する前記カラー補正量の割り付け処理は、任意のコマ画像の全ての画素に対して、前記構成画素のうち少なくとも一個の画素が存在する前記空間領域に、前記補正量演算手段１０で算出された当該コマ画像のカラー補正量を割り当てる。同様の割り付け処理を、他の全てのコマ画像に対しても行う。
【００６１】
例えば、ＲＧＢ成分のカラー補正量が夫々ΔＲＡ、ΔＧＡ、ΔＢＡのコマ画像Ａとカラー補正量が夫々ΔＲＢ、ΔＧＢ、ΔＢＢのコマ画像Ｂに対して上述の割り付け処理を実行する場合、前記コマ画像Ａの画素が少なくとも１個存在しており前記コマ画像Ｂの画素が存在していない空間領域には、前記カラー補正量ΔＲＡ、ΔＧＡ、ΔＢＡが割り付けられ、前記コマ画像Ｂの画素が少なくとも１個存在しており前記コマ画像Ａの画素が存在していない空間領域には、前記カラー補正量ΔＲＢ、ΔＧＢ、ΔＢＢが割り付けられ、前記コマ画像Ａの画素が少なくとも１個存在しており且つ前記コマ画像Ｂの画素が少なくとも１個存在している空間領域には、前記カラー補正量ΔＲＡ、ΔＧＡ、ΔＢＡとΔＲＢ、ΔＧＢ、ΔＢＢが割り付けられ、前記コマ画像Ａと前記コマ画像Ｂの何れにも画素が存在していない空間領域には、前記カラー補正量が割り付けられない。
【００６２】
また、前記補正量割付手段１１は、一連のコマ画像の空間領域への分布状態に基づいて領域サイズを調整する領域サイズ調整手段１１０を備えて構成されており、例えば、前記領域サイズ調整手段１１０は、図９では前記三原色の１２ビットを夫々８４個の領域に分割していたが、該領域の分割を夫々４２個とすることで合計４２×４２×４２個の空間領域に領域サイズを調整する構成や、図１０（ａ）に示すように、前記三原色の各色成分で領域サイズが異なるように調整する構成や、図１０（ｂ）に示すように、前記三原色の１２ビットの低彩度領域と高彩度領域で領域サイズが異なるように調整する構成や、上述のサイズ調整を適宜組み合わせた構成が挙げられる。
【００６３】
尚、領域サイズの調整は、入力された各コマ画像の色相、彩度、カラー補正量等により自動で調整する方法であってもよく、オペレータが手動で調整する方法であってもよい。
【００６４】
前記領域補正量演算手段１２は、空間領域毎に割り付けられた前記カラー補正量に基づいて領域に対する補正量を求めるように構成されており、例えば、前記領域に対する補正量として、補正量を平均した値である領域平均補正量ａｖｇ（ΔＲ）、ａｖｇ（ΔＧ）、ａｖｇ（ΔＢ）を求めるように構成されている。詳述すると、任意の空間領域における前記領域補正量を算出するためには、当該空間領域に割り付けられた全てのカラー補正量を合算してから、当該空間領域にカラー補正量を割り付けたコマ画素の数で除する。
【００６５】
例えば、コマ画像ＡのＲＧＢ成分のカラー補正量ΔＲＡ、ΔＧＡ、ΔＢＡとコマ画像Ｂのカラー補正量ΔＲＢ、ΔＧＢ、ΔＢＢが割り付けられている空間領域に対して前記領域平均補正量ａｖｇ（ΔＲ）、ａｖｇ（ΔＧ）、ａｖｇ（ΔＢ）を求める場合に、前記カラー補正量ΔＲＡ、ΔＧＡ、ΔＢＡのみが割り付けられている空間領域の前記領域平均補正量ａｖｇ（ΔＲ）、ａｖｇ（ΔＧ）、ａｖｇ（ΔＢ）は夫々ΔＲＡ、ΔＧＡ、ΔＢＡであり、前記カラー補正量ΔＲＢ、ΔＧＢ、ΔＢＢのみが割り付けられている空間領域の前記領域平均補正量ａｖｇ（ΔＲ）、ａｖｇ（ΔＧ）、ａｖｇ（ΔＢ）は夫々ΔＲＢ、ΔＧＢ、ΔＢＢであり、前記カラー補正量ΔＲＡ、ΔＧＡ、ΔＢＡとΔＲＢ、ΔＧＢ、ΔＢＢのみが割り付けられている空間領域の前記領域平均補正量ａｖｇ（ΔＲ）、ａｖｇ（ΔＧ）、ａｖｇ（ΔＢ）は夫々（ΔＲＡ＋ΔＲＢ）／２、（ΔＧＡ＋ΔＧＢ）／２、（ΔＢＡ＋ΔＢＢ）／２である。
【００６６】
前記揃い補正量演算手段１３は、特定のコマ画像の構成画素の属する空間領域から抽出した前記領域補正量を当該特定のコマ画像の構成要素に対する補正量、或いは補正量算出の要因として求めるように構成されており、例えば、特定のコマ画像の構成画素の属する空間領域から抽出した前記領域補正量の平均を当該特定のコマ画像の揃い補正量として求めるように構成されている。具体的には、前記空間領域から抽出した前記領域補正量を合算してから、前記領域補正量が抽出された空間領域の数で除する。
【００６７】
例えば、３画素より構成されるコマ画像Ａ（カラー補正量はΔＲＡ、ΔＧＡ、ΔＢＡ）の構成画素が、夫々領域Ｐ（１，１，１）、領域Ｑ（２，１０，１５）、領域Ｒ（ｌ，ｍ，ｎ）に属しており、夫々３画素より構成されるコマ画像Ｂ（カラー補正量はΔＲＢ、ΔＧＢ、ΔＢＢ）、コマ画像Ｃ（カラー補正量はΔＲＣ、ΔＧＣ、ΔＢＣ）と共に各コマ画像の全構成画素が空間領域に割り付けられて領域補正量、例えば、領域平均補正量が求められた結果、前記領域Ｐ（１，１，１）の前記領域平均補正量ａｖｇ（ΔＲＰ）、ａｖｇ（ΔＧＰ）、ａｖｇ（ΔＢＰ）と、前記領域Ｑ（２，１０，１５）の前記領域平均補正量ａｖｇ（ΔＲＱ）、ａｖｇ（ΔＧＱ）、ａｖｇ（ΔＢＱ）と、前記領域Ｒ（ｌ，ｍ，ｎ）の前記領域平均補正量ａｖｇ（ΔＲＲ）、ａｖｇ（ΔＧＲ）、ａｖｇ（ΔＢＲ）が〔数６〕で示されたとすると、前記コマ画像Ａの揃い補正量ΔＳＲ、ΔＳＧ、ΔＳＢは〔数７〕のように示される。
【００６８】
【数６】
0005168809.tif 000011
【００６９】
【数７】
0005168809.tif 000012
【００７０】
また、前記揃い補正量演算手段１３は、各空間領域に属する画素数が当該特定のコマ画像の所定比率以上となる空間領域から前記領域補正量を抽出するように構成されている。
【００７１】
例えば、画素数が１０００個より構成されるコマ画像があり、前記コマ画像の前記所定比率が５％である場合は、属する画素数が５０個以上の空間領域からは前記揃い補正量を抽出するが、属する画素数が４９個以下の空間領域からは前記領域揃い補正量を抽出しない。つまり、コマ画像全体における占有率が低い空間領域は、揃い補正量の抽出においては当該空間領域に当該コマ画像の画素が存在しないものとして扱われるのである。
【００７２】
前記適正補正量演算手段１４は、前記揃い補正量演算手段１３により求められた揃い補正量ΔＳＲ、ΔＳＧ、ΔＳＢと、対応するコマ画像のカラー補正量ΔＲ、ΔＧ、ΔＢとを所定の重み係数で融合した適正カラー補正量ΔＳＣＲ、ΔＳＣＧ、ΔＳＣＢを求めるように構成されており、例えば、所定の重み係数をｒとした場合、前記適正カラー補正量ΔＳＣＲ、ΔＳＣＧ、ΔＳＣＢを〔数８〕に示す演算式に基づいて算出する。〔数８〕において前記重み係数ｒは、前記揃い補正量と前記カラー補正量のうち、前記カラー補正量の使用割合を示しているが、前記重み係数ｒは、前記揃い補正量の使用割合を示していてもよい。
【００７３】
【数８】
0005168809.tif 000013
【００７４】
〔数８〕における前記重み係数ｒの具体的な値としては、前記揃い補正量と前記カラー補正量を６：４の割合とする値、つまり０．４が好ましい値として挙げられる。これは、図８に示すニューラルネットワークを３段構成としたものに多数のサンプル画像を適用することで得られた値である。
【００７５】
前記カラー補正手段１５は、前記適正補正量演算手段１４により求められた適正カラー補正量に基づいて対応するコマ画像をカラー補正するように構成されており、詳述すると、任意のコマ画像の全ての画素について、前記適正補正量演算手段１４で求められた適正カラー補正量ΔＳＣＲ、ΔＳＣＧ、ΔＳＣＢを加算することによって前記コマ画像を補正する。
【００７６】
以上説明した本発明によるコマ画像の補正に際して、マニュアル補正で適正と評価された画素値を正解画素値とした場合の、キー補正の程度による前記正解画素値を得る確率および不正解であった場合の誤差の標準偏差を、図１１に示すように、従来のＬＡＴＤ補正を用いた方法、前記カラー補正量演算手段１０３の出力であるカラー補正量を用いた方法、前記揃い補正量演算手段１３の出力である揃い補正量を用いた方法、そして前記適正補正量演算手段１４の出力である適正カラー補正量を前記重み係数０．５で用いた方法について示す。
【００７７】
図１１において、前記揃い補正量演算手段１３の出力である揃い補正量を用いた方法は、従来のＬＡＴＤ補正を用いた方法や前記カラー補正量演算手段１０３の出力であるカラー補正量を用いた方法よりも向上しており、前記適正補正量演算手段１４を用いた方法は、前記揃い補正量演算手段１３を用いた方法よりも向上している。
【００７８】
以上説明した本発明による写真画像処理装置は、一連のコマ画像に対してコマ画像単位にカラー補正量を求める補正量演算手段１０と、三原色のカラー成分を各座標軸とする色空間を所定間隔で画定して複数の空間領域を構成し、任意のコマ画像の構成画素の属する空間領域に当該コマ画像に対するカラー補正量を割り付ける処理を全コマ画像に対して実行する補正量割付手段１１と、空間領域毎に割り付けられたカラー補正量に基づいて領域に対する補正量を求める領域補正量演算手段１２と、特定のコマ画像の構成画素の属する空間領域から抽出した前記領域補正量を当該特定のコマ画像の構成要素に対する補正量、或いは補正量算出の要因として求める揃い補正量演算手段１３としてコンピュータを機能させるための画像処理プログラムがインストールされることにより容易に実現される。
【００７９】
同様に、入力されたコマ画像単位に全構成画素に対するカラー成分毎の平均画素値を求める平均画素値演算手段１００と、前記コマ画像の各画素を、中心からの距離で彩度が規定され中心周りの角度で色相が規定される色相彩度座標系の画素に変換する座標変換手段１０１と、前記色相彩度座標系に前記平均画素値を中心とする同心円状の複数の円領域を画定し、各円領域に含まれる画素の分布特性を示す円領域特性値を領域毎に求める円領域特性値演算手段１０２と、前記円領域特性値に基づいてカラー補正量を求めるカラー補正量演算手段１０３としてコンピュータを機能させるための画像処理プログラムがインストールされることにより容易に実現される。
【００８０】
以下、別実施形態について説明する。上述の実施形態では、補正量演算手段１０の一例として、コマ画像単位に全構成画素に対するカラー成分毎の平均画素値を求める平均画素値演算手段１００と、前記コマ画像の各画素を、中心からの距離で彩度が規定され中心周りの角度で色相が規定される色相彩度座標系の画素に変換する座標変換手段１０１と、前記色相彩度座標系に前記平均画素値を中心とする同心円状の複数の円領域を画定し、各円領域に含まれる画素の分布特性を示す円領域特性値を領域毎に求める円領域特性値演算手段１０２と、前記円領域特性値に基づいてカラー補正量を求めるカラー補正量演算手段１０３とを備えた構成について説明したが、下記に示す構成であってもよい。
【００８１】
つまり、写真プリンタに搭載されたＹＭＣ等のカラー補正キーを用いたオペレータによるマニュアル補正や、予め記憶されている関数に従って入力画素値を変換して出力するルックアップテーブルに基づく補正を行う構成であってもよい。これらの構成によれば、コマ画像単位にカラー補正量を算出する際に、前記カラー補正量は、前記コマ画像を構成する各画素について個々に算出されることになる。
【００８２】
上述の実施形態では、揃い補正量演算手段１３は、各空間領域に属する画素数が当該特定のコマ画像の所定比率以上となる空間領域から前記領域補正量を抽出する構成について説明したが、前記領域補正量の抽出基準は比率ではなくて所定画素数であってもよい。
【００８３】
例えば、画素数が１０００個より構成されるコマ画像があり、前記コマ画像の前記所定画素数が１０である場合は、属する画素数が１０個以上の空間領域からは前記揃い補正量を抽出するが、属する画素数が９個以下の空間領域からは前記領域揃い補正量を抽出しない構成が挙げられる。
【００８４】
また、前記揃い補正量演算手段１３において、所定の比率や画素数といった前記領域補正量の抽出基準を設けない構成であってもよい。つまり、少なくとも１個の画素が属する全ての空間領域から揃い補正量を抽出するので、前記揃い補正量は、当該コマ画像を構成する全ての画素を反映したものとなる。
【００８５】
上述の実施形態では、前記補正量演算手段１０は、平均画素値演算手段１００、座標変換手段１０１、円領域特性値演算手段１０２、及びカラー補正量演算手段１０３を備えており、円領域特性値に基づいてカラー補正量を求める構成について説明したが、前記補正量演算手段１０は、これらに追加して、前記カラー補正量で表される画素を中心として、所定角度間隔で周方向に複数の放射状領域を画定し、各放射状領域に含まれる画素の分布特性を示す放射状領域特性値を領域毎に求める放射状領域特性値演算手段１０４と、前記放射状領域特性値に基づいてカラー補正量を求める第二カラー補正量演算手段１０５を備えて構成されていてもよい。
【００８６】
前記放射状領域特性値演算手段１０４は、前記座標変換手段１０１においてコマ画像の各画素が変換させられた色相彩度座標系に、前記カラー補正量演算手段１０３により求められたカラー補正量で表される画素を中心として、所定角度間隔で周方向に複数の放射状領域を画定し、各放射状領域に含まれる画素の分布特性を示す放射状領域特性値を領域毎に求めるように構成されている。
【００８７】
詳述すると、図１２に示すように、前記カラー補正量演算手段１０３で導出したカラー補正量Ｒ１、Ｇ１、Ｂ１を夫々〔数２〕における画素ｒｎ、ｇｎ、ｂｎに代入して算出したｘｎとｙｎを夫々中心座標Ｘ、Ｙとして前記色相彩度座標系にプロットし、前記中心座標を中心として、所定角度間隔で周方向に複数の放射状領域を画定する。尚、放射状領域の数は適宜設定される。例えば、図１２では、３０度毎に１２の放射状領域を設定されており、この設定に基づいて以下の説明を行う。
【００８８】
前記放射状領域特性値は、各放射状領域に含まれる画素の画素数と、各放射状領域に含まれる画素に対して低彩度画素ほど大きな重み係数を掛けて求まるカラー成分毎の重み平均画素値で構成される。前記放射状領域特性値の算出は２ステップで構成されており、以下に詳述する。
【００８９】
第一ステップでは、コマ画像の全構成画素を前記色相彩度座標系にプロットして、放射状の複数の領域Ｈｊ（ｊは１から１２の整数）毎にプロットされた画素数ｌＨｊをカウントする。
【００９０】
第二ステップでは、前記領域Ｈｊにおいてプロットされた画素のうちｎ番目の画素のＲＧＢ成分を夫々ｒＨｊｎ、ｇＨｊｎ、ｂＨｊｎとして、各画素の中心からの距離（彩度）に基づいて算出した重み係数を夫々ｚｒＨｊｎ、ｚｇＨｊｎ、ｚｂＨｊｎとして前記放射状領域特性値ΔｒＨｊ、ΔｇＨｊ、ΔｂＨｊを〔数９〕に示す演算式に基づいて算出する。尚、〔数９〕において、分母に加算されている１と分子に加算されているカラー補正量Ｒ１、Ｇ１、Ｂ１については後述する。
【００９１】
【数９】
0005168809.tif 000014
0005168809.tif 000015
0005168809.tif 000016
【００９２】
つまり、前記領域Ｈｊに存在する各画素に対して、前記各画素に対応する重み係数を乗じた値を算出し、算出した値を合計してから前記領域Ｈｊに存在する画素数ｌＨｊで割ることで、前記放射状領域特性値ΔｒＨｊ、ΔｇＨｊ、ΔｂＨｊは算出される。さらに前記領域Ｈｊに占める画素数ｌＨｊも前記放射状領域特性値とする。つまり、前記放射状領域特性値は、ｌＨｊ、ΔｒＨｊ、ΔｇＨｊ、ΔｂＨｊの４値である。
【００９３】
尚、各画素（例えば画素値をｒＨｊｎ、ｇＨｊｎ、ｂＨｊｎとする）についての前記重み係数ｚｒＨｊｎ、ｚｇＨｊｎ、ｚｂＨｊｎの算出は、各画素の濃度値を以下に説明する方法で正規化した値である。正規化の方法は、前記原画像データにおける濃度設定値の範囲（例えば、原画像データが８ビットで構成される場合には、０から２５５の範囲）を０から４の範囲へと正規化して、正規化した値を１から減じることで正規化する。尚、正規化した値を１から減じた際に負の値となった場合には、重み係数は０とする。このように算出した重み係数は、彩度の低い画素の重みが大きくなり、全濃度範囲の高彩度側７５％に相当する領域では０となる。
【００９４】
また、〔数９〕の演算において、前記領域Ｈｊに存在する画素に追加して、前記カラー補正量Ｒ１、Ｇ１、Ｂ１を、前記重み係数を１として加算している。その結果、画素数が１個増えることから、〔数９〕においては、画素数はｌＨｊではなくｌＨｊ＋１で除している。これは、前記領域Ｈｊに画素が存在しているが、当該画素の全てについて重み係数が０であるために前記放射状領域特性値が０となることと、前記領域Ｈｊに画素が全く存在しないために前記放射状領域特性値が０となることを区別するためである。
【００９５】
前記第二カラー補正量演算手段１０５は、前記放射状領域特性値ｌＨｊ、ΔｒＨｊ、ΔｇＨｊ、ΔｂＨｊに基づいてカラー補正量Ｒ２、Ｇ２、Ｂ２を求める。前記カラー補正量Ｒ２、Ｇ２、Ｂ２は、元の画素値と新たな画素値との差分であり、前記第二補正量演算手段１６は、例えば、前記第一補正量演算手段１０３と同様に、誤差逆伝搬学習法により学習されるニューラルネットワーク、または、重回帰分析手段で構成される。
【００９６】
誤差逆伝搬学習法により学習されるニューラルネットワークで構成される場合は、図８に示す入力層には、前記放射状領域特性値ｌＨｊ、ΔｒＨｊ、ΔｇＨｊ、ΔｂＨｊが入力され、出力層からＲＧＢ各成分に対するカラー補正量Ｒ２、Ｇ２、Ｂ２が出力される。当該ニューラルネットワークも先と同様に、撮影条件の異なる９０００枚のサンプル画像に対してなされたマニュアル補正によるカラー補正量及びサンプル画像の放射状領域特性値を教師データとして、サンプル画像毎に各々一億回の繰り返し学習を行なって最もばらつきが少なかった結果を得たものを設定している。
【００９７】
重回帰分析手段で構成される場合は、例えば、〔数５〕に示すような重回帰式に、従属変数として前記放射状領域特性値ｌＨｊ、ΔｒＨｊ、ΔｇＨｊ、ΔｂＨｊを夫々適用することによって、独立変数としてのＲＧＢ成分毎のカラー補正量Ｒ２、Ｇ２、Ｂ２を導出する構成が挙げられる。
【００９８】
尚、前記放射状領域特性値演算手段１０４において、前記カラー補正量演算手段１０３で導出したカラー補正量Ｒ１、Ｇ１、Ｂ１を、夫々画素ｒｎ、ｇｎ、ｂｎとして〔数２〕に適用して算出したｘｎとｙｎを夫々中心座標Ｘ、Ｙとして前記色相彩度座標系にプロットする構成について説明したが、前記中心座標が図６に示すａｖｇ（ｘ）とａｖｇ（ｙ）であってもよい。
【００９９】
上述した実施形態では、前記円領域特性値が、各円領域に含まれる画素の画素数と、各円領域に含まれる画素のカラー成分毎の平均値である領域平均画素値と前記平均画素値との差分値で構成されるものを説明したが、前記円領域特性値としては、各円領域に含まれる画素の画素数と、各円領域に含まれる画素のカラー成分毎の平均値である領域平均画素値を採用することもできる。
【０１００】
この場合、カラー補正量演算手段１０３において、誤差逆伝搬学習法により学習されるニューラルネットワークで構成される場合は、例えば、図１３に示すように、入力層、中間層、出力層の三層で構成されたパーセプトロン形で構成されており、平均画素値演算手段により求められた各画素のＲＧＢ成分の平均値ａｖｇ（ｒ）、ａｖｇ（ｇ）、ａｖｇ（ｂ）を正規化した値ａｖｇｎ（ｒ）、ａｖｇｎ（ｇ）、ａｖｇｎ（ｂ）と、各領域Ｓｊの円領域特性値ｍｊ、ａｖｇ（ｒｊｎ）、ａｖｇ（ｇｊｎ）、ａｖｇ（ｂｊｎ）（これらは〔数３〕による平均画素値ａｖｇ（ｒＳｊｎ）、ａｖｇ（ｇＳｊｎ）、ａｖｇ（ｂＳｊｎ）を夫々正規化した値である）が入力される入力層と、入力データが所定の結合加重で重み演算された結果が入力される中間層と、中間層の出力が所定の結合加重で重み演算されてカラー補正量Ｒ１、Ｇ１、Ｂ１として出力される出力層の三層から構成され、各結合加重が誤差逆伝搬学習法により予め設定される。
【０１０１】
上述の実施形態では、カラー補正量演算手段１０３や第二カラー補正量演算手段１０５がニューラルネットワークで構成される場合に、前記ニューラルネットワークの中間層が一層で構成されるものについて説明したが、前記中間層が複数層で構成されるものであってもよい。
【０１０２】
上述の実施形態では、補正量演算手段１０として、コマ画像単位に全構成画素に対するカラー成分毎の平均画素値を求める平均画素値演算手段１００と、前記コマ画像の各画素を、中心からの距離で彩度が規定され中心周りの角度で色相が規定される色相彩度座標系の画素に変換する座標変換手段１０１と、前記色相彩度座標系に前記平均画素値を中心とする同心円状の複数の円領域を画定し、各円領域に含まれる画素の分布特性を示す円領域特性値を領域毎に求める円領域特性値演算手段１０２と、前記円領域特性値に基づいてカラー補正量を求めるカラー補正量演算手段１０３とを備えた構成について説明したが、前記補正量演算手段１０は、上述の構成に限定されるものではなく、一連のコマ画像に対してコマ画像単位にカラー補正量を求めるものであれば、上記以外の構成であってもよい。
【０１０３】
例えば、従来のＬＡＴＤ露光方式、つまり、エバンスの説に基づいて、色に偏りが見られる場合は、ネガフィルムを透過したＲＧＢの積算光が印画紙上でグレーに再現されるようにＲＧＢの各露光量を調節して露光を行う方式に基づいてカラー補正量を求める構成や、従来のＬＡＴＤ露光方式において高彩度画素を除去したり、彩度によって重み付けをした条件付平均値を求める方法によってカラー補正量を求める構成等が挙げられる。
【０１０４】
尚、上述した実施形態は、本発明の一例に過ぎず、本発明の作用効果を奏する範囲において各ブロックの具体的構成等を適宜変更設計できることは言うまでもない。
【図面の簡単な説明】
【０１０５】
【図１】画像処理部の機能ブロック構成図
【図２】写真画像処理装置の外観構成の説明図
【図３】写真プリンタの説明図
【図４】写真画像処理装置の機能ブロック構成の説明図
【図５】色相彩度座標系の説明図
【図６】色相彩度座標系における複数の同心円上領域の説明図
【図７】各同心円上領域における画素数の累積ヒストグラムの説明図
【図８】ニューラルネットワークで構成される第二補正量演算手段の説明図
【図９】三原色のカラー成分を各座標軸とする空間領域の説明図
【図１０】（ａ）は、各色成分で領域サイズが異なる空間領域を示し、（ｂ）は、低彩度側と高彩度領域で領域サイズが異なる空間領域を示す説明図
【図１１】本発明によるカラー補正の正解率の説明表
【図１２】色相彩度座標系における複数の放射状領域の説明図
【図１３】別実施形態を示すニューラルネットワークで構成される第一補正量演算手段の説明図
【符号の説明】
【０１０６】
１０：補正量演算手段
１１：補正量割付手段
１２：領域補正量演算手段
１３：揃い補正量演算手段
１４：適正補正量演算手段
１５：カラー補正手段
１００：平均画素値演算手段
１０１：座標変換手段
１０２：円領域特性値演算手段
１０３：カラー補正量演算手段
１１０：領域サイズ調整手段
【図１】
0005168809.tif 000017
【図２】
0005168809.tif 000018
【図３】
0005168809.tif 000019
【図４】
0005168809.tif 000020
【図５】
0005168809.tif 000021
【図６】
0005168809.tif 000022
【図７】
0005168809.tif 000023
【図８】
0005168809.tif 000024
【図９】
0005168809.tif 000025
【図１０】
0005168809.tif 000026
【図１１】
0005168809.tif 000027
【図１２】
0005168809.tif 000028
【図１３】
0005168809.tif 000029