(57)【要約】【課題】  明るい無彩色を撮像した際に、高輝度部分に色がつくことを防止する。【解決手段】  本発明の画像処理装置は、飽和線分３２０、３３０上の信号値を無彩色軸３１０上の対応する信号にマッピングする。また、併せて、飽和線分３２０、３３０の周囲に存在する色信号を、マッピング前の連続性を保持しつつマッピングを行う。画像処理装置は、飽和線分の周囲の色信号のマッピングにかかる移動量を、飽和線分までの距離に反比例するように決定する。飽和線分３２０、３３０上の信号の彩度が、マッピング前後で急峻に変化しているため、このように周辺の信号も原信号の連続性を保持したマッピングを行うことにより、マッピング前の飽和線分３２０、３３０上の信号と、その周辺の信号との彩度の階調を高い精度で保持することができる。【選択図】  図３
2005318499.tif 000002
【特許請求の範囲】
【請求項１】
３色の色成分を表す色空間によって表現された画像データを処理する画像処理装置であって、
前記３色の色成分のうち、１色が飽和している色信号および２色が飽和している色信号を、任意の光源に基づき特定される無彩色軸にマッピングする第１のマッピング処理部と、
前記飽和している色信号の周囲の色信号を、前記無彩色軸方向へ、該色信号の階調を保持するようマッピングを行う第２のマッピング処理部とを備える画像処理装置。
【請求項２】
３色の色成分の表色座標系によって表された画像データを処理する画像処理装置であって、
前記３色の色成分のうち、第１の色成分が飽和する第１飽和点と、前記第１の色成分および第２の色成分が飽和する第２飽和点と、前記第１の色成分、前記第２の色成分および第３の色成分が飽和する第３飽和点とを特定する飽和点特定部と、
前記第１飽和点、前記第２飽和点、および、前記第３飽和点を連結して構成される飽和線分上の色信号を、任意の光源に基づき特定される無彩色軸上にマッピングする第１のマッピング処理部と、
前記飽和線分から所定の範囲の領域を、マッピング候補領域として設定する候補領域設定と、
前記飽和線分から前記マッピング候補領域内の色信号までの距離に応じて、マッピングに伴う変位量を変化させて、該マッピング候補領域内の色信号を、前記無彩色軸方向へマッピングする第２のマッピング処理部とを備える画像処理装置。
【請求項３】
前記変位量は、前記飽和線分からの距離に反比例して定められる請求項２記載の画像処理装置。
【請求項４】
請求項１ないし請求項３いずれか記載の画像処理装置であって、
前記任意の光源とは、前記画像データの撮像時の環境に基づき設定された光源である画像処理装置。
【請求項５】
請求項１ないし請求項３いずれか記載の画像処理装置であって、
更に、
前記無彩色軸を、指定された光源に基づき特定される無彩色軸である指定光源無彩色軸に設定する無彩色軸設定部を備える画像処理装置。
【請求項６】
請求項５記載の画像処理装置であって、
更に、
前記３色の各色成分値を調整するための調整値を設定する調整値設定部を備え、
前記無彩色軸設定部は、前記調整値に基づき、前記無彩色軸上の座標点の各色成分値を調整して、前記無彩色軸を前記指定光源無彩色軸と同一の軸に設定し、
前記第１のマッピング処理部および前記第２のマッピング処理部は、前記設定された無彩色軸に基づき、前記マッピング処理を実行する画像処理装置。
【請求項７】
請求項６記載の画像処理装置であって、
更に、
前記調整値を特定するための調整値特定情報を入力する入力部と、
前記調整値特定情報と、前記調整値とを対応付ける調整値管理情報を格納する管理情報格納部とを備え、
前記調整値設定部は、前記入力された調整値特定情報および前記調整値管理情報に基づき、前記調整値を設定する画像処理装置。
【請求項８】
請求項６記載の画像処理装置であって、
前記調整値設定部は、前記画像データを解析し、前記解析の結果に基づき前記調整値を設定する画像処理装置。
【請求項９】
請求項８記載の画像処理装置であって、
前記調整値設定部は、前記画像データの無彩色近傍の画素群における前記３色の色成分の各平均値を算出し、前記各平均値を等量とする前記調整値を設定する画像処理装置。
【請求項１０】
請求項８記載の画像処理装置であって、
更に、
前記調整値を特定するための情報である前記調整値特定情報と、前記調整値とを対応付ける調整値管理情報を格納する管理情報格納部とを備え、
前記調整値設定部は、前記解析により前記調整値特定情報を取得し、前記取得した調整値特定情報および前記調整値管理情報に基づき、前記調整値を設定する画像処理装置。
【請求項１１】
請求項１０記載の画像処理装置であって、
前記調整値特定情報は、予め、前記画像データに設定されている情報である画像処理装置。
【請求項１２】
３色の色成分を表す色空間で表現された画像データを画像処理装置により処理する画像処理方法であって、
（ａ）前記３色の色成分のうち、１色が飽和している色信号および２色が飽和している色信号を、任意の光源に基づき特定される無彩色軸にマッピングする工程と、
（ｂ）前記飽和している色信号の周囲の色信号を、前記無彩色軸方向へ、該色信号の彩度の階調を保持するようマッピングする工程とを備える画像処理方法。
【請求項１３】
３色の色信号の表色座標系において、明るい無彩色を撮像した際に画像処理装置が行う色抑圧処理方法であって、
（ａ）前記３色の色成分のうち、第１の色成分が飽和する第１飽和点と、前記第１の色成分および第２の色成分が飽和する第２飽和点と、前記第１の色成分、前記第２の色成分および第３の色成分が飽和する第３飽和点とを特定する工程と、
（ｂ）前記第１飽和点、前記第２飽和点、および、前記第３飽和点を連結して構成される飽和線分上の信号を、任意の光源に基づき特定される無彩色軸上にマッピングする工程と、
（ｃ）前記飽和線分から所定の範囲の領域を、マッピング候補領域として設定する工程と、
（ｄ）前記飽和線分から前記マッピング候補領域内の色信号までの距離に応じてマッピングに伴う変位量を変化させて、該マッピング候補領域内の色信号を、前記無彩色軸方向へマッピングする工程とを備える色抑制処理方法。
【請求項１４】
コンピュータに、画像処理を実行させるためのコンピュータプログラムであって、
前記３色の色成分のうち、１色が飽和している色信号および２色が飽和している色信号を、任意の光源に基づき特定される無彩色軸にマッピングする機能と、
前記飽和している色信号の周囲の色信号を、前記無彩色軸方向へ、該色信号の彩度の階調を保持するようマッピングする機能とをコンピュータに実行させるためのコンピュータプログラム。
【請求項１５】
コンピュータに、色抑圧処理を実行させるためのコンピュータプログラムであって、
前記３色の色成分のうち、第１の色成分が飽和する第１飽和点と、前記第１の色成分および第２の色成分が飽和する第２飽和点と、前記第１の色成分、前記第２の色成分および第３の色成分が飽和する第３飽和点とを特定する機能と、
前記第１飽和点、前記第２飽和点、および、前記第３飽和点を連結して構成される飽和線分上の信号を、任意の光源に基づき特定される無彩色軸上にマッピングする機能と、
前記飽和線分から所定の範囲の領域を、マッピング候補領域として設定する機能と、
前記飽和線分から前記マッピング候補領域内の色信号までの距離に応じてマッピングに伴う変位量を変化させて、該マッピング候補領域内の色信号を、前記無彩色軸方向へマッピングする機能とをコンピュータに実行させるためのコンピュータプログラム。
【請求項１６】
請求項１４または請求項１５記載のコンピュータプログラムを、コンピュータ読み取り可能に記録した記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理方法に関し、特に、撮像した画像の色抑制処理に関するものである。
【背景技術】
【０００２】
近年、ＣＣＤ（Charge Coupled Devices）などの撮像素子を使用して画像を撮影するデジタルカメラが普及している。撮像素子が有する透過フィルタの透過率は各色で異なり、各色の飽和する信号値が異なる。そのため、撮像した画像では、各色の信号値が同時ではなく順に飽和する結果、高輝度部分で、本来その部分には存在しない色が付く、いわゆる「偽色」と呼ばれる現象が発生し、正確な色が再現されないという問題があった。
【０００３】
かかる問題点を改善するために、特許文献１では、撮像された画像において、輝度情報を用いて、高輝度部分の色信号のゲインを抑圧することにより、彩度を低下させ、偽色の発生を抑制するという技術を開示している。
【０００４】
【特許文献１】特開２０００−３２４８２号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１記載の技術では、広い領域でゲインを抑圧するため、飽和していない色信号に対しても彩度を低下させることとなる。従って、輝度信号が一定以上の高輝度部分では、彩度が完全に抑圧され階調が保持されず、色とびが発生するという問題があった。
【０００６】
本発明はかかる課題を解決するためになされたものであり、偽色の発生を抑制すると共に、原信号の階調情報を保持することを目的とする。
【課題を解決するための手段】
【０００７】
上述した課題の少なくとも一部を解決するために、本発明は第１の構成として以下の構成をとることとした。すなわち、
３色の色成分を表す色空間で表現された画像データを処理する画像処理装置において、
前記３色の色成分のうち、１色が飽和している色信号および２色が飽和している色信号を、任意の光源に基づき特定される無彩色軸にマッピングする第１のマッピング処理部と、
前記飽和している色信号の周囲の色信号を、前記無彩色軸方向へ、該色信号の階調を保持するようマッピングする第２のマッピング処理部とを備えることを要旨とする。
【０００８】
本発明によれば、偽色の発生を抑制するとともに、原信号の彩度の階調を保持した階調再現を行うことができる。
【０００９】
本発明の第２の構成として、３色の色信号の表色座標系において、画像データを処理する画像処理装置において、
前記３色の色信号のうち、第１の色信号が飽和する点を第１飽和点、第２の色信号が飽和する点を第２飽和点、および、第３の色信号が飽和する点を第３飽和点を特定する飽和点特定部と、
前記第１飽和点、前記第２飽和点、および、前記第３飽和点を連結して構成される飽和線分上の信号を、任意の光源に基づき特定される無彩色軸上にマッピングする第１のマッピング処理部と、
前記飽和線分から所定の範囲の領域を、マッピング候補領域として設定する候補領域設定と、
前記飽和線分から前記マッピング候補領域内の色信号までの距離に応じてマッピングに伴う変位量を変化させて、該マッピング候補領域内の色信号を、前記無彩色軸方向へマッピングする第２のマッピング処理部とを備えることを要旨とする。
【００１０】
本発明の構成をとることとすれば、飽和した色信号だけでなく、その周囲の色信号についても、彩度の連続性を保持しつつマッピングを行うことができる。従って、明るい無彩色を撮像した場合に、高輝度部分に色がつく、いわゆる偽色の発生を抑制するとともに、原信号の彩度の階調を高い精度で再現することができる。
【００１１】
本発明の画像処理装置において、前記変位量は、前記飽和線分からの距離に反比例して定められることとしてもよい。こうすれば、簡易な構成で彩度の連続性を保持しつつマッピングを行うことができる。従って、原信号の彩度の階調を、高い精度で、簡易に、再現することが可能となる。
【００１２】
本発明の画像処理装置において、前記任意の光源とは、前記画像データの撮像時の環境に基づき設定された光源であることとしてもよい。こうすれば、撮影時の環境を反映した無彩色軸に基づき、マッピング処理を行うことができる。
【００１３】
本発明の画像処理装置において、
更に、
前記無彩色軸を、指定された光源に基づき特定される無彩色軸である指定光源無彩色軸に設定する無彩色軸設定部を備えることとしてもよい。
【００１４】
本発明によれば、光源によって異なる無彩色軸を考慮したマッピング処理を行うことができるため、種々の光源について、適切な右色抑制処理を行うことができる。
【００１５】
本発明の画像処理装置において、
更に、
前記３色の各色成分値を調整するための調整値を設定する調整値設定部を備え、
前記無彩色軸設定部は、前記調整値に基づき、前記無彩色軸上の座標点の各色成分値を調整して、前記無彩色軸を前記指定光源無彩色軸と同一の軸に設定し、
前記第１のマッピング処理部および前記第２のマッピング処理部は、前記設定された無彩色軸に基づき、前記マッピング処理を実行することとしてもよい。
【００１６】
本発明の構成によれば、無彩色軸上の座標点の各色成分値を同一となるよう調整し、調整された無彩色軸に基づきマッピング処理を行うことができる。従って、色度のバランスを考慮した階調再現を行うことができる。
【００１７】
本発明の画像処理装置において、
更に、
前記調整値を特定するための調整値特定情報を入力する入力部と、
前記調整値特定情報と、前記調整値とを対応付ける調整値管理情報を格納する管理情報格納部とを備え、
前記調整値設定部は、前記入力された調整値特定情報および前記調整値管理情報に基づき、前記調整値を設定することとしてもよい。
【００１８】
本発明によれば、予め、調整値管理情報を有しているため、調整値管理情報と、入力された調整値特定情報とから、簡易に、調整値を特定することができる。
【００１９】
本発明の画像処理装置において、前記調整値設定部は、前記画像データを解析し、前記解析の結果に基づき前記調整値を設定することとしてもよい。
【００２０】
本発明によれば、画像データごとに、適切に調整値を設定することができ、階調再現の精度を向上することができる。
【００２１】
本発明の画像処理装置において、
前記調整値設定部における前記解析は、前記画像データの無彩色近傍の画素群における前記３色の色成分の各平均値を算出し、前記各平均値を等量とする前記調整値を設定することとしてもよい。
【００２２】
本発明によれば、画像データの、無彩色の被写体が撮像されている領域の画素値に基づき調整値を設定するため、種々の光源環境下で撮影された場合にも、適切な階調再現を行うことができる。
【００２３】
本発明の画像処理装置において、
更に、
前記調整値を特定するための情報である前記調整値特定情報と、前記調整値とを対応付ける調整値管理情報を格納する管理情報格納部とを備え、
前記調整値設定部は、前記解析により前記調整値特定情報を取得し、前記取得した調整値特定情報および前記調整値管理情報に基づき、前記調整値を設定することとしてもよい。
【００２４】
本発明によれば、解析により取得した調整値特定情報と、予め格納している調整値管理情報とに基づき、調整値を設定するため、簡易に、かつ、高い精度で、調整値を設定することができる。
【００２５】
本発明の画像処理装置において、前記調整値特定情報は、予め、前記画像データに設定されている情報であることとしてもよい。
【００２６】
本発明によれば、画像データを参照することにより、簡易に調整値特定情報を取得することができ、画像処理装置の処理負荷を軽減することができる。
【００２７】
本発明において、上述した種々の態様は、適宜、組み合わせたり、一部を省略したりして適用することができる。また、本発明は、上述した画像処理装置としての構成の他に、画像処理装置により画像処理、色抑制処理を行う画像処理方法、色抑制処理方法、コンピュータに画像処理、色抑制処理を実行させるためのコンピュータプログラム、かかるコンピュータプログラムをコンピュータ読み取り可能に記録した記録媒体等としても構成できる。いずれの構成においても、上述した各態様を適宜適用可能である。コンピュータが読み取り可能な記録媒体としては、例えば、フレキシブルディスクや、ＣＤ−ＲＯＭ、ＤＶＤ−ＲＯＭ、光磁気ディスク、ＩＣカード、ハードディスク等種々の媒体を利用することが可能である。
【発明を実施するための最良の形態】
【００２８】
以下、本発明の実施の形態について、実施例に基づき、次の順序で説明する。
Ａ．第１実施例：
Ａ１．システム構成：
Ａ２．画像処理：
Ａ３．色抑制処理：
Ａ４．マッピング処理：
Ａ５．効果：
Ｂ．第２実施例：
Ｂ１．光源情報：
Ｂ２．機能ブロック：
Ｂ３．色抑制処理：
Ｂ４．第２実施例の変形例：
Ｃ．変形例：
【００２９】
Ａ．実施例：
Ａ１．システム構成：
図１は、本実施例の画像処理システム１０００のシステム構成を例示する説明図である。画像処理システム１０００は、画像処理装置１００と、デジタルカメラ２００とを含み構成されている。デジタルカメラ２００は、撮像素子として、ＣＣＤを備えており、画像を撮像する。撮像された画像はＲＡＷ形式で生成され、画像データ１０として、図示しないメモリカードに保存される。本実施例では、デジタルカメラ２００の撮像素子としてＣＣＤを備えることとしたが、例えば、ＣＭＯＳ（Complementary Metal Oxide Semiconductor）などとしてもよい。
【００３０】
画像処理装置１００は、一般的なコンピュータであり、デジタルカメラ２００から画像データ１０の入力を受け付けると、かかる画像データ１０に対して、色抑制処理を含む画像処理を行い、ビットマップ形式の画像データ２０に変換する。
【００３１】
図に併せて画像処理装置１００の機能ブロックを示した。画像処理装置１００は、主制御部１０１と、画像入力部１０２と、画像処理部１０３と、画像格納部１０４とから構成されている。画像格納部１０４は、画像処理装置１００のハードディスクドライブ１０７の所定の領域に構成された領域である。他の機能ブロックはソフトウェア的に構成され、主制御部１０１により制御される。各機能ブロックは、ハードウェア的に構成することとしてもよい。
【００３２】
画像入力部１０２は、画像処理対象となる画像データ１０の入力を受け付け、画像処理部１０３に受け渡す機能を奏する。画像処理部１０３は、受け渡された画像データに対して、色抑制処理を含む画像処理を行い、画像データ２０に変換する機能を備える。かかる画像処理の詳細は、後述する。主制御部１０１は、変換された画像データ２０を、画像処理部１０３から受け取り、画像格納部１０４に格納する。
【００３３】
Ａ２．画像処理：
図２は、本実施例における画像処理を説明するブロックチャートである。画像処理部１０３は、図示するように、機能ブロックとして、欠陥画素補正処理部１５０と、色補間処理部１５１と、色抑制処理部１５２と、色変換処理部１５３と、逆γ補正処理部１５４とから構成されている。画像処理部１０３による画像処理は、画像入力部１０２から画像データ１０を受け取ることにより開始される。
【００３４】
欠陥画素補正処理部１５０は、画像データ１０の入力を検出すると、画像データ１０に対して、欠陥画素を補正する処理を行う。欠陥画素とは、撮像素子を構成する画素に欠陥画素が含まれている場合に、撮像された画像データ１０の欠陥画素に対応する画素を指しており、かかる画素データには、不正確な信号情報が設定された状態となっている。欠陥画素補正処理部１５０は、欠陥画素の周囲の画素の信号情報に基づき、かかる欠陥画素の画素値を補正する。
【００３５】
次に、色補間処理部１５１は、欠陥補正後の画像データに対して、色補間処理を実行する。色補間処理は、画像データ１０を構成する各画素データのフィルタ色以外の色成分を算出し、画素データごとにＲ，Ｇ，Ｂの各色成分の画素値を有するように、欠けているデータを補完する処理である。色補間処理によって、デジタルカメラ２００の機器特性に依存した色空間におけるＲＧＢデータが生成される。本実施例では、撮像素子としてＣＣＤを用いているため、生成される色空間をＣＣＤＲＧＢ色空間と呼ぶこととする。
【００３６】
色抑制処理部１５２は、色補間処理後の画像データ１０に対して色抑制処理を実行する。色抑制処理は、Ｒ，Ｇ，Ｂの各色成分を変化させ、彩度を抑圧する処理である。色抑制処理によって、高輝度の無彩色を撮像した際に、各色成分のフィルタ透過率の相異による偽色の発生や色飛びが抑制される。かかる処理は、本発明の要部であり、詳細は後述する。
【００３７】
次に、色変換処理部１５３は、色抑制処理後の画像データ１０に対して、色変換処理を行う。色変換処理とは、画像データの色空間を、目的に応じて適切な色空間へ変換する処理である。本実施例では、ＣＣＤＲＧＢ色空間から、ＲＧＢ色空間へ変換する。
【００３８】
次に、逆γ補正処理部１５４は、色変換後の画像データ１０に対して、出力機器の特性を示すガンマ値を使用して逆ガンマ補正を実行する。
【００３９】
画像処理部１０３は、以上の処理を行い、ＲＡＷ形式の画像データ１０を、ビットマップ形式の画像データ２０に変換する。画像データ２０は、画像格納部１０４に格納される。
【００４０】
画像処理部１０３は、上述の処理と併せて、例えば、カラーバランス（ホワイトバランス）、コントラスト、色彩、シャープネスの補正等、一般的な画質調整処理を実行することとしてもよい。
【００４１】
Ａ３．色抑制処理：
図３は、色抑制処理を模式的に示す説明図である。図３（ａ）は、ＣＣＤＲＧＢ色空間を模式的に示した図である。ＣＣＤＲＧＢ色空間３００（以降、色空間３００と呼ぶこととする）とは、ＣＣＤにより撮像した際に表現する事のできる信号の集合を表している。色空間３００は、図示するように、色信号の３色の成分であるＲ成分、Ｇ成分、Ｂ成分を軸とした表色座標系において、３次元的に略立方体となるよう構成されている。色空間３００において、全ての色信号の値が「０」のとき、すなわち、彩度「０」、かつ、輝度「０」のときには「黒色（Ｂｋ）」となり、かかる表色座標系の原点を表す。無彩色軸３１０は、ユーザにより指定された光源に基づき特定される彩度「０」となる軸である。ユーザによる指定に限らず、例えば、画像データに含まれるＥｘｉｆ情報から特定した光源情報を用いて無彩色軸を設定することとしてもよい。
【００４２】
本実施例における撮像素子であるＣＣＤには、ＲＧＢの各色成分を透過するフィルタが配置されている。色成分ごとにフィルタを透過する透過率が異なるため、各色成分は、同量のときに飽和せず、別々に飽和する。本実施例の色空間３００では、白色などの明るい無彩色を撮像した際に、Ｇ成分、Ｂ成分、Ｒ成分の順に飽和するものとする。図示するように、Ｇ成分が飽和する点を飽和点Ｐ１、Ｂ成分が飽和する点を飽和点Ｐ２、Ｒ成分が飽和する点を飽和点Ｐ３とよぶこととする。飽和点Ｐ１およびＰ２を結ぶ線分を飽和線分３２０と呼ぶこととし、同様に、飽和点Ｐ２およびＰ３を結ぶ線分を飽和線分３３０と呼ぶこととする。
【００４３】
飽和線分３２０、３３０上の信号値を、無彩色軸３１０上の対応する信号値にマッピングする。具体的には、飽和線分３２０上の信号の各色成分比を、無彩色軸上の信号の色成分比と同一とすることにより実現する。飽和線分３２０は、Ｇ成分が飽和しているため、Ｇ成分の値を変化させ、３色の色成分比を無彩色軸上の信号が有する信号の色成分比と同一に設定する。飽和線分３３０は、Ｇ成分およびＢ成分が飽和しているため、同様に、Ｇ成分およびＢ成分を変化させ、各色成分比を無彩色軸上の信号が有する色成分比と同一に設定する。
【００４４】
図３（ｂ）は、図３（ａ）に示す色空間３００を斜め上方向から見た斜視図である。飽和線分３２０上の信号は、Ｇ成分を変化させるため、Ｇ軸方向に垂直に移動して無彩色軸３１０上にマッピングされる。飽和線分３３０上の信号は、Ｇ成分およびＢ成分を増加するため、Ｇ成分の変化分と、Ｂ成分の増加分との合成方向に移動して、無彩色軸３１０上にマッピングされる。本実施例におけるマッピングに伴う「移動量」とは、本発明の「変位量」を示すものとする。
【００４５】
本実施例の色抑制処理部１５２は、飽和線分３２０、３３０上の信号値と、かかる飽和線分３２０、３３０の周囲に存在する色信号についても、併せてマッピングする。色抑制処理部１５２は、飽和線分３２０、３３０の周囲に存在する色信号のマッピングにかかる移動量を、飽和線分までの距離に反比例するように決定する。すなわち、飽和線分３２０、３３０の周辺領域に存在する信号についても、無彩色軸３１０方向へ、原信号の連続性を保持したマッピングを行う。飽和線分３２０、３３０上の信号の彩度が、マッピング処理の前後で急峻に変化しているため、このように周辺の信号もマッピングすることにより、飽和線分３２０、３３０上の原信号と、その周辺の信号との彩度の階調を高い精度で保持することができる。以下、図４〜図９に、本実施例における色抑制処理の詳細を示す。
【００４６】
Ａ４．マッピング処理：
図４は、本実施例における色抑制処理を説明するフローチャートである。色抑制処理部１５２が行う処理である。図中右上に示すように、色空間３００において、任意の色信号値をＤ（ｒ，ｇ，ｂ）、色抑制処理を行った変換後の色信号値をＤ'（ｒ'，ｇ'，ｂ'）、無彩色軸３１０のＲ信号、Ｇ信号、Ｂ信号の信号比を、Ｒｇ，Ｇｇ，Ｂｇとし、また、Ｒ成分の飽和値をＲｓ、Ｇ成分の飽和値をＧｓ、Ｂ成分の飽和値をＢｓとすることとする。
【００４７】
色抑制処理部１５２は、まず、計算量を低減するために、色空間３００において、３色の信号値が等量のときに無彩色となるよう正規化を行う（ステップＳ１０）。かかる計算式は、以下に示すとおりである。
【００４８】
式（１）
ｒａ＝ｒ／Ｒｇ；
ｇａ＝ｇ／Ｇｇ；
ｂａ＝ｂ／Ｂｇ；
Ｒｓａ＝Ｒｓ／Ｒｇ；
Ｇｓａ＝Ｇｓ／Ｇｇ；
Ｂｓａ＝Ｂｓ／Ｂｇ；
【００４９】
図５は、色空間３００を正規化して得られる色空間４００である。図４のステップＳ１０における処理により得られる空間である。無彩色軸４１０は、色空間４００における彩度「０」となる信号値の集合を表している。飽和点Ｐ１では、３色の色成分が等量のときに、Ｇ成分が飽和する。図に一点鎖線で示す飽和線分４２０、４３０は、飽和線分３２０、３３０に対応する線分である。
【００５０】
図４に戻り、説明を続ける。色抑制処理部１５２は、飽和線分４３０の周囲の領域のうち、マッピング対象とすべき領域を選択する（ステップＳ１１）。色抑制処理部１５２は、式（２）を適用してかかる処理を実行する。
【００５１】
式（２）
ｒａ＋ｂａ＞２＊Ｂｓａ  ａｎｄ
ｒａ      ≦Ｒｓａ      ａｎｄ
ｂａ      ≦Ｂｓａ；
【００５２】
次に、色抑制処理部１５２は、式（２）により選択された領域内の色信号のＢ成分を伸張する（ステップＳ１２）。色抑制処理部１５２は、式（３）を適用してかかる処理を実行する。
【００５３】
式（３）
ｒａ＝ｒａ；
ｂａ＝ｒａ＋２＊ｂａ−２＊Ｂｓａ；
【００５４】
図６は、色抑制処理部１５２が行うマッピング処理を模式的に示す説明図であり、図４のステップＳ１１〜Ｓ１２において行われる処理である。図６（ａ）は、本実施例におけるＢｋを原点とするＢ−Ｒ座標系５００を示している。無彩色線４４０は、飽和線分４２０を延長した直線である。また、直線４５０は、無彩色線４４０が、和線分４３０に対して対称となる直線である。かかる直線４５０は、以下の式で表すことができる。
【００５５】
ｒａ＋ｂａ＝２Ｂｓａ；
【００５６】
すなわち、式（２）が表す領域は、図にハッチングを付した領域５１０であり、Ｒ，Ｂ成分が、Ｐ２−Ｐ３−Ｑ２'に含まれる値を有する信号からなる領域である。式（３）は、Ｇ成分に関係なく、Ｂ成分を変化させ、領域５１０内の任意の信号Ｄ１（ｒａ，ｂａ）を、Ｄ１'（ｒａ、ｂａ'）に移動する式である。
【００５７】
図６（ｂ）は、式（３）を行った後の色空間４００を示す。図示するように、Ｑ２'−Ｐ３−Ｐ２−Ｑ５−Ｑ６−Ｑ４'から構成されるマッピング候補領域８００内の信号は、矢印Ａ方向に移動し、Ｑ２−Ｐ３−Ｐ２−Ｑ５−Ｑ６−Ｑ４から構成される領域にマッピングされる。
【００５８】
図４に戻り、説明を続ける。色抑制処理部１５２は、次に、飽和線分４２０の周囲の領域のうち、マッピング対象とすべき領域を選択する（ステップＳ１３）。色抑制処理部１５２は、式（４）を適用してかかる処理を実行する。
【００５９】
式（４）
ｒａ＋ｇａ＞２＊Ｇｓａ  ａｎｄ
ｒａ−ｂａ≦０          ａｎｄ
ｒａ      ≦Ｒｓａ      ａｎｄ
ｇａ      ≦Ｇｓａ；
【００６０】
次に、色抑制処理部１５２は、式（４）により選択された領域内の色信号のＧ成分を伸張する（ステップＳ１７）。色抑制処理部１５２は、式（８）を適用してかかる処理を実行する。
【００６１】
式（８）
ｒａ'＝ｒａ；
ｇａ'＝ｒａ＋２＊ｇａ−２＊Ｇｓａ；
ｂａ'＝ｂａ；
【００６２】
図７は、本実施例におけるマッピング処理を模式的に示す説明図であり、図４のステップＳ１３、Ｓ１７において行われる処理を示している。図７（ａ）は、本実施例におけるＢｋを原点とするＢ−Ｒ座標系５００および、Ｒ−Ｇ座標系６００を示している。直線４６０は、以下の式により表される。
【００６３】
ｒａ＋ｇａ＝２＊Ｇｓａ；
【００６４】
すなわち、式（４）が表す領域は、図７（ａ）にハッチングで示す領域５２０および６１０であり、Ｒ，Ｂ成分が、Ｐ１−Ｐ２−Ｑ７に含まれる値を有する信号で、かつ、Ｒ，Ｇ成分が、Ｐ２−Ｑ７−Ｑ１０の領域に含まれる値を有する信号の集合からなる領域である。かかる領域は、図７（ｂ）において、ハッチングを付したマッピング候補領域８１０である。
【００６５】
図７（ｂ）に示すように、Ｐ１−Ｐ２−Ｑ７−Ｑ１０ら構成されるマッピング候補領域８１０内の信号は、矢印Ｂ方向に移動され、Ｐ１−Ｐ２−Ｑ７−Ｑ３から構成される領域にマッピングされる。
【００６６】
図４に戻り、説明を続ける。色抑制処理部１５２は、更に、飽和線分４３０の周囲の領域のうち、マッピング対象とすべき領域を選択する（ステップＳ１４）。色抑制処理部１５２は、式（５）を適用してかかる処理を実行する。
式（５）の判断を行う（ステップＳ１４）。
【００６７】
式（５）
ｇａ＋ｂａ＞２＊Ｇｓａ  ａｎｄ
ｒａ−ｂａ＞０          ａｎｄ
ｇａ      ≦Ｇｓａ      ａｎｄ
ｂａ      ≦Ｒｓａ；
【００６８】
次に、色抑制処理部１５２は、式（５）により選択された領域内の色信号のＧ成分を伸張する（ステップＳ１８）。色抑制処理部１５２は、式（６）を適用してかかる処理を実行する。
【００６９】
式（６）
ｒａ'＝ｒａ；
ｇａ'＝ｒａ＋２＊ｇａ−２＊Ｇｓａ；
ｂａ'＝ｂａ；
【００７０】
図８は、本実施例におけるマッピング処理を模式的に示す説明図であり、図４のステップＳ１４〜Ｓ１５において行われる処理を示している。図８（ａ）は、本実施例におけるＢｋを原点とするＢ−Ｒ座標系５００および、Ｂ−Ｇ座標系７００を示している。直線４８０は、以下の式により表される。
【００７１】
ｇａ＋ｂａ＝２＊Ｇｓａ；
【００７２】
すなわち、式（５）が表す領域は、図８（ａ）にハッチングで示す領域５３０および７１０であり、Ｒ，Ｂ成分が、Ｐ１−Ｑ２−Ｑ７に含まれる値を有する信号で、かつ、Ｇ，Ｂ成分が、Ｑ７−Ｑ２−Ｑ９の領域に含まれる値を有する信号の集合から構成される領域である。かかる領域は、図８（ｂ）において、ハッチングを付したマッピング候補領域８２０である。
【００７３】
図８（ｂ）に、図示するように、Ｐ１−Ｑ７−Ｑ２−Ｑ９から構成されるマッピング候補領域８２０内の信号は、Ｇ成分が変更され、矢印Ｃ方向に移動し、Ｐ１−Ｑ７−Ｑ２−Ｑ１から構成される領域にマッピングされる。
【００７４】
図４に戻り説明を続ける。色抑制処理部１５２は、マッピング候補領域以外の領域内の色信号をそのまま保持する（ステップＳ１６）。色抑制処理部１５２は、式（７）を適用してかかる処理を実行する。
【００７５】
式（７）
ｒａ'＝ｒａ；
ｇａ'＝ｇａ；
ｂａ'＝ｂａ；
【００７６】
そして、色抑制処理部１５２は、正規化した色空間４００から元の色空間３００に戻す処理を行う（ステップＳ１８）。色抑制処理部１５２は、式（９）を適用してかかる処理を実行する。
【００７７】
式（９）
ｒ'＝ｒａ'＊Ｒｇ；
ｇ'＝ｇａ'＊Ｇｇ；
ｂ'＝ｂａ'＊Ｂｇ；
【００７８】
以上の処理を行った色空間３００を図９ａおよび図９ｂに示した。図９ａは、マッピング処理後の色空間３００を、Ｒ成分側から見た状態を表す説明図である。図に太線で示す領域（Ｐ１−Ｑ７−Ｑ２−Ｑ１とＰ２−Ｑ２−Ｐ３−Ｑ６−Ｑ５−Ｑ４）が、飽和線分３３０上の信号およびその周囲の信号をマッピングした領域であり、図４のステップＳ１２，Ｓ１５によりマッピングされた領域である。図９ｂは、マッピング処理後の色空間３００をＢ成分側から見た状態を表す説明図である。図に太線で示す領域（Ｐ１−Ｐ２−Ｑ７−Ｑ３）が、飽和線分３２０上の信号およびその周囲の信号をマッピングした領域である。
【００７９】
色抑制処理部１５２は、かかる処理を行った画像データ１０を、色変換処理部１５３において変換処理を行うことができるよう圧縮して受け渡す。
【００８０】
Ａ５．効果：
以上説明した実施例の画像処理装置によれば、明るい無彩色を撮像した際にも、彩度の階調を保持しつつ、偽色の発生を抑制することができる。図１０は、本実施例における色抑制処理の効果を示す説明図である。図１０には、従来の技術におけるマッピング模式図９００と、本発明におけるマッピング模式図９１０とを示した。
【００８１】
マッピング模式図９００は広い領域で彩度を抑制するした結果を示している。図に示すように、輝度ｋ２以上の領域９０２（ハッチング部分）と高輝度領域９０１内の信号は全て、一定量、彩度が低下される。そのため、輝度ｋ１以上の高輝度領域９０１（グレー部分）の信号は全て彩度が「０」となり無彩色軸Ｙにマッピングされる。従って、彩度が急峻に変化し、色飛びを発生することがあった。
【００８２】
本発明では、マッピング模式図９１０に示すように、マッピングにかかる移動量を連続的に変化させている。従って、マッピング前後での彩度の階調を高い精度で保持して自然な階調再現を行うことができ、また、偽色の発生を抑制することができる。
【００８３】
また、本発明では、無彩色軸および飽和線分に基づいてマッピング候補領域を設定することとしているため、簡易に、かかる領域を設定することができ、計算負荷を軽減することが可能となる。
【００８４】
Ｂ．第２実施例：
第１実施例では、無彩色の被写体を撮像した際のＲ，Ｇ，Ｂの色成分比であるＲｇ，Ｇｇ，Ｂｇを用いて色抑制処理を実行した。第２実施例では、ホワイトバランス処理において使用するゲイン量を考慮したうえで、色抑制処理を行う。ホワイトバランス処理とは、光源が異なる場合に、撮像した画像に実際の被写体の色とは異なる色が付くことを防止する処理であり、各色成分値に対して、光源によって異なるゲイン量を乗算することにより実行される。
【００８５】
本実施例では、以下に示す処理を行う。すなわち、任意に指定する光源に対応して設定されているゲイン量を用いて、無彩色軸上の座標点の、Ｒ，Ｇ，Ｂの色成分値を同一とするよう調整を行い、無彩色軸を、任意に指定する光源に基づき特定される無彩色軸と同一の軸に設定する。設定した無彩色軸に基づき色抑制処理を行う。ゲイン量は、光源によって異なる値が設定されている。
【００８６】
本実施例において、デジタルカメラ２００は、ユーザの入力に基づき、撮影時の光源を特定するための光源情報を画像データ３０に設定する機能を備える。画像処理装置１００は、画像データ３０に設定された環境条件に基づき、ホワイトバランス処理を行う。環境条件の詳細について、図１１を用いて説明する。
【００８７】
Ｂ１．光源情報：
図１１（ａ）〜（ｃ）は、本実施例における環境条件を説明する説明図である。図１１（ａ）は、撮影時における環境条件の設定方法を説明する説明図である。デジタルカメラ２００の背面には、液晶パネルＤＳＰと操作部２１０とが備えられている。操作部２１０は、操作ボタン２０１，２０２，２０３，２０４，２０５から構成される。デジタルカメラ２００を使用して撮影を行うユーザは、撮影前に操作ボタン２０１〜２０５を操作して液晶パネルＤＳＰに環境条件の設定メニューを表示する。図の液晶パネルＤＳＰに環境条件の設定メニューを示した。
【００８８】
図示するように、環境条件として「晴天」、「曇天」等が表示され、各環境条件と共にラジオボタン２２１，２２２，２２３，２２４が表示される。ユーザは、操作ボタン２０１〜２０５を操作して所望の環境条件を選択する。本実施例では、環境条件として「晴天」を選択することとした。
【００８９】
図１１（ｂ）は、画像データ３０のデータ構成を例示する説明図である。画像データ３０は、Ｅｘｉｆ情報３１と、データ３２とから構成されている。撮影時に設定した環境条件は、Ｅｘｉｆ情報３１に設定されている。Ｅｘｉｆ情報の内容の一例を、図１１（ｃ）を用いて説明する。
【００９０】
図１１（ｃ）は、Ｅｘｉｆ情報３１の設定内容を例示する説明図である。Ｅｘｉｆ情報３１には、図示するように、「画像タイトル」、「撮影日時」などの画像データ３０に関する情報とともに、「光源情報」「色空間情報」、「ＩＳＯ感度」、「露出時間」、「被写体距離」、「焦点距離」等の画像データ３０の撮影時の撮影環境が記録されている。図に破線で示すように、画像データ３０の「光源情報」は、「晴天」と設定されている。Ｅｘｉｆ情報３１は、簡易に参照することができ、種々の値を取得することができるため、光源情報の判断材料として好適である。
【００９１】
Ｂ２．機能ブロック：
図１２は、本実施例における画像処理の詳細を説明するブロックチャートである。色補間処理部１５１の後に、ゲイン量を設定する処理を行うゲイン量設定部１５５による処理が行われる。色抑制処理部１５２は、設定されたゲイン量を使用して無彩色軸の色成分値を調整し色抑制処理を行う。ゲイン量設定部１５５の詳細な構成を、図１３を用いて説明する。
【００９２】
図１３（ａ）は、本実施例におけるゲイン量設定部１５５の詳細な機能ブロックを説明する説明図である。ゲイン量設定部１５５は、ゲイン量特定情報取得部１６０と、ゲイン量特定部１６１とから構成される。
【００９３】
ゲイン量特定情報取得部１６０は、画像データ３０のＥｘｉｆ情報３１を参照して、「光源情報」を取得する。本実施例は、Ｅｘｉｆ情報３１の「光源情報」には、「晴天」と記録されている。
【００９４】
ゲイン量特定部１６１は、ゲイン表１７０を有している。ゲイン表１７０の詳細を図１３（ｂ）にしめした。図１３（ｂ）に示すように、ゲイン表１７０は、各種光源によって特定される各色信号のゲイン量を表している。ゲイン表５０は、「光源種」と、「Ｒｇａｉｎ」と、「Ｇｇａｉｎ」と、「Ｂｇａｉｎ」とから構成されている。「光源種」は、光源の種類を表す。「Ｒｇａｉｎ」はＲ信号のゲイン量を表し、「Ｇｇａｉｎ」はＧ信号のゲイン量を表し、「Ｂｇａｉｎ」は、Ｂ信号のゲイン量を表す。例えば、光源種が「白熱灯」の場合のゲイン量は、「Ｒｇａｉｎ」が「１．２０」、「Ｇｇａｉｎ」が「１．００」、「Ｂｇａｉｎ」が「２．４４」である。
【００９５】
ゲイン量特定部１６１は、ゲイン量特定情報取得部１６０から、「光源情報」を取得し、取得した「光源情報」に対応するゲイン量を特定する。本実施例では、ゲイン量特定部１６１は、ゲイン量特定情報取得部１６０から光源情報として「晴天」という情報を取得しているため、光源種が「晴天」であるゲイン量を用いる。図示するように、光源種が「晴天」と設定されている画像データに対して、ホワイトバランス処理を行う際のゲイン量は、「Ｒｇａｉｎ」は「１．９０」、「Ｇｇａｉｎ」は「１．００」、「Ｂｇａｉｎ」は「１．７１」である。このように特定したゲイン量を使用して、色抑制処理の処理内容を図１４を用いて説明する。
【００９６】
Ｂ３．色抑制処理：
図１４は、本実施例における色抑制処理を説明するフローチャートである。色抑制処理部１５２は、ゲイン量設定部１５５からゲイン量を取得することにより、本処理を開始する。
【００９７】
色抑制処理部１５２は、以下の式（１０）を適用して、無彩色軸上の座標点のＲ，Ｇ，Ｂの３色の色成分比を同一とする（ステップＳ２０）。
【００９８】
式（１０）
Ｒｇ＝１／Ｒｇａｉｎ；
Ｇｇ＝１／Ｇｇａｉｎ；
Ｂｇ＝１／Ｂｇａｉｎ；
【００９９】
色抑制処理部１５２は、式（１０）によって算出されたＲｇ，Ｇｇ，Ｂｇを使用して、色抑制処理を実行する（ステップＳ２１〜Ｓ２９）。ステップＳ２１〜ステップＳ２８までのしょりは、第１実施例の図４のステップＳ１０〜Ｓ１７間での処理と同様の処理であるため、詳細な説明を省略する。
【０１００】
色抑制処理部１５２は、以下の式（１１）を適用して、色抑制処理後の色空間の座標点を求める（ステップＳ２９）。
【０１０１】
式（１１）
ｒ'＝ｒａ'；
ｇ'＝ｇａ'；
ｂ'＝ｂａ'；
【０１０２】
以上説明した第２実施例の画像処理システムによれば、ホワイトバランス処理において使用されるゲイン量を用いて、無彩色軸上の座標点の３色の色成分値を同一となるよう調整し、調整された無彩色軸に基づいて色抑制処理を実行することができる。従って、種々の光源下で撮像された画像データにおいて、適切な無彩色軸に基づいて色抑制処理を行うことができ、精度の高い階調再現を行うことができる。
【０１０３】
また、本実施例では、画像データに記録されているＥｘｉｆ情報から光源情報を取得することとしているため、簡易に光源情報を取得でき、画像処理装置の処理負荷を軽減することができる。
【０１０４】
Ｂ４．第２実施例の変形例：
Ｂ４−１．変形例（１）：
本実施例では、画像データの撮像時設定した光源情報をＥｘｉｆ情報に記録
し、これを参照して光源情報を取得しゲイン量を設定することとしたがこれに限られない。例えば、画像処理装置１００の画像処理時に、任意に光源を設定することとしても良い。このような場合には、例えば、例えば、図２０に示すように、画像処理装置１００のディスプレイに、撮影条件をユーザに選択させる画面ＷＤを表示し、ユーザに選択させることとしてもよい。
【０１０５】
図１５は、本変形例における光源情報の設定画面を例示する画面例である。画面ＷＤは、光源情報選択部１８０と、決定ボタン１９０とキャンセルボタン１９１とから構成される。光源情報選択部１８０は、図示するように、光源情報として「晴天」、「曇天」等が表示され、各光源情報と共にラジオボタン１８１，１８２，１８３，１８４が表示される。ユーザは、所望する光源のラジオボタンを選択し、決定ボタン１９０を押下する。光源情報の選択処理をキャンセルする場合には、キャンセルボタン１９１を押下する。本図は、光源情報として「晴天」が選択された状態を示している。
【０１０６】
こうすれば、撮影時に光源情報を設定せずとも、後から簡易に光源情報を設定しゲイン量を設定することができるため、利便性の向上を図ることができる。
【０１０７】
Ｂ４−１．変形例（１）：
第２実施例では、Ｅｘｉｆ情報に記録された光源情報を参照し、予め有していたゲイン表１７０に基づきゲイン量を設定することとしたがこれに限られない。例えば、画像データの画素値を解析し、解析結果に基づきゲイン量を算出することとしてもよい。かかる場合には、以下のような処理により実現することができる。
【０１０８】
図１６は、本変形例におけるゲイン量算出処理のフローチャートである。ゲイン量設定部１５５が、色補間処理部１５１によって色補間された画像データを受け取ると、本処理は開始される。
【０１０９】
ゲイン量設定部１５５は、画像データ３０において、無彩色の被写体が撮像された付近の画素をサンプリングし（ステップＳ３０）、かかる画素群のＲ，Ｇ，Ｂの各色成分値の平均値を算出する（ステップＳ３１）。次に、ゲイン量設定部１５５は算出した各平均値を等量とする、Ｒ，Ｇ，Ｂ色成分についての各ゲイン量を算出する（ステップＳ３２）。
【０１１０】
以上の処理によって算出されたゲイン量を使用して、色抑制処理部１５２は、色抑制処理を実行する。
【０１１１】
本変形例によれば、画像データごとに、適切なゲイン量を柔軟に算出することができるため、階調再現の精度を向上することができる。
【０１１２】
Ｃ．変形例：
以上、本発明の種々の実施例について説明したが、本発明はこれらの実施例に限定されず、その趣旨を逸脱しない範囲で種々の構成をとることができることは言うまでもない。例えば、以下のような態様が挙げられる。
【０１１３】
Ｃ１．変形例１：
第１実施例では、Ｇ−Ｂ−Ｒの順に飽和することとしたがこれに限られない。どの順に飽和しても、本発明の色抑制処理方法によれば、自然な階調再現を行うことができる。
【０１１４】
Ｃ２．変形例２：
第１実施例では、３色の色信号の表色座標系として、ＲＧＢ表色座標系を使用することとしたがこれに限られない他の表色座標系を使用することとしてもよい。また、本実施例では、飽和線分および無彩色軸に基づき、マッピング候補とする領域を設定することとしたが、例えば、ユーザが任意に指定した領域としてもよい。また、飽和線分からの距離が一定以下の領域等としてもよい。
【図面の簡単な説明】
【０１１５】
【図１】第１実施例における画像処理システムのシステム構成を例示する説明図である。
【図２】第１実施例における画像処理を説明するブロックチャートである。
【図３】第１実施例における色抑制処理を模式的に示す説明図である。
【図４】第１実施例における色抑制処理を説明するフローチャートである。
【図５】第１実施例における色抑制処理中の色空間を例示する説明図である。
【図６】第１実施例におけるマッピング処理を模式的に示す説明図である。
【図７】第１実施例におけるマッピング処理を模式的に示す説明図である。
【図８】第１実施例におけるマッピング処理を模式的に示す説明図である。
【図９ａ】第１実施例における色抑制処理後の色空間を例示する説明図である。
【図９ｂ】第１実施例における色抑制処理後の色空間を例示する説明図である。
【図１０】第１実施例における色抑制処理の効果を示す説明図である。
【図１１】第２実施例における環境条件を説明する説明図である。
【図１２】第２実施例における画像処理の詳細を説明するブロックチャートである。
【図１３】第２実施例におけるゲイン量設定部１５５の詳細な機能ブロックを説明する説明図である。
【図１４】第２実施例における色抑制処理を説明するフローチャートである。
【図１５】変形例における光源情報の設定画面を例示する画面例である。
【図１６】変形例におけるゲイン量算出処理のフローチャートである。
【符号の説明】
【０１１６】
１０００...画像処理システム
１０，２０...画像データ
１００...画像処理装置
１０１...主制御部
１０２...画像入力部
１０３...画像処理部
１０４...画像格納部
１５０...欠陥画素補正処理部
１５１...色補間処理部
１５２...色抑制処理部
１５３...色変換処理部
１５４...逆γ補正処理部
２００...デジタルカメラ
３００...色空間
３１０...無彩色軸
３２０...飽和線分
３３０...飽和線分
４００...色空間
４１０...無彩色軸
４２０...飽和線分
４３０...飽和線分
４４０...無彩色線
８００，８１０，８２０...マッピング候補領域
９００，９１０...マッピング模式図
９０１，９０２、９１１...高輝度領域
【図１】
2005318499.tif 000003
【図２】
2005318499.tif 000004
【図３】
2005318499.tif 000005
【図４】
2005318499.tif 000006
【図５】
2005318499.tif 000007
【図６】
2005318499.tif 000008
【図７】
2005318499.tif 000009
【図８】
2005318499.tif 000010
【図９ａ】
2005318499.tif 000011
【図９ｂ】
2005318499.tif 000012
【図１０】
2005318499.tif 000013
【図１１】
2005318499.tif 000014
【図１２】
2005318499.tif 000015
【図１３】
2005318499.tif 000016
【図１４】
2005318499.tif 000017
【図１５】
2005318499.tif 000018
【図１６】
2005318499.tif 000019