(57)【特許請求の範囲】
【請求項１】
再現対象物の上に形成する画像を示す画像データを印刷手段に出力する色変換装置であって、
前記再現対象物の高さを示す高さ情報を各画素が有する形状データを取得する形状データ取得手段と、
前記形状データに基づいて、前記画像の各画素の色を変換するための色域データを導出する色域導出手段と、
前記色域データを用いて前記画像の各画素の色を変換する色変換手段と、を有し、
前記色域導出手段は、前記形状データが示す形状が第一の形状である場合は、前記形状データが前記示す形状が前記第一の形状よりも振幅が小さい形状である第二の形状である場合よりも小さい色域を示す前記色域データを導出することを特徴とする色変換装置。
【請求項２】
前記形状データの解析結果に基づいて、前記印刷手段が前記再現対象物の上に前記画像を形成する際の、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出する縮小率導出手段をさらに有し、
前記色域導出手段は、前記印刷手段の色域を示す色域データを、前記縮小率に基づいて変更することによって前記色域データを導出することを特徴とする請求項１に記載の色変換装置。
【請求項３】
前記縮小率導出手段は、前記形状データに対して周波数分解を行った結果に基づいて、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出することを特徴とする請求項２に記載の色変換装置。
【請求項４】
前記縮小率導出手段は、前記形状データに対して周波数分解を行うことによって得られる周波数に基づいて、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出することを特徴とする請求項２又は請求項３に記載の色変換装置。
【請求項５】
前記縮小率導出手段は、形状の周波数と縮小率とを対応付けたテーブルを参照して、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出することを特徴とする請求項３又は請求項４に記載の色変換装置。
【請求項６】
前記縮小率導出手段は、形状の各周波数成分についての縮小率の線形和を求めることによって、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出することを特徴とする請求項５に記載の色変換装置。
【請求項７】
前記縮小率導出手段は、前記形状データに対して周波数分解を行うことによって得られる周波数と振幅とに基づいて、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出することを特徴とする請求項２又は請求項３に記載の色変換装置。
【請求項８】
前記縮小率導出手段は、形状の周波数と振幅と縮小率とをそれぞれ対応付けたテーブルを参照して、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出することを特徴とする請求項７に記載の色変換装置。
【請求項９】
前記縮小率導出手段は、縮小率の特性が異なる周波数成分ごとに形状の各周波数成分についての縮小率の線形和を求めることによって、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出することを特徴とする請求項８に記載の色変換装置。
【請求項１０】
前記色域導出手段は、縮小率の特性が低明度部の色域縮小に対応する周波数成分の縮小率に基づいて、前記印刷手段の色域を示す色域データを変更し、変更した色域データに対して、さらに、縮小率の特性が光量の減少による全体明度の低下に対応する周波数成分の縮小率に基づいて変更を行なうことによって、前記画像の各画素の色を変換するための色域データを導出することを特徴とする請求項９に記載の色変換装置。
【請求項１１】
前記縮小率導出手段は、最大の振幅を有する方向の周波数成分に基づいて、前記画像の各画素に対応する前記印刷手段の色域の縮小率を導出することを特徴とする請求項２乃至請求項１０のいずれか一項に記載の色変換装置。
【請求項１２】
前記高さ情報は、基準面から前記再現対象物の最表面までの高さを示すことを特徴とする請求項１乃至請求項１１のいずれか一項に記載の色変換装置。
【請求項１３】
各画素の色が変換された前記画像を表示する表示手段をさらに有し、
前記表示手段は、前記画像の各画素のうち、色域の変更前後の隣接画素との色差の比が所定の閾値より大きい画素を強調した画像を表示することを特徴とする請求項１乃至請求項１２のいずれか一項に記載の色変換装置。
【請求項１４】
前記閾値の変更をユーザから受け付ける第一の受け付け手段をさらに有し、
前記表示手段は、前記閾値の変更がされた場合、前記画像の各画素のうち、色域の変更前後の隣接画素との色差の比が変更後の閾値より大きい画素を強調した画像を表示することを特徴とする請求項１３に記載の色変換装置。
【請求項１５】
前記色域導出手段で導出された前記色域データが示す色域のうち、最大色域を前記色変換手段における色変換に用いる第一のモードか、最小色域を前記色変換手段における色変換に用いる第二のモードか、の設定をユーザから受け付ける第二の受け付け手段をさらに有し、
前記色変換手段は、
前記第二の受け付け手段が前記第一のモードの設定を受け付けた場合、前記画像の各画素の色を、前記色域導出手段で導出された各画素の色域データを用いてそれぞれ変換し、
前記第二の受け付け手段が前記第二のモードの設定を受け付けた場合、前記画像の各画素の色を、前記色域導出手段で導出された各画素の色域データのうち、共通の色域を示す色域データを用いてそれぞれ変換することを特徴とする請求項１３又は請求項１４に記載の色変換装置。
【請求項１６】
コンピュータを請求項１乃至請求項１５のいずれか一項に記載の色変換装置の各手段として機能させるためのプログラム。
【請求項１７】
再現対象物の上に形成する画像を示す画像データを印刷手段に出力する色変換方法であって、
前記再現対象物の高さを示す高さ情報を各画素が有する形状データを取得する取得ステップと、
前記形状データに基づいて、前記画像の各画素の色を変換するための色域データを導出する導出ステップと、
前記色域データを用いて前記画像の各画素の色を変換する色変換ステップと、を有し、
前記導出ステップにおいては、前記形状データが示す形状が第一の形状である場合は、前記形状データが前記示す形状が前記第一の形状よりも振幅が小さい形状である第二の形状である場合よりも小さい色域を示す前記色域データを導出することを特徴とする色変換方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像処理装置で用いられる色変換技術に関する。
【背景技術】
【０００２】
近年、隆起を有する印刷物を得る技術が広がり、その色合わせの要求が高まっている。このような印刷物では、印刷層の重ね方や凹凸形状により色域が拡大縮小（例えば、印刷層を重ねるほど色材濃度が向上し色域は拡大）することが知られている。ここでユーザは、入力色データに基づく印刷物を確認して所望の色域が得られているかを判断し、不足の場合には印刷層を重ねることで所望の色域を得ることがある。このような場合、印刷結果が所望の色域を得られているかを、印刷物の測定やユーザによる目視によって確認する必要がある。また、所望の色域が得られていない場合は色パラメータを変更して再度印刷を行い、上記の確認作業を繰り返し行うことになるため作業負荷が大きくなる。
【０００３】
特許文献１では、所望の色域が得られているかを確認するために、略平面状に印刷層を重ねた場合の印刷層毎の色域目安をユーザに提示する方法が提案されている。特許文献１では、入力色データと各印刷層に使用されるインク量との対応関係を規定したプロファイルを用いて、入力色データから各印刷層を印刷するためのインク量を決定する。この際、プロファイルには、各印刷層を重ねて印刷する毎に色域が拡大するようインク量が規定されているため、ユーザは試行錯誤を重ねることなく印刷層の重ね数から簡易に印刷層の数と色域の関係を確認することができる。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１３−１４０５７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
上述のように、略平面での色予測を可能とする技術は提案されているが、隆起印刷のように凹凸のある起伏面においては画素に応じて凹凸の形状が異なるので単純な色予測が成立しない。特許文献１の技術は略平面上に印刷層を重ねることを前提としているため、凹凸のある起伏面における色域目安をユーザに提示することができない。  このような凹凸のある起伏面での色域の変化に対応して印刷色を確認するには、凹凸のある起伏面それぞれについて色域を示す情報を保持しておく方法が考えられる。そして、入力された凹凸のある起伏面を示す情報と保持している色域の情報とを利用してその起伏面の印刷色を確認する方法が考えられる。このような場合、全ての凹凸のある起伏面の色域の情報に対応するプロファイルの保持が必須であり、メモリコスト増大やプロファイル作成のための大量印字／測定が必要という問題が発生する。
【課題を解決するための手段】
【０００６】
本発明に係る色変換装置は、再現対象物の上に形成する画像を示す画像データを印刷手段に出力する色変換装置であって、前記再現対象物の高さを示す高さ情報を各画素が有する形状データを取得する形状データ取得手段と、前記形状データに基づいて、前記画像の各画素の色を変換するための色域データを導出する色域導出手段と、前記色域データを用いて前記画像の各画素の色を変換する色変換手段と、を有し、前記色域導出手段は、前記形状データが示す形状が第一の形状である場合は、前記形状データが前記示す形状が前記第一の形状よりも振幅が小さい形状である第二の形状である場合よりも小さい色域を示す前記色域データを導出することを特徴とする。
【発明の効果】
【０００７】
本発明によれば、多数のプロファイル保持によるメモリコスト増大やプロファイルの作成にかかる作業を増大させることなく、凹凸のある起伏面に対応してユーザが簡易に印刷色の確認を行うことが可能となる。
【図面の簡単な説明】
【０００８】
【図１】凹凸のある起伏面による色域縮小を説明するための模式図である。
【図２】実施例１における色変換処理装置の構成図である。
【図３】実施例１における色変換処理装置の論理構成を示すブロック図である。
【図４】実施例１における色変換処理を示すフローチャート図である。
【図５】実施例１における色変換ＬＵＴの一例を示す図である。
【図６】実施例１におけるＲＧＢ座標とＬａｂ座標とを模式した図である。
【図７】実施例１における四面体補間の例を示す図である。
【図８】実施例１におけるユーザインターフェースの一例を示す図である。
【図９】実施例１におけるＵＩの操作に応じた処理を示すフローチャートである。
【図１０】実施例１における色域縮小率算出を示すフローチャート図である。
【図１１】実施例１における色域縮小率と周波数の関係を示す模式図である。
【図１２】実施襟１におけるウェーブレット変換における周波数分解結果を示す模式図である。
【図１３】実施例１における周波数の方向判定を示すフローチャート図である。
【図１４】実施例１における色域導出処理を示すフローチャート図である。
【図１５】実施例２における凹凸のある起伏面による色域縮小を説明するための模式図である。
【図１６】実施例２における周波数と振幅による色域特性の分布を示す模式図である。
【図１７】実施例２における色域縮小率算出を示すフローチャート図である。
【図１８】実施例２における色域縮小率と周波数の関係を示す模式図である。
【図１９】実施例２における色域算出処理を示すフローチャート図である。
【図２０】実施例３における色域縮小率の算出処理を示すフローチャート図である。
【発明を実施するための形態】
【０００９】
以下に凹凸のある起伏面によって色域が縮小する概念を図１の模式図を用いて説明する。図１（ａ）及び図１（ｂ）は、色計測の幾何条件を基材断面方向から観察した模式図である。また図１（ｃ）は図１（ａ）及び（ｂ）におけるＬＣ平面での色域分布を模式的に表す図である。基材上に凹凸形状層が形成され、凹凸形状層の上部にさらに色材層が形成されているものとして説明を行う。図１（ａ）では、光源から発する光は入射角４５°で進入し、色材上で反射する。この時、受光部には主として表面拡散成分が入射し、図１（ａ）のような凹凸形状層においては図１（ｃ）のＧａに相当する色域となる。次に図１（ｂ）では、表面拡散成分だけでなく凹凸形状により入射角４５°以外の光における鏡面反射成分が受光部に入射する。図１（ｂ）のような凹凸形状層においては、低明度部の色域が縮小した図１（ｃ）のＧｂに相当する色域となる。なお、上述した現象は、印字表層の法線ベクトルの分布が不均一である場合に発生することが知られている。この際、色域体積の縮小率は起伏面の周波数特性と起伏の高さ（振幅）にもよるが、数十μメートルの振幅において数％〜数十％程度色域が縮小する場合がある。
【００１０】
以下で説明する実施形態においては、このように凹凸形状層において色域が縮小する場合に、凹凸のある起伏面に対応してユーザが簡易に印刷色の確認を行うことを可能とする例を説明する。また、以下で説明する実施形態では、メモリコストを増大させず、またプロファイルの作成にかかる作業を増大させずに済む。
【００１１】
以下、添付の図面を参照して、本発明を好適な実施例に基づいて詳細に説明する。なお以下の実施例において示す構成は一例にすぎず、本発明は図示された構成に限定されるものではない。
【００１２】
［実施例１］
＜システム構成例＞
図２は本実施例における色変換処理装置のハードウェア構成の一例を示す図である。色変換処理装置は、ＣＰＵ２０１、メインメモリ２０２、ＨＤＤ２０３、汎用インタフェース２０４、ディスプレイ２０５、メインバス２０６、キーボードやマウス等の指示入力部２０７、外部記憶装置２０８、プリンタ２０９を備える。汎用インタフェース２０４は、指示入力部２０７や外部記憶装置２０８などをメインバス２０６に接続する。
【００１３】
例えば、ＣＰＵ２０１は、指示入力部２０７に対するユーザの指示により、ＨＤＤ２０３や外部記憶装置２０８に格納されている色変換アプリケーションを起動する。そして、色変換アプリケーションをメインメモリ２０２に展開すると共にディスプレイ２０５にユーザインターフェースを表示する。続いて、ＨＤＤ２０３や外部記憶装置２０８に格納されている各種データが、ＣＰＵ２０１からの指令に基づきメインバス２０６経由によりメインメモリ２０２に転送される。メインメモリ２０２に転送された各種データは、ＣＰＵ２０１からの指令により所定の演算処理が行われ、演算処理の結果がメインバス２０６経由によってディスプレイ２０５上に表示あるいは、ＨＤＤ２０３や外部記憶装置２０８に格納される。
【００１４】
＜処理の構成＞
図３は、色変換処理装置３００のソフトウェアの機能構成例を示すブロック図である。なお、図３に示す構成は、ＣＰＵ２０１が色変換アプリケーションを実行することにより実現される。以下では、ＣＰＵ２０１がＨＤＤ２０３に格納された各種ソフトウェア（コンピュータプログラム）を動作させることで実現する各種処理について述べる。
【００１５】
色変換処理装置３００は、画像データ取得部３０１と、形状データ取得部３０２と、凹凸形状解析部３０３と、色域導出部３０４と、色域縮小率データ保持部３０５と、色域特性データ保持部３０６と、色変換部３０７と、出力部３０８とを有する。
【００１６】
画像データ取得部３０１は、汎用Ｉ／Ｆ２０４を介してＨＤＤ２０３や外部記憶装置２０８から入力された画像データを取得する。形状データ取得部３０２は、画像データ取得部３０１で取得した画像データの各画素に対応する凹凸形状データを、汎用Ｉ／Ｆ２０４を介してＨＤＤ２０３や外部記憶装置２０８から取得する。
【００１７】
凹凸形状解析部３０３は、形状データ取得部３０２で取得された凹凸形状データから凹凸形状特性を解析する。そして、解析の結果得られた凹凸形状特性と色域縮小率データ保持部に保持されているデータとに基づいて、各画素の色域の縮小率を決定する。そして、凹凸形状解析部３０３は、決定した各画素の色域の縮小率を色域導出部３０４に送信する。ここで、本実施例においては、凹凸形状特性は周波数特性を用いるものとし、凹凸形状解析部３０３ではウェーブレット変換により周波数分解を行い、その結果抽出された周波数情報を凹凸形状特性とする。ウェーブレット変換は公知の離散ウェーブレット変換を用いることができる。また、色域縮小率データ保持部３０５に保持されているデータは、例えば凹凸形状特性と縮小率とを対応付けたルックアップテーブル（ＬＵＴ）である。
【００１８】
色域導出部３０４は、色域特性データ保持部３０６に保持されている色域特性データと、凹凸形状解析部３０３から送られる縮小率とを用いて各画素の色域データを算出する。色域特性データ保持部３０６に保持されている色域特性データは、例えばプリンタ２０９に対応する平坦部の色域を示すデータである。
【００１９】
色変換部３０７は、画像データ取得部３０１で取得された画像データの各画素について、色域導出部３０４で算出された各画素に対応する色域データを利用してガマットマッピングを行い、予測される色データである出力結果を出力部３０８へ出力する。出力部３０８に出力された出力結果はディスプレイ２０５のユーザインターフェース上に表示される。ユーザは、ディスプレイ２０５上に表示された色データの出力結果（色変換結果）を見て、適切な色変換が成されるのか否かを確認することができる。
【００２０】
＜色変換処理動作＞
図４は色変換アプリケーションにて実行される色変換処理のフローチャートである。図４に示す処理は、ＣＰＵ２０１がＨＤＤ２０３に格納された色変換アプリケーションプログラムをメインメモリ２０２に読み出して実行することにより実現される。
【００２１】
ステップＳ４０１で画像データ取得部３０１は、色変換処理装置３００に入力された画像データを取得する。
【００２２】
ステップＳ４０２で形状データ取得部３０２は、色変換処理装置３００に入力された凹凸形状データを取得する。この凹凸形状データは、ステップＳ４０１で取得した画像データに対応したデータであり、例えば各画素の値が基準面からの高さを示すデータである。
【００２３】
ステップＳ４０３からステップＳ４０６の処理は各画素ごとに行なわれる処理である。ステップＳ４０３で凹凸形状解析部３０３は、i番目の画素に対応する縮小率を、ステップＳ４０２で取得した形状データを解析して得られた結果に基づいて導出する。詳細については後述する。
【００２４】
ステップＳ４０４で色域導出部３０４は、色域特性データ保持部３０６に格納された色域特性データと、ステップＳ４０３で導出された縮小率とから、画素ｉに対応する色域データを導出する。詳細については後述するが、ステップＳ４０４の処理の結果、縮小率に応じた変更後の色域データが出力される。
【００２５】
ステップＳ４０５で色変換部３０７は、ステップＳ４０４で導出された、変更後の色域データを用いて、画素ｉのガマットマッピング（色変換処理）を行う。ここで、ガマットマッピングはＣＩＥＬＡＢ空間のＬａｂ色空間にて行うものとし、マッピング手法については測色的色再現を用いるものとする。入力画像データの各画素はＲＧＢの値を有するデータであるものとする。入力画像のi番目の画素値をＲiＧiＢiとする。色変換部３０７は、ＲiＧiＢiから、Ｌａｂ色空間の色座標値への変換を、色変換ＬＵＴを用いた四面体補間演算により行い、ｉ番目の画素値のＲＧＢ→Ｌａｂ変換を施す。色変換ＬＵＴはＲＧＢ色空間での格子点の色座標値と、格子点が再現するＬａｂ色空間の座標値との対応を記したルックアップテーブルである。図５は色変換ＬＵＴのデータ構造の一例を示す図である。データ構造の先頭には、Ｒ／Ｇ／Ｂ値のステップが記述され、この後に、各格子点に対応するＬａｂ色座標が、Ｒ、Ｇ、Ｂの順でネストされて記述される。このデータ構造をＲＧＢ色空間で模式的に表すと、図６の様に表される。また、四面体補間は、隣接する８つの格子点で囲まれる六面体領域を、図７の模式図に示す様な６つに分割された四面体領域毎に、補間演算を行う。例えば、色変換部３０７は、ｉ番目の画素のＲiＧiＢi値が含まれている六面体領域をまず探索する。そして、探索された六面体領域を図７のように分割した四面体領域のうちの、ＲＧＢ値が含まれている四面体を特定する。そして、その特定した四面体の頂点の値とＲiＧiＢi値との距離に応じて、その特定した四面体の頂点のＬａｂ値に補間演算を行なうことでｉ番目の画素のＬａｂ変換された値であるＬiａiｂi値を導出する。ＲiＧiＢi値が頂点に対応する値そのものであれば、色変換ＬＵＴを参照して得られた値に変換すればよい。
【００２６】
また、このステップＳ４０５では、色変換部３０７は色域特性データ保持部３０６に保持されている、プリンタ２０９に対応した色域特性データを用いてｉ番目の画素の色変換処理も行なう。すなわち、ステップＳ４０５では、ｉ番目の画素について、変更前の色域データと変更後の色域データとを用いた色変換処理後の値がそれぞれ求められる。なお、色変換部３０７での色変換処理自体はいずれも上記で説明したような処理を行なえばよく、このとき用いる色変換ＬＵＴが異なる。
【００２７】
次に、ステップＳ４０６で色変換部３０７は、色域の変更前後の隣接画素との色差の比が閾値より大きいか否かを式１により判定する。すなわち、色域の変更前のｉ番目の画素とｉ−１番目の画素との色差と、変更後のｉ番目の画素とｉ−１番目の画素との色差と、の比が閾値より大きいか否かを判定する。式１において、色域の変更前のｉ番目の画素の色変換結果をＬiａiｂiとし、色域の変更後のｉ番目の画素の色変換結果をＬi!ａi!ｂi!とする。
【００２８】
【数１】
0006512884.tif 000002
【００２９】
本実施例では、画素ごとに凹凸形状に応じた色域を用いて色変換処理を行なう。このとき、凹凸形状によっては色域の変更前の色と変更後の色とが大きく異なる画素が生じる可能性がある。そこで、後述するように、ディスプレイ２０５に表示されるユーザインターフェースにおいてユーザに対して色域の変更前の色と変更後の色とが大きく異なる箇所（画素）を明示する。これによりユーザが印刷色の確認をする際に注視する部分をユーザに認識させることができる。ステップＳ４０６では、このような処理を行なう前処理として、画素ごとに、色域の変更前後の隣接画素との色差の比と閾値との比較判定処理を行なっている。
【００３０】
ステップＳ４０６で色変換部３０７は、各画素における判定結果を変数Ｐｉに格納する。変数Ｐは例えば画像データ（凹凸形状データ）の画素と同じ画素数を有し、Ｐｉはｉ番目の画素に対応するものとする。色変換部３０７は、例えば、色域の変更前後の隣接画素との色差の比が閾値より大きければＰｉに１を、小さければＰｉに０を割り当てる。なお、閾値の一例としてＬａｂ色空間における色差を用いるが、例えばＸＹＺ色空間のユークリッド距離やＲＧＢ色空間のユークリッド距離を利用しても良い。
【００３１】
次に、色変換処理装置３００は、ステップＳ４０１及びステップＳ４０２で取得した画像データ及び凹凸形状データの全ての画素についてステップＳ４０３からステップＳ４０５の処理を行ったかを判定する。処理を行っていなければ画素の番号を表す変数ｉに１を加えてステップＳ４０３へジャンプする。そうでなければ出力部３０８が色変換された画像データとステップＳ４０６の結果得られる変数Ｐとを出力し、終了に関する動作を行う。
【００３２】
＜ユーザインターフェース＞
図８はディスプレイ２０５に色変換アプリケーションの実行結果が表示されるユーザインターフェースの例を示す図である。ウインドウ８０１は色変換を行なった結果得られた画像データに基づく画像が表示される。ウインドウ８０２にも同様に、色変換を行なった結果得られる画像が表示される。なお、ウインドウ８０２の画像には、各画素位置に対してステップＳ４０６の結果得られる変数Ｐを参照し、変数Ｐｉに１が割り当てられた画素をウィンドウ８０２上に明示する。例えば図８のウィンドウ８０２においては変数Ｐｉが１が割り当てられた領域を黒塗りで示している、この黒塗り部分は、色域の変更前後で隣接画像との色差が異なる箇所（画素）である。なお、図８の例では色域の変更前後で隣接画像との色差が異なる箇所（画素）を黒塗りして表示する例を示しているが、これに限られることはない。例えば、ウィンドウを１つだけ用意し、色域の変更前後で隣接画像との色差が異なる箇所（画素）のみを点滅させてもよい。ユーザが色域変更前後で隣接画像との色差が異なる箇所（画素）を確認できるような形態であればどのような形態でもよい。例えば、閾値にΔＥ＝１（隣接レベルで色の差異がわかるレベル）を設定するとＵＩ上に閾値を越えるトーンジャンプの領域が表示される。ユーザはこのトーンジャンプの領域を確認し、画像障害の有無によって、コンボボックス８０４を操作して最大色域を使用するか最小色域を使用するか選択することができる。
【００３３】
スライドバー８０３は前述の色域の変更前後で隣接画像との色差が異なる箇所（画素）を黒塗りで表示する際の基準となる閾値をユーザが設定する際に用いられるスライドバーである。コンボボックス８０４は色変換のインテントを選択するコンボボックスであり、最小色域あるいは最大色域を利用するかを選択する。ボタン８０５は印刷実行ボタンであり、ユーザが押下すると選択されたインテントで印刷を行う。ボタン８０６はキャンセルボタンであり、ユーザが押下すると色変換アプリケーションを終了する。
【００３４】
ここで、コンボボックス８０４で選択される「最大色域」と「最小色域」について説明する。本実施例では、各画素ごとに形状に応じた縮小率を決定し、その縮小率に応じた色域データを導出する。そして、その導出した色域データを用いて各画素ごとにそれぞれ色変換が行なわる。「最大色域で印刷」とは、各画素ごとに導出された色域データをそれぞれ用いて求められた色変換後の画素値（色値）で印刷を行う処理である。「最小色域で印刷」とは、各画素の共通の色域データを用いて各画素をそれぞれ色変換した画素値（色値）で印刷を行なう処理である。各画素の共通の色域データを用いた色変換処理とは、例えば、図１（ｃ）の最小色域Ｇｂに対応する画素が１つでもあればその色域Ｇｂの色域データを用いて全ての画素を色変換処理することである。最大色域を用いる場合はプリンタの再現能力を最大限発揮して高彩度色の出力ができるようにするモードといえる。最大色域を用いると、画素ごとに異なる色域情報を使って処理を行うため、階調性が破たんする可能性がある。一方、最小色域を用いる場合は階調性を重視し、プリントの再現能力をあまり使わないモードといえる。最小色域を用いると、階調性の破たんがない代わりに高彩度の色が出力できない可能性がある。図８のコンボボックス８０４はこのようなトレードオフに対してユーザがどちらを使うかを視覚的に選択することをサポートするためのＵＩである。
【００３５】
次に、図９を用いて図８で示すＵＩの操作に応じた処理を説明する。図９は、図４に示す処理が終った後の色変換処理装置３００の処理の一例を示す図である。
【００３６】
ステップＳ９０１において色変換処理装置３００は、色変換結果の画像をディスプレイ２０５に表示する。色変換処理装置３００は、例えば図８に示すようなユーザインターフェース画面を表示する。次に、ステップＳ９０２において、色変換処理装置３００は、色域の変更指示があるかを判定する。本実施例では、デフォルトの設定が「最大色域で印刷」となっているものとする。ユーザは前述のように図８のようなＵＩを参照して、画像障害などを考慮して、色域の変更をコンボボックス８０４を介して指示する。
【００３７】
色域の変更指示があった場合、ステップＳ９０３で色変換処理装置３００は、変更後の色域でガマットマッピングを行なう。例えば、「最大色域で印刷」から「最小色域で印刷」に設定が変更された場合、色域導出部３０４は、縮小率に応じて求められた各画素の色域データを取得し、全ての画素で共通となる色域データを決定する。そして、色変換部３０７は決定された全ての画素で共通となる色域データを用いて各画素のガマットマッピングを行なう。その後、ステップＳ９０１に処理を進める。
【００３８】
一方、色域の変更指示がなかった場合、ステップＳ９０４で色変換処理装置３００は、閾値の変更指示があるかを判定する。例えばスライドバー８０３の設定が変更されたかを判定する。閾値の変更指示がある場合、ステップＳ９０５で色変換部３０７は、変更された閾値と、色域の変更前後の各画素の隣接画素との色差とを比較する。ステップＳ９０４はステップＳ４０６と同様の処理を全ての画素について行なう処理である。その後、ステップＳ９０１に処理を進める。
【００３９】
一方、閾値の変更指示がない場合、ステップＳ９０６で色変換処理装置３００は、印刷指示があるかを判定する。印刷指示がある場合、ステップＳ９０７に進み、設定された変更後の色域を用いて色変換された画素値（色値）を用いた印刷を行ない、印刷指示がない場合、ステップＳ９０８で印刷処理を中止する。
【００４０】
＜縮小率導出処理＞
次に、ステップＳ４０３で実行する、画素ｉでの色域縮小率を導出する処理を説明する。図１０は凹凸形状解析部３０３における処理の例を示すフローチャートである。
【００４１】
ステップＳ１００１で凹凸形状解析部３０３は、凹凸形状データ、及び、色域縮小率データ保持部３０５に保持されている凹凸形状と色域縮小率との関係を示す縮小率ＬＵＴを取得する。図１１は色域縮小率データ保持部３０５に保持されている縮小率ＬＵＴの一例の模式図である。色域縮小率データ保持部３０５には、図１１（ａ）に示すように凹凸形状特徴量の一例として凹凸形状データの周波数と縮小率との対応関係が記載されている。また、図１１（ｂ）に凹凸形状データの周波数と縮小率との１ＤＬＵＴ対応関係の一例を示す。凹凸形状解析部３０３は画素ｉの凹凸形状特徴量（周波数）を求め、縮小率ＬＵＴを参照して、その求められた凹凸形状特徴量（周波数）に対応する色域の縮小率を導出することになる。
【００４２】
ステップＳ１００２で凹凸形状解析部３０３は、ウェーブレット変換により、凹凸形状データの周波数分解を行う。凹凸形状データは、前述のように画素ごとに凹凸の高さ情報を有する一種の画像データである。周波数分解の方法としては、公知の離散ウェーブレット変換を用いる。ウェーブレット変換を用いることで凹凸形状データを複数の周波数帯の信号に変換することができる。なお、本実施例では説明のため離散ウェーブレット変換のレベルを２として説明する。もちろん、３以上のレベルでもよい。離散ウェーブレット変換の手法については公知の技術であるのでここでは詳細な説明は省略する。以下では、ウェーブレット変換を用いて得られた結果を用いる処理について説明する。
【００４３】
図１２は、離散ウェーブレットの結果を表す模式図である。凹凸形状データが示す形状画像Ｉは、高周波成分ＩH、低周波成分ＩLに分解され、それぞれの成分は波形の方向（縦方向、横方向、斜め方向）に応じてＩH1，ＩH2，ＩH3およびＩL1，ＩL2，ＩL3に分けられる。なお、ここでの元画像Ｉの縦横画素数Ｍ，Ｎ、画素ｉ（ｘ、ｙ）とＩL1（ｘL1、ｙL1），ＩL2（ｘL2、ｙL2），ＩL3（ｘL3、ｙL3）の関係は以下の式２のようになる。また、画素ｉ（ｘ、ｙ）とＩH1（ｘH1、ｙH1），ＩH2（ｘH2、ｙH2），ＩH3（ｘH3、ｙH3）の関係は以下式３および式４のようになる。すなわち、画素ｉに対応する位相の周波数成分は離散ウェーブレット変換の結果に基づいて、以下の式２から式４から求められる。
【００４４】
【数２】
0006512884.tif 000003
【００４５】
【数３】
0006512884.tif 000004
【００４６】
【数４】
0006512884.tif 000005
【００４７】
ステップＳ１００３で凹凸形状解析部３０３は、後述の方法により、パワーの最も大きい方向を抽出する。パワーの最も大きい方向が視覚的に最も影響すると考えられるからである。
【００４８】
ステップＳ１００４で凹凸形状解析部３０３は、ステップＳ６０３で抽出された方向の全周波数分解結果について、周波数の成分比から縮小率の線形和を算出し終了に関する動作を行う。各周波数の係数をｐkとすると、縮小率の線形和ｓは、レベル２の離散ウェーブレット変換の結果を用いているので、式５のように表せる。
【００４９】
【数５】
0006512884.tif 000006
【００５０】
すなわち、パワーの最も大きい方向の各周波数の係数に対応する縮小率を図１１に示すような縮小ＬＵＴから取得し、その得られた縮小率の線形和ｓを、画素ｉの縮小率として決定する。
【００５１】
＜方向抽出動作＞
図１３はステップＳ１００３で実行する、パワーの最も大きい方向を抽出する処理のフローチャートである。ステップＳ１３０１で凹凸形状解析部３０３は、ステップＳ１００２の周波数成分の分解結果を取得する。ステップＳ１３０２で凹凸形状解析部３０３は、画素ｉに対する各周波数成分ＩH1，ＩH2，ＩH3およびＩL1，ＩL2，ＩL3の係数を上記の式２から４に従って取得する。
【００５２】
ステップＳ１３０３で凹凸形状解析部３０３は、ＩH1とＩL1とを比較しＩH1の方が大きければＩ1maxにＩH1を設定してステップＳ１３０４へ、そうでなければＩ1maxにＩL1を設定してステップＳ１３０４へ進む。ステップＳ１３０４で凹凸形状解析部３０３は、ＩH2とＩL2とを比較しＩH2の方が大きければＩ2maxにＩH2を設定してステップＳ１３０５へ、そうでなければＩ2maxにＩL2を設定してステップＳ１３０５へ進む。ステップＳ１３０５で凹凸形状解析部３０３は、ＩH3とＩL3とを比較しＩH3の方が大きければＩ3maxにＩH3を設定してステップＳ１３０６へ、そうでなければＩ3maxにＩL3を設定してステップＳ１３０６へ進む。
【００５３】
ステップＳ１３０６で凹凸形状解析部３０３は、Ｉ1maxとＩ2maxとを比較しＩ1maxの方が大きければＳ１３０７へ、そうでなければステップＳ１３０８へ進む。ステップＳ１３０７で凹凸形状解析部３０３は、Ｉ1maxとＩ3maxとを比較しＩ1maxの方が大きければＳ１３０９へ、そうでなければステップＳ１３１０へ進む。ステップＳ１３０８で凹凸形状解析部３０３は、Ｉ2maxとＩ3maxとを比較しＩ2maxの方が大きければＳ１３１１へ、そうでなければステップＳ１３１０へ進む。
【００５４】
ステップＳ１３０９で凹凸形状解析部３０３は、画素ｉの最も大きいパワーの方向として縦方向を選択し、終了に関する動作を行う。ステップＳ１３１０で凹凸形状解析部３０３は、画素ｉの最も大きいパワーの方向として斜め方向を選択し、終了に関する動作を行う。ステップＳ１３１１で凹凸形状解析部３０３は、画素ｉの最も大きいパワーの方向として横方向を選択し、終了に関する動作を行う。
【００５５】
＜色域導出動作＞
図１４はステップＳ４０４で実行する、色域導出処理のフローチャートである。ステップＳ１４０１で色域導出部３０４は、色域特性データ保持部３０６に保持されている色域データを取得する。この色域データはプリンタ２０９を用いて平滑面上に印字した場合の色域特性を示すデータであり、例えば前述のような色変換ＬＵＴである。ステップＳ１００２で色域導出部３０４は、ステップＳ４０３で得られた縮小率を取得する。
【００５６】
ステップＳ１００３で色域導出部３０４は、色域データを縮小率に基づいて変換する。ここで、平滑面上に印字した場合の色域データ（色変換ＬＵＴ）の格子点ｊの色値をＬjａjｂj、図１（ｄ）の縮小収束点ＣＰ１の色値をＬcp1ａcp1ｂcp1、縮小率をｓとすると、色域を変更後の格子点ｊの色値Ｌj’ａj’ｂj’は式６のようになる。なお、縮小収束点とは、色域の縮小が収束する点を示しており、図１（ｄ）では彩度が０で明度が共通の点を示している。
【００５７】
【数６】
0006512884.tif 000007
【００５８】
ステップＳ１４０４で色域導出部３０４は、色変換ＬＵＴの全ての格子点について縮小率に基づく色域の変更処理を実行したかを判定し、実行していれば終了に関する動作を行い、そうでなければ処理対象の格子点ｊを変更しステップＳ１４０３に戻る。
【００５９】
本実施例によれば、画素ごとの凹凸形状に対応した色域の縮小率に応じて平滑面上に印刷層を形成した場合の色域データを変更する。そして各画素に対応する、変更した色域データを各画素の色変換を行うことにより、凹凸形状に対応した印字色を算出する。これにより、ユーザは印字色を簡易に確認できる。本実施例では、各凹凸形状に対応する色域データを直接用いないため、各凹凸形状に対応する色域データを持つ場合や多数のプロファイルを保持する場合に比べてメモリコスト削減やプロファイル作成の効率化が可能である。
【００６０】
［実施例２］
実施例１では、図１（ｂ）に示すように、凹凸形状特性を示すデータとして凹凸形状データの周波数を用いて、各周波数の係数と縮小率との対応関係を示すＬＵＴを参照して縮小率を導出して、変更後の色域を導出する例を説明した。一方、凹凸形状による色域縮小は周波数だけでなく、凹凸形状の起伏の大きさにより発生する。実施例２では、周波数及び凹凸形状の起伏の大きさ（以後、振幅と呼ぶ）と色域縮小との関係を考慮した処理について説明する。
【００６１】
図１５（ｃ）に凹凸の起伏によって、照明光が遮光される様子を模式的に示す。図１５（ｃ）では、起伏の大きさ（すなわち、周波数成分に分解した際の振幅）が図１５（ｂ）よりも大きい例を示している。図１５（ｃ）では、凹凸形状によって、入射角４５°以外の光における鏡面反射成分の一部が受光部に到達することがあることを示している。このように、図１５（ｃ）の凹凸形状では、入射光が凹凸形状に遮光されることにより表面拡散成分ならびに鏡面反射成分の一部が受光部に到達しないため、全体的な明度低下が発生し、図１５（ｄ）のＧｃに相当する色域となる。なお、図１５（ａ）及び（ｂ）と、色域ＧａとＧｂとは、図１に示すものと同様の例である。
【００６２】
上述したように、凹凸形状に応じてどのような色域縮小が発生するかは、周波数及び振幅の組み合わせに応じて分類することができる。例えば、図１６に凹凸形状の形状特性（周波数と振幅）による縮小率の分布パターンを模式的に示す。図１６のＧａ´は図１１（ａ）の略平面の凹凸形状特性における縮小特性に相当する。また、図１６のＧｂ´は図１５（ｂ）の凹凸形状特性の縮小特性（低明度部の色域縮小）を表し、図１６のＧｃ´は図１５（ｃ）の凹凸形状特性の縮小特性（光量の減少による全体明度の低下）に相当する。
【００６３】
このように、周波数とパワーの関係により、特に、形状特性Ｇｂ´〜Ｇｃ´間で色域縮小の傾向が異なる。このため、例えば画素ｉにおいて、形状特性Ｇｂ´に相当する周波数成分と形状特性Ｇｃ´に相当する周波数成分とが含まれる場合は単純な線形結合では縮小率の予測が困難となる。
【００６４】
そこで、実施例２では、周波数と振幅との関係に応じて、色域縮小の傾向が異なる周波数成分を含んでいても、縮小率を予測して色域データを導出する手法について説明する。なお、実施例１と異なる点を中心に簡潔に説明する。
【００６５】
＜縮小率算出動作＞
図１７は、実施例２においてステップＳ４０３で凹凸形状解析部３０３が実行する、画素ｉでの色域縮小率を導出する処理のフローチャートである。
【００６６】
ステップＳ１７０１で凹凸形状解析部３０３は、凹凸形状データ、及び、周波数と振幅との対応に応じた縮小率を示す縮小率ＬＵＴを取得する。ここで、図１８に各周波数の係数と振幅との対応に応じた縮小率を示す縮小率ＬＵＴの模式図を示す。図１８の模式図において、左端の列は領域を表し、左から２番目の列は振幅を表し、左から３番目の列は各周波数の係数を表し、左から４番目の列は縮小率を表す。このように、各領域にて振幅（パワー）と周波数の組み合わせについて縮小率が記されている。また、領域Ｇｂ´に相当する縮小率を表す変数をｓｂ、領域Ｇｃ´に相当する縮小率を表す変数をｓｃとし、それぞれ０に初期化する。
【００６７】
ステップＳ１７０２で凹凸形状解析部３０３は、ウェーブレット変換により、周波数分解を行う。そして、ステップＳ１７０３で凹凸形状解析部３０３は、パワーの最も大きい方向を抽出する。ステップＳ１７０２及びＳ１７０３とは、ステップＳ１００２及びステップＳ１００３と同じである。
【００６８】
ステップＳ１７０４で凹凸形状解析部３０３は、ステップＳ１７０３で抽出された方向の周波数分解結果から、ｋ（ｋ＝１（低周波成分）、２（高周波成分））番目の周波数成分について、図１６の周波数と振幅との分布を判定する。すなわち、図１６の周波数と振幅の分布のうち、領域Ｇｂ´に相当するか、領域Ｇｃ´に相当するかを判定する。領域Ｇｂ´に相当する場合はステップＳ１７０５へ、領域Ｇｃ´に相当する場合はステップＳ１７０６へ進む。
【００６９】
ステップＳ１７０５で凹凸形状解析部３０３は、領域Ｇｂ´に相当するカウントを表すＮｂに１を加えるとともに、領域Ｇｂ´に相当する縮小率の和を算出する。具体的には、ステップＳ１７０３で選択された周波数成分に対応する縮小率Ｓｋを後述の処理により算出する。そして、領域Ｇｂ´に相当する縮小率を表す変数ｓｂに、選択された周波数成分に対応する縮小率Ｓｋを加算する。すなわち、式７により縮小率の和を算出する。
【００７０】
【数７】
0006512884.tif 000008
【００７１】
ステップＳ１７０６で凹凸形状解析部３０３は、領域Ｇｃ´に相当するカウントを表すＮｃに１を加えるとともに、領域Ｇｃ´に相当する縮小率の和を算出する。具体的には、ステップＳ１７０３で選択された周波数成分に対応する縮小率ＳｋをステップＳ１７０５と同様の処理により算出する。そして、領域Ｇｃ´に相当する縮小率を表す変数ｓｃに、選択された周波数成分に対応する縮小率Ｓｋを加算する。すなわち、式８により縮小率の和を算出する。
【００７２】
【数８】
0006512884.tif 000009
【００７３】
ステップＳ１７０７で凹凸形状解析部３０３は、全ての周波数成分について分布判定を行ったかを判定し、行っていればステップＳ１７０８へ、そうでなければステップＳ１７０４へジャンプする。
【００７４】
ステップＳ１７０８で凹凸形状解析部３０３は、式９により、領域Ｇｂ´、領域Ｇｃ´にそれぞれ相当する縮小率を算出し、終了に関する動作を行う。
【００７５】
【数９】
0006512884.tif 000010
【００７６】
＜色域算出動作＞
図１９はステップＳ４０４において色域導出部３０４が実行する色域導出処理のフローチャートである。ステップＳ１９０１で色域導出部３０４は、色域データを取得する。ここで取得する色域データは、実施例１で説明したように、プリンタ２０９を用いて平滑面上に印字した場合の色域特性を示すデータである。ステップＳ１９０２で色域導出部３０４は、縮小率を取得する。ここで取得する縮小率は、ステップＳ１７０８の処理で得られた、領域Ｇｂ´、領域Ｇｃ´にそれぞれ相当する縮小率ｓｂ、ｓｃである。
【００７７】
ステップＳ１９０３で色域導出部３０４は、ステップＳ１９０１で取得した色域データとステップＳ１９０２で取得した縮小率sｂとを用いて、式１０の処理を行い色域データを変換する。式１０は実施例１で説明した式６と同様の処理である。
【００７８】
【数１０】
0006512884.tif 000011
【００７９】
ステップＳ１９０４で色域導出部３０４は、ステップＳ１９０３で変換された色域データとステップＳ１９０２で取得した縮小率ｓｃとを用いて、式１１の処理を行い色域データを変換する。ここで、ステップＳ１９０３で算出した色値Ｌj’ａj’ｂj’とし、彩度及び明度が０の点ＣＰ２の色値をＬcp2ａcp2ｂcp2とすると、変換後の色値Ｌj“ａj”ｂj“は式１１のようになる。
【００８０】
【数１１】
0006512884.tif 000012
【００８１】
ステップＳ１９０５で色域導出部３０４は、全ての格子点情報について実行したかを判定し、実行していれば終了に関する動作を行い、そうでなければステップＳ１９０３にジャンプする。
【００８２】
本実施例によれば、画像の位相に含まれる周波数と振幅の関係から色域縮小率を導出し、平滑面上に印刷層を形成した場合の色特性データとの掛け合わせによって、各凹凸形状に対応した色域を算出する。さらに、周波数と振幅により色域の縮小傾向の異なる特性を含む場合であっても、色域の縮小率を算出可能である。これにより各凹凸形状における色域情報の保持が不要、多数のプロファイルの保持する不要なため、メモリコスト削減やプロファイル作成が効率化する。
【００８３】
［実施例３］
実施例１、実施例２では、凹凸形状に対応した色域縮小率の導出について、複数の周波数成分の係数や振幅を基に色域の縮小率を導出する方法について説明した。しかしながら、実施例１、実施例２の方法では、処理に時間がかかる場合がある。以下では、単一の周波数成分の情報から色域縮小率を導出する簡易な色域導出方法について説明する。なお、実施例１、実施例２と異なる点を中心に簡潔に説明する。
【００８４】
＜縮小率導出動作＞
図２０はステップＳ４０３において凹凸形状解析部３０３が実行する画素ｉの色域の縮小率を導出する処理のフローチャートである。
【００８５】
ステップＳ２００１及びステップＳ２００２は図１０のＳ１００１及びＳ１００２と同じである。
【００８６】
ステップＳ２００３で凹凸形状解析部３０３は、ステップＳ２００２により分解した周波数成分の中からパワーの最も大きい周波数成分を抽出する。ステップＳ２００４で凹凸形状解析部３０３は、ステップＳ２００３で抽出されたパワーが最大である周波数成分に対応する縮小率を縮小率ＬＵＴから求め終了に関する動作を行う。
【００８７】
本実施例によれば、最もパワーの大きい周波数から色域縮小率を導出し、平滑面上に印刷層を形成した場合の色特性データとの掛け合わせによって、各凹凸形状に対応した色域を導出する。単一の周波数成分のみから色域縮小率を算出するため、処理時間の高速化が可能である。
【００８８】
＜他の実施例＞
なお、凹凸形状データは前述のように各画素位置における基材面からの高さに相当する情報のほか、３次元空間のＸＹＺ座標や単にＺ情報や法線情報であってもよい。
【００８９】
また、ユーザインタフェースで強調表示する画素の求め方として、色域の変更前後の隣接画素との色差の比と、閾値と、の比較を行ったが、色差を用いても良い。また、離散ウェーブレット変換のレベルを２としたが、３以上でもよい。周波数分解の方法は離散ウェーブレット変換が好適であるが、これに限定されない。高周波成分と低周波成分を分離する任意の手法が適用可能である。また、周波数と縮小率の関係はＬＵＴ形式に限定されるものではなく、γや、多項式で定義してもよい。
【００９０】
また、各実施例においては、ＬＡＢ色空間上の値で評価値を算出していたが、例えば、ＣＩＥＸＹＺ空間や、ＣＩＥＬＵＶ空間やＣＩＥＣＡＭ９７、ＣＩＥＣＡＭ０２空間等の値に変換しても良いことは言うまでもない。さらに、測色的一致を例に説明したが、色域内にマッピングする色変換であればこれに限定されるものではないことは言うまでもない。
【００９１】
また、各実施例では、測定幾何条件の一例として、４５度入射０°受光の測定系を用いた説明をしたが、ユーザの任意の観察条件に合わせた測定系を用いても良いことは言うまでもない。更に、ユーザの任意の観察条件に合わせて、図１（ｄ）に記載した色域縮小の収束点ｃｐ１ならびに、図１５（ｄ）に記載した色域縮小の収束点ｃｐ１、ｃｐ２については任意の点を設定しても構わないことは言うまでもない。
【００９２】
また、各実施例では、縮小率の線形和の算出において、ウェーブレット変換による周波数分解結果を基に最もパワーを持つ方向の情報だけを用いているが、方向判定を行わず全ての方向の全ての周波数のパワーの線形和を用いても構わない。
【００９３】
（その他の実施例）
本発明は、上述の実施形態の１以上の機能を実現するプログラムを、ネットワーク又は記憶媒体を介してシステム又は装置に供給し、そのシステム又は装置のコンピュータにおける１つ以上のプロセッサーがプログラムを読出し実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＡＳＩＣ）によっても実現可能である。
【符号の説明】
【００９４】
３０１  画像データ取得部
３０２  形状データ取得部
３０３  凹凸形状解析部
３０４  色域導出部
３０５  色域縮小率データ保持部
３０６  色域特性データ保持部
３０７  色変換部
３０８  出力部
【図１】
0006512884.tif 000013
【図２】
0006512884.tif 000014
【図３】
0006512884.tif 000015
【図４】
0006512884.tif 000016
【図５】
0006512884.tif 000017
【図６】
0006512884.tif 000018
【図７】
0006512884.tif 000019
【図８】
0006512884.tif 000020
【図９】
0006512884.tif 000021
【図１０】
0006512884.tif 000022
【図１１】
0006512884.tif 000023
【図１２】
0006512884.tif 000024
【図１３】
0006512884.tif 000025
【図１４】
0006512884.tif 000026
【図１５】
0006512884.tif 000027
【図１６】
0006512884.tif 000028
【図１７】
0006512884.tif 000029
【図１８】
0006512884.tif 000030
【図１９】
0006512884.tif 000031
【図２０】
0006512884.tif 000032