JP2000113215A - 画像処理装置およびその処理を実行するためのプログラムを記録した記録媒体 - Google Patents

Abstract

(57)【要約】 【課題】 画像に含まれる鏡面反射領域を容易に決定す
る技術を提供する。 【解決手段】 画像処理装置は、画像に含まれる鏡面反
射領域を決定する領域決定部を備えている。領域決定部
は、画像に含まれる画素の画素値の等高線を示す等高線
データを生成する。次に、生成された等高線データで示
される等高線の凸部の頂上となる領域を、鏡面反射領域
の候補領域として選出する。候補領域が選出されると、
候補領域を囲む所定の複数の等高線を選択し、選択され
た複数の等高線でそれぞれ規定される等高線領域の面積
の変化率を検査する。複数の等高線領域の面積の変化率
が所定の範囲内にある場合には、候補領域を鏡面反射領
域として決定することができる。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】この発明は、画像処理の技術
に関し、特に、画像に含まれる鏡面反射領域を決定する
ための技術に関する。

【０００２】

【従来の技術】照明光を反射する物体を撮像した画像に
は、鏡面反射領域（光沢のある部分）が含まれている。
鏡面反射領域の色は物体の通常の色とは異なるので、画
像中の鏡面反射領域を自動的に認識できれば、これを種
々の画像処理に利用することができる。例えば、鏡面反
射領域の色を用いて、画像内の物体を照明する光源の色
（光源色）を推定することができ、推定された光源色を
画像の色順応処理などに利用できる。

【０００３】

【発明が解決しようとする課題】しかし、従来は、画像
中の鏡面反射領域を自動的に認識する技術は知られてお
らず、ユーザが、鏡面反射領域を指定しなければならな
かった。

【０００４】この発明は、従来技術における上述の課題
を解決するためになされたものであり、画像に含まれる
鏡面反射領域を容易に決定する技術を提供することを目
的とする。

【０００５】

【課題を解決するための手段およびその作用・効果】上
述の課題の少なくとも一部を解決するため、本発明の装
置は、画像処理装置であって、画像に含まれる鏡面反射
領域を決定する領域決定部を備え、前記領域決定部は、
前記画像に含まれる画素の画素値の等高線を示す等高線
データを生成する等高線データ生成部と、前記等高線デ
ータで示される等高線の凸部の頂上となる領域を、前記
鏡面反射領域の候補領域として選出する候補領域選出部
と、前記候補領域を囲む所定の複数の等高線を選択し、
選択された複数の等高線でそれぞれ規定される等高線領
域の面積の変化率を検査し、前記面積の変化率が所定の
範囲内にある場合に、前記候補領域を前記鏡面反射領域
として決定する鏡面反射領域決定部と、を備えることを
特徴とする。

【０００６】上記の画像処理装置においては、画像から
画素値の等高線を示す等高線データが生成される。この
等高線データを用いれば、候補領域選出部において鏡面
反射領域の候補となる候補領域を容易に選出することが
できる。さらに、候補領域とを囲む複数の等高線領域の
面積の変化率を解析すれば、鏡面反射領域の特徴を有す
る候補領域を特定することができるので、鏡面反射領域
を容易に決定することが可能となる。

【０００７】上記の装置において、前記面積の変化率の
前記所定の範囲は、２つの等高線でそれぞれ示される２
つの画素値が画素値のダイナミックレンジの１％減少す
る毎に、前記面積が約５％〜約１００％の間で増加する
範囲であることが好ましい。

【０００８】２つの等高線領域の面積の変化率が上記の
所定の範囲である場合には、その候補領域は鏡面反射領
域とみなすことができる。したがって、このような面積
の増加率を示す候補領域を検出すれば、確実に鏡面反射
領域を決定することが可能となる。

【０００９】また、上記の装置において、前記鏡面反射
領域決定部は、前記候補領域を含む複数の等高線領域の
うち、最も面積の小さい所定数の等高線領域の前記面積
の変化率を無視し、無視した等高線領域よりも面積の大
きい所定数の等高線領域についての前記面積の変化率を
検査するようにしてもよい。

【００１０】鏡面反射領域を含む複数の等高線領域のう
ち最も面積の小さい所定数の等高線領域の面積は大きく
変化する傾向がある。したがって、これらの面積変化率
を無視すれば、うまく鏡面反射領域を決定することが可
能となる。

【００１１】また、上記の装置において、前記鏡面反射
領域決定部は、さらに、前記候補領域における前記画像
の第１の色と、前記候補領域の周辺領域における前記画
像の第２の色とを比較し、前記第１の色と前記第２の色
との差が所定量以上の場合に、前記候補領域を前記鏡面
反射領域として決定することが好ましい。

【００１２】鏡面反射領域における画像の色は光源の色
に近く、鏡面反射領域とその周辺領域の色が近い色であ
る場合には、鏡面反射領域を用いて鏡面反射領域特有の
情報を得ることができない。したがって、上記のように
すれば、鏡面反射領域特有の情報が得られる候補領域を
鏡面反射領域として決定することが可能となる。

【００１３】さらに、上記の装置において、前記等高線
データを生成する前に前記画像を平滑化する平滑化処理
部を備えることが好ましい。

【００１４】このようにすれば、画像内のノイズが除去
されるので、等高線データをうまく生成することができ
る。

【００１５】また、上記の装置において、さらに、前記
画像が複数の色成分によって構成されている場合に、前
記複数の色成分を合成して単色の画像を生成する単色画
像生成部を備え、前記領域決定部は、前記単色画像を用
いて前記鏡面反射領域の決定を行うことが好ましい。

【００１６】このようにすれば、鏡面反射領域を迅速に
決定することが可能となる。

【００１７】さらに、上記の装置において、前記鏡面反
射領域における前記画像の所定の表色系の３刺激値と、
前記鏡面反射領域の周辺領域における前記画像の前記所
定の表色系の３刺激値とに基づいて、前記画像の撮像時
に用いられた光源の光源色を推定する光源色推定部と、
推定された前記光源色に応じて、前記画像の色順応処理
を行う色順応処理部と、を備え、前記色順応処理部は、
前記光源の光源色が経時的に変化する場合に、変化前に
推定された光源色と変化後に推定された光源色とを用い
て、時間経過に応じた色順応処理を行うようにしてもよ
い。

【００１８】上記のように推定された光源色の時間変化
に応じた色順応処理を行えば、実際に人間の目で見るよ
うな色順応処理された画像を得ることができる。

【００１９】本発明の記録媒体は、画像に含まれる鏡面
反射領域を決定するためのコンピュータプログラムを記
録したコンピュータ読み取り可能な記録媒体であって、
前記画像に含まれる画素の画素値の等高線を示す等高線
データを生成する機能と、前記等高線データで示される
等高線の凸部の頂上となる領域を、前記鏡面反射領域の
候補領域として選出する機能と、前記候補領域を囲む所
定の複数の等高線を選択し、選択された複数の等高線で
それぞれ規定される等高線領域の面積の変化率を検査
し、前記面積の変化率が所定の範囲内にある場合に、前
記候補領域を前記鏡面反射領域として決定する機能と、
をコンピュータに実現させるためのコンピュータプログ
ラムを記録したコンピュを記録する。

【００２０】このような記録媒体に記録されたコンピュ
ータプログラムがコンピュータシステムによって実行さ
れた場合にも、上記の画像処理装置と同様の作用・効果
を有し、画像に含まれる鏡面反射領域を容易に決定する
ことができる。

【００２１】

【発明の他の態様】この発明は、以下のような態様も含
んでいる。第１の態様は、画像処理方法であって、画像
に含まれる鏡面反射領域を決定する工程を備え、前記工
程は、前記画像に含まれる画素の画素値の等高線を示す
等高線データを生成する工程と、前記等高線データで示
される等高線の凸部の頂上となる領域を、前記鏡面反射
領域の候補領域として選出する工程と、前記候補領域を
囲む所定の複数の等高線を選択し、選択された複数の等
高線でそれぞれ規定される等高線領域の面積の変化率を
検査し、前記面積の変化率が所定の範囲内にある場合
に、前記候補領域を前記鏡面反射領域として決定する工
程と、を備えることを特徴とする。

【００２２】第２の態様は、コンピュータに上記の発明
の各工程または各部の機能を実行させるコンピュータプ
ログラムを通信経路を介して供給するプログラム供給装
置である。

【００２３】なお、上記の等高線領域の面積の変化率
は、単純に等高線領域の面積そのものから算出される変
化率のみならず、例えば、等高線領域の面積の平方根に
応じて算出される変化率でもよい。本明細書において、
面積の変化率とは、一般的に面積を関数として表現した
ものから算出される種々の変化率を含む広義の意味とす
る。

【００２４】

【発明の実施の形態】Ａ．第１実施例：次に、本発明の
実施の形態を実施例に基づき説明する。図１は、本発明
の第１実施例としての画像処理装置を示すブロック図で
ある。この装置は、画像から光源色を推定し、推定され
た光源色を用いて画像の色順応処理を行うためのもので
ある。この装置は、ＣＰＵ１００と、バスライン１０２
とを備えている。バスライン１０２には、画像データメ
モリ１１０と、表示部１２０と、ＲＯＭ１３０と、磁気
ディスク１４０とが接続されている。また、バスライン
１０２には、画像内に含まれる鏡面反射領域を決定する
ための領域決定部１５０と、画像内の物体を照明する光
源の光源色を推定する光源色推定部１６０と、推定され
た光源色に基づいて画像の色順応処理を行う色順応処理
部１７０とが接続されている。領域決定部１５０は、画
像調整部１５１と、等高線データ生成部１５２と、等高
線データリンク部１５３と、候補領域選出部１５４と、
鏡面反射領域決定部１５５とを備えている。この画像処
理装置は、伝送路に接続して一般的なネットワークシス
テムを構成するようにしてもよい。なお、本実施例にお
いて領域決定部１５０と光源色推定部１６０と色順応処
理部１７０との機能は、ソフトウェアで実現されてい
る。

【００２５】上記の領域決定部１５０と光源色推定部１
６０と色順応処理部１７０との機能を実現するコンピュ
ータプログラムは、フレキシブルディスクやＣＤ−ＲＯ
Ｍ等の、コンピュータ読み取り可能な記録媒体に記録さ
れた形態で提供される。コンピュータは、その記録媒体
からコンピュータプログラムを読み取って内部記憶装置
または外部記憶装置に転送する。あるいは、通信経路を
介してコンピュータにコンピュータプログラムを供給す
るようにしてもよい。コンピュータプログラムの機能を
実現するときには、内部記憶装置に格納されたコンピュ
ータプログラムがコンピュータのマイクロプロセッサに
よって実行される。また、記録媒体に記録されたコンピ
ュータプログラムをコンピュータが読み取って直接実行
するようにしてもよい。

【００２６】この明細書において、コンピュータとは、
ハードウェア装置とオペレーションシステムとを含む概
念であり、オペレーションシステムの制御の下で動作す
るハードウェア装置を意味している。また、オペレーシ
ョンシステムが不要でアプリケーションプログラム単独
でハードウェア装置を動作させるような場合には、その
ハードウェア装置自体がコンピュータに相当する。ハー
ドウェア装置は、ＣＰＵ等のマイクロプロセッサと、記
録媒体に記録されたコンピュータプログラムを読み取る
ための手段とを少なくとも備えている。コンピュータプ
ログラムは、このようなコンピュータに、上述の各部の
機能を実現させるプログラムコードを含んでいる。な
お、上述の機能の一部は、アプリケーションプログラム
でなく、オペレーションシステムによって実現されてい
ても良い。

【００２７】なお、この発明における「記録媒体」とし
ては、フレキシブルディスクやＣＤ−ＲＯＭ、光磁気デ
ィスク、ＩＣカード、ＲＯＭカートリッジ、パンチカー
ド、バーコードなどの符号が印刷された印刷物、コンピ
ュータの内部記憶装置（ＲＡＭやＲＯＭなどのメモリ）
および外部記憶装置等の、コンピュータが読取り可能な
種々の媒体を利用できる。

【００２８】図２は、画像の色順応処理の手順を示すフ
ローチャートである。ステップＳ１０１においては、図
１の画像調整部１５１が画像データの調整を行う。な
お、画像データは、磁気ディスク１４０（図１）等から
読み出され、画像データメモリ１１０に記憶されてい
る。本実施例において、画像データは、Ｒ，Ｇ，Ｂの３
色の画像データで構成されている。

【００２９】ステップＳ１０１では、まず、Ｒ，Ｇ，Ｂ
の３色で表された画像データを合成することにより、グ
レー画像、すなわち、単色の画像データを生成する。グ
レー画像は、例えば、数式１によって生成される。

【００３０】

【数１】

【００３１】ここで、ｄ_R ，ｄ_G ，ｄ_B は、Ｒ，Ｇ，Ｂ
の画像データのうち、それぞれ対応する位置（画素）に
おける画素値を示している。ｄは、グレー画像におい
て、Ｒ，Ｇ，Ｂの各画像データに対応する位置（画素）
の画素値を示している。ｄ_R ，ｄ_G ，ｄ_B に乗じられて
いる各係数は、グレー画像においてカラー画像の輝度を
うまく表現するための係数である。

【００３２】このように、複数の色成分によって構成さ
れている画像データを合成して、単色の画像データに生
成することにより、後述するステップＳ１０２〜Ｓ１０
４（図２）の処理を迅速に行うことが可能となる。

【００３３】また、ステップＳ１０１では、生成された
単色の画像データについて、平滑化フィルタを用いて平
滑化処理を行う。平滑化フィルタとしては、例えば、ｍ
×ｍ画素のサイズを有し、すべてのフィルタ要素が「１
／ｍ² 」である平均値フィルタなどを用いることができ
る。整数ｍとしては、例えば、３〜９の値が用いられ
る。このように画像データを平滑化処理することによ
り、画像データ内のノイズを除去することができる。こ
れにより、後述するステップＳ１０２の処理をうまく行
うことが可能となる。

【００３４】なお、本実施例においては、複数色の画像
データから単色の画像データを生成した後に、平滑化処
理を行っているが、この順序は逆にしてもよい。すなわ
ち、複数色の画像データのそれぞれを平滑化処理した後
に、単色の画像データを生成してもよい。また、本実施
例のステップＳ１０１における処理は、後述のステップ
Ｓ１０２〜Ｓ１０４の処理を迅速に行うための処理であ
るため、省略可能である。

【００３５】上記の説明から分かるように、本実施例の
画像調整部１５１は、本発明における平滑化処理部と単
色画像生成部との機能を有している。

【００３６】図２のステップＳ１０２においては、図１
の等高線データ生成部１５２が、画像データから等高線
データを生成する。

【００３７】図３は、画像データの等高線を示す説明図
である。図３（Ａ）は、ステップＳ１０１において調整
された画像データを示しており、この画像内には、球状
の物体Ｍ１と円錐の形状を有する物体Ｍ２とが配置され
てる。図３（Ａ）から分かるように、２つの物体Ｍ１，
Ｍ２を照明する光源は、画像に対し右手上方にある。図
３（Ｂ）には、図３（Ａ）の画像データについての等高
線が示されている。

【００３８】図４は、図３（Ａ）の物体Ｍ１についての
画素値分布と等高線との関係を示す説明図である。図４
（Ａ）には、図３（Ａ）の物体Ｍ１を横切る直線Ｐ−Ｑ
上の画素についての画素値分布が示されている。図中、
横軸はＰ−Ｑ間における画素の位置を示している。縦軸
は、各画素における画素値を示している。Ｐ−Ｑ間の画
素値分布が図４（Ａ）に示すような場合、物体Ｍ１内に
おいて、画素値が同じ値となる画素を結ぶ等高線は図４
（Ｂ）に示すようになる。図４（Ｂ）の等高線データ
は、図４（Ａ）に示す画素値分布を複数の閾値で２値化
し、さらに、その２値画像の輪郭をベクトル化すること
によって生成される。なお、本実施例においては、それ
ぞれの等高線データは後述するように関連付けられてい
る。

【００３９】図５は、種々の等高線の例を示す説明図で
ある。図５（Ａ）は、画素値のピークが１つだけ存在す
る場合の等高線の例を示している。図５（Ｂ）は、画素
値のピークが２つ存在する場合の等高線の例を示してい
る。なお、図５（Ａ），（Ｂ）の右側の画素値分布にお
ける符号Ｈ１，Ｈ２，Ｈ３…は２値化の際の閾値（画素
値）を示しており、大きい順に、Ｈ１，Ｈ２，Ｈ３…と
表されている。なお、この明細書では、各等高線を求め
る際に使用された閾値を、その等高線の「階層」と呼
ぶ。階層（閾値）の差分は均等であることが好ましい
が、必ずしも均等にする必要はなく、任意の差分で複数
の階層を設定してもよい。

【００４０】図６は、図５の各階層の等高線データを示
す説明図である。図６（Ａ−１）〜（Ａ−４），（Ｂ−
１）〜（Ｂ−４）は、それぞれ図５（Ａ），（Ｂ）に示
す各階層における等高線データの一部を示している。な
お、以下の説明では、各等高線データによって決定され
る等高線の内側の領域を「等高線領域」と呼ぶ。

【００４１】図６（Ａ−１），（Ａ−２）の第１および
第２の階層Ｈ１，Ｈ２には、それぞれ第１および第２の
等高線データＣＡ_H1，ＣＡ_H2が含まれている。第２の階
層Ｈ２に含まれる等高線データＣＡ_H2によって決定され
る第２の等高線領域は、第１の階層Ｈ１に含まれる等高
線データＣＡ_H1によって決定される第１の等高線領域よ
り大きく、第２の等高線領域内に第１の等高線領域が含
まれる。すなわち、第２の等高線領域は、第１の等高線
領域を囲む。このとき、第１の等高線データＣＡ_H1と第
２の等高線データＣＡ_H2とは互いに関連付けられる。同
様にして、第２の等高線データＣＡ_H2と図示しない第３
の等高線データＣＡ_H3とも互いに関連付けられる。この
ようにして、第１の等高線データＣＡ_H1から第９の等高
線データＣＡ_H9まで、順次、関連付けられていく。

【００４２】図６（Ｂ−１）では、図５（Ｂ）に示すよ
うに第１の階層Ｈ１に等高線がないので、第１の階層Ｈ
１には等高線データが存在しない。図６（Ｂ−２）に示
す第２の階層Ｈ２には、２つの等高線データＣＢ_H2a ，
ＣＢ_H2a が含まれており、図示しない第３の階層Ｈ３に
も、２つの等高線データＣＢ_H3a ，ＣＢ_H3b が含まれて
いる。このとき、上記のようにして、第２の階層Ｈ２に
含まれる２つの等高線データＣＢ_H2a ，ＣＢ_H2b は、そ
れぞれ第３の階層Ｈ３に含まれる２つの等高線データＣ
Ｂ_H3a ，ＣＢ_H3b と関連付けられる。同様に、第１ない
し第５の階層Ｈ１〜Ｈ５に２つずつ含まれる等高線デー
タは、それぞれ互いに関連付けられる。

【００４３】図６（Ｂ−４）に示す第６の階層Ｈ６に
は、１つの等高線データＣＢ_H6が含まれている。第６の
階層Ｈ６に含まれる等高線データＣＢ_H6によって決定さ
れる第６の等高線領域は、図６（Ｂ−３）の第５の階層
Ｈ５に含まれる２つの等高線データＣＢ_H5a ，ＣＢ_H5b
によって決定される２つの等高線領域のそれぞれを囲
む。このとき、第５の一方の等高線データＣＢ_H5a と第
６の等高線データＣＢ_H6とは、互いに関連付けられる。
また、第５の他方の等高線データＣＢ_H5b と第６の等高
線データＣＢ_H6とも、互いに関連付けられる。このよう
にして、各階層Ｈ１〜Ｈ９に含まれる等高線データは、
順次、関連付けられていく。

【００４４】なお、各階層に含まれる等高線データの関
連付けは、図１の等高線データリンク部１５３によって
行われる。

【００４５】このようにステップＳ１０２（図２）にお
いて等高線データが生成されると、ステップＳ１０３に
おいて鏡面反射領域の候補となる候補領域を選出する。
なお、候補領域は、図１の候補領域選出部１５４によっ
て選出される。

【００４６】候補領域の選出は、関連付けられた複数の
等高線データのうち、最も画素値の大きい階層に含まれ
る等高線データによって決定される等高線領域を選出す
ることによって行われる。換言すれば、等高線で表され
る画素値の分布において、凸部の頂上となる領域が鏡面
反射領域の候補領域として選出される。例えば、図５
（Ａ）のような画素値分布を有する場合には、図６
（Ａ）の第１の階層Ｈ１に含まれる第１の等高線データ
ＣＡ_H1によって決定される等高線領域を候補領域として
選出する。また、図５（Ｂ）のような画素値分布を有す
る場合には、図６（Ｂ）に示す第２の階層Ｈ２に含まれ
る２つの等高線データＣＢ_H2a ，ＣＢ_H2b によって決定
される等高線領域をそれぞれ候補領域として選出する。

【００４７】上記のように候補領域を選出する場合に
は、画素値があまり大きくない領域、すなわち、鏡面反
射領域となり得ない領域も候補領域として選出される恐
れがある。そこで、本実施例においては、所定の画素値
以上の階層に含まれている等高線データによって決定さ
れる等高線領域を候補領域として選出するようにしてい
る。例えば、画素値がＨ１以上となる階層に含まれる等
高線データによって決定される等高線領域を候補領域と
して選出する場合には、図６（Ａ−１）では、第１の階
層Ｈ１に含まれる第１の等高線データＣＡ_H1によって決
定される等高線領域が候補領域として選出される。一
方、図６（Ｂ−１）では、第１の階層Ｈ１には等高線デ
ータは含まれていないので候補領域は存在しない。ま
た、画素値がＨ２以上となる階層に含まれる等高線デー
タによって決定される等高線領域を候補領域として選出
する場合には、図６（Ａ−１），（Ａ−２）では、第１
および第２の階層Ｈ１，Ｈ２に等高線データが含まれて
いるが、最も画素値の大きい階層Ｈ１に含まれる第１の
等高線データＣＡ_H1によって決定される等高線領域が候
補領域として選出される。一方、図６（Ｂ−１），（Ｂ
−２）では、第１の階層Ｈ１に等高線データが含まれて
おらず、第２の階層Ｈ２に２つの等高線データＣ
Ｂ _H2a ，ＣＢ_H2b が含まれているので、２つの等高線デ
ータＣＢ_H2a ，ＣＢ_H2b によって決定される等高線領域
が候補領域として選出される。

【００４８】なお、候補領域を選出するための階層の所
定の画素値は、実験により経験的に求めることができ
る。本実施例においては、画像データの各画素値は０〜
２５５の値（８ビット）を取り、画素値「１７０」以上
となる階層に含まれる等高線領域が候補領域として選出
されるように設定されている。

【００４９】図３（Ｃ）には、図３（Ｂ）に示す等高線
データについて、上記のように選出された鏡面反射領域
の候補領域が図示されている。

【００５０】次に、ステップＳ１０４（図２）において
は、ステップＳ１０３において選出された候補領域の中
から鏡面反射領域を決定する。なお、鏡面反射領域は、
図１の鏡面反射領域決定部１５５によって決定される。
鏡面反射領域は、候補領域およびその周辺の複数の等高
線領域の面積の変化率から決定される。

【００５１】本実施例においては、候補領域を含む複数
の等高線領域を検査対象領域として、これらの検査対象
領域のうち、隣り合う２つの階層における２つの等高線
領域の面積の増加率を順次解析することによって鏡面反
射領域が決定される。

【００５２】図７は、図３（Ａ）に示す画像のうち、物
体Ｍ１についての等高線領域の面積増加率を示す説明図
である。図７の左欄は画素値（階層）を示している。図
７の中欄は、左欄の階層の等高線領域の面積（図７では
画素数）を示している。また、右欄は、隣接する２つの
階層に含まれる等高線領域の面積増加率を示している。
なお、面積増加率Ｒは、数式２によって決定されてい
る。

【００５３】

【数２】

【００５４】ここで、Ｓ_kは、階層「ｋ」における等高
線領域の面積を示しており、Ｓ_k+1は、階層「ｋ」と隣
り合い、階層「ｋ」より画素値の小さい階層「ｋ＋１」
における等高線領域の面積を示している。

【００５５】図７では、各階層の画素値間隔が「５」で
ある。前述のように、画素値「１７０」以上となる階層
の等高線領域を候補領域とする場合には、図７には画素
値「１７０」の等高線領域が存在するので、これが候補
領域となる。

【００５６】本実施例では、候補領域を含む複数の等高
線領域についての面積増加率の最初の２つを無視し、か
つ、１５％〜８０％の範囲の面積増加率が５つ以上存在
する場合に、その候補領域を鏡面反射領域として決定す
る。ここで、面積増加率の範囲を限定している理由は、
面積増加率が過度に大きくなったり過度に小さくなった
りする場合には、候補領域が真の鏡面反射領域とは考え
られないからである。例えば、面積増加率が過度に大き
いとき（例えば５００％のとき）には、候補領域は単に
なだらかな画像部分であると考えることができる。ま
た、面積増加率が過度に小さいとき（例えば１％のと
き）には、候補領域は画像のエッジ部分であると考える
ことができる。したがって、候補領域が真の鏡面反射領
域であるときには、等高線領域の面積増加率がある程度
の範囲内にあるはずである。この面積増加率の好ましい
範囲は、実験により経験的に決定されるものであり、本
実施例では１５％〜８０％の範囲に設定されている。面
積増加率の最初の２つの値を無視する理由は、候補領域
が真の鏡面反射領域であるときにも、最初の２つの面積
増加率は、好ましい範囲を超えてしまう場合があること
が経験的に知られているからである。なお、候補領域が
真の鏡面反射領域ではないときにも、等高線領域の面積
増加率がたまたま好ましい範囲に入る場合がある。この
ような偶然を避けるために、本実施例では、好ましい範
囲の面積増加率が５つ以上存在するときにのみ、候補領
域を鏡面反射領域として認識するようにしている。

【００５７】図７において、候補領域（画素値「１７
０」）の面積（画素数）は「５８」であり、候補領域を
囲む等高線領域の面積は、画素値が小さくなるにつれ
「７１２」、「１２４７」、「１８３２」…と大きくな
っている。このとき、隣り合う等高線領域の面積増加率
は、順に、「１１２８％」、「７５％」、「４７％」…
となっている。図７の候補領域とその次の等高線領域と
の間での面積増加率は「１１２８％」であり、かなり大
きくなっているが、この面積増加率「１１２８％」と次
の面積増加率「７５％」との２つの面積増加率は無視さ
れる。このとき、図７の例では、１５％〜８０％の範囲
を満たす面積増加率の数は７つあるため、図７の候補領
域は、鏡面反射領域として決定される。

【００５８】なお、鏡面反射領域の決定に際しては、上
記の条件以外の他の条件を用いてもよい。例えば、１５
％〜８０％の範囲を満たす面積増加率の数が７つ以上あ
る場合にのみ、候補領域を鏡面反射領域として決定する
ようにしてもよい。また、面積増加率の範囲としては、
等高線領域の階層が画素値のダイナミックレンジ（８ビ
ットの場合、「２５６」）の１％減少する毎に、等高線
領域の面積が約５％〜約１００％増加する範囲を使用し
てもよい。

【００５９】ところで、図２のステップＳ１０２〜Ｓ１
０４までの処理は、以下に説明するような前処理を行っ
た後に実行してもよい。図８は、図３（Ａ）の物体Ｍ１
についての画素値分布とエッジ強度との関係を示す説明
図である。図８（Ａ）には、図３（Ａ）に示す物体Ｍ１
を横切る直線Ｐ−Ｑ上の画素についての画素値分布が示
されており、図４（Ａ）の画素値分布と同じである。図
８（Ｂ）には、各画素におけるエッジ強度が示されてい
る。なお、図８（Ｂ）に示すエッジ強度は、周知の微分
フィルタを用いてフィルタ処理することにより容易に求
めることができる。微分フィルタとしては、例えば、ラ
プラシアンフィルタや、１次微分フィルタなどを用いる
ことができる。

【００６０】図８（Ｂ）に示すように、エッジ強度は、
物体Ｍ１の外縁の画素Ｅ１，Ｅ２の近傍において他の画
素に比べ大きな値となる。すなわち、エッジ強度を求め
ることにより、画像内の物体の輪郭を知ることができ
る。したがって、ステップＳ１０２における等高線デー
タの関連付けを、エッジが明らかとなった領域内（すな
わち、物体の輪郭内）においてのみ行うことが可能とな
る。このとき、ステップＳ１０３では、物体の輪郭内に
おいて関連付けられた等高線データのみを用いて候補領
域を選出すればよい。また、ステップＳ１０４では、物
体の輪郭内において関連付けられた等高線データのみを
用いて面積変化率を検査し、鏡面反射領域を決定すれば
よい。このようにすれば、図２のステップＳ１０２〜Ｓ
１０４の処理を比較的迅速に行うことが可能となる。

【００６１】上記のように鏡面反射領域が決定される
と、図２のステップＳ１０５において、光源色推定部１
６０（図１）によって光源色が推定される。図９は、撮
像される画像内に含まれ得る鏡面反射領域およびその周
辺領域を示す説明図である。なお、「周辺領域」とは、
鏡面反射領域を除いた鏡面反射領域の周辺の領域を意味
する。図１０は、図９において撮像された画像内の鏡面
反射領域およびその周辺領域のＣＩＥ−ＸＹＺ表色系に
おける３刺激値Ｘ，Ｙ，Ｚの分布を示す説明図である。
図１０では、図示の便宜上、Ｚ軸を省略している。

【００６２】図９に示す光源２００から射出された光の
うち、球体２１０上において入射角と反射角とが等しい
角度θで反射される鏡面反射光ＬｓについてのＣＩＥ−
ＸＹＺ表色系における３刺激値は、図１０の点（Ｘｓ，
Ｙｓ，Ｚｓ）である。すなわち、図９の鏡面反射領域Ａ
ｓにおける色が、図１０の３刺激値（Ｘｓ，Ｙｓ，Ｚ
ｓ）に対応する。図１０の点（Ｘｂ，Ｙｂ，Ｚｂ）は、
図９の鏡面反射領域Ａｓと同じ領域において、光源２０
０からの光が鏡面反射しない場合、すなわち、光源２０
０からの光が球体２１０上において完全拡散反射する場
合の３刺激値を示している。なお、この３刺激値（Ｘ
ｂ，Ｙｂ，Ｚｂ）は、物体の色を示している。図１０に
おいて、光源の色を示す光源色ベクトル（Ｘｄ，Ｙｄ，
Ｚｄ）は、次の数式３によって求められる。

【００６３】

【数３】

【００６４】ところで、図１０の３刺激値（Ｘｂ，Ｙ
ｂ，Ｚｂ）は、図９の周辺領域Ａｅにおける色からほぼ
正確に求めることができる。すなわち、図９の鏡面反射
領域の近傍の周辺領域Ａｅの色は、図１０における点
（Ｘｂ，Ｙｂ，Ｚｂ）の近傍の点（Ｘｅ，Ｙｅ，Ｚｅ）
に対応する。したがって、周辺領域Ａｅの色を平均して
３刺激値を求めれば、図１０の点（Ｘｂ，Ｙｂ，Ｚｂ）
の値をほぼ正確に求めることができる。

【００６５】図１１は、図３（Ａ）の物体Ｍ１について
の鏡面反射領域とその周辺領域との３刺激値を示す説明
図である。図１１では、物体Ｍ１についての鏡面反射領
域における色の３刺激値（Ｘｓ１，Ｙｓ１，Ｚｓ１）
と、鏡面反射領域の近傍の周辺領域における色の３刺激
値（Ｘｂ１，Ｙｂ１，Ｚｂ１）とが示されてる。同様に
して、図３（Ａ）の物体Ｍ２についても、鏡面反射領域
における色の３刺激値（Ｘｓ２，Ｙｓ２，Ｚｓ２）と、
鏡面反射領域の近傍の周辺領域における色の３刺激値
（Ｘｂ２，Ｙｂ２，Ｚｂ２）とが求められる。なお、鏡
面反射領域およびその周辺領域における各３刺激値は、
それぞれの領域に含まれる画素の画素値を平均して求め
られる値である。

【００６６】図１２は、図３（Ａ）の画像から求められ
る光源色ベクトルを示す説明図である。このように、図
３（Ａ）の画像からは、物体Ｍ１，Ｍ２における鏡面反
射領域およびその周辺領域から２つの光源色ベクトルが
求められる。１つの画像から複数の光源色ベクトルが求
められる場合には、数式４に従って、平均の光源色ベク
トル（Ｘｗ，Ｙｗ，Ｚｗ）を求める。

【００６７】

【数４】

【００６８】ここで、（Ｘｄｉ，Ｙｄｉ，Ｚｄｉ）は画
像内のｉ番目の鏡面反射領域を用いて求められる光源色
ベクトルを示しており、ｎは画像内に含まれる鏡面反射
領域の総数を示している。

【００６９】このように複数の光源色ベクトルを平均化
することにより、より正確な光源色を推定することが可
能となる。

【００７０】ところで、鏡面反射領域の色および周辺領
域の色の関係によっては、光源色を推定することが困難
となる場合があるため、ステップＳ１０５（図２）にお
ける処理においては、以下の処理を併用することが好ま
しい。

【００７１】図１３は、光源色ベクトルと物体色ベクト
ルとの関係を示す説明図である。図１３（Ａ）は、光源
色ベクトル（Ｘｄ，Ｙｄ，Ｚｄ）と物体色ベクトル（Ｘ
ｂ，Ｙｂ，Ｚｂ）との方向が大きく異なる場合を示して
いる。図１３（Ｂ）は、光源色ベクトル（Ｘｄ，Ｙｄ，
Ｚｄ）と物体色ベクトル（Ｘｂ，Ｙｂ，Ｚｂ）との方向
がほぼ同じであるがわずかに異なる場合を示している。
図１３（Ａ）に示すような場合には、前述のように鏡面
反射領域の近傍の周辺領域における色を平均化すること
により、光源色ベクトルの始点（Ｘｂ，Ｙｂ，Ｚｂ）を
容易に求めることができる。したがって、ほぼ正確な光
源色ベクトルを推定することが可能である。一方、図１
３（Ｂ）に示すような場合には、光源色ベクトルの始点
（Ｘｂ，Ｙｂ，Ｚｂ）を正確に求めるのが困難となり、
この場合には、光源色ベクトルをあまり正確に推定する
ことができない。

【００７２】したがって、物体色ベクトルと光源色ベク
トルとが図１３（Ａ）に示すような関係にある場合にの
み、候補領域を鏡面反射領域として決定することが好ま
しい。このことは、候補領域とその周辺領域の色を比較
することによって判別可能である。例えば、候補領域に
おけるＲ，Ｇ，Ｂの画素値の比と、周辺領域における
Ｒ，Ｇ，Ｂの画素値の比とを比較し、その比が異なると
判断できる場合にのみ、その候補領域を鏡面反射領域と
して決定すればよい。あるいは、光源色ベクトル（Ｘ
ｄ，Ｙｄ，Ｚｄ）と物体色ベクトル（Ｘｂ，Ｙｂ，Ｚ
ｂ）との内積から、２つのベクトルによって形成される
角度を調べることによっても判断することが可能であ
る。

【００７３】なお、光源色ベクトルと物体色ベクトルと
の方向がほぼ完全に一致する場合には、光源色ベクトル
の始点をうまく特定できない可能性があるが、光源色ベ
クトルの方向は正確であるため、これを光源色ベクトル
として用いるようにしてもよい。

【００７４】上記のようにして、光源色が推定される
と、図２のステップＳ１０６において、色順応処理部１
７０（図１）は、画像データの色順応処理を実行する。
色順応処理は、ステップＳ１０５において推定された光
源色を用いて、画像内の各画素の色（３刺激値）を人間
の目の色順応によって認識される色（３刺激値）に変換
することによって行われる。本実施例の色順応処理は、
ＣＩＥ規格「ＣＩＥＣＡＭ９７ｓ」に基づいて行われ
る。すなわち、推定された光源色（Ｘｗ，Ｙｗ，Ｚｗ）
を用いて、画像内の各画素のＣＩＥ−ＸＹＺ表色系にお
ける３刺激値（Ｘ，Ｙ，Ｚ）を、上記の規格に基づいて
色順応処理し、Ｈ，Ｑ，Ｍ値を求める。ここで、Ｈ値は
色相に関する値、Ｑ値は明るさに関する値、Ｍ値は色彩
に関する値である。次に、求められたＨ，Ｑ，Ｍ値を、
通常用いられる所定の変換式を用いることによって、Ｃ
ＩＥ−ＸＹＺ表色系における３刺激値（Ｘ’，Ｙ’，
Ｚ’）に変換する。このように変換された３刺激値
（Ｘ’，Ｙ’，Ｚ’）を、さらに、表示部１２０（図
１）に適した表色系における３刺激値（例えば、Ｒ，
Ｇ，Ｂ）に変換することによって、色順応処理された画
像を表現することが可能となる。

【００７５】Ｂ．第２実施例：図１４は、本発明の第２
実施例としての画像処理装置を示すブロック図である。
この装置（デジタルカメラ）は、図１に示す画像処理装
置とほぼ同じ構成であるが、図１の装置に加え、画像入
力部３００と、光源色設定部３１０と、光源色センサ３
２０とが備えられている。

【００７６】画像入力部３００は、ＣＣＤなどの受光素
子を有する撮像系を備えている。光源色設定部３１０
は、画像を撮像する環境に応じてユーザが任意に光源色
を設定するために用いられる。光源色センサ３２０は、
光源色を直接測定するために用いられる。光源色設定部
３１０と光源色センサ３２０は、撮像された画像に鏡面
反射領域が含まれていない場合や、含まれていても光源
色の推定が困難である場合に用いることができる。

【００７７】すわなち、本実施例においては、光源色推
定部１６０と、光源色設定部３１０と、光源色センサ３
２０とから、３つの光源色を求めることができるので、
３つの光源色のうち、少なくとも１つを用いることによ
って光源色を決定することができる。例えば、３つの光
源色のそれぞれを、信頼度に応じて加重平均して求めら
れた光源色を用いてもよい。

【００７８】画像入力部３００において得られた画像デ
ータは、画像データメモリ１１０に、順次、記憶され
る。この装置では、撮像された画像に対して、略リアル
タイムに前述の色順応処理を行うことが可能である。

【００７９】このとき、撮像される物体を照明する光源
の光源色が変化する場合には、時間経過に応じた色順応
処理を行うことができる。すなわち、物体を照明する光
源色が変化した場合には、通常、人間の目の色順応は、
数分程度の時間を要すると言われている。したがって、
時間経過を考慮した光源色を用いて色順応処理を行え
ば、実際に人間の目で見ているような画像を得ることが
できる。なお、以下の説明では、実際に撮像する際の照
明光の光源色を「撮像光源色」と呼び、時間経過を考慮
した光源色を「擬似光源色」と呼ぶ。

【００８０】図１５は、時間経過に応じた色順応処理を
模式的に示す説明図である。図１５（Ａ），（Ｂ）は、
それぞれ、撮像光源色と擬似光源色との時間変化を示し
ている。図１５（Ａ）に示すように、撮像光源色は、時
刻taにおいて、（Ｘｗ_t1，Ｙｗ_t1，Ｚｗ_t1）から（Ｘｗ
_t2，Ｙｗ_t2，Ｚｗ_t2）に変化している。このとき、擬似
光源色として、図１５（Ｂ）に示すような時刻taから時
刻tbの間で、光源色が（Ｘｗ_t1，Ｙｗ_t1，Ｚｗ_t1）から
（Ｘｗ_t2，Ｙｗ_t2，Ｚｗ_t2）に徐々に変化するものを用
いる。擬似光源色（Ｘｗ(t) ，Ｙｗ(t) ，Ｚｗ(t) ）
は、数式５で表される。

【００８１】

【数５】

【００８２】ここで、Ｋ(t) は、Ｋ(ta)＝１、Ｋ(tb)＝
０を満たし、時間t とともに単調に減少する関数であ
る。

【００８３】このように求められた擬似光源色を用いて
前述の色順応処理を行えば、人間の目の色順応をうまく
再現した画像を得ることができる。なお、本実施例にお
いては、数式５に示すように、ＣＩＥ−ＸＹＺ表色系に
おける３刺激値を直接用いて時間変化に応じた色順応処
理を行っているが、ＣＩＥＬＡＢ色空間における色差Δ
Ｅの時間変化率が一定となるように変化させてもよい。

【００８４】以上、説明したように、本発明において
は、画像データから等高線データを生成し、等高線デー
タを解析することによって、画像内の鏡面反射領域を決
定している。すなわち、等高線データから鏡面反射領域
の候補となる候補領域を選出し、候補領域を囲む複数の
等高線領域の面積の変化率から、鏡面反射領域の特徴を
有する候補領域を鏡面反射領域として決定している。こ
れにより、画像に含まれる鏡面反射領域を容易に決定す
ることが可能となる。

【００８５】なお、この発明は上記の実施例や実施形態
に限られるものではなく、その要旨を逸脱しない範囲に
おいて種々の態様において実施することが可能であり、
例えば以下のような変形も可能である。

【００８６】（１）上記実施例では、画像データは、
Ｒ，Ｇ，Ｂの３色の画像データで構成されているが、他
の色で構成された画像データを用いてもよい。すなわ
ち、画像データを構成する色は、通常、画像を撮像した
際の受光素子に設置された光学フィルタによって決定さ
れている。例えば、光学フィルタとして、Ｒ，Ｇ，Ｂの
３色のフィルタを用いれば、画像データはＲ，Ｇ，Ｂの
３色から構成される。また、ＣＩＥ−ＸＹＺ表色系にお
ける等色関数にほぼ従うような光学フィルタを用いる場
合には、画像データは、Ｘ，Ｙ，Ｚの３刺激値で表され
る画像データとなる。このような他の色で構成された画
像データを用いても、上述のように鏡面反射領域を決定
することが可能である。

【００８７】なお、Ｘ，Ｙ，Ｚの３色で構成された画像
データを用いる場合には、図２のステップＳ１０１にお
いて調整される単色の画像データとしては、各色の画素
値を平均したものを用いてもよいし、輝度を表す刺激値
Ｙのみをそのまま用いてもよい。

【００８８】（２）上記実施例においては、図２のステ
ップＳ１０２において生成された等高線データは、異な
る階層間で関連付けられているが、この等高線データの
関連付けは、省略してもよい。ただし、等高線データが
関連付けられている場合には、候補領域が決定されれ
ば、候補領域を囲む複数の等高線領域も自動的に決定さ
れるので、図２のステップＳ１０３およびステップＳ１
０４における処理を迅速に行うことができるという利点
がある。

【００８９】（３）上記実施例においては、等高線デー
タはベクトル化されたデータであるが、等高線データは
２値化データでもよい。一般に、等高線データは画素値
の等高線を表すデータであればよい。

【００９０】（４）上記実施例においては、数式２の等
高線領域の面積変化率に基づいて鏡面反射領域を決定し
ているが、次の数式６に従って求められる変化率Ｒ’に
基づいて鏡面反射領域を決定してもよい。

【００９１】

【数６】

【００９２】ここで、ｒ_k+1 ，ｒ_k は、それぞれ数式２
の等高線領域の面積Ｓ_k+1 ，Ｓ_k と同じ面積を有する円
を考えた場合の半径を示している。すなわち、半径ｒ
_k+1 ，ｒ_k は、それぞれ、（Ｓ_k+1 ／π）^1/2 ，（Ｓ_k
／π）^1/2 である。

【００９３】数式６を用いた場合にも、変化率Ｒ’の範
囲を予め所定の範囲に設定しておくことにより、鏡面反
射領域を決定することが可能である。すなわち、本明細
書における面積の変化率は、面積に関する値の変化率を
含む広い意味である。

【図面の簡単な説明】

【図１】本発明の第１実施例としての画像処理装置を示
すブロック図。

【図２】画像の色順応処理の手順を示すフローチャー
ト。

【図３】画像データの等高線を示す説明図。

【図４】図３（Ａ）の物体Ｍ１についての画素値分布と
等高線との関係を示す説明図。

【図５】種々の等高線の例を示す説明図。

【図６】図５の各階層の等高線データを示す説明図。

【図７】図３（Ａ）に示す画像のうち、物体Ｍ１につい
ての面積増加率を示す説明図。

【図８】図３（Ａ）の物体Ｍ１についての画素値分布と
エッジ強度との関係を示す説明図。

【図９】撮像される画像内に含まれ得る鏡面反射領域お
よびその周辺領域を示す説明図。

【図１０】図９において撮像された画像内の鏡面反射領
域および周辺領域のＣＩＥ−ＸＹＺ表色系における３刺
激値Ｘ，Ｙ，Ｚの分布を示す説明図。

【図１１】図３（Ａ）の物体Ｍ１についての鏡面反射領
域とその周辺領域との３刺激値を示す説明図。

【図１２】図３（Ａ）の画像から求められる光源色ベク
トルを示す説明図。

【図１３】光源色ベクトルと物体色ベクトルとの関係を
示す説明図。

【図１４】本発明の第２実施例としての画像処理装置を
示すブロック図。

【図１５】時間経過に応じた色順応処理を模式的に示す
説明図。

【符号の説明】

１００…ＣＰＵ １０２…バスライン １１０…画像データメモリ １２０…表示部 １３０…ＲＯＭ １４０…磁気ディスク １５０…領域決定部 １５１…画像調整部 １５２…等高線データ生成部 １５３…等高線データリンク部 １５４…候補領域選出部 １５５…鏡面反射領域決定部 １６０…光源色推定部 １７０…色順応処理部 ２００…光源 ２１０…球体 ３００…画像入力部 ３１０…光源色設定部 ３２０…光源色センサ

───────────────────────────────────────────────────── フロントページの続き (72)発明者 嶋津 茂昭 京都市上京区堀川通寺之内上る４丁目天神 北町１番地の１ 大日本スクリーン製造株 式会社内 Ｆターム(参考） 5B050 BA01 DA02 EA02 EA08 5B080 GA11

Claims (8)

【特許請求の範囲】
1. 【請求項１】 画像処理装置であって、 画像に含まれる鏡面反射領域を決定する領域決定部を備
   え、 前記領域決定部は、 前記画像に含まれる画素の画素値の等高線を示す等高線
   データを生成する等高線データ生成部と、 前記等高線データで示される等高線の凸部の頂上となる
   領域を、前記鏡面反射領域の候補領域として選出する候
   補領域選出部と、 前記候補領域を囲む所定の複数の等高線を選択し、選択
   された複数の等高線でそれぞれ規定される等高線領域の
   面積の変化率を検査し、前記面積の変化率が所定の範囲
   内にある場合に、前記候補領域を前記鏡面反射領域とし
   て決定する鏡面反射領域決定部と、を備えることを特徴
   とする画像処理装置。
2. 【請求項２】 請求項１記載の画像処理装置であって、 前記面積の変化率の前記所定の範囲は、２つの等高線で
   それぞれ示される２つの画素値が画素値のダイナミック
   レンジの１％減少する毎に、前記面積が約５％〜約１０
   ０％の間で増加する範囲である、画像処理装置。
3. 【請求項３】 請求項１または２記載の画像処理装置で
   あって、 前記鏡面反射領域決定部は、 前記候補領域を含む複数の等高線領域のうち、最も面積
   の小さい所定数の等高線領域の前記面積の変化率を無視
   し、無視した等高線領域よりも面積の大きい所定数の等
   高線領域についての前記面積の変化率を検査する、画像
   処理装置。
4. 【請求項４】 請求項２または３記載の画像処理装置で
   あって、 前記鏡面反射領域決定部は、さらに、前記候補領域にお
   ける前記画像の第１の色と、前記候補領域の周辺領域に
   おける前記画像の第２の色とを比較し、前記第１の色と
   前記第２の色との差が所定量以上の場合に、前記候補領
   域を前記鏡面反射領域として決定する、画像処理装置。
5. 【請求項５】 請求項１ないし４のいずれかに記載の画
   像処理装置であって、さらに、 前記等高線データを生成する前に前記画像を平滑化する
   平滑化処理部を備える、画像処理装置。
6. 【請求項６】 請求項１ないし５のいずれかに記載の画
   像処理装置であって、さらに、 前記画像が複数の色成分によって構成されている場合
   に、前記複数の色成分を合成して単色の画像を生成する
   単色画像生成部を備え、 前記領域決定部は、前記単色画像を用いて前記鏡面反射
   領域の決定を行う、画像処理装置。
7. 【請求項７】 請求項１ないし６のいずれかに記載の画
   像処理装置であって、さらに、 前記鏡面反射領域における前記画像の所定の表色系の３
   刺激値と、前記鏡面反射領域の周辺領域における前記画
   像の前記所定の表色系の３刺激値とに基づいて、前記画
   像の撮像時に用いられた光源の光源色を推定する光源色
   推定部と、 推定された前記光源色に応じて、前記画像の色順応処理
   を行う色順応処理部と、を備え、 前記色順応処理部は、 前記光源の光源色が経時的に変化する場合に、変化前に
   推定された光源色と変化後に推定された光源色とを用い
   て、時間経過に応じた色順応処理を行う、画像処理装
   置。
8. 【請求項８】 画像に含まれる鏡面反射領域を決定する
   ためのコンピュータプログラムを記録したコンピュータ
   読み取り可能な記録媒体であって、 前記画像に含まれる画素の画素値の等高線を示す等高線
   データを生成する機能と、 前記等高線データで示される等高線の凸部の頂上となる
   領域を、前記鏡面反射領域の候補領域として選出する機
   能と、 前記候補領域を囲む所定の複数の等高線を選択し、選択
   された複数の等高線でそれぞれ規定される等高線領域の
   面積の変化率を検査し、前記面積の変化率が所定の範囲
   内にある場合に、前記候補領域を前記鏡面反射領域とし
   て決定する機能と、をコンピュータに実現させるための
   コンピュータプログラムを記録したコンピュータ読み取
   り可能な記録媒体。