JP3861808B2 - 撮像装置 - Google Patents
撮像装置 Download PDFInfo
- Publication number
- JP3861808B2 JP3861808B2 JP2002375423A JP2002375423A JP3861808B2 JP 3861808 B2 JP3861808 B2 JP 3861808B2 JP 2002375423 A JP2002375423 A JP 2002375423A JP 2002375423 A JP2002375423 A JP 2002375423A JP 3861808 B2 JP3861808 B2 JP 3861808B2
- Authority
- JP
- Japan
- Prior art keywords
- color
- signal
- processing unit
- filter
- linear matrix
- 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.)
- Expired - Lifetime
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/80—Camera processing pipelines; Components thereof
- H04N23/84—Camera processing pipelines; Components thereof for processing colour signals
- H04N23/843—Demosaicing, e.g. interpolating colour pixel values
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N23/00—Cameras or camera modules comprising electronic image sensors; Control thereof
- H04N23/10—Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths
- H04N23/12—Cameras or camera modules comprising electronic image sensors; Control thereof for generating image signals from different wavelengths with one sensor only
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Color Television Image Signal Generators (AREA)
- Facsimile Image Signal Circuits (AREA)
- Image Processing (AREA)
- Processing Of Color Television Signals (AREA)
- Color Image Communication Systems (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、被写体を撮像する撮像装置及び方法に関し、特に撮像感度に対して好適に撮像を行なう撮像装置及び方法に関する。
【０００２】
【従来の技術】
近年、コンシューマ向けの画像撮像装置（ディジタルカメラ、カラースキャナ等）および画像処理ソフトウェアが普及しており、当該画像撮像装置により被写体を撮像して生成した画像を画像処理ソフトウェアにより自ら編集するユーザが増加している。また、画像撮像装置により撮像された画像の画質に対する要求も高くなってきており、画像撮像装置を購入する際の主要な条件として、画質がよいことを挙げるユーザの割合が増えてきている。ここで、画像撮像装置の一般的な構成について以下に述べる。
【０００３】
画像撮像装置は、例えば、図１に示すようなＲＧＢの３原色のカラーフィルタ１が用いられている。この例では、図１の一点鎖線で示すように、緑色（Ｇ）の光のみを透過するＧフィルタが２個、赤（Ｒ）の光のみを透過するＲフィルタが１個、および青（Ｂ）の光のみを透過するＢフィルタが１個の、合計４個を最小単位として、いわゆるベイヤー配列（Bayer配列）により、カラーフィルタ１が構成される。
【０００４】
図２は、ＲＧＢカラーフィルタ１を有するＣＣＤ(Charge Coupled Device)撮像素子により取得されたＲＧＢ信号に対して各種の処理を施す信号処理部１１の構成例を示すブロック図である。
【０００５】
オフセット補正処理部２１は、ＣＣＤ撮像素子により取得された信号に所定の処理を施すフロントエンド１３から供給されてきた画像信号に含まれるオフセット成分を除去し、得られた画像信号をホワイトバランス補正処理部２２に出力する。ホワイトバランス補正処理部２２は、オフセット補正処理部２１から供給されてきた画像信号の色温度、およびカラーフィルタ１の各フィルタの感度の違いに基づいて、各色のバランスを補正する。ホワイトバランス補正処理部２２により補正が施され、取得された色信号はガンマ補正処理部２３に出力される。ガンマ補正処理部２３は、ホワイトバランス補正処理部２２から供給されてきた信号に対してガンマ補正を行い、取得した信号を垂直方向同時化処理部２４に出力する。垂直方向同時化処理部２４にはディレイ素子が設けられており、ガンマ補正処理部２３から供給されてきた信号の垂直方向の時間のずれが同時化される。
【０００６】
ＲＧＢ信号生成処理部２５は、垂直方向同時化処理部２４から供給されてきた色信号を、同一空間の位相に補間する補間処理、信号のノイズ成分を除去するノイズ除去処理、信号帯域を制限するフィルタリング処理、および信号帯域の高域成分を補正する高域周波数補正処理等を行い、得られたＲＧＢ信号を、輝度信号生成処理部２６、および色差信号生成処理部２７に出力する。
【０００７】
輝度信号生成処理部２６は、ＲＧＢ信号生成処理部２５から供給されてきたＲＧＢ信号を所定の合成比で合成し、輝度信号（Ｙ）を生成する。色差信号生成処理部２７も同様に、ＲＧＢ信号生成処理部２５から供給されてきたＲＧＢ信号を所定の合成比で合成し、色差信号（Ｃｂ，Ｃｒ）を生成する。輝度信号生成処理部２６により生成された輝度信号（Ｙ）、および色差信号生成処理部２７により生成された色差信号（Ｃｂ，Ｃｒ）は、例えば、信号処理部１１の外部に設けられているモニタに出力される。
【０００８】
このように、元信号に対してガンマ処理を施してから線形変換によって画像処理（リニアマトリクス処理）を行うことが一般的に行われている。
【０００９】
上述したような画像撮像装置は、被写体を撮像し画像を生成したとき、観察時の視環境によって見え方が異なるために、所望の色で再現されない場合がある。これは、撮像時の光源（以下、撮像光源という。）と観察時の光源（以下、観察光源という。）の演色性が大きく異なる場合に生じる現象である。そこで、撮像光源と演色性の異なる観察光源で画像が再現される場合でも、良好に色再現する方法が提案されている（例えば、特許文献１参照。）。
【００１０】
また、図３に分光感度の特性を示す。図３の曲線Ｌ１は、Ｒの分光感度を示し、曲線Ｌ２は、Ｇの分光感度を示し、曲線Ｌ３は、Ｂの分光感度を示している。
【００１１】
【特許文献１】
特開２００２−１４２２３１号公報
【００１２】
【発明が解決しようとする課題】
一方、上述したような画像撮像装置に、リニアマトリクス処理の追加、更に撮像素子のカラーフィルタの追加を行うことにより、色再現性を劇的に向上させることができる。その際、使用するリニアマトリクスの係数を単に色差が最小になるように決定すると、ノイズを従来よりも増長させてしまうことがある。
【００１３】
本発明は、上述したような実情に鑑みて提案されたものであり、撮像する環境や条件等に応じて色再現性とノイズ低減性を考慮した係数を用いたリニアマトリクス処理を行なうことが可能な撮像装置及び撮像方法を提供することを目的とする。
【００１４】
【課題を解決するための手段】
本発明に係る撮像装置は、上述の課題を解決するために、分光特性の異なるカラーフィルタによりなっており、被写体を撮像する撮像素子部を有し、該撮像素子部で撮像した画像にマトリクス変換処理を行って色補正を行う撮像装置において、色再現性指標ΔＥ（Ｍ）とノイズ低減性指標σＮ（Ｍ）とから評価指数ＥＥＶ（Ｍ）を以下の条件により導出し、導出された上記評価指数ＥＥＶ（Ｍ）に基づいてマトリクス係数を決定し、決定された上記マトリクス係数に基づいて上記撮像素子部で撮像した画像にマトリクス変換処理を行うマトリクス変換処理手段を備える。
（条件）
ＥＥＶ（Ｍ）＝√｛（ｗｃ・ΔＥ（Ｍ）） ２ ＋（ｗｎ・σＮ（Ｍ）） ２ ｝
但し、ＥＥＶ（Ｍ）：評価指数、ΔＥ（Ｍ）：色再現性指標、σＮ（Ｍ）：ノイズ低減性指標、ｗｃ：色再現性に対する重み係数、ｗｎ：ノイズ低減性に対する重み係数
【００１６】
【発明の実施の形態】
図５は、本発明を適用した撮像装置の構成例を示すブロック図である。
【００１７】
図５に示す撮像装置は、CCD(Charge Coupled Device)等よりなる画像センサ４５の前面（レンズ４２に対向する面）に、４種類の色（光）を識別するカラーフィルタが設けられている。なお、撮像装置では、図５の画像センサ４５に設けられるカラーフィルタを図６に示す４色カラーフィルタ６１とする。
【００１８】
４色カラーフィルタ６１は、図６の一点鎖線で示すように、赤の光のみを透過するＲフィルタ、青の光のみを透過するＢフィルタ、第１の波長帯域の緑色の光のみを透過するＧ１フィルタ、およびＧ１フィルタと相関が高い、第２の波長帯域の緑色の光のみを透過するＧ２フィルタの、合計４個のフィルタを最小単位として構成される。また、Ｇ１フィルタとＧ２フィルタは、その最小単位内において、互いに対角する位置に配置されている。
【００１９】
後に詳述するように、画像センサ４５により取得される画像の色を４種類とし、取得される色情報を増やすことにより、３種類の色（RGB）のみが取得される場合と較べて、より正確に色を表現することができ、眼で違うものに見える色は違う色に、同じものに見える色は同じ色に、それぞれ再現すること（「色の判別性」）を向上させることができる。
【００２０】
なお、図７に示される視感度曲線から分かるように、人間の眼は輝度に敏感である。従って、図６に示す４色カラーフィルタ６１においては、より正確な輝度情報を取得することにより、輝度の階調を上げることができ、かつ、眼の見え方に近い画像を再現することができるように、視感度曲線に近い分光感度特性を有するＧ２のカラーフィルタが追加されている（図１のＲ，Ｇ，Ｂに対応するＲ，Ｇ１，Ｂのフィルタに対して、新たに決定された緑Ｇ２フィルタが追加されている）。
【００２１】
また、４色カラーフィルタ６１を決定する際に使用されるフィルタ評価係数として、例えば、「色再現性」と「ノイズ低減性」の両方を考慮した係数であるUMG(Unified Measure of Goodness)が用いられている。
【００２２】
UMGを用いた評価においては、評価対象のフィルタが単にルータ条件を満たしているだけでは、その評価値は高くならず、それぞれのフィルタの分光感度分布の重なりも考慮される。従って、ｑファクタ、μファクタ、またはFOM(Figure of Merit)を利用して評価されたカラーフィルタの場合と比較して、ノイズをより低減させることができる。すなわち、UMGを用いた評価により、それぞれのフィルタの分光感度特性が、ある程度の重なりを有するが、図４のＲの特性とＧの特性のようにほぼ全てが重なるものではないフィルタが選択されるため、色の分離のため、それぞれの色信号を増幅する場合であっても、増幅率をそれ程大きくする必要がなく、それに伴ってノイズ成分が増幅されることが抑制される。
【００２３】
図８は、各フィルタ評価係数の特徴を示す図である。なお、図８は、各評価係数に対して、一度に評価できるフィルタの数、物体の分光反射率が考慮されているか否か、およびノイズの低減が考慮されているか否かに関する事項を示すものである。
【００２４】
図８に示すように、ｑファクタ（q-factor）は、一度に評価できるフィルタの数が「１個」のみであり、物体の分光反射率と、ノイズの低減が考慮されていない。また、μファクタ（μ-factor）は、一度に複数のフィルタを評価することができるものの、物体の分光反射率とノイズの低減は考慮されていない。さらに、FOMは、一度に複数のフィルタを評価することができ、物体の分光反射率が考慮されているものの、ノイズの低減が考慮されていない。
【００２５】
これに対して、４色カラーフィルタ６１を決定する際に使用されるUMGは、一度に複数のフィルタを評価することができ、物体の分光反射率が考慮され、かつ、ノイズの低減も考慮されている。
【００２６】
なお、ｑファクタについては「H.E.J.Neugebauer "Quality Factor for Filters Whose Spectral Transmittances are Different from Color Mixture Curves, and Its Application to Color Photography" JOURNAL OF THE OPTICAL SOCIETY OF AMERICA, VOLUME 46, NUMBER 10」にその詳細が開示されており、μファクタについては「P.L.Vora and H.J.Trussell, "Measure of Goodness of a set of color-scanning filters", JOURNAL OF THE OPTICAL SOCIETY OF AMERICA, VOLUME 10, NUMBER 7」にその詳細が開示されている。また、FOMについては「G.Sharma and H.J.Trussell, "Figures of Merit for Color Scanners", IEEE TRANSACTION ON IMAGE PROCESSING, VOLUME 6」にその詳細が開示されており、UMGについては「S.Quan, N.Ohta, and N.Katoh, "Optimal Design of Camera Spectral Sensitivity Functions Based on Practical Filter Components", CIC, 2001」にその詳細が開示されている。
【００２７】
図５の説明に戻る。マイクロコンピュータ４１は、所定の制御プログラムに従って全体の動作を制御する。例えば、マイクロコンピュータ４１は、絞り４３による露光制御、シャッター４４の開閉制御、TG(Timing Generator)４６の電子シャッターの制御、フロントエンド４７でのゲインコントロール、カメラシステムLSI(Large Scale Integrated Circuit)４８のモード制御、パラメータ制御等を行う。
【００２８】
絞り４３は、レンズ４２により集光された光の通過（絞り）を調整し、画像センサ４５により取り込まれる光量を制御する。シャッター４４は、マイクロコンピュータ４１の指示に基づいて、レンズ４２により集光された光の通過を制御する。
【００２９】
画像センサ４５は、さらに、CCDやCMOS(Complementary Metal Oxide Semiconductor)により構成される撮像素子を有し、その撮像素子の前面に形成する４色カラーフィルタ６１を介して入射された光を電気信号に変換し、４種類の色信号（Ｒ信号，Ｇ１信号，Ｇ２信号，Ｂ信号）をフロントエンド４７に出力する。画像センサ４５には、図６の４色カラーフィルタ６１が設けられ、レンズ４２を介して入射された光から、Ｒ，Ｇ１，Ｇ２，Ｂのそれぞれの帯域の波長の成分が抽出される。なお、その詳細については、図１５を参照して後述する。
【００３０】
フロントエンド４７は、画像センサ４５から供給されてきた色信号に対して、ノイズ成分を除去するための相関二重サンプリング処理、ゲインコントロール処理、およびディジタル変換処理等を施す。フロントエンド４７により各種の処理が施され、得られた画像データは、カメラシステムLSI４８に出力される。
【００３１】
カメラシステムLSI４８は、後に詳述するように、フロントエンド４７から供給されてきた画像データに対して各種の処理を行い、例えば、輝度信号および色信号を生成して画像モニタ５０に出力し、信号に対応する画像を表示させる。
【００３２】
画像メモリ４９は、例えば、DRAM(Dynamic Random Access Memory)やSDRAM(Synchronous Dynamic Random Access Memory)などにより構成され、カメラシステムLSI４８が各種の処理を行う際に適宜利用される。半導体メモリ、ディスク等により構成される外部記憶媒体５１は、例えば、図５の撮像装置に対して着脱可能に構成され、カメラシステムLSI４８によりJPEG(Joint Photographic Expert Group)フォーマットで圧縮された画像データが記憶される。
【００３３】
画像モニタ５０は、例えば、LCD(Liquid Crystal Display)などにより構成され、撮影された画像や各種のメニュー画面等を表示する。
【００３４】
図９は、図５に示すカメラシステムLSI４８の構成例を示すブロック図である。カメラシステムLSI４８を構成する各ブロックは、マイクロコンピュータインタフェース(I/F)７３を介して、図５に示すマイクロコンピュータ４１により制御される。
【００３５】
信号処理部７１は、フロントエンド４７から供給されてきた４種類の色情報に対して、補間処理、フィルタリング処理、マトリクス演算処理、輝度信号生成処理、色差信号生成処理等の各種の処理を行い、例えば、生成した画像信号をモニタインタフェース７７を介して画像モニタ５０に出力する。
【００３６】
画像検波部７２は、フロントエンド４７の出力に基づいて、オートフォーカス、オートエキスポージャー、オートホワイトバランス等の検波処理を行い、その結果を、適宜、マイクロコンピュータ４１に出力する。
【００３７】
メモリコントローラ７５は、処理ブロック間同士のデータの送受信、または、所定の処理ブロックと画像メモリ４９との間のデータの送受信を制御し、例えば、信号処理部７１から供給されてきた画像データをメモリインタフェース７４を介して画像メモリ４９に出力し、記憶させる。
【００３８】
画像圧縮・解凍部７６は、例えば、信号処理部７１から供給されてきた画像データをJPEGフォーマットで圧縮し、得られたデータをマイクロコンピュータインタフェース７３を介して外部記憶媒体５１に出力し、記憶させる。画像圧縮・解凍部７６は、また、外部記憶媒体５１から読み出された圧縮データを解凍（伸張）し、モニタインタフェース７７を介して画像モニタ５０に出力する。
【００３９】
図１０は、図９に示す信号処理部７１の詳細な構成例を示すブロック図である。信号処理部７１を構成する各ブロックは、マイクロコンピュータインタフェース７３を介して、マイクロコンピュータ４１により制御される。
【００４０】
オフセット補正処理部９１は、フロントエンド４７から供給されてきた画像信号に含まれるノイズ成分（オフセット成分）を除去し、得られた画像信号をホワイトバランス補正処理部９２に出力する。ホワイトバランス補正処理部９２は、オフセット補正処理部９１から供給されてきた画像信号の色温度、および４色カラーフィルタ６１の各フィルタの感度の違いに基づいて、各色のバランスを補正する。ホワイトバランス補正処理部９２により補正が施され、取得された色信号は、垂直方向同時化処理部９３に出力される。垂直方向同時化処理部９３にはディレイ素子が設けられており、ホワイトバランス補正処理部９２から出力されてきた信号（以下、RG1G2B信号という。）の垂直方向の時間のずれが同時化（補正）される。
【００４１】
信号生成処理部９４は、垂直方向同時化処理部９３から供給されてきたRG1G2B信号の最小単位の２×２画素の色信号を、同一空間の位相に補間する補間処理、信号のノイズ成分を除去するノイズ除去処理、信号帯域を制限するフィルタリング処理、および信号帯域の高域成分を補正する高域周波数補正処理等を行い、得られたRG1G2B信号を、リニアマトリクス処理部９５に出力する。
【００４２】
リニアマトリクス処理部９５は、所定のリニアマトリクス係数Ｍ（３×４の行列）に基づいて式（１）によるRG1G2B信号の演算を行ない、３色のRGB信号を生成する。
【００４３】
【数１】
【００４４】
リニアマトリクス処理部９５により生成されたＲ信号は、ガンマ補正処理部９６−１に出力され、Ｇ信号は、ガンマ補正処理部９６−２に出力され、Ｂ信号は、ガンマ補正処理部９６−３に出力される。
【００４５】
ガンマ補正処理部９６−１乃至９６−３は、リニアマトリクス処理部９５から出力されてきたRGB信号のそれぞれの信号に対してガンマ補正を行い、取得されたRGB信号を輝度（Ｙ）信号生成処理部９７および色差（Ｃ）信号生成処理部９８に出力する。
【００４６】
輝度信号生成処理部９７は、ガンマ補正処理部９６−１乃至９６−３から供給されてきたRGB信号を、例えば、式（２）に従って、所定の合成比で合成し、輝度信号（Ｙ）を生成する。
Ｙ＝０．２１２６Ｒ＋０．７１５２Ｇ＋０．０７２２Ｂ・・・（２）
色差信号生成処理部９８も同様に、ガンマ補正処理部９６−１乃至９６−３から供給されてきたRGB信号を所定の合成比で合成して色差信号（Ｃ）を生成し、帯域制限間引き処理部９９に出力する。帯域制限間引き処理部９９は、色差信号（Ｃ）に基づき色差信号（Ｃｂ，Ｃｒ）を生成する。なお、単板２×２色コーディングによって得られた信号には、一般的に色情報の帯域が輝度信号ほど存在しない。したがって、帯域制限間引き処理部９９は、色差信号生成処理部９８から供給された色差信号（Ｃ）を帯域制限処理と間引き処理を行なうことで色情報データを削減することで色差信号（Ｃｂ，Ｃｒ）を生成する。
【００４７】
輝度信号生成処理部９７により生成された輝度信号（Ｙ）、および色差信号生成処理部９８により生成された色差信号（Ｃ）又は帯域制限間引き処理部９９により生成された色差信号（Ｃｂ，Ｃｒ）は、例えば、図９に示すモニタインタフェース７７を介して画像モニタ５０に出力される。
【００４８】
以上のような構成を有する撮像装置において、画像の撮影が指示された場合、マイクロコンピュータ４１は、TG４６を制御し、画像センサ４５により画像を取り込ませる。すなわち、画像センサ４５を構成するCCD等の撮像素子（以下、CCD撮像素子という。）の前面に形成されている４色カラーフィルタ６１により４色の光が透過され、透過された光がCCD撮像素子により取り込まれる。CCD撮像素子により取り込まれた光は、４色の色信号に変換され、それがフロントエンド４７に出力される。
【００４９】
フロントエンド４７は、画像センサ４５から供給されてきた色信号に対して、ノイズ成分を除去するための相関二重サンプリング処理、ゲインコントロール処理、およびディジタル変換処理等を施し、得られた画像データをカメラシステムLSI４８に出力する。
【００５０】
カメラシステムLSI４８の信号処理部７１においては、オフセット補正処理部９１により色信号のオフセット成分が除去され、ホワイトバランス補正処理部９２により、画像信号の色温度、および４色カラーフィルタ６１の各フィルタの感度の違いに基づいて、各色のバランスが補正される。
【００５１】
また、垂直方向同時化処理部９３により、ホワイトバランス補正処理部９２により補正された信号の垂直方向の時間のずれが同時化（補正）され、信号生成処理部９４により、垂直方向同時化処理部９３から供給されてきたRG1G2B信号の最小単位の２×２画素の色信号を、同一空間の位相に補間する補間処理、信号のノイズ成分を除去するノイズ除去処理、信号帯域を制限するフィルタリング処理、および信号帯域の高域成分を補正する高域周波数補正処理等が行われる。
【００５２】
さらに、リニアマトリクス処理部９５においては、信号生成処理部９４により生成された信号（RG1G2B信号）が、所定のリニアマトリクス係数Ｍ（３×４の行列）に基づいて変換され、３色のRGB信号が生成される。リニアマトリクス処理部９５により生成されたＲ信号は、ガンマ補正処理部９６−１に出力され、Ｇ信号は、ガンマ補正処理部９６−２に出力され、Ｂ信号は、ガンマ補正処理部９６−３に出力される。
【００５３】
ガンマ補正処理部９６−１乃至９６−３により、リニアマトリクス処理部９５の処理により得られたRGB信号のそれぞれの信号に対してガンマ補正が行われ、取得されたRGB信号が輝度信号生成処理部９７および色差信号生成処理部９８に出力される。輝度信号生成処理部９７、および色差信号生成処理部９８においては、ガンマ補正処理部９６−１乃至９６−３から供給されてきたＲ信号、Ｇ信号、Ｂ信号のそれぞれの信号が所定の合成比で合成され、輝度信号（Ｙ）および色差信号（Ｃ）が生成される。輝度信号生成処理部９７により生成された輝度信号（Ｙ）、および色差信号生成処理部９８により生成された色差信号（Ｃ）は、図９の画像圧縮・解凍部７６に出力され、例えば、JPEGフォーマットで圧縮される。圧縮され、得られた画像データは、マイクロコンピュータインタフェース７３を介して外部記憶媒体５１に出力され、記憶される。
【００５４】
以上のように、１つの画像データが４種類の色信号に基づいて形成されるため、その再現性は、人間の眼の見えにより近いものとなる。
【００５５】
一方、外部記憶媒体５１に記憶されている画像データの再生（表示）が指示されたとき、マイクロコンピュータ４１により、外部記憶媒体５１に記憶されている画像データが読み出され、それがカメラシステムLSI４８の画像圧縮・解凍部７６に出力される。画像圧縮・解凍部７６においては、圧縮されている画像データが伸張され、モニタインタフェース７７を介して、得られたデータに対応する画像が画像モニタ５０に表示される。
【００５６】
次に、図１１に示すフローチャートを参照して、以上のような構成を有する撮像装置を作成する処理（手順）について説明する。
【００５７】
ステップＳ１において、図５に示す画像センサ４５に設けられる４色カラーフィルタ６１の分光感度特性を決定する４色カラーフィルタ決定処理が行われ、ステップＳ２において、図１０に示すリニアマトリクス処理部９５に設定されるマトリクス係数Ｍを決定するリニアマトリクス係数Ｍ決定処理が行われる。ステップＳ１において実行される４色カラーフィルタ決定処理の詳細については、図１２に示すフローチャートを参照して、また、ステップＳ２において実行されるリニアマトリクス係数Ｍ決定処理の詳細については、図１６に示すフローチャートを参照して、それぞれ後述する。
【００５８】
４色カラーフィルタ６１が決定され、マトリクス係数が決定された後、ステップＳ３において、図１０に示す信号処理部７１が作成され、ステップＳ４に進み、図９に示すカメラシステムLSI４８が作成される。また、ステップＳ５において、図５に示すような撮像装置（例えば、ディジタルカメラ）の全体が作成される。ステップＳ６において、ステップＳ５で作成された撮像装置の画質（「色再現性」、「色判別性」）の評価が行われ、処理が終了される。
【００５９】
ここで、「色再現性」、「色判別性」などを評価する際に参照される物体色について説明する。物体色は、「物体の分光反射率」、「標準照明の分光エネルギー分布」、および「物体を感知するセンサ（カラーフィルタ）の分光感度分布（特性）」の積を可視光領域（例えば、４００乃至７００nm）の範囲で積分した値によって算出される。すなわち、式（３）により物体色が算出される。
【００６０】
【数２】
【００６１】
例えば、所定の物体を眼で観察する場合、式（３）の「センサの分光感度特性」が等色関数で表され、その物体の物体色は、Ｘ，Ｙ，Ｚの三刺激値で表される。具体的には、Ｘの値は式（４−１）で算出され、Ｙの値は式（４−２）で算出され、Ｚの値は式（４−３）で算出される。なお、式（４−１）乃至（４−３）における定数ｋの値は式（４−４）で算出される。
【００６２】
【数３】
【００６３】
また、ディジタルカメラなどの撮像装置により、所定の物体の画像を取り込む場合、式（３）の「センサの分光感度特性」がカラーフィルタの分光感度特性で表され、その物体の物体色は、フィルタの数の色値（例えば、RGBフィルタ（３種類）の場合はRGB値（３値））の物体色が算出される。撮像装置に、３種類の色を検出するRGBフィルタが設けられている場合、具体的には、Ｒの値は式（５−１）で算出され、Ｇの値は式（５−２）で算出され、Ｂの値は式（５−３）で算出される。また、式（５−１）における定数ｋrの値は式（５−４）で算出され、式（５−２）における定数ｋgの値は式（５−５）で算出され、式（５−３）における定数ｋbの値は式（５−６）で算出される。
【００６４】
【数４】
【００６５】
次に、図１２に示すフローチャートを参照して、図１１に示すステップＳ１で行われる４色カラーフィルタ決定処理について説明する。
【００６６】
なお、４色カラーフィルタの決定方法としては様々な方法があるが、例えば、RGBフィルタを基調として（既存の（図１）Ｇフィルタの一方をＧ１フィルタとして）、Ｇ１フィルタを透過する色と相関の高い色を透過するＧ２フィルタを選択し、それを追加して４色カラーフィルタを決定する処理について説明する。
【００６７】
ステップＳ２１において、UMG値を算出するために使用されるカラーターゲットが選択される。例えば、ステップＳ２１において、現存する色を代表するカラーパッチを多く含み、かつ、人間の記憶色（肌色、植物の緑、空の青等）を重視したカラーパッチを多く含むカラーターゲットが選択される。カラーターゲットとしては、例えば、IT8.7，Macbeth Color Checker，GretagMacbeth DigitalCamera Color Checker，CIE，Color Bar等がある。
【００６８】
また、目的に応じて、SOCS(Standard Object Color Spectra Database)などのデータから、標準となりうるカラーパッチを作成し、それを使用するようにしてもよい。なお、SOCSについては「田島譲二, "標準物体色分光データベース(SOCS)による統計的色再現評価", カラーフォーラムJAPAN 99」にその詳細が開示されている。以下、Macbeth Color Checkerがカラーターゲットとして選択された場合について説明する。
【００６９】
ステップＳ２２において、Ｇ２フィルタの分光感度特性が決定される。分光感度特性として、実存する材料から作成可能なものを使用するようにしてもよいし、図１３に示すようなcubic spline曲線（３次スプライン関数）で仮想曲線Ｃ（λ）を想定し、仮想曲線Ｃ（λ）のピーク値λ0、値ｗ（ｗ1とｗ2の和を２で除算した値）、値Δｗ（ｗ1からｗ2を減算した値を２で除算した値）を、図中に示す範囲で変化させたものを使用するようにしてもよい。なお、ｗ，Δｗの値は、半値幅の値に基づく値とされる。λ0，ｗ，Δｗの変化の方法は、例えば、５nm刻みとする。仮想曲線Ｃ（λ）は、それぞれの範囲において、下式（６−１）乃至（６−５）で表される。
【００７０】
【数５】
【００７１】
なお、この例では、フィルタＧ２のみが追加されるが、図１に示すフィルタ（Ｒ，Ｇ，Ｇ，Ｂ）のＲフィルタとＢフィルタのみを使用し、残りの２つのＧ１，Ｇ２フィルタを、緑色付近の上式（６−１）乃至（６−５）の仮想曲線として定義することも可能である。また、同様に、ＲとＧのみ、ＧとＢのみを、図１に示すフィルタから使用するようにしてもよい。さらに、４色のフィルタのうち、３色を仮想曲線、または４色とも仮想曲線として定義することも可能である。
【００７２】
ステップＳ２３において、追加するフィルタ（Ｇ２フィルタ）と、現存のフィルタ（Ｒフィルタ，Ｇ１フィルタ，Ｂフィルタ）が組み合わされ、４色カラーフィルタの最小単位（セット）が作成される。また、ステップＳ２４において、ステップＳ２３で作成された４色カラーフィルタに対して、フィルタ評価係数としてUMGが用いられ、UMG値が算出される。
【００７３】
図８を参照して説明したように、UMGを用いた場合、４色のそれぞれのカラーフィルタに対して一度で評価を行うことができる。また、物体の分光反射率を考慮して評価が行われるだけでなく、ノイズの低減性をも考慮して評価が行われる。UMGを用いた評価においては、それぞれのフィルタの分光感度特性に、適度な重なりがあるフィルタに対して高い評価が示されるため、例えば、Ｒの特性とＧの特性が広い波長帯域にわたって重複する特性を有するフィルタ（それぞれの色信号を分離した際、ノイズが増幅されるフィルタ）に対して、高い評価が示されるということを抑制することができる。
【００７４】
図１４は、３色カラーフィルタにおいて算出されるUMG値の例を示す図である。例えば、RGBの特性がそれぞれ重複しない図１４Ａに示すような特性のフィルタにおいては、「０．７９４２」のUMG値が算出され、Ｒの特性とＧの特性が広い波長帯域にわたって重複する図１４Ｂに示すような特性のフィルタにおいては、「０．８２１１」のUMG値が算出される。また、RGBのそれぞれの特性が適度に重なる図１４Ｃに示すような特性を有するフィルタにおいては、「０．８８７９」のUMG値が算出される。すなわち、RGBのそれぞれの特性が適度に重なる、図１４Ｃに示すような特性を有するフィルタに対して最も高い評価が示される。このことは、４色カラーフィルタにおいても同様である。なお、図１４Ａに示す曲線Ｌ３１、図１４Ｂに示す曲線Ｌ４１、および図１４Ｃに示す曲線Ｌ５１はＲの分光感度を表し、図１４Ａに示す曲線Ｌ３２、図１４Ｂに示す曲線Ｌ４２、および図１４Ｃに示す曲線Ｌ５２はＧの分光感度を表し、図１４Ａに示す曲線Ｌ３３、図１４Ｂに示す曲線Ｌ４３、および図１４Ｃに示す曲線Ｌ５３はＢの分光感度を表している。
【００７５】
ステップＳ２５において、ステップＳ２４で算出されたUMG値が、所定の閾値である「０．９５」以上であるか否かが判定され、「０．９５」未満であると判定された場合、ステップＳ２６に進み、作成された４色カラーフィルタが却下される（使用されない）。ステップＳ２６において、４色カラーフィルタが却下された場合、その後、処理は終了される（図１１に示すステップＳ２以降の処理は実行されない）。
【００７６】
一方、ステップＳ２５において、ステップＳ２４で算出されたUMG値が「０．９５」以上であると判定された場合、ステップＳ２７において、その４色カラーフィルタが、ディジタルカメラで使用される候補のフィルタとされる。
【００７７】
ステップＳ２８において、ステップＳ２７で候補のフィルタとされた４色カラーフィルタが、現存する材料、染料で実現可能であるか否かが判定される。材料、染料等が取得困難である場合、実現不可能であると判定され、ステップＳ２６に進み、その４色カラーフィルタが却下される。
【００７８】
一方、ステップＳ２８において、材料、染料等が取得可能であり、実現可能であると判定された場合、ステップＳ２９に進み、作成された４色カラーフィルタが、ディジタルカメラにおいて使用されるフィルタとして決定される。その後、図１１に示すステップＳ２以降の処理が実行される。
【００７９】
図１５は、ステップＳ２９において決定された４色カラーフィルタの分光感度特性の例を示す図である。
【００８０】
図１５において、曲線Ｌ６１はＲの分光感度を表し、曲線Ｌ６２はＧ１の分光感度を表している。また、曲線Ｌ６３はＧ２の分光感度を表し、曲線Ｌ６４はＢの分光感度を表している。図１５に示すように、Ｇ２の分光感度曲線（曲線Ｌ６３）は、Ｇ１の分光感度曲線（曲線Ｌ６２）に対して相関が高い。また、Ｒの分光感度、Ｇ（Ｇ１，Ｇ２）の分光感度、およびＢの分光感度は、それぞれ適度な範囲において重複している。
【００８１】
以上のようにして決定された４色カラーフィルタを利用することにより、特に、「色再現性」のうちの「色の判別性」を改善させることができる。
【００８２】
なお、光の利用効率の観点から、以上のように、現存するRGBフィルタのGフィルタと相関の高いフィルタを、追加するフィルタ（Ｇ２フィルタ）とすることが好ましい。この場合、追加するフィルタの分光感度曲線のピーク値は、経験的に４９５乃至５３５nmの範囲（現存するＧフィルタの分光感度曲線のピーク値の近傍）にあることが望ましい。
【００８３】
また、現存するＧフィルタと相関の高いフィルタを追加する場合、図１に示す最小単位（Ｒ，Ｇ，Ｇ，Ｂ）を構成する２つのＧフィルタのいずれか一方を追加色のフィルタにするだけで４色カラーフィルタを作成することができるため、作成の行程に大きな変更を加える必要がない。
【００８４】
以上のようにして４色カラーフィルタが作成され、それがディジタルカメラに設けられた場合、図１０に示す信号処理装置７１には４種類の色信号が信号生成処理部９４から供給されるため、リニアマトリクス処理部９５において、４色（Ｒ，Ｇ１，Ｇ２，Ｂ）の信号から３色（Ｒ，Ｇ，Ｂ）の信号を生成する変換処理が行われる。この変換処理は、輝度リニアな（輝度値を線形的な変換により表すことができる）入力信号値に対するマトリクス処理であるため、以下、リニアマトリクス処理部９５において行われる変換処理を、適宜、リニアマトリクス処理と称する。
【００８５】
次に、図１６に示すフローチャートを参照して、図１１に示すステップＳ２において実行される、リニアマトリクス係数Ｍ決定処理について説明する。なお、リニアマトリクス係数Ｍ決定処理において使用されるカラーターゲットはMacbeth Color Checkerとし、また、使用される４色カラーフィルタは図１５に示す分光感度特性を有するものとする。
【００８６】
ステップＳ４１において、例えば、CIE(Commision Internationale del'Eclairange)において標準光源とされている一般的な昼光Ｄ６５（照明光Ｌ（λ））が照明光として選択される。なお、照明光は、画像処理装置が頻繁に使用されることが予想される環境の照明光等に変更してもよい。また、想定される照明環境が複数ある場合、リニアマトリクスを複数用意することも考えられる。以下、照明光として昼光Ｄ６５が選択された場合について説明する。
【００８７】
ステップＳ４２において、リファレンス値（参照値）Ｘｒ，Ｙｒ，Ｚｒが算出される。具体的には、リファレンス値Ｘｒは式（７−１）により算出され、Ｙｒは式（７−２）により算出され、Ｚｒは式（７−３）により算出される。
【００８８】
【数６】
【００８９】
また、定数ｋは式（８）により算出される。
【００９０】
【数７】
【００９１】
例えば、カラーターゲットがMacbeth Color Checkerの場合、２４色分のリファレンス値が算出される。
【００９２】
次に、ステップＳ４３において、４色カラーフィルタの出力値Ｒf，Ｇ１f，Ｇ２f，Ｂfが算出される。具体的には、Ｒfは式（９−１）により算出され、Ｇ１fは式（９−２）により算出され、Ｇ２fは式（９−３）により算出され、Ｂfは式（９−４）により算出される。
【００９３】
【数８】
【００９４】
また、定数ｋrは式（１０−１）により算出され、定数ｋg1は式（１０−２）により算出され、定数ｋg2は式（１０−３）により算出され、定数ｋbは式（１０−４）により算出される
【００９５】
【数９】
【００９６】
例えば、カラーターゲットがMacbeth Color Checkerの場合は、２４色分の出力値Ｒf，Ｇ１f，Ｇ２f，Ｂfが算出される。
【００９７】
ステップＳ４４において、ステップＳ４３で算出されたフィルタ出力値を、ステップＳ４２で算出されたリファレンス値（XYZref）に近似させる変換を行うマトリクスが、例えば、XYZ色空間における誤差最小二乗法により算出される。
【００９８】
例えば、算出される、３×４のマトリクスを式（１１）で表されるＡとした場合、マトリクス変換（XYZexp）は、次の式（１２）で表される。
【００９９】
【数１０】
【０１００】
また、リファレンス値に対するマトリクス変換（式（１２））の誤差の２乗（Ｅ2）は、次の式（１３）で表され、これに基づいてリファレンス値に対するマトリクス変換の誤差を最小にするマトリクスＡが算出される。
【０１０１】
【数１１】
【０１０２】
また、誤差最小二乗法で使用する色空間を、XYZ色空間以外のものに変更するようにしてもよい。例えば、人間の知覚に対して均等なLab，Luv，Lch色空間（知覚均等色空間）に変換した後に同様の演算を行うことにより、知覚的な誤差の少ない色の再現を可能にするリニアマトリクスを算出することができる。なお、これらの色空間の値は、XYZ値から非線形な変換によって算出されるため、誤差最小二乗法においても非線形な計算アルゴリズムが使用される。
【０１０３】
上述したような演算により、例えば、図１５に示す分光感度特性を有するフィルタに対するマトリクス係数として式（１４）で表されるものが算出される。
【０１０４】
【数１２】
【０１０５】
ステップＳ４５において、リニアマトリクスが決定される。例えば、作成される最終のRGB画像データが次の式（１５）で表されるとした場合、リニアマトリクス（LinearＭ）は、以下のようにして算出される。
【０１０６】
RGBout＝［Ｒ０,Ｇ０,Ｂ０］ｔ・・・(15)
すなわち、照明光がＤ６５である場合、sRGB色空間をXYZ色空間に変換する変換式はITU-R709.BTマトリクスを含む式（１６）で表され、そのITU-R709.BTマトリクスの逆マトリクスにより式（１７）が算出される。
【０１０７】
【数１３】
【０１０８】
式（１２）のマトリクス変換式、式（１５）および式（１７）のITU-R709.BTマトリクスの逆マトリクスにより、式（１８）が算出される。式（１８）の右辺には、ITU-R709.BTマトリクスの逆マトリクスと、上述したマトリクスＡを乗算した値としてのリニアマトリクスが含まれる。
【０１０９】
【数１４】
【０１１０】
すなわち、３×４のリニアマトリクス（LinearＭ）は式（１９−１）により表され、例えば、式（１４）のマトリクス係数が用いられる、図１５に示す分光分布特性を有する４色カラーフィルタに対するリニアマトリクスは、式（１９−２）により表される。
【０１１１】
【数１５】
【０１１２】
以上のようにして算出されたリニアマトリクスが、図１０に示すリニアマトリクス処理部９５に与えられる。これにより、輝度を線形変換により表すことができる信号（Ｒ，Ｇ１，Ｇ２，Ｂ）に対してマトリクス処理を行うことができるので、図２に示す信号処理部１１における処理のように、ガンマ処理を施した後に得られる信号に対してマトリクス処理を行う場合に較べて、色彩工学的に、より忠実な色を再現させることができる。
【０１１３】
次に、図１１に示すステップＳ６において行われる評価について説明する。
【０１１４】
以上のようにして作成された、例えば、図１５に示す分光感度特性を有する４色カラーフィルタが設けられた撮像装置の色再現性と、図１に示す３色カラーフィルタが設けられる画像処理装置の色再現性を比較した場合、以下のような差が現れる。
【０１１５】
例えば、マクベスチャートを２種類の画像入力装置（４色カラーフィルタが設けられる撮像装置と、３色カラーフィルタが設けられる撮像装置）で撮像したときの出力値とリファレンス値とのLab色空間における色差が次の式（２０）によりそれぞれ算出される。
【０１１６】
【数１６】
【０１１７】
Ｌ1−Ｌ2は２つの試料の明度差であり、ａ1−ａ2、ｂ1−ｂ2は２つの試料の色相・彩度の成分差を表している。
【０１１８】
図１７は、式（２０）による算出結果を示す図である。図１７に示すように、３色カラーフィルタが設けられる撮像装置の場合は色差が「３．３２」であるのに対して、４色カラーフィルタが設けられる撮像装置の場合は「１．３９」であり、「色の見え方」は４色カラーフィルタが設けられた撮像装置の方が優れている（色差が小さい）。
【０１１９】
図１８においては、物体Ｒ１のＲ値が「４９．４」、Ｇ値が「６４．１」、Ｂ値が「１４９．５」とされ、物体Ｒ２のＲ値が「６６．０」、Ｇ値が「６３．７」、Ｂ値が「１５５．６」とされている。従って、４色カラーフィルタにおいては、物体Ｒ１と物体Ｒ２のRGB値はそれぞれ異なる値となり、眼で見た場合と同様に、それぞれの物体の色が識別されている。すなわち、４種類の色を識別できるフィルタを設けることにより、「色の判別性」が改善されている。
【０１２０】
以上においては、４色カラーフィルタ６１は、図６に示すように、Ｇ１フィルタの左右にＢフィルタが設けられるとともに、Ｇ２フィルタの左右にＲフィルタが設けられるような配列により構成されるとしたが、図１９に示すような配列により構成されるようにしてもよい。図１９に示す４色カラーフィルタ６１においては、Ｇ１フィルタの左右にＲフィルタが設けられるとともに、Ｇ２フィルタの左右にＢフィルタが設けられている。４色カラーフィルタ６１をこのように構成することによっても、図６に示すものと同様に、「色の判別性」、「色の再現性」および「ノイズの低減性」を向上させることができる。
【０１２１】
ところで、色差（ΔＥ値）を最小にするように、リニアマトリクス係数Ｍを決定した場合、もしも、画像センサ４５の前段部に形成されているカラーフィルタの分光感度が図４に示すように重なっているときには、式（２１）のようにリニアマトリクス係数Ｍの差が大きくなってしまう。
【０１２２】
【数１７】
【０１２３】
撮像素子の出力信号には、微小なノイズ成分が含まれているため、このようなリニアマトリクス係数Ｍを使用して色分離処理を行なうと、微小なノイズ成分も増幅されてしまう。したがって、リニアマトリクス係数Ｍの差が大きくならないように、色再現性よりもノイズ低減性を考慮する必要が生じる。しかし、当該撮像装置により実際に被写体を撮像するときに、撮像するシーンや環境により、色再現性よりもノイズ低減性を重視してリニアマトリクス係数Ｍを決定し、適応的にリニアマトリクス処理を行なった方が画質の向上を図れる場合や、逆に、ノイズ低減性よりも色再現性を重視してリニアマトリクス係数Ｍを決定し、適応的にリニアマトリクス処理を行なった方が画質の向上を図れる場合がある。また、撮像装置の使用の用途は、ユーザごとに異なるため、リニアマトリクス係数Ｍの決定をユーザが任意に行ないたい場合がある。
【０１２４】
そこで、本願発明に係る撮像装置では、上述したような問題点を解決するために、図２０に示すフローチャートにしたがってリニアマトリクス係数Ｍを決定する。
【０１２５】
まず、使用するチャート及び照明光の決定を行ない（ステップＳ５０）、つぎに色再現性指標ΔＥ（Ｍ）の定義を行ない（ステップＳ５１）、つぎにノイズ低減性指標σＮ（Ｍ）の定義を行ない（ステップＳ５２）、つぎに評価指数ＥＥＶ（Ｍ）（Error Evaluation Value）の定義を行ない（ステップＳ５３）、そして、評価指数ＥＥＶ（Ｍ）に基づきリニアマトリクス係数Ｍを決定する（ステップＳ５４）。なお、ステップ５では、評価指数ＥＥＶ（Ｍ）の係数を撮像条件等により適応的に変化させ、対応するリニアマトリクス係数Ｍを決定する。以下に、各ステップの詳細について説明する。
【０１２６】
使用するチャート及び照明光の決定（ステップＳ５０）について説明する。リニアマトリクス係数Ｍを決定するためには、カラーチャートと、当該カラーチャートを照明する光源を決める必要がある。カラーチャートとしては、Macbeth Color Checker, Digital Camera Color Checker, IT8.7等の複数の均一な色平面を持つカラーパッチから成る様々な反射チャート又は透過チャートが考えられる。照明光としては、撮像装置がよく使用される環境の光に近い分光感度を持つもの（例えばＤ５５や光源等）が考えられる。なお、撮像装置は、ユーザの使用目的により様々な光源下で用いられることが考えられるので、照明光としては、よく使用される環境の光のみに限定されるものではない。
【０１２７】
つぎに、色再現性指標ΔＥ（Ｍ）の定義（ステップＳ５１）について説明する。色再現性は、目標色と撮像装置のリニアマトリクス処理部９５でリニアマトリクス処理が行なわれた信号値が示す色（以下、出力色という。）との差で定義される。なお、色の値としては、RGB値、YCｂCr値又はXYZ値等様々考えられるが、人間の眼の見えに対して知覚が均等な色空間の値（L*a*b*値、L*u*v*値等）を用いて定義する方が、より色の差を正確に示すことが可能となる。例えば、カラーチャート中のｋ番目のカラーパッチの目標色をLabref_k(L*ref_k, a*ref_k, b*ref_k)とし、撮像装置の出力色をL*a*b*shot_k(L*shot_k, a*shot_k, b*shot_k)とすると、このパッチの色差ΔＥｋは式（２２）のようになる。
【０１２８】
【数１８】
【０１２９】
また、色再現性指標ΔＥ（Ｍ）としては、カラーチャートの各パッチの平均ΔＥ値や、各パッチに重み付けを行ってある特定色の色再現性を重視した値等が考えられる。
【０１３０】
【数１９】
【０１３１】
ただし、ｗｋは、各パッチに対する重み付け係数を示し、TotalPatchNumは、カラーパッチの総数を示す。
【０１３２】
また、実際には撮像装置の出力色には、リニアマトリクス処理が施されているので、L*a*b*shot_kは、リニアマトリクス係数Ｍの関数値となり、したがってΔＥk、ΔＥもＭの関数値となる。
【０１３３】
【数２０】
【０１３４】
つぎに、ノイズ低減性指標σＮ（Ｍ）の定義（ステップＳ５２）について説明する。ノイズ低減性指標σＮ（Ｍ）は、撮像装置のリニアマトリクス処理部９５でリニアマトリクス処理が行なわれた信号値の標準偏差で定義する。信号値としては、RGB値、YCｂCr値又はXYZ値等様々考えられるが、人間の眼の見えに対して知覚が均等な色空間の値（L*a*b*値、L*u*v*値）等を用いて信号値を定義する方が、より人が感じるノイズ感と相関のあるノイズ値σＮｋを得ることが可能となる。なお、ノイズ値は、信号値の色空間に対応して、色空間の各成分の標準偏差となる。例えば、RGB空間の場合には、ノイズ値は、σR,σG,σBとなり、XYZ空間の場合には、ノイズ値は、σX,σY,σZとなる。ノイズ低減性指標σＮ（Ｍ）の定義では、これらのノイズ値を使って１つのノイズ指標を決定する。例えば、あるカラーパッチを照明光下で撮像した場合、L*a*b*空間のノイズ値σＮｋは、明度ノイズをσL*k、色ノイズをσa*k,σb*kとし、これらの値から例えば式（２５）のような明度及び色ノイズを考慮したものとして定義される。
【０１３５】
【数２１】
【０１３６】
ただし、ｗＬ＊ ｋ、ｗａ＊ ｋ、ｗｂ＊ ｋは、各標準偏差値への重み付け係数を示し、人間の眼が感じるノイズ感との相関で適切に設定されるものである。なお、ノイズ値σＮｋとしては、他の色空間の分散値を使ったもの等様々考えられる。ノイズ低減性指標σＮ（Ｍ）としては、カラーチャートの各パッチの平均σＮ値や、各パッチに重み付けを行ってある特定色のノイズ低減性を重視した値等が考えられる。
【０１３７】
【数２２】
【０１３８】
実際には、撮像装置の信号値には、リニアマトリクス処理が施されているので、σＮｋ及びσＮはリニアマトリクス係数Ｍの関数値となる。
【０１３９】
【数２３】
【０１４０】
つぎに、評価指数ＥＥＶ（Ｍ）の定義（ステップＳ５３）について説明する。上述したステップＳ５１及びステップＳ５２の定義により、リニアマトリクス係数Ｍの関数値である色再現性指標ΔＥ（Ｍ）及びノイズ低減性指標σＮ（Ｍ）の二値を考慮した評価指数ＥＥＶ(Error Evaluation Value)（Ｍ）を式（２８）の様に定義する。
【０１４１】
【数２４】
【０１４２】
ただし、h, I, j, k, l, は、関数を示し、ｗｃは、色差に対する重み付け係数を示し, ｗｎは、ノイズ値に対する重み付け係数を示す。ｗｃとｗｎを変化させて、評価指数ＥＥＶ（Ｍ）が最小となるようにリニアマトリクス係数Ｍを決定することにより、色再現性及びノイズ低減性の双方を考慮したリニアマトリクス係数Ｍの決定が可能となる。なお、色再現性を重視したい場合は、ｗｃ>ｗｎとし、ノイズ低減性を重視したい場合は、ｗｃ < ｗｎと重み付けを設定すれば良い。
【０１４３】
つぎに、リニアマトリクス係数Ｍを決定する手段（ステップＳ５４）について説明する。ステップＳ５３により定義した評価指数ＥＥＶ（Ｍ）に誤差最小二乗法を適用してリニアマトリクス係数Ｍを決定する。なお、ステップＳ５４では、ｗｃ及びｗｎを適当に定めて、例えばNewton法、Steepest Descent法又はConjugate Gradient法等を回帰アルゴリズムとして使用して誤差最小二乗法を適用し、リニアマトリクス係数Ｍを決定する。
【０１４４】
また、ステップＳ５４では、撮像装置により被写体を撮像するときの環境や条件等によりステップＳ５３で定義した評価指数ＥＥＶ（Ｍ）の色差に対する重み付け係数ｗｃ及びノイズ値に対する重み付け係数 ｗｎを適応的に変化させ、誤差最小二乗法によりリニアマトリクス係数Ｍを決定する。図２１に、色再現性指標ΔＥ（Ｍ）を変化させたときのノイズ値低減性指標σＮ（Ｍ）の変化の様子を示す。図２１に示すように、一つの撮像素子においても、リニアマトリクス係数Ｍによって色再現性指標ΔＥ（Ｍ）とノイズ低減性指標σＮ（Ｍ）のトレードオフが存在する。この結果を用いて、様々な撮像環境及び条件等に応じて適応的にリニアマトリクス係数Ｍを決定する。また、あらかじめいくつかのリニアマトリクス係数Ｍのセットを用意し、ユーザが必要に応じてリニアマトリクス係数Ｍを選択して色再現性指標ΔＥ（Ｍ）及びノイズ低減性指標σＮ（Ｍ）を調節することも可能である。
【０１４５】
ここで、撮像装置が図２２に示すような特性の４色カラーフィルターよりなるCCD撮像素子を有している場合に、上述したステップＳ５０〜ステップＳ５４にしたがってリニアマトリクス係数Ｍを決定する具体例について述べる。
【０１４６】
まず、使用するチャート及び照明光を決定する（ステップＳ５０）。カラーチャートは、Macbeth Color Checker（２４色のカラーパッチを含む）を使用し、照明光は、Ｄ５５光源（ＣＩＥで定義される５５００ｋの標準昼光）を使用する。なお、それぞれの分光データは、例えば、分光放射輝度計を用いて測定されているものとする。
【０１４７】
つぎに、色再現性指標ΔＥ（Ｍ）の定義を行なう（ステップＳ５１）。目標色を人間の眼の見えとし、Lab空間での色差ΔＥを指標とする。一般的に物体色は、「物体の分光反射率」と「照明の分光エネルギー分布」と「物体を感知するセンサの分光感度分布」との積を可視光領域vis（通常400nm〜700nm）の範囲で積分した値によって式（２９）のように定義される。
物体色＝∫vis（物体の分光反射率）・（照明の分光輝度）・（物体を感知するセンサの分光感度）・・・（２９）
さらに、人間の眼の分光感度を表す等色関数を用いると、カラーチャートのｋ番目のパッチの目標色XYZref_k(Xref_k, Yref_k, Zref_k)は、式（２９）を用いて式（３０）のように表すことができる。
【０１４８】
【数２５】
【０１４９】
また、通常は、式（３１）を使用してXYZ空間の色をL*a*b*空間の色に変換する。
【０１５０】
【数２６】
【０１５１】
ただし、（Ｘｎ、Ｙｎ、Ｚｎ）は、完全拡散反射面（白点）のＸＹＺ値を示している。
【０１５２】
そして、目標色XYZref_kを式（３１）を用いてL*a*b*ref_k(L*ref_k, a*ref_k, b*ref_k)に変換する。
【０１５３】
また、CCD撮像素子から出力される信号値であるｒａｗデータRGBXraw_k(Rraw_k, Graw_k, Braw_k, Xraw_k)は、式（２９）を用いると式（３２）のようになる。
【０１５４】
【数２７】
【０１５５】
撮像装置は、リニアマトリクス処理部９５でｒａｗデータRGBXraw_k(Rraw_k, Graw_k, Braw_k, Xraw_k)をリニアマトリクス係数Ｍ(ｍ０〜ｍ１１)を用いてリニアマトリクス処理をするので、リニアマトリクス処理後の撮像データは式（３３）のようになる。
【０１５６】
【数２８】
【０１５７】
そして、上記RGBcam_kをXYZ値に変換する。なお、変換には、式（３４）に示すように、一般的に使用されている７０９系マトリクスＭ７０９を使用する。
【０１５８】
【数２９】
【０１５９】
つぎに、式（３４）を用いて、Xcam_k、Ycam_k、及びZcam_kを求める。
【０１６０】
【数３０】
【０１６１】
また、式（３１）を用いてL*a*b*値(L*a*b*cam_k(L*cam_k, a*cam_k, b*cam_k))に変換し、カラーチャートのｋ番目のパッチの色差ΔEkを式（３６）のように定義する。
【０１６２】
【数３１】
【０１６３】
なお、ΔEkは、L*a*b*cam_kの値を使用しているので、リニアマトリクス係数Ｍの関数値でもあるので、ΔEk(M)と表すことができる。色再現性指標ΔE(M)は、式（３７）に示すように、各カラーパッチの色差の平均値として定義する。
【０１６４】
【数３２】
【０１６５】
つぎに、ノイズ低減性指標σＮ（Ｍ）の定義を行なう（ステップＳ５２）。なお、ここでは、撮像装置のリニアマトリクス処理部９５によるリニアマトリクス処理後の信号値に含まれているσＬ（Ｍ）成分に基づきノイズ低減性指標σＮ（Ｍ）を定義する。
【０１６６】
一般的にCCD撮像素子自体が出力する信号CV_CCDに含まれているノイズNoise_rawは、式（３８）のように定義される。
【０１６７】
【数３３】
【０１６８】
なお、ShotNoiseCoefとDarkNoiseは、CCD撮像素子のデバイス特性によって決まる値である。DarkNoiseは、信号値に依存しないノイズ成分(Fixed Pattern Noise等)を表し、ShotNoiseは、信号値に依存するノイズ成分（Sensor Dark Noise, Photon Shot Noise等）を表す。
【０１６９】
式（３１）を使用すると、評価する撮像装置のｋ番目のカラーパッチのｒａｗデータに含まれているノイズ成分は式（３９）のように定義される。
【０１７０】
【数３４】
【０１７１】
また、文献（P.D.Burns and R.S.Berns, "Error Propagation Analysis in Color Measurement and Imaging", Color Research and Application, 1997）には、以下に示すようなノイズ伝搬理論が記載されている。
【０１７２】
所定の入力信号Ｘ'＝［ｘ１，ｘ２，・・・，ｘｎ］が、（ｍ×ｎ）マトリクスＡによって、Ｙ'＝［ｙ１，ｙ２，・・・，ｙｍ］に線形変換されるとすると、式（４０）のように表される。
Ｙ'＝Ａ・Ｘ'・・・（４０）
いま、入力信号Ｘ'の分散共分散行列Σｘが、式（４１）に示すようになると、対角成分が入力信号のNoise分散値となる。
【０１７３】
【数３５】
【０１７４】
もし、入力信号値同士の相関がなければ、行列成分中の共分散成分（すなわち非対角成分）は０となる。このとき、出力信号Ｙの分散共分散行列Σｙは、式（４２）のように定義される。
Σｙ＝Ａ・ΣｘＡ'・・・（４２）
なお、式（４２）は、線形変換で変換することができる色空間同士での、Noise分散値の伝搬理論式となる。
【０１７５】
また、撮像装置のリニアマトリクス処理部９５によるリニアマトリクス処理後の信号RGBcam_kをL*a*b*cam_kに変換するためには、RGB空間からXYZ空間への変換処理（以下、RGB→XYZ変換処理という。）を行ない、その後、XYZ空間からL*a*b*空間への変換処理（以下、XYZ→L*a*b*変換処理という。）を行なう必要がある。RGB→XYZ変換処理は、式（３４）に示す７０９系マトリクスＭ７０９を使用して線形変換を行なうことができるが、XYZ→L*a*b*変換処理は、式（３１）に示すように非線形変換を行なう必要がある。しかし、XYZ→L*a*b*変換処理は、ノイズ量が微小量であることから、JacobianマトリクスJL*a*b*_kを使用して線形近似することができ、RGB→XYZ変換処理と同様に線形変換により行なうことができる。
【０１７６】
リニアマトリクス後の信号値をXYZ値に変換した値をXYZcam_k(Xcam_k,Ycam_k,Zcam_k)とすると、式（４３）のように表すことができる。
JL*a*b*_k＝J0D(XYZcam_k)・・・（４３）
ただし、Ｊ０は、式（４４）とする。
【０１７７】
【数３６】
【０１７８】
したがって、CCD撮像素子から出力されるｒａｗデータをLab値に線形変換するための近似マトリクスMtotal_kは式（４６）のようになる。
【０１７９】
【数３７】
【０１８０】
式（４６）のマトリクスと式（４２）に示したノイズ伝搬理論の式を適用すると、ｋ番目のカラーパッチにおける明度ノイズσLkは、式（４７）で計算することができる。
【０１８１】
【数３８】
【０１８２】
したがって、式（４７）より式（４８）を導き出すことができる。
【０１８３】
【数３９】
【０１８４】
なお、式（４８）は、リニアマトリクス係数Ｍの関数なのでσL*k(M)で表すことができる。ノイズ低減性指標σN(M)は、カラーパッチの各明度ノイズの平均値なので式（４９）のように定義することができる。
【０１８５】
【数４０】
【０１８６】
つぎに、上述のように定義した色再現性指標ΔE(M)及びノイズ低減性指標σN(M)を考慮した評価指数EEV(M)を式（５０）のように定義する（ステップＳ５３）。なお、EEVは、Error Evaluation Valueの略である。
【０１８７】
【数４１】
【０１８８】
つぎに、式（５０）を誤差最小二乗法によって解き、リニアマトリクス係数Ｍを決定する。例えば、ｗｃ＝１としｗｎ＝２とすると評価指数EEV(M)は、式（５１）となる。
【０１８９】
【数４２】
【０１９０】
式（５１）を誤差最小二乗法を用いて解き、リニアマトリクス係数Ｍを式（５２）のように決定する（ステップＳ５４）。
【０１９１】
【数４３】
【０１９２】
一方、ｗｃ＝１としｗｎ＝0とした、すなわち色再現性だけを追い込んだ色差最小マトリクスは式（５３）のようになる。
【０１９３】
【数４４】
【０１９４】
式（５２）と、式（５３）を比較すると、明らかに式（５３）の係数の方が係数間の大小差が大きく、ノイズを増長するマトリクスであることが確認できる。
【０１９５】
ここで、撮像装置のISO感度の設定に基づきリニアマトリクス係数Ｍを適応的に決定する具体例について述べる。撮像装置では、ISO感度の設定に基づいてCCD撮像素子から入力される信号（以下、入力信号という。）を増幅又は減衰させる。撮像装置のISO感度の設定をISO100からISO200に変更したときには、入力信号は、例えば、ISO100のときの２倍に増幅されて撮像装置に入力される。しかし、撮像装置では、ISO感度の設定状態にかかわらず全ての入力信号に対して同一のリニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なうため、ISO感度の設定が高い場合には、入力信号の増幅に伴い、入力信号に含まれているノイズ成分も一緒に増幅されてしまう。したがって、ISO感度の設定を高くして高解像度の画像を得ようとしても、増幅されたノイズ成分が含まれている画像が生成されてしまう。
【０１９６】
そこで、本願発明に係る撮像装置では、ISO感度の設定に基づき増幅又は減衰される入力信号に含まれているノイズ成分の低減性を考慮してリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なう。例えば、表１に示すように、ノイズ低減性に対する重み係数（ｗｎ）をISO感度の設定に応じて変化できるようにし、ｗｃ及びｗｎを式（５０）に代入してISO感度の設定ごとのリニアマトリクス係数Ｍを決定する。したがって、撮像装置は、ISO感度の設定状態に基づいて決定されたリニアマトリクス係数Ｍを用いてリニアマトリクス処理を行なうことができるので、ISO感度の設定を高くしても、それに伴いノイズ成分が増幅することがなく、高解像度の画像を得ることが可能となる。
【０１９７】
【表１】
【０１９８】
ここで、撮像装置により被写体を撮像する環境に基づきリニアマトリクス係数Ｍを適応的に決定する例について述べる。例えば、夜景等の被写体を撮像装置により撮像した場合には、生成された画像の大半を暗い部分（暗部）が占めることがあり、その暗部にノイズが生じていると非常に目立ってしまう。このような場合には、色再現性よりもノイズ成分の低減性を考慮した方が良い。
【０１９９】
そこで、本発明に係る撮像装置では、被写体を撮像するシーンに基づきノイズ低減性と色再現性を考慮してリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なう。例えば、図２３のヒストグラムに示すように、撮像装置の輝度ダイナミックレンジの半分以下の領域に、生成した画像の輝度成分の７０％以上が含まれているようなときには、ノイズ低減性を重視してリニアマトリクス係数Ｍを決定し、それ以外のときには色再現性及びノイズ低減性を考慮してリニアマトリクス係数Ｍを決定する。具体的には、表２に示すように、ノイズ低減性に対する重み係数（ｗｎ）を撮像シーンに応じて変化できるようにし、ｗｃ及びｗｎを式（５０）に代入して、撮像シーンごとにリニアマトリクス係数Ｍを決定する。したがって、撮像装置は、撮像シーンに基づいて決定されたリニアマトリクス係数Ｍを用いてリニアマトリクス処理を行なうことができるので、生成された画像の大半を暗部が占めても、ノイズ成分を目立たなくすることが可能である。
【０２００】
【表２】
【０２０１】
ここで、撮像装置を使用するユーザの要求に基づきリニアマトリクス係数Ｍを適応的に決定する具体例について述べる。撮像装置で被写体を撮影して生成した画像は、ユーザの使用用途によって、色再現性よりもノイズが少ないことが要求されることが多々ある。使用用途は撮像装置の製作会社には分からない事で、ユーザだけが知っている事実である。
【０２０２】
そこで、本願発明に係る撮像装置では、ユーザが意図的する条件に基づきリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なう。例えば、表３に示すように、ノイズ低減性に対する重み係数（ｗｎ）をノイズ量調節変数に応じて変化できるようにし、ｗｎ及びｗｃを式（５０）に代入して、ノイズ量調整変数ごとにリニアマトリクス係数Ｍを決定し、決定したリニアマトリクス係数Ｍを記憶しておき、ユーザがノイズ量調節変数を当該撮像装置のインターフェースを通して変化させたときに、所定のリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なう。したがって、撮像装置は、ユーザの要求に応じて決定されたリニアマトリクス係数Ｍを用いてリニアマトリクス処理を行なうことができるので、ユーザの使用状況に応じたノイズ量調節を行なうことが可能である。
【０２０３】
【表３】
【０２０４】
なお、表３では、ノイズ量調節変数の数値が大きくなればなるほどノイズ量が減ることを意味している。また、撮像前にノイズ量調節変数の設定を行なわず、撮像装置の撮像開始ボタンを押圧したときに、いくつかのリニアマトリクス係数Ｍを使用したリニアマトリクス処理後の画像が複数個メモリに保存されるようにしても良い。
【０２０５】
このように構成された撮像装置は、図１２に示すフローチャートに従って画像センサ４５の前段部に形成される４色カラーフィルタ６１が決定されるので、「色再現性」のうちの「色の判別性」を改善させることができ、また、信号処理部７１で輝度を線形変換により表すことができる信号（Ｒ，Ｇ１，Ｇ２，Ｂ）に対してマトリクス処理を行うので、図２に示す信号処理部１１における処理のように、ガンマ処理を施した後に得られる信号に対してマトリクス処理を行う場合に較べて、色彩工学的に、より忠実な色を再現することができ、さらに、信号処理部７１のリニアマトリクス処理部９５で行なうリニアマトリクス処理において、リニアマトリクス係数Ｍの決定を撮像条件等に応じて決定するので、撮像画像の色再現性及びノイズの低減性を向上させることができる。
【０２０６】
【発明の効果】
以上詳細に説明したように、本発明に係る撮像装置は、人間の眼の見え方に対して忠実な色の再現を表す色再現値と人間が感じるノイズ感を表すノイズ値を調整する調整部と、上記調整部の調整に基づきマトリクス係数を決定するマトリクス係数決定部と、上記マトリクス係数に基づき当該撮像装置に備えられている撮像素子部で撮像した画像にマトリクス変換処理を行なうマトリクス変換処理部とを備えるので、撮像環境及び条件に応じて適応的にリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なうことができる。
【０２０７】
また、本発明に係る撮像方法は、分光特性の異なるカラーフィルタによりなっており、人間の眼の見え方に対して忠実な色の再現を表す色再現値と人間が感じるノイズ感を表すノイズ値を調整する第１のステップと、上記第１のステップの調整に基づきマトリクス係数を決定する第２のステップと、上記マトリクス係数に基づき被写体を撮像する撮像素子部で撮像した画像にマトリクス変換処理を行なう第３のステップとを有するので、撮像環境及び条件に応じて適応的にリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なうことができる。
【図面の簡単な説明】
【図１】従来のRGBカラーフィルタの例を示す図である。
【図２】従来の撮像装置に設けられる信号処理部の構成例を示すブロック図である。
【図３】分光感度特性の例を示す図である。
【図４】分光感度特性の他の例を示す図である。
【図５】本発明を適用した撮像装置の構成例を示すブロック図である。
【図６】本発明を適用した撮像装置に設けられる４色カラーフィルタの例を示す図である。
【図７】視感度曲線の例を示す図である。
【図８】評価係数の特徴を示す図である。
【図９】本発明を適用した撮像装置が備えるカメラシステムLSIの構成例を示すブロック図である。
【図１０】図９の信号処理部の構成例を示すブロック図である。
【図１１】画像処理装置の作成処理を説明するフローチャートである。
【図１２】図１１のステップＳ１の４色カラーフィルタ決定処理の詳細を説明するフローチャートである。
【図１３】仮想曲線の例を示す図である。
【図１４】フィルタ毎のUMG値の例を示す図である。
【図１５】４色カラーフィルタの分光感度特性の例を示す図である。
【図１６】図１１のステップＳ２のリニアマトリクス決定処理の詳細を説明するフローチャートである。
【図１７】色差の評価結果の例を示す図である。
【図１８】４色カラーフィルタによる所定の物体の色度を示す図である。
【図１９】本発明を適用した撮像装置に設けられる４色カラーフィルタの他の例を示す図である。
【図２０】適応的にリニアマトリクス係数Ｍの決定を示すフローチャートである。
【図２１】色再現性指標を変化させたときのノイズ低減性指標の変化の様子を示す図である。
【図２２】４色カラーフィルタの分光感度特性の例を示す図である。
【図２３】画像のヒストグラムを示す図である。
【符号の説明】
４５ 画像センサ、６１ ４色カラーフィルタ、７１ 信号処理部、９１ オフセット補正処理部、９２ ホワイトバランス補正処理部、９３ 垂直方向同時化処理部、９４ 信号生成処理部、９５ リニアマトリクス処理部、９６−１，９６−２，９６−３ ガンマ補正処理部、９７ 輝度信号生成処理部、９８ 色差信号生成処理部、９９ 帯域制限間引き処理部
【発明の属する技術分野】
本発明は、被写体を撮像する撮像装置及び方法に関し、特に撮像感度に対して好適に撮像を行なう撮像装置及び方法に関する。
【０００２】
【従来の技術】
近年、コンシューマ向けの画像撮像装置（ディジタルカメラ、カラースキャナ等）および画像処理ソフトウェアが普及しており、当該画像撮像装置により被写体を撮像して生成した画像を画像処理ソフトウェアにより自ら編集するユーザが増加している。また、画像撮像装置により撮像された画像の画質に対する要求も高くなってきており、画像撮像装置を購入する際の主要な条件として、画質がよいことを挙げるユーザの割合が増えてきている。ここで、画像撮像装置の一般的な構成について以下に述べる。
【０００３】
画像撮像装置は、例えば、図１に示すようなＲＧＢの３原色のカラーフィルタ１が用いられている。この例では、図１の一点鎖線で示すように、緑色（Ｇ）の光のみを透過するＧフィルタが２個、赤（Ｒ）の光のみを透過するＲフィルタが１個、および青（Ｂ）の光のみを透過するＢフィルタが１個の、合計４個を最小単位として、いわゆるベイヤー配列（Bayer配列）により、カラーフィルタ１が構成される。
【０００４】
図２は、ＲＧＢカラーフィルタ１を有するＣＣＤ(Charge Coupled Device)撮像素子により取得されたＲＧＢ信号に対して各種の処理を施す信号処理部１１の構成例を示すブロック図である。
【０００５】
オフセット補正処理部２１は、ＣＣＤ撮像素子により取得された信号に所定の処理を施すフロントエンド１３から供給されてきた画像信号に含まれるオフセット成分を除去し、得られた画像信号をホワイトバランス補正処理部２２に出力する。ホワイトバランス補正処理部２２は、オフセット補正処理部２１から供給されてきた画像信号の色温度、およびカラーフィルタ１の各フィルタの感度の違いに基づいて、各色のバランスを補正する。ホワイトバランス補正処理部２２により補正が施され、取得された色信号はガンマ補正処理部２３に出力される。ガンマ補正処理部２３は、ホワイトバランス補正処理部２２から供給されてきた信号に対してガンマ補正を行い、取得した信号を垂直方向同時化処理部２４に出力する。垂直方向同時化処理部２４にはディレイ素子が設けられており、ガンマ補正処理部２３から供給されてきた信号の垂直方向の時間のずれが同時化される。
【０００６】
ＲＧＢ信号生成処理部２５は、垂直方向同時化処理部２４から供給されてきた色信号を、同一空間の位相に補間する補間処理、信号のノイズ成分を除去するノイズ除去処理、信号帯域を制限するフィルタリング処理、および信号帯域の高域成分を補正する高域周波数補正処理等を行い、得られたＲＧＢ信号を、輝度信号生成処理部２６、および色差信号生成処理部２７に出力する。
【０００７】
輝度信号生成処理部２６は、ＲＧＢ信号生成処理部２５から供給されてきたＲＧＢ信号を所定の合成比で合成し、輝度信号（Ｙ）を生成する。色差信号生成処理部２７も同様に、ＲＧＢ信号生成処理部２５から供給されてきたＲＧＢ信号を所定の合成比で合成し、色差信号（Ｃｂ，Ｃｒ）を生成する。輝度信号生成処理部２６により生成された輝度信号（Ｙ）、および色差信号生成処理部２７により生成された色差信号（Ｃｂ，Ｃｒ）は、例えば、信号処理部１１の外部に設けられているモニタに出力される。
【０００８】
このように、元信号に対してガンマ処理を施してから線形変換によって画像処理（リニアマトリクス処理）を行うことが一般的に行われている。
【０００９】
上述したような画像撮像装置は、被写体を撮像し画像を生成したとき、観察時の視環境によって見え方が異なるために、所望の色で再現されない場合がある。これは、撮像時の光源（以下、撮像光源という。）と観察時の光源（以下、観察光源という。）の演色性が大きく異なる場合に生じる現象である。そこで、撮像光源と演色性の異なる観察光源で画像が再現される場合でも、良好に色再現する方法が提案されている（例えば、特許文献１参照。）。
【００１０】
また、図３に分光感度の特性を示す。図３の曲線Ｌ１は、Ｒの分光感度を示し、曲線Ｌ２は、Ｇの分光感度を示し、曲線Ｌ３は、Ｂの分光感度を示している。
【００１１】
【特許文献１】
特開２００２−１４２２３１号公報
【００１２】
【発明が解決しようとする課題】
一方、上述したような画像撮像装置に、リニアマトリクス処理の追加、更に撮像素子のカラーフィルタの追加を行うことにより、色再現性を劇的に向上させることができる。その際、使用するリニアマトリクスの係数を単に色差が最小になるように決定すると、ノイズを従来よりも増長させてしまうことがある。
【００１３】
本発明は、上述したような実情に鑑みて提案されたものであり、撮像する環境や条件等に応じて色再現性とノイズ低減性を考慮した係数を用いたリニアマトリクス処理を行なうことが可能な撮像装置及び撮像方法を提供することを目的とする。
【００１４】
【課題を解決するための手段】
本発明に係る撮像装置は、上述の課題を解決するために、分光特性の異なるカラーフィルタによりなっており、被写体を撮像する撮像素子部を有し、該撮像素子部で撮像した画像にマトリクス変換処理を行って色補正を行う撮像装置において、色再現性指標ΔＥ（Ｍ）とノイズ低減性指標σＮ（Ｍ）とから評価指数ＥＥＶ（Ｍ）を以下の条件により導出し、導出された上記評価指数ＥＥＶ（Ｍ）に基づいてマトリクス係数を決定し、決定された上記マトリクス係数に基づいて上記撮像素子部で撮像した画像にマトリクス変換処理を行うマトリクス変換処理手段を備える。
（条件）
ＥＥＶ（Ｍ）＝√｛（ｗｃ・ΔＥ（Ｍ）） ２ ＋（ｗｎ・σＮ（Ｍ）） ２ ｝
但し、ＥＥＶ（Ｍ）：評価指数、ΔＥ（Ｍ）：色再現性指標、σＮ（Ｍ）：ノイズ低減性指標、ｗｃ：色再現性に対する重み係数、ｗｎ：ノイズ低減性に対する重み係数
【００１６】
【発明の実施の形態】
図５は、本発明を適用した撮像装置の構成例を示すブロック図である。
【００１７】
図５に示す撮像装置は、CCD(Charge Coupled Device)等よりなる画像センサ４５の前面（レンズ４２に対向する面）に、４種類の色（光）を識別するカラーフィルタが設けられている。なお、撮像装置では、図５の画像センサ４５に設けられるカラーフィルタを図６に示す４色カラーフィルタ６１とする。
【００１８】
４色カラーフィルタ６１は、図６の一点鎖線で示すように、赤の光のみを透過するＲフィルタ、青の光のみを透過するＢフィルタ、第１の波長帯域の緑色の光のみを透過するＧ１フィルタ、およびＧ１フィルタと相関が高い、第２の波長帯域の緑色の光のみを透過するＧ２フィルタの、合計４個のフィルタを最小単位として構成される。また、Ｇ１フィルタとＧ２フィルタは、その最小単位内において、互いに対角する位置に配置されている。
【００１９】
後に詳述するように、画像センサ４５により取得される画像の色を４種類とし、取得される色情報を増やすことにより、３種類の色（RGB）のみが取得される場合と較べて、より正確に色を表現することができ、眼で違うものに見える色は違う色に、同じものに見える色は同じ色に、それぞれ再現すること（「色の判別性」）を向上させることができる。
【００２０】
なお、図７に示される視感度曲線から分かるように、人間の眼は輝度に敏感である。従って、図６に示す４色カラーフィルタ６１においては、より正確な輝度情報を取得することにより、輝度の階調を上げることができ、かつ、眼の見え方に近い画像を再現することができるように、視感度曲線に近い分光感度特性を有するＧ２のカラーフィルタが追加されている（図１のＲ，Ｇ，Ｂに対応するＲ，Ｇ１，Ｂのフィルタに対して、新たに決定された緑Ｇ２フィルタが追加されている）。
【００２１】
また、４色カラーフィルタ６１を決定する際に使用されるフィルタ評価係数として、例えば、「色再現性」と「ノイズ低減性」の両方を考慮した係数であるUMG(Unified Measure of Goodness)が用いられている。
【００２２】
UMGを用いた評価においては、評価対象のフィルタが単にルータ条件を満たしているだけでは、その評価値は高くならず、それぞれのフィルタの分光感度分布の重なりも考慮される。従って、ｑファクタ、μファクタ、またはFOM(Figure of Merit)を利用して評価されたカラーフィルタの場合と比較して、ノイズをより低減させることができる。すなわち、UMGを用いた評価により、それぞれのフィルタの分光感度特性が、ある程度の重なりを有するが、図４のＲの特性とＧの特性のようにほぼ全てが重なるものではないフィルタが選択されるため、色の分離のため、それぞれの色信号を増幅する場合であっても、増幅率をそれ程大きくする必要がなく、それに伴ってノイズ成分が増幅されることが抑制される。
【００２３】
図８は、各フィルタ評価係数の特徴を示す図である。なお、図８は、各評価係数に対して、一度に評価できるフィルタの数、物体の分光反射率が考慮されているか否か、およびノイズの低減が考慮されているか否かに関する事項を示すものである。
【００２４】
図８に示すように、ｑファクタ（q-factor）は、一度に評価できるフィルタの数が「１個」のみであり、物体の分光反射率と、ノイズの低減が考慮されていない。また、μファクタ（μ-factor）は、一度に複数のフィルタを評価することができるものの、物体の分光反射率とノイズの低減は考慮されていない。さらに、FOMは、一度に複数のフィルタを評価することができ、物体の分光反射率が考慮されているものの、ノイズの低減が考慮されていない。
【００２５】
これに対して、４色カラーフィルタ６１を決定する際に使用されるUMGは、一度に複数のフィルタを評価することができ、物体の分光反射率が考慮され、かつ、ノイズの低減も考慮されている。
【００２６】
なお、ｑファクタについては「H.E.J.Neugebauer "Quality Factor for Filters Whose Spectral Transmittances are Different from Color Mixture Curves, and Its Application to Color Photography" JOURNAL OF THE OPTICAL SOCIETY OF AMERICA, VOLUME 46, NUMBER 10」にその詳細が開示されており、μファクタについては「P.L.Vora and H.J.Trussell, "Measure of Goodness of a set of color-scanning filters", JOURNAL OF THE OPTICAL SOCIETY OF AMERICA, VOLUME 10, NUMBER 7」にその詳細が開示されている。また、FOMについては「G.Sharma and H.J.Trussell, "Figures of Merit for Color Scanners", IEEE TRANSACTION ON IMAGE PROCESSING, VOLUME 6」にその詳細が開示されており、UMGについては「S.Quan, N.Ohta, and N.Katoh, "Optimal Design of Camera Spectral Sensitivity Functions Based on Practical Filter Components", CIC, 2001」にその詳細が開示されている。
【００２７】
図５の説明に戻る。マイクロコンピュータ４１は、所定の制御プログラムに従って全体の動作を制御する。例えば、マイクロコンピュータ４１は、絞り４３による露光制御、シャッター４４の開閉制御、TG(Timing Generator)４６の電子シャッターの制御、フロントエンド４７でのゲインコントロール、カメラシステムLSI(Large Scale Integrated Circuit)４８のモード制御、パラメータ制御等を行う。
【００２８】
絞り４３は、レンズ４２により集光された光の通過（絞り）を調整し、画像センサ４５により取り込まれる光量を制御する。シャッター４４は、マイクロコンピュータ４１の指示に基づいて、レンズ４２により集光された光の通過を制御する。
【００２９】
画像センサ４５は、さらに、CCDやCMOS(Complementary Metal Oxide Semiconductor)により構成される撮像素子を有し、その撮像素子の前面に形成する４色カラーフィルタ６１を介して入射された光を電気信号に変換し、４種類の色信号（Ｒ信号，Ｇ１信号，Ｇ２信号，Ｂ信号）をフロントエンド４７に出力する。画像センサ４５には、図６の４色カラーフィルタ６１が設けられ、レンズ４２を介して入射された光から、Ｒ，Ｇ１，Ｇ２，Ｂのそれぞれの帯域の波長の成分が抽出される。なお、その詳細については、図１５を参照して後述する。
【００３０】
フロントエンド４７は、画像センサ４５から供給されてきた色信号に対して、ノイズ成分を除去するための相関二重サンプリング処理、ゲインコントロール処理、およびディジタル変換処理等を施す。フロントエンド４７により各種の処理が施され、得られた画像データは、カメラシステムLSI４８に出力される。
【００３１】
カメラシステムLSI４８は、後に詳述するように、フロントエンド４７から供給されてきた画像データに対して各種の処理を行い、例えば、輝度信号および色信号を生成して画像モニタ５０に出力し、信号に対応する画像を表示させる。
【００３２】
画像メモリ４９は、例えば、DRAM(Dynamic Random Access Memory)やSDRAM(Synchronous Dynamic Random Access Memory)などにより構成され、カメラシステムLSI４８が各種の処理を行う際に適宜利用される。半導体メモリ、ディスク等により構成される外部記憶媒体５１は、例えば、図５の撮像装置に対して着脱可能に構成され、カメラシステムLSI４８によりJPEG(Joint Photographic Expert Group)フォーマットで圧縮された画像データが記憶される。
【００３３】
画像モニタ５０は、例えば、LCD(Liquid Crystal Display)などにより構成され、撮影された画像や各種のメニュー画面等を表示する。
【００３４】
図９は、図５に示すカメラシステムLSI４８の構成例を示すブロック図である。カメラシステムLSI４８を構成する各ブロックは、マイクロコンピュータインタフェース(I/F)７３を介して、図５に示すマイクロコンピュータ４１により制御される。
【００３５】
信号処理部７１は、フロントエンド４７から供給されてきた４種類の色情報に対して、補間処理、フィルタリング処理、マトリクス演算処理、輝度信号生成処理、色差信号生成処理等の各種の処理を行い、例えば、生成した画像信号をモニタインタフェース７７を介して画像モニタ５０に出力する。
【００３６】
画像検波部７２は、フロントエンド４７の出力に基づいて、オートフォーカス、オートエキスポージャー、オートホワイトバランス等の検波処理を行い、その結果を、適宜、マイクロコンピュータ４１に出力する。
【００３７】
メモリコントローラ７５は、処理ブロック間同士のデータの送受信、または、所定の処理ブロックと画像メモリ４９との間のデータの送受信を制御し、例えば、信号処理部７１から供給されてきた画像データをメモリインタフェース７４を介して画像メモリ４９に出力し、記憶させる。
【００３８】
画像圧縮・解凍部７６は、例えば、信号処理部７１から供給されてきた画像データをJPEGフォーマットで圧縮し、得られたデータをマイクロコンピュータインタフェース７３を介して外部記憶媒体５１に出力し、記憶させる。画像圧縮・解凍部７６は、また、外部記憶媒体５１から読み出された圧縮データを解凍（伸張）し、モニタインタフェース７７を介して画像モニタ５０に出力する。
【００３９】
図１０は、図９に示す信号処理部７１の詳細な構成例を示すブロック図である。信号処理部７１を構成する各ブロックは、マイクロコンピュータインタフェース７３を介して、マイクロコンピュータ４１により制御される。
【００４０】
オフセット補正処理部９１は、フロントエンド４７から供給されてきた画像信号に含まれるノイズ成分（オフセット成分）を除去し、得られた画像信号をホワイトバランス補正処理部９２に出力する。ホワイトバランス補正処理部９２は、オフセット補正処理部９１から供給されてきた画像信号の色温度、および４色カラーフィルタ６１の各フィルタの感度の違いに基づいて、各色のバランスを補正する。ホワイトバランス補正処理部９２により補正が施され、取得された色信号は、垂直方向同時化処理部９３に出力される。垂直方向同時化処理部９３にはディレイ素子が設けられており、ホワイトバランス補正処理部９２から出力されてきた信号（以下、RG1G2B信号という。）の垂直方向の時間のずれが同時化（補正）される。
【００４１】
信号生成処理部９４は、垂直方向同時化処理部９３から供給されてきたRG1G2B信号の最小単位の２×２画素の色信号を、同一空間の位相に補間する補間処理、信号のノイズ成分を除去するノイズ除去処理、信号帯域を制限するフィルタリング処理、および信号帯域の高域成分を補正する高域周波数補正処理等を行い、得られたRG1G2B信号を、リニアマトリクス処理部９５に出力する。
【００４２】
リニアマトリクス処理部９５は、所定のリニアマトリクス係数Ｍ（３×４の行列）に基づいて式（１）によるRG1G2B信号の演算を行ない、３色のRGB信号を生成する。
【００４３】
【数１】
リニアマトリクス処理部９５により生成されたＲ信号は、ガンマ補正処理部９６−１に出力され、Ｇ信号は、ガンマ補正処理部９６−２に出力され、Ｂ信号は、ガンマ補正処理部９６−３に出力される。
【００４５】
ガンマ補正処理部９６−１乃至９６−３は、リニアマトリクス処理部９５から出力されてきたRGB信号のそれぞれの信号に対してガンマ補正を行い、取得されたRGB信号を輝度（Ｙ）信号生成処理部９７および色差（Ｃ）信号生成処理部９８に出力する。
【００４６】
輝度信号生成処理部９７は、ガンマ補正処理部９６−１乃至９６−３から供給されてきたRGB信号を、例えば、式（２）に従って、所定の合成比で合成し、輝度信号（Ｙ）を生成する。
Ｙ＝０．２１２６Ｒ＋０．７１５２Ｇ＋０．０７２２Ｂ・・・（２）
色差信号生成処理部９８も同様に、ガンマ補正処理部９６−１乃至９６−３から供給されてきたRGB信号を所定の合成比で合成して色差信号（Ｃ）を生成し、帯域制限間引き処理部９９に出力する。帯域制限間引き処理部９９は、色差信号（Ｃ）に基づき色差信号（Ｃｂ，Ｃｒ）を生成する。なお、単板２×２色コーディングによって得られた信号には、一般的に色情報の帯域が輝度信号ほど存在しない。したがって、帯域制限間引き処理部９９は、色差信号生成処理部９８から供給された色差信号（Ｃ）を帯域制限処理と間引き処理を行なうことで色情報データを削減することで色差信号（Ｃｂ，Ｃｒ）を生成する。
【００４７】
輝度信号生成処理部９７により生成された輝度信号（Ｙ）、および色差信号生成処理部９８により生成された色差信号（Ｃ）又は帯域制限間引き処理部９９により生成された色差信号（Ｃｂ，Ｃｒ）は、例えば、図９に示すモニタインタフェース７７を介して画像モニタ５０に出力される。
【００４８】
以上のような構成を有する撮像装置において、画像の撮影が指示された場合、マイクロコンピュータ４１は、TG４６を制御し、画像センサ４５により画像を取り込ませる。すなわち、画像センサ４５を構成するCCD等の撮像素子（以下、CCD撮像素子という。）の前面に形成されている４色カラーフィルタ６１により４色の光が透過され、透過された光がCCD撮像素子により取り込まれる。CCD撮像素子により取り込まれた光は、４色の色信号に変換され、それがフロントエンド４７に出力される。
【００４９】
フロントエンド４７は、画像センサ４５から供給されてきた色信号に対して、ノイズ成分を除去するための相関二重サンプリング処理、ゲインコントロール処理、およびディジタル変換処理等を施し、得られた画像データをカメラシステムLSI４８に出力する。
【００５０】
カメラシステムLSI４８の信号処理部７１においては、オフセット補正処理部９１により色信号のオフセット成分が除去され、ホワイトバランス補正処理部９２により、画像信号の色温度、および４色カラーフィルタ６１の各フィルタの感度の違いに基づいて、各色のバランスが補正される。
【００５１】
また、垂直方向同時化処理部９３により、ホワイトバランス補正処理部９２により補正された信号の垂直方向の時間のずれが同時化（補正）され、信号生成処理部９４により、垂直方向同時化処理部９３から供給されてきたRG1G2B信号の最小単位の２×２画素の色信号を、同一空間の位相に補間する補間処理、信号のノイズ成分を除去するノイズ除去処理、信号帯域を制限するフィルタリング処理、および信号帯域の高域成分を補正する高域周波数補正処理等が行われる。
【００５２】
さらに、リニアマトリクス処理部９５においては、信号生成処理部９４により生成された信号（RG1G2B信号）が、所定のリニアマトリクス係数Ｍ（３×４の行列）に基づいて変換され、３色のRGB信号が生成される。リニアマトリクス処理部９５により生成されたＲ信号は、ガンマ補正処理部９６−１に出力され、Ｇ信号は、ガンマ補正処理部９６−２に出力され、Ｂ信号は、ガンマ補正処理部９６−３に出力される。
【００５３】
ガンマ補正処理部９６−１乃至９６−３により、リニアマトリクス処理部９５の処理により得られたRGB信号のそれぞれの信号に対してガンマ補正が行われ、取得されたRGB信号が輝度信号生成処理部９７および色差信号生成処理部９８に出力される。輝度信号生成処理部９７、および色差信号生成処理部９８においては、ガンマ補正処理部９６−１乃至９６−３から供給されてきたＲ信号、Ｇ信号、Ｂ信号のそれぞれの信号が所定の合成比で合成され、輝度信号（Ｙ）および色差信号（Ｃ）が生成される。輝度信号生成処理部９７により生成された輝度信号（Ｙ）、および色差信号生成処理部９８により生成された色差信号（Ｃ）は、図９の画像圧縮・解凍部７６に出力され、例えば、JPEGフォーマットで圧縮される。圧縮され、得られた画像データは、マイクロコンピュータインタフェース７３を介して外部記憶媒体５１に出力され、記憶される。
【００５４】
以上のように、１つの画像データが４種類の色信号に基づいて形成されるため、その再現性は、人間の眼の見えにより近いものとなる。
【００５５】
一方、外部記憶媒体５１に記憶されている画像データの再生（表示）が指示されたとき、マイクロコンピュータ４１により、外部記憶媒体５１に記憶されている画像データが読み出され、それがカメラシステムLSI４８の画像圧縮・解凍部７６に出力される。画像圧縮・解凍部７６においては、圧縮されている画像データが伸張され、モニタインタフェース７７を介して、得られたデータに対応する画像が画像モニタ５０に表示される。
【００５６】
次に、図１１に示すフローチャートを参照して、以上のような構成を有する撮像装置を作成する処理（手順）について説明する。
【００５７】
ステップＳ１において、図５に示す画像センサ４５に設けられる４色カラーフィルタ６１の分光感度特性を決定する４色カラーフィルタ決定処理が行われ、ステップＳ２において、図１０に示すリニアマトリクス処理部９５に設定されるマトリクス係数Ｍを決定するリニアマトリクス係数Ｍ決定処理が行われる。ステップＳ１において実行される４色カラーフィルタ決定処理の詳細については、図１２に示すフローチャートを参照して、また、ステップＳ２において実行されるリニアマトリクス係数Ｍ決定処理の詳細については、図１６に示すフローチャートを参照して、それぞれ後述する。
【００５８】
４色カラーフィルタ６１が決定され、マトリクス係数が決定された後、ステップＳ３において、図１０に示す信号処理部７１が作成され、ステップＳ４に進み、図９に示すカメラシステムLSI４８が作成される。また、ステップＳ５において、図５に示すような撮像装置（例えば、ディジタルカメラ）の全体が作成される。ステップＳ６において、ステップＳ５で作成された撮像装置の画質（「色再現性」、「色判別性」）の評価が行われ、処理が終了される。
【００５９】
ここで、「色再現性」、「色判別性」などを評価する際に参照される物体色について説明する。物体色は、「物体の分光反射率」、「標準照明の分光エネルギー分布」、および「物体を感知するセンサ（カラーフィルタ）の分光感度分布（特性）」の積を可視光領域（例えば、４００乃至７００nm）の範囲で積分した値によって算出される。すなわち、式（３）により物体色が算出される。
【００６０】
【数２】
例えば、所定の物体を眼で観察する場合、式（３）の「センサの分光感度特性」が等色関数で表され、その物体の物体色は、Ｘ，Ｙ，Ｚの三刺激値で表される。具体的には、Ｘの値は式（４−１）で算出され、Ｙの値は式（４−２）で算出され、Ｚの値は式（４−３）で算出される。なお、式（４−１）乃至（４−３）における定数ｋの値は式（４−４）で算出される。
【００６２】
【数３】
また、ディジタルカメラなどの撮像装置により、所定の物体の画像を取り込む場合、式（３）の「センサの分光感度特性」がカラーフィルタの分光感度特性で表され、その物体の物体色は、フィルタの数の色値（例えば、RGBフィルタ（３種類）の場合はRGB値（３値））の物体色が算出される。撮像装置に、３種類の色を検出するRGBフィルタが設けられている場合、具体的には、Ｒの値は式（５−１）で算出され、Ｇの値は式（５−２）で算出され、Ｂの値は式（５−３）で算出される。また、式（５−１）における定数ｋrの値は式（５−４）で算出され、式（５−２）における定数ｋgの値は式（５−５）で算出され、式（５−３）における定数ｋbの値は式（５−６）で算出される。
【００６４】
【数４】
次に、図１２に示すフローチャートを参照して、図１１に示すステップＳ１で行われる４色カラーフィルタ決定処理について説明する。
【００６６】
なお、４色カラーフィルタの決定方法としては様々な方法があるが、例えば、RGBフィルタを基調として（既存の（図１）Ｇフィルタの一方をＧ１フィルタとして）、Ｇ１フィルタを透過する色と相関の高い色を透過するＧ２フィルタを選択し、それを追加して４色カラーフィルタを決定する処理について説明する。
【００６７】
ステップＳ２１において、UMG値を算出するために使用されるカラーターゲットが選択される。例えば、ステップＳ２１において、現存する色を代表するカラーパッチを多く含み、かつ、人間の記憶色（肌色、植物の緑、空の青等）を重視したカラーパッチを多く含むカラーターゲットが選択される。カラーターゲットとしては、例えば、IT8.7，Macbeth Color Checker，GretagMacbeth DigitalCamera Color Checker，CIE，Color Bar等がある。
【００６８】
また、目的に応じて、SOCS(Standard Object Color Spectra Database)などのデータから、標準となりうるカラーパッチを作成し、それを使用するようにしてもよい。なお、SOCSについては「田島譲二, "標準物体色分光データベース(SOCS)による統計的色再現評価", カラーフォーラムJAPAN 99」にその詳細が開示されている。以下、Macbeth Color Checkerがカラーターゲットとして選択された場合について説明する。
【００６９】
ステップＳ２２において、Ｇ２フィルタの分光感度特性が決定される。分光感度特性として、実存する材料から作成可能なものを使用するようにしてもよいし、図１３に示すようなcubic spline曲線（３次スプライン関数）で仮想曲線Ｃ（λ）を想定し、仮想曲線Ｃ（λ）のピーク値λ0、値ｗ（ｗ1とｗ2の和を２で除算した値）、値Δｗ（ｗ1からｗ2を減算した値を２で除算した値）を、図中に示す範囲で変化させたものを使用するようにしてもよい。なお、ｗ，Δｗの値は、半値幅の値に基づく値とされる。λ0，ｗ，Δｗの変化の方法は、例えば、５nm刻みとする。仮想曲線Ｃ（λ）は、それぞれの範囲において、下式（６−１）乃至（６−５）で表される。
【００７０】
【数５】
なお、この例では、フィルタＧ２のみが追加されるが、図１に示すフィルタ（Ｒ，Ｇ，Ｇ，Ｂ）のＲフィルタとＢフィルタのみを使用し、残りの２つのＧ１，Ｇ２フィルタを、緑色付近の上式（６−１）乃至（６−５）の仮想曲線として定義することも可能である。また、同様に、ＲとＧのみ、ＧとＢのみを、図１に示すフィルタから使用するようにしてもよい。さらに、４色のフィルタのうち、３色を仮想曲線、または４色とも仮想曲線として定義することも可能である。
【００７２】
ステップＳ２３において、追加するフィルタ（Ｇ２フィルタ）と、現存のフィルタ（Ｒフィルタ，Ｇ１フィルタ，Ｂフィルタ）が組み合わされ、４色カラーフィルタの最小単位（セット）が作成される。また、ステップＳ２４において、ステップＳ２３で作成された４色カラーフィルタに対して、フィルタ評価係数としてUMGが用いられ、UMG値が算出される。
【００７３】
図８を参照して説明したように、UMGを用いた場合、４色のそれぞれのカラーフィルタに対して一度で評価を行うことができる。また、物体の分光反射率を考慮して評価が行われるだけでなく、ノイズの低減性をも考慮して評価が行われる。UMGを用いた評価においては、それぞれのフィルタの分光感度特性に、適度な重なりがあるフィルタに対して高い評価が示されるため、例えば、Ｒの特性とＧの特性が広い波長帯域にわたって重複する特性を有するフィルタ（それぞれの色信号を分離した際、ノイズが増幅されるフィルタ）に対して、高い評価が示されるということを抑制することができる。
【００７４】
図１４は、３色カラーフィルタにおいて算出されるUMG値の例を示す図である。例えば、RGBの特性がそれぞれ重複しない図１４Ａに示すような特性のフィルタにおいては、「０．７９４２」のUMG値が算出され、Ｒの特性とＧの特性が広い波長帯域にわたって重複する図１４Ｂに示すような特性のフィルタにおいては、「０．８２１１」のUMG値が算出される。また、RGBのそれぞれの特性が適度に重なる図１４Ｃに示すような特性を有するフィルタにおいては、「０．８８７９」のUMG値が算出される。すなわち、RGBのそれぞれの特性が適度に重なる、図１４Ｃに示すような特性を有するフィルタに対して最も高い評価が示される。このことは、４色カラーフィルタにおいても同様である。なお、図１４Ａに示す曲線Ｌ３１、図１４Ｂに示す曲線Ｌ４１、および図１４Ｃに示す曲線Ｌ５１はＲの分光感度を表し、図１４Ａに示す曲線Ｌ３２、図１４Ｂに示す曲線Ｌ４２、および図１４Ｃに示す曲線Ｌ５２はＧの分光感度を表し、図１４Ａに示す曲線Ｌ３３、図１４Ｂに示す曲線Ｌ４３、および図１４Ｃに示す曲線Ｌ５３はＢの分光感度を表している。
【００７５】
ステップＳ２５において、ステップＳ２４で算出されたUMG値が、所定の閾値である「０．９５」以上であるか否かが判定され、「０．９５」未満であると判定された場合、ステップＳ２６に進み、作成された４色カラーフィルタが却下される（使用されない）。ステップＳ２６において、４色カラーフィルタが却下された場合、その後、処理は終了される（図１１に示すステップＳ２以降の処理は実行されない）。
【００７６】
一方、ステップＳ２５において、ステップＳ２４で算出されたUMG値が「０．９５」以上であると判定された場合、ステップＳ２７において、その４色カラーフィルタが、ディジタルカメラで使用される候補のフィルタとされる。
【００７７】
ステップＳ２８において、ステップＳ２７で候補のフィルタとされた４色カラーフィルタが、現存する材料、染料で実現可能であるか否かが判定される。材料、染料等が取得困難である場合、実現不可能であると判定され、ステップＳ２６に進み、その４色カラーフィルタが却下される。
【００７８】
一方、ステップＳ２８において、材料、染料等が取得可能であり、実現可能であると判定された場合、ステップＳ２９に進み、作成された４色カラーフィルタが、ディジタルカメラにおいて使用されるフィルタとして決定される。その後、図１１に示すステップＳ２以降の処理が実行される。
【００７９】
図１５は、ステップＳ２９において決定された４色カラーフィルタの分光感度特性の例を示す図である。
【００８０】
図１５において、曲線Ｌ６１はＲの分光感度を表し、曲線Ｌ６２はＧ１の分光感度を表している。また、曲線Ｌ６３はＧ２の分光感度を表し、曲線Ｌ６４はＢの分光感度を表している。図１５に示すように、Ｇ２の分光感度曲線（曲線Ｌ６３）は、Ｇ１の分光感度曲線（曲線Ｌ６２）に対して相関が高い。また、Ｒの分光感度、Ｇ（Ｇ１，Ｇ２）の分光感度、およびＢの分光感度は、それぞれ適度な範囲において重複している。
【００８１】
以上のようにして決定された４色カラーフィルタを利用することにより、特に、「色再現性」のうちの「色の判別性」を改善させることができる。
【００８２】
なお、光の利用効率の観点から、以上のように、現存するRGBフィルタのGフィルタと相関の高いフィルタを、追加するフィルタ（Ｇ２フィルタ）とすることが好ましい。この場合、追加するフィルタの分光感度曲線のピーク値は、経験的に４９５乃至５３５nmの範囲（現存するＧフィルタの分光感度曲線のピーク値の近傍）にあることが望ましい。
【００８３】
また、現存するＧフィルタと相関の高いフィルタを追加する場合、図１に示す最小単位（Ｒ，Ｇ，Ｇ，Ｂ）を構成する２つのＧフィルタのいずれか一方を追加色のフィルタにするだけで４色カラーフィルタを作成することができるため、作成の行程に大きな変更を加える必要がない。
【００８４】
以上のようにして４色カラーフィルタが作成され、それがディジタルカメラに設けられた場合、図１０に示す信号処理装置７１には４種類の色信号が信号生成処理部９４から供給されるため、リニアマトリクス処理部９５において、４色（Ｒ，Ｇ１，Ｇ２，Ｂ）の信号から３色（Ｒ，Ｇ，Ｂ）の信号を生成する変換処理が行われる。この変換処理は、輝度リニアな（輝度値を線形的な変換により表すことができる）入力信号値に対するマトリクス処理であるため、以下、リニアマトリクス処理部９５において行われる変換処理を、適宜、リニアマトリクス処理と称する。
【００８５】
次に、図１６に示すフローチャートを参照して、図１１に示すステップＳ２において実行される、リニアマトリクス係数Ｍ決定処理について説明する。なお、リニアマトリクス係数Ｍ決定処理において使用されるカラーターゲットはMacbeth Color Checkerとし、また、使用される４色カラーフィルタは図１５に示す分光感度特性を有するものとする。
【００８６】
ステップＳ４１において、例えば、CIE(Commision Internationale del'Eclairange)において標準光源とされている一般的な昼光Ｄ６５（照明光Ｌ（λ））が照明光として選択される。なお、照明光は、画像処理装置が頻繁に使用されることが予想される環境の照明光等に変更してもよい。また、想定される照明環境が複数ある場合、リニアマトリクスを複数用意することも考えられる。以下、照明光として昼光Ｄ６５が選択された場合について説明する。
【００８７】
ステップＳ４２において、リファレンス値（参照値）Ｘｒ，Ｙｒ，Ｚｒが算出される。具体的には、リファレンス値Ｘｒは式（７−１）により算出され、Ｙｒは式（７−２）により算出され、Ｚｒは式（７−３）により算出される。
【００８８】
【数６】
また、定数ｋは式（８）により算出される。
【００９０】
【数７】
例えば、カラーターゲットがMacbeth Color Checkerの場合、２４色分のリファレンス値が算出される。
【００９２】
次に、ステップＳ４３において、４色カラーフィルタの出力値Ｒf，Ｇ１f，Ｇ２f，Ｂfが算出される。具体的には、Ｒfは式（９−１）により算出され、Ｇ１fは式（９−２）により算出され、Ｇ２fは式（９−３）により算出され、Ｂfは式（９−４）により算出される。
【００９３】
【数８】
また、定数ｋrは式（１０−１）により算出され、定数ｋg1は式（１０−２）により算出され、定数ｋg2は式（１０−３）により算出され、定数ｋbは式（１０−４）により算出される
【００９５】
【数９】
例えば、カラーターゲットがMacbeth Color Checkerの場合は、２４色分の出力値Ｒf，Ｇ１f，Ｇ２f，Ｂfが算出される。
【００９７】
ステップＳ４４において、ステップＳ４３で算出されたフィルタ出力値を、ステップＳ４２で算出されたリファレンス値（XYZref）に近似させる変換を行うマトリクスが、例えば、XYZ色空間における誤差最小二乗法により算出される。
【００９８】
例えば、算出される、３×４のマトリクスを式（１１）で表されるＡとした場合、マトリクス変換（XYZexp）は、次の式（１２）で表される。
【００９９】
【数１０】
また、リファレンス値に対するマトリクス変換（式（１２））の誤差の２乗（Ｅ2）は、次の式（１３）で表され、これに基づいてリファレンス値に対するマトリクス変換の誤差を最小にするマトリクスＡが算出される。
【０１０１】
【数１１】
また、誤差最小二乗法で使用する色空間を、XYZ色空間以外のものに変更するようにしてもよい。例えば、人間の知覚に対して均等なLab，Luv，Lch色空間（知覚均等色空間）に変換した後に同様の演算を行うことにより、知覚的な誤差の少ない色の再現を可能にするリニアマトリクスを算出することができる。なお、これらの色空間の値は、XYZ値から非線形な変換によって算出されるため、誤差最小二乗法においても非線形な計算アルゴリズムが使用される。
【０１０３】
上述したような演算により、例えば、図１５に示す分光感度特性を有するフィルタに対するマトリクス係数として式（１４）で表されるものが算出される。
【０１０４】
【数１２】
ステップＳ４５において、リニアマトリクスが決定される。例えば、作成される最終のRGB画像データが次の式（１５）で表されるとした場合、リニアマトリクス（LinearＭ）は、以下のようにして算出される。
【０１０６】
RGBout＝［Ｒ０,Ｇ０,Ｂ０］ｔ・・・(15)
すなわち、照明光がＤ６５である場合、sRGB色空間をXYZ色空間に変換する変換式はITU-R709.BTマトリクスを含む式（１６）で表され、そのITU-R709.BTマトリクスの逆マトリクスにより式（１７）が算出される。
【０１０７】
【数１３】
式（１２）のマトリクス変換式、式（１５）および式（１７）のITU-R709.BTマトリクスの逆マトリクスにより、式（１８）が算出される。式（１８）の右辺には、ITU-R709.BTマトリクスの逆マトリクスと、上述したマトリクスＡを乗算した値としてのリニアマトリクスが含まれる。
【０１０９】
【数１４】
すなわち、３×４のリニアマトリクス（LinearＭ）は式（１９−１）により表され、例えば、式（１４）のマトリクス係数が用いられる、図１５に示す分光分布特性を有する４色カラーフィルタに対するリニアマトリクスは、式（１９−２）により表される。
【０１１１】
【数１５】
以上のようにして算出されたリニアマトリクスが、図１０に示すリニアマトリクス処理部９５に与えられる。これにより、輝度を線形変換により表すことができる信号（Ｒ，Ｇ１，Ｇ２，Ｂ）に対してマトリクス処理を行うことができるので、図２に示す信号処理部１１における処理のように、ガンマ処理を施した後に得られる信号に対してマトリクス処理を行う場合に較べて、色彩工学的に、より忠実な色を再現させることができる。
【０１１３】
次に、図１１に示すステップＳ６において行われる評価について説明する。
【０１１４】
以上のようにして作成された、例えば、図１５に示す分光感度特性を有する４色カラーフィルタが設けられた撮像装置の色再現性と、図１に示す３色カラーフィルタが設けられる画像処理装置の色再現性を比較した場合、以下のような差が現れる。
【０１１５】
例えば、マクベスチャートを２種類の画像入力装置（４色カラーフィルタが設けられる撮像装置と、３色カラーフィルタが設けられる撮像装置）で撮像したときの出力値とリファレンス値とのLab色空間における色差が次の式（２０）によりそれぞれ算出される。
【０１１６】
【数１６】
Ｌ1−Ｌ2は２つの試料の明度差であり、ａ1−ａ2、ｂ1−ｂ2は２つの試料の色相・彩度の成分差を表している。
【０１１８】
図１７は、式（２０）による算出結果を示す図である。図１７に示すように、３色カラーフィルタが設けられる撮像装置の場合は色差が「３．３２」であるのに対して、４色カラーフィルタが設けられる撮像装置の場合は「１．３９」であり、「色の見え方」は４色カラーフィルタが設けられた撮像装置の方が優れている（色差が小さい）。
【０１１９】
図１８においては、物体Ｒ１のＲ値が「４９．４」、Ｇ値が「６４．１」、Ｂ値が「１４９．５」とされ、物体Ｒ２のＲ値が「６６．０」、Ｇ値が「６３．７」、Ｂ値が「１５５．６」とされている。従って、４色カラーフィルタにおいては、物体Ｒ１と物体Ｒ２のRGB値はそれぞれ異なる値となり、眼で見た場合と同様に、それぞれの物体の色が識別されている。すなわち、４種類の色を識別できるフィルタを設けることにより、「色の判別性」が改善されている。
【０１２０】
以上においては、４色カラーフィルタ６１は、図６に示すように、Ｇ１フィルタの左右にＢフィルタが設けられるとともに、Ｇ２フィルタの左右にＲフィルタが設けられるような配列により構成されるとしたが、図１９に示すような配列により構成されるようにしてもよい。図１９に示す４色カラーフィルタ６１においては、Ｇ１フィルタの左右にＲフィルタが設けられるとともに、Ｇ２フィルタの左右にＢフィルタが設けられている。４色カラーフィルタ６１をこのように構成することによっても、図６に示すものと同様に、「色の判別性」、「色の再現性」および「ノイズの低減性」を向上させることができる。
【０１２１】
ところで、色差（ΔＥ値）を最小にするように、リニアマトリクス係数Ｍを決定した場合、もしも、画像センサ４５の前段部に形成されているカラーフィルタの分光感度が図４に示すように重なっているときには、式（２１）のようにリニアマトリクス係数Ｍの差が大きくなってしまう。
【０１２２】
【数１７】
撮像素子の出力信号には、微小なノイズ成分が含まれているため、このようなリニアマトリクス係数Ｍを使用して色分離処理を行なうと、微小なノイズ成分も増幅されてしまう。したがって、リニアマトリクス係数Ｍの差が大きくならないように、色再現性よりもノイズ低減性を考慮する必要が生じる。しかし、当該撮像装置により実際に被写体を撮像するときに、撮像するシーンや環境により、色再現性よりもノイズ低減性を重視してリニアマトリクス係数Ｍを決定し、適応的にリニアマトリクス処理を行なった方が画質の向上を図れる場合や、逆に、ノイズ低減性よりも色再現性を重視してリニアマトリクス係数Ｍを決定し、適応的にリニアマトリクス処理を行なった方が画質の向上を図れる場合がある。また、撮像装置の使用の用途は、ユーザごとに異なるため、リニアマトリクス係数Ｍの決定をユーザが任意に行ないたい場合がある。
【０１２４】
そこで、本願発明に係る撮像装置では、上述したような問題点を解決するために、図２０に示すフローチャートにしたがってリニアマトリクス係数Ｍを決定する。
【０１２５】
まず、使用するチャート及び照明光の決定を行ない（ステップＳ５０）、つぎに色再現性指標ΔＥ（Ｍ）の定義を行ない（ステップＳ５１）、つぎにノイズ低減性指標σＮ（Ｍ）の定義を行ない（ステップＳ５２）、つぎに評価指数ＥＥＶ（Ｍ）（Error Evaluation Value）の定義を行ない（ステップＳ５３）、そして、評価指数ＥＥＶ（Ｍ）に基づきリニアマトリクス係数Ｍを決定する（ステップＳ５４）。なお、ステップ５では、評価指数ＥＥＶ（Ｍ）の係数を撮像条件等により適応的に変化させ、対応するリニアマトリクス係数Ｍを決定する。以下に、各ステップの詳細について説明する。
【０１２６】
使用するチャート及び照明光の決定（ステップＳ５０）について説明する。リニアマトリクス係数Ｍを決定するためには、カラーチャートと、当該カラーチャートを照明する光源を決める必要がある。カラーチャートとしては、Macbeth Color Checker, Digital Camera Color Checker, IT8.7等の複数の均一な色平面を持つカラーパッチから成る様々な反射チャート又は透過チャートが考えられる。照明光としては、撮像装置がよく使用される環境の光に近い分光感度を持つもの（例えばＤ５５や光源等）が考えられる。なお、撮像装置は、ユーザの使用目的により様々な光源下で用いられることが考えられるので、照明光としては、よく使用される環境の光のみに限定されるものではない。
【０１２７】
つぎに、色再現性指標ΔＥ（Ｍ）の定義（ステップＳ５１）について説明する。色再現性は、目標色と撮像装置のリニアマトリクス処理部９５でリニアマトリクス処理が行なわれた信号値が示す色（以下、出力色という。）との差で定義される。なお、色の値としては、RGB値、YCｂCr値又はXYZ値等様々考えられるが、人間の眼の見えに対して知覚が均等な色空間の値（L*a*b*値、L*u*v*値等）を用いて定義する方が、より色の差を正確に示すことが可能となる。例えば、カラーチャート中のｋ番目のカラーパッチの目標色をLabref_k(L*ref_k, a*ref_k, b*ref_k)とし、撮像装置の出力色をL*a*b*shot_k(L*shot_k, a*shot_k, b*shot_k)とすると、このパッチの色差ΔＥｋは式（２２）のようになる。
【０１２８】
【数１８】
また、色再現性指標ΔＥ（Ｍ）としては、カラーチャートの各パッチの平均ΔＥ値や、各パッチに重み付けを行ってある特定色の色再現性を重視した値等が考えられる。
【０１３０】
【数１９】
ただし、ｗｋは、各パッチに対する重み付け係数を示し、TotalPatchNumは、カラーパッチの総数を示す。
【０１３２】
また、実際には撮像装置の出力色には、リニアマトリクス処理が施されているので、L*a*b*shot_kは、リニアマトリクス係数Ｍの関数値となり、したがってΔＥk、ΔＥもＭの関数値となる。
【０１３３】
【数２０】
つぎに、ノイズ低減性指標σＮ（Ｍ）の定義（ステップＳ５２）について説明する。ノイズ低減性指標σＮ（Ｍ）は、撮像装置のリニアマトリクス処理部９５でリニアマトリクス処理が行なわれた信号値の標準偏差で定義する。信号値としては、RGB値、YCｂCr値又はXYZ値等様々考えられるが、人間の眼の見えに対して知覚が均等な色空間の値（L*a*b*値、L*u*v*値）等を用いて信号値を定義する方が、より人が感じるノイズ感と相関のあるノイズ値σＮｋを得ることが可能となる。なお、ノイズ値は、信号値の色空間に対応して、色空間の各成分の標準偏差となる。例えば、RGB空間の場合には、ノイズ値は、σR,σG,σBとなり、XYZ空間の場合には、ノイズ値は、σX,σY,σZとなる。ノイズ低減性指標σＮ（Ｍ）の定義では、これらのノイズ値を使って１つのノイズ指標を決定する。例えば、あるカラーパッチを照明光下で撮像した場合、L*a*b*空間のノイズ値σＮｋは、明度ノイズをσL*k、色ノイズをσa*k,σb*kとし、これらの値から例えば式（２５）のような明度及び色ノイズを考慮したものとして定義される。
【０１３５】
【数２１】
ただし、ｗＬ＊ ｋ、ｗａ＊ ｋ、ｗｂ＊ ｋは、各標準偏差値への重み付け係数を示し、人間の眼が感じるノイズ感との相関で適切に設定されるものである。なお、ノイズ値σＮｋとしては、他の色空間の分散値を使ったもの等様々考えられる。ノイズ低減性指標σＮ（Ｍ）としては、カラーチャートの各パッチの平均σＮ値や、各パッチに重み付けを行ってある特定色のノイズ低減性を重視した値等が考えられる。
【０１３７】
【数２２】
実際には、撮像装置の信号値には、リニアマトリクス処理が施されているので、σＮｋ及びσＮはリニアマトリクス係数Ｍの関数値となる。
【０１３９】
【数２３】
つぎに、評価指数ＥＥＶ（Ｍ）の定義（ステップＳ５３）について説明する。上述したステップＳ５１及びステップＳ５２の定義により、リニアマトリクス係数Ｍの関数値である色再現性指標ΔＥ（Ｍ）及びノイズ低減性指標σＮ（Ｍ）の二値を考慮した評価指数ＥＥＶ(Error Evaluation Value)（Ｍ）を式（２８）の様に定義する。
【０１４１】
【数２４】
ただし、h, I, j, k, l, は、関数を示し、ｗｃは、色差に対する重み付け係数を示し, ｗｎは、ノイズ値に対する重み付け係数を示す。ｗｃとｗｎを変化させて、評価指数ＥＥＶ（Ｍ）が最小となるようにリニアマトリクス係数Ｍを決定することにより、色再現性及びノイズ低減性の双方を考慮したリニアマトリクス係数Ｍの決定が可能となる。なお、色再現性を重視したい場合は、ｗｃ>ｗｎとし、ノイズ低減性を重視したい場合は、ｗｃ < ｗｎと重み付けを設定すれば良い。
【０１４３】
つぎに、リニアマトリクス係数Ｍを決定する手段（ステップＳ５４）について説明する。ステップＳ５３により定義した評価指数ＥＥＶ（Ｍ）に誤差最小二乗法を適用してリニアマトリクス係数Ｍを決定する。なお、ステップＳ５４では、ｗｃ及びｗｎを適当に定めて、例えばNewton法、Steepest Descent法又はConjugate Gradient法等を回帰アルゴリズムとして使用して誤差最小二乗法を適用し、リニアマトリクス係数Ｍを決定する。
【０１４４】
また、ステップＳ５４では、撮像装置により被写体を撮像するときの環境や条件等によりステップＳ５３で定義した評価指数ＥＥＶ（Ｍ）の色差に対する重み付け係数ｗｃ及びノイズ値に対する重み付け係数 ｗｎを適応的に変化させ、誤差最小二乗法によりリニアマトリクス係数Ｍを決定する。図２１に、色再現性指標ΔＥ（Ｍ）を変化させたときのノイズ値低減性指標σＮ（Ｍ）の変化の様子を示す。図２１に示すように、一つの撮像素子においても、リニアマトリクス係数Ｍによって色再現性指標ΔＥ（Ｍ）とノイズ低減性指標σＮ（Ｍ）のトレードオフが存在する。この結果を用いて、様々な撮像環境及び条件等に応じて適応的にリニアマトリクス係数Ｍを決定する。また、あらかじめいくつかのリニアマトリクス係数Ｍのセットを用意し、ユーザが必要に応じてリニアマトリクス係数Ｍを選択して色再現性指標ΔＥ（Ｍ）及びノイズ低減性指標σＮ（Ｍ）を調節することも可能である。
【０１４５】
ここで、撮像装置が図２２に示すような特性の４色カラーフィルターよりなるCCD撮像素子を有している場合に、上述したステップＳ５０〜ステップＳ５４にしたがってリニアマトリクス係数Ｍを決定する具体例について述べる。
【０１４６】
まず、使用するチャート及び照明光を決定する（ステップＳ５０）。カラーチャートは、Macbeth Color Checker（２４色のカラーパッチを含む）を使用し、照明光は、Ｄ５５光源（ＣＩＥで定義される５５００ｋの標準昼光）を使用する。なお、それぞれの分光データは、例えば、分光放射輝度計を用いて測定されているものとする。
【０１４７】
つぎに、色再現性指標ΔＥ（Ｍ）の定義を行なう（ステップＳ５１）。目標色を人間の眼の見えとし、Lab空間での色差ΔＥを指標とする。一般的に物体色は、「物体の分光反射率」と「照明の分光エネルギー分布」と「物体を感知するセンサの分光感度分布」との積を可視光領域vis（通常400nm〜700nm）の範囲で積分した値によって式（２９）のように定義される。
物体色＝∫vis（物体の分光反射率）・（照明の分光輝度）・（物体を感知するセンサの分光感度）・・・（２９）
さらに、人間の眼の分光感度を表す等色関数を用いると、カラーチャートのｋ番目のパッチの目標色XYZref_k(Xref_k, Yref_k, Zref_k)は、式（２９）を用いて式（３０）のように表すことができる。
【０１４８】
【数２５】
また、通常は、式（３１）を使用してXYZ空間の色をL*a*b*空間の色に変換する。
【０１５０】
【数２６】
ただし、（Ｘｎ、Ｙｎ、Ｚｎ）は、完全拡散反射面（白点）のＸＹＺ値を示している。
【０１５２】
そして、目標色XYZref_kを式（３１）を用いてL*a*b*ref_k(L*ref_k, a*ref_k, b*ref_k)に変換する。
【０１５３】
また、CCD撮像素子から出力される信号値であるｒａｗデータRGBXraw_k(Rraw_k, Graw_k, Braw_k, Xraw_k)は、式（２９）を用いると式（３２）のようになる。
【０１５４】
【数２７】
撮像装置は、リニアマトリクス処理部９５でｒａｗデータRGBXraw_k(Rraw_k, Graw_k, Braw_k, Xraw_k)をリニアマトリクス係数Ｍ(ｍ０〜ｍ１１)を用いてリニアマトリクス処理をするので、リニアマトリクス処理後の撮像データは式（３３）のようになる。
【０１５６】
【数２８】
そして、上記RGBcam_kをXYZ値に変換する。なお、変換には、式（３４）に示すように、一般的に使用されている７０９系マトリクスＭ７０９を使用する。
【０１５８】
【数２９】
つぎに、式（３４）を用いて、Xcam_k、Ycam_k、及びZcam_kを求める。
【０１６０】
【数３０】
また、式（３１）を用いてL*a*b*値(L*a*b*cam_k(L*cam_k, a*cam_k, b*cam_k))に変換し、カラーチャートのｋ番目のパッチの色差ΔEkを式（３６）のように定義する。
【０１６２】
【数３１】
なお、ΔEkは、L*a*b*cam_kの値を使用しているので、リニアマトリクス係数Ｍの関数値でもあるので、ΔEk(M)と表すことができる。色再現性指標ΔE(M)は、式（３７）に示すように、各カラーパッチの色差の平均値として定義する。
【０１６４】
【数３２】
つぎに、ノイズ低減性指標σＮ（Ｍ）の定義を行なう（ステップＳ５２）。なお、ここでは、撮像装置のリニアマトリクス処理部９５によるリニアマトリクス処理後の信号値に含まれているσＬ（Ｍ）成分に基づきノイズ低減性指標σＮ（Ｍ）を定義する。
【０１６６】
一般的にCCD撮像素子自体が出力する信号CV_CCDに含まれているノイズNoise_rawは、式（３８）のように定義される。
【０１６７】
【数３３】
なお、ShotNoiseCoefとDarkNoiseは、CCD撮像素子のデバイス特性によって決まる値である。DarkNoiseは、信号値に依存しないノイズ成分(Fixed Pattern Noise等)を表し、ShotNoiseは、信号値に依存するノイズ成分（Sensor Dark Noise, Photon Shot Noise等）を表す。
【０１６９】
式（３１）を使用すると、評価する撮像装置のｋ番目のカラーパッチのｒａｗデータに含まれているノイズ成分は式（３９）のように定義される。
【０１７０】
【数３４】
また、文献（P.D.Burns and R.S.Berns, "Error Propagation Analysis in Color Measurement and Imaging", Color Research and Application, 1997）には、以下に示すようなノイズ伝搬理論が記載されている。
【０１７２】
所定の入力信号Ｘ'＝［ｘ１，ｘ２，・・・，ｘｎ］が、（ｍ×ｎ）マトリクスＡによって、Ｙ'＝［ｙ１，ｙ２，・・・，ｙｍ］に線形変換されるとすると、式（４０）のように表される。
Ｙ'＝Ａ・Ｘ'・・・（４０）
いま、入力信号Ｘ'の分散共分散行列Σｘが、式（４１）に示すようになると、対角成分が入力信号のNoise分散値となる。
【０１７３】
【数３５】
もし、入力信号値同士の相関がなければ、行列成分中の共分散成分（すなわち非対角成分）は０となる。このとき、出力信号Ｙの分散共分散行列Σｙは、式（４２）のように定義される。
Σｙ＝Ａ・ΣｘＡ'・・・（４２）
なお、式（４２）は、線形変換で変換することができる色空間同士での、Noise分散値の伝搬理論式となる。
【０１７５】
また、撮像装置のリニアマトリクス処理部９５によるリニアマトリクス処理後の信号RGBcam_kをL*a*b*cam_kに変換するためには、RGB空間からXYZ空間への変換処理（以下、RGB→XYZ変換処理という。）を行ない、その後、XYZ空間からL*a*b*空間への変換処理（以下、XYZ→L*a*b*変換処理という。）を行なう必要がある。RGB→XYZ変換処理は、式（３４）に示す７０９系マトリクスＭ７０９を使用して線形変換を行なうことができるが、XYZ→L*a*b*変換処理は、式（３１）に示すように非線形変換を行なう必要がある。しかし、XYZ→L*a*b*変換処理は、ノイズ量が微小量であることから、JacobianマトリクスJL*a*b*_kを使用して線形近似することができ、RGB→XYZ変換処理と同様に線形変換により行なうことができる。
【０１７６】
リニアマトリクス後の信号値をXYZ値に変換した値をXYZcam_k(Xcam_k,Ycam_k,Zcam_k)とすると、式（４３）のように表すことができる。
JL*a*b*_k＝J0D(XYZcam_k)・・・（４３）
ただし、Ｊ０は、式（４４）とする。
【０１７７】
【数３６】
したがって、CCD撮像素子から出力されるｒａｗデータをLab値に線形変換するための近似マトリクスMtotal_kは式（４６）のようになる。
【０１７９】
【数３７】
式（４６）のマトリクスと式（４２）に示したノイズ伝搬理論の式を適用すると、ｋ番目のカラーパッチにおける明度ノイズσLkは、式（４７）で計算することができる。
【０１８１】
【数３８】
したがって、式（４７）より式（４８）を導き出すことができる。
【０１８３】
【数３９】
なお、式（４８）は、リニアマトリクス係数Ｍの関数なのでσL*k(M)で表すことができる。ノイズ低減性指標σN(M)は、カラーパッチの各明度ノイズの平均値なので式（４９）のように定義することができる。
【０１８５】
【数４０】
つぎに、上述のように定義した色再現性指標ΔE(M)及びノイズ低減性指標σN(M)を考慮した評価指数EEV(M)を式（５０）のように定義する（ステップＳ５３）。なお、EEVは、Error Evaluation Valueの略である。
【０１８７】
【数４１】
つぎに、式（５０）を誤差最小二乗法によって解き、リニアマトリクス係数Ｍを決定する。例えば、ｗｃ＝１としｗｎ＝２とすると評価指数EEV(M)は、式（５１）となる。
【０１８９】
【数４２】
式（５１）を誤差最小二乗法を用いて解き、リニアマトリクス係数Ｍを式（５２）のように決定する（ステップＳ５４）。
【０１９１】
【数４３】
一方、ｗｃ＝１としｗｎ＝0とした、すなわち色再現性だけを追い込んだ色差最小マトリクスは式（５３）のようになる。
【０１９３】
【数４４】
式（５２）と、式（５３）を比較すると、明らかに式（５３）の係数の方が係数間の大小差が大きく、ノイズを増長するマトリクスであることが確認できる。
【０１９５】
ここで、撮像装置のISO感度の設定に基づきリニアマトリクス係数Ｍを適応的に決定する具体例について述べる。撮像装置では、ISO感度の設定に基づいてCCD撮像素子から入力される信号（以下、入力信号という。）を増幅又は減衰させる。撮像装置のISO感度の設定をISO100からISO200に変更したときには、入力信号は、例えば、ISO100のときの２倍に増幅されて撮像装置に入力される。しかし、撮像装置では、ISO感度の設定状態にかかわらず全ての入力信号に対して同一のリニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なうため、ISO感度の設定が高い場合には、入力信号の増幅に伴い、入力信号に含まれているノイズ成分も一緒に増幅されてしまう。したがって、ISO感度の設定を高くして高解像度の画像を得ようとしても、増幅されたノイズ成分が含まれている画像が生成されてしまう。
【０１９６】
そこで、本願発明に係る撮像装置では、ISO感度の設定に基づき増幅又は減衰される入力信号に含まれているノイズ成分の低減性を考慮してリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なう。例えば、表１に示すように、ノイズ低減性に対する重み係数（ｗｎ）をISO感度の設定に応じて変化できるようにし、ｗｃ及びｗｎを式（５０）に代入してISO感度の設定ごとのリニアマトリクス係数Ｍを決定する。したがって、撮像装置は、ISO感度の設定状態に基づいて決定されたリニアマトリクス係数Ｍを用いてリニアマトリクス処理を行なうことができるので、ISO感度の設定を高くしても、それに伴いノイズ成分が増幅することがなく、高解像度の画像を得ることが可能となる。
【０１９７】
【表１】
ここで、撮像装置により被写体を撮像する環境に基づきリニアマトリクス係数Ｍを適応的に決定する例について述べる。例えば、夜景等の被写体を撮像装置により撮像した場合には、生成された画像の大半を暗い部分（暗部）が占めることがあり、その暗部にノイズが生じていると非常に目立ってしまう。このような場合には、色再現性よりもノイズ成分の低減性を考慮した方が良い。
【０１９９】
そこで、本発明に係る撮像装置では、被写体を撮像するシーンに基づきノイズ低減性と色再現性を考慮してリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なう。例えば、図２３のヒストグラムに示すように、撮像装置の輝度ダイナミックレンジの半分以下の領域に、生成した画像の輝度成分の７０％以上が含まれているようなときには、ノイズ低減性を重視してリニアマトリクス係数Ｍを決定し、それ以外のときには色再現性及びノイズ低減性を考慮してリニアマトリクス係数Ｍを決定する。具体的には、表２に示すように、ノイズ低減性に対する重み係数（ｗｎ）を撮像シーンに応じて変化できるようにし、ｗｃ及びｗｎを式（５０）に代入して、撮像シーンごとにリニアマトリクス係数Ｍを決定する。したがって、撮像装置は、撮像シーンに基づいて決定されたリニアマトリクス係数Ｍを用いてリニアマトリクス処理を行なうことができるので、生成された画像の大半を暗部が占めても、ノイズ成分を目立たなくすることが可能である。
【０２００】
【表２】
ここで、撮像装置を使用するユーザの要求に基づきリニアマトリクス係数Ｍを適応的に決定する具体例について述べる。撮像装置で被写体を撮影して生成した画像は、ユーザの使用用途によって、色再現性よりもノイズが少ないことが要求されることが多々ある。使用用途は撮像装置の製作会社には分からない事で、ユーザだけが知っている事実である。
【０２０２】
そこで、本願発明に係る撮像装置では、ユーザが意図的する条件に基づきリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なう。例えば、表３に示すように、ノイズ低減性に対する重み係数（ｗｎ）をノイズ量調節変数に応じて変化できるようにし、ｗｎ及びｗｃを式（５０）に代入して、ノイズ量調整変数ごとにリニアマトリクス係数Ｍを決定し、決定したリニアマトリクス係数Ｍを記憶しておき、ユーザがノイズ量調節変数を当該撮像装置のインターフェースを通して変化させたときに、所定のリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なう。したがって、撮像装置は、ユーザの要求に応じて決定されたリニアマトリクス係数Ｍを用いてリニアマトリクス処理を行なうことができるので、ユーザの使用状況に応じたノイズ量調節を行なうことが可能である。
【０２０３】
【表３】
なお、表３では、ノイズ量調節変数の数値が大きくなればなるほどノイズ量が減ることを意味している。また、撮像前にノイズ量調節変数の設定を行なわず、撮像装置の撮像開始ボタンを押圧したときに、いくつかのリニアマトリクス係数Ｍを使用したリニアマトリクス処理後の画像が複数個メモリに保存されるようにしても良い。
【０２０５】
このように構成された撮像装置は、図１２に示すフローチャートに従って画像センサ４５の前段部に形成される４色カラーフィルタ６１が決定されるので、「色再現性」のうちの「色の判別性」を改善させることができ、また、信号処理部７１で輝度を線形変換により表すことができる信号（Ｒ，Ｇ１，Ｇ２，Ｂ）に対してマトリクス処理を行うので、図２に示す信号処理部１１における処理のように、ガンマ処理を施した後に得られる信号に対してマトリクス処理を行う場合に較べて、色彩工学的に、より忠実な色を再現することができ、さらに、信号処理部７１のリニアマトリクス処理部９５で行なうリニアマトリクス処理において、リニアマトリクス係数Ｍの決定を撮像条件等に応じて決定するので、撮像画像の色再現性及びノイズの低減性を向上させることができる。
【０２０６】
【発明の効果】
以上詳細に説明したように、本発明に係る撮像装置は、人間の眼の見え方に対して忠実な色の再現を表す色再現値と人間が感じるノイズ感を表すノイズ値を調整する調整部と、上記調整部の調整に基づきマトリクス係数を決定するマトリクス係数決定部と、上記マトリクス係数に基づき当該撮像装置に備えられている撮像素子部で撮像した画像にマトリクス変換処理を行なうマトリクス変換処理部とを備えるので、撮像環境及び条件に応じて適応的にリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なうことができる。
【０２０７】
また、本発明に係る撮像方法は、分光特性の異なるカラーフィルタによりなっており、人間の眼の見え方に対して忠実な色の再現を表す色再現値と人間が感じるノイズ感を表すノイズ値を調整する第１のステップと、上記第１のステップの調整に基づきマトリクス係数を決定する第２のステップと、上記マトリクス係数に基づき被写体を撮像する撮像素子部で撮像した画像にマトリクス変換処理を行なう第３のステップとを有するので、撮像環境及び条件に応じて適応的にリニアマトリクス係数Ｍを決定し、当該リニアマトリクス係数Ｍを使用してリニアマトリクス処理を行なうことができる。
【図面の簡単な説明】
【図１】従来のRGBカラーフィルタの例を示す図である。
【図２】従来の撮像装置に設けられる信号処理部の構成例を示すブロック図である。
【図３】分光感度特性の例を示す図である。
【図４】分光感度特性の他の例を示す図である。
【図５】本発明を適用した撮像装置の構成例を示すブロック図である。
【図６】本発明を適用した撮像装置に設けられる４色カラーフィルタの例を示す図である。
【図７】視感度曲線の例を示す図である。
【図８】評価係数の特徴を示す図である。
【図９】本発明を適用した撮像装置が備えるカメラシステムLSIの構成例を示すブロック図である。
【図１０】図９の信号処理部の構成例を示すブロック図である。
【図１１】画像処理装置の作成処理を説明するフローチャートである。
【図１２】図１１のステップＳ１の４色カラーフィルタ決定処理の詳細を説明するフローチャートである。
【図１３】仮想曲線の例を示す図である。
【図１４】フィルタ毎のUMG値の例を示す図である。
【図１５】４色カラーフィルタの分光感度特性の例を示す図である。
【図１６】図１１のステップＳ２のリニアマトリクス決定処理の詳細を説明するフローチャートである。
【図１７】色差の評価結果の例を示す図である。
【図１８】４色カラーフィルタによる所定の物体の色度を示す図である。
【図１９】本発明を適用した撮像装置に設けられる４色カラーフィルタの他の例を示す図である。
【図２０】適応的にリニアマトリクス係数Ｍの決定を示すフローチャートである。
【図２１】色再現性指標を変化させたときのノイズ低減性指標の変化の様子を示す図である。
【図２２】４色カラーフィルタの分光感度特性の例を示す図である。
【図２３】画像のヒストグラムを示す図である。
【符号の説明】
４５ 画像センサ、６１ ４色カラーフィルタ、７１ 信号処理部、９１ オフセット補正処理部、９２ ホワイトバランス補正処理部、９３ 垂直方向同時化処理部、９４ 信号生成処理部、９５ リニアマトリクス処理部、９６−１，９６−２，９６−３ ガンマ補正処理部、９７ 輝度信号生成処理部、９８ 色差信号生成処理部、９９ 帯域制限間引き処理部
Claims (1)
- 分光特性の異なるカラーフィルタによりなっており、被写体を撮像する撮像素子部を有し、該撮像素子部で撮像した画像にマトリクス変換処理を行って色補正を行う撮像装置において、
色再現性指標ΔＥ（Ｍ）とノイズ低減性指標σＮ（Ｍ）とから評価指数ＥＥＶ（Ｍ）を以下の条件により導出し、導出された上記評価指数ＥＥＶ（Ｍ）に基づいてマトリクス係数を決定し、決定された上記マトリクス係数に基づいて上記撮像素子部で撮像した画像にマトリクス変換処理を行うマトリクス変換処理手段を備えることを特徴とする撮像装置。
（条件）
ＥＥＶ（Ｍ）＝√｛（ｗｃ・ΔＥ（Ｍ）） ２ ＋（ｗｎ・σＮ（Ｍ）） ２ ｝
但し、ＥＥＶ（Ｍ）：評価指数、ΔＥ（Ｍ）：色再現性指標、σＮ（Ｍ）：ノイズ低減性指標、ｗｃ：色再現性に対する重み係数、ｗｎ：ノイズ低減性に対する重み係数
Priority Applications (7)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002375423A JP3861808B2 (ja) | 2002-12-25 | 2002-12-25 | 撮像装置 |
| CN2003801076614A CN1732694B (zh) | 2002-12-25 | 2003-12-02 | 图像拾取装置和方法 |
| PCT/JP2003/015437 WO2004059987A1 (ja) | 2002-12-25 | 2003-12-02 | 撮像装置及び方法 |
| KR1020057012037A KR20050088335A (ko) | 2002-12-25 | 2003-12-02 | 촬상 장치 및 방법 |
| EP03813977A EP1578139A4 (en) | 2002-12-25 | 2003-12-02 | PICTURE DEVICE AND METHOD |
| US10/540,404 US7489346B2 (en) | 2002-12-25 | 2003-12-02 | Image pick-up device and image pick-up method adapted with image pick-up sensitivity |
| TW092134391A TWI225940B (en) | 2002-12-25 | 2003-12-05 | Image taking apparatus and method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2002375423A JP3861808B2 (ja) | 2002-12-25 | 2002-12-25 | 撮像装置 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2004208079A JP2004208079A (ja) | 2004-07-22 |
| JP2004208079A5 JP2004208079A5 (ja) | 2005-03-17 |
| JP3861808B2 true JP3861808B2 (ja) | 2006-12-27 |
Family
ID=32677336
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002375423A Expired - Lifetime JP3861808B2 (ja) | 2002-12-25 | 2002-12-25 | 撮像装置 |
Country Status (7)
| Country | Link |
|---|---|
| US (1) | US7489346B2 (ja) |
| EP (1) | EP1578139A4 (ja) |
| JP (1) | JP3861808B2 (ja) |
| KR (1) | KR20050088335A (ja) |
| CN (1) | CN1732694B (ja) |
| TW (1) | TWI225940B (ja) |
| WO (1) | WO2004059987A1 (ja) |
Families Citing this family (38)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060164533A1 (en) * | 2002-08-27 | 2006-07-27 | E-Phocus, Inc | Electronic image sensor |
| JP4501634B2 (ja) | 2004-10-29 | 2010-07-14 | 富士フイルム株式会社 | マトリクス係数決定方法及び画像入力装置 |
| EP1833264B1 (en) | 2004-12-17 | 2012-11-07 | Nikon Corporation | Image processing method |
| JP4859502B2 (ja) * | 2005-03-30 | 2012-01-25 | Ｈｏｙａ株式会社 | 撮像装置 |
| US7864235B2 (en) | 2005-03-30 | 2011-01-04 | Hoya Corporation | Imaging device and imaging method including generation of primary color signals |
| JP4293174B2 (ja) * | 2005-09-28 | 2009-07-08 | ソニー株式会社 | 撮像装置および画像処理装置 |
| JP2007103401A (ja) * | 2005-09-30 | 2007-04-19 | Matsushita Electric Ind Co Ltd | 撮像装置及び画像処理装置 |
| KR100721543B1 (ko) * | 2005-10-11 | 2007-05-23 | (주) 넥스트칩 | 통계적 정보를 이용하여 노이즈를 제거하는 영상 처리 방법및 시스템 |
| KR101147084B1 (ko) | 2005-12-20 | 2012-05-17 | 엘지디스플레이 주식회사 | 액정 표시장치의 구동장치 및 구동방법 |
| JP4983093B2 (ja) * | 2006-05-15 | 2012-07-25 | ソニー株式会社 | 撮像装置および方法 |
| TWI400667B (zh) * | 2006-09-06 | 2013-07-01 | Magic Pixel Inc | 色彩補插方法及應用其之影像處理裝置 |
| JP4874752B2 (ja) * | 2006-09-27 | 2012-02-15 | Ｈｏｙａ株式会社 | デジタルカメラ |
| TWI318536B (en) * | 2006-09-29 | 2009-12-11 | Sony Taiwan Ltd | A method of color matching applied on the image apparatus |
| US7847830B2 (en) * | 2006-11-21 | 2010-12-07 | Sony Ericsson Mobile Communications Ab | System and method for camera metering based on flesh tone detection |
| US8532428B2 (en) * | 2007-02-27 | 2013-09-10 | Nec Corporation | Noise reducing apparatus, noise reducing method, and noise reducing program for improving image quality |
| US8108211B2 (en) | 2007-03-29 | 2012-01-31 | Sony Corporation | Method of and apparatus for analyzing noise in a signal processing system |
| US8711249B2 (en) | 2007-03-29 | 2014-04-29 | Sony Corporation | Method of and apparatus for image denoising |
| US9285670B2 (en) * | 2007-09-14 | 2016-03-15 | Capso Vision, Inc. | Data communication between capsulated camera and its external environments |
| JP5080934B2 (ja) * | 2007-10-22 | 2012-11-21 | キヤノン株式会社 | 画像処理装置及び方法、及び撮像装置 |
| JP5178226B2 (ja) * | 2008-02-08 | 2013-04-10 | オリンパス株式会社 | 画像処理装置および画像処理プログラム |
| JP4735994B2 (ja) * | 2008-08-27 | 2011-07-27 | ソニー株式会社 | 撮像装置及び方法、プログラム、並びに記録媒体 |
| US20110075935A1 (en) | 2009-09-25 | 2011-03-31 | Sony Corporation | Method to measure local image similarity based on the l1 distance measure |
| KR101129220B1 (ko) * | 2009-11-03 | 2012-03-26 | 중앙대학교 산학협력단 | 레인지 영상의 노이즈 제거장치 및 방법 |
| JP5460276B2 (ja) * | 2009-12-04 | 2014-04-02 | キヤノン株式会社 | 撮像装置及び撮像システム |
| JP5308375B2 (ja) * | 2010-02-25 | 2013-10-09 | シャープ株式会社 | 信号処理装置、固体撮像装置、電子情報機器、信号処理方法、制御プログラムおよび記録媒体 |
| JP2012010276A (ja) * | 2010-06-28 | 2012-01-12 | Sony Corp | 画像処理装置、画像処理方法及び画像処理プログラム |
| US8731281B2 (en) * | 2011-03-29 | 2014-05-20 | Sony Corporation | Wavelet transform on incomplete image data and its applications in image processing |
| JP2012253727A (ja) * | 2011-06-07 | 2012-12-20 | Toshiba Corp | 固体撮像装置及びカメラモジュール |
| JP5954661B2 (ja) * | 2011-08-26 | 2016-07-20 | パナソニックＩｐマネジメント株式会社 | 撮像素子、及び撮像装置 |
| US8780225B2 (en) * | 2011-10-12 | 2014-07-15 | Apple Inc. | Use of noise-optimized selection criteria to calculate scene white points |
| US8654210B2 (en) | 2011-10-13 | 2014-02-18 | Canon Kabushiki Kaisha | Adaptive color imaging |
| JP5895569B2 (ja) * | 2012-02-08 | 2016-03-30 | ソニー株式会社 | 情報処理装置、情報処理方法およびコンピュータプログラム |
| US9531922B2 (en) * | 2013-10-30 | 2016-12-27 | Canon Kabushiki Kaisha | Image processing apparatus, image processing method, and computer readable storage medium for improving sharpness |
| JP2015180045A (ja) * | 2014-02-26 | 2015-10-08 | キヤノン株式会社 | 画像処理装置、画像処理方法及びプログラム |
| JP6013382B2 (ja) * | 2014-02-27 | 2016-10-25 | 富士フイルム株式会社 | 内視鏡システム及びその作動方法 |
| CN104902153B (zh) * | 2015-05-25 | 2017-11-28 | 北京空间机电研究所 | 一种多光谱相机色彩校正方法 |
| WO2020221731A1 (en) * | 2019-04-30 | 2020-11-05 | Signify Holding B.V. | Luminance distribution determination |
| EP3896967A1 (en) | 2020-04-17 | 2021-10-20 | Leica Microsystems CMS GmbH | Digital imaging device and method for generating a digital color image |
Family Cites Families (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0815338B2 (ja) * | 1984-08-13 | 1996-02-14 | 株式会社日立製作所 | カラービデオカメラ |
| JP3263924B2 (ja) | 1990-09-14 | 2002-03-11 | ソニー株式会社 | カラー撮像装置 |
| US5668596A (en) * | 1996-02-29 | 1997-09-16 | Eastman Kodak Company | Digital imaging device optimized for color performance |
| JP3750830B2 (ja) * | 1996-08-30 | 2006-03-01 | ソニー株式会社 | 撮像装置における色補正装置 |
| JP3880276B2 (ja) * | 2000-03-03 | 2007-02-14 | 三菱電機株式会社 | 撮像装置 |
| JP4064038B2 (ja) * | 2000-06-09 | 2008-03-19 | 富士フイルム株式会社 | 固体撮像素子を用いた画像取得装置および画像取得方法並びにその方法を実行するためのプログラムを記録した記録媒体 |
| JP4534340B2 (ja) | 2000-10-31 | 2010-09-01 | ソニー株式会社 | 色再現補正装置 |
| US6822657B2 (en) * | 2001-08-10 | 2004-11-23 | Agilent Technologies, Inc. | Method and apparatus for improving image quality in digital cameras |
| JP2003235050A (ja) * | 2002-02-12 | 2003-08-22 | Nikon Corp | 画像処理装置、画像処理プログラム、および画像処理方法 |
| US7053935B2 (en) * | 2002-02-21 | 2006-05-30 | Eastman Kodak Company | Apparatus and method for accurate electronic color capture and reproduction |
| JP2003284084A (ja) * | 2002-03-20 | 2003-10-03 | Sony Corp | 画像処理装置および方法、並びに画像処理装置の製造方法 |
-
2002
- 2002-12-25 JP JP2002375423A patent/JP3861808B2/ja not_active Expired - Lifetime
-
2003
- 2003-12-02 KR KR1020057012037A patent/KR20050088335A/ko not_active Application Discontinuation
- 2003-12-02 EP EP03813977A patent/EP1578139A4/en not_active Ceased
- 2003-12-02 WO PCT/JP2003/015437 patent/WO2004059987A1/ja active Application Filing
- 2003-12-02 CN CN2003801076614A patent/CN1732694B/zh not_active Expired - Fee Related
- 2003-12-02 US US10/540,404 patent/US7489346B2/en not_active Expired - Fee Related
- 2003-12-05 TW TW092134391A patent/TWI225940B/zh not_active IP Right Cessation
Also Published As
| Publication number | Publication date |
|---|---|
| JP2004208079A (ja) | 2004-07-22 |
| TWI225940B (en) | 2005-01-01 |
| EP1578139A4 (en) | 2006-05-03 |
| CN1732694A (zh) | 2006-02-08 |
| EP1578139A1 (en) | 2005-09-21 |
| US7489346B2 (en) | 2009-02-10 |
| KR20050088335A (ko) | 2005-09-05 |
| US20060082665A1 (en) | 2006-04-20 |
| WO2004059987A1 (ja) | 2004-07-15 |
| CN1732694B (zh) | 2012-11-14 |
| TW200424569A (en) | 2004-11-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3861808B2 (ja) | 撮像装置 | |
| EP1487219A1 (en) | Image processing device, image processing method, and image processing device manufacturing method | |
| TWI255645B (en) | Image pickup device | |
| JP4265546B2 (ja) | 撮像装置、画像処理装置および画像処理方法 | |
| Ramanath et al. | Color image processing pipeline | |
| US8803994B2 (en) | Adaptive spatial sampling using an imaging assembly having a tunable spectral response | |
| WO2010053029A1 (ja) | 画像入力装置 | |
| JP2007259344A (ja) | 撮像装置および画像処理方法 | |
| US7053935B2 (en) | Apparatus and method for accurate electronic color capture and reproduction | |
| JP2009105576A (ja) | 画像処理装置及び方法、及び撮像装置 | |
| US9060159B2 (en) | Image processing device and method, and imaging device | |
| JP4874752B2 (ja) | デジタルカメラ | |
| US20080225135A1 (en) | Imaging Device Element | |
| WO2007148576A1 (ja) | 撮像システム及び撮像プログラム | |
| JP2005278213A (ja) | 製造方法 | |
| JP4285587B2 (ja) | 撮像素子 | |
| JP3863773B2 (ja) | 画像撮影方法および装置 | |
| JP2007097202A (ja) | 画像処理装置 | |
| JP2014042138A (ja) | 画像処理装置、コンピュータプログラム、およびデジタルカメラ | |
| JP4298595B2 (ja) | 撮像装置及びその信号処理方法 | |
| JP2007267404A (ja) | 画像処理装置の製造方法 | |
| JP2005323141A (ja) | 撮像装置、カメラ、及び信号処理方法 | |
| JP7203441B2 (ja) | 画像生成装置及び撮像装置 | |
| JP5104452B2 (ja) | 画像処理装置、画像処理方法および画像処理プログラム | |
| KR20050006424A (ko) | 색신호 처리방법 및 그 방법을 사용하는 촬상장치 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20040419 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20040419 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20051108 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060106 |
|
| 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: 20060905 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060918 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091006 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101006 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111006 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121006 Year of fee payment: 6 |