(57)【要約】  一列に並ぶ画素のうち注目画素を中心として一方向及び他方向に並ぶ画素の画素値を夫々色毎に合計する合計画素値算出部（１０）と、合計画素値算出部による算出結果又は注目画素の画素値の少なくとも一方に基づいて、注目画素が軸上色収差の発生起点であるかを判定し発生起点を検出する発生起点検出部（１１）と、発生起点を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定部（１２）と、発生起点に対して合計画素値算出部の算出結果に基づいて、一方向及び他方向につき特定色空間上の色空間情報を夫々算出する色空間情報算出部（１５）と、一方向の色空間情報と他方向の色空間情報との差分を算出する色空間差分算出部（１６）と、差分に応じて軸上色収差補正の補正量を算出する補正量算出部（１７）と、補正量を用いて軸上色収差領域について補正する補正部（４）とを備える画像処理装置である。
2016181476.tif 000002
【特許請求の範囲】
【請求項１】
入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出部と、
該合計画素値算出部による算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素が軸上色収差の発生起点であるか否かを判定し、発生起点を検出する発生起点検出部と、
前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定部と、
前記発生起点として検出された前記注目画素に対し、前記合計画素値算出部による算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出部と、
該色空間情報算出部により算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出部と、
該色空間差分算出部により算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出部と、
前記軸上色収差領域について前記補正量を用いて補正する補正部と、
を備える画像処理装置。
【請求項２】
前記色空間情報算出部が、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
前記色空間差分算出部が、前記色空間情報算出部による算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
前記補正量算出部が、前記色空間差分算出部による算出結果に基づいて軸上色収差補正の補正量を算出する請求項１記載の画像処理装置。
【請求項３】
前記補正量算出部が、前記色の角度の差分が１８０度に近いほど大きな補正量を算出し、１８０度から離れている程小さな補正量を算出する請求項２記載の画像処理装置。
【請求項４】
前記補正量算出部が、前記色の強度の差分の絶対値が大きいほど大きな補正量を算出し、絶対値が小さいほど小さな補正量を算出する請求項２記載の画像処理装置。
【請求項５】
前記補正量算出部が、前記色の明るさの差分の絶対値が大きいほど大きな補正量を算出し、絶対値が小さいほど小さな補正量を算出する請求項２記載の画像処理装置。
【請求項６】
前記補正量算出部が、前記色の角度が前記一方向と前記他方向で同方向となる補正量を算出する請求項２記載の画像処理装置。
【請求項７】
前記補正量算出部が、前記色の強度が前記一方向と前記他方向とで同一となる補正量を算出する請求項２記載の画像処理装置。
【請求項８】
前記補正量算出部が、前記色の明るさが前記一方向と前記他方向とで同一となる補正量を算出する請求項２記載の画像処理装置。
【請求項９】
前記特定色空間が、ＨＳＶ空間である請求項１乃至請求項７の何れか１項記載の画像処理装置。
【請求項１０】
前記特定色空間が、Ｌ＊ａ＊ｂ＊空間である請求項１乃至請求項７の何れか１項記載の画像処理装置。
【請求項１１】
前記発生起点検出部が、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項１乃至請求項１０の何れか１項記載の画像処理装置。
【請求項１２】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ２＜ＰｒＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計である。
【請求項１３】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｂＤ２
ここで、
ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
【請求項１４】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ２＜ＰｒＤ２
ＰｇＤ２＜ＰｂＤ２
ここで、
ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
【請求項１５】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｒＤ２
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
【請求項１６】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｒＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、である。
【請求項１７】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計
である。
【請求項１８】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＋ＴＨＲ１＞ＰｒＤ１
ＰｇＤ２−−ＴＨＲ２＜ＰｒＤ２
ＰｇＤ１＋ＴＨＢ１＞ＰｂＤ１
ＰｇＤ２−−ＴＨＢ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２は任意の定数である。
【請求項１９】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＋ＴＨＲ＞ＰｒＤ１
ＰｇＤ２＋ＴＨＲ＜ＰｒＤ２
ＰｇＤ１＋ＴＨＢ＞ＰｂＤ１
ＰｇＤ２＋ＴＨＢ＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ，ＴＨＢは任意の定数である。
【請求項２０】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１から請求項１０の何れか１項記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１×ＫＲ
ＰｇＤ２＜ＰｒＤ２×ＫＲ
ＰｇＤ１＞ＰｂＤ１×ＫＢ
ＰｇＤ２＜ＰｂＤ２×ＫＢ
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＫＲ，ＫＢは任意の定数である。
【請求項２１】
前記画像が、各画素にＲＧＢ全ての画素値を有する請求項１から２０の何れか１項記載の画像処理装置。
【請求項２２】
前記画像が、各画素においてＲＧＢの少なくとも１つの画素値が欠落した画素値を有する請求項１から請求項２０の何れか１項記載の画像処理装置。
【請求項２３】
入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、
該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、
前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、
前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、
該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、
該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、
前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、
を備える画像処理方法。
【請求項２４】
前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出する請求項２３記載の画像処理方法。
【請求項２５】
前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項２３又は請求項２４記載の画像処理方法。
【請求項２６】
入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、
該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、
前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、
前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、
該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、
該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、
前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、
をコンピュータに実行させる画像処理プログラム。
【請求項２７】
前記コンピュータが、
前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出する請求項２７記載の画像処理プログラム。
【請求項２８】
前記コンピュータが、
前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項２６又は請求項２７記載の画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置、画像処理方法及び画像処理プログラムに関するものである。
【背景技術】
【０００２】
通常の光学系を通して撮影された画像には軸上色収差が発生することが知られている。ここで、軸上色収差とは、光学系、撮像系及び画像処理系に起因し、光軸上の焦点位置が波長によって異なるために生じる偽色、すなわち本来の被写体には生じない色付きをいう。そして、このような軸上色収差を画像処理で補正するための技術が提案されている（例えば、特許文献１参照。）。
特許文献１では、画像データ中の注目画素について、偽色可能性を示すパラメータとして、白とびへの近さ、彩度の大きさ及び特定色相への近さを評価し、これらの少なくともいずれかを適用して軸上色収差を補正している。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特許第４５３９２７８号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、軸上色収差を補正するに際し、特許文献１のように注目画素について評価し、これに基づいて補正を行うだけでは不十分である。すなわち、いくつかの注目画素が本来同一色であったとしても、当該注目画素に隣接する画素との関係によって、人間にはその視覚特性上それら注目画素が互いに異なる色であるように見える。このため、人間の視覚特性に鑑みると、注目画素のみを評価して軸上色収差を補正するだけでは不十分となる。【０００５】
本発明は、上述した事情に鑑みてなされたものであって、人間の視覚特性に基づいて軸上色収差を精度よく補正することを目的としている。
【課題を解決するための手段】
【０００６】
上記目的を達成するため、本発明は以下の手段を提供する。
本発明の一態様は、入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出部と、該合計画素値算出部による算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素が軸上色収差の発生起点であるか否かを判定し、発生起点を検出する発生起点検出部と、前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定部と、前記発生起点として検出された前記注目画素に対し、前記合計画素値算出部による算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出部と、該色空間情報算出部により算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出部と、該色空間差分算出部により算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出部と、前記軸上色収差領域について前記補正量を用いて補正する補正部と、を備える画像処理装置である。【０００７】
本態様によれば、注目画素を挟む両側に一直線上に配列されている画素の画素値の合計を合計画素値算出部により算出し、この算出結果又は注目画素の画素値に基づいて発生起点検出部により注目画素が軸上色収差の発生起点であるか判定し、発生起点を検出する。注目画素が軸上色収差の発生起点であるとして検出された場合、色空間情報算出部により、注目画素について、特定色空間上の色空間情報を一方向及び他方向について夫々算出し、色空間差分算出部によりこれらの差分を算出する。そして、この差分に基づいて算出された補正量を用いて、注目画素を中心とする周辺領域からなる軸上色収差領域について補正処理を行う。このように、注目画素だけでなく、注目画素を挟んで両側に一直線上に配列されている画素の画素値をも評価した補正量を算出することで、人間の視覚特性に基づいて軸上色収差補正を行うことができる。【０００８】
上記態様においては、前記色空間情報算出部が、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、前記色空間差分算出部が、前記色空間情報算出部による算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、前記補正量算出部が、前記色空間差分算出部による算出結果に基づいて軸上色収差補正の補正量を算出することとしてもよい。
このようにすることで、人間がその視覚特性上、軸上色収差を知覚しやすい隣接する領域の色の角度の差、色の強度の差、明るさの差に応じた最適な補正量を算出することができる。【０００９】
また上記様態においては、前記補正量算出部が、前記色の角度の差分が１８０度に近いほど大きな補正量を算出し、１８０度から離れている程小さな補正量を算出することとしてもよい。
このようにすることで、色空間の角度を色相環のように０度から３６０度で表した場合に、２つの色が対極となる１８０度のように、人間の視覚的に色付きが検知されやすい場合に適切な補正量を算出することができる。【００１０】
また上記様態においては、前記補正量算出部が、前記色の強度の差分の絶対値が大きいほど大きな補正量を算出し、絶対値が小さいほど小さな補正量を算出することとしてもよい。
このようにすることで、色の強度の差分が大きく、人間の視覚的に色付きが検知されやすい場合に適切な補正量を算出することができる。【００１１】
また上記様態においては、前記補正量算出部が、前記色の明るさの差分の絶対値が大きいほど大きな補正量を算出し、絶対値が小さいほど小さな補正量を算出することとしてもよい。
このようにすることで、色の明るさの差分が大きく、人間の視覚的に色付きが検知されやすい場合にも適切な補正量を算出することができる。【００１２】
また上記様態においては、前記補正量算出部が、前記色の角度が前記一方向と前記他方向で同方向となる補正量を算出することとしてもよい。
このようにすることで、補正後に色の角度の差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。【００１３】
また上記様態においては、前記補正量算出部が、前記色の強度が前記一方向と前記他方向とで同一となる補正量を算出することとしてもよい。
このようにすることで、補正後に色の強度の差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。【００１４】
また上記様態においては、前記補正量算出部が、前記色の明るさが前記一方向と前記他方向とで同一となる補正量を算出することとしてもよい。
このようにすることで、補正後に色の明るさの差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。【００１５】
また上記様態においては、前記特定色空間が、ＨＳＶ空間であることとしてもよい。
このようにすることで、比較的に低い計算コストで色空間の算出を行うことができる。【００１６】
また上記様態においては、前記特定色空間が、Ｌ＊ａ＊ｂ＊空間であることとしてもよい。
このようにすることで、人間の視覚的な感覚に近い色空間の算出を行うことができる。【００１７】
また、上記態様においては、前記発生起点検出部が、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出することとしてもよい。
このようにすることで、近飽和領域、近飽和以外の領域の少なくとも一方を軸上色収差の発生起点として検出することができる。【００１８】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ２＜ＰｒＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計である。
このようにすることで、フォーカス位置近傍における軸上色収差の発生起点を検出することができる。【００１９】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
このようにすることによっても、フォーカス位置近傍における軸上色収差の発生起点を検出することができる。【００２０】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ２＜ＰｒＤ２
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
このようにすることで、他方向にパープルフリンジが発生し、一方向のＲの画素値の合計が、Ｇの画素値の合計より小さく、色相環で反対色に近い色となってパープルフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。【００２１】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｒＤ２
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
このようにすることによっても、他方向にパープルフリンジが発生し、一方向のＢの画素値の合計が、Ｇの画素値の合計より小さく、色相環で反対色に近い色となってパープルフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。【００２２】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｒＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計である。【００２３】
このようにすることで、一方向にグリーンフリンジが発生し、他方向のＲの画素値の合計がＧの画素値の合計より大きく、色相環で反対色に近い色となってグリーンフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。【００２４】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。【００２５】
このようにすることによっても、一方向にグリーンフリンジが発生し、他方向のＢの画素値の合計がＧの画素値の合計より大きく、色相環で反対色に近い色となってグリーンフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。【００２６】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＋ＴＨＲ１＞ＰｒＤ１
ＰｇＤ２−ＴＨＲ２＜ＰｒＤ２
ＰｇＤ１＋ＴＨＢ１＞ＰｂＤ１
ＰｇＤ２−ＴＨＢ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２は任意の定数である。
このようにすることで、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２を適当に選択することにより、ノイズの影響を低減し、軸上色収差の発生起点を安定して検出することができる。【００２７】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＋ＴＨＲ＞ＰｒＤ１
ＰｇＤ２＋ＴＨＲ＜ＰｒＤ２
ＰｇＤ１＋ＴＨＢ＞ＰｂＤ１
ＰｇＤ２＋ＴＨＢ＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ，ＴＨＢは任意の定数である。
このようにすることで、ＴＨＲ，ＴＨＢを適当に選択することにより、無彩色ではなく、色付き被写体上に発生した軸上色収差の発生起点を検出することができる。【００２８】
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１×ＫＲ
ＰｇＤ２＜ＰｒＤ２×ＫＲ
ＰｇＤ１＞ＰｂＤ１×ＫＢ
ＰｇＤ２＜ＰｂＤ２×ＫＢ
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＫＲ，ＫＢは任意の定数である。
このようにすることによっても、ＫＲ，ＫＢを適当に選択することにより、無彩色ではなく、色付き被写体上に発生した軸上色収差の発生起点を検出することができる。【００２９】
また、上記態様においては、前記画像が、各画素にＲＧＢ全ての画素値を有していてもよい。
このようにすることで、各画素においてＲＧＢ全ての画素値を同時に取得する同時化された画像に対しても軸上色収差の発生起点を検出する処理を行うことができる。【００３０】
また、上記態様においては、前記画像が、各画素においてＲＧＢの少なくとも１つの画素値が欠落した画素値を有していてもよい。
このようにすることで、各画素においてＲＧＢ全ての画素値が同時に取得されない、同時化されていない画像に対しても軸上色収差の発生起点を検出することができ、検出時に同時化のための補間エラーの影響をなくすことができる。【００３１】
また、本発明の他の態様は、入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、を備える画像処理方法である。【００３２】
上記他の態様において、前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出することとしてもよい。【００３３】
また、上記他の態様において、前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出することとしてもよい。【００３４】
また、本発明の他の態様は、入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、をコンピュータに実行させる画像処理プログラムである。【００３５】
上記他の態様において、前記コンピュータが、前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出することとしてもよい。【００３６】
また、上記他の態様において、前記コンピュータが、前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出することとしてもよい。
【発明の効果】
【００３７】
本発明によれば、近飽和領域もしくは近飽和領域以外の領域少なくとも一方においても軸上色収差を検出して人間の視覚特性に基づいて補正することができるという効果を奏する。
【図面の簡単な説明】
【００３８】
【図１】本発明の一実施形態に係る画像処理装置を示すブロック図である。
【図２】図１の画像処理装置を適用する光学系と軸上色収差との関係を模式的に示す図である。
【図３】図１の画像処理装置を適用する光学系と軸上色収差との関係を模式的に示す図である。
【図４】図１の画像処理装置を用いた画像処理方法を説明するフローチャートである。
【図５】図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図である。
【図６】図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図であり、特に、補正前の注目画素とその周辺画素との座標関係を示す。
【図７】図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図であり、特に、補正後の注目画素とその周辺画素との座標関係を示す。
【図８】図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図であり、特に、注目画素の画素位置がＤ１方向として定まっている場合の補正前の注目画素とその周辺画素との座標関係を示す。
【図９】図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図であり、特に、注目画素の画素位置がＤ１方向として定まっている場合の補正後の注目画素とその周辺画素との座標関係を示す。
【図１０】図２とは異なる態様で軸上色収差が発生している場合のステップエッジの第１の例を示す図である。
【図１１】図２とは異なる態様で軸上色収差が発生している場合のステップエッジの第２の例を示す図である。
【図１２】画像信号にノイズが含まれる場合の対応例をノイズのない場合及び対応前の状態との比較で示す図である。
【図１３】有彩色のステップエッジの場合の対応例を無彩色の場合及び対応前の状態との比較で示す図である。
【図１４】有彩色のステップエッジの場合の他の対応例を無彩色の場合及び対応前の状態との比較で示す図である。
【図１５】ベイヤ配列された撮像素子により取得された画像のＲを注目画素とする場合の画素配列を示す図である。
【図１６】ベイヤ配列された撮像素子により取得された画像のＢを注目画素とする場合の画素配列を示す図である。
【図１７】ベイヤ配列された撮像素子により取得された画像のＧＲを注目画素とする場合の画素配列を示す図である。
【図１８】ベイヤ配列された撮像素子により取得された画像のＧＢを注目画素とする場合の画素配列を示す図である。
【発明を実施するための形態】
【００３９】
本発明の一実施形態に係る画像処理装置及び画像処理方法について、図面を参照して以下に説明する。
ここで、本実施形態に係る画像処理装置は、カメラにより取得された入力画像に対して後述する画像処理を施すものである。そして、画像処理装置１の処理対象である入力画像は、図示しない撮像光学系を通して入射された光線が固体撮像素子において電気信号に変換され、固定ノイズパターン除去及びアナログゲイン調整等が行われた後に、Ａ／Ｄ変換器でディジタル信号に変換され、各種画像処理が施された画像信号である。この入力画像は、各画素においてＲＧＢ全ての信号値が得られるフルカラー画像であり、三板撮像素子により取得された三板カラー画像、もしくは、単板撮像素子により取得された単板カラー画像を各画素でＲＧＢ値が得られるように補間処理された画像である。
【００４０】
本実施形態に係る画像処理装置は、図１に示すように、入力画像において軸上色収差が発生している軸上色収差領域を抽出する領域抽出部２と、領域抽出部２により抽出された軸上色収差領域に対する補正量を算出する軸上色収差補正量算出部３と、軸上色収差領域の軸上色収差を補正した出力画像を出力する収差補正部４とを備えている。
【００４１】
領域抽出部２は、合計画素値算出部１０、発生起点検出部１１、及び領域判定部１２を備えている。
合計画素値算出部１０は、入力画像中の注目画素を設定し、設定された注目画素を含んで一列に配列している複数の画素について、注目画素を中心として一方向（Ｄ１）及び他方向（Ｄ２）に配列している画素の画素値を夫々色毎に合計する。
【００４２】
発生起点検出部１１は、注目画素の画素値、又は合計画素値算出部１１により算出された注目画素を挟んで一方向に並ぶ画素の合計画素値若しくは他方向に並ぶ画素の合計画素値の少なくとも一方に基づいて、注目画素が軸上色収差の発生起点であるかを判定し、発生起点となる注目画素を検出する。
【００４３】
発生起点検出部５により検出する軸上色収差の発生起点とは、図２に示すように、例えば、注目画素のＧ画素値の大きさが飽和レベルになる注目画素、又は図３に示すように、注目画素を挟んで一直線に並ぶ一方向及び他方向の色が異なる場合の注目画素の少なくとも一方を意味している。
すなわち、図２及び図３に示すように、通常の光学系７では、波長の異なる光は光軸の奥行き方向に結像位置が異なることが知られている。一般にＧに対して長波長のＲは光軸方向の奥側に結像し、Ｇに対して短波長のＢは光軸方向の手前側に結像する。このＲＧＢの結像位置のずれが軸上色収差となって現れる。
【００４４】
無彩色のステップエッジを撮影した場合、像面ではＧがステップエッジの形状を有するが、Ｒ，Ｂは結像せずにステップエッジの形状が鈍っている。一方、像面奥では、Ｇは結像しないためステップエッジの形状が鈍っているが、Ｒが結像してステップエッジの形状を保持している。ＢはＧよりさらに鈍った形状となる。
反対に、像面手前では、Ｇは結像しないためステップエッジの形状が鈍り、Ｂが結像してステップエッジの形状を保持している。ＲはＧよりさらに鈍った形状となる。
【００４５】
領域判定部１２は、発生起点検出部５により検出された発生起点である注目画素の所定の周辺領域を軸上色収差領域として判定する。すなわち、注目画素が軸上色収差の発生起点として検出された場合には、軸上色収差の補正対象となるので、注目画素と注目画素を中心とする所定の周辺画素からなる領域を軸上色収差領域として判定する。軸上色収差領域としては、例えば、注目画素を中心とする５×５画素の領域や７×７画素の領域等、適宜設定することができる。
【００４６】
軸上色収差補正量算出部３は、色空間情報算出部１５と、色空間差分算出部１６と、補正量算出部１７とを備えている。
色空間情報算出部１５は、発生起点検出部５により発生起点であると検出された注目画素について、特定色空間上の色空間情報を、注目画素に対する一方向及び他方向の夫々について算出する。特定色空間としては、例えば、ＨＳＶ空間やｌ＊ａ＊ｂ空間などを適用することができる。そして、色空間情報は、合計画素値算出部１０により算出された、注目画素を中心として一方向（Ｄ１）及び他方向（Ｄ２）に配列している画素についての色毎の画素値の合計に基づいて算出する。
【００４７】
色空間差分算出部１６は、色空間情報算出部１５により算出された一方向の色空間情報と他方向の色空間情報との差分を算出する。
補正量算出部１７は、色空間差分算出部１６により算出された差分に応じて軸上色収差領域に対する補正量を算出する。
収差補正部４は、補正量算出部１７により算出された補正量を用いて軸上色収差領域を補正する。
【００４８】
以下に、このように構成された本実施形態に係る画像処理装置の作用としての画像処理方法について、図４のフローチャートを参照して説明する。
ステップＳ１において、合計画素値算出部１０は、図５に示すように、入力画像中に注目画素（０，０）を設定し、次のステップＳ２（合計画素値算出ステップ）において、この注目画素を中心として、図５中の縦方向、横方向、右斜め方向及び左斜め方向の４方向について、注目画素を挟んで両側に並ぶ複数画素の画素値の合計を夫々算出する（ステップＳ２）。
具体的には、数１に示す式に従って、各方向の色毎の画素値の合計を算出する。
【００４９】
【数１】
2016181476.tif 000003
【００５０】
ここで、Ｐｃ（ｘ，ｙ）は、座標（ｘ，ｙ）に位置する画素の画素値、ＰｃＬは注目画素の左側の画素列の画素値の合計、ＰｃＲは注目画素の右側の画素列の画素値の合計、ＰｃＴは注目画素の上側の画素列の画素値の合計、ＰｃＢは注目画素の下側の画素列の画素値の合計、ＰｃＬＴは注目画素の左上側の画素列の画素値の合計、ＰｃＬＢは注目画素の左下側の画素列の画素値の合計、ＰｃＲＴは注目画素の右上側の画素列の画素値の合計、ＰｃＲＢは注目画素の右下側の画素列の画素値の合計である。また、添字ｃはカラーを意味し、ｒ、ｇ、ｂに置き換えることができ、置き換えることにより色毎の画素値の合計を算出する。また、ｉ，ｍは正の整数である。
【００５１】
ステップＳ３（発生起点検出ステップ）において、発生起点検出部１１は、以下の条件式（１）又は、合計画素値算出部１０により算出された画素値の合計が、条件式（２）〜（５）の全てを満たす場合に、注目画素を軸上色収差の発生起点として検出する。
条件式（１）は、注目画素のＧ画素値の大きさが任意の飽和レベル閾値ＴＨｓａｔより大きい場合に、注目画素を軸上色収差の発生起点として検出することを示す。
【００５２】
Ｐｇ（０，０）＞ＴＨｓａｔ      …（１）
ＰｇＤ１＞ＰｒＤ１                  …（２）
ＰｇＤ２＜ＰｒＤ２                  …（３）
ＰｇＤ１＞ＰｂＤ１                  …（４）
ＰｇＤ２＜ＰｂＤ２                  …（５）
【００５３】
ここで、ＰｇＤ１は注目画素から一方向Ｄ１に配列している画素のＧの画素値の合計、ＰｇＤ２は注目画素から他方向Ｄ２に配列している画素のＧの画素値の合計、ＰｒＤ１は注目画素から一方向Ｄ１に配列している画素のＲの画素値の合計、ＰｒＤ２は注目画素から他方向Ｄ２に配列している画素のＲの画素値の合計、ＰｂＤ１は注目画素から一方向Ｄ１に配列している画素のＢの画素値の合計、ＰｂＤ２は注目画素から他方向Ｄ２に配列している画素のＢの画素値の合計である。
【００５４】
添字のＤ１は注目画素からの一方向であることを示し、例えば、数１におけるＬ、Ｔ、ＬＴ、ＬＢが該当する。また、添字のＤ２は注目画素からの他方向であることを示し、例えば、Ｒ、Ｂ、ＲＴ、ＲＢが該当する。方向の定義は任意であるが、Ｄ１としてＬ、Ｔ、ＬＴ、ＬＢを選択した場合にはＤ２はＲ、Ｂ、ＲＴ、ＲＢとなる。
【００５５】
なお、条件式（１）において、注目画素と任意の飽和レベル閾値ＴＨｓａｔとを比較する際に、注目画素のＧ画素のみを用いることとしたが、必ずしもこれに限られない。例えば、飽和レベル閾値と輝度値とを比較することもできる他、その他公知の方法を用いることもできる。また、近飽和以外の領域の軸上色収差の発生起点を算出する方法についても上記した例に限られず、公知の方法を用いることができる。
【００５６】
ステップＳ４（検出情報生成ステップ）において、発生起点検出部１１は、軸上色収差の発生起点であると検出された注目画素について「０」、発生起点ではないと検出された画素について「１」とする検出情報を生成する。発生起点検出部１１は、注目画素が、軸上色収差の発生起点であると検出された場合、その画素位置に対応してＰｒＤ１，ＰｒＤ２，ＰｇＤ１，ＰｇＤ２，ＰｂＤ１，ＰｂＤ２を保持する。
【００５７】
また、注目画素が、軸上色収差の発生起点であると検出された場合には、その画素位置と対応させて、図５中の縦方向、横方向、右斜め方向及び左斜め方向の４方向の何れの方向が上記条件式（２）〜（５）を満たしたのかを保持し、その方向についての一方向Ｄ１及び他方向Ｄ２を保持する。
つまり、発生起点検出部１１は、検出情報に関連する情報として、Ｄ１が、Ｌ、Ｔ、ＬＴ、ＬＢの何れであるのか、及びＤ２がＲ、Ｂ、ＲＴ、ＲＢの何れであるのかを保持しておく。
【００５８】
ステップＳ５（領域判定ステップ）において、領域判定部１２は、注目画素が軸上色収差の発生起点であると検出された場合、すなわち、注目画素について発生起点検出部１１により生成された検出情報が「０」であった場合、当該注目画素とその所定の周辺画素を含む領域を軸上色収差領域と判定する。軸上色収差領域として、例えば、検出情報が「０」である画素を中心とした７×７画素を含む領域の等、適宜設定することができる。
【００５９】
ステップＳ６（色空間情報算出ステップ）において、色空間情報算出部１５は、検出情報が「０」であった画素位置に対応するＰｒＤ１，ＰｇＤ１，ＰｂＤ１及びＰｒＤ２，ＰｇＤ２，ＰｂＤ２から特定色空間の色空間情報を算出する。
例えば、特定色空間をＨＳＶ空間とすると、ＰｒＤ１，ＰｇＤ１，ＰｂＤ１からＰｈＤ１，ＰｓＤ１，ＰｖＤ１の少なくとも１つを、ＰｒＤ２，ＰｇＤ２，ＰｂＤ２からＰｈＤ２，ＰｓＤ２，ＰｖＤ２の内Ｄ１方向と同じ色空間情報を算出する。
【００６０】
なお、ＰｈＤ１，ＰｈＤ２は、発生起点であると検出された注目画素におけるＤ１方向、Ｄ２方向についての角度（色相）に関する色空間情報を示し、ＰｓＤ１，ＰｓＤ２は、発生起点であると検出された注目画素におけるＤ１方向、Ｄ２方向についての強度（彩度）に関する色空間情報を示し、ＰｖＤ１，ＰｖＤ２は、発生起点であると検出された注目画素におけるＤ１方向、Ｄ２方向についての明るさ（明度）に関する色空間情報を示す。
【００６１】
また、特定色空間をＬ＊ａ＊ｂ＊空間とすると、ＰｒＤ１，ＰｇＤ１，ＰｂＤ１からＰｌ＊Ｄ１，Ｐａ＊Ｄ１，Ｐｂ＊Ｄ１の少なくとも１つを、ＰｒＤ２，ＰｇＤ２，ＰｂＤ２からＰｌ＊Ｄ２，Ｐａ＊Ｄ２，Ｐｂ＊Ｄ２の内Ｄ１方向と同じ色空間情報を算出する。なお、Ｐｌ＊Ｄ１，Ｐｌ＊Ｄ２は、発生起点であると検出された注目画素におけるＤ１方向、Ｄ２方向についての明るさ（明度）に関する色空間情報を示し、Ｐａ＊Ｄ１，Ｐｂ＊Ｄ１，Ｐａ＊Ｄ２，Ｐｂ＊Ｄ２は、補色に関する情報であって、角度及び強度に関する色空間情報として利用することができる。
【００６２】
ステップＳ７（色空間差分算出ステップ）において、色空間差分算出部１６により、２つの方向の色空間情報から色空間の差分を算出する。色区間情報算出部１６は、差分として、色の角度の差分、色の強度の差分、明るさの差分の少なくとも１つを算出する。
【００６３】
ステップＳ８（補正量算出ステップ）では、色空間差分算出部１６により算出された差分に基づいて、補正量を算出し、次のステップＳ９において、全ての画素について処理が行われたか判定し、全ての画素についての処理が完了していない場合には、ステップＳ１に戻り、処理が未だ行われていない画素を注目画素としてステップＳ２〜Ｓ７の工程が繰り返される。
【００６４】
全ての画素について処理が行われた場合には、ステップＳ１０（収差補正ステップ）において、領域判定部１２により設定された軸上色収差領域について、補正量算出部１７により算出された補正量を用いて収差補正部４により補正処理を実施し、処理を終了する。
【００６５】
上述したように、色空間情報として、色についての、角度、強度及び明るさの少なくとも一つを算出し、これらの色空間情報から角度の差分、強度の差分、明るさの差分の少なくとも一つを算出し、これら差分のうち何れか一又は複数を組み合わせた情報に基づいて補正量を算出し、軸上色収差領域の補正を行うことができる。
以下、角度の差分、強度の差分及び明るさの差分のうち、何れか１つまたは複数を組み合わせて補正量を算出する場合の具体例について説明する。
【００６６】
（角度に関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、角度に関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した角度に関する色空間情報からその差分を算出する（ステップＳ７）。
【００６７】
具体的には、特定色空間がＨＳＶ空間である場合、角度に関する色空間情報として、ＰｈＤ１，ＰｈＤ２を算出し、これらの差分の絶対値Ｐａｎｇｌｅを以下の式（６）に従って算出する。
Pangle=|PhD1-PhD2|      …（６）
【００６８】
また、特定色空間がＬ＊ａ＊ｂ＊空間である場合、角度に関する色空間情報として、Ｐａ＊Ｄ１，Ｐｂ＊Ｄ１，Ｐａ＊Ｄ２，Ｐｂ＊Ｄ２を算出し、これらの差分の絶対値Ｐａｎｇｌｅを以下の式（７）に従って算出する。
Pangle=|tan-1(Pb*D1/Pa*D1)-
tan-1(Pb*D2/Pa*D2)|      …（７）
【００６９】
ここで、発生起点検出部１１において、条件式（１）を満たすことにより注目画素のＧ画素値の大きさが任意の飽和レベル閾値ＴＨｓａｔより大きいとして、当該注目画素が発生起点と判定された場合は、特定方向のＤ１，Ｄ２が定まっていない。この場合には、例えば、上述した縦、横、斜めの計４つの全ての方向のＤ１，Ｄ２についてＰａｎｇｌｅを算出し、Ｐａｎｇｌｅが最も大きい値を採用することができる。また、全ての方向のＤ１，Ｄ２についてＰａｎｇｌｅを算出し、これらの平均値を採用することもできる。
【００７０】
補正量算出部１７は、Ｐａｎｇｌｅから補正量Ｃａを算出する。補正量Ｃａは、例えば、Ｐａｎｇｌｅが０から３６０度の値の場合、Ｐａｎｇｌｅが１８０度に近い程補正量を大きくし、Ｐａｎｇｌｅが１８０度から離れている程補正量を小さくする。色空間の角度を色相環のように０度から３６０度で表した場合に、２つの色が対極となる１８０度のように、人間の視覚的に色付きが検知されやすい場合に適切な補正量を算出することができる。色の角度が一方向と他方向で同方向となるような補正量を算出することで、補正後に色の角度の差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。
【００７１】
そして、全ての画素について補正量が算出されると、角度に関する色空間情報から算出された補正量に基づく収差補正部の補正処理は、例えば、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’として以下の式（８）〜（１０）に従って行われる。
【００７２】
Pr’=Pr-Ca       …（８）
Pg’=Pg          …（９）
Pb’=Pb-Ca      …（１０）
【００７３】
また、補正前の画素値の角度をＰａ、補正後の画素値の角度をＰａ’として、以下の式（１１）に従って補正を行うこともできる。
Pa’=Pa-Ca    …（１１）
【００７４】
以下、角度に関する色空間情報に基づいて補正を行う場合の他の例について説明する。
色空間情報算出部１５は、例えば、特定色空間がＨＳＶ空間である場合には、角度に関する色空間情報としてのＰａｎｇｌｅ１，Ｐａｎｇｌｅ２を以下の式（１２）及び式（１３）に従って算出する。
Pangle1=PhD1      …（１２）
Pangle2=PhD2      …（１３）
【００７５】
また、特定色空間がＬ＊ａ＊ｂ＊空間である場合、角度に関する色空間情報としてＰａｎｇｌｅ１，Ｐａｎｇｌｅ２を以下の式（１４）及び式（１５）に従って算出する。
Pangle1=tan-1(Pb*D1/Pa*D1)      …（１４）
Pangle2=tan-1(Pb*D2/Pa*D2)      …（１５）
【００７６】
ここで、発生起点検出部１１において、条件式（１）を満たすことにより注目画素が発生起点と判定された場合には、特定方向のＤ１，Ｄ２が定まっていない。この場合には、色空間情報算出部１５は、例えば、上述した縦、横、斜めの計４つの全ての方向のＤ１，Ｄ２でＰａｎｇｌｅ１，Ｐａｎｇｌｅ２を算出し、色空間差分算出部１６は、Ｐａｎｇｌｅ１とＰａｎｇｌｅ２の差の絶対値の最大値を、角度に関する色空間情報の差分として採用し、最大の絶対値を示すＤ１，Ｄ２を特定方向とすることができる。
例えば、角度に関する色空間情報の差分の絶対値Ｐａｎｇｌｅは式（１６）により算出する。
Pangle=|Pangle1-Pangle2|      …（１６）
【００７７】
補正量算出部１７はＰａｎｇｌｅ，Ｐａｎｇｌｅ１，Ｐａｎｇｌｅ２から補正量Ｃａ１，Ｃａ２を、例えば、以下の式（１７）〜（２０）に従って算出する。
Ｐａｎｇｌｅ１がＰａｎｇｌｅ２より大きい場合
Ca1=Pangle      …（１７）
Ca2=0           …（１８）
Ｐａｎｇｌｅ２がＰａｎｇｌｅ１より大きい場合
Ca1=0           …（１９）
Ca2=Pangle      …（２０）
【００７８】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度に基づく補正量を用いて収差補正部４により補正処理を実施する。
収差補正部４による角度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（２１）〜（２３）で表すことができる。
Pr1’=Pr-Ca1      …（２１）
Pg1’=Pg          …（２２）
Pb1’=Pb-Ca1      …（２３）
【００７９】
なお、補正前の画素値の角度をＰａｎｇｌｅ１、補正後の画素値の角度をＰａｎｇｌｅ１’とし、次の式（２４）に従って補正処理を行うこともできる。。
Pangle1’=Pangle1-Ca1      …（２４）
同様に、補正対象の画素位置がＤ２方向の場合、補正前の画素値をＰｒ２，Ｐｇ２，Ｐｂ２、補正後の画素値をＰｒ２’，Ｐｇ２’，Ｐｂ２’とし、次の式（２５）〜（２７）に従って補正処理を行うことができる。
Pr2’=Pr-Ca2      …（２５）
Pg2’=Pg          …（２６）
Pb2’=Pb-Ca2      …（２７）
【００８０】
また、補正前の画素値の角度をＰａｎｇｌｅ２、補正後の画素値の角度をＰａａｎｇｌｅ２’として、次の式（２８）で補正処理を行うこともできる。
Pangle2’=Pangle2-Ca2      …（２８）
【００８１】
角度に基づく色空間情報に従って軸上色収差を補正する場合、画像処理等により色空間上で特定の色について角度の差が広がるような処理が行われる場合に特に有効である。
【００８２】
（強度に関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、強度に関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した強度に関する色空間情報からその差分を算出する（ステップＳ７）。
具体的には、特定色空間がＨＳＶ空間である場合、強度に関する色空間情報として、ＰｓＤ１，ＰｓＤ２を算出し、これらの差分の絶対値Ｐｓｔｒｅｎｇｔｈを以下の式（２９）に従って算出する。
Pstrength=|PsD1-PsD2|      …（２９）
また、特定色空間がＬ＊ａ＊ｂ＊空間である場合、強度に関する色空間情報として、Ｐａ＊Ｄ１，Ｐｂ＊Ｄ１，Ｐａ＊Ｄ２，Ｐｂ＊Ｄ２を算出し、これらの差分の絶対値Ｐｓｔｒｅｎｇｔｈを以下の式（３０）に従って算出する。
Pstrength=|√(Pa*D1^2+Pb*D1^2)-√(Pa*D2^2+Pb*D2^2)|…（３０）
【００８３】
発生起点検出部１１において、条件式（１）を満たすことにより注目画素が発生起点と判定され、特定方向のＤ１，Ｄ２が定まっていない場合には、例えば、上述した縦、横、斜めの計４つの全ての方向のＤ１，Ｄ２でＰｓｔｒｅｎｇｔｈを算出し、Ｐｓｔｒｅｎｇｔｈの最大値を採用し、その最大値を示すＤ１，Ｄ２を特定方向とすることができる。
補正量算出部１７は、Ｐｓｔｒｅｎｇｔｈから以下の式（３１）に従って、補正量Ｃｂを算出する。Ｋは任意の定数とする。
Cb=K*Pstrength      …（３１）
【００８４】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部により算出された強度に基づく補正量を用いて収差補正部４により補正処理を実施する。収差補正部４の強度に基づく補正処理は、例えば、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’とし、以下の式（３２）〜（３４）で表すことができる。
Pr’=Pr-Cb      …（３２）
Pg’=Pg         …（３３）
Pb’=Pb-Cb      …（３４）
【００８５】
また、補正前の画素値の強度をＰｓ、補正後の画素値の強度をＰｓ’とし、次の式（３５）に従って補正処理を行うこともできる。
Ps’=Ps-Cb      …（３５）
【００８６】
色の強度の差分の絶対値が大きいほど大きな補正量を算出し、絶対値が小さいほど小さな補正量を算出することで、色の強度の差分が大きく、人間の視覚的に色付きが検知されやすい場合に適切な補正量を算出することができる。
また、色の強度が一方向と他方向とで同一となる補正量を算出することで、補正後に色の強度の差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。
【００８７】
続いて、強度に関する色空間情報に基づいて補正を行う場合の他の例について説明する。色空間情報算出部１５は、例えば、特定色空間がＨＳＶ空間である場合には、強度に関する色空間情報としてのＰｓｔｒｅｎｇｔｈ１，Ｐｓｔｒｅｎｇｔｈ２を以下の式（３６）及び式（３７）に従って算出する
Pstrength1=PsD1      …（３６）
Pstrength2=PsD2      …（３７）
また、特定色空間がＬ＊ａ＊ｂ＊空間である場合、強度に関する色空間情報としてＰｓｔｒｅｎｇｔｈ１，Ｐｓｔｒｅｎｇｔｈ２を以下の式（３８）及び式（３９）に従って算出する。
Pstrength1=√(Pa*D1^2+Pb*D1^2)      …（３８）
Pstrength2=√(Pa*D2^2+Pb*D2^2)      …（３９）
【００８８】
ここで、発生起点検出部１１において、条件式（１）を満たすことにより注目画素が発生起点と判定され、特定方向のＤ１，Ｄ２が定まっていない場合、例えば、色空間情報算出部１５が上述した縦、横、斜めの計４つの全ての方向のＤ１，Ｄ２でＰｓｔｒｅｎｇｔｈ１，Ｐｓｔｒｅｎｇｔｈ２を算出し、色空間差分算出部１６がＰｓｔｒｅｎｇｔｈ１とＰｓｔｒｅｎｇｔｈ２の差の絶対値の最大値を強度に関する色空間情報の差分として採用し、最大の絶対値を示すＤ１，Ｄ２を特定方向とすることができる。
例えば、強度に関する色空間情報の差分の絶対値Ｐｓｔｒｅｎｇｔｈは次の式（４０）により算出する。
Pstrength=|Pstrength1-Pstrength2|      …（４０）
【００８９】
補正量算出部１７はＰｓｔｒｅｎｇｔｈ，Ｐｓｔｒｅｎｇｔｈ１，Ｐｓｔｒｅｎｇｔｈ２から補正量Ｃｂ１，Ｃｂ２を、例えば、以下の式（４１）〜（４４）に従って算出する。
Ｐｓｔｒｅｎｇｔｈ１がＰｓｔｒｅｎｇｔｈ２より大きい場合
Cb1=Pstrength    …（４１）
Cb2=0            …（４２）
Ｐｓｔｒｅｎｇｔｈ２がＰｓｔｒｅｎｇｔｈ１より大きい場合
Cb1=0            …（４３）
Cb2=Pstrength    …（４４）
【００９０】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された強度に基づく補正量を用いて収差補正部４により補正処理を実施する。
【００９１】
収差補正部１７による強度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式次の式（４５）〜（４７）で表すことができる。
Pr1’=Pr-Cb1      …（４５）
Pg1’=Pg          …（４６）
Pb1’=Pb-Cb1      …（４７）
【００９２】
なお、補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ１、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ１’とし、次の式（４８）に従って補正処理を行うこともできる。
Pstrength1’=Pstrength1-Cb1    …（４８）
【００９３】
同様に、補正対象の画素位置がＤ２方向の場合、補正前の画素値をＰｒ２，Ｐｇ２，Ｐｂ２、補正後の画素値をＰｒ２’，Ｐｇ２’，Ｐｂ２’とし、次の式（４９）〜（５１）に従って補正処理を行うことができる。
Pr2’=Pr-Cb2    …（４９）
Pg2’=Pg        …（５０）
Pb2’=Pb-Cb2    …（５１）
【００９４】
また、補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ２、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ２’として、次の式（５２）で補正処理を行うこともできる。
Pstrength2’=Pstrength2-Cb2    …（５２）
【００９５】
強度に基づく色空間情報に従って軸上色収差を補正する場合、画像処理等により色空間上で特定の色の強度の差が広がるような処理が行われる場合に特に有効である。
【００９６】
（明度に関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、明るさに関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した明るさに関する色空間情報からその差分を算出する（ステップＳ７）。
【００９７】
具体的には、特定色空間がＨＳＶ空間である場合、明るさに関する色空間情報として、ＰｖＤ１，ＰｖＤ２を算出し、これらの差分の絶対値Ｐｌｉｇｈｔｎｅｓｓを以下の式（５３）に従って算出する。
色空間の２つの明るさの差の絶対値Ｐｌｉｇｈｔｎｅｓｓは次式で表される。
Plightness=|PvD1-PvD2|      …（５３）
【００９８】
また、特定色空間がＬ＊ａ＊ｂ＊空間である場合、強度に関する色空間情報として、Ｐｌ＊Ｄ１，Ｐｌ＊Ｄ２を算出し、これらの差分の絶対値Ｐｌｉｇｈｔｎｅｓｓを次の式（５４）に従って算出する。
Plightness=|Pl*D1-Pl*D2|      …（５４）
【００９９】
発生起点検出部１１において、条件式（１）を満たすことにより注目画素が発生起点と判定され、特定方向のＤ１，Ｄ２が定まっていない場合には、例えば、上述した縦、横、斜めの計４つの全ての方向のＤ１，Ｄ２でＰｌｉｇｈｔｎｅｓｓを算出し、Ｐｌｉｇｈｔｎｅｓｓの最大値を採用し、その最大値を示すＤ１，Ｄ２を特定方向とすることができる。
【０１００】
補正量算出部１７は、Ｐｌｉｇｈｔｎｅｓｓから以下の式に従って、補正量Ｃｃを算出する。Ｋは任意の定数とする。
Cc=K*Plightness      …（５５）
【０１０１】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部により算出された明るさに基づく補正量を用いて収差補正部４により補正処理を実施する。
収差補正部の補正処理は、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’とし、以下の式（５６）〜（５８）で表すことができる。
Pr’=Pr-Cc      …（５６）
Pg’=Pg         …（５７）
Pb’=Pb-Cc      …（５８）
【０１０２】
また、補正前の画素値の明度をＰｌ、補正後の画素値の強度をＰｌ’とし、次の式（５９）に従って補正処理を行うこともできる。
Pl’=Pl-Cc      …（５９）
【０１０３】
色の明るさの差分の絶対値が大きいほど大きな補正量を算出し、絶対値が小さいほど小さな補正量を算出することで、色の明るさの差分が大きく、人間の視覚的に色付きが検知されやすい場合にも適切な補正量を算出することができる。
また、色の明るさが一方向と他方向とで同一となるような補正量を算出することで、補正後に色の明るさの差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。
【０１０４】
続いて、明るさに関する色空間情報に基づいて補正を行う場合の他の例について説明する。色空間情報算出部１５は、例えば、特定色空間がＨＳＶ空間である場合には、明るさに関する色空間情報としてのＰｌｉｇｈｔｎｅｓｓ１，Ｐｌｉｇｈｔｎｅｓｓ２を以下の式（６０）及び式（６１）に従って算出する。
Plightness1=PvD1      …（６０）
Plightness2=PvD2      …（６１）
【０１０５】
また、特定色空間がＬ＊ａ＊ｂ＊空間である場合、明るさに関する色空間情報としてＰｌｉｇｈｔｎｅｓｓ１，Ｐｌｉｇｈｔｎｅｓｓ２を以下の式（６２）及び式（６３）に従って算出する。
Plightness1=Pl*D1      …（６２）
Plightnss2=Pl*D2       …（６３）
【０１０６】
ここで、発生起点検出部１１において、条件式（１）を満たすことにより注目画素が発生起点と判定され、特定方向のＤ１，Ｄ２が定まっていない場合、例えば、色空間情報算出部１５が上述した縦、横、斜めの計４つの全ての方向のＤ１，Ｄ２でＰｌｉｇｈｔｎｅｓｓ１，Ｐｌｉｇｈｔｎｅｓｓ２を算出し、色空間差分算出部１６がＰｌｉｇｈｔｎｅｓｓ１とＰｌｉｇｈｔｎｅｓｓ２の差分の絶対値の最大値を明るさに関する色空間情報の差分として採用し、最大の絶対値を示すＤ１，Ｄ２を特定方向とすることができる。
例えば、明るさに関する色空間情報差分の絶対値Ｐｌｉｇｈｔｎｅｓｓは次の式（６４）で表される。
Plightness=|Plightness1-Plightness2|      …（６４）
【０１０７】
補正量算出部１７はＰｌｉｇｈｔｎｅｓｓ，Ｐｌｉｇｈｔｎｅｓｓ１，Ｐｌｉｇｈｔｎｅｓｓ２から補正量Ｃｃ１，Ｃｃ２を、例えば、以下の式（６５）〜（６８）に従って算出する。
Ｐｌｉｇｈｔｎｅｓｓ１がＰｌｉｇｈｔｎｅｓｓ２より大きい場合
Cc1=Plightness      …（６５）
Cc2=0               …（６６）
Ｐｌｉｇｈｔｎｅｓｓ２がＰｌｉｇｈｔｎｅｓｓ１より大きい場合
Cc1=0               …（６７）
Cc2=Plightness      …（６８）
【０１０８】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された明るさに基づく補正量を用いて収差補正部４により補正処理を実施する。収差補正部１７による明るさに基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（６９）〜（７１）で表すことができる。
Pr1’=Pr-Cc1      …（６９）
Pg1’=Pg          …（７０）
Pb1’=Pb-Cc1      …（７１）
【０１０９】
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１、補正後の画素値の強度をＰｌｉｇｈｔｎｅｓｓ１’とし、次の式（７２）に従って補正処理を行うこともできる。
Plightness1’=Plightness1-Cc1      …（７２）
【０１１０】
同様に、補正対象の画素位置がＤ２方向の場合、補正前の画素値をＰｒ２，Ｐｇ２，Ｐｂ２、補正後の画素値をＰｒ２’，Ｐｇ２’，Ｐｂ２’とし、次の式（７３）〜（７５）に従って補正処理を行うことができる。
Pr2’=Pr-Cc2      …（７３）
Pg2’=Pg          …（７４）
Pb2’=Pb-Cc2      …（７５）
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の強度をＰｌｉｇｈｔｎｅｓｓ２’として、次の式（７６）で補正処理を行うこともできる。
Plightness2’=Plightness2-Cc2      …（７６）
【０１１１】
明るさに基づく色空間情報に従って軸上色収差を補正する場合、画像処理などにより色空間上で特定の色の明るさの差が広がるような処理が行われる場合に特に有効である。
【０１１２】
（角度及び強度に関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、角度及び強度に関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した角度及び強度に関する色空間情報からそれらの差分を夫々算出する（ステップＳ７）。
【０１１３】
補正量算出部１７は、上記した角度に基づく補正量Ｃａ及び強度に基づく補正量Ｃｂを算出し、これら２つの補正量に基づいて角度及び強度に基づく補正量Ｃａｂを、例えば、以下の式（７７）に従って算出する。
Cab=Ca+Cb      …（７７）
【０１１４】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度及び強度に基づく補正量Ｃａｂを用いて収差補正部４により補正処理を実施する。収差補正部４の補正処理は、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’として以下の式（７８）〜（８０）に従って行われる。
Pr’=Pr-Cab    …（７８）
Pg’=Pg        …（７９）
Pb’=Pb-Cab    …（８０）
【０１１５】
なお、補正前の画素値の角度をＰａ、補正後の画素値の角度をＰａ’として、次の式（８１）に従って補正を行うこともできる。
Pa’=Pa-Ca      …（８１）
また、補正前の画素値の強度をＰｓ、補正後の画素値の強度をＰｓ’として、以下の式（８２）に従って補正を行うこともできる。
Ps’=Ps-Cb      …（８２）
【０１１６】
続いて、角度に関する色空間情報に基づいて補正を行う場合の他の例について説明する。
補正量算出部１７により、例えば式（１７）〜（２０）に従って補正量Ｃａ１，Ｃａ２及びＣｂ１，Ｃｂ２を算出し、これらを用いて角度及び強度に基づく補正量Ｃａｂ１，Ｃａｂ２を、例えば、以下の式（８３）及び（８４）に従って算出することができる。算出方法は例えば次のようにする。
Cab1=Ca1+Cb1      …（８３）
Cab2=Ca1+Cb2      …（８４）
【０１１７】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度に基づく補正量を用いて収差補正部４により補正処理を実施する。
収差補正部１７による角度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（８５）〜（８７）で表すことができる。
Pr1’=Pr-Cab1      …（８５）
Pg1’=Pg           …（８６）
Pb1’=Pb-Cab1      …（８７）
【０１１８】
なお、補正前の画素値の角度をＰａｎｇｌｅ１、補正後の画素値の角度をＰａｎｇｌｅ１’とし、次の式（８８）に従って補正処理を行うこともできる。
Pangle1’=Pangle1-Cb1      …（８８）
また補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ１、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ１’とし、次の式（８９）に従って補正処理を行うこともできる。
Pstrength1’=Pstrength1-Cb1      …（８９）
【０１１９】
同様に、補正対象の画素位置がＤ２方向の場合、補正前の画素値をＰｒ２，Ｐｇ２，Ｐｂ２、補正後の画素値をＰｒ２’，Ｐｇ２’，Ｐｂ２’とし、次の式（９０）〜（９２）に従って補正処理を行うことができる。
Pr2’=Pr-Cab2      …（９０）
Pg2’=Pg           …（９１）
Pb2’=Pb-Cab2      …（９２）
【０１２０】
また、補正前の画素値の角度をＰａｎｇｌｅ２、補正後の画素値の角度をＰａｎｇｌｅ２’として、次の式（９３）で補正処理を行うこともできる。
Pangle2’=Pangle2-Cb2      …（９３）
さらに、補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ２、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ２’として、次の式（９４）で補正処理を行うこともできる。
Pstrength2’=Pstrength2-Cb2      …（９４）
【０１２１】
角度及び強度に基づく色空間情報に従って軸上色収差を補正する場合、、その他画像処理などにより色空間上で特定の色の角度の差及び強度の差が広がるような処理が行われる場合に特に有効である。
【０１２２】
（角度及び明るさに関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、角度及び明るさに関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した角度及び明るさに関する色空間情報からそれらの差分を夫々算出する（ステップＳ７）。
【０１２３】
補正量算出部１７で算出した角度に基づく補正量Ｃａ及び明るさに基づく補正量Ｃｃを用いて角度及び明るさに基づく補正量Ｃａｃを、例えば、以下の式（９５）に従って算出する。
Cac=Ca+Cc      …（９５）
【０１２４】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度及び強度に基づく補正量Ｃａｂを用いて収差補正部４により補正処理を実施する。収差補正部４の補正処理は、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’として以下の式（７８）〜（８０）に従って行われる。
Pr’=Pr-Cac      …（９６）
Pg’=Pg          …（９７）
Pb’=Pb-Cac      …（９８）
【０１２５】
なお、補正前の画素値の角度をＰａ、補正後の画素値の角度をＰａ’として、以下の式（９９）に従って補正を行うこともできる。
Pa’=Pa-Ca      …（９９）
また、補正前の画素値の明度をＰｌ、補正後の画素値の強度をＰｌ’として、以下の式（１００）に従って補正を行うこともできる。
Pl’=Pl-Cc      …（１００）
【０１２６】
続いて、角度及び明るさに関する色空間情報に基づいて補正を行う場合の他の例について説明する。補正量算出部１７により、角度に基づく補正量Ｃａ１、Ｃａ２及び明るさに基づく補正量Ｃｃ１，Ｃｃ２を算出し、これらを用いて角度及び明るさに基づく補正量Ｃａｃ１，Ｃａｃ２を、例えば、以下の式に従って算出することができる。
Cac1=Ca1+Cc1      …（１０１）
Cac2=Ca2+Cc2      …（１０２）
【０１２７】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度に基づく補正量を用いて収差補正部４により補正処理を実施する。収差補正部１７による角度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（１０３）〜（１０５）で表すことができる。
Pr1’=Pr-Cac1       …（１０３）
Pg1’=Pg            …（１０４）
Pb1’=Pb- Cac1      …（１０５）
【０１２８】
なお、補正前の画素値の角度をＰａｎｇｌｅ１、補正後の画素値の角度をＰａｎｇｌｅ１’とし、次の式（１０６）に従って補正処理を行うこともできる。
Pangle1’=Pangle1-Cb1      …（１０６）
【０１２９】
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１、補正後の画素値の強度をＰｌｉｇｈｔｎｅｓｓ１’とすると次式で補正を行っても良い。
Plightness1’=Plightness1-Cc1      …（１０７）
【０１３０】
同様に、補正対象の画素位置がＤ２方向の場合、補正前の画素値をＰｒ２，Ｐｇ２，Ｐｂ２、補正後の画素値をＰｒ２’，Ｐｇ２’，Ｐｂ２’とし、次の式（１０８）〜（１１０）に従って補正処理を行うことができる。
Pr2’=Pr-Cac2      …（１０８）
Pg2’=Pg           …（１０９）
Pb2’=Pb-Cac2      …（１１０）
【０１３１】
また、補正前の画素値の角度をＰａｎｇｌｅ２、補正後の画素値の角度をＰａｎｇｌｅ２’とし、次の式（１１１）に従って補正処理を行うことができる。
Pangle2’=Pangle2-Cb2      …（１１１）
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の強度をＰｌｉｇｈｔｎｅｓｓ２’とし、次の式（１１２）に従って補正処理を行うことができる。
Plightness2’=Plightness2-Cc2      …（１１２）
【０１３２】
角度及び明るさに基づく色空間情報に従って軸上色収差を補正する場合、画像処理などにより色空間上で特定の色の角度の差及び明るさの差が広がるような処理が行われる場合に特に有効である。
【０１３３】
（強度及び明るさに関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、強度及び明るさに関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した強度及び明るさに関する色空間情報からそれらの差分を夫々算出する（ステップＳ７）。
【０１３４】
補正量算出部１７で算出した強度に基づく補正量Ｃｂ及び明るさに基づく補正量Ｃｃを用いて強度及び明るさに基づく補正量Ｃｂｃを、例えば、以下の式（９５）に従って算出する。
Cbc=Cb+Cc      …（１１３）
【０１３５】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度及び強度に基づく補正量Ｃａｂを用いて収差補正部４により補正処理を実施する。
収差補正部４の補正処理は、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’として以下の式（１１４）〜（１１６）に従って行われる。
Pr’=Pr-Cbc      …（１１４）
Pg’=Pg          …（１１５）
Pb’=Pb-Cbc      …（１１６）
【０１３６】
なお、補正前の画素値の強度をＰｓ、補正後の画素値の強度をＰｓ’として、以下の式（１１７）に従って補正を行うこともできる。
Ps’=Ps-Cb      …（１１７）
また、補正前の画素値の明度をＰｌ、補正後の画素値の強度をＰｌ’として、以下の式（１１８）に従って補正を行うこともできる。
Pl’=Pl-Cc      …（１１８）
【０１３７】
続いて、強度及び明るさに関する色空間情報に基づいて補正を行う場合の他の例について説明する。補正量算出部１７により、強度に基づく補正量Ｃｂ１，Ｃｂ２及び明るさに基づく補正量Ｃｃ１，Ｃｃ２を算出し、これらを用いて強度及び明るさに基づく補正量Ｃｂｃ１，Ｃｂｃ２を例えば、以下の式（１１９）及び（１２０）に従って算出する。
Cbc1=Cb1+Cc1      …（１１９）
Cbc2=Cb2+Cc2      …（１２０）
【０１３８】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度に基づく補正量を用いて収差補正部４により補正処理を実施する。収差補正部１７による角度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（１２１）〜（１２３）で表すことができる。
Pr1’=Pr-Cbc1       …（１２１）
Pg1’=Pg            …（１２２）
Pb1’=Pb-Cbc1       …（１２３）
【０１３９】
なお、補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ１、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ１’とし、次の式（１２４）に従って補正処理を行うこともできる。
Pstrength1’=Pstrength1-Cb1      …（１２４）
【０１４０】
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１、補正後の画素値の強度をＰｌｉｇｈｔｎｅｓｓ１’とし、次の式（１２５）に従って補正処理を行うこともできる。
Plightness1’=Plightness1-Cc1      …（１２５）
【０１４１】
同様に、補正対象の画素位置がＤ２方向の場合、補正前の画素値をＰｒ２，Ｐｇ２，Ｐｂ２、補正後の画素値をＰｒ２’，Ｐｇ２’，Ｐｂ２’とし、次の式（１２６）〜（１２８）に従って補正処理を行うことができる。
Pr2’=Pr-Cbc2      …（１２６）
Pg2’=Pg           …（１２７）
Pb2’=Pb-Cbc2      …（１２８）
【０１４２】
なお、補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ２、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ２’とし、次の式（１２９）に従って補正処理を行うことができる。
Pstrength2’=Pstrength2-Cb2      …（１２９）
さらに、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の強度をＰｌｉｇｈｔｎｅｓｓ２’とし、次の式（１３０）に従って補正処理を行うことができる。
Plightness2’=Plightness2-Cc2      …（１３０）
【０１４３】
強度及び明るさに基づく色空間情報に従って軸上色収差を補正する場合、画像処理等により色空間上で特定の色の強度の差及び明るさの差が広がるような処理が行われる場合に特に有効である。
【０１４４】
（強度及び明るさに関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、角度、強度及び明るさに関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した角度、強度及び明るさに関する色空間情報からそれらの差分を夫々算出する（ステップＳ７）。
【０１４５】
補正量算出部で算出した角度に基づく補正量Ｃａ、強度に基づく補正量Ｃｂ、明るさに基づく補正量Ｃｃを用いて角度、強度及び明るさ補正量Ｃａｂｃ、例えば、以下の式（１３１）に従って算出する。
Cabc=Ca+Cb+Cc      …（１３１）
【０１４６】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度及び強度に基づく補正量Ｃａｂを用いて収差補正部４により補正処理を実施する。収差補正部４の補正処理は、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’として以下の式（１３２）〜（１３４）に従って行われる。
Pr’=Pr-Cabc        …（１３２）
Pg’=Pg             …（１３３）
Pb’=Pb-Cabc        …（１３４）
【０１４７】
なお、補正前の画素値の角度をＰａ、補正後の画素値の強度をＰａ’として、以下の式（１３５）に従って補正を行うこともできる。
Pa’=Pa-Ca      …（１３５）
また、補正前の画素値の強度をＰｓ、補正後の画素値の強度をＰｓ’として、以下の式（１３６）に従って補正を行うこともできる。
Ps’=Ps-Cb      …（１３６）
さらに、補正前の画素値の明度をＰｌ、補正後の画素値の強度をＰｌ’として、以下の式（１３７）に従って補正を行うこともできる。
Pl’=Pl-Cc      …（１３７）
【０１４８】
続いて、角度、強度及び明るさに関する色空間情報に基づいて補正を行う場合の他の例について説明する。  補正量算出部１７で算出した明るさに基づく補正量Ｃａ１，Ｃａ２、強度に基づく補正量Ｃｂ１，Ｃｂ２、明るさに基づく補正量Ｃｃ１，Ｃｃ２を用いて角度、強度及び明るさに基づく補正量Ｃａｂｃ１，Ｃａｂｃ２を、例えば、以下の式（１３８）及び（１３９）に従って算出する。
Cabc1=Ca1+Cb1+Cc1      …（１３８）
Cabc2=Ca2+Cb2+Cc2      …（１３９）
【０１４９】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度に基づく補正量を用いて収差補正部４により補正処理を実施する。収差補正部１７による角度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（１４０）〜（１４２）で表すことができる。
Pr1’=Pr-Cabc1      …（１４０）
Pg1’=Pg            …（１４１）
Pb1’=Pb-Cabc1      …（１４２）
【０１５０】
なお、補正前の画素値の角度をＰａｎｇｌｅ１、補正後の画素値の強度をＰａｎｇｌｅ１’とし、次の式（１４３）に従って補正処理を行うこともできる。
Pangle1’=Pangle1-Ca1      …（１４３）
【０１５１】
また、補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ１、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ１’とし、次の式（１４４）に従って補正処理を行うこともできる。
Pstrength1’=Pstrength1-Cb1      …（１４４）
【０１５２】
さらに、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１、補正後の画素値の強度をＰｌｉｇｈｔｎｅｓｓ１’とし、次の式（１４５）に従って補正処理を行うこともできる。
Plightness1’=Plightness1-Cc1      …（１４５）
【０１５３】
同様に、補正対象の画素位置がＤ２方向の場合、補正前の画素値をＰｒ２，Ｐｇ２，Ｐｂ２、補正後の画素値をＰｒ２’，Ｐｇ２’，Ｐｂ２’と
し、次の式（１４６）〜（１４８）に従って補正処理を行うことができる。
Pr2’=Pr-Cabc2        …（１４６）
Pg2’=Pg              …（１４７）
Pb2’=Pb-Cabc2        …（１４８）
【０１５４】
なお、補正前の画素値の角度をＰａｎｇｌｅ２、補正後の画素値の強度をＰａａｎｇｌｅ２’とし、次の式（１４９）に従って補正処理を行うことができる。
。
【０１５５】
Paangle2’=Paangle2-Ca2      …（１４９）
また補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ２、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ２’とし、次の式（１５０）に従って補正処理を行うことができる。
Pstrength2’=Pstrength2-Cb2      …（１５０）
【０１５６】
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の強度をＰｌｉｇｈｔｎｅｓｓ２’とし、次の式（１５１）に従って補正処理を行うことができる。
Plightness2’=Plightness2-Cc2      …（１５１）
【０１５７】
角度、強度及び明るさに基づく色空間情報に従って軸上色収差を補正する場合、画像処理などにより色空間上で特定の色の角度の差、強度の差及び明るさの差が広がるような処理が行われる場合に特に有効である。
【０１５８】
このように、本実施形態に係る画像処理装置及び画像処理装置によって実施される画像処理方法によれば、近飽和領域もしくは近飽和領域以外の領域少なくとも一方においても軸上色収差を検出して人間の視覚特性に基づいて補正することができる。
【０１５９】
なお、本実施形態に係る画像処理装置１においては、上述した条件式（２）〜（５）全てを満たす場合に軸上色収差の発生起点であるとして検出することとしたが、これに代えて、図１０及び図１１に示されるように、条件式（２）及び条件式（３）を満たす場合、または、条件式（４）及び条件式（５）を満たす場合に、注目画素を軸上色収差の発生起点であるとして検出することにしてもよい。
このようにすることで、フォーカス位置より広い奥行き範囲で軸上色収差の発生起点を検出することができるという利点がある。
【０１６０】
また、条件式（２）、（３）、（５）を全て満たす場合、または、条件式（３）、（４）、（５）を全て満たす場合に、注目画素を軸上色収差の発生起点として検出してもよい。
このようにすることで、他方向Ｄ２にパープルフリンジが発生し、一方向Ｄ１のＲの画素値の合計またはＢの画素値の合計が、Ｇの画素値の合計より小さく、色相環で反対色に近い色となってパープルフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。
【０１６１】
また、条件式（２）、（３）、（４）を全て満たす場合、または、条件式（２）、（４）、（５）を全て満たす場合に、注目画素を軸上色収差の発生起点として検出してもよい。
このようにすることで、一方向Ｄ１にグリーンフリンジが発生し、他方向Ｄ２のＲの画素値の合計またはＢの画素値の合計が、Ｇの画素値の合計より大きく、色相環で反対色に近い色となってグリーンフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。
【０１６２】
また、数１で算出された画素値の合計が、次の条件式（１５２）〜（１５５）を全て満たす場合に注目画素を軸上色収差の発生起点として検出することにしてもよい。
ＰｇＤ１＋ＴＨＲ１＞ＰｒＤ１                …（１５２）
ＰｇＤ２−ＴＨＲ２＜ＰｒＤ２                …（１５３）
ＰｇＤ１＋ＴＨＢ１＞ＰｂＤ１                …（１５４）
ＰｇＤ２−ＴＨＢ２＜ＰｂＤ２                …（１５５）
ここで、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２は任意の定数である。
【０１６３】
ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２としては、図１２に示すように、例えば、処理対象の画素位置におけるＧとＲ、または、ＧとＢのノイズ量を示す標準偏差の差を例示することができる。一方向Ｄ１及び他方向Ｄ２のＲ，Ｇ，Ｂ夫々のノイズ量をＮｒＤ１、ＮｇＤ１、ＮｂＤ１、ＮｒＤ２、ＮｇＤ２、ＮｂＤ２として、以下の式（１５６）〜式（１５９）で示すことができる。
ＴＨＲ１＝ＮｒＤ１−ＮｇＤ１      …（１５６）
ＴＨＲ２＝ＮｒＤ２−ＮｇＤ２      …（１５７）
ＴＨＢ１＝ＮｂＤ１−ＮｇＤ１      …（１５８）
ＴＨＢ２＝ＮｂＤ２−ＮｇＤ２      …（１５９）
このようにすることで、ノイズの影響を低減し、安定して軸上色収差の発生起点を検出することができるという利点がある。
【０１６４】
また、数１で算出された画素値の合計が、次の条件式（１６０）〜（１６３）を全て満たす場合に注目画素を軸上色収差の発生起点として検出することにしてもよい。
ＰｇＤ１＋ＴＨＲ＞ＰｒＤ１                …（１６０）
ＰｇＤ２＋ＴＨＲ＜ＰｒＤ２                …（１６１）
ＰｇＤ１＋ＴＨＢ＞ＰｂＤ１                …（１６２）
ＰｇＤ２＋ＴＨＢ＜ＰｂＤ２                …（１６３）
ここで、ＴＨＲ，ＴＨＢは任意の定数である。
【０１６５】
定数ＴＨＲ，ＴＨＢとしては、図１３に示されるように、例えば、処理対象の画素位置におけるＧとＲ、または、ＧとＢの画素値の差分を例示することができる。注目画素のＲ，Ｇ，Ｂの夫々の画素値をＰｒ（０，０）、Ｐｇ（０，０）、Ｐｂ（０，０）とすると、以下の式（１６４）及び式（１６５）となる。
ＴＨＲ＝Ｐｇ（０，０）−−Ｐｒ（０，０）      …（１６４）
ＴＨＢ＝Ｐｇ（０，０）−−Ｐｂ（０，０）      …（１６５）
このようにすることで、無彩色ではなく、色つきの被写体上に発生した軸上色収差の発生起点を検出することができるという利点がある。
【０１６６】
また、数１で算出された画素値の合計が、次の条件式（１６６）〜（１６９）を全て満たす場合に注目画素を軸上色収差の発生起点として検出することにしてもよい。
ＰｇＤ１＞ＰｒＤ１×ＫＲ                …（１６６）
ＰｇＤ２＜ＰｒＤ２×ＫＲ                …（１６７）
ＰｇＤ１＞ＰｂＤ１×ＫＢ                …（１６８）
ＰｇＤ２＜ＰｂＤ２×ＫＢ                …（１６９）
ここで、ＫＲ，ＫＢは任意の定数である。
【０１６７】
定数ＫＲ，ＫＢとしては、図１４に示されるように、例えば、処理対象の画素位置におけるＧとＲ、または、ＧとＢの画素値の比率を例示することができる。注目画素のＲ，Ｇ，Ｂの夫々の画素値をＰｒ（０，０）、Ｐｇ（０，０）、Ｐｂ（０，０）とすると、以下の式（１７０）及び式（１７１）となる。
ＫＲ＝Ｐｇ（０，０）／Ｐｒ（０，０）      …（１７０）
ＫＢ＝Ｐｇ（０，０）／Ｐｂ（０，０）      …（１７１）
このようにすることによっても、無彩色ではなく、色つきの被写体上に発生した軸上色収差の発生起点を検出することができるという利点がある。
【０１６８】
また、本実施形態においては、各画素においてＲＧＢ全ての信号値が得られるフルカラー画像が入力画像として入力される場合について例示したが、これに代えて、いわゆるベイヤ画像である場合に適用することにしてもよい。この場合、一方向Ｄ１の合計画素値算出部４における合計画素値の算出方法は数２〜数４の通りとなる。
【０１６９】
図１５に示すように、注目画素がＲ（斜線部分）である場合には、数２及び数３に示す数式により合計画素値が算出される。
【数２】
2016181476.tif 000004
【数３】
2016181476.tif 000005
【０１７０】
図１６に示されるように、注目画素がＢ（斜線部分）である場合には、数４及び数５に示される数式により合計画素値が算出される。
【数４】
2016181476.tif 000006
【数５】
2016181476.tif 000007
【０１７１】
図１７に示されるように、注目画素がＧＲ（斜線部分）である場合には、数６及び数７に示される数式により合計画素値が算出される。
【数６】
2016181476.tif 000008
【数７】
2016181476.tif 000009
【０１７２】
図１８に示されるように、注目画素がＧＢ（斜線部分）である場合には、数８及び数９に示される数式により合計画素値が算出される。
【数８】
2016181476.tif 000010
【数９】
2016181476.tif 000011
【０１７３】
このようにすることで、同時化されていない画像に対しても軸上色収差の発生起点を検出することができ、検出時に同時化のための補間エラーの影響をなくすことができる。
本実施形態に係る画像処理装置１は、放送用据え置き型カメラ、ＥＮＧカメラ、民生用ハンディカメラ、デジタルカメラ等の製品に搭載して使用される。また、動画を扱う画像信号補正プログラム（ＣＧプログラム）や画像編集装置等にも用いることができる。また、本実施形態に係る画像処理方法をコンピュータ（例えば、パーソナルコンピュータ）に実行させる画像処理プログラムとして用いることもできる。
【符号の説明】
【０１７４】
２  領域抽出部
３  軸上色収差補正量算出部
４  収差補正部
１０  合計画素値算出部
１１  発生起点検出部
１２  領域判定部
１５  色空間情報算出部
１６  色空間差分算出部
１７  補正量算出部
【図１】
2016181476.tif 000012
【図２】
2016181476.tif 000013
【図３】
2016181476.tif 000014
【図４】
2016181476.tif 000015
【図５】
2016181476.tif 000016
【図６】
2016181476.tif 000017
【図７】
2016181476.tif 000018
【図８】
2016181476.tif 000019
【図９】
2016181476.tif 000020
【図１０】
2016181476.tif 000021
【図１１】
2016181476.tif 000022
【図１２】
2016181476.tif 000023
【図１３】
2016181476.tif 000024
【図１４】
2016181476.tif 000025
【図１５】
2016181476.tif 000026
【図１６】
2016181476.tif 000027
【図１７】
2016181476.tif 000028
【図１８】
2016181476.tif 000029
【手続補正書】
【提出日】平成29年11月8日(2017.11.8)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出部と、
該合計画素値算出部による算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素が軸上色収差の発生起点であるか否かを判定し、発生起点を検出する発生起点検出部と、
前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定部と、
前記発生起点として検出された前記注目画素に対し、前記合計画素値算出部による算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出部と、
該色空間情報算出部により算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出部と、
該色空間差分算出部により算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出部と、
前記軸上色収差領域について前記補正量を用いて補正する補正部と、
を備える画像処理装置。
【請求項２】
前記色空間情報算出部が、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
前記色空間差分算出部が、前記色空間情報算出部による算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
前記補正量算出部が、前記色空間差分算出部による算出結果に基づいて軸上色収差補正の補正量を算出する請求項１記載の画像処理装置。
【請求項３】
前記発生起点検出部が、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項１または２記載の画像処理装置。
【請求項４】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ２＜ＰｒＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計である。
【請求項５】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｂＤ２
ここで、
ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
【請求項６】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ２＜ＰｒＤ２
ＰｇＤ２＜ＰｂＤ２
ここで、
ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
【請求項７】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｒＤ２
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
【請求項８】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｒＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、である。
【請求項９】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計
である。
【請求項１０】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＋ＴＨＲ１＞ＰｒＤ１
ＰｇＤ２−−ＴＨＲ２＜ＰｒＤ２
ＰｇＤ１＋ＴＨＢ１＞ＰｂＤ１
ＰｇＤ２−−ＴＨＢ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２は任意の定数である。
【請求項１１】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＋ＴＨＲ＞ＰｒＤ１
ＰｇＤ２＋ＴＨＲ＜ＰｒＤ２
ＰｇＤ１＋ＴＨＢ＞ＰｂＤ１
ＰｇＤ２＋ＴＨＢ＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ，ＴＨＢは任意の定数である。
【請求項１２】
前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
ＰｇＤ１＞ＰｒＤ１×ＫＲ
ＰｇＤ２＜ＰｒＤ２×ＫＲ
ＰｇＤ１＞ＰｂＤ１×ＫＢ
ＰｇＤ２＜ＰｂＤ２×ＫＢ
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＫＲ，ＫＢは任意の定数である。
【請求項１３】
前記画像が、各画素にＲＧＢ全ての画素値を有する請求項１または２記載の画像処理装置。
【請求項１４】
前記画像が、各画素においてＲＧＢの少なくとも１つの画素値が欠落した画素値を有する請求項１または２記載の画像処理装置。
【請求項１５】
入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、
該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、
前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、
前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、
該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、
該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、
前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、
を備える画像処理方法。
【請求項１６】
前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出する請求項１５記載の画像処理方法。
【請求項１７】
前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出ステップにおいて算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項１５又は１６記載の画像処理方法。
【請求項１８】
入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、
該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、
前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、
前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、
該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、
該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、
前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、
をコンピュータに実行させる画像処理プログラム。
【請求項１９】
前記コンピュータが、
前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出する請求項１８記載の画像処理プログラム。
【請求項２０】
前記コンピュータが、
前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出ステップにおいて算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項１７又は１８記載の画像処理プログラム。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００３３
【補正方法】変更
【補正の内容】
【００３３】
また、上記他の態様において、前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出ステップにおいて算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出することとしてもよい。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００３６
【補正方法】変更
【補正の内容】
【００３６】
また、上記他の態様において、前記コンピュータが、前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出ステップにおいて算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出することとしてもよい。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００９８
【補正方法】変更
【補正の内容】
【００９８】
また、特定色空間がＬ＊ａ＊ｂ＊空間である場合、明るさに関する色空間情報として、Ｐｌ＊Ｄ１，Ｐｌ＊Ｄ２を算出し、これらの差分の絶対値Ｐｌｉｇｈｔｎｅｓｓを次の式（５４）に従って算出する。
Plightness=|Pl*D1-Pl*D2|      …（５４）
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】０１０２
【補正方法】変更
【補正の内容】
【０１０２】
また、補正前の画素値の明度をＰｌ、補正後の画素値の明度をＰｌ’とし、次の式（５９）に従って補正処理を行うこともできる。
Pl’=Pl-Cc      …（５９）
【手続補正６】
【補正対象書類名】明細書
【補正対象項目名】０１１０
【補正方法】変更
【補正の内容】
【０１１０】
同様に、補正対象の画素位置がＤ２方向の場合、補正前の画素値をＰｒ２，Ｐｇ２，Ｐｂ２、補正後の画素値をＰｒ２’，Ｐｇ２’，Ｐｂ２’とし、次の式（７３）〜（７５）に従って補正処理を行うことができる。
Pr2’=Pr-Cc2      …（７３）
Pg2’=Pg          …（７４）
Pb2’=Pb-Cc2      …（７５）
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２’として、次の式（７６）で補正処理を行うこともできる。
Plightness2’=Plightness2-Cc2      …（７６）
【手続補正７】
【補正対象書類名】明細書
【補正対象項目名】０１１６
【補正方法】変更
【補正の内容】
【０１１６】
続いて、角度及び強度に関する色空間情報に基づいて補正を行う場合の他の例について説明する。
補正量算出部１７により、例えば式（１７）〜（２０）に従って補正量Ｃａ１，Ｃａ２及びＣｂ１，Ｃｂ２を算出し、これらを用いて角度及び強度に基づく補正量Ｃａｂ１，Ｃａｂ２を、例えば、以下の式（８３）及び（８４）に従って算出することができる。算出方法は例えば次のようにする。
Cab1=Ca1+Cb1      …（８３）
Cab2=Ca2+Cb2      …（８４）
【手続補正８】
【補正対象書類名】明細書
【補正対象項目名】０１２４
【補正方法】変更
【補正の内容】
【０１２４】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度及び明度に基づく補正量Ｃａｃを用いて収差補正部４により補正処理を実施する。収差補正部４の補正処理は、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’として以下の式（９６）〜（９８）に従って行われる。
Pr’=Pr-Cac      …（９６）
Pg’=Pg          …（９７）
Pb’=Pb-Cac      …（９８）
【手続補正９】
【補正対象書類名】明細書
【補正対象項目名】０１２５
【補正方法】変更
【補正の内容】
【０１２５】
なお、補正前の画素値の角度をＰａ、補正後の画素値の角度をＰａ’として、以下の式（９９）に従って補正を行うこともできる。
Pa’=Pa-Ca      …（９９）
また、補正前の画素値の明度をＰｌ、補正後の画素値の明度をＰｌ’として、以下の式（１００）に従って補正を行うこともできる。
Pl’=Pl-Cc      …（１００）
【手続補正１０】
【補正対象書類名】明細書
【補正対象項目名】０１２７
【補正方法】変更
【補正の内容】
【０１２７】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度及び明度に基づく補正量を用いて収差補正部４により補正処理を実施する。収差補正部１７による角度及び明度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（１０３）〜（１０５）で表すことができる。
Pr1’=Pr-Cac1       …（１０３）
Pg1’=Pg            …（１０４）
Pb1’=Pb- Cac1      …（１０５）
【手続補正１１】
【補正対象書類名】明細書
【補正対象項目名】０１２９
【補正方法】変更
【補正の内容】
【０１２９】
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１、補正後の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１’とすると次式で補正を行っても良い。
Plightness1’=Plightness1-Cc1      …（１０７）
【手続補正１２】
【補正対象書類名】明細書
【補正対象項目名】０１３１
【補正方法】変更
【補正の内容】
【０１３１】
また、補正前の画素値の角度をＰａｎｇｌｅ２、補正後の画素値の角度をＰａｎｇｌｅ２’とし、次の式（１１１）に従って補正処理を行うことができる。
Pangle2’=Pangle2-Cb2      …（１１１）
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２’とし、次の式（１１２）に従って補正処理を行うことができる。
Plightness2’=Plightness2-Cc2      …（１１２）
【手続補正１３】
【補正対象書類名】明細書
【補正対象項目名】０１３４
【補正方法】変更
【補正の内容】
【０１３４】
補正量算出部１７で算出した強度に基づく補正量Ｃｂ及び明るさに基づく補正量Ｃｃを用いて強度及び明るさに基づく補正量Ｃｂｃを、例えば、以下の式（１１３）に従って算出する。
Cbc=Cb+Cc      …（１１３）
【手続補正１４】
【補正対象書類名】明細書
【補正対象項目名】０１３５
【補正方法】変更
【補正の内容】
【０１３５】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された強度及び明度に基づく補正量Ｃａｃを用いて収差補正部４により補正処理を実施する。
収差補正部４の補正処理は、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’として以下の式（１１４）〜（１１６）に従って行われる。
Pr’=Pr-Cbc      …（１１４）
Pg’=Pg          …（１１５）
Pb’=Pb-Cbc      …（１１６）
【手続補正１５】
【補正対象書類名】明細書
【補正対象項目名】０１３６
【補正方法】変更
【補正の内容】
【０１３６】
なお、補正前の画素値の強度をＰｓ、補正後の画素値の強度をＰｓ’として、以下の式（１１７）に従って補正を行うこともできる。
Ps’=Ps-Cb      …（１１７）
また、補正前の画素値の明度をＰｌ、補正後の画素値の明度をＰｌ’として、以下の式（１１８）に従って補正を行うこともできる。
Pl’=Pl-Cc      …（１１８）
【手続補正１６】
【補正対象書類名】明細書
【補正対象項目名】０１３８
【補正方法】変更
【補正の内容】
【０１３８】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された強度及び明度に基づく補正量を用いて収差補正部４により補正処理を実施する。収差補正部１７による強度及び明度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（１２１）〜（１２３）で表すことができる。
Pr1’=Pr-Cbc1       …（１２１）
Pg1’=Pg            …（１２２）
Pb1’=Pb-Cbc1       …（１２３）
【手続補正１７】
【補正対象書類名】明細書
【補正対象項目名】０１４０
【補正方法】変更
【補正の内容】
【０１４０】
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１、補正後の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１’とし、次の式（１２５）に従って補正処理を行うこともできる。
Plightness1’=Plightness1-Cc1      …（１２５）
【手続補正１８】
【補正対象書類名】明細書
【補正対象項目名】０１４２
【補正方法】変更
【補正の内容】
【０１４２】
なお、補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ２、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ２’とし、次の式（１２９）に従って補正処理を行うことができる。
Pstrength2’=Pstrength2-Cb2      …（１２９）
さらに、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２’とし、次の式（１３０）に従って補正処理を行うことができる。
Plightness2’=Plightness2-Cc2      …（１３０）
【手続補正１９】
【補正対象書類名】明細書
【補正対象項目名】０１４４
【補正方法】変更
【補正の内容】
【０１４４】
（角度、強度及び明るさに関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、角度、強度及び明るさに関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した角度、強度及び明るさに関する色空間情報からそれらの差分を夫々算出する（ステップＳ７）。
【手続補正２０】
【補正対象書類名】明細書
【補正対象項目名】０１４５
【補正方法】変更
【補正の内容】
【０１４５】
補正量算出部で算出した角度に基づく補正量Ｃａ、強度に基づく補正量Ｃｂ、明るさに基づく補正量Ｃｃを用いて角度、強度及び明るさに基づく補正量Ｃａｂｃ、例えば、以下の式（１３１）に従って算出する。
Cabc=Ca+Cb+Cc      …（１３１）
【手続補正２１】
【補正対象書類名】明細書
【補正対象項目名】０１４６
【補正方法】変更
【補正の内容】
【０１４６】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度、強度及び明度に基づく補正量Ｃａｂｃを用いて収差補正部４により補正処理を実施する。収差補正部４の補正処理は、図６に示すように補正前の画素値をＰｒ，Ｐｇ，Ｐｂ、図７に示すように補正後の画素値をＰｒ’，Ｐｇ’，Ｐｂ’として以下の式（１３２）〜（１３４）に従って行われる。
Pr’=Pr-Cabc        …（１３２）
Pg’=Pg             …（１３３）
Pb’=Pb-Cabc        …（１３４）
【手続補正２２】
【補正対象書類名】明細書
【補正対象項目名】０１４７
【補正方法】変更
【補正の内容】
【０１４７】
なお、補正前の画素値の角度をＰａ、補正後の画素値の角度をＰａ’として、以下の式（１３５）に従って補正を行うこともできる。
Pa’=Pa-Ca      …（１３５）
また、補正前の画素値の強度をＰｓ、補正後の画素値の強度をＰｓ’として、以下の式（１３６）に従って補正を行うこともできる。
Ps’=Ps-Cb      …（１３６）
さらに、補正前の画素値の明度をＰｌ、補正後の画素値の明度をＰｌ’として、以下の式（１３７）に従って補正を行うこともできる。
Pl’=Pl-Cc      …（１３７）
【手続補正２３】
【補正対象書類名】明細書
【補正対象項目名】０１４８
【補正方法】変更
【補正の内容】
【０１４８】
続いて、角度、強度及び明るさに関する色空間情報に基づいて補正を行う場合の他の例について説明する。  補正量算出部１７で算出した角度に基づく補正量Ｃａ１，Ｃａ２、強度に基づく補正量Ｃｂ１，Ｃｂ２、明るさに基づく補正量Ｃｃ１，Ｃｃ２を用いて角度、強度及び明るさに基づく補正量Ｃａｂｃ１，Ｃａｂｃ２を、例えば、以下の式（１３８）及び（１３９）に従って算出する。
Cabc1=Ca1+Cb1+Cc1      …（１３８）
Cabc2=Ca2+Cb2+Cc2      …（１３９）
【手続補正２４】
【補正対象書類名】明細書
【補正対象項目名】０１４９
【補正方法】変更
【補正の内容】
【０１４９】
入力画像中の全ての画素について処理が完了したら、軸上色収差領域について、補正量算出部１７により算出された角度、強度及び明度に基づく補正量を用いて収差補正部４により補正処理を実施する。収差補正部１７による角度に基づく補正処理は、例えば、図８に示すように補正対象の画素位置がＤ１方向の場合、補正前の画素値をＰｒ１，Ｐｇ１，Ｐｂ１（図８）、補正後の画素値をＰｒ１’，Ｐｇ１’，Ｐｂ１’（図９）として、以下の式（１４０）〜（１４２）で表すことができる。
Pr1’=Pr-Cabc1      …（１４０）
Pg1’=Pg            …（１４１）
Pb1’=Pb-Cabc1      …（１４２）
【手続補正２５】
【補正対象書類名】明細書
【補正対象項目名】０１５０
【補正方法】変更
【補正の内容】
【０１５０】
なお、補正前の画素値の角度をＰａｎｇｌｅ１、補正後の画素値の角度をＰａｎｇｌｅ１’とし、次の式（１４３）に従って補正処理を行うこともできる。
Pangle1’=Pangle1-Ca1      …（１４３）
【手続補正２６】
【補正対象書類名】明細書
【補正対象項目名】０１５２
【補正方法】変更
【補正の内容】
【０１５２】
さらに、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１、補正後の画素値の明度をＰｌｉｇｈｔｎｅｓｓ１’とし、次の式（１４５）に従って補正処理を行うこともできる。
Plightness1’=Plightness1-Cc1      …（１４５）
【手続補正２７】
【補正対象書類名】明細書
【補正対象項目名】０１５４
【補正方法】変更
【補正の内容】
【０１５４】
なお、補正前の画素値の角度をＰａｎｇｌｅ２、補正後の画素値の角度をＰａａｎｇｌｅ２’とし、次の式（１４９）に従って補正処理を行うことができる。
Pangle2’=Pangle2-Ca2      …（１４９）
【手続補正２８】
【補正対象書類名】明細書
【補正対象項目名】０１５５
【補正方法】変更
【補正の内容】
【０１５５】
また補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ２、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ２’とし、次の式（１５０）に従って補正処理を行うことができる。
Pstrength2’=Pstrength2-Cb2      …（１５０）
【手続補正２９】
【補正対象書類名】明細書
【補正対象項目名】０１５６
【補正方法】変更
【補正の内容】
【０１５６】
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２’とし、次の式（１５１）に従って補正処理を行うことができる。
Plightness2’=Plightness2-Cc2      …（１５１）
【国際調査報告】
2016181476.tif 000030
2016181476.tif 000031