JP6633624B2 - 画像処理装置、画像処理方法及び画像処理プログラム - Google Patents

画像処理装置、画像処理方法及び画像処理プログラム Download PDF

Info

Publication number
JP6633624B2
JP6633624B2 JP2017517499A JP2017517499A JP6633624B2 JP 6633624 B2 JP6633624 B2 JP 6633624B2 JP 2017517499 A JP2017517499 A JP 2017517499A JP 2017517499 A JP2017517499 A JP 2017517499A JP 6633624 B2 JP6633624 B2 JP 6633624B2
Authority
JP
Japan
Prior art keywords
pixel
sum
pixel values
pixels arranged
color space
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
JP2017517499A
Other languages
English (en)
Other versions
JPWO2016181476A1 (ja
Inventor
武史 福冨
武史 福冨
理子 朝倉
理子 朝倉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Olympus Corp
Original Assignee
Olympus Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Olympus Corp filed Critical Olympus Corp
Publication of JPWO2016181476A1 publication Critical patent/JPWO2016181476A1/ja
Application granted granted Critical
Publication of JP6633624B2 publication Critical patent/JP6633624B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00—Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/60—Noise processing, e.g. detecting, correcting, reducing or removing noise
    • H04N25/61—Noise processing, e.g. detecting, correcting, reducing or removing noise the noise originating only from the lens unit, e.g. flare, shading, vignetting or "cos4"
    • H04N25/611—Correction of chromatic aberration
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N25/00—Circuitry of solid-state image sensors [SSIS]; Control thereof
    • H04N25/10—Circuitry of solid-state image sensors [SSIS]; Control thereof for transforming different wavelengths into image signals
    • H04N25/11—Arrangement of colour filter arrays [CFA]; Filter mosaics
    • H04N25/13—Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements
    • H04N25/134—Arrangement of colour filter arrays [CFA]; Filter mosaics characterised by the spectral characteristics of the filter elements based on three different wavelength filter elements
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N9/00—Details of colour television systems
    • H04N9/64—Circuits for processing colour signals
    • H04N9/646—Circuits for processing colour signals for image enhancement, e.g. vertical detail restoration, cross-colour elimination, contour correction, chrominance trapping filters

Landscapes

  • Engineering & Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Image Processing (AREA)
  • Color Television Image Signal Generators (AREA)

Description

本発明は、画像処理装置、画像処理方法及び画像処理プログラムに関するものである。
通常の光学系を通して撮影された画像には軸上色収差が発生することが知られている。ここで、軸上色収差とは、光学系、撮像系及び画像処理系に起因し、光軸上の焦点位置が波長によって異なるために生じる偽色、すなわち本来の被写体には生じない色付きをいう。そして、このような軸上色収差を画像処理で補正するための技術が提案されている（例えば、特許文献１参照。）。
特許文献１では、画像データ中の注目画素について、偽色可能性を示すパラメータとして、白とびへの近さ、彩度の大きさ及び特定色相への近さを評価し、これらの少なくともいずれかを適用して軸上色収差を補正している。
特許第４５３９２７８号公報
しかしながら、軸上色収差を補正するに際し、特許文献１のように注目画素について評価し、これに基づいて補正を行うだけでは不十分である。すなわち、いくつかの注目画素が本来同一色であったとしても、当該注目画素に隣接する画素との関係によって、人間にはその視覚特性上それら注目画素が互いに異なる色であるように見える。このため、人間の視覚特性に鑑みると、注目画素のみを評価して軸上色収差を補正するだけでは不十分となる。
本発明は、上述した事情に鑑みてなされたものであって、人間の視覚特性に基づいて軸上色収差を精度よく補正することを目的としている。
上記目的を達成するため、本発明は以下の手段を提供する。
本発明の一態様は、入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出部と、該合計画素値算出部による算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素が軸上色収差の発生起点であるか否かを判定し、発生起点を検出する発生起点検出部と、前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定部と、前記発生起点として検出された前記注目画素に対し、前記合計画素値算出部による算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出部と、該色空間情報算出部により算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出部と、該色空間差分算出部により算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出部と、前記軸上色収差領域について前記補正量を用いて補正する補正部と、を備える画像処理装置である。
本態様によれば、注目画素を挟む両側に一直線上に配列されている画素の画素値の合計を合計画素値算出部により算出し、この算出結果又は注目画素の画素値に基づいて発生起点検出部により注目画素が軸上色収差の発生起点であるか判定し、発生起点を検出する。注目画素が軸上色収差の発生起点であるとして検出された場合、色空間情報算出部により、注目画素について、特定色空間上の色空間情報を一方向及び他方向について夫々算出し、色空間差分算出部によりこれらの差分を算出する。そして、この差分に基づいて算出された補正量を用いて、注目画素を中心とする周辺領域からなる軸上色収差領域について補正処理を行う。このように、注目画素だけでなく、注目画素を挟んで両側に一直線上に配列されている画素の画素値をも評価した補正量を算出することで、人間の視覚特性に基づいて軸上色収差補正を行うことができる。
上記態様においては、前記色空間情報算出部が、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、前記色空間差分算出部が、前記色空間情報算出部による算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、前記補正量算出部が、前記色空間差分算出部による算出結果に基づいて軸上色収差補正の補正量を算出することとしてもよい。
このようにすることで、人間がその視覚特性上、軸上色収差を知覚しやすい隣接する領域の色の角度の差、色の強度の差、明るさの差に応じた最適な補正量を算出することができる。
また上記様態においては、前記補正量算出部が、前記色の角度の差分が１８０度に近いほど大きな補正量を算出し、１８０度から離れている程小さな補正量を算出することとしてもよい。
このようにすることで、色空間の角度を色相環のように０度から３６０度で表した場合に、２つの色が対極となる１８０度のように、人間の視覚的に色付きが検知されやすい場合に適切な補正量を算出することができる。
また上記様態においては、前記補正量算出部が、前記色の強度の差分の絶対値が大きいほど大きな補正量を算出し、絶対値が小さいほど小さな補正量を算出することとしてもよい。
このようにすることで、色の強度の差分が大きく、人間の視覚的に色付きが検知されやすい場合に適切な補正量を算出することができる。
また上記様態においては、前記補正量算出部が、前記色の明るさの差分の絶対値が大きいほど大きな補正量を算出し、絶対値が小さいほど小さな補正量を算出することとしてもよい。
このようにすることで、色の明るさの差分が大きく、人間の視覚的に色付きが検知されやすい場合にも適切な補正量を算出することができる。
また上記様態においては、前記補正量算出部が、前記色の角度が前記一方向と前記他方向で同方向となる補正量を算出することとしてもよい。
このようにすることで、補正後に色の角度の差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。
また上記様態においては、前記補正量算出部が、前記色の強度が前記一方向と前記他方向とで同一となる補正量を算出することとしてもよい。
このようにすることで、補正後に色の強度の差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。
また上記様態においては、前記補正量算出部が、前記色の明るさが前記一方向と前記他方向とで同一となる補正量を算出することとしてもよい。
このようにすることで、補正後に色の明るさの差分がなくなり、人間の視覚的に色付きに対する望ましい補正量を算出することができる。
また上記様態においては、前記特定色空間が、ＨＳＶ空間であることとしてもよい。
このようにすることで、比較的に低い計算コストで色空間の算出を行うことができる。
また上記様態においては、前記特定色空間が、Ｌ＊ａ＊ｂ＊空間であることとしてもよい。
このようにすることで、人間の視覚的な感覚に近い色空間の算出を行うことができる。
また、上記態様においては、前記発生起点検出部が、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出することとしてもよい。
このようにすることで、近飽和領域、近飽和以外の領域の少なくとも一方を軸上色収差の発生起点として検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ２＜ＰｒＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計である。
このようにすることで、フォーカス位置近傍における軸上色収差の発生起点を検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
このようにすることによっても、フォーカス位置近傍における軸上色収差の発生起点を検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ２＜ＰｒＤ２
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
このようにすることで、他方向にパープルフリンジが発生し、一方向のＲの画素値の合計が、Ｇの画素値の合計より小さく、色相環で反対色に近い色となってパープルフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｒＤ２
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
このようにすることによっても、他方向にパープルフリンジが発生し、一方向のＢの画素値の合計が、Ｇの画素値の合計より小さく、色相環で反対色に近い色となってパープルフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｒＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計である。
このようにすることで、一方向にグリーンフリンジが発生し、他方向のＲの画素値の合計がＧの画素値の合計より大きく、色相環で反対色に近い色となってグリーンフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１
ＰｇＤ１＞ＰｂＤ１
ＰｇＤ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
このようにすることによっても、一方向にグリーンフリンジが発生し、他方向のＢの画素値の合計がＧの画素値の合計より大きく、色相環で反対色に近い色となってグリーンフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＋ＴＨＲ１＞ＰｒＤ１
ＰｇＤ２−ＴＨＲ２＜ＰｒＤ２
ＰｇＤ１＋ＴＨＢ１＞ＰｂＤ１
ＰｇＤ２−ＴＨＢ２＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２は任意の定数である。
このようにすることで、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２を適当に選択することにより、ノイズの影響を低減し、軸上色収差の発生起点を安定して検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＋ＴＨＲ＞ＰｒＤ１
ＰｇＤ２＋ＴＨＲ＜ＰｒＤ２
ＰｇＤ１＋ＴＨＢ＞ＰｂＤ１
ＰｇＤ２＋ＴＨＢ＜ＰｂＤ２
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ，ＴＨＢは任意の定数である。
このようにすることで、ＴＨＲ，ＴＨＢを適当に選択することにより、無彩色ではなく、色付き被写体上に発生した軸上色収差の発生起点を検出することができる。
また、上記態様においては、前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定してもよい。
ＰｇＤ１＞ＰｒＤ１×ＫＲ
ＰｇＤ２＜ＰｒＤ２×ＫＲ
ＰｇＤ１＞ＰｂＤ１×ＫＢ
ＰｇＤ２＜ＰｂＤ２×ＫＢ
ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＫＲ，ＫＢは任意の定数である。
このようにすることによっても、ＫＲ，ＫＢを適当に選択することにより、無彩色ではなく、色付き被写体上に発生した軸上色収差の発生起点を検出することができる。
また、上記態様においては、前記画像が、各画素にＲＧＢ全ての画素値を有していてもよい。
このようにすることで、各画素においてＲＧＢ全ての画素値を同時に取得する同時化された画像に対しても軸上色収差の発生起点を検出する処理を行うことができる。
また、上記態様においては、前記画像が、各画素においてＲＧＢの少なくとも１つの画素値が欠落した画素値を有していてもよい。
このようにすることで、各画素においてＲＧＢ全ての画素値が同時に取得されない、同時化されていない画像に対しても軸上色収差の発生起点を検出することができ、検出時に同時化のための補間エラーの影響をなくすことができる。
また、本発明の他の態様は、入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、を備える画像処理方法である。
上記他の態様において、前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出することとしてもよい。
また、上記他の態様において、前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出ステップにおいて算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出することとしてもよい。
また、本発明の他の態様は、入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、をコンピュータに実行させる画像処理プログラムである。
上記他の態様において、前記コンピュータが、前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出することとしてもよい。
また、上記他の態様において、前記コンピュータが、前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出ステップにおいて算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出することとしてもよい。
本発明によれば、近飽和領域もしくは近飽和領域以外の領域少なくとも一方においても軸上色収差を検出して人間の視覚特性に基づいて補正することができるという効果を奏する。
本発明の一実施形態に係る画像処理装置を示すブロック図である。 図１の画像処理装置を適用する光学系と軸上色収差との関係を模式的に示す図である。 図１の画像処理装置を適用する光学系と軸上色収差との関係を模式的に示す図である。 図１の画像処理装置を用いた画像処理方法を説明するフローチャートである。 図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図である。 図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図であり、特に、補正前の注目画素とその周辺画素との座標関係を示す。 図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図であり、特に、補正後の注目画素とその周辺画素との座標関係を示す。 図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図であり、特に、注目画素の画素位置がＤ１方向として定まっている場合の補正前の注目画素とその周辺画素との座標関係を示す。 図１の画像処理装置おいて設定される注目画素とその周辺画素との座標関係を示す図であり、特に、注目画素の画素位置がＤ１方向として定まっている場合の補正後の注目画素とその周辺画素との座標関係を示す。 図２とは異なる態様で軸上色収差が発生している場合のステップエッジの第１の例を示す図である。 図２とは異なる態様で軸上色収差が発生している場合のステップエッジの第２の例を示す図である。 画像信号にノイズが含まれる場合の対応例をノイズのない場合及び対応前の状態との比較で示す図である。 有彩色のステップエッジの場合の対応例を無彩色の場合及び対応前の状態との比較で示す図である。 有彩色のステップエッジの場合の他の対応例を無彩色の場合及び対応前の状態との比較で示す図である。 ベイヤ配列された撮像素子により取得された画像のＲを注目画素とする場合の画素配列を示す図である。 ベイヤ配列された撮像素子により取得された画像のＢを注目画素とする場合の画素配列を示す図である。 ベイヤ配列された撮像素子により取得された画像のＧＲを注目画素とする場合の画素配列を示す図である。 ベイヤ配列された撮像素子により取得された画像のＧＢを注目画素とする場合の画素配列を示す図である。
本発明の一実施形態に係る画像処理装置及び画像処理方法について、図面を参照して以下に説明する。
ここで、本実施形態に係る画像処理装置は、カメラにより取得された入力画像に対して後述する画像処理を施すものである。そして、画像処理装置１の処理対象である入力画像は、図示しない撮像光学系を通して入射された光線が固体撮像素子において電気信号に変換され、固定ノイズパターン除去及びアナログゲイン調整等が行われた後に、Ａ／Ｄ変換器でディジタル信号に変換され、各種画像処理が施された画像信号である。この入力画像は、各画素においてＲＧＢ全ての信号値が得られるフルカラー画像であり、三板撮像素子により取得された三板カラー画像、もしくは、単板撮像素子により取得された単板カラー画像を各画素でＲＧＢ値が得られるように補間処理された画像である。
本実施形態に係る画像処理装置は、図１に示すように、入力画像において軸上色収差が発生している軸上色収差領域を抽出する領域抽出部２と、領域抽出部２により抽出された軸上色収差領域に対する補正量を算出する軸上色収差補正量算出部３と、軸上色収差領域の軸上色収差を補正した出力画像を出力する収差補正部４とを備えている。
領域抽出部２は、合計画素値算出部１０、発生起点検出部１１、及び領域判定部１２を備えている。
合計画素値算出部１０は、入力画像中の注目画素を設定し、設定された注目画素を含んで一列に配列している複数の画素について、注目画素を中心として一方向（Ｄ１）及び他方向（Ｄ２）に配列している画素の画素値を夫々色毎に合計する。
発生起点検出部１１は、注目画素の画素値、又は合計画素値算出部１１により算出された注目画素を挟んで一方向に並ぶ画素の合計画素値若しくは他方向に並ぶ画素の合計画素値の少なくとも一方に基づいて、注目画素が軸上色収差の発生起点であるかを判定し、発生起点となる注目画素を検出する。
発生起点検出部５により検出する軸上色収差の発生起点とは、図２に示すように、例えば、注目画素のＧ画素値の大きさが飽和レベルになる注目画素、又は図３に示すように、注目画素を挟んで一直線に並ぶ一方向及び他方向の色が異なる場合の注目画素の少なくとも一方を意味している。
すなわち、図２及び図３に示すように、通常の光学系７では、波長の異なる光は光軸の奥行き方向に結像位置が異なることが知られている。一般にＧに対して長波長のＲは光軸方向の奥側に結像し、Ｇに対して短波長のＢは光軸方向の手前側に結像する。このＲＧＢの結像位置のずれが軸上色収差となって現れる。
無彩色のステップエッジを撮影した場合、像面ではＧがステップエッジの形状を有するが、Ｒ，Ｂは結像せずにステップエッジの形状が鈍っている。一方、像面奥では、Ｇは結像しないためステップエッジの形状が鈍っているが、Ｒが結像してステップエッジの形状を保持している。ＢはＧよりさらに鈍った形状となる。
反対に、像面手前では、Ｇは結像しないためステップエッジの形状が鈍り、Ｂが結像してステップエッジの形状を保持している。ＲはＧよりさらに鈍った形状となる。
領域判定部１２は、発生起点検出部５により検出された発生起点である注目画素の所定の周辺領域を軸上色収差領域として判定する。すなわち、注目画素が軸上色収差の発生起点として検出された場合には、軸上色収差の補正対象となるので、注目画素と注目画素を中心とする所定の周辺画素からなる領域を軸上色収差領域として判定する。軸上色収差領域としては、例えば、注目画素を中心とする５×５画素の領域や７×７画素の領域等、適宜設定することができる。
軸上色収差補正量算出部３は、色空間情報算出部１５と、色空間差分算出部１６と、補正量算出部１７とを備えている。
色空間情報算出部１５は、発生起点検出部５により発生起点であると検出された注目画素について、特定色空間上の色空間情報を、注目画素に対する一方向及び他方向の夫々について算出する。特定色空間としては、例えば、ＨＳＶ空間やｌ＊ａ＊ｂ空間などを適用することができる。そして、色空間情報は、合計画素値算出部１０により算出された、注目画素を中心として一方向（Ｄ１）及び他方向（Ｄ２）に配列している画素についての色毎の画素値の合計に基づいて算出する。
色空間差分算出部１６は、色空間情報算出部１５により算出された一方向の色空間情報と他方向の色空間情報との差分を算出する。
補正量算出部１７は、色空間差分算出部１６により算出された差分に応じて軸上色収差領域に対する補正量を算出する。
収差補正部４は、補正量算出部１７により算出された補正量を用いて軸上色収差領域を補正する。
以下に、このように構成された本実施形態に係る画像処理装置の作用としての画像処理方法について、図４のフローチャートを参照して説明する。
ステップＳ１において、合計画素値算出部１０は、図５に示すように、入力画像中に注目画素（０，０）を設定し、次のステップＳ２（合計画素値算出ステップ）において、この注目画素を中心として、図５中の縦方向、横方向、右斜め方向及び左斜め方向の４方向について、注目画素を挟んで両側に並ぶ複数画素の画素値の合計を夫々算出する（ステップＳ２）。
具体的には、数１に示す式に従って、各方向の色毎の画素値の合計を算出する。
Figure 0006633624
ここで、Ｐｃ（ｘ，ｙ）は、座標（ｘ，ｙ）に位置する画素の画素値、ＰｃＬは注目画素の左側の画素列の画素値の合計、ＰｃＲは注目画素の右側の画素列の画素値の合計、ＰｃＴは注目画素の上側の画素列の画素値の合計、ＰｃＢは注目画素の下側の画素列の画素値の合計、ＰｃＬＴは注目画素の左上側の画素列の画素値の合計、ＰｃＬＢは注目画素の左下側の画素列の画素値の合計、ＰｃＲＴは注目画素の右上側の画素列の画素値の合計、ＰｃＲＢは注目画素の右下側の画素列の画素値の合計である。また、添字ｃはカラーを意味し、ｒ、ｇ、ｂに置き換えることができ、置き換えることにより色毎の画素値の合計を算出する。また、ｉ，ｍは正の整数である。
ステップＳ３（発生起点検出ステップ）において、発生起点検出部１１は、以下の条件式（１）又は、合計画素値算出部１０により算出された画素値の合計が、条件式（２）〜（５）の全てを満たす場合に、注目画素を軸上色収差の発生起点として検出する。
条件式（１）は、注目画素のＧ画素値の大きさが任意の飽和レベル閾値ＴＨｓａｔより大きい場合に、注目画素を軸上色収差の発生起点として検出することを示す。
Ｐｇ（０，０）＞ＴＨｓａｔ …（１）
ＰｇＤ１＞ＰｒＤ１ …（２）
ＰｇＤ２＜ＰｒＤ２ …（３）
ＰｇＤ１＞ＰｂＤ１ …（４）
ＰｇＤ２＜ＰｂＤ２ …（５）
ここで、ＰｇＤ１は注目画素から一方向Ｄ１に配列している画素のＧの画素値の合計、ＰｇＤ２は注目画素から他方向Ｄ２に配列している画素のＧの画素値の合計、ＰｒＤ１は注目画素から一方向Ｄ１に配列している画素のＲの画素値の合計、ＰｒＤ２は注目画素から他方向Ｄ２に配列している画素のＲの画素値の合計、ＰｂＤ１は注目画素から一方向Ｄ１に配列している画素のＢの画素値の合計、ＰｂＤ２は注目画素から他方向Ｄ２に配列している画素のＢの画素値の合計である。
添字のＤ１は注目画素からの一方向であることを示し、例えば、数１におけるＬ、Ｔ、ＬＴ、ＬＢが該当する。また、添字のＤ２は注目画素からの他方向であることを示し、例えば、Ｒ、Ｂ、ＲＴ、ＲＢが該当する。方向の定義は任意であるが、Ｄ１としてＬ、Ｔ、ＬＴ、ＬＢを選択した場合にはＤ２はＲ、Ｂ、ＲＴ、ＲＢとなる。
なお、条件式（１）において、注目画素と任意の飽和レベル閾値ＴＨｓａｔとを比較する際に、注目画素のＧ画素のみを用いることとしたが、必ずしもこれに限られない。例えば、飽和レベル閾値と輝度値とを比較することもできる他、その他公知の方法を用いることもできる。また、近飽和以外の領域の軸上色収差の発生起点を算出する方法についても上記した例に限られず、公知の方法を用いることができる。
ステップＳ４（検出情報生成ステップ）において、発生起点検出部１１は、軸上色収差の発生起点であると検出された注目画素について「０」、発生起点ではないと検出された画素について「１」とする検出情報を生成する。発生起点検出部１１は、注目画素が、軸上色収差の発生起点であると検出された場合、その画素位置に対応してＰｒＤ１，ＰｒＤ２，ＰｇＤ１，ＰｇＤ２，ＰｂＤ１，ＰｂＤ２を保持する。
また、注目画素が、軸上色収差の発生起点であると検出された場合には、その画素位置と対応させて、図５中の縦方向、横方向、右斜め方向及び左斜め方向の４方向の何れの方向が上記条件式（２）〜（５）を満たしたのかを保持し、その方向についての一方向Ｄ１及び他方向Ｄ２を保持する。
つまり、発生起点検出部１１は、検出情報に関連する情報として、Ｄ１が、Ｌ、Ｔ、ＬＴ、ＬＢの何れであるのか、及びＤ２がＲ、Ｂ、ＲＴ、ＲＢの何れであるのかを保持しておく。
ステップＳ５（領域判定ステップ）において、領域判定部１２は、注目画素が軸上色収差の発生起点であると検出された場合、すなわち、注目画素について発生起点検出部１１により生成された検出情報が「０」であった場合、当該注目画素とその所定の周辺画素を含む領域を軸上色収差領域と判定する。軸上色収差領域として、例えば、検出情報が「０」である画素を中心とした７×７画素を含む領域の等、適宜設定することができる。
ステップＳ６（色空間情報算出ステップ）において、色空間情報算出部１５は、検出情報が「０」であった画素位置に対応するＰｒＤ１，ＰｇＤ１，ＰｂＤ１及びＰｒＤ２，ＰｇＤ２，ＰｂＤ２から特定色空間の色空間情報を算出する。
例えば、特定色空間をＨＳＶ空間とすると、ＰｒＤ１，ＰｇＤ１，ＰｂＤ１からＰｈＤ１，ＰｓＤ１，ＰｖＤ１の少なくとも１つを、ＰｒＤ２，ＰｇＤ２，ＰｂＤ２からＰｈＤ２，ＰｓＤ２，ＰｖＤ２の内Ｄ１方向と同じ色空間情報を算出する。
なお、ＰｈＤ１，ＰｈＤ２は、発生起点であると検出された注目画素におけるＤ１方向、Ｄ２方向についての角度（色相）に関する色空間情報を示し、ＰｓＤ１，ＰｓＤ２は、発生起点であると検出された注目画素におけるＤ１方向、Ｄ２方向についての強度（彩度）に関する色空間情報を示し、ＰｖＤ１，ＰｖＤ２は、発生起点であると検出された注目画素におけるＤ１方向、Ｄ２方向についての明るさ（明度）に関する色空間情報を示す。
また、特定色空間をＬ＊ａ＊ｂ＊空間とすると、ＰｒＤ１，ＰｇＤ１，ＰｂＤ１からＰｌ＊Ｄ１，Ｐａ＊Ｄ１，Ｐｂ＊Ｄ１の少なくとも１つを、ＰｒＤ２，ＰｇＤ２，ＰｂＤ２からＰｌ＊Ｄ２，Ｐａ＊Ｄ２，Ｐｂ＊Ｄ２の内Ｄ１方向と同じ色空間情報を算出する。なお、Ｐｌ＊Ｄ１，Ｐｌ＊Ｄ２は、発生起点であると検出された注目画素におけるＤ１方向、Ｄ２方向についての明るさ（明度）に関する色空間情報を示し、Ｐａ＊Ｄ１，Ｐｂ＊Ｄ１，Ｐａ＊Ｄ２，Ｐｂ＊Ｄ２は、補色に関する情報であって、角度及び強度に関する色空間情報として利用することができる。
ステップＳ７（色空間差分算出ステップ）において、色空間差分算出部１６により、２つの方向の色空間情報から色空間の差分を算出する。色区間情報算出部１６は、差分として、色の角度の差分、色の強度の差分、明るさの差分の少なくとも１つを算出する。
ステップＳ８（補正量算出ステップ）では、色空間差分算出部１６により算出された差分に基づいて、補正量を算出し、次のステップＳ９において、全ての画素について処理が行われたか判定し、全ての画素についての処理が完了していない場合には、ステップＳ１に戻り、処理が未だ行われていない画素を注目画素としてステップＳ２〜Ｓ７の工程が繰り返される。
全ての画素について処理が行われた場合には、ステップＳ１０（収差補正ステップ）において、領域判定部１２により設定された軸上色収差領域について、補正量算出部１７により算出された補正量を用いて収差補正部４により補正処理を実施し、処理を終了する。
上述したように、色空間情報として、色についての、角度、強度及び明るさの少なくとも一つを算出し、これらの色空間情報から角度の差分、強度の差分、明るさの差分の少なくとも一つを算出し、これら差分のうち何れか一又は複数を組み合わせた情報に基づいて補正量を算出し、軸上色収差領域の補正を行うことができる。
以下、角度の差分、強度の差分及び明るさの差分のうち、何れか１つまたは複数を組み合わせて補正量を算出する場合の具体例について説明する。
（角度に関する色空間情報に基づいて補正を行う場合）
色空間角度情報算出部１５は、発生起点検出部１１が検出情報に関連する情報として保持している２つの方向（Ｄ１，Ｄ２）について、角度に関する色空間情報を夫々算出する（ステップＳ６）。
色空間角度差分算出部１６は、算出した角度に関する色空間情報からその差分を算出する（ステップＳ７）。
具体的には、特定色空間がＨＳＶ空間である場合、角度に関する色空間情報として、ＰｈＤ１，ＰｈＤ２を算出し、これらの差分の絶対値Ｐａｎｇｌｅを以下の式（６）に従って算出する。
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=Ca2+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 …（１４８）
なお、補正前の画素値の角度をＰａｎｇｌｅ２、補正後の画素値の角度をＰａａｎｇｌｅ２’とし、次の式（１４９）に従って補正処理を行うことができる。
Pangle2’=Pangle2-Ca2 …（１４９）
また補正前の画素値の強度をＰｓｔｒｅｎｇｔｈ２、補正後の画素値の強度をＰｓｔｒｅｎｇｔｈ２’とし、次の式（１５０）に従って補正処理を行うことができる。
Pstrength2’=Pstrength2-Cb2 …（１５０）
また、補正前の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２、補正後の画素値の明度をＰｌｉｇｈｔｎｅｓｓ２’とし、次の式（１５１）に従って補正処理を行うことができる。
Plightness2’=Plightness2-Cc2 …（１５１）
角度、強度及び明るさに基づく色空間情報に従って軸上色収差を補正する場合、画像処理などにより色空間上で特定の色の角度の差、強度の差及び明るさの差が広がるような処理が行われる場合に特に有効である。
このように、本実施形態に係る画像処理装置及び画像処理装置によって実施される画像処理方法によれば、近飽和領域もしくは近飽和領域以外の領域少なくとも一方においても軸上色収差を検出して人間の視覚特性に基づいて補正することができる。
なお、本実施形態に係る画像処理装置１においては、上述した条件式（２）〜（５）全てを満たす場合に軸上色収差の発生起点であるとして検出することとしたが、これに代えて、図１０及び図１１に示されるように、条件式（２）及び条件式（３）を満たす場合、または、条件式（４）及び条件式（５）を満たす場合に、注目画素を軸上色収差の発生起点であるとして検出することにしてもよい。
このようにすることで、フォーカス位置より広い奥行き範囲で軸上色収差の発生起点を検出することができるという利点がある。
また、条件式（２）、（３）、（５）を全て満たす場合、または、条件式（３）、（４）、（５）を全て満たす場合に、注目画素を軸上色収差の発生起点として検出してもよい。
このようにすることで、他方向Ｄ２にパープルフリンジが発生し、一方向Ｄ１のＲの画素値の合計またはＢの画素値の合計が、Ｇの画素値の合計より小さく、色相環で反対色に近い色となってパープルフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。
また、条件式（２）、（３）、（４）を全て満たす場合、または、条件式（２）、（４）、（５）を全て満たす場合に、注目画素を軸上色収差の発生起点として検出してもよい。
このようにすることで、一方向Ｄ１にグリーンフリンジが発生し、他方向Ｄ２のＲの画素値の合計またはＢの画素値の合計が、Ｇの画素値の合計より大きく、色相環で反対色に近い色となってグリーンフリンジが目立つ場所を軸上色収差の発生起点として検出することができる。
また、数１で算出された画素値の合計が、次の条件式（１５２）〜（１５５）を全て満たす場合に注目画素を軸上色収差の発生起点として検出することにしてもよい。
ＰｇＤ１＋ＴＨＲ１＞ＰｒＤ１ …（１５２）
ＰｇＤ２−ＴＨＲ２＜ＰｒＤ２ …（１５３）
ＰｇＤ１＋ＴＨＢ１＞ＰｂＤ１ …（１５４）
ＰｇＤ２−ＴＨＢ２＜ＰｂＤ２ …（１５５）
ここで、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２は任意の定数である。
ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２としては、図１２に示すように、例えば、処理対象の画素位置におけるＧとＲ、または、ＧとＢのノイズ量を示す標準偏差の差を例示することができる。一方向Ｄ１及び他方向Ｄ２のＲ，Ｇ，Ｂ夫々のノイズ量をＮｒＤ１、ＮｇＤ１、ＮｂＤ１、ＮｒＤ２、ＮｇＤ２、ＮｂＤ２として、以下の式（１５６）〜式（１５９）で示すことができる。
ＴＨＲ１＝ＮｒＤ１−ＮｇＤ１ …（１５６）
ＴＨＲ２＝ＮｒＤ２−ＮｇＤ２ …（１５７）
ＴＨＢ１＝ＮｂＤ１−ＮｇＤ１ …（１５８）
ＴＨＢ２＝ＮｂＤ２−ＮｇＤ２ …（１５９）
このようにすることで、ノイズの影響を低減し、安定して軸上色収差の発生起点を検出することができるという利点がある。
また、数１で算出された画素値の合計が、次の条件式（１６０）〜（１６３）を全て満たす場合に注目画素を軸上色収差の発生起点として検出することにしてもよい。
ＰｇＤ１＋ＴＨＲ＞ＰｒＤ１ …（１６０）
ＰｇＤ２＋ＴＨＲ＜ＰｒＤ２ …（１６１）
ＰｇＤ１＋ＴＨＢ＞ＰｂＤ１ …（１６２）
ＰｇＤ２＋ＴＨＢ＜ＰｂＤ２ …（１６３）
ここで、ＴＨＲ，ＴＨＢは任意の定数である。
定数ＴＨＲ，ＴＨＢとしては、図１３に示されるように、例えば、処理対象の画素位置におけるＧとＲ、または、ＧとＢの画素値の差分を例示することができる。注目画素のＲ，Ｇ，Ｂの夫々の画素値をＰｒ（０，０）、Ｐｇ（０，０）、Ｐｂ（０，０）とすると、以下の式（１６４）及び式（１６５）となる。
ＴＨＲ＝Ｐｇ（０，０）−−Ｐｒ（０，０） …（１６４）
ＴＨＢ＝Ｐｇ（０，０）−−Ｐｂ（０，０） …（１６５）
このようにすることで、無彩色ではなく、色つきの被写体上に発生した軸上色収差の発生起点を検出することができるという利点がある。
また、数１で算出された画素値の合計が、次の条件式（１６６）〜（１６９）を全て満たす場合に注目画素を軸上色収差の発生起点として検出することにしてもよい。
ＰｇＤ１＞ＰｒＤ１×ＫＲ …（１６６）
ＰｇＤ２＜ＰｒＤ２×ＫＲ …（１６７）
ＰｇＤ１＞ＰｂＤ１×ＫＢ …（１６８）
ＰｇＤ２＜ＰｂＤ２×ＫＢ …（１６９）
ここで、ＫＲ，ＫＢは任意の定数である。
定数ＫＲ，ＫＢとしては、図１４に示されるように、例えば、処理対象の画素位置におけるＧとＲ、または、ＧとＢの画素値の比率を例示することができる。注目画素のＲ，Ｇ，Ｂの夫々の画素値をＰｒ（０，０）、Ｐｇ（０，０）、Ｐｂ（０，０）とすると、以下の式（１７０）及び式（１７１）となる。
ＫＲ＝Ｐｇ（０，０）／Ｐｒ（０，０） …（１７０）
ＫＢ＝Ｐｇ（０，０）／Ｐｂ（０，０） …（１７１）
このようにすることによっても、無彩色ではなく、色つきの被写体上に発生した軸上色収差の発生起点を検出することができるという利点がある。
また、本実施形態においては、各画素においてＲＧＢ全ての信号値が得られるフルカラー画像が入力画像として入力される場合について例示したが、これに代えて、いわゆるベイヤ画像である場合に適用することにしてもよい。この場合、一方向Ｄ１の合計画素値算出部４における合計画素値の算出方法は数２〜数４の通りとなる。
図１５に示すように、注目画素がＲ（斜線部分）である場合には、数２及び数３に示す数式により合計画素値が算出される。
Figure 0006633624
Figure 0006633624
図１６に示されるように、注目画素がＢ（斜線部分）である場合には、数４及び数５に示される数式により合計画素値が算出される。
Figure 0006633624
Figure 0006633624
図１７に示されるように、注目画素がＧＲ（斜線部分）である場合には、数６及び数７に示される数式により合計画素値が算出される。
Figure 0006633624
Figure 0006633624
図１８に示されるように、注目画素がＧＢ（斜線部分）である場合には、数８及び数９に示される数式により合計画素値が算出される。
Figure 0006633624
Figure 0006633624
このようにすることで、同時化されていない画像に対しても軸上色収差の発生起点を検出することができ、検出時に同時化のための補間エラーの影響をなくすことができる。
本実施形態に係る画像処理装置１は、放送用据え置き型カメラ、ＥＮＧカメラ、民生用ハンディカメラ、デジタルカメラ等の製品に搭載して使用される。また、動画を扱う画像信号補正プログラム（ＣＧプログラム）や画像編集装置等にも用いることができる。また、本実施形態に係る画像処理方法をコンピュータ（例えば、パーソナルコンピュータ）に実行させる画像処理プログラムとして用いることもできる。
２ 領域抽出部
３ 軸上色収差補正量算出部
４ 収差補正部
１０ 合計画素値算出部
１１ 発生起点検出部
１２ 領域判定部
１５ 色空間情報算出部
１６ 色空間差分算出部
１７ 補正量算出部

Claims (20)

  1. 入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出部と、
    該合計画素値算出部による算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素が軸上色収差の発生起点であるか否かを判定し、発生起点を検出する発生起点検出部と、
    前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定部と、
    前記発生起点として検出された前記注目画素に対し、前記合計画素値算出部による算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出部と、
    該色空間情報算出部により算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出部と、
    該色空間差分算出部により算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出部と、
    前記軸上色収差領域について前記補正量を用いて補正する補正部と、
    を備える画像処理装置。
  2. 前記色空間情報算出部が、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
    前記色空間差分算出部が、前記色空間情報算出部による算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
    前記補正量算出部が、前記色空間差分算出部による算出結果に基づいて軸上色収差補正の補正量を算出する請求項１記載の画像処理装置。
  3. 前記発生起点検出部が、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出部により算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項１または２記載の画像処理装置。
  4. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＞ＰｒＤ１
    ＰｇＤ２＜ＰｒＤ２
    ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計である。
  5. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＞ＰｂＤ１
    ＰｇＤ２＜ＰｂＤ２
    ここで、
    ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
  6. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＞ＰｒＤ１
    ＰｇＤ２＜ＰｒＤ２
    ＰｇＤ２＜ＰｂＤ２
    ここで、
    ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
  7. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＞ＰｂＤ１
    ＰｇＤ２＜ＰｒＤ２
    ＰｇＤ２＜ＰｂＤ２
    ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計である。
  8. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＞ＰｒＤ１
    ＰｇＤ１＞ＰｂＤ１
    ＰｇＤ２＜ＰｒＤ２
    ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、である。
  9. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＞ＰｒＤ１
    ＰｇＤ１＞ＰｂＤ１
    ＰｇＤ２＜ＰｂＤ２
    ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計
    である。
  10. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＋ＴＨＲ１＞ＰｒＤ１
    ＰｇＤ２−−ＴＨＲ２＜ＰｒＤ２
    ＰｇＤ１＋ＴＨＢ１＞ＰｂＤ１
    ＰｇＤ２−−ＴＨＢ２＜ＰｂＤ２
    ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ１，ＴＨＲ２，ＴＨＢ１，ＴＨＢ２は任意の定数である。
  11. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＋ＴＨＲ＞ＰｒＤ１
    ＰｇＤ２＋ＴＨＲ＜ＰｒＤ２
    ＰｇＤ１＋ＴＨＢ＞ＰｂＤ１
    ＰｇＤ２＋ＴＨＢ＜ＰｂＤ２
    ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＴＨＲ，ＴＨＢは任意の定数である。
  12. 前記発生起点検出部が、以下の条件式を満足する場合に、前記注目画素を軸上色収差の発生起点であると判定する請求項１または２記載の画像処理装置。
    ＰｇＤ１＞ＰｒＤ１×ＫＲ
    ＰｇＤ２＜ＰｒＤ２×ＫＲ
    ＰｇＤ１＞ＰｂＤ１×ＫＢ
    ＰｇＤ２＜ＰｂＤ２×ＫＢ
    ここで、ＰｇＤ１は前記注目画素から前記一方向に配列している画素のＧの画素値の合計、ＰｇＤ２は前記注目画素から前記他方向に配列している画素のＧの画素値の合計、ＰｒＤ１は前記注目画素から前記一方向に配列している画素のＲの画素値の合計、ＰｒＤ２は前記注目画素から前記他方向に配列している画素のＲの画素値の合計、ＰｂＤ１は前記注目画素から前記一方向に配列している画素のＢの画素値の合計、ＰｂＤ２は前記注目画素から前記他方向に配列している画素のＢの画素値の合計、ＫＲ，ＫＢは任意の定数である。
  13. 前記画像が、各画素にＲＧＢ全ての画素値を有する請求項１または２記載の画像処理装置。
  14. 前記画像が、各画素においてＲＧＢの少なくとも１つの画素値が欠落した画素値を有する請求項１または２記載の画像処理装置。
  15. 入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、
    該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、
    前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、
    前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、
    該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、
    該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、
    前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、
    を備える画像処理方法。
  16. 前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
    前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
    前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出する請求項１５記載の画像処理方法。
  17. 前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出ステップにおいて算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項１５又は１６記載の画像処理方法。
  18. 入力画像中の注目画素を含み一列に配列する複数の画素について、該注目画素を挟んで一方向及び他方向に配列する画素の画素値を夫々色毎に合計する合計画素値算出ステップと、
    該合計画素値算出ステップにおける算出結果又は前記注目画素の画素値の少なくとも一方に基づいて、前記注目画素を軸上色収差の発生起点として検出する発生起点検出ステップと、
    前記発生起点として検出された注目画素を中心とする所定の周辺領域を軸上色収差領域と判定する領域判定ステップと、
    前記発生起点として検出された前記注目画素に対し、前記合計画素値算出ステップにおける算出結果に基づいて、前記一方向及び前記他方向について特定色空間上の色空間情報を夫々算出する色空間情報算出ステップと、
    該色空間情報算出ステップにおいて算出された前記一方向の色空間情報と前記他方向の色空間情報との差分を算出する色空間差分算出ステップと、
    該色空間差分算出ステップにおいて算出された差分に応じて軸上色収差補正の補正量を算出する補正量算出ステップと、
    前記軸上色収差領域について前記補正量を用いて補正する補正ステップと、
    をコンピュータに実行させる画像処理プログラム。
  19. 前記コンピュータが、
    前記色空間情報算出ステップにおいて、特定色空間上の色空間情報として、色の角度に関する情報、色の強度に関する情報、及び色の明るさに関する情報の少なくとも１つを算出し、
    前記色空間差分算出ステップにおいて、前記色空間情報算出ステップによる算出結果に基づいて、前記一方向の色空間情報と前記他方向の色空間情報との差分として、前記色の角度の差分、前記色の強度の差分、及び前記色の明るさの差分のうち少なくとも１つを算出し、
    前記補正量算出ステップにおいて、前記色空間差分算出ステップによる算出結果に基づいて軸上色収差補正の補正量を算出する請求項１８記載の画像処理プログラム。
  20. 前記コンピュータが、
    前記発生起点検出ステップにおいて、前記注目画素の画素値が所定の閾値よりも大きい場合、又は、前記合計画素値算出ステップにおいて算出された前記一方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＲの画素値の合計との大小関係と逆である場合、若しくは、前記一方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係が、前記他方向に配列している画素のＧの画素値の合計とＢの画素値の合計との大小関係と逆である場合に、前記注目画素を軸上色収差の発生起点として検出する請求項１８記載の画像処理プログラム。
JP2017517499A 2015-05-12 2015-05-12 画像処理装置、画像処理方法及び画像処理プログラム Active JP6633624B2 (ja)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2015/063540 WO2016181476A1 (ja) 2015-05-12 2015-05-12 画像処理装置、画像処理方法及び画像処理プログラム

Publications (2)

Publication Number Publication Date
JPWO2016181476A1 JPWO2016181476A1 (ja) 2018-03-01
JP6633624B2 true JP6633624B2 (ja) 2020-01-22

Family

ID=57247848

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2017517499A Active JP6633624B2 (ja) 2015-05-12 2015-05-12 画像処理装置、画像処理方法及び画像処理プログラム

Country Status (4)

Country Link
US (1) US10659738B2 (ja)
JP (1) JP6633624B2 (ja)
CN (1) CN107534758B (ja)
WO (1) WO2016181476A1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106454054B (zh) 2016-11-29 2019-03-19 Oppo广东移动通信有限公司 控制方法、控制装置及电子装置
CN106341670B (zh) 2016-11-29 2017-09-22 广东欧珀移动通信有限公司 控制方法、控制装置及电子装置
CN106454288B (zh) 2016-11-29 2018-01-19 广东欧珀移动通信有限公司 控制方法、控制装置、成像装置及电子装置
CN106504218B (zh) 2016-11-29 2019-03-12 Oppo广东移动通信有限公司 控制方法、控制装置及电子装置
GB2568261B (en) * 2017-11-08 2022-01-26 Displaylink Uk Ltd System and method for presenting data at variable quality
US11442266B1 (en) * 2019-09-09 2022-09-13 Apple Inc. Method and device for correcting chromatic aberration in multiple bands

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4539278B2 (ja) 2004-10-13 2010-09-08 ソニー株式会社 画像処理装置、および画像処理方法、並びにコンピュータ・プログラム
WO2007007878A1 (ja) 2005-07-14 2007-01-18 Nikon Corporation 画像処理装置および画像処理方法
JP4487902B2 (ja) 2005-11-09 2010-06-23 ソニー株式会社 画像処理装置および方法、並びにプログラム
US7577292B2 (en) 2005-12-30 2009-08-18 Microsoft Corporation Automatic removal of purple fringing from images
JP4192975B2 (ja) * 2006-07-26 2008-12-10 ソニー株式会社 画像処理装置および方法、プログラム並びに記録媒体
JP5200645B2 (ja) 2008-04-24 2013-06-05 株式会社ニコン 画像処理装置、撮像装置、画像処理方法およびプログラム
US8958009B2 (en) 2010-01-12 2015-02-17 Nikon Corporation Image-capturing device
JP5163733B2 (ja) 2010-06-17 2013-03-13 株式会社ニコン 撮像装置、および焦点判定プログラム
JP5324508B2 (ja) 2010-03-26 2013-10-23 富士フイルム株式会社 画像処理装置および方法，ならびに画像処理プログラム
JP5552956B2 (ja) 2010-08-12 2014-07-16 株式会社ニコン 画像処理装置、画像処理方法およびプログラム
US8553978B2 (en) * 2010-11-03 2013-10-08 Csr Technology Inc. System and method for providing multi resolution purple fringing detection and correction
JP5840008B2 (ja) 2011-02-21 2016-01-06 キヤノン株式会社 画像処理装置、画像処理方法およびプログラム
JP2012209709A (ja) * 2011-03-29 2012-10-25 Sony Corp 画像処理装置および方法、記録媒体並びにプログラム
JP5847471B2 (ja) 2011-07-20 2016-01-20 キヤノン株式会社 画像処理装置、撮像装置、画像処理方法および画像処理プログラム
JP2013172402A (ja) * 2012-02-22 2013-09-02 Nikon Corp 画像処理装置、撮像装置および画像処理プログラム
JP2014086932A (ja) 2012-10-25 2014-05-12 Sony Corp 画像処理装置、画像処理方法、および画像処理プログラム
JP6441379B2 (ja) 2014-11-28 2018-12-19 オリンパス株式会社 画像処理装置、画像処理方法および画像処理プログラム

Also Published As

Publication number Publication date
US10659738B2 (en) 2020-05-19
CN107534758B (zh) 2019-05-28
WO2016181476A1 (ja) 2016-11-17
US20180077392A1 (en) 2018-03-15
JPWO2016181476A1 (ja) 2018-03-01
CN107534758A (zh) 2018-01-02

Similar Documents

Publication Publication Date Title
JP6633624B2 (ja) 画像処理装置、画像処理方法及び画像処理プログラム
KR102170686B1 (ko) 화이트 밸런스 보정 장치 및 방법
JP6600936B2 (ja) 画像処理装置、画像処理方法、画像処理システム、プログラム、及び記録媒体
JP4346634B2 (ja) 目標物検出装置
JP2016510108A5 (ja)
JP2015211319A5 (ja) 画像処理装置、画像処理方法、およびプログラム
JP6239358B2 (ja) 画素補間装置、撮像装置、プログラムおよび集積回路
JP2014086932A (ja) 画像処理装置、画像処理方法、および画像処理プログラム
JP6282123B2 (ja) 画像処理装置、画像処理方法及びプログラム
JP5718138B2 (ja) 画像信号処理装置及びプログラム
WO2009098827A1 (ja) 画像処理装置及び画像処理プログラム
US10489890B2 (en) Image-processing apparatus, image-processing method, and image-processing program
JP6415093B2 (ja) 画像処理装置、画像処理方法、及びプログラム
KR101810393B1 (ko) 화상처리장치, 촬상장치, 화상처리방법 및 기억매체
JP2015211343A5 (ja)
CN107340296A (zh) 变化程度导出设备、变化程度导出系统和已知颜色体
JP2013038655A (ja) 画像処理装置及びその制御方法
JP6552248B2 (ja) 画像処理装置及び方法、撮像装置、並びにプログラム
JP4334150B2 (ja) 画像補間装置
JP5836878B2 (ja) 画像処理装置、方法、及びプログラム
JP2014158165A (ja) 画像処理装置、画像処理方法およびプログラム
JP6725060B2 (ja) 画像処理装置、画像処理システム、画像処理方法、及びプログラム
JP4334152B2 (ja) 画像補間装置
JP2021170767A (ja) 画像処理装置、画像処理方法、および画像処理システム
JP6033182B2 (ja) 画像読取装置

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20171108

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180425

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20190611

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20190718

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20191119

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20191212

R151 Written notification of patent or utility model registration

Ref document number: 6633624

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250