JP2001292454A - 画像処理方法および装置並びに記録媒体 - Google Patents

画像処理方法および装置並びに記録媒体

Info

Publication number
JP2001292454A
JP2001292454A JP2000104977A JP2000104977A JP2001292454A JP 2001292454 A JP2001292454 A JP 2001292454A JP 2000104977 A JP2000104977 A JP 2000104977A JP 2000104977 A JP2000104977 A JP 2000104977A JP 2001292454 A JP2001292454 A JP 2001292454A
Authority
JP
Japan
Prior art keywords
value
signal
predetermined direction
orthogonal direction
weighting
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.)
Granted
Application number
JP2000104977A
Other languages
English (en)
Other versions
JP4302854B2 (ja
Inventor
Tatsuya Aoyama
達也 青山
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.)
Fujifilm Holdings Corp
Original Assignee
Fuji Photo Film Co Ltd
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 Fuji Photo Film Co Ltd filed Critical Fuji Photo Film Co Ltd
Priority to JP2000104977A priority Critical patent/JP4302854B2/ja
Publication of JP2001292454A publication Critical patent/JP2001292454A/ja
Application granted granted Critical
Publication of JP4302854B2 publication Critical patent/JP4302854B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Color Television Image Signal Generators (AREA)
  • Editing Of Facsimile Originals (AREA)

Abstract

(57)【要約】 【課題】 単板CCDのような撮像デバイスにおいて得
られる画像データについて、色エッジ、グレーエッジに
拘わらず偽色を生じさせることなく、全画素位置におけ
る信号値を推定する。 【解決手段】 判定手段10において、信号値を算出す
る画素位置がグレーエッジにあるか色エッジにあるか否
かを判定する。第1補間手段11では所定方向について
の信号値Q1が算出され、第2補間手段12ではこれに
直交する方向について信号値Q2が算出される。変動値
算出手段13では、判定結果を表す判定情報Hに基づ
き、グレーエッジ、色エッジに応じて、その画素位置で
の所定方向および直交方向の信号値の変動量を変動値S
A,SBとして算出する。重み付け係数算出手段14は
変動値SA,SBから重み付け係数Sjを算出し、重み
付け加算手段15において、重み付け係数Sjにより信
号Q1,Q2を重み付け加算してR信号QR、B信号Q
Bを求める。

Description

【発明の詳細な説明】
【0001】
【発明の属する技術分野】本発明は、単板CCDのよう
な撮像デバイスにおいて得られた信号値を用いて、撮像
デバイスの全画素位置における信号値を推定する画像処
理方法および装置並びに画像処理方法をコンピュータに
実行させるためのプログラムを記録したコンピュータ読
取り可能な記録媒体に関するものである。
【0002】
【従来の技術】デジタルカメラに用いられるCCD等の
撮像デバイスとしては、分光感度が異なる複数種類の光
電変換素子を同一平面上に交互に配置して構成されてい
るものが知られている(以下単板CCDと称する)。こ
こで、R,G,Bのそれぞれに分光感度を有する光電変
換素子、すなわちR,G,Bの各チャンネルの光電変換
素子を交互に配置した単板CCDの場合、連続したR,
G,Bチャンネルの3個の光電変換素子の組が1つの画
素を構成することとなる。しかしながら、このような単
板CCDにおいては各画素のR,G,B信号値を同一画
素位置において得ることができないため、色ずれや偽色
が生じることがある。また、各チャンネルの光電変換素
子数は単板CCDを構成する全素子数よりも少ないた
め、高解像度の画像を得ることができない。例えばR,
G,B各チャンネルの光電変換素子を交互に配置した単
板CCDにおいては、各チャンネルの光電変換素子数は
全素子数の1/3しかないため、同一素子数のモノクロ
撮像装置に比べて解像度が1/3となってしまう。この
ため、R,G,B各チャンネルの光電変換素子が存在し
ない部分における信号値を補間処理により求める方法が
提案されているが、単に補間処理を行うのみでは、信号
値が大きく変化する部分において偽色が発生することが
ある。この場合、撮像系に光学ローパスフィルタを使用
したり、撮像信号に対してローパスフィルタによる平滑
化処理を行うことにより偽色の発生を防止することがで
きるが、この場合解像度が悪化するという問題がある。
【0003】ここで、人間の視覚特性は色よりも輝度に
対して感度が高いものである。このため、単板CCDに
おいて得られたカラー撮像信号から、各画素の輝度を表
す高周波の輝度信号と、上述した補間処理およびローパ
スフィルタによる平滑化処理によって得られた低周波の
色信号とを生成し、輝度信号および色信号を用いてカラ
ー画像信号を再構成するようにした方法が提案されてい
る(特開平10−200906号等)。この方法によれ
ば、人間の視覚特性において感度が高い輝度成分に対し
てより多くの情報が与えられることとなるため、見かけ
上解像度が高い画像を再現可能なカラー画像信号を得る
ことができる。
【0004】ところで、単板CCDとして、例えば図3
3に示すように市松状に画素が配置されたハニカム配列
のアレイ構造を有するCCDが知られている(例えば特
開平10−136391号)。なお、これを市松状の画
素配列と称することもある。また、図34に示すように
正方状に画素が配置されたベイヤー配列のアレイ構造を
有するCCDも知られている。なお、これを正方状の画
素配列と称することもある。このようなアレイ構造すな
わち画素配列を有する単板CCDにおいても上記と同様
に偽色の問題が生じる。このため、上述したベイヤー配
列のCCDにおいて得られた光量ベースの信号について
偽色を除去するために、画像中の局所的な領域において
はr,g,bの信号の比は略一定であるという仮定に基
づいて、ベイヤー配列のCCDにおける垂直方向または
水平方向のラインにおいて、隣接するラインにおけるr
信号とg信号との比をg信号に乗算することにより、そ
のラインにおけるr信号を算出するようにした方法が知
られている(特開平9−214989号)。この方法
は、具体的には、図35に示す画素配列において、g1
2画素位置におけるr信号r12を求めるには、まずr
11画素位置におけるg11信号を(g6+g16)/
2の演算により算出し、r11:g11=r12:g1
2という仮定に基づいて、r12=g12×r11/g
11の演算によりr12信号を算出するものである。
【0005】しかしながら、上記特開平10−2009
06号等に記載された方法においては、単板CCDにお
いて得られた撮像信号に対してどのようなローパスフィ
ルタにより平滑化処理を施しても、実際の画像における
高周波成分が画像中にすでに折り返されてしまっている
ため、折り返しひずみによるモアレを取り除くことがで
きず、その結果偽色を十分に除去することができない。
【0006】一方、上記特開平9−214989号に記
載された方法によれば、効果的に偽色を除去することが
できる。とくにこの方法は、画像の局所的な領域におい
ては、r:g:bの光量の比が一定であるという仮定に
基づいており、得られたRGB信号の比が光量に比例し
ているアナログ信号である場合には、ベイヤー配列のC
CDにおいて得られる画像信号における偽色を効率よく
除去することができる。しかしながら、デジタルカメラ
により得られる画像信号は、A/D変換を行って光量r
gbをデジタルのRGB信号に変換する場合に、量子化
誤差を低減するため、およびコンピュータ系のビデオ回
路へ信号を入力するために、例えばR=r0.45、R
=log(r)のように、光量に対する指数値、対数値
となるように信号値が表されていることから、r:g:
b=R:G:Bとはならないものである。このため、上
記特開平9−214989号に記載された方法は、信号
値が光量に比例するアナログ信号に対しては偽色を除去
することができるが、信号値が光量の指数値や対数値に
より表されている場合には偽色を除去することはできな
い。また、偽色はベイヤー配列の単板CCDのみならず
ハニカム配列の単板CCDにおいても発生する。
【0007】このため、撮像デバイスにおいて得られた
信号値が、光量に対する指数値や対数値により表されて
いる場合、画像の局所的な領域においてはRGB各信号
値の差は一定であるという前提に基づき、信号値を算出
する画素位置における信号値が変化する方向を考慮し
て、各画素位置における信号値を算出する方法が提案さ
れている(特願平11−212202号)。以下、この
方法について具体的に説明する。
【0008】まず、撮像デバイスであるCCDの画素配
列を図33に示すハニカム配列とし、図36に示すよう
に各画素位置に参照番号を付与する。なお、この画素配
列を正方配列と見なした場合に、信号値を有さない画素
位置を図36において*を用いて示す。また、各チャン
ネルにおいて得られる信号値は8ビットであるものと
し、紙面左上から右下に延在する方向を矢印A方向、紙
面右上から左下に延在する方向を矢印B方向とし、以下
本明細書中においては矢印A方向および矢印B方向は不
変とする。なお、ここではB24画素位置におけるR信
号(以下R(B24)信号と表記する)の算出について
説明する。まず、矢印A方向について、下記の式(1)
によりB24画素位置のR(B24)A信号を求める。
【0009】 R(B24)A=G(B24)A+((R(G33)A-G33)+(R(G15)A-G15))/2 (1) 但し、G(B24)A信号は、画素がB02,G13,
B24,G35…と並ぶライン上のG画素位置における
G信号に対して、矢印A方向に1次元Cubicスプライン
補間演算等の1次元補間演算を施して求めたB24画素
位置のG信号である。また、R(G33)A信号は、画
素が…R22,G33,R44,G55…と並ぶライン
上のR画素位置におけるR信号に対して、矢印A方向に
1次元補間演算を施して求めたG33画素位置のR信号
である。さらに、R(G15)A信号は、画素がR0
4,G15,R26,G37…と並ぶライン上のR画素
位置におけるR信号に対して、矢印A方向に1次元補間
演算を施して求めたG15画素位置のR信号である。
【0010】次に、矢印B方向について、下記の式
(2)によりB24画素位置のR(B24)B信号を求
める。
【0011】 R(B24)B=G(B24)B+((R(G13)B-G13)+(R(G35)B-G35))/2 (2) 但し、G(B24)B信号は、画素がB06,G15,
B24,G33…と並ぶライン上のG画素位置における
G信号に対して、矢印B方向に1次元補間演算を施して
求めたB24画素位置のG信号である。また、R(G1
3)B信号は、画素がR04,G13,R22,G31
…と並ぶライン上のR画素位置におけるR信号に対し
て、矢印B方向に1次元補間演算を施して求めたG13
画素位置のR信号である。さらに、R(G35)B信号
は、画素が…R26,G35,R44,G53…と並ぶ
ライン上のR画素位置におけるR信号に対して、矢印B
方向に1次元補間演算を施して求めたG35画素位置の
R信号である。
【0012】そして、下記の式(3a)、(3b)によ
り、B24画素位置における矢印A方向および矢印B方
向についての信号値の変動量を表す変動値SA,SBを
算出し、さらに、変動値SA,SBを用いて式(3c)
によりR(B24)A信号とR(B24)B信号とを重
み付け加算するための重み付け係数Sjを算出する。
【0013】 SA=|B24-G(B24)A|/(B24+G(B24)A) (3a) SB=|B24-G(B24)B|/(B24+G(B24)B) (3b) Sj=SA/(SA+SB)(SA+SB=0ならSj=0.5) (3c) このように重み付け係数Sjが求められると、下記の式
(4)により、R(B24)信号を算出する。
【0014】 R(B24)=(1-Sj)*R(B24)A+Sj*R(B24)B (4) この方法によれば、撮像デバイスにおいて得られた信号
値が指数値や対数値により表されている場合であって
も、偽色を発生させることなく全画素位置における信号
値を求めることができる。また、この方法において、信
号値が変化する方向に応じて、撮像デバイス上の所定方
向およびこれに直交する直交方向について算出された信
号値を上記式(4)に示すように重み付け加算して、各
画素位置における信号値を算出することにより、信号値
の変化方向に拘わらず偽色の発生を防止できるととも
に、アーチファクトの発生をも防止することができる。
【0015】
【発明が解決しようとする課題】しかしながら、特願平
11−212202号の方法では、信号値を推定する画
素位置が画像の色エッジ上にある場合には、偽色を発生
させることなく信号値を算出することができるが、信号
値を推定する画素がグレーエッジ上にある場合には、偽
色を十分に低減することができない。以下、この問題点
について説明する。まず、上記図36に示す画素配列に
おいて、図16に示すように信号値が得られて白(信号
値が255)および黒(信号値が0または0に近い値)
のグレーエッジが表されているものとする。なお、ここ
ではR44画素位置におけるB(R44)信号の算出に
ついて説明する。まず、矢印A方向について、下記の式
(5)によりR44画素位置のB(R44)A信号を求
める。
【0016】 B(R44)A=G(R44)A+((B(G35)A-G35)+(B(G53)A-G53))/2 =0+((2-0)+(255-255))/2=1 (5) 但し、G(R44)A=0,B(G35)A=2,B(G53)A=255 また、矢印B方向について、下記の式(6)によりR4
4画素位置のB(R44)B信号を求める。
【0017】 B(R44)B=G(R44)B+((B(G33)B-G33)+(B(G55)B-G55))/2 =103+((105-0)+(105-0))/2=208 (6) 但し、G(R44)B≒103,B(G33)B≒105,B(G55)B≒105 そして、下記の式(7)によりB(R44)A信号とB
(R44)B信号とを重み付け加算するための重み付け
係数を算出する。
【0018】 SA=|R44-G(R44)A|/(R44+G(R44)A)=|1-0|/(1+0)=1 SB=|R44-G(R44)B|/(R44+G(R44)B)=|1-103|/(1+103)≒0.98 (7) Sj=SA/(SA+SB)=1/(1+0.98)≒0.5 このように重み付け係数Sjが求められると、下記の式
(8)により、B(R44)信号を算出する。
【0019】 B(R44)=(1-Sj)*B(R44)A+Sj*B(R44)B=0.5*1+0.5*208≒105 (8) ここで、R44画素位置は黒いエッジ上にあるため、B
(R44)信号の値は0に近い値であることが好まし
い。しかしながら、上記式(7)に示すように重み付け
係数Sjを算出して式(8)によりB(R44)信号の
値を求めると、その値は105となるため、ここに偽色
が発生することとなる。
【0020】本発明は上記事情に鑑みなされたものであ
り、グレーエッジであっても偽色を生じさせることなく
信号値を推定できる画像処理方法および装置並びに画像
処理方法をコンピュータに実行させるためのプログラム
を記録したコンピュータ読取り可能な記録媒体を提供す
ることを目的とするものである。
【0021】
【課題を解決するための手段】本発明による画像処理方
法は、異なる分光感度を有する第1から第3の光電変換
素子を単一面上に配置して撮像面が形成された撮像デバ
イスであって、前記第1および前記第2の光電変換素子
を所定方向に交互に配置して第1ラインを形成し、前記
第1および前記第3の光電変換素子を前記所定方向に交
互に配置して第2ラインを形成し、前記所定方向と直交
する直交方向において前記第1および前記第2の光電変
換素子が交互に、かつ前記第1および前記第3の光電変
換素子が交互に配置されるよう、前記第1および前記第
2ラインを前記直交方向に交互に配置することにより、
第1から第3画素からなる前記撮像面が形成された撮像
デバイスにおいて得られた第1から第3信号値に基づい
て、全画素位置における前記第1から第3信号値のうち
少なくとも1つの信号値を推定する画像処理方法におい
て、信号値を推定する画素位置が、グレーエッジにある
か色エッジにあるかを判定し、前記所定方向についての
信号値の推定値である所定方向推定値および前記直交方
向についての信号値の推定値である直交方向推定値を算
出し、前記判定の結果に応じた異なる演算方法により、
前記所定方向についての信号値の変動量を表す所定方向
変動値および前記直交方向についての信号値の変動量を
表す直交方向変動値を算出し、前記所定方向変動値およ
び前記直交方向変動値に基づいて、前記所定方向推定値
および前記直交方向推定値を重み付け加算する際の重み
付け係数を設定し、該重み付け係数に基づいて、前記所
定方向推定値および前記直交方向推定値を重み付け加算
して、前記第1から第3信号値のうち少なくとも1つの
信号値を推定することを特徴とするものである。
【0022】ここで、図33に示すようなハニカム配列
の単板CCDにおいて、第1から第3画素、第1信号値
から第3信号値をそれぞれG,R,Bに対応させるとす
ると、第1ラインは図33上左上から右下方向(これを
所定方向とする)にGR画素を交互に配置したGRライ
ンとなり、第2ラインは第1ラインと同様の方向にGB
画素を交互に配置したGBラインとなる。また、所定方
向と直交する方向にもGRラインおよびGBラインが交
互に配置されている。本発明による画像処理方法は、単
板CCDのような撮像デバイスにおいて得られる信号値
が光量に対する指数値や対数値により表されている場合
には、画像上の局所的な領域において、各信号値の差分
値が等しくなるという前提に基づいて、このような撮像
デバイスにおいて得られた第1から第3信号値から、各
画素位置における全ての信号値を推定するようにしたも
のである。なお、本発明においては撮像デバイスの画素
配列はこれに限定されるものではない。
【0023】「所定方向推定値」および「直交方向推定
値」とは、所定方向および直交方向についてそれぞれ推
定した第1から第3信号値のうち少なくとも1つの信号
値のことである。
【0024】「判定の結果に応じた異なる演算方法によ
り、所定方向変動値および直交方向変動値を算出する」
とは、色エッジであると判定された場合には、例えば上
記式(3a)、(3b)に示す演算方法により変動値S
Aおよび変動値SBを所定方向変動値および直交方向変
動値として求め、グレーエッジであると判定された場合
には、例えば下記の式(9a)、(9b)に示す演算方
法により変動値SAおよび変動値SBを所定方向変動値
および直交方向変動値として求めることをいう。
【0025】 SA=|B24-G(B24)A| (9a) SB=|B24-G(B24)B| (9b) 「重み付け係数を設定する」とは、グレーエッジと判定
された場合には、偽色が発生しないように所定方向推定
値および直交方向推定値を重み付ける重み付け係数を設
定し、色エッジと判定された場合には、偽色を低減させ
つつも細かな色味を残すことができるような重み付け係
数を設定することをいう。
【0026】また、本発明による画像処理方法において
は、前記所定方向変動値と前記直交方向変動値との加算
値を算出し、該加算値に対する前記所定方向変動値また
は前記直交方向変動値の比率を算出し、該比率を前記所
定方向推定値および前記直交方向推定値を重み付け加算
する際の重み付け係数として設定することが好ましい。
【0027】さらに、本発明による画像処理方法におい
ては、前記所定方向変動値および前記直交方向変動値の
加算値に対する前記所定方向変動値または前記直交方向
変動値の比率を、前記所定方向推定値および前記直交方
向推定値を重み付け加算する際の重み付け係数として設
定するに際し、前記重み付け係数が一定であるとした場
合の前記所定方向変動値と前記直交方向変動値との関係
に基づいて、複数の重み付け係数に対する所定方向変動
値および直交方向変動値の対応関係を予め算出して記憶
した記憶手段を参照して、前記所定方向変動値および前
記直交方向変動値に基づいて、前記複数の重み付け係数
から一の重み付け係数を選択し、該選択された一の重み
付け係数を、前記所定方向推定値および前記直交方向推
定値を重み付け加算する際の重み付け係数として設定す
ることが好ましい。
【0028】また、この場合、前記記憶手段に記憶され
た前記対応関係における前記複数の重み付け係数を、前
記加算値に対する前記所定方向変動値または前記直交方
向変動値の比率が1/2(nは自然数)となるように
設定し、前記重み付け加算を、前記重み付け係数に基づ
くビットシフトにより行うことが好ましい。
【0029】ここで、加算値に対する所定方向変動値ま
たは直交方向変動値の比率が1/2 (nは自然数)と
なるように重み付け係数が設定されていると、所定方向
推定値および直交方向推定値を2進数で表した場合、こ
れらに1/2を乗算することは、これらの推定値を右
側にnビットシフトさせることと等価な演算を行ってい
ることとなる。したがって、「重み付け係数に基づくビ
ットシフト」とは、選択された一の重み付け係数の値1
/2に基づいて、所定方向推定値および直交方向推定
値を右側にnビットシフトさせることにより、重み付け
加算を行うことをいう。
【0030】また、本発明による画像処理方法において
は、前記信号値を推定する画素位置が、前記第1から第
3信号値のいずれか1つの信号値のみからなる単色平坦
部にあるか否かを判定し、前記単色平坦部にあると判定
された場合には、該単色平坦部内における前記第1から
第3信号値に基づく線形補間演算により、前記第1から
第3信号値のうち少なくとも1つの信号値を推定し、前
記単色平坦部にないと判定された場合に、前記グレーエ
ッジにあるか色エッジにあるか否かの判定、前記所定方
向推定値および前記直交方向推定値の算出、前記所定方
向変動値および前記直交方向変動値の算出、前記重み付
け係数の設定、および該重み付け係数に基づく信号値の
推定を行うことが好ましい。
【0031】また、本発明による画像処理方法において
は、前記第1から前記第3の光電変換素子が、それぞれ
G(グリーン)、B(ブルー)、R(レッド)のいずれ
かの色、あるいはY(イエロー)、G(グリーン)、C
(シアン)のいずれかの色に分光感度を有するものであ
ることが好ましい。
【0032】本発明による画像処理装置は、異なる分光
感度を有する第1から第3の光電変換素子を単一面上に
配置して撮像面が形成された撮像デバイスであって、前
記第1および前記第2の光電変換素子を所定方向に交互
に配置して第1ラインを形成し、前記第1および前記第
3の光電変換素子を前記所定方向に交互に配置して第2
ラインを形成し、前記所定方向と直交する直交方向にお
いて前記第1および前記第2の光電変換素子が交互に、
かつ前記第1および前記第3の光電変換素子が交互に配
置されるよう、前記第1および前記第2ラインを前記直
交方向に交互に配置することにより、第1から第3画素
からなる前記撮像面が形成された撮像デバイスにおいて
得られた第1から第3信号値に基づいて、全画素位置に
おける前記第1から第3信号値のうち少なくとも1つの
信号値を推定する画像処理装置において、信号値を推定
する画素位置が、グレーエッジにあるか色エッジにある
かを判定する判定手段と、前記所定方向についての信号
値の推定値である所定方向推定値および前記直交方向に
ついての信号値の推定値である直交方向推定値を算出す
る推定値算出手段と、前記判定の結果に応じた異なる演
算方法により、前記所定方向についての信号値の変動量
を表す所定方向変動値および前記直交方向についての信
号値の変動量を表す直交方向変動値を算出する変動値算
出手段と、前記所定方向変動値および前記直交方向変動
値に基づいて、前記所定方向推定値および前記直交方向
推定値を重み付け加算する際の重み付け係数を設定する
重み付け係数設定手段と、該重み付け係数に基づいて、
前記所定方向推定値および前記直交方向推定値を重み付
け加算して、前記第1から第3信号値のうち少なくとも
1つの信号値を推定する重み付け加算手段とを備えたこ
とを特徴とするものである。
【0033】なお、本発明による画像処理装置において
は、前記重み付け係数設定手段は、前記所定方向変動値
と前記直交方向変動値との加算値を算出し、該加算値に
対する前記所定方向変動値または前記直交方向変動値の
比率を算出する比率算出手段を備え、該比率を、前記所
定方向推定値および前記直交方向推定値を重み付け加算
する際の重み付け係数として設定する手段であることが
好ましい。
【0034】また、本発明による画像処理装置において
は、前記重み付け係数設定手段は、前記所定方向変動値
と前記直交方向変動値との加算値に対する前記所定方向
変動値または前記直交方向変動値の比率を、前記所定方
向推定値および前記直交方向推定値を重み付け加算する
際の重み付け係数として設定するに際し、前記重み付け
係数が一定であるとした場合の前記所定方向変動値と前
記直交方向変動値との関係に基づいて、複数の重み付け
係数に対する所定方向変動値および直交方向変動値の対
応関係を予め算出して、該対応関係を記憶した記憶手段
と、前記所定方向変動値および前記直交方向変動値に基
づいて、前記記憶手段を参照して前記複数の重み付け係
数から一の重み付け係数を選択する重み付け係数選択手
段とを備え、該選択された一の重み付け係数を、前記所
定方向推定値および前記直交方向推定値を重み付け加算
する際の重み付け係数として設定する手段であることが
好ましい。
【0035】また、この場合、前記記憶手段に記憶され
た前記対応関係における前記複数の重み付け係数を、前
記加算値に対する前記所定方向変動値または前記直交方
向変動値の比率が1/2(nは自然数)となるように
設定し、前記重み付け加算手段は、前記重み付け加算
を、前記重み付け係数に基づくビットシフトにより行う
手段であることが好ましい。
【0036】また、本発明による画像処理装置において
は、前記信号値を推定する画素位置が、前記第1から第
3信号値のいずれか1つの信号値のみからなる単色平坦
部にあるか否かを判定する平坦部判定手段と、前記単色
平坦部にあると判定された場合に、該単色平坦部内にお
ける前記第1から第3信号値に基づく線形補間演算によ
り、前記第1から第3信号値のうち少なくとも1つの信
号値を推定する線形補間手段とをさらに備え、前記単色
平坦部にないと判定された場合に、前記判定手段による
グレーエッジにあるか色エッジにあるかの判定、前記推
定値算出手段による前記所定方向推定値および前記直交
方向推定値の算出、前記変動値算出手段による前記所定
方向変動値および前記直交方向変動値の算出、前記重み
付け係数設定手段による前記重み付け係数の設定、およ
び前記重み付け加算手段による前記重み付け係数に基づ
く信号値の推定を行うことが好ましい。
【0037】さらに、本発明による画像処理装置におい
ては、前記第1から前記第3の光電変換素子が、それぞ
れG、B、Rのいずれかの色に、あるいはY、G、Cの
いずれかの色に分光感度を有するものであることが好ま
しい。
【0038】なお、本発明による画像処理方法をコンピ
ュータに実行させるためのプログラムとして、コンピュ
ータ読取り可能な記録媒体に記録して提供してもよい。
【0039】また、本発明による画像処理装置をデジタ
ルカメラ等の撮像装置、あるいはスキャナ等の画像読取
装置に搭載して提供してもよい。
【0040】
【発明の効果】本発明によれば、信号値を推定する画素
位置がグレーエッジにあるか色エッジにあるかを判定
し、その判定の結果に応じた異なる演算方法により、所
定方向についての信号値の変動量を表す所定方向変動値
および直交方向についての信号値の変動量を表す直交方
向変動値を算出し、この所定方向変動値および直交方向
変動値に基づいて所定方向推定値および直交方向推定値
を重み付け加算する際の重み付け係数を設定するように
したものである。ここで、仮に所定方向にエッジが延在
するとした場合、エッジ上のある画素位置の所定方向推
定値は、エッジ上の他の画素位置の信号値と略一致する
値となるが、直交方向推定値は、エッジに直交する方向
に存在する画素位置における信号値の影響を受けて、エ
ッジ上の他の画素位置の信号値とは異なる値となる。エ
ッジがグレーエッジである場合には、直交方向推定値の
影響がなくなるように重み付け係数を設定することが好
ましい。一方、エッジが色エッジである場合、直交方向
推定値の影響をなくすことが好ましいが、直交方向推定
値の影響を完全になくして信号値を推定すると、エッジ
付近における細かい色味がなくなってしまうため、却っ
て不自然な画像となってしまう。
【0041】ここで、所定方向にエッジが延在するとし
た場合、上記式(3a)、(3b)および式(9a)、
(9b)により算出したSA,SBをそれぞれ所定方向
変動値、直交方向変動値として、これらを比較すると、
式(9a)、(9b)の場合除算を行っていないため、
所定方向変動値に対する直交方向変動値の比率は、式
(9a)、(9b)により算出した場合の方が大きくな
る。したがって、式(9a)、(9b)により算出した
所定方向変動値および直交方向変動値により求められた
重み付け係数は、直交方向推定値に対する重み付けが所
定方向推定値に対する重み付けと比較して、非常に小さ
くなる。一方、式(3a)、(3b)により算出した所
定方向変動値および直交方向変動値により求められた重
み付け係数は、直交方向推定値に対する重み付けが所定
方向推定値に対する重み付けと比較して小さくなるもの
の、式(9a)、(9b)の場合ほど小さくはならな
い。
【0042】したがって、所定方向変動値および直交方
向変動値の演算方法を判定結果に応じて異なるものとし
て重み付け係数を設定することにより、グレーエッジの
場合には偽色をほとんど発生させることなく信号値を推
定でき、色エッジの場合には偽色を低減させつつも細か
な色味を残すことができるように、信号値を推定するこ
とができる。
【0043】なお、上記式(3c)に示すように重み付
け係数を算出してもよいが、この式(3c)において
は、SAを(SA+SB)にて除算しているため、この
演算を実施するための装置の構成が大型化し、また演算
量も多くなる。このため、請求項3,8,13の発明の
ように、所定方向変動値と直交方向変動値との加算値に
対する所定方向変動値または直交方向変動値の比率を重
み付け係数として設定するに際し、重み付け係数が一定
であるとした場合の所定方向変動値と直交方向変動値と
の関係に基づいて、複数の重み付け係数に対する所定方
向変動値および直交方向変動値の対応関係を予め算出し
て、これを記憶手段に記憶しておき、算出された所定方
向変動値および直交方向変動値に基づいて、上記記憶手
段を参照して一の重み付け係数を選択し、この一の重み
付け係数に基づいて、所定方向推定値と直交方向推定値
との重み付け加算を行うことが好ましい。
【0044】ここで、重み付け係数を上記式(3c)に
より算出する場合、これを変動量SBについて解くと、 SB=((1-Sj)/Sj)*SA (3c′) となる。これは重み付け係数Sjが一定であると、変動
値SAと変動値SBとの関係は比例関係にあることを示
している。したがって、変動値SAを所定方向変動値、
変動値SBを直交方向変動値と対応させ、x軸を所定方
向変動値、y軸を直交方向変動値とする2次元平面を設
定すると、この2次元平面上においては、複数の重み付
け係数はその重み付け係数に応じた傾きを有する原点を
通る(厳密には原点を含まないが)直線として表される
ものとなる。
【0045】よって、例えば、重み付け係数に応じた傾
きを有する複数の直線により、2次元平面を複数の領域
に分割して、これを複数の重み付け係数に対する所定方
向変動値および直交方向変動値の対応関係として求めて
これを記憶手段に記憶しておき、算出された所定方向変
動値および直交方向変動値が上記2次元平面上のいずれ
の領域に含まれるかを判断し、これらが位置すると判断
された領域に対応する重み付け係数を一の重み付け係数
として選択することができる。
【0046】したがって、請求項3,8,13の発明に
よれば、重み付け係数を求めるために、式(3c)に示
すように、所定方向変動値または直交方向変動値を所定
方向変動値および直交方向変動値の加算値により除算す
る必要がなくなり、これにより、演算量を低減でき、演
算を行うための回路の構成も簡易なものとすることがで
きる。したがって、本発明を実施するための装置をデジ
タルカメラ等の小型の撮像装置や、スキャナ等の画像読
取装置に搭載することが容易となる。
【0047】また、重み付け係数が加算値に対する所定
方向変動値または直交方向変動値の値が1/2(nは
自然数)となるように設定されている場合、所定方向推
定値および直交方向推定値を2進数で表すと、これらに
1/2を乗算することは、値を右側にnビットシフト
させることと等価な演算を行っていることとなる。した
がって、ビットシフトにより所定方向推定値および直交
方向推定値の重み付け加算を行うことにより、演算を簡
易に行って処理時間を短縮することができる。
【0048】また、信号値を推定する画素位置が、単色
の信号値のみからなる単色平坦部にある場合には、単純
な線形補間演算を行っても偽色を発生させることなく信
号値を推定できる。したがって、信号値を推定する画素
位置が単色平坦部にあるか否かを判定し、単色平坦部に
ある場合には、線形補間演算により信号値を推定するこ
とにより、信号値推定のための演算時間を短縮すること
ができる。
【0049】
【発明の実施の形態】以下図面を参照して本発明の実施
形態について説明する。図1は本発明の第1の実施形態
による画像処理装置の構成を示す概略ブロック図であ
る。図1に示すように本発明の第1の実施形態による画
像処理装置は、単板CCD1を構成する各光電変換素子
において得られた信号値に対して補間処理を施して、全
画素位置における信号値を求めるものであり、各信号値
により構成される画像データS0に対して補間処理を施
して、補間済み画像データS1を得る補間手段2を備え
る。補間手段2は、後述するようにして補間演算により
G信号QGを算出するG補間手段3と、R信号QRおよ
びB信号QBを算出するRB補間手段4と、G信号Q
G、R信号QRおよびB信号QBから後述する図3に示
す空孔画素位置*の信号値を算出する市松正方補間手段
5とを備える。なお、図1に示す単板CCD1において
は、これを構成する光電変換素子からはアナログ信号が
得られるものであるが、本実施形態における画像データ
S0は、アナログ信号をA/D変換したデジタル信号で
あり、なおかつこのデジタル信号は単板CCD1に入力
する光量の0.45乗あるいは対数値となるように変換
されているものとする。
【0050】なお、本実施形態による画像処理装置は、
デジタルカメラのような撮像装置やフイルムから画像を
読み取るスキャナ等の読取装置に設けられるものであっ
てもよく、これらの装置において得られた画像信号を再
生するモニタやプリンタ等の再生装置に設けられるもの
であってもよい。また、画像処理装置単体として用いて
もよい。
【0051】図2は本実施形態において用いられる単板
CCD1の画素配列を示す図である。図2(a)に示す
画素配列は、y方向にR,Gチャンネルに対応する画素
を交互に配列した第1のラインと、y方向にG,Bチャ
ンネルに対応する画素を交互に配列した第2のラインと
を、x方向に交互に配列したものであり、x方向におけ
る各ラインにおいてもR,GチャンネルおよびG,Bチ
ャンネルが交互に配列された正方状の画素配列となって
いる。このような正方状の画素配列としては、図34に
示すベイヤー配列が知られており、本実施形態において
は図2(a)に示す画素配列をベイヤー配列と称するも
のとする。
【0052】また、図2(b)に示す画素配列は、y方
向にR,Bチャンネルに対応する画素を交互に配設した
ラインと、y方向にGチャンネルに対応する画素を配列
したラインとを、各ラインの画素の配列間隔が他のライ
ンの画素配列に対してy方向に略1/2相対的にずれる
ように配列したものであり、市松状の画素配列となって
いる。このような市松状の画素配列としては、上記特開
平10−136391号に記載されたハニカム配列(図
33)が知られており、本実施形態においては図2
(b)に示す画素配列をハニカム配列と称するものとす
る。このハニカム配列は、x方向に対して45度傾斜し
た方向についてみれば、R,Gチャンネルの画素を交互
に配列したラインと、同様に45度傾斜した方向にG,
Bチャンネルの画素を交互に配列したラインとを、この
方向に直交する方向に交互に配列したものとなってい
る。なお、ハニカム配列は図2(a)に示すようなベイ
ヤー配列と画素の配列を45度回転させた関係となって
いる。また、ハニカム配列は上述したように市松状に画
素が配列されてなるものであり、信号値を有さない空孔
画素位置*を用いて図3に示すように正方状に表現する
ことも可能である。本実施形態においては、図2(b)
すなわち図3に示すハニカム配列の単板CCD1におい
て得られた画像データS0に対して処理を行うものとし
て説明する。
【0053】図4は、ハニカム配列の単板CCD1の画
素位置を示す図であり、各画素位置に参照番号を付して
示すものである。ここでは、図4の左上から右下に向か
う方向、すなわち図4において画素がR00,G11,
R22,G33,R44…と並ぶ方向(矢印A方向)を
所定方向とし、この所定方向すなわち矢印A方向に直交
する方向を矢印B方向とする。なお、本実施形態では画
素位置に付与された参照番号を信号値の参照番号として
も使用する。
【0054】図5は、RB補間手段4の構成を示す概略
ブロック図である。図5に示すように、RB補間手段4
は、信号値を推定する画素位置がグレーエッジにあるか
色エッジにあるか否かを判定する判定手段10と、図4
における矢印A方向について全画素位置でのR信号およ
びB信号(総称して信号Q1とする)を算出する第1補
間手段11と、矢印B方向について全画素位置でのR信
号およびB信号(総称して信号Q2とする)を算出する
第2補間手段12と、判定手段10における判定結果を
表す判定情報Hに基づいて、矢印A方向および矢印B方
向における信号値の変動量を表す値を変動値SA,SB
として算出する変動値算出手段13と、変動値算出手段
13において算出された変動値SA,SBに基づいて、
第1補間手段11および第2補間手段12において得ら
れた信号Q1,Q2を重み付け補間する際の重み付け係
数Sjを算出する重み付け係数算出手段14と、重み付
け係数算出手段14において算出された重み付け係数S
jにより、第1補間手段11および第2補間手段におい
て得られた信号Q1,Q2を重み付け加算してR信号Q
RおよびB信号QBを得る重み付け加算手段15とを備
える。なお、図5における重み付け係数算出手段14が
本発明における重み付け係数設定手段と観念される。
【0055】次いで、本実施形態のG補間手段3におけ
るG信号QGおよびRB補間手段4におけるR信号Q
R、B信号QBの算出について説明する。
【0056】(1)まず、RおよびB画素位置における
G信号の算出処理について説明する。このG信号の算出
はG補間手段3において行われる。RおよびB画素位置
におけるG信号は、その画素位置周辺のG画素位置にお
いて得られるG信号に対してスプライン補間演算を施す
ことにより算出することができる。図6はスプライン補
間演算を行う2次元Cubicスプライン補間フィルタの例
を示す図である。図6に示す2次元Cubicスプライン補
間フィルタは、G信号を算出する画素位置の近傍16画
素のG画素位置におけるG信号に対して補間演算を施す
ものである。したがって、図7において実線で囲んだR
44画素位置でのG(R44)信号は、破線で囲んだそ
の周囲の16のG画素位置(G11,G13,G15,
G17,G31,G33,G35,G37,G51,G
53,G55,G57,G71,G73,G75,G7
7)のG信号に対して、図6に示す2次元Cubicスプラ
イン補間フィルタによりフィルタリング処理を施すこと
により算出される。一方、B46画素位置でのG(B4
6)信号は、破線で囲んだその周囲のG画素位置(G1
3,G15,G17,G19,G33,G35,G3
7,G39,G53,G55,G57,G59,G7
3,G75,G77,G79)のG信号に対して、同様
にフィルタリング処理を施すことにより算出される。こ
れにより、補間前にG信号が有する周波数成分を損なう
ことなく補間を行うことができ、この結果、全ての画素
位置において元の周波数情報を保持したG信号を得るこ
とができる。
【0057】ここで、補間演算の手法としては図4にお
ける縦横2次元方向における補間演算であれば、いかな
る手法をも適用することができる。なお、G信号を算出
する画素位置周辺のG画素位置において得られるG信号
に対して単純な線形補間演算を施すことによりR、G画
素位置におけるG信号を算出してもよい。例えば図7に
おいて実線で囲んだR44画素位置でのG(R44)信
号およびB46画素位置でのG(B46)信号は、その
周囲の4つのG画素位置におけるG信号を用いて下記の
式(10)、(11)により算出してもよい。
【0058】 G(R44)=(G33+G35+G53+G55)/4 (10) G(B46)=(G35+G37+G55+G57)/4 (11) (2)次に、RB補間手段4におけるR信号QRおよび
B信号QBの算出について説明する。まず、第1および
第2補間手段11,12における信号Q1,Q2の算出
について説明する。
【0059】1まず、矢印A方向においてR,G画素が
並ぶライン(以下第1のラインとする)におけるG画素
位置での、矢印A方向についてのR信号の算出処理につ
いて説明する。この矢印A方向についての信号の算出は
第1補間手段11において行われる。なお、以降の説明
において、矢印A方向について算出された信号には参照
符号に(A)を付し、矢印B方向について算出された信
号には参照符号に(B)を付するものとする。また、矢
印A方向にR,G画素が並ぶラインを第1のライン、矢
印A方向にB,G画素が並ぶラインを第2のライン、矢
印B方向にR,G画素が並ぶラインを第3のライン、矢
印B方向にB,G画素が並ぶラインを第4のラインとす
る。
【0060】第1のラインにおける矢印A方向について
のG画素位置におけるR信号の算出処理は、そのG画素
位置が含まれる第1のライン上におけるR信号に対して
1次元補間演算を施すことにより行われる。例えば、図
8における実線で囲んだG33画素位置でのR(G3
3)A信号は、G33画素位置が存在する第1のライン
上のG33画素位置周辺の4点のR画素(R00,R2
2,R44,R66)において得られるR信号に対し
て、下記の式(12)に示す1次元Cubicスプライン補
間演算等の1次元補間演算を施すことにより算出され
る。
【0061】 R(G33)A=(-3*ROO+19*R22+19*R44-3*R66)/32 (12) なお、第1のラインにおいて矢印A方向に隣接する2つ
のR画素位置におけるR22,R44信号を用いて下記
の式(13)により算出してもよい。
【0062】 R(G33)A=(R22+R44)/2 (13) また、補間演算の手法としてはこの第1のライン上の矢
印A方向における1次元補間演算であれば、いかなる手
法をも適用することができる。そしてこれにより第1の
ライン上のG画素位置におけるR信号を算出することが
できる。
【0063】2次に、矢印A方向においてR,G画素が
並ぶ第1のラインにおけるG画素位置での、矢印B方向
についてのR信号の算出処理について説明する。この矢
印B方向についての信号の算出は第2補間手段12にお
いて以下のようにして行われる。まず、R信号を算出す
るG画素位置に矢印A方向において隣接する画素位置を
求める。例えば、ここでは図9に示す実線で囲んだG3
3画素位置でのR(G33)B信号を算出するものとす
ると、G33画素位置に隣接する画素位置は破線で囲ん
だR22画素位置およびR44画素位置となる。そして
このR22画素位置におけるR22信号とG(R22)
信号との差およびR44画素位置におけるR44信号と
G(R44)信号との差の平均値を、G33画素位置に
おけるG33信号に加算し、その加算結果をG33画素
位置におけるR(G33)B信号とする。
【0064】ここで、R22信号のみを用いた場合、G
33画素位置でのR(G33)B信号は下記の式(1
4)により算出される。
【0065】 R(G33)B=G33+(R22-G(R22)) (14) なお、式(14)は画像の局所的な領域におけるR信号
とG信号との差は等しいという前提により定められるも
のである。例えば、R22画素位置およびG33画素位
置でのR信号とG信号との差が等しいとは、 R(G33)B-G33=R22-G(R22) (15) の関係を満たすということであり、式(15)をR(G
33)B信号について解いたものが式(14)となる。
なお、式(14)においてG(R22)信号は、上記
(1)の処理により算出されたものではなく、R22画
素位置が存在する第3のライン上のG信号に対して、矢
印B方向に下記の式(16)に示す線形補間演算や式
(12)に示す1次元Cubicスプライン補間演算等の1
次元補間演算を施すことにより算出したものである。し
たがって、G(R22)信号には参照符号Bが付される
こととなる。
【0066】G(R22)B=(G31+G13)/2 (16) すなわち、式(14)では、画像の局所的な領域におけ
るR信号とG信号との差は等しいという前提により信号
値を算出しており、信号値の算出に使用するのは矢印B
方向においてG33画素位置が存在する第4のライン
と、R22画素位置が存在する第3のラインである。こ
れは、第3および第4のライン上にある画素値の相関関
係に基づいて信号値を算出しているということであり、
不明な信号値について矢印B方向における信号値の変化
を反映させて信号値を算出する必要がある。2の処理に
おいて、R22信号のみを用いた場合、G33画素位置
でのR(G33)B信号は上記式(15)に示す関係に
基づいて算出されるが、G(R22)信号が不明である
ため、これを推定する必要がある。ここで、式(15)
は第3および第4のライン間における相関関係を表して
いるため、G(R22)信号を推定するには、第3のラ
イン上において矢印B方向の信号値の変化を反映させる
必要がある。この場合、上記(1)の方法により算出し
たG信号を使用することが考えられるが、このG信号は
式(10)、(11)に示すように、別のラインの信号
値を用いて算出しているものであるため、第3のライン
上における矢印B方向の信号値の変化を反映したもので
はない。このため、矢印B方向の信号値の変化を反映さ
せるために、式(16)のように第3のライン上におい
て、矢印B方向についての1次元補間演算を行ってG
(R22)B信号を算出しているものである。
【0067】一方、R44信号のみを用いると、上記式
(15)の関係に基づいて、G33画素位置でのR(G
44)B信号は下記の式(17)により算出される。
【0068】 R(G33)B=G33+(R44-G(R44)B) (17) なお、G(R44)B信号は式(16)と同様に算出さ
れる。ここで、本実施形態においては、式(14)ある
いは(17)のいずれかによりR(G33)B信号を求
めてもよいが、G画素位置の片側にのみ隣接する画素位
置の信号値を用いたのでは、画像の位相がその方向に近
寄ってしまう。これを防止するために、下記の式(1
8)に示すように、R信号を算出するG画素位置(例え
ばG33画素位置)の両隣の画素位置(R22,R44
画素位置)におけるR信号とG信号との差の平均値を、
R信号を算出するG画素位置におけるG信号に加算した
ものを、そのG画素位置におけるR信号としたものであ
る。
【0069】 R(G33)B=G33+((R22-G(R22)B)+(R44-G(R44)B))/2 (18) そしてこれにより、矢印B方向についてのG画素位置で
のR(G33)B信号を算出することができる。
【0070】3次に、矢印A方向においてB,G画素が
並ぶ第2のラインにおけるG画素位置での、矢印A方向
についてのR信号の算出処理について説明する。この矢
印A方向についての信号の算出は第1補間手段11にお
いて行われる。まず、R信号を算出するG画素位置に矢
印B方向において隣接する画素位置を求める。例えば、
ここでは図10に示す実線で囲んだG35画素位置での
R(G35)A信号を算出するものとすると、G35画
素位置に隣接する画素位置は破線で囲んだR26画素位
置およびR44画素位置となる。そして上記2の処理と
同様に、下記の式(19)に示すように、R26画素位
置におけるR26信号とG(R26)A信号との差およ
びR44画素位置におけるR44信号とG(R44)A
信号との差の平均値を、G35画素位置におけるG35
信号に加算し、その加算結果をG35画素位置における
R(G35)A信号とする。
【0071】 R(G35)A=G35+((R26-G(R26)A)+(R44-G(R44)A))/2 (19) なお、G(R26)A信号およびG(R44)A信号
は、R26およびR44画素位置が存在する第1のライ
ン上におけるG信号に対して、矢印A方向に1次元補間
演算を施すことにより算出したものである。そしてこれ
により、矢印A方向についてのG画素位置でのR(G3
5)A信号を算出することができる。
【0072】4次に、矢印A方向においてB,G画素が
並ぶ第2のラインにおけるG画素位置での、矢印B方向
についてのR信号の算出処理について説明する。この矢
印B方向についての信号の算出は、第2補間手段12に
おいて、そのG画素位置が含まれる第3のライン上にお
けるR信号に対して、1次元補間演算を施すことにより
行われる。例えば、図11における実線で囲んだG35
画素位置でのR(G35)B信号は、G35画素位置が
存在する第3のライン上におけるG35画素位置周辺の
4点のR画素(R08,R26,R44,R62)にお
いて得られるR信号に対して、1の処理と同様に1次元
補間演算を施すことにより算出される。1次元Cubicス
プライン補間演算による算出の例を下記の式(20)に
示す。
【0073】 R(G35)B=(-3*RO8+19*R26+19*R44-3*R62)/32 (20) 5次に、矢印A方向においてB,G画素が並ぶ第2のラ
インにおけるB画素位置での矢印A方向についてのR信
号の算出処理について説明する。この矢印A方向につい
ての信号の算出は第1補間手段11において行われる。
まず、R信号を算出するB画素位置に矢印B方向におい
て隣接する画素位置を求める。例えば、ここでは図12
に示す実線で囲んだB24画素位置でのR(B24)A
信号を算出するものとすると、B24画素位置に隣接す
る画素位置は破線で囲んだG15画素位置およびG33
画素位置となる。そして上記2の処理と同様に、下記の
式(21)に示すように、G15画素位置におけるR
(G15)A信号とG15信号との差およびG33画素
位置におけるR(G33)A信号とG33信号との差の
平均値を、B24画素位置におけるG(B24)A信号
に加算し、その加算結果をB24画素位置におけるR
(B24)A信号とする。
【0074】 R(B24)A=G(B24)A+((R(G33)A-G33)+(R(G15)A-G15))/2 (21) なお、R(G15)A信号およびR(G33)A信号
は、G15およびG33画素位置が存在する第1のライ
ン上におけるR信号に対して、矢印A方向に1次元補間
演算を施すことにより算出したものであり、1の処理に
おいて算出したR信号そのものである。また、G(B2
4)A信号は、B24画素位置が存在する第2のライン
上におけるG信号に対して矢印A方向に1次元補間演算
を施すことにより算出したものである。そしてこれによ
り、矢印A方向についてのB画素位置でのR(B24)
A信号を算出することができる。
【0075】6次に、矢印A方向においてB,G画素が
並ぶ第2のラインにおけるB画素位置での矢印B方向に
ついてのR信号の算出処理について説明する。この矢印
B方向についての信号の算出は第2補間手段12におい
て行われる。まず、R信号を算出するB画素位置に矢印
A方向において隣接する画素位置を求める。例えば、こ
こでは図13に示す実線で囲んだB24画素位置でのR
(B24)B信号を算出するものとすると、B24画素
位置に隣接する画素位置は破線で囲んだG13画素位置
およびG35画素位置となる。そして上記2の処理と同
様に、下記の式(22)に示すように、G13画素位置
におけるR(G13)B信号とG13信号との差および
G35画素位置におけるR(G35)B信号とG35信
号との差の平均値を、B24画素位置におけるG(B2
4)B信号に加算し、その加算結果をB24画素位置に
おけるR(B24)B信号とする。
【0076】 R(B24)B=G(B24)B+((R(G35)B-G35)+(R(G13)B-G13))/2 (22) なお、R(G13)B信号およびR(G35)B信号
は、G13およびG35画素位置が存在する第3のライ
ン上におけるR信号に対して、矢印B方向に1次元補間
演算を施すことにより算出したものであり、4の処理に
おいて算出したR信号そのものである。また、G(B2
4)B信号は、B24画素位置が存在する第4のライン
上におけるG信号に対して、矢印B方向に1次元補間演
算を施すことにより算出したものである。そしてこれに
より、矢印B方向についてのB画素位置でのR(B2
4)B信号を算出することができる。
【0077】以上、1の処理により、矢印A方向におい
てR,G画素が並ぶ第1のラインにおけるG画素位置で
の矢印A方向についてのR信号が算出され、2の処理に
より、矢印A方向においてR,G画素が並ぶ第1のライ
ンにおけるG画素位置での矢印B方向についてのR信号
が算出される。
【0078】また、3の処理により、矢印A方向におい
てB,G画素が並ぶ第2のラインにおけるG画素位置で
の矢印A方向についてのR信号が算出され、4の処理に
より、矢印A方向においてB,G画素が並ぶ第2のライ
ンにおけるG画素位置での矢印B方向についてのR信号
が算出される。
【0079】さらに、5の処理により、矢印A方向にお
いてB,G画素が並ぶ第2のラインにおけるB画素位置
での矢印A方向についてのR信号が算出され、6の処理
により、矢印A方向においてB,G画素が並ぶ第2のラ
インにおけるB画素位置での矢印B方向についてのR信
号が算出される。
【0080】以上、各画素位置における矢印A,B両方
向でのR信号の算出について説明したが、B信号につい
ても上記1から6の処理において説明したようにR信号
と同様にして算出することができる。
【0081】判定手段10は、以下のようにして信号値
を推定する画素位置がグレーエッジにあるか色エッジに
あるかを判定する。図14は赤(R)および緑(G)か
らなる色エッジの例を示す図、図15は白および黒から
なるグレーエッジの例を示す図である。なお、図14,
15においては、信号値が0の画素位置を白抜きで示
し、それ以外の画素は信号値が255であるものとして
説明する。また、ここではR44画素位置での判定につ
いて説明する。まず、判定手段10は、下記の式(2
3)から(26)により、エッジの判断値Dat1〜D
at4を求める。
【0082】 Dat1=|R44-G33|+|R44-G55|+|G35-B24|+|G35-B46| (23) Dat2=|R44-G33|+|R44-G55|+|G53-B42|+|G53-B64| (24) Dat3=|R44-G35|+|R44-G53|+|G33-B24|+|G33-B42| (25) Dat4=|R44-G35|+|R44-G53|+|G55-B46|+|G55-B64| (26) ここで、式(23)は矢印A方向において、R44画素
位置の右上部分のエッジの色を判断するための判断値D
at1を算出するものであり、式(24)は矢印A方向
において、R44画素位置の左下部分のエッジの色を判
断するための判断値Dat2を算出するものである。ま
た、式(25)は矢印B方向において、R44画素位置
の左上部分のエッジの色を判断するための判断値Dat
3を算出するものであり、式(26)は矢印B方向にお
いて、R44画素位置の右下部分のエッジの色を判断す
るための判断値Dat4を算出するものである。
【0083】信号値を推定する画素位置がグレーエッジ
にある場合は、上記判断値Dat1〜Dat4のうちい
ずれか1つは必ず値が小さくなる。例えば、図15に示
すグレーエッジの場合、上記式(23)〜(26)にお
いて算出される判断値Dat1〜Dat4において、判
断値Dat1およびDat2が0となる。一方、信号値
を推定する画素位置が図14に示すような色エッジにあ
る場合は、上記判断値Dat1〜Dat4のいずれの値
も大きくなる。したがって、判断値Dat1〜Dat4
の最小値Dminを予め定められた閾値Tと比較し、最
小値Dminが閾値T以下である場合、すなわち、 Dmin=min(Dat1,Dat2,Dat3,Dat4)≦T (27) である場合には、その画素位置をグレーエッジであると
判定し、それ以外の場合、すなわち、 Dmin=min(Dat1,Dat2,Dat3,Dat4)>T (28) である場合には色エッジであると判定する。ここで、閾
値Tの値としては、信号値が8ビットで表される場合1
28程度の値とする。そして、判定手段10はこの判定
結果を判定情報Hとして変動値算出手段13に入力す
る。
【0084】変動値算出手段13は、後述する重み付け
係数算出手段14において、1と2、3と4および5と
6の処理により算出された矢印A,B方向の信号値を重
み付け加算するための重み付け係数の算出に必要な変動
値SA,SBを算出するものである。ここで、変動値S
A,SBとは信号値を推定する画素位置において、信号
値の変動量を表す値のことをいう。まず、1および2の
処理により算出された矢印A,B方向のR信号の重み付
け係数を算出するための変動値SA,SBの算出につい
て説明する。
【0085】まず、判定情報Hにより信号値を算出する
画素位置が色エッジにあると判定された場合、矢印A方
向についての変動値SAおよび矢印B方向についての変
動値SBを下記の式(29)、(30)により算出す
る。なお、ここでは1および2の処理と同様に、図8に
示すG33画素位置での信号値の変動値SA,SBを算
出するものとする。
【0086】 SA=|R(G33)A-G33|/(R(G33)A+G33) (29) SB=|B(G33)B-G33|/(B(G33)B+G33) (30) ここで、式(29)におけるR(G33)Aは1の処理
において式(12)により算出された信号値、式(3
0)におけるB(G33)Bは第4のライン上のG33
画素位置周辺の4点のB画素位置(B06,B24,B
42,B60)において得られるB信号に対して、1次
元補間演算を施すことにより算出された信号値である。
このように算出された変動値SAは、G33画素位置の
矢印A方向における信号値の変化を表し、変動値SBは
G33画素位置の矢印B方向における信号値の変化を表
す。
【0087】一方、判定情報Hにより信号値を算出する
画素位置がグレーエッジにあると判定された場合、変動
値SAおよび変動値SBを式(29)、(30)の分母
を省略した下記の式(31)、(32)により算出す
る。
【0088】 SA=|R(G33)A-G33| (31) SB=|B(G33)B-G33| (32) 次に、3および4の処理により算出された矢印A,B方
向のR信号の重み付け係数を算出するための変動値の算
出について説明する。
【0089】まず、判定情報Hにより信号値を算出する
画素位置が色エッジにあると判定された場合、矢印A方
向についての変動値SAおよび矢印B方向についての変
動値SBを下記の式(33)、(34)により算出す
る。なお、ここでは3および4の処理と同様に、図10
に示すG35画素位置での信号値の変動値を算出するも
のとする。
【0090】 SA=|B(G35)A-G35|/(B(G35)A+G35) (33) SB=|R(G35)B-G35|/(R(G35)B+G35) (34) ここで、式(33)におけるB(G35)Aは第2のラ
イン上のG35画素位置周辺の4点のB画素位置(B0
2,B24,B46,B68)において得られるB信号
に対して、1次元補間演算を施すことにより算出された
信号値、式(34)におけるR(G35)Bは上記4の
処理において式(20)により算出された信号値であ
る。このように算出された変動値SAは、G35画素位
置の矢印A方向における信号値の変化を表し、変動値S
BはG35画素位置の矢印B方向における信号値の変化
を表す。
【0091】一方、判定情報Hにより信号値を算出する
画素位置がグレーエッジにあると判定された場合、変動
値SAおよび変動値SBを式(33)、(34)の分母
を省略した下記の式(35)、(36)により算出す
る。
【0092】 SA=|B(G35)A-G35| (35) SB=|R(G35)B-G35| (36) 次に、5および6の処理により算出された矢印A,B方
向のR信号の重み付け係数を算出するための変動値の算
出について説明する。
【0093】まず、判定情報Hにより信号値を算出する
画素位置が色エッジにあると判定された場合、矢印A方
向についての変動値SAおよび矢印B方向についての変
動値SBを下記の式(37)、(38)により算出す
る。なお、ここでは5および6の処理と同様に、図1
2、図13に示すB24画素位置での信号値の変動値を
算出するものとする。
【0094】 SA=|B24-G(B24)A|/(B24+G(B24)A) (37) SB=|B24-G(B24)B|/(B24+G(B24)B) (38) ここで、式(37)におけるG(B24)Aは、上記5
の処理においてB24画素位置が存在する第2のライン
上のG信号に対して、矢印A方向に1次元補間演算を施
すことにより算出した信号値、式(38)におけるG
(B24)Bは、上記6の処理においてB24画素位置
が存在する第4のライン上のG信号に対して、矢印B方
向に1次元補間演算を施すことにより算出した信号値で
ある。このように算出された変動値SAは、B24画素
位置の矢印A方向における信号値の変化を表し、変動値
SBはB24画素位置の矢印B方向における信号値の変
化を表す。
【0095】一方、判定情報Hにより信号値を算出する
画素位置がグレーエッジにあると判定された場合、変動
値SAおよび変動値SBを式(37)、(38)の分母
を省略した下記の式(39)、(40)により算出す
る。
【0096】 SA=|B24-G(B24)A| (39) SB=|B24-G(B24)B| (40) 以上のように算出された変動値SA,SBは重み付け係
数算出手段14に入力され、下記の式(41)により、
重み付け係数Sjが算出される。重み付け係数Sjは、
重み付け加算手段15に入力される。
【0097】 Sj=SA/(SA+SB) (if SA+SB=0 then Sj=0.5) (41) なお、式(41)においてSA+SB=0となるのは、
SA,SB≧0であることから、SA=SB=0の場合
である。
【0098】重み付け加算手段15においては、1と
2、3と4および5と6のそれぞれの場合について、下
記の式(42)〜(44)により信号値が算出される。
なお、式(42)〜(44)は、それぞれ矢印A方向に
おいてR,G画素が並ぶ第1のラインにおけるG33画
素位置でのR(G33)信号、矢印A方向においてB,
G画素が並ぶ第2のラインにおけるG35画素位置での
R(G35)信号、矢印A方向においてB,G画素が並
ぶ第2のラインにおけるB24画素位置でのR(B2
4)信号を算出するものである。
【0099】 R(G33)=(1-Sj)*R(G33)A+Sj*R(G33)B (42) R(G35)=(1-Sj)*R(G35)A+Sj*R(G35)B (43) R(B24)=(1-Sj)*R(B24)A+Sj*R(B24)B (44) なお、信号Q1、Q2を用いた一般式を式(45)に示
す。
【0100】 QR,QB=(1-Sj)*Q1+Sj*Q2 (45) これにより、信号値の変化が少ない方向ほど大きな重み
付けがなされて、R、B信号QR,QBが算出されるこ
ととなる。
【0101】ここで、色エッジおよびグレーエッジの場
合における変動値SA,SBの算出方法の差異について
説明する。図16は変動値SA,SBの算出方法の差異
を説明するための図である。図16に示す画素配列にお
いては白(信号値が255)および黒(信号値が0また
は0に近い値)のグレーエッジが表されているものとす
る。なお、ここではR44画素におけるB(R44)信
号の算出について説明する。まず、矢印A方向につい
て、下記の式(46)によりR44画素位置のB(R4
4)A信号を求める。
【0102】 B(R44)A=G(R44)A+((B(G35)A-G35)+(B(G53)A-G53))/2 =0+((2-0)+(255-255))/2=1 (46) 但し、G(R44)A=0,B(G35)A=2,B(G53)A=255であり、こ
れらの値の算出には、式(12)、(20)に示すよう
な1次元Cubicスプライン補間演算を用いた。
【0103】また、矢印B方向について、下記の式(4
7)によりR44画素位置のB(R44)B信号を求め
る。
【0104】 B(R44B)=G(R44)B+((B(G33)B-G33)+(B(G55)B-G55))/2 =103+((105-0)+(105-0))/2=208 (47) 但し、G(R44)B≒103,B(G33)B≒105,B(G55)B≒105であ
り、これらの値の算出には、式(12)、(20)に示
すような1次元Cubicスプライン補間演算を用いた。
【0105】そして、下記の式(48)により変動値S
A,SBを求め、式(49)によりB(R44)A信号
とB(R44)B信号とを重み付け加算するための重み
付け係数Sjを算出する。
【0106】 SA=|R44-G(R44)A|/(R44+G(R44)A)=|1-0|/(1+0)=1 SB=|R44-G(R44)B|/(R44+G(R44)B)=|1-103|/(1+103)≒0.98 (48) Sj=SA/(SA+SB)=1/(1+0.98)≒0.5 (49) このように重み付け係数Sjが求められると、下記の式
(50)によりB(R44)信号を算出する。
【0107】 B(R44)=(1-Sj)*B(R44)A+Sj*B(R44)B=0.5*1+0.5*208≒105 (50) ここで、R44画素位置は黒いエッジ上にあるため、B
(R44)信号の値は0に近い値であることが好まし
い。しかしながら、上記式(48)、(49)に示すよ
うに重み付け係数Sjを算出して式(50)によりB
(R44)信号の値を求めるとその値は105となるた
め、ここに偽色が発生することとなる。これは、矢印A
方向において信号値の変動が少ないほど重み付け係数S
jの値は0に近くなるが、信号値を算出する画素位置周
辺の画素位置における信号値によっては、式(50)に
示すように、重み付け係数Sjが0に近い値とならない
場合があるからである。
【0108】一方、本実施形態においては、信号値を算
出する画素位置がグレーエッジにある場合には、式(4
8)において分母を省略して変動値SA,SBを算出す
るものであり、R44画素位置においてそのように変動
値SA,SBを算出すると、 SA=|R44-G(R44)A|=|1-0|=1 SB=|R44-G(R44)B|=|1-103|=102 (51) となり、重み付け係数Sjは、 Sj=SA/(SA+SB)=1/(1+102)≒0.01 (52) となる。このように算出された重み付け係数Sjにより
B(R44)信号を算出すると、 B(R44)=(1-Sj)*B(R44)A+Sj*B(R44)B=0.99*1+0.01*208≒2.99 (53) となる。これにより、B(R44)信号は0に近い値と
なり、その結果偽色を低減することができる。
【0109】一方、色エッジであると判定された場合
に、式(51)により変動値SA,SBを算出しても偽
色を低減できるが、エッジ付近における細かな色味がな
くなってしまうため、画像が不自然なものとなる。した
がって、本実施形態においては、信号値を算出する画素
位置がグレーエッジであるか色エッジであるかに応じ
て、変動値SA,SBの算出方法を変更するようにした
ものである。
【0110】なお、図2(a)に示すベイヤー配列の単
板CCD1においては画素配列を45度回転させると、
その画素配列は図2(b)に示すハニカム配列と空孔画
素位置*を除けば同様のものとなる。したがって、上記
の説明において、矢印A方向を図2における紙面上下方
向、矢印B方向を図2における紙面左右方向とすること
により、ベイヤー配列の場合もハニカム配列の場合と同
様に、全画素位置におけるRGB信号を算出することが
できる。
【0111】このようにして、全ての画素位置において
RGB全ての信号値QR,QG,QBが得られると、市
松正方補間手段5において、図3における空孔画素位置
*でのRGB信号値を補間演算により算出して、画素が
正方状に配列されてなる補間済み画像データS1を得
る。この補間演算は、図17に示すような空孔画素位置
*周辺の4画素位置の信号値を用いた補間フィルタや、
図18に示すように4×4画素についての2次元Cubic
スプライン補間演算を行う補間フィルタの補間係数の配
置を45度傾斜させた補間係数の配置を有する補間フィ
ルタを用いて補間演算を行うことにより求めることがで
きる。なお、この空孔画素位置*における信号値を算出
するための補間演算を市松正方補間演算とする。また、
補間演算についてはこれらに限定されるものではなく、
上述したように求められた各画素位置におけるRGB信
号をYCC輝度色差空間に変換し、YCC毎に異なる補
間フィルタによる補間演算を施す等、空孔画素位置*に
おける信号値を算出するための補間演算であればいかな
る方法をも採用可能である。
【0112】次いで、第1の実施形態の動作について説
明する。図19は第1の実施形態の動作を示すフローチ
ャートである。まず、被写体を撮影して単板CCD1に
おいて画像データS0を得る(ステップS1)。次い
で、補間手段2のG補間手段3において、上記(1)の
処理により所定方向におけるRまたはB画素位置におけ
るG信号QGが算出される(ステップS2)。そして、
補間手段2において、算出する信号をR信号にセットし
(ステップS3)、RB補間手段4においてR信号QR
が算出される(ステップS4)。
【0113】図20はステップS4の処理を示すフロー
チャートである。まず、上記処理1および2により信号
値を算出し(ステップS11)、ステップS3および4
により信号値を算出し(ステップS12)、処理5およ
び6により信号値を算出する(ステップS13)。一
方、判定手段10においては、信号値を算出する画素位
置がグレーエッジにあるか否かが判定され(ステップS
14)、その判定結果を表す判定情報Hが変動値算出手
段13に入力される。変動値算出手段13においては、
信号値を算出する画素位置がグレーエッジにある場合に
は、上記式(31)、(32)等に示すようにグレーエ
ッジ用の変動値SA,SBを算出する(ステップS1
5)。逆に色エッジである場合には、上記式(29)、
(30)等に示すように色エッジ用の変動値SA,SB
を算出する(ステップS16)。そして、変動値SA,
SBは重み付け係数算出手段14に入力されて、ここで
重み付け係数Sjが算出され(ステップS17)、ステ
ップS11,S12,S13の処理により算出された信
号値が重み付け加算されて、R信号QRが算出される
(ステップS18)。
【0114】図19に戻り、R,B全ての信号を算出し
たか否かが判断され(ステップS5)、ステップS5が
否定された場合には、ステップS6において算出する信
号をB信号にセットし、ステップS4に戻ってステップ
S4,S5の処理を繰り返してB信号QBが算出され
る。
【0115】ステップS5が肯定されると、空孔画素位
置以外の全ての画素位置においてRGB信号QR,Q
G,QBが算出されたとして、市松正方補間手段5にお
いて、各画素の信号値に対して図17あるいは図18に
示す補間フィルタによる市松正方補間演算を行って、空
孔画素位置における信号値を算出し(ステップS7)、
処理を終了する。
【0116】そしてこれにより、図2(b)に示すハニ
カム配列による単板CCD1の空孔画素位置を含む全画
素位置におけるRGB信号を得ることができ、全画素位
置においてRGB信号を有する補間済み画像データS1
を得ることができる。
【0117】なお、図19に示すフローチャートのステ
ップS3においては、算出する信号をR信号にセットし
ているが、先にB信号にセットしてもよい。この場合、
ステップS5が否定されるとステップS6において、算
出する信号がB信号からR信号に切り替えられることと
なる。
【0118】このように、第1の実施形態によれば、信
号値を算出する画素位置がグレーエッジにあるか色エッ
ジにあるかを判定し、判定結果に応じて重み付け係数S
jを算出するための変動値SA,SBの算出方法を変更
するようにしたため、信号値を算出する画素位置がグレ
ーエッジにある場合であっても、偽色を低減できる。ま
た、信号値を算出する画素位置が色エッジにある場合に
は、エッジ付近における細かい色味を残して自然な感じ
となる画像を得ることができる。
【0119】次いで、本発明の第2の実施形態について
説明する。
【0120】図21は本発明の第2の実施形態による画
像処理装置のRB補間手段の構成を示す概略ブロック図
である。なお、図21において図5と同一の構成につい
ては同一の参照番号を付し、詳細な説明は省略する。第
2の実施形態においては、判定手段10に代えて、信号
値を推定する画素位置がグレーエッジにあるか色エッジ
にあるかの判定に加えて、単色の信号値からなる単色平
坦部にあるか否かを判定する判定手段21を備え、さら
に、画素位置が単色平坦部にあると判定された場合に、
線形補間演算によりその画素位置におけるR信号QRお
よびB信号QBを算出する線形補間手段22を備えた点
が第1の実施形態と異なるものである。
【0121】判定手段21におけるグレーエッジにある
か色エッジにあるかの判定処理は、第1の実施形態にお
ける判定手段10における判定処理と同一であり、グレ
ーエッジまたは色エッジにあると判定された場合は、そ
の判定結果を表す判定情報Hが変動値算出手段13に入
力され、第1の実施形態と同様に変動値SA,SBが算
出され、さらに重み付け係数算出手段14において重み
付け係数Sjが算出される。この場合、第1の実施形態
と同様に第1補間手段11および第2補間手段12にお
いて、各画素位置における矢印A方向および矢印B方向
の信号Q1,Q2が算出され、その信号値が重み付け加
算手段15に入力されて、重み付け係数Sjにより重み
付け加算されてR信号QRおよびB信号QBが算出され
る。
【0122】以下、判定手段21における単色平坦部に
あるか否かの判定について説明する。図22は緑(G)
の単色平坦部の信号値を示す図、図23は青(B)の単
色平坦部の信号値を示す図、図24は赤(R)の単色平
坦部の信号値を示す図である。図22においてはG画素
位置のみが255の信号値を有し、図23においてはB
画素位置のみが255の信号値を有し、図24において
はR画素位置のみが255の信号値を有するものとなっ
ている。判定手段21は、信号値を算出する画素位置お
よびその画素位置周辺のR信号の平均値をDat11、
G信号の平均値をDat12、B信号の平均値をDat
13として求める。例えば、図22から図24に示す画
素配列において、R44画素位置についてみると、 Dat11=(R22+R26+R44+R62+R66)/5 Dat12=(G33+G35+G53+G55)/4 (54) Dat13=(B24+B42+B46+B64)/4 により、平均値Dat11〜Dat13を算出する。
【0123】また、B42画素位置については、 Dat11=(R22+R40+R44+R62)/4 Dat12=(G31+G33+G51+G53)/4 (55) Dat13=(B20+B24+B42+B60+B64)/5 により、平均値Dat11〜Dat13を算出する。
【0124】さらに、G33画素位置については、 Dat11=(R04+R22+R44+R62)/4 Dat12=(G11+G13+G15+G31+G33+G35+G51+G53+G55)/9 (56) Dat13=(B02+B24+B42+B64)/4 により、平均値Dat11〜Dat13を算出する。
【0125】そして、各平均値Dat11〜Dat13
の差分値の絶対値|Dat11-Dat12|、|Dat12-Dat13|、|Dat
13-Dat11|を算出し、これらの和が予め定められた閾値
T1より大きいか否かを判断し、大きい場合、すなわち
下記の式(57)の条件を満たす場合にはその画素位置
は単色平坦部にあると見なすものである。
【0126】 |Dat11-Dat12|+|Dat12-Dat13|+|Dat13-Dat11|>T1 (57) すなわち、図22に示すGの単色平坦部の場合、G画素
位置のみが信号値を有するため、式(54)〜(56)
においては、平均値Dat12のみが値を有し、Dat
11,Dat13は値を有さない。したがって、|Dat11
-Dat12|、|Dat12-Dat13|の値は非常に大きい値となる。
また、図23に示すBの単色平坦部の場合、B画素位置
のみが信号値を有するため、式(54)〜(56)にお
いては、平均値Dat13のみが値を有し、Dat1
1,Dat12は値を有さない。したがって、|Dat12-D
at13|、|Dat13-Dat11|の値は非常に大きい値となる。さ
らに、図24に示すRの単色平坦部の場合、R画素位置
のみが信号値を有するため、式(54)〜(56)にお
いては、平均値Dat11のみが値を有し、Dat1
2,Dat13は値を有さない。したがって、|Dat11-D
at12|、|Dat13-Dat11|の値は非常に大きい値となる。
【0127】このように、信号値を算出する画素位置が
単色平坦部にあると、|Dat11-Dat12|、|Dat12-Dat13|、
|Dat13-Dat11|のうち2つは大きな値となるため、これ
らの値の和も大きな値となる。したがって、閾値T1を
適切に設定しておき、上記式(57)に示す関係を満た
すか否かを判断することにより、画素位置が単色平坦部
にあるか否かを判定することができる。
【0128】ここで、信号値を算出する画素位置が単色
平坦部にある場合には、上記第1補間演算11および第
2補間演算12における信号値の算出を行わなくとも、
単純に線形補間演算を行っても偽色は発生しないもので
ある。したがって、判定手段21は、信号値を算出する
画素位置が単色平坦部にあると判定された場合には、画
像データS0を線形補間手段22に入力する。そして、
線形補間手段22においては、信号値を算出する画素位
置とその周囲の画素位置との信号値を用いて線形補間演
算を行って、信号値を算出する。
【0129】例えば、G画素位置におけるR信号を算出
する場合、信号値を算出する画素位置をG33画素位置
とすると、 R(G33)=(R22+R24)/2 (58) によりR(G33)信号を算出する。また、G33画素
位置におけるB(G33)信号を算出する場合、 B(G33)=(B24+B42)/2 (59) によりB(G33)信号を算出する。
【0130】また、R画素位置におけるB信号を算出す
る場合、信号値を算出する画素位置をR44とすると、 B(R44)=(B24+B42+B46+B64)/4 (60) によりB(R44)信号を算出する。
【0131】さらに、B画素位置におけるR信号を算出
する場合、信号値を算出する画素位置をB24とする
と、 R(B24)=(R04+R22+R26+R44)/4 (61) によりR(B24)信号を算出する。
【0132】次いで、第2の実施形態の動作について説
明する。なお、第2の実施形態の動作は基本的には図1
9に示す第1の実施形態の動作と同一であり、図19に
おけるステップS4における処理のみが異なるものであ
るため、ここではステップS4における信号算出処理に
ついてのみ説明する。図25は第2の実施形態における
信号算出処理の動作を示すフローチャートである。ま
ず、判定手段21において、上記式(57)により、信
号値を算出する画素位置が単色平坦部にあるか否かが判
断される(ステップS21)。単色平坦部にある場合に
はステップS21が肯定され、線形補間手段22におい
て上記式(58)および(61)に示す線形補間演算に
よりR信号QRが算出され(ステップS30)、処理を
終了する。なお、B信号QBについては、上記式(5
9)および(60)に示す線形補間演算により算出され
る。
【0133】ステップS21が否定された場合は、ステ
ップS22に進み、第1の実施形態と同様に、処理1お
よび2による信号値の算出(ステップS22)、処理3
および4による信号値の算出(ステップS23)、処理
5および6による信号値の算出(ステップS24)、グ
レーエッジか否かの判断(ステップS25)、グレーエ
ッジ用変動値算出(ステップS26)、色エッジ用変動
値算出(ステップS27)、重み付け係数算出(ステッ
プS28)および重み付け加算(ステップS29)が行
われて、R信号QR、B信号QBが算出される。
【0134】このように、第2の実施形態においては、
R、B信号QR、QBを算出する際に、まず信号値を算
出する画素位置が単色平坦部にあるか否かを判断し、単
色平坦部にある場合には、線形補間演算により信号値を
算出するようにしたため、全ての画素について第1の実
施形態の方法により信号値を算出する場合と比較して、
単色平坦部の画素位置については演算量を低減すること
ができ、これにより高速に処理を行うことができること
となる。
【0135】ここで、単板CCD1を使用した撮像装置
においては、偽色を低減するために、画像データを得る
ための撮像系において、光学ローパスフィルタを用いる
必要がある。図26はハニカム配列の単板CCDを用い
たデジタルカメラの撮像系を示す図、図27はハニカム
配列の単板CCDにおいて得られた信号値の周波数特性
を示す図である。なお、図27ではfs/2がナイキス
ト周波数であり、図27(a)における実線は従来の処
理により再現可能な周波数帯域、図27(b)における
実線は本実施形態の処理により再現可能な周波数帯域を
それぞれ示している。図26に示すように、デジタルカ
メラの撮像系は被写体から順に撮影レンズ31、光学ロ
ーパスフィルタ32およびハニカム配列の単板CCD3
3が配列されてなる。なお、光軸をXで示す。上記特開
平10−20906号、同10−136391号等に記
載された方法においては、偽色の低減効果が不十分であ
るため、図27(a)に示す実線より外側の周波数帯域
をカットするための帯域の狭い光学ローパスフィルタ3
2を使用する必要がある。しかしながら、このような帯
域の狭い光学ローパスフィルタを使用すると、得られる
画像の解像度が低下する。また、このような周波数帯域
の狭い光学ローパスフィルタ32を作成するには、中心
波長630〜785nm程度の1/4波長板を図28
(a)に示すように3枚張り合わせる必要があるため、
撮像系の構成が大型化する。
【0136】本発明によれば、従来の方法と比較して偽
色を低減することができるため、光学ローパスフィルタ
32の周波数帯域を図27(b)に示すように広げるこ
とができ、これにより、単板CCD33において得られ
る画像データS0により表される画像の解像度を向上さ
せることができる。また、図28(b)に示すように2
枚の1/4波長板を張り合わせることにより、図27
(b)に示す特性を得ることができるため、撮像系を小
型化することができる。
【0137】なお、上記第1および第2の実施形態にお
いては、重み付け係数算出手段14において、式(4
1)に示すように、変動値SAに対して(SA+SB)
による除算を行うことにより、重み付け係数Sjを算出
している。
【0138】 Sj=SA/(SA+SB) (if SA+SB=0 then Sj=0.5) (41) ここで、種々の演算を行う演算回路においては、演算が
加算、減算および乗算のみである場合には、比較的少な
い演算素子により回路を構成することができる。しかし
ながら、上記式(41)のように除算が含まれると、演
算が複雑となるため多数の演算素子が必要となり回路の
構成が複雑化するとともに、演算に長時間を要するもの
となる。
【0139】ここで、重み付け係数Sjを式(41)に
より算出する場合、これを変動量SBについて解くと、 SB=((1-Sj)/Sj)*SA (41′) となる。これは重み付け係数Sjが一定であると、変動
値SAと変動値SBとの関係は、比例関係にあることを
示している。したがって、x軸を変動値SA、y軸を変
動値SBとする2次元平面を設定すると、この2次元平
面上においては、図29に示すように、重み付け係数S
jはその重み付け係数に応じた傾きを有する原点を通る
(厳密には原点を含まないが)直線として表されるもの
となる。以下、この関係を用いて式(41)における除
算処理を省略する実施形態を第3の実施形態として説明
する。
【0140】図30は第3の実施形態の処理を説明する
ための図、図31は第3の実施形態による画像処理装置
のRB補間手段の構成を示す概略ブロック図である。な
お、図31において図5と同一の構成については同一の
参照番号を付し、詳細な説明は省略する。第3の実施形
態においては、重み付け係数算出手段14に代えて、第
1補間手段11および第2補間手段12において得られ
た信号Q1,Q2を重み付け補間する際の重み付け係数
Sjに対する変動値SA,SBの対応関係を、複数の重
み付け係数について記憶したメモリ14Aと、変動値算
出手段13において算出された変動値SA,SBに基づ
いて、メモリ14Aに記憶された対応関係を参照して複
数の重み付け係数から一の重み付け係数Sjを選択する
重み付け係数選択手段24とを備え、重み付け加算手段
15においては、重み付け係数選択手段24において選
択された重み付け係数Sjにより、重み付け加算を行う
ようにしたものである。なお、第3の実施形態において
は、重み付け係数選択手段24およびメモリ14Aが本
発明における重み付け係数設定手段と観念される。
【0141】まず、重み付け係数Sjに対する変動値S
A,SBの対応関係の求め方について説明する。第3の
実施形態においては、重み付け係数Sjが一定であると
変動値SAと変動値SBとは比例関係にあり、SA−S
B平面においては重み付け係数Sjはその値に応じた傾
きを有する直線となる、ということを利用して、まず重
み付け係数Sjを図29に示すようにSj=0,0.2
5,0.5,0.75,1.0の5段階に量子化する。
そして、図30に示すようにSA−SB平面上において
Sj=0,0.125,0.375,0.625,0.
875,1を境界として5つの領域A1〜A5を設定
し、各領域A1〜A5をそれぞれSj=0,0.25,
0.5,0.75,1.0に対応付ける。ここで、Sj
=0.125の場合SB=7SA、Sj=0.375の
場合SB=5/3SA、Sj=0.625の場合SB=
3/5SA、Sj=0.875の場合SB=1/7SA
となる。
【0142】したがって、図31に示すように、メモリ
14Aを設け、このメモリ14Aに図30に示す関係を
テーブルとして記憶させておき、重み付け係数選択手段
24において変動値SA,SBに基づいてこのテーブル
を参照することにより、式(41)に示すように除算を
行うことなく重み付け係数Sjを求めることができる。
例えば、算出された変動値SA,SBの値が3/5SA
<SB≦5/3SAの関係、すなわち3SA<5SBか
つ5SA≧3SBの関係にある場合には、変動値SA,
SBは図30に示すグラフの領域A3にプロットされる
ことから、重み付け係数をSj=0.5として求めるこ
とができる。また、7SA<SBの関係にある場合には
プロットは領域A1にあるためSj=0、5/3SA<
SB≦7SAの関係、すなわち5SA<3SBかつ7S
A≧SBの関係にある場合にはプロットは領域A2にあ
るためSj=0.25、1/7SA<SB≦3/5SA
の関係、すなわちSA<7SBかつ3SA≧5SBの関
係にある場合にはプロットは領域A4にあるためSj=
0.75、SB≦1/7SAの関係、すなわちSA≧7
SBの関係にある場合にはプロットは領域A5にあるた
めSj=1として求めることができる。
【0143】なお、重み付け係数Sjはより細かく量子
化した方が好ましいが、変動値SA,SBが存在する領
域の判断のための演算時間を考慮すると、5段階程度で
あることが好ましい。本出願人の実験によれば、5段階
に量子化したものであっても、十分に偽色低減効果が得
られるものである。
【0144】そして、重み付け加算手段15において
は、この重み付け係数Sjを用いて上記式(42)〜
(44)に示す重み付け加算が行われて信号値が算出さ
れる。例えばR44画素におけるB(R44)信号の算
出は、重み付け係数Sj=0.5とした場合、 B(R44)=(1-Sj)*B(R44)A+Sj*B(R44)B=0.5*B(R44)A+0.5*B(R44)B (62) として算出されるが、第3の実施形態においては、Sj
または(1−Sj)=1/2(nは自然数)となるよ
うに重み付け係数が設定されているため、下記の式(6
3)のように、ビットシフトを用いて整数化して演算を
行うことも可能である。
【0145】 B(R44)=(B(R44)A+B(R44)B+1)>>1 (63) ここで、>>1は信号値を2進数で表した際に、右側に
1ビットシフトさせることを表す。このように、右側に
1ビットシフトさせることにより、信号値は1/2倍さ
れるため、式(62)と同様の演算を行っていることと
なる。なお、式(63)において1を加算しているの
は、演算結果を四捨五入するためである。
【0146】さらに、上記第1および第2の実施形態に
おいて、変動値算出手段13において、色エッジである
場合の変動値SA,SBの算出は、上記式(29)、
(30)等により行っているが、これらの式においては
除算が含まれるため、演算量を低減させるためには、除
算を行わないことが好ましい。この場合、変動値SA,
SBの比率が変わらなければ重み付け係数Sjの算出に
支障はない。したがって、上記式(29)、(30)の
場合、右辺に(R(G33)A+G33)*(B(G33)B+G33)を乗算し、
下記の式(29′)、(30′)により、変動値SA,
SBを算出することが好ましい。
【0147】 SA=|R(G33)A-G33|*(B(G33)B+G33) (29′) SB=|B(G33)B-G33|*(R(G33)A+G33) (30′) なお、上記式(33)、(34)については、(B(G35)A
+G35)*(R(G35)B+G35)を乗算して、下記の式(3
3′)、(34′)により変動値SA,SBを算出すれ
ばよい。
【0148】 SA=|B(G35)A-G35|*(R(G35)B+G35) (33′) SB=|R(G35)B-G35|*(B(G35)A+G35) (34′) また、上記式(37)、(38)については、(B24+G(B
24)A)*(B24+G(B24)B)を乗算して、下記の式(3
7′)、(38′)により変動値SA,SBを算出すれ
ばよい。
【0149】 SA=|B24-G(B24)A|*(B24+G(B24)B) (37′) SB=|B24-G(B24)B|*(B24+G(B24)A) (38′) 以下、第3の実施形態においてビットシフトを用いて信
号値を求める動作について説明する。図32は第3の実
施形態の動作を示すフローチャートである。なお、図3
2においては、上記第1および第2の実施形態と変動値
SA,SBの算出までの処理は同一であるため詳細な説
明は省略し、変動値SA,SBが算出されたものとし
て、重み付け係数Sjの算出および重み付け加算の処理
についてのみ説明する。また、図32においては、R4
4画素位置におけるB(R44)信号の算出を表すもの
である。まず、変動値SA,SBの和が0であるか否か
が判断される(ステップS41)。ステップS41が肯
定されると、SA=SB=0.5としてビットシフトを
用いて、式(63)によりB(R44)信号が算出され
る(ステップS42)。
【0150】 B(R44)=(B(R44)A+B(R44)B+1)>>1 (63) ステップS41が否定されると、SA≧7SBであるか
否かが判断される(ステップS43)。ステップS43
が肯定されると、変動値SA,SBの関係は図30の領
域A5にあることからSj=1とされ、下記の式(6
4)によりB(R44)信号が算出される(ステップS
44)。
【0151】 B(R44)=B(R44)B (64) ステップS43が否定されると、3SA≧5SBである
か否かが判断される(ステップS45)。ステップS4
5が肯定されると、変動値SA,SBの関係は図30の
領域A4にあることからSj=0.75とされ、下記の
式(65)によりB(R44)信号が算出される(ステ
ップS46)。
【0152】 B(R44)=(B(R44)A+3B(R44)B+2)>>2 (65) ステップS45が否定されると、5SA≧3SBである
か否かが判断される(ステップS47)。ステップS4
7が肯定されると、変動値SA,SBの関係は図30の
領域A3にあることからSj=0.5とされ、式(6
3)によりB(R44)信号が算出される(ステップS
48)。
【0153】ステップS47が否定されると、7SA≧
SBであるか否かが判断される(ステップS49)。ス
テップS49が肯定されると、変動値SA,SBの関係
は図30の領域A2にあることからSj=0.25とさ
れ、下記の式(66)によりB(R44)信号が算出さ
れる(ステップS50)。
【0154】 B(R44)=(3*B(R44)A+B(R44)B+2)>>2 (66) そして、ステップS49が否定されると、変動値SA,
SBの関係は図30の領域A1にあることからSj=0
とされ、下記の式(67)によりB(R44)信号が算
出される(ステップS51)。
【0155】B(R44)=B(R44)A (67) このように、第3の実施形態によれば、重み付け係数S
jを算出する際に除算を行わないため、重み付け係数S
jの算出を行う回路の構成を簡易なものとすることがで
き、その結果、本実施形態を実施するための装置の構成
を簡易なものとすることができるとともに、また演算時
間も短縮することができる。また、ビットシフトを行う
ことにより、整数値にて演算を行うことができるため、
演算を簡易に行って処理時間を短縮することができる。
したがって、本発明を実施するための装置をデジタルカ
メラ等の撮像装置やスキャナ等の画像読取装置、あるい
はプリンタ等に搭載することが容易となり、これにより
高機能の装置を作成することができる。
【0156】なお、上記各実施形態においては単板CC
D1として、G画素がR,B画素に対して倍の密度を有
するものについて説明したが、R画素がG,B画素に対
して、あるいはB画素がR,G画素に対して倍の密度を
有するものであってもよい。また、R,G,Bに対して
分光感度を有するものについて説明したが、単板CCD
1としてはこれに限定されるものではなく、Y(イエロ
ー)、G(グリーン)、C(シアン)に対して分光感度
を有するもの、あるいはY、W(ホワイト)、Cに対し
て分光感度を有するものであってもよい。
【0157】また、上記各実施形態においては、3色全
ての信号値を算出しているが、1色または2色のみの信
号値を算出するようにしてもよい。
【0158】さらに、上記第1および第2の実施形態に
おいては、式(41)により重み付け係数Sjを算出し
ているが、下記の式(68)により重み付け係数Sjを
算出してもよい。
【0159】 Sj=SB/(SA+SB) (if SA+SB=0 then Sj=0.5) (68) この場合、R、B信号QR、QBはは、下記の式(6
9)により算出される。
【0160】 QR,QB=Sj*Q1+(1-Sj)*Q2 (69)
【図面の簡単な説明】
【図1】本発明の第1の実施形態による画像処理装置の
構成を示す概略ブロック図
【図2】単板CCDにおける画素配列を示す図
【図3】ハニカム配列の単板CCDの画素配列を示す図
【図4】ハニカム配列の単板CCDの各画素位置を参照
番号を付して示す図
【図5】RB補間手段の構成を示す概略ブロック図
【図6】Cubicスプライン補間フィルタの例を示す図
【図7】信号値の算出を説明するための図(その1)
【図8】信号値の算出を説明するための図(その2)
【図9】信号値の算出を説明するための図(その3)
【図10】信号値の算出を説明するための図(その4)
【図11】信号値の算出を説明するための図(その5)
【図12】信号値の算出を説明するための図(その6)
【図13】信号値の算出を説明するための図(その7)
【図14】赤(R)および緑(G)からなる色エッジの
例を示す図
【図15】グレーエッジの例を示す図
【図16】変動値SA,SBの算出方法の差異を説明す
るための図
【図17】市松正方フィルタの例を示す図(その1)
【図18】市松正方フィルタの例を示す図(その2)
【図19】第1の実施形態の動作を示すフローチャート
【図20】信号算出処理の動作を示すフローチャート
【図21】本発明の第2の実施形態による画像処理装置
におけるRB補間手段の構成を示す概略ブロック図
【図22】緑(G)の単色平坦部の信号値を示す図
【図23】青(B)の単色平坦部の信号値を示す図
【図24】赤(R)の単色平坦部の信号値を示す図
【図25】第2の実施形態における信号算出処理の動作
を示すフローチャート
【図26】ハニカム配列の単板CCDを用いたデジタル
カメラの撮像系を示す図
【図27】ハニカム配列の単板CCDにおいて得られた
信号値の周波数特性を示す図
【図28】光学ローパスフィルタの配置を示す図
【図29】変動値SA,SBと重み付け係数Sjとの関
係を示す図
【図30】第3の実施形態の処理を説明するための図
【図31】本発明の第3の実施形態による画像処理装置
におけるRB補間手段の構成を示す概略ブロック図
【図32】第3の実施形態の動作を示すフローチャート
【図33】ハニカム配列を示す図
【図34】ベイヤー配列を示す図
【図35】ベイヤー配列の画素位置を参照番号を付して
示す図
【図36】ハニカム配列の画素位置を参照番号を付して
示す図
【符号の説明】
1,33 単板CCD 2 補間手段 3 G補間手段 4 RB補間手段 5 市松正方補間手段 10,21 判定手段 11 第1補間手段 12 第2補間手段 13 変動値算出手段 14 重み付け係数算出手段 14A メモリ 15 重み付け加算手段 22 線形補間手段 24 重み付け係数選択手段 31 撮影レンズ 32 光学ローパスフィルタ

Claims (17)

    【特許請求の範囲】
  1. 【請求項1】 異なる分光感度を有する第1から第3
    の光電変換素子を単一面上に配置して撮像面が形成され
    た撮像デバイスであって、前記第1および前記第2の光
    電変換素子を所定方向に交互に配置して第1ラインを形
    成し、前記第1および前記第3の光電変換素子を前記所
    定方向に交互に配置して第2ラインを形成し、前記所定
    方向と直交する直交方向において前記第1および前記第
    2の光電変換素子が交互に、かつ前記第1および前記第
    3の光電変換素子が交互に配置されるよう、前記第1お
    よび前記第2ラインを前記直交方向に交互に配置するこ
    とにより、第1から第3画素からなる前記撮像面が形成
    された撮像デバイスにおいて得られた第1から第3信号
    値に基づいて、全画素位置における前記第1から第3信
    号値のうち少なくとも1つの信号値を推定する画像処理
    方法において、 信号値を推定する画素位置が、グレーエッジにあるか色
    エッジにあるかを判定し、 前記所定方向についての信号値の推定値である所定方向
    推定値および前記直交方向についての信号値の推定値で
    ある直交方向推定値を算出し、 前記判定の結果に応じた異なる演算方法により、前記所
    定方向についての信号値の変動量を表す所定方向変動値
    および前記直交方向についての信号値の変動量を表す直
    交方向変動値を算出し、 前記所定方向変動値および前記直交方向変動値に基づい
    て、前記所定方向推定値および前記直交方向推定値を重
    み付け加算する際の重み付け係数を設定し、 該重み付け係数に基づいて、前記所定方向推定値および
    前記直交方向推定値を重み付け加算して、前記第1から
    第3信号値のうち少なくとも1つの信号値を推定するこ
    とを特徴とする画像処理方法。
  2. 【請求項2】 前記所定方向変動値と前記直交方向変
    動値との加算値を算出し、該加算値に対する前記所定方
    向変動値または前記直交方向変動値の比率を算出し、該
    比率を前記所定方向推定値および前記直交方向推定値を
    重み付け加算する際の重み付け係数として設定すること
    を特徴とする請求項1記載の画像処理方法。
  3. 【請求項3】 前記所定方向変動値および前記直交方
    向変動値の加算値に対する前記所定方向変動値または前
    記直交方向変動値の比率を、前記所定方向推定値および
    前記直交方向推定値を重み付け加算する際の重み付け係
    数として設定するに際し、 前記重み付け係数が一定であるとした場合の前記所定方
    向変動値と前記直交方向変動値との関係に基づいて、複
    数の重み付け係数に対する所定方向変動値および直交方
    向変動値の対応関係を予め算出して記憶した記憶手段を
    参照して、前記所定方向変動値および前記直交方向変動
    値に基づいて、前記複数の重み付け係数から一の重み付
    け係数を選択し、該選択された一の重み付け係数を、前
    記所定方向推定値および前記直交方向推定値を重み付け
    加算する際の重み付け係数として設定することを特徴と
    する請求項1記載の画像処理方法。
  4. 【請求項4】 前記記憶手段に記憶された前記対応関
    係における前記複数の重み付け係数が、前記加算値に対
    する前記所定方向変動値または前記直交方向変動値の比
    率が1/2(nは自然数)となるように設定されてな
    り、 前記重み付け加算を、前記重み付け係数に基づくビット
    シフトにより行うことを特徴とする請求項3記載の画像
    処理方法。
  5. 【請求項5】 前記信号値を推定する画素位置が、前
    記第1から第3信号値のいずれか1つの信号値のみから
    なる単色平坦部にあるか否かを判定し、 前記単色平坦部にあると判定された場合には、該単色平
    坦部内における前記第1から第3信号値に基づく線形補
    間演算により、前記第1から第3信号値のうち少なくと
    も1つの信号値を推定し、 前記単色平坦部にないと判定された場合に、前記グレー
    エッジにあるか色エッジにあるか否かの判定、前記所定
    方向推定値および前記直交方向推定値の算出、前記所定
    方向変動値および前記直交方向変動値の算出、前記重み
    付け係数の設定、および該重み付け係数に基づく信号値
    の推定を行うことを特徴とする請求項1から4のいずれ
    か1項記載の画像処理方法。
  6. 【請求項6】 異なる分光感度を有する第1から第3
    の光電変換素子を単一面上に配置して撮像面が形成され
    た撮像デバイスであって、前記第1および前記第2の光
    電変換素子を所定方向に交互に配置して第1ラインを形
    成し、前記第1および前記第3の光電変換素子を前記所
    定方向に交互に配置して第2ラインを形成し、前記所定
    方向と直交する直交方向において前記第1および前記第
    2の光電変換素子が交互に、かつ前記第1および前記第
    3の光電変換素子が交互に配置されるよう、前記第1お
    よび前記第2ラインを前記直交方向に交互に配置するこ
    とにより、第1から第3画素からなる前記撮像面が形成
    された撮像デバイスにおいて得られた第1から第3信号
    値に基づいて、全画素位置における前記第1から第3信
    号値のうち少なくとも1つの信号値を推定する画像処理
    装置において、 信号値を推定する画素位置が、グレーエッジにあるか色
    エッジにあるかを判定する判定手段と、 前記所定方向についての信号値の推定値である所定方向
    推定値および前記直交方向についての信号値の推定値で
    ある直交方向推定値を算出する推定値算出手段と、 前記判定の結果に応じた異なる演算方法により、前記所
    定方向についての信号値の変動量を表す所定方向変動値
    および前記直交方向についての信号値の変動量を表す直
    交方向変動値を算出する変動値算出手段と、 前記所定方向変動値および前記直交方向変動値に基づい
    て、前記所定方向推定値および前記直交方向推定値を重
    み付け加算する際の重み付け係数を設定する重み付け係
    数設定手段と、 該重み付け係数に基づいて、前記所定方向推定値および
    前記直交方向推定値を重み付け加算して、前記第1から
    第3信号値のうち少なくとも1つの信号値を推定する重
    み付け加算手段とを備えたことを特徴とする画像処理装
    置。
  7. 【請求項7】 前記重み付け係数設定手段は、前記所
    定方向変動値と前記直交方向変動値との加算値を算出
    し、該加算値に対する前記所定方向変動値または前記直
    交方向変動値の比率を算出する比率算出手段を備え、 該比率を、前記所定方向推定値および前記直交方向推定
    値を重み付け加算する際の重み付け係数として設定する
    手段であることを特徴とする請求項6記載の画像処理装
    置。
  8. 【請求項8】 前記重み付け係数設定手段は、前記所
    定方向変動値と前記直交方向変動値との加算値に対する
    前記所定方向変動値または前記直交方向変動値の比率
    を、前記所定方向推定値および前記直交方向推定値を重
    み付け加算する際の重み付け係数として設定するに際
    し、前記重み付け係数が一定であるとした場合の前記所
    定方向変動値と前記直交方向変動値との関係に基づい
    て、複数の重み付け係数に対する所定方向変動値および
    直交方向変動値の対応関係を予め算出して、該対応関係
    を記憶した記憶手段と、 前記所定方向変動値および前記直交方向変動値に基づい
    て、前記記憶手段を参照して前記複数の重み付け係数か
    ら一の重み付け係数を選択する重み付け係数選択手段と
    を備え、 該選択された一の重み付け係数を、前記所定方向推定値
    および前記直交方向推定値を重み付け加算する際の重み
    付け係数として設定する手段であることを特徴とする請
    求項6記載の画像処理装置。
  9. 【請求項9】 前記記憶手段に記憶された前記対応関
    係における前記複数の重み付け係数が、前記加算値に対
    する前記所定方向変動値または前記直交方向変動値の比
    率が1/2(nは自然数)となるように設定されてな
    り、 前記重み付け加算手段は、前記重み付け加算を、前記重
    み付け係数に基づくビットシフトにより行う手段である
    ことを特徴とする請求項8記載の画像処理装置。
  10. 【請求項10】 前記信号値を推定する画素位置が、
    前記第1から第3信号値のいずれか1つの信号値のみか
    らなる単色平坦部にあるか否かを判定する平坦部判定手
    段と、 前記単色平坦部にあると判定された場合に、該単色平坦
    部内における前記第1から第3信号値に基づく線形補間
    演算により、前記第1から第3信号値のうち少なくとも
    1つの信号値を推定する線形補間手段とをさらに備え、 前記単色平坦部にないと判定された場合に、前記判定手
    段によるグレーエッジにあるか色エッジにあるかの判
    定、前記推定値算出手段による前記所定方向推定値およ
    び前記直交方向推定値の算出、前記変動値算出手段によ
    る前記所定方向変動値および前記直交方向変動値の算
    出、前記重み付け係数設定手段による前記重み付け係数
    の設定、および前記重み付け加算手段による前記重み付
    け係数に基づく信号値の推定を行うことを特徴とする請
    求項6から9のいずれか1項記載の画像処理装置。
  11. 【請求項11】 異なる分光感度を有する第1から第
    3の光電変換素子を単一面上に配置して撮像面が形成さ
    れた撮像デバイスであって、前記第1および前記第2の
    光電変換素子を所定方向に交互に配置して第1ラインを
    形成し、前記第1および前記第3の光電変換素子を前記
    所定方向に交互に配置して第2ラインを形成し、前記所
    定方向と直交する直交方向において前記第1および前記
    第2の光電変換素子が交互に、かつ前記第1および前記
    第3の光電変換素子が交互に配置されるよう、前記第1
    および前記第2ラインを前記直交方向に交互に配置する
    ことにより、第1から第3画素からなる前記撮像面が形
    成された撮像デバイスにおいて得られた第1から第3信
    号値に基づいて、全画素位置における前記第1から第3
    信号値のうち少なくとも1つの信号値を推定する画像処
    理方法をコンピュータに実行させるためのプログラムを
    記録したコンピュータ読取り可能な記録媒体において、 前記プログラムは、信号値を推定する画素位置が、グレ
    ーエッジにあるか色エッジにあるかを判定する手順と、 前記所定方向についての信号値の推定値である所定方向
    推定値および前記直交方向についての信号値の推定値で
    ある直交方向推定値を算出する手順と、 前記判定の結果に応じた異なる演算方法により、前記所
    定方向についての信号値の変動量を表す所定方向変動値
    および前記直交方向についての信号値の変動量を表す直
    交方向変動値を算出する手順と、 前記所定方向変動値および前記直交方向変動値に基づい
    て、前記所定方向推定値および前記直交方向推定値を重
    み付け加算する際の重み付け係数を設定する手順と、 該重み付け係数に基づいて、前記所定方向推定値および
    前記直交方向推定値を重み付け加算して、前記第1から
    第3信号値のうち少なくとも1つの信号値を推定する手
    順とを有することを特徴とするコンピュータ読取り可能
    な記録媒体。
  12. 【請求項12】 前記重み付け係数を設定する手順
    は、前記所定方向変動値と前記直交方向変動値との加算
    値を算出し、該加算値に対する前記所定方向変動値また
    は前記直交方向変動値の比率を算出する手順を有し、 該比率を、前記所定方向推定値および前記直交方向推定
    値を重み付け加算する際の重み付け係数として設定する
    手順であることを特徴とする請求項11記載のコンピュ
    ータ読取り可能な記録媒体。
  13. 【請求項13】 前記重み付け係数を設定する手順
    は、前記所定方向変動値および前記直交方向変動値の加
    算値に対する前記所定方向変動値または前記直交方向変
    動値の比率を、前記所定方向推定値および前記直交方向
    推定値を重み付け加算する際の重み付け係数として設定
    するに際し、 前記重み付け係数が一定であるとした場合の前記所定方
    向変動値と前記直交方向変動値との関係に基づいて、複
    数の重み付け係数に対する所定方向変動値および直交方
    向変動値の対応関係を予め算出して記憶した記憶手段を
    参照して、前記所定方向変動値および前記直交方向変動
    値に基づいて、前記複数の重み付け係数から一の重み付
    け係数を選択する手順を有し、 該選択された一の重み付け係数を、前記所定方向推定値
    および前記直交方向推定値を重み付け加算する際の重み
    付け係数として設定する手順であることを特徴とする請
    求項11記載のコンピュータ読取り可能な記録媒体。
  14. 【請求項14】 前記記憶手段に記憶された前記対応
    関係における前記複数の重み付け係数が、前記加算値に
    対する前記所定方向変動値または前記直交方向変動値の
    比率が1/2(nは自然数)となるように設定されて
    なり、 前記重み付け加算を行う手順は、前記重み付け加算を、
    前記重み付け係数に基づくビットシフトにより行う手順
    であることを特徴とする請求項13記載のコンピュータ
    読取り可能な記録媒体。
  15. 【請求項15】 前記信号値を推定する画素位置が、
    前記第1から第3信号値のいずれか1つの信号値のみか
    らなる単色平坦部にあるか否かを判定する手順と、 前記単色平坦部にあると判定された場合には、該単色平
    坦部内における前記第1から第3信号値に基づく線形補
    間演算により、前記第1から第3信号値のうち少なくと
    も1つの信号値を推定する手順と、 前記単色平坦部にないと判定された場合に、前記グレー
    エッジにあるか色エッジにあるか否かの判定、前記所定
    方向推定値および前記直交方向推定値の算出、前記所定
    方向変動値および前記直交方向変動値の算出、前記重み
    付け係数の設定、および該重み付け係数に基づく信号値
    の推定を行う手順とをさらに有することを特徴とする請
    求項11から14のいずれか1項記載のコンピュータ読
    取り可能な記録媒体。
  16. 【請求項16】 請求項6から10のいずれか1項記
    載の画像処理装置を備えたことを特徴とする撮像装置。
  17. 【請求項17】 請求項6から10のいずれか1項記
    載の画像処理装置を備えたことを特徴とする画像読取装
    置。
JP2000104977A 2000-04-06 2000-04-06 画像処理方法および装置並びに記録媒体 Expired - Fee Related JP4302854B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2000104977A JP4302854B2 (ja) 2000-04-06 2000-04-06 画像処理方法および装置並びに記録媒体

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2000104977A JP4302854B2 (ja) 2000-04-06 2000-04-06 画像処理方法および装置並びに記録媒体

Publications (2)

Publication Number Publication Date
JP2001292454A true JP2001292454A (ja) 2001-10-19
JP4302854B2 JP4302854B2 (ja) 2009-07-29

Family

ID=18618431

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2000104977A Expired - Fee Related JP4302854B2 (ja) 2000-04-06 2000-04-06 画像処理方法および装置並びに記録媒体

Country Status (1)

Country Link
JP (1) JP4302854B2 (ja)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006135919A (ja) * 2004-11-09 2006-05-25 Eastman Kodak Co カラー撮像画像データの補間方法およびプログラム
KR100782812B1 (ko) 2005-06-15 2007-12-06 삼성전자주식회사 에지 적응적 컬러 보간 방법 및 장치
US7894110B2 (en) 2004-12-01 2011-02-22 Megachips Corporation Pixel interpolation method and image distinction method
CN111383552A (zh) * 2020-05-09 2020-07-07 广州视源电子科技股份有限公司 一种led灯珠及led显示屏

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103164702B (zh) * 2011-12-13 2016-09-28 北京慧眼智行科技有限公司 标记中心点提取方法、装置及图像处理系统

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006135919A (ja) * 2004-11-09 2006-05-25 Eastman Kodak Co カラー撮像画像データの補間方法およびプログラム
US7894110B2 (en) 2004-12-01 2011-02-22 Megachips Corporation Pixel interpolation method and image distinction method
KR100782812B1 (ko) 2005-06-15 2007-12-06 삼성전자주식회사 에지 적응적 컬러 보간 방법 및 장치
CN111383552A (zh) * 2020-05-09 2020-07-07 广州视源电子科技股份有限公司 一种led灯珠及led显示屏

Also Published As

Publication number Publication date
JP4302854B2 (ja) 2009-07-29

Similar Documents

Publication Publication Date Title
JP4162111B2 (ja) 画像処理方法および装置並びに記録媒体
KR101342806B1 (ko) 최종 디지털 컬러 이미지를 형성하는 방법
US7502505B2 (en) High-quality gradient-corrected linear interpolation for demosaicing of color images
JP4127638B2 (ja) イメージセンサ出力データ処理装置及び処理方法
US8314863B2 (en) Image processing device, image processing method, and program pertaining to image correction
US20080123997A1 (en) Providing a desired resolution color image
JP3771054B2 (ja) 画像処理装置及び画像処理方法
JP2000078597A (ja) デジタルカラ―イメ―ジ値の処理装置
US8320714B2 (en) Image processing apparatus, computer-readable recording medium for recording image processing program, and image processing method
US20080018750A1 (en) Method and apparatus for color interpolation in digital photographing device
KR20090050460A (ko) 이미지 처리 장치 및 방법
US7489822B2 (en) Image processing apparatus and method for detecting a direction of an edge in the vicinity of a pixel of interest and generating all color signals for each pixel by interpolation using color signals of a pixel of interest and its neighbor pixels, and a recording medium having a program recorded thereon for causing the apparatus to perform the method
JP4302855B2 (ja) 画像処理方法および装置並びに記録媒体
JP4302854B2 (ja) 画像処理方法および装置並びに記録媒体
US7986361B2 (en) Color processing apparatus and method thereof for generating an interpolated image
JP4333544B2 (ja) 画像処理方法、及び画像処理装置、半導体装置、電子機器、画像処理プログラム、並びにコンピュータ読み取り可能な記録媒体
KR101327790B1 (ko) 영상 보간 방법 및 장치
JP5484015B2 (ja) 撮像装置、撮像方法、及びプログラム
JP4239483B2 (ja) 画像処理方法、画像処理プログラム、画像処理装置
JP4122082B2 (ja) 信号処理装置およびその処理方法
JP3946913B2 (ja) 信号生成方法および装置並びに記録媒体
JP4239480B2 (ja) 画像処理方法、画像処理プログラム、画像処理装置
US20050031222A1 (en) Filter kernel generation by treating algorithms as block-shift invariant
JP2003087808A (ja) 画像処理方法および装置並びにプログラム
JP2007325199A (ja) 撮像装置および信号処理装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20050912

A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A712

Effective date: 20061201

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20080201

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20080318

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20080516

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: 20090421

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20090423

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120501

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130501

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140501

Year of fee payment: 5

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees