(57)【特許請求の範囲】
【請求項１】
印刷機で複数のインキにより印刷された印刷物を複数の異なる波長領域で撮像して、各波長領域の実測画像データを画素単位で出力する撮像手段と、
前記印刷物の目標状態における前記各波長領域の基準画像データを生成する画像処理手段と、
前記実測画像データと、前記基準画像データと、前記実測画像データと前記基準画像データの差分を、前記実測画像データと前記基準画像データに対応する印刷物のベタ濃度偏差、すなわち、前記両画像データに対応する印刷物のベタ濃度の差、またはその部分がベタ部でない場合には、ベタ部であると仮定して実現するベタ濃度の差、に写像する変換関数と、に基づいて、前記実測画像データと前記基準画像データの対応するベタ濃度偏差を各インキについて１画素または複数画素単位で算出する演算手段と、
前記演算手段が算出したベタ濃度偏差を集計して各インキのベタ濃度偏差を一定範囲毎に算出し、このベタ濃度偏差に基づいて、前記印刷機の各インキ量を一定範囲毎に制御する制御信号を出力する制御信号出力手段と、
を備えた色調監視装置。
【請求項２】
前記画像処理手段は、ＩＣＣプロファイルまたは多次元ＬＵＴに基づいて、上流の製版データから、前記撮像手段の実測画像データを推定した基準画像データを生成する請求項１に記載の色調監視装置。
【請求項３】
前記画像処理手段は、前記撮像手段が撮像した基準印刷物の画像に基づいて、前記基準画像データを生成する請求項１に記載の色調監視装置。
【請求項４】
前記演算手段は、前記変換関数が一次式の場合、前記実測画像データと前記基準画像データの差分ΔＣｏｌｏｒ１，ΔＣｏｌｏｒ２，・・・，ΔＣｏｌｏｒＭに対して、各インキにおける前記一次式の係数すなわち変化率係数の行列である下記のＮ×Ｍ行列
∂ink1/∂color1,∂ink1/∂Color2,・・・,∂ink1/∂ColorM
∂ink2/∂Color1,∂ink2/∂Color2,・・・,∂ink2/∂ColorM
・・・・・
∂inkN/∂Color1,∂inkN/∂Color2,・・・,∂inkN/∂ColorM
を乗算して、各インキのベタ濃度偏差Δｉｎｋ１，Δｉｎｋ２，・・・，ΔｉｎｋＮを算出して出力する請求項１乃至３のいずれかに記載の色調監視装置。
【請求項５】
前記演算手段は、前記変換関数が一次式の場合、前記実測画像データと前記基準画像データの差分ΔＣｏｌｏｒ１，ΔＣｏｌｏｒ２，・・・，ΔＣｏｌｏｒＭに対して、各インキにおける前記一次式の係数すなわち変化率係数の行列である下記のＮ×Ｍ行列
∂ink1/∂color1,∂ink1/∂Color2,・・・,∂ink1/∂ColorM
∂ink2/∂Color1,∂ink2/∂Color2,・・・,∂ink2/∂ColorM
・・・・・
∂inkN/∂Color1,∂inkN/∂Color2,・・・,∂inkN/∂ColorM
を乗算して、各インキのベタ濃度偏差Δｉｎｋ１，Δｉｎｋ２，・・・，ΔｉｎｋＮを算出し、これに基準画像に対応する印刷物のベタ濃度値、基準ｉｎｋ１、基準ｉｎｋ２、・・・、基準ｉｎｋＮを加算したものから、別途設定された目標濃度値、目標ｉｎｋ１、目標ｉｎｋ２、・・・、目標ｉｎｋＮを減算して出力する請求項１乃至３のいずれかに記載の色調監視装置。
【請求項６】
前記変換関数が一次式の場合、印刷物における各インキの網点面積率と、各インキのベタ濃度値と、に基づいて、印刷物の画像データを算出する印刷機モデルを用いて、各インキにおける前記一次式の係数すなわち変化率係数の行列を求める変化率係数取得手段を備えた請求項１乃至５のいずれかに記載の色調監視装置。
【請求項７】
前記変換関数が一次式の場合、複数の格子点を有し、画像の網点面積率を入力されると、各インキにおける前記一次式の係数すなわち変化率係数の行列の各変化率係数を出力するＬＵＴと、画像の網点面積率と、に基づいて、前記変化率係数の行列を求める変化率係数取得手段を備えた請求項１乃至５のいずれかに記載の色調監視装置。
【請求項８】
前記変換関数が一次式の場合、複数の格子点を有し、前記基準画像データまたは前記実測画像データを入力されると、各インキにおける前記一次式の係数すなわち変化率係数の行列の各変化率係数を出力するＬＵＴと、前記基準画像データまたは前記実測画像データと、に基づいて、前記変化率係数の行列を求める変化率係数取得手段を備えた請求項１乃至５のいずれかに記載の色調監視装置。
【請求項９】
前記ＬＵＴは、印刷物における各インキの網点面積率と、各インキのベタ濃度値と、に基づいて、印刷物の画像データを算出する印刷機モデルを用いて作成された請求項７または８に記載の色調監視装置。
【請求項１０】
前記ＬＵＴは、複数のインキの濃度を独立して変動させて印刷した複数のカラーチャートにおいて、各カラーチャートの各パッチにおける画像データの組に、行列の要素を、回帰分析手法を用いて求めた各パッチの変化率係数を用いて作成された請求項７または８に記載の色調監視装置。
【請求項１１】
前記撮像手段は、前記実測画像データとして、波長領域毎の反射光量または反射光量に基づく別の値を画素単位で出力する請求項１乃至１０のいずれかに記載の色調監視装置。
【請求項１２】
前記撮像手段が波長領域毎の反射光量値を出力すると、この反射光量値に基づく別の値に変換する変換手段を備えた請求項１乃至１１のいずれかに記載の色調監視装置。
【請求項１３】
前記印刷物に複数の監視点を設定して、その座標情報を前記画像処理手段及び前記演算手段に出力する監視点設定手段を備え、
前記画像処理手段及び前記演算手段は、監視点単位で処理を行う請求項１乃至１２のいずれかに記載の色調監視装置。
【請求項１４】
前記演算手段が算出した各インキの１画素または複数画素単位のベタ濃度偏差、または前記制御信号出力手段が算出した各インキの一定範囲毎のベタ濃度偏差の少なくとも一方を報知する報知手段を備えた請求項１乃至１３のいずれかに記載の色調監視装置。
【請求項１５】
前記制御信号出力手段に代えて、前記演算手段が出力した各インキの１画素または複数画素単位のベタ濃度偏差を報知する報知手段を備えた請求項１乃至１３のいずれかに記載の色調監視装置。
【請求項１６】
前記制御信号出力手段に代えて、
前記演算手段が算出した各インキの１画素または複数画素単位のベタ濃度偏差を集計して各インキのベタ濃度偏差を一定範囲毎に算出する集計手段と、
前記演算手段が演算した各インキの１画素または複数画素単位のベタ濃度偏差と、前記集計手段が算出した各インキの一定範囲毎のベタ濃度偏差と、の少なくとも一方を報知する報知手段と、
を備えた請求項１乃至１３のいずれかに記載の色調監視装置。
【請求項１７】
請求項１乃至１４のいずれかに記載の色調監視装置と、
前記色調監視装置の制御信号出力手段が出力した制御信号に基づいて、前記印刷機が印刷に用いるインキ量を調整するインキ量制御装置と、
を備えた色調管理システム。
【請求項１８】
印刷物の色調を監視する色調監視装置の色調監視方法であって、
撮像手段が、印刷機で複数のインキにより印刷された印刷物を複数の異なる波長領域で撮像して、各波長領域の実測画像データを画素単位で出力する手順と、
画像処理手段が、前記印刷物の目標状態における前記各波長領域の基準画像データを生成する手順と、
演算手段が、前記実測画像データと、前記基準画像データと、前記実測画像データと前記基準画像データの差分を前記実測画像データと前記基準画像データのベタ濃度偏差、すなわち、前記両画像データに対応する印刷物のベタ濃度の差、またはその部分がベタ部でない場合には、ベタ部であると仮定して実現するベタ濃度の差、に写像する変換関数と、に基づいて前記実測画像データと前記基準画像データの対応するベタ濃度偏差を各インキについて１画素または複数画素単位で算出する手順と、
報知手段が、前記演算手段が算出したベタ濃度偏差を報知する手順と、
を備えた色調監視方法。
【請求項１９】
印刷物の色調を監視する色調監視装置と、印刷機のインキ量を調整するインキ量制御装置を含む色調管理システムの色調管理方法であって、
撮像手段が、印刷機で複数のインキにより印刷された印刷物を複数の異なる波長領域で撮像して、各波長領域の実測画像データを画素単位で出力する手順と、
画像処理手段が、前記印刷物の目標状態における前記各波長領域の基準画像データを生成する手順と、
演算手段が、前記実測画像データと、前記基準画像データと、前記実測画像データと前記基準画像データの差分を前記実測画像データと前記基準画像データのベタ濃度偏差、すなわち、前記両画像データに対応する印刷物のベタ濃度の差、またはその部分がベタ部でない場合には、ベタ部であると仮定して実現するベタ濃度の差、に写像する変換関数と、に基づいて前記実測画像データと前記基準画像データの台頭するベタ濃度偏差を各インキについて１画素または複数画素単位で算出する手順と、
制御信号出力手段が、前記演算手段が算出したベタ濃度偏差を集計して一定範囲毎のベタ濃度偏差を算出し、このベタ濃度偏差に基づいて、前記印刷機の各インキ量を一定範囲毎に制御する制御信号を出力する手順と、
インキ量制御手段が、前記制御信号出力手段が出力した制御信号に基づいて、前記印刷機が印刷に用いるインキ量を一定範囲毎に調整する手順と、
を備えた色調管理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、印刷機で連続的に印刷した印刷物の色調の変動を監視する色調監視装置及び色調監視方法、並びに色調の変動の監視結果に基づき、インキ量を補正して一定の色調に保つ色調管理システム及び色調管理方法に関する。
【背景技術】
【０００２】
印刷機で印刷した印刷面の色調を目的の色調に合わせるための手段として、従来はカラーパッチ等の単色のベタ部や特定の組み合わせのパターンを印刷物の絵柄以外の部分に印刷し、その部分の濃度や分光特性または色座標値などの表色系で管理する方法が多く提案されている。
【特許文献１】特開２００１−１８３６４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００３】
しかし、カラーパッチを用いた場合、その部分の色を監視することは容易であるが、実際の印刷物は場所によってインキ毎の消費量が異なり、同じ網点面積率の組み合わせであっても同じようには印刷されないので、本当に合わせたい部分の色を合わせることにはならないという問題があった。
【０００４】
また、これらは選択部分を濃度や分光特性または色座標値などの表色系で管理する方法であるため、監視している部分の現在の色の状態はわかっても、そのためにはどの場所のどのインキをどれくらい増減させれば良いかが明確ではないという問題があった。
【０００５】
さらに、新聞やチラシ印刷のように、断裁しろが無く、カラーパッチを設けることができない印刷物があり、この場合には上記の方法が使用できないという問題があった。
【０００６】
そこで、本発明は、色調の監視及び色合わせを容易に行うことができ、また画像の画素単位でインキ濃度の変化を検出でき、さらにカラーパッチ等を用いなくてもインキ濃度を正確に推定できる色調監視装置、色調管理システム、色調監視方法、及び色調管理方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
（１）印刷機で複数のインキにより印刷された印刷物を複数の異なる波長領域で撮像して、各波長領域の実測画像データを画素単位で出力する撮像手段と、
前記印刷物の目標状態における前記各波長領域の基準画像データを生成する画像処理手段と、
前記実測画像データと、前記基準画像データと、前記実測画像データと前記基準画像データの差分を、前記実測画像データと前記基準画像データに対応する印刷物のベタ濃度偏差、すなわち、前記両画像データに対応する印刷物のベタ濃度の差、またはその部分がベタ部でない場合には、ベタ部であると仮定して実現するベタ濃度の差、に写像する変換関数と、に基づいて、前記実測画像データと前記基準画像データの対応するベタ濃度偏差を各インキについて１画素または複数画素単位で算出する演算手段と、
前記演算手段が算出したベタ濃度偏差を集計して各インキのベタ濃度偏差を一定範囲毎に算出し、このベタ濃度偏差に基づいて、前記印刷機の各インキ量を一定範囲毎に制御する制御信号を出力する制御信号出力手段と、
を備えたことを特徴とする。
【０００８】
この構成においては、色調監視装置は、撮像手段が印刷物の画像を撮像して出力した実測画像データと、画像処理手段が生成した基準画像データと、
変換関数と、に基づいて、前記実測画像データと前記基準画像データの対応するベタ濃度偏差を各インキについて１画素または複数画素単位で算出する。そして、算出したベタ濃度偏差を集計して各インキのベタ濃度偏差を一定範囲毎に算出し、このベタ濃度偏差に基づいて、前記印刷機の各インキ量を一定範囲毎に制御する制御信号を出力する。したがって、特定のパターンではなく、印刷された面全体を確認できるため、実際に色調を合わせたい部分を見ることが可能となるとともに、実測画像データと基準画像データの差は常に各インキの単色ベタ濃度の変動に結びついているため、どの部分のどのインキをどれくらい変動させれば目標とする画像データのベタ濃度に一致させられるかが明確となる。また、カラーバッチのようなベタ部または色毎に独立した既知の網点面積率の部分が無くても、ベタ部の濃度値を推定することができる。さらに、どの部分のどのインキをどれくらい変動させれば良いかを、インキ供給機構の関係する領域毎に集計し、インキ供給量を制御することにより自動的に色調を目的とするものに合わせることが容易となる。加えて、印刷機が例えばオフセット輪転機の場合には、インキキーの幅単位で各インキのベタ濃度偏差を集計し、その結果に基づく制御信号をインキキー制御装置に出力することで、印刷機のインキ量をフィードバック制御することができ、印刷物の色調を常に一定に管理することができる。
【０００９】
（２）前記画像処理手段は、ＩＣＣプロファイルまたは多次元ＬＵＴに基づいて、上流の製版データから、前記撮像手段の実測画像データを推定した基準画像データを生成することを特徴とする。
【００１０】
この構成においては、画像処理手段は、基準画像データとして、ＩＣＣプロファイルまたは多次元ＬＵＴに基づいて、上流の製版データから生成された基準画像データを生成する。したがって、上流の製版データから容易に基準画像データを取得できる。また、上流の製版データから基準画像データを生成・取得することにより、基準画像をセンサから手動で取り込む必要がなくなり、印刷開始直後から全自動で色調の監視・制御が可能となり、損紙の低減・品質の向上・印刷時間の短縮等を実現できる。さらに、上流の製版データから基準画像データを生成・取得することにより、色調を合わせる基準が統一され、事業所単位、会社単位での一貫したカラーマネジメントが可能になる。
【００１１】
（３）前記画像処理手段は、前記撮像手段が撮像した基準印刷物の画像に基づいて、前記基準画像データを生成することを特徴とする。
【００１２】
この構成においては、画像処理手段は、撮像手段で撮像した基準印刷物の画像から基準画像データを生成する。したがって、上流の製版データが取得できない場合でも、基準画像データを容易に取得できる。また、印刷物の色調を人が望ましい色調に調整を行った後に基準画像データを取得し、それを基準に色調の監視・制御を行うことにより、人が判断した望ましい色調を維持することができる。
【００１３】
（４）前記演算手段は、前記変換関数が一次式の場合、前記実測画像データと前記基準画像データの差分ΔＣｏｌｏｒ１，ΔＣｏｌｏｒ２，・・・，ΔＣｏｌｏｒＭに対して、各インキにおける前記一次式の係数すなわち変化率係数の行列である下記のＮ×Ｍ行列
∂ink1/∂color1,∂ink1/∂Color2,・・・,∂ink1/∂ColorM
∂ink2/∂Color1,∂ink2/∂Color2,・・・,∂ink2/∂ColorM
・・・・・
∂inkN/∂Color1,∂inkN/∂Color2,・・・,∂inkN/∂ColorM
を乗算して、各インキのベタ濃度偏差Δｉｎｋ１，Δｉｎｋ２，・・・，ΔｉｎｋＮを算出して出力することを特徴とする。
【００１４】
この構成においては、演算手段は、変換関数が一次式の場合、各インキにおける前記一次式の係数すなわち変化率係数の行列であるＮ×Ｍ行列を、実測画像データと基準画像データの差分ΔＣｏｌｏｒ１，ΔＣｏｌｏｒ２，・・・，ΔＣｏｌｏｒＭに乗算して、各インキのベタ濃度偏差Δｉｎｋ１，Δｉｎｋ２，・・・，ΔｉｎｋＮを算出する。したがって、各インキのベタ濃度偏差を容易に算出できる。また、インキのベタ濃度偏差算出に必要な演算量が少ないため、高速に処理できる。さらに、係数はあらかじめ求めておけばよいので、係数を求める際には、必要な精度に応じて複雑で時間のかかる処理を採用することができる。
【００１５】
（５）前記演算手段は、前記変換関数が一次式の場合、前記実測画像データと前記基準画像データの差分ΔＣｏｌｏｒ１，ΔＣｏｌｏｒ２，・・・，ΔＣｏｌｏｒＭに対して、各インキにおける前記一次式の係数すなわち変化率係数の行列である下記のＮ×Ｍ行列
∂ink1/∂color1,∂ink1/∂Color2,・・・,∂ink1/∂ColorM
∂ink2/∂Color1,∂ink2/∂Color2,・・・,∂ink2/∂ColorM
・・・・・
∂inkN/∂Color1,∂inkN/∂Color2,・・・,∂inkN/∂ColorM
を乗算して、各インキのベタ濃度偏差Δｉｎｋ１，Δｉｎｋ２，・・・，ΔｉｎｋＮを算出し、これに基準画像に対応する印刷物のベタ濃度値、基準ｉｎｋ１、基準ｉｎｋ２、・・・、基準ｉｎｋＮを加算したものから、別途設定された目標濃度値、目標ｉｎｋ１、目標ｉｎｋ２、・・・、目標ｉｎｋＮを減算して出力することを特徴とする。
【００１６】
この構成においては、演算手段は、変換関数が一次式の場合、各インキにおける前記一次式の係数すなわち変化率係数の行列であるＮ×Ｍ行列を、実測画像データと基準画像データの差分ΔＣｏｌｏｒ１，ΔＣｏｌｏｒ２，・・・，ΔＣｏｌｏｒＭに乗算して、各インキのベタ濃度偏差Δｉｎｋ１，Δｉｎｋ２，・・・，ΔｉｎｋＮを算出する。そして、これに基準画像に対応する印刷物のベタ濃度値、基準ｉｎｋ１、基準ｉｎｋ２、・・・、基準ｉｎｋＮを加算したものから、別途設定された目標濃度値、目標ｉｎｋ１、目標ｉｎｋ２、・・・、目標ｉｎｋＮを減算して出力する。したがって、監視画像と基準画像の濃度差ΔＩｎｋに基準画像の濃度を加算することにより、監視画像の濃度を求めることができ、監視画像の濃度と任意に指定された目標濃度との差分を新たにΔｉｎｋとして出力することにより、基準画像の濃度とは異なる目標濃度に合わせた印刷が可能となる。
【００１７】
（６）前記変換関数が一次式の場合、印刷物における各インキの網点面積率と、各インキのベタ濃度値と、に基づいて、印刷物の画像データを算出する印刷機モデルを用いて、各インキにおける前記一次式の係数すなわち変化率係数の行列を求める変化率係数取得手段を備えたことを特徴とする。
【００１８】
この構成においては、変換関数が一次式の場合、印刷機モデルを用いて、各インキにおける一次式の係数すなわち変化率係数の行列を求める変化率係数作成手段を備えている。したがって、印刷画像における各インキの網点面積率と、各インキのベタ濃度値と、に基づいて、変化率係数を作成できる。
【００１９】
（７）前記変換関数が一次式の場合、複数の格子点を有し、画像の網点面積率を入力されると、各インキにおける前記一次式の係数すなわち変化率係数の行列の各変化率係数を出力するＬＵＴと、画像の網点面積率と、に基づいて、前記変化率係数の行列を求める変化率係数取得手段を備えたことを特徴とする。
【００２０】
この構成においては、変換関数が一次式の場合、複数の格子点を有し、画像の網点面積率を入力されると、各インキにおける一次式の係数すなわち変化率係数の行列の各変化率係数を出力するＬＵＴに対して、画像の網点面積率を入力して、変化率係数の行列の各変化率係数を得ることができる。したがって、予めＬＵＴを作成しておくことで、画像の網点面積率を入力することで、短時間で容易に各変化率係数を得られる。また、ＬＵＴを用いることにより変化率係数を求める演算量が減少し、高速に処理できる。さらに、ＬＵＴはあらかじめ作成しておけばよいので、ＬＵＴを作成する際には、必要な精度に応じて複雑で時間のかかる処理を採用することができる。
【００２１】
（８）前記変換関数が一次式の場合、複数の格子点を有し、前記基準画像データまたは前記実測画像データを入力されると、各インキにおける一次式の係数すなわち変化率係数の行列の各変化率係数を出力するＬＵＴと、基準画像データまたは実測画像データと、に基づいて、前記変化率係数の行列を求める変化率係数取得手段を備えたことを特徴とする。
【００２２】
この構成においては、変換関数が一次式の場合、各インキにおける前記一次式の係数すなわち変化率係数の行列を、複数の格子点を有し、実測画像データを入力されると変化率係数の行列の各変化率係数を出力するＬＵＴと、実測画像データと、に基づいて求める。ＬＵＴの入力は、実測画像データである撮像手段が出力した画像（センサ画像）でもよいので、上流からインキのベタ濃度値画像および基準画像をもらわないシステムに用いることができる。
【００２３】
（９）前記ＬＵＴは、印刷物における各インキの網点面積率と、各インキのベタ濃度値と、に基づいて、印刷物の画像データを算出する印刷機モデルを用いて作成されたことを特徴とする。
【００２４】
この構成においては、印刷機モデルを用いてＬＵＴが作成される。したがって、印刷画像における各インキの網点面積率と、各インキのベタ濃度値と、を与えることで、ＬＵＴを作成することができる。また、印刷機モデルからＬＵＴを作成するため、回帰分析を実施しなくても良い。
【００２５】
（１０）前記ＬＵＴは、複数のインキの濃度を独立して変動させて印刷した複数のカラーチャートにおいて、各カラーチャートの各パッチにおける画像データの組に、行列の要素を、回帰分析手法を用いて求めた各パッチの変化率係数を用いて作成されたことを特徴とする。
【００２６】
この構成においては、複数のインキの濃度を独立して変動させて印刷した複数のカラーチャートの各格子点における画像データに対して、回帰分析手法を用いてＬＵＴを作成する。したがって、精度の良いＬＵＴを作成することが可能となる。また、条件がよければ印刷機モデルを用いる場合よりもはるかに高い精度に到達可能である。さらに、印刷機モデルを構築する必要がないため、いかなる印刷方式に対しても適用が可能である。
【００２７】
（１１）前記撮像手段は、前記実測画像データとして、波長領域毎の反射光量または反射光量に基づく別の値を画素単位で出力することを特徴とする。
【００２８】
この構成においては、撮像手段は、波長領域毎の反射光量値または反射光量値に基づく別の値を実測画像データとして画素単位で出力する。したがって、例えば、反射光量値の変換を対数関数のようなものにすることにより、濃度の濃い部分と薄い部分で要求する精度を変えることができ、濃度の階調に近い、適正な精度を全体に渡って得ることができる。
【００２９】
（１２）前記撮像手段が波長領域毎の反射光量値を出力すると、この反射光量値に基づく別の値に変換する変換手段を備えたことを特徴とする。
【００３０】
この構成においては、変換手段は、撮像手段が出力した波長領域毎の反射光量値を別の値に変換する。したがって、撮像手段が波長領域毎の反射光量値を出力する場合には、別の値に変換することで、（１１）と同様の効果が得られる。
【００３１】
（１３）前記印刷物に複数の監視点を設定して、その座標情報を前記画像処理手段及び前記演算手段に出力する監視点設定手段を備え、
前記画像処理手段及び前記演算手段は、監視点単位で処理を行うことを特徴とする。
【００３２】
この構成においては、監視点設定手段により設定された監視点単位で、画像処理手段及び演算手段は処理を行う。したがって、適切に監視点を設定することで監視する画素数を削減できるので、色調の監視処理を高速化できる。
【００３３】
（１４）前記演算手段が算出した各インキの１画素または複数画素単位のベタ濃度偏差、または前記制御信号出力手段が算出した各インキの一定範囲毎のベタ濃度偏差の少なくとも一方を報知する報知手段を備えたことを特徴とする。
【００３４】
この構成においては、報知手段は、演算手段が算出した各インキのベタ濃度偏差を報知する。したがって、印刷機のオペレータは、報知手段の報知内容に基づいてインキ量を調整することができる。
【００３５】
（１５）前記制御信号出力手段に代えて、前記演算手段が出力した各インキの１画素または複数画素単位のベタ濃度偏差を報知する報知手段を備えたことを特徴とする。
【００３６】
この構成においては、演算手段が出力した各インキのベタ濃度偏差を報知手段で報知する。したがって、印刷された面全体を比較できるため、実際に色調を合わせたい部分を見ることが可能となるとともに、推定画像データの差は常に各インキの単色ベタ濃度の変動に結びついているため、どの部分のどのインキをどれくらい変動させれば目標とする画像データのベタ濃度に一致させられるかが明確となり、オペレータは、報知手段の報知内容に基づいて印刷機のインキ量を制御可能となる。また、カラーバッチが無くても、ベタ部の濃度値を推定することができる。
【００３７】
（１６）前記制御信号出力手段に代えて、
前記演算手段が算出した各インキの１画素または複数画素単位のベタ濃度偏差を集計して各インキのベタ濃度偏差を一定範囲毎に算出する集計手段と、
前記演算手段が演算した各インキの１画素または複数画素単位のベタ濃度偏差と、前記集計手段が算出した各インキの一定範囲毎のベタ濃度偏差と、の少なくとも一方を報知する報知手段と、
を備えたことを特徴とする。
【００３８】
この構成においては、報知手段は、演算手段が演算した各インキの１画素または複数画素単位のベタ濃度偏差と、集計手段が算出した各インキの一定範囲毎のベタ濃度偏差と、の少なくとも一方を報知する。したがって、印刷機のオペレータは、報知手段の報知内容に基づいてインキ量を調整することができる。
【００３９】
（１７）（１）乃至（１４）のいずれかに記載の色調監視装置と、
前記色調監視装置の制御信号出力手段が出力した制御信号に基づいて、前記印刷機が印刷に用いるインキ量を調整するインキ量制御装置と、
を備えたことを特徴とする。
【００４０】
この構成においては、色調監視装置は、印刷機が印刷した印刷物の色調を監視し、制御信号出力手段で、印刷機のインキ量を調整する。したがって、印刷物の色調を常に一定範囲内に精度良く保つことができる。
【００４１】
（１８）印刷物の色調を監視する色調監視装置の色調監視方法であって、
撮像手段が、印刷機で複数のインキにより印刷された印刷物を複数の異なる波長領域で撮像して、各波長領域の実測画像データを画素単位で出力する手順と、
画像処理手段が、前記印刷物の目標状態における前記各波長領域の基準画像データを生成する手順と、
演算手段が、前記実測画像データと、前記基準画像データと、前記実測画像データと前記基準画像データの差分を前記実測画像データと前記基準画像データのベタ濃度偏差、すなわち、前記両画像データに対応する印刷物のベタ濃度の差、またはその部分がベタ部でない場合には、ベタ部であると仮定して実現するベタ濃度の差、に写像する変換関数と、に基づいて前記実測画像データと前記基準画像データの対応するベタ濃度偏差を各インキについて１画素または複数画素単位で算出する手順と、
報知手段が、前記演算手段が算出したベタ濃度偏差を報知する手順と、
を備えたことを特徴とする。
【００４２】
この構成においては、（１５）と同様の効果を得ることができる。
【００４３】
（１９）  印刷物の色調を監視する色調監視装置と、印刷機のインキ量を調整するインキ量制御装置を含む色調管理システムの色調管理方法であって、
撮像手段が、印刷機で複数のインキにより印刷された印刷物を複数の異なる波長領域で撮像して、各波長領域の実測画像データを画素単位で出力する手順と、
画像処理手段が、前記印刷物の目標状態における前記各波長領域の基準画像データを生成する手順と、
演算手段が、前記実測画像データと、前記基準画像データと、前記実測画像データと前記基準画像データの差分を前記実測画像データと前記基準画像データのベタ濃度偏差、すなわち、前記両画像データに対応する印刷物のベタ濃度の差、またはその部分がベタ部でない場合には、ベタ部であると仮定して実現するベタ濃度の差、に写像する変換関数と、に基づいて前記実測画像データと前記基準画像データの台頭するベタ濃度偏差を各インキについて１画素または複数画素単位で算出する手順と、
制御信号出力手段が、前記演算手段が算出したベタ濃度偏差を集計して一定範囲毎のベタ濃度偏差を算出し、このベタ濃度偏差に基づいて、前記印刷機の各インキ量を一定範囲毎に制御する制御信号を出力する手順と、
インキ量制御手段が、前記制御信号出力手段が出力した制御信号に基づいて、前記印刷機が印刷に用いるインキ量を一定範囲毎に調整する手順と、
を備えたことを特徴とする。
【００４４】
この構成においては、（１）と同様の効果を得ることができる。
【発明の効果】
【００４５】
本発明によれば、印刷物を取り込んだ光学センサの実測画像データを、センサの画像と等価な基準画像データと比較し、その差から各インキの単色ベタ部の濃度変動と結びつける方法を採用する。これにより、カラーパッチなどの特定のパターンではなく、印刷された面全体を比較できるため、実際に色調を合わせたい部分を見ることが可能となるとともに、実測画像データと基準画像データの差は常に各インキの単色ベタ濃度の変動に結びついているため、どの部分のどのインキをどれくらい変動させれば目標とする基準画像データに一致させられるかが明確となる。また、どの部分のどのインキをどれくらい変動させれば良いかを、インキ供給機構の関係する領域毎に集計し、インキ供給量を制御することにより自動的に色調を目的とするものに合わせることが容易となる。また、新聞、ちらし等のように断裁しろが無く、カラーパッチを印刷することができない印刷においても色調を合わせることができる。
【発明を実施するための最良の形態】
【００４６】
図１は、本発明の実施形態に係る色調管理システムを含む印刷システムの概略構成を示したブロック図である。図１に示す色調管理システム３は、印刷物の一例として新聞紙の色調を管理するシステムである。なお、本構成の新聞社側の構成については一般的なものを記述しているが、これに限るものではない。
【００４７】
印刷システム１は、モニタ１１３が接続された紙面編集装置（コンピュータ）１１、ＲＩＰ(Raster Image Processor)１２、画像ファイルを蓄積・提供する画像サーバ１３、校正刷りデータの画像処理を行う画像変換装置（コンピュータ）１４、校正紙を印刷する校正用プリンタ１５、刷版を製作するＣＴＰ（Computer to Plate)１７、印刷物（新聞紙）を印刷するオフセット輪転機（以下、単に輪転機と称する。）１９、印刷物の色調を監視する色調監視装置２１、色調監視装置２１が推定したインキの状態を表示するモニタ２２、インキキー７２Ａ〜７２Ｈの開度を制御するインキキー制御装置２３、及び輪転機１９で印刷された印刷物の画像を撮像する光学センサ２５Ｆ，２５Ｒを備えた構成である。また、本発明の色調管理システム３は、色調監視装置２１と、インキキー制御装置２３と、光学センサ２５Ｆ・２５Ｒから成る。なお、輪転機１９のインキの状態を報知するために、モニタ２２により画像により色調の情報を報知するだけでなくスピーカやブサー等を設けて音声により色調の情報を報知するようにしても良い。
【００４８】
ここで、輪転機１９では、印刷物（新聞紙）の両面に文字や絵柄が印刷され、これらを光学センサ２５Ｆ・２５Ｒで撮像して、印刷物の両面の色調を管理するが、以下の説明では、説明を簡略化するために印刷物の表面を光学センサ２５Ｆで撮像して色調監視・管理を行う説明のみ記載するが、印刷物の裏面についても同様に光学センサ２５Ｒで撮像して色調監視・管理を行うものとする。
【００４９】
紙面編集装置１１は、編集者Ｈが、記事やカラー写真、クライアントＫから依頼された新聞掲載用のカラー広告等をモニタ１１３に表示させて紙面編集を行うためのものである。紙面編集装置１１は、標準的な輪転機に合った色調特性となるように設定されたプロファイルを記憶部１１１で記憶している。紙面編集装置１１では、制御部１１２が記憶部１１１からこのプロファイルを読み出して、編集されたカラー広告等の画像をプロファイル変換してＲＩＰ１２へ出力する。ＲＩＰ１２は、プロファイル変換された画像から複数のインキに対応する２値網点画像データを作成し、画像サーバ１３にこれらの２値網点画像データを出力する。ＲＩＰ１２は、２値網点画像データとして、例えば１ｂｉｔ  Ｔｉｆｆ形式の画像データを作成する。
【００５０】
画像サーバ１３は、これらの画像データを蓄積する。
【００５１】
画像変換装置１４は、校正用プリンタ１５の色調特性となるように設定されたプロファイルを記憶部１４１で記憶している。画像変換装置１４では、制御部１４２が記憶部１４１から、このプロファイルを読み出すとともに、画像サーバ１３に蓄積された各２値網点画像データを読み出して、校正用プリンタ１５の色調特性となるように多値化及びプロファイル変換を行い、変換した画像データを校正用プリンタ１５に出力する。
【００５２】
校正用プリンタ１５は、画像変換装置１４から出力された画像データに基づく画像を印刷して校正紙３１を出力する。クライアントＫは、校正紙３１の内容を確認して、校正紙３１の内容に問題があれば編集者Ｈに問題点を伝える。この場合、編集者Ｈは、紙面編集装置１１で画像を修正し、校正用プリンタ１５で印刷した校正紙３１を再度クライアントＫに確認してもらう。クライアントＫは、校正紙３１の内容に問題が無ければ（または無くなれば）、その旨を編集者Ｈに伝える。
【００５３】
紙面編集装置１１は、編集者Ｈにより校正紙３１の内容に問題が無い旨の操作を検出すると、この画像データを正規の画像データとしてＲＩＰ１２に出力し、画像サーバ１３に２値網点画像データを蓄積させる。
【００５４】
ＣＴＰ１７は、画像サーバ１３から正規の画像データ（２値網点画像データ）を読み出して、輪転機１９の各インキ用の刷版を制作する。
【００５５】
輪転機１９は、ＣＴＰ１７で製作された刷版がセットされており、ロール紙Ｒからウェブ状の連続紙Ｐを供給し、印刷ユニット７３でブラック(Black：以下、Ｋと称する。）・シアン（Cyan：以下、Ｃと称する。）・マゼンタ(Magenta：以下、Ｍと称する。）・イエロー（Yellow：以下、Ｙと称する。）のインキが複数のインキキー７２Ａ〜７２Ｈから供給され、連続紙Ｐの両面に文字や絵柄を印刷する。なお、インキの供給順はこの順番に限定されるものではなく、新聞社や印刷工場により異なる順番でインキが供給される。インキキー７２Ａ〜７２Ｈは、一定幅であり、連続紙Ｐの紙面に対して複数配列されている。また、輪転機１９では、図示していないが、例えば、ＣＭＹＫとオレンジとグリーンの６色によるヘキサクローム印刷や、ＣＭＹＫと赤・緑・青（ＲＧＢ）の７色によるハイファイ印刷など５色以上のインキによる印刷が可能である。また、輪転機１９では、ＣＭＹＫと１つまたは複数の特色（スポットカラー、特練色とも言う。）インキによる特色印刷も可能である。
【００５６】
印刷ユニット７３の後段（下流側）には、連続紙Ｐの表側に対向して光学センサ２５Ｆが、連続紙Ｐの裏面側に対向して光学センサ２５Ｒが設けられている。光学センサ２５Ｆ・２５Ｒは、印刷ユニット７３でＣＭＹＫや他の色のインキによる印刷が完了した走行中の印刷紙面の表面（表胴側）及び裏面（裏胴側）を撮像して、複数の波長領域の画像データ（センサ値）を画素毎に出力する。
【００５７】
光学センサ２５Ｆ・２５Ｒは、一例として、不図示の光源から照射されて連続紙Ｐの表面で反射した光を受光して、赤（Red:以下、Ｒと称する。）・緑（Green:以下、Ｇと称する。）・青（Blue:以下、Ｂと称する。）・赤外（Infrared rays:以下、Ｉｒと称する。）の４つの波長領域で撮像を行い、Ｒ・Ｇ・Ｂ・Ｉｒ（以下、ＲＧＢＩｒと記載する場合もある。）の４つの波長領域の画像データ（センサ値）を画素毎に出力する。なお、光学センサ２５Ｆ・２５Ｒは、一例として反射光量値をデジタルデータとして出力する。
【００５８】
なお、不図示の光源は、可視光の全波長を含む光を照射するものでも、Ｒ・Ｇ・Ｂ・Ｉｒの４つの波長領域のみの光を照射するものでも良く、また、１つの光源でも複数の光源でも良い。なお、ＣＭＹＫに加えて他の色のインキを使用して印刷を行う場合には、印刷した絵柄をＲＧＢＩｒの波長領域で撮像しても、さらに他の波長領域で撮像しても良い。
【００５９】
輪転機１９には、複数のインキキー７２Ａ〜７２Ｈの開度を制御して紙面に供給するインキ量を調整するインキキー制御装置２３が接続されて、色調監視装置２１には、モニタ２２が接続されている。インキキー制御装置２３は、色調監視装置２１からの制御信号により、またはオペレータがモニタ２２の表示内容に応じて行う操作に応じて、各インキキー７２Ａ〜７２Ｈの開度を調整する。
【００６０】
色調管理システム３の色調監視装置２１は、輪転機１９が印刷したドライダウン前の新聞紙画像を光学センサ２５Ｆで監視し、インキキー制御装置２３に制御信号を出力して、輪転機１９の各インキキー７２Ａ〜７２Ｈの開度を調整する。これにより、輪転機１９で印刷した連続紙Ｐをドライダウンした後の新聞紙３３の色調（インキ濃度値）が、新聞社で決めた適正印刷状態の色調に、結果的に、ほぼ一致するようになる。
【００６１】
次に、本願発明の色調監視の考え方について説明する。
【００６２】
ここで、画像データの変化分と各インキの変化分との関係を表す変換関数により解を求める方法としては、次のようなものがある。
【００６３】
変換関数は、ΔR,ΔG,ΔB,ΔIrをΔC,ΔM,ΔY,ΔKに写像するもので、一般にΔR,ΔG,ΔB,ΔIrに関する多項式等が当てはまる。本実施例ではΔRGBIｒの1次式で説明しているが、例えば、変換関数が2次式の場合には、下記の式のようになる（但し、ΔCのみ記載)。
【００６４】
ΔC＝∂C/∂R×ΔR＋∂C/∂G×ΔG＋∂C/∂B×ΔB＋∂C/∂Ir×ΔIr
＋1/2{∂2C/∂2R×ΔR2＋∂2C/∂2G×ΔG2＋∂2C/∂2B×ΔB2＋∂2C/∂2Ir×ΔIr2
＋2×(∂2C/∂R∂G×ΔR×ΔG＋∂2C/∂R∂B×ΔR×ΔB＋∂2C/∂R∂Ir×ΔR×ΔIr
＋∂2C/∂G∂B×ΔG×ΔB＋∂2C/∂G∂Ir×ΔG×ΔIr＋∂2C/∂B∂Ir×ΔB×ΔIr)}
また、変換関数が３次式の場合には下記のようになる(但し、ΔCのみ記載)。
【００６５】
【数１】
0004275716.tif 000002
【００６６】
また、変換関数としては、さらに高次の式を用いることも可能である。
【００６７】
本発明では、上記の多項式による近似や、その他の関数による近似により変換関数の解を求める。
【００６８】
なお、以下の説明では、説明を簡略化するために、変換関数が一次式の場合を例に挙げて説明する。また、この一次式の係数を変化率係数と称する。
【００６９】
ここで、用語及び記号を以下のように定義する。
センサ値：光学センサが検出した反射光量そのもの、または、光学センサが検出した反射光量と白紙部の反射光量との関係を以下のように演算した値とする。
【００７０】
センサ値＝-log10(光量／白紙部の光量)
但し、logの底は10に限るものではなく、他の値でも良い。なお、センサ値は、これらの値に限定するものではなく、その他、任意の１価関数でも良く、光学センサ２５Ｆ・２５Ｒが反射光量値を出力するのか、それとも反射光量値を変換した所定の値（反射光量値とは別の値）を出力するのか、に応じて定義すればよい。
【００７１】
C,M,Y,K         ：  印刷物の単色ベタ部の濃度
R,G,B,Ir        ：  印刷物の監視点のセンサ値
ΔC,ΔM,ΔY,ΔK ：  基準印刷物に対する監視印刷物の単色ベタ部の濃度偏差
ΔR,ΔG,ΔB,ΔIr：  基準印刷物に対する監視印刷物のセンサ値の偏差
C0,M0,Y0,K0     ：  基準印刷物の単色ベタ部の濃度値
R0,G0,B0,Ir0    ：  基準印刷物の監視点のセンサ値
Ac,Am,Ay,Ak     ：  監視点のＣＭＹＫ網点面積率
とする。
【００７２】
なお、光学センサ２５Ｆ・２５Ｒは、上記のようにＲＧＢＩｒの４つの波長領域で撮像して、各波長領域のセンサ値（実測画像データ）を画素毎に出力するので、以下の説明では、４つの波長領域のセンサ値をまとめてＲＧＢＩｒ値と称する。
【００７３】
色調の管理のために必要なものは、基準印刷物に対する監視印刷物の単色ベタ部の濃度偏差ΔC,ΔM,ΔY,ΔKであり、これを基準印刷物の監視点のセンサ値（R0,G0,B0,Ir0）と、印刷物の監視点のセンサ値（R,G,B,Ir）から精度良く求めることができれば、色調の管理を問題なく実施できる。
【００７４】
なお、色調監視・管理の対象となる印刷物には、実際には単色ベタ部がない場合が多いので、この場合には上記のC,M,Y,Kは、単色ベタ部があったと仮定して単色ベタ部の濃度を表すことになる。
【００７５】
仮に、
C  =  C(R,G,B,Ir)
M  =  M(R,G,B,Ir)
Y  =  Y(R,G,B,Ir)
K  =  K(R,G,B,Ir)
を満たす関数C( ),M( ),Y( ),K()を精度良く求めることができれば、基準印刷物に対する監視印刷物の単色ベタ部の濃度偏差ΔC,ΔM,ΔY,ΔKは、
ΔC =  C(R,G,B,Ir)  −  C(R0,G0,B0,Ir0)
ΔM =  M(R,G,B,Ir)  −  M(R0,G0,B0,Ir0)
ΔY =  Y(R,G,B,Ir)  −  Y(R0,G0,B0,Ir0)
ΔK =  K(R,G,B,Ir)  −  K(R0,G0,B0,Ir0)
であり課題が解決できる。しかしながら、変数（変化の要因）が複数有るため、C( ),M( ),Y( )，K( )を精度良く求めることは困難である。そこで、以下のような方法を用いる。
【００７６】
ＣＭＹＫ面積率空間を複数の十分小さいブロックに分割すれば、任意のｊ番目のブロックに対して、以下のような関数を精度良く定義することが可能である。
【００７７】
ΔC   =   Cj(ΔR,ΔG,ΔB,ΔIr)
ΔM   =   Mj(ΔR,ΔG,ΔB,ΔIr)
ΔY   =   Yj(ΔR,ΔG,ΔB,ΔIr)
ΔK   =   Kj(ΔR,ΔG,ΔB,ΔIr)
各Cj( ),Mj( ),Yj( ),Kj()は、ΔR,ΔG,ΔB,ΔIrの、例えば１〜３次の多項式で表される。同様にRGBIｒ空間で複数の十分小さいブロックに分割することも可能である。
【００７８】
本発明では、上記のアルゴリズムを監視点毎に適用した結果を、コラム（インキキー幅）毎に平均してインキキー７２Ａ〜７２Ｈの制御を行う。なお、監視点毎の結果は、必ず平均化するものではなく、装置の構成等に応じて行えばよい。
【００７９】
各Cj( ),Mj( ),Yj( ),Kj()が、ΔR,ΔG,ΔB,ΔIrの１次式である場合の具体的な計算フローは、以下の通りである。
【００８０】
ΔC = Cj(ΔR,ΔG,ΔB,ΔIr) = ∂C/∂R×ΔR +∂C/∂G×ΔG +∂C/∂B×ΔB +∂C/∂Ir×ΔIr…（式１）
ΔM = Mj(ΔR,ΔG,ΔB,ΔIr) = ∂M/∂R×ΔR +∂M/∂G×ΔG +∂M/∂B×ΔB +∂M/∂Ir×ΔIr…（式２）
ΔY = Yj(ΔR,ΔG,ΔB,ΔIr) = ∂Y/∂R×ΔR +∂Y/∂G×ΔG +∂Y/∂B×ΔB +∂Y/∂Ir×ΔIr…（式３）
ΔK = Kj(ΔR,ΔG,ΔB,ΔIr) = ∂K/∂R×ΔR +∂K/∂G×ΔG +∂K/∂B×ΔB +∂K/∂Ir×ΔIr…（式４）
ここで、∂C/∂R等の変化率係数は監視点に依存する。つまり、変化率係数は、監視点毎に異なる値となる。
【００８１】
上記の式１〜式４を行列式で表すと、以下のようになる。
【００８２】
【数２】
0004275716.tif 000003
【００８３】
次に、色調監視装置２１の具体的な構成について説明する。図２は、色調監視装置の構成を示す図である。
【００８４】
色調監視装置２１は、多値化処理部４１、画像処理部４３、プロファイル記憶部４４、変化率係数画像作成部４５、変化率係数決定部４６、監視点座標出力部４７、補間付変化率ＬＵＴ（Look Up Table ）４８、差分演算部４９、積演算部５０、集計部５１、センサ値変換部５２、及び画像メモリ５３を備えている。
【００８５】
ここで、色調監視装置２１では、光学センサ２５Ｆ・２５Ｒが撮像して画素単位で出力した画像データを、１画素単位、または複数画素をまとめたものを１要素として、これを単位として、多値化処理部４１、画像処理部４３、及びセンサ値変換部５２は処理を行うことが可能である。この取りまとめ方法には、以下のようなものがある。
【００８６】
１．複数画素を平均化して１処理単位とする。例えば、縦３画素×横３画素のエリアの平均値を１処理単位として扱う。これにより、１画素単位のノイズの影響を緩和したり、処理速度を上げたりするメリットがある。
【００８７】
２．RGBIｒまたはCMYKの網点面積率が、ある一定の範囲内のものを１単位としてこの平均値を１処理単位とする。例えば、RGBIｒをそれぞれ０〜１００の値で表現した場合、これを５刻みで１つの範囲としてRGBIｒをそれぞれ２０の範囲に分割し、RGBIｒの組み合わせで２０×２０×２０×２０=１６００００に区分して、この同じ区分に入る画素を平均化して１処理単位として扱う。この方法により、均一な絵柄(平網部)がある絵柄に対しては処理速度及び精度を上げることができる。
【００８８】
また、以下の説明では、複数の監視点を、１画素または１要素単位に設定し、各監視点について、色調の監視・管理を行う場合を例に挙げて説明する。
【００８９】
多値化処理部４１は、画像サーバ１３から２値網点画像データを読み出して、この２値網点画像データから各画素の刷版上の網点面積率（Ａｃ，Ａｍ，Ａｙ，Ａｋ＝デジタルデータでの網点面積率）を作成して、変化率係数画像作成部４５へ出力する。
【００９０】
画像処理部４３は、画像サーバ１３からテストチャートの２値網点画像データを読み出して、この２値網点画像データと、画像メモリ５３から読み出した、光学センサ２５Ｆで読み取った色調調節済テストチャートの画像データ（ＲＧＢＩｒ値）と、の関係を分析して、プロファイルＰ３を作成する。そして、このプロファイルＰ３をプロファイル記憶部４４に記憶させる。また、画像処理部４３は、記憶部４４からプロファイルＰ３を読み出して、画像サーバ１３から読み出した２値網点画像データを多値化・プロファイル変換して作成したＲＧＢＩｒ値の基準値（基準ＲＧＢＩｒ値）を、差分演算部４９へ出力する。
【００９１】
プロファイル記憶部４４は、プロファイルＰ３を記憶する。
【００９２】
変化率係数画像作成部４５は、多値化処理部４１が出力した印刷物の網点面積率のデータと、補間付変化率ＬＵＴ４８から読み出してデータと、に基づいて変化率係数画像を作成し、この変化率係数画像を変化率係数決定部４６へ出力する。変化率係数画像は、基準画像の各画素と１対１対応で変化率係数が入っている画像である。
【００９３】
変化率係数決定部４６は、変化率係数画像を記憶する不図示のメモリを備えており、このメモリに記憶させた変化率係数画像から監視点の座標の変化率係数を読み出して、積演算部５０へ出力する。
【００９４】
監視点座標出力部４７は、画素毎の画像データを解析して、ＲＧＢＩｒの各色の複数の監視点を設定する。また、監視点座標出力部４７は、前記のように設定した複数の監視点の座標データを記憶しており、監視点の座標データを画像処理部４３、変化率係数決定部４６、及び差分演算部４９に出力する。なお、監視点としては、すべての画素を設定しても、一部の画素を選択して設定しても良い。
【００９５】
補間付変化率ＬＵＴ（Look UpTable ）４８は、Ｎ×Ｎ×Ｎ×Ｎの格子点を持ち、網点面積率の４つの値Ac,Am,Ay,Akを入力とし、それに対応した以下の１６個の値を出力する。
【００９６】
∂C/∂R    ∂C/∂G    ∂C/∂B    ∂C/∂Ir
∂M/∂R    ∂M/∂G    ∂M/∂B    ∂M/∂Ir
∂Y/∂R    ∂Y/∂G    ∂Y/∂B    ∂Y/∂Ir
∂K/∂R    ∂K/∂G    ∂K/∂B    ∂K/∂Ir
なお、補間付変化率ＬＵＴの詳細は後述する。
【００９７】
差分演算部４９は、画像メモリ５３から読み出した実測RGBIr値と、画像処理部４３が出力した基準RGBIｒ値と、の差を監視点毎に算出して、積演算部５０に出力する。
【００９８】
積演算部５０は、差分演算部４９が出力した実測ＲＧＢＩｒ値と基準ＲＧＢＩｒ値との差に、変化率係数決定部４６が出力した変化率係数を乗算することで、基準印刷物に対する監視印刷物の単色ベタ濃度差（ベタ濃度偏差）を算出して、集計部５１へ出力する。
【００９９】
集計部５１は、積演算部５０が出力した演算結果（ベタ濃度偏差）を、コラム（インキ量の制御単位幅であるインキキー幅）で平均して、コラム毎のベタ濃度偏差を算出し、その平均値に基づいた制御信号をインキキー制御装置２３に出力する。また、集計部５１は、差分演算部４９が出力した監視点毎のＲＧＢＩｒ値の実測値と基準値の差分、積演算部５０が出力したベタ濃度偏差、またはコラム毎のベタ濃度偏差の少なくともいずれかをモニタ２２が表示するように、これらの信号をモニタ２２に出力する。
【０１００】
モニタ２２は、集計部５１から受け取った信号に基づいた画像を表示する。
【０１０１】
センサ値変換部５２は、光学センサ２５Ｆが印刷物の画像を読み取って出力した光量値を上記用語定義によるセンサ値に変換し、画像メモリ５３に出力する。なお、色調監視装置２１において、センサ値として光学センサ２５Ｆ・２５Ｒが出力した反射光量値に対して対数変換処理などを行わずに用いる場合や、光学センサ２５Ｆ・２５Ｒが反射光量に基づく別の値（例えば、反射光量値を対物変換した値）を出力する場合には、センサ値変換部５４は設けなくても良い。また、センサ値変換部５２は、画像メモリ５３の下流側に設けても良い。
【０１０２】
画像メモリ５３は、光学センサ２５Ｆが印刷物の画像を読み取って出力したＲＧＢＩｒ値を記憶する。
【０１０３】
次に、色調監視装置２１の概略の動作について、図２〜図４に基づいて説明する。図３は、前処理の手順を示すフローチャートである。図４は、色調監視処理の手順を示すフローチャートである。
【０１０４】
＜前処理＞
色調監視装置２１は、初期設定時に、前処理として以下の処理を行う。
【０１０５】
まず、色調監視装置２１は、図２，図３（Ａ）に示すように、テストチャートを用いて色合わせを行い、色合わせが正常に完了した状態で印刷したテストチャート（印刷物）の画像を光学センサ（ＲＧＢＩｒセンサ）２５Ｆで読み取って（ｓ１）、ＲＧＢＩｒの各反射光量値をセンサ値変換部５２へ出力する（ｓ２）。センサ値変換部５２は、光学センサ２５Ｆが印刷物の画像を読み取って出力した反射光量値をセンサ値（ＲＧＢＩｒ値）に変換し、画像メモリ５３に出力し、画像メモリ５３は、この画像のセンサ値（ＲＧＢＩｒ値）を蓄積する（ｓ３）。
【０１０６】
画像処理部４３は、画像メモリ５３からセンサ値（ＲＧＢＩｒ値）を読み出すとともに、画像サーバ１３に蓄積されているテストチャートの画像データを読み出して（ｓ４）、双方の関係を画像処理部４３で分析して、プロファイル（ＩＣＣプロファイル）Ｐ３を作成する（ｓ５）。そして、画像処理部４３は、このプロファイルＰ３をプロファイル記憶部４４に記憶させる（ｓ６）。
【０１０７】
色調監視装置２１では、図２、図３（Ｂ）に示すように、多値化処理部４１が、画像サーバ１３に蓄積された画像ファイルのデジタルデータ（２値網点画像データ）を取得する（ｓ１１）。そして、多値化処理部４１は、２値網点画像データから各画素の刷版上の網点面積率（Ａｃ，Ａｍ，Ａｙ，Ａｋ＝デジタルデータでの網点面積率）を作成して、監視点座標出力部４７及び変化率係数画像作成部４５へ出力する（ｓ１２）。
【０１０８】
監視点座標出力部４７は、監視点の設定を行う。すなわち、画像データに対して、画素単位の大きさで、各コラム（インキキー幅）にＣＭＹＫの各色それぞれ複数個の監視点を設定する（ｓ１３）。
【０１０９】
一方、変化率係数画像作成部４５は、多値化処理部４１が出力した印刷物の網点面積率のデータと、補間付変化率ＬＵＴ４８から読み出したデータと、に基づいて、監視点の面積率に対応した変化率係数を求めて、変化率係数画像を作成し、変化率係数決定部４６へ変化率係数画像を出力する（ｓ１４）
変化率係数決定部４６は、変化率係数画像作成部４５から変化率係数画像を受け取ると、この画像を不図示のメモリに記憶させる（ｓ１５）。
【０１１０】
なお、補間付変化率ＬＵＴ４８は、後述する手順で別途作成する。
【０１１１】
また、変化率係数画像は、印刷物毎に必要であり、監視する印刷物が変われば、変化率係数画像を作り直す必要がある。また、監視点も同様に印刷物が変われば、作り直す必要がある。
【０１１２】
また、画像処理部４３が作成するプロファイルＰ３は、ＩＣＣプロファイルでも良いが、このプロファイルＰ３に代えて、４入力４出力の４次元ＬＵＴをプロファイル記憶部４４に記憶させても良い。
【０１１３】
また、色調監視装置２１では、変化率係数画像を作成せずに、印刷物の網点面積率のデータと、補間付変化率ＬＵＴ４８から読み出したデータと、に基づいて、監視点の面積率に対応した変化率係数を求めるようにすることも可能である。
【０１１４】
＜色調管理処理＞
輪転機１９は、オペレータによって印刷が開始されると、ロール紙Ｒから供給されるウェブ状の連続紙Ｐに、所定の画像として例えば新聞８ページ分の画像を連続して印刷する。
【０１１５】
図４に示すように、監視点座標出力部４７は、輪転機１９でこれから印刷する画像の監視点の座標を、画像処理部４３、変化率係数決定部４６、及び差分演算部４９に出力する。また、輪転機１９で印刷が開始されると、光学センサ２５Ｆは、連続紙Ｐに印刷された画像を、Ｒ・Ｇ・Ｂ・Ｉｒの４つの波長領域に分けて連続して撮像し、Ｒ・Ｇ・Ｂ・Ｉｒの各波長領域の画像データを画素単位で、センサ値変換部５２へ出力する（ｓ３１）。
【０１１６】
センサ値変換部５２は、光学センサ２５Ｆが印刷物の画像を読み取って出力した光量値を、上記用語定義によるセンサ値（ＲＧＢＩｒ値）に変換し、画像メモリ５３に出力する（ｓ３２）。
【０１１７】
画像メモリ５３は、センサ値変換部５２から送られてきたセンサ値（ＲＧＢＩｒ値）を、印刷の１周期分（例えば新聞８ページ分）記憶する（ｓ３３）。
【０１１８】
一方、画像処理部４３は、画像サーバ１３から輪転機１９で印刷中の画像の２値網点画像データを読み出して、記憶部４４からプロファイルＰ３（または４次元ＬＵＴ）に基づいて基準ＲＧＢＩｒ値を生成する（ｓ３４）。そして、画像処理部４３は、各監視点のＲＧＢＩｒ値（基準値）を差分演算部４９へ出力する（ｓ３５）。
【０１１９】
差分演算部４９は、画像処理部４３が出力したＲＧＢＩｒ値（基準値）と、画像メモリ５３から読み出したＲＧＢＩｒ値（実測値）と、を監視点毎に比較し、両ＲＧＢＩｒ値の差を積演算部５０へ出力する（ｓ３６）。
【０１２０】
変化率係数決定部４６では、監視点座標出力部４７が出力した監視点の座標情報に基づいて、変化率係数画像作成部４５から、各監視点の座標の変化率係数を読み出す（ｓ３７）。続いて、そして、変化率係数決定部４６は、監視点毎の変化率係数を積演算部５０へ出力する（ｓ３８）。
【０１２１】
積演算部５０は、差分演算部４９が出力したＲＧＢＩｒ値の差の値に対して、変化率係数決定部４６が出力した変化率係数を乗算することで、基準印刷物に対する監視印刷物の単色ベタ部の濃度偏差を監視点毎に算出する。そして、積演算部５０は、監視点毎のベタ濃度偏差を集計部５１に出力する（ｓ３９）。
【０１２２】
集計部５１は、積演算部５０が出力した監視点毎のベタ濃度偏差（演算結果）をコラム（インキ量の制御単位幅であるインキキー幅）毎に平均化し、コラム毎の単色ベタ濃度偏差を、インキキー制御装置２３に出力する（ｓ４０）。
【０１２３】
このとき、集計部５１は、コラム毎の平均ベタ濃度偏差をモニタ２２が表示するように、信号を出力させるように構成してもよい。また、集計部５１では、ベタ濃度偏差の平均化をコラム毎に行わずに、監視点毎のベタ濃度偏差をモニタ２２が表示するように、信号を出力させるように構成してもよい。さらに、前記の推定ベタ濃度値とベタ濃度偏差とをモニタ２２に表示するように集計部５１を構成してもよい。
【０１２４】
インキキー制御装置２３は、積演算部５０が出力したコラム毎の単色ベタ濃度の偏差に基づいて、各色のインキキー７２Ａ〜７２Ｈを制御して、輪転機１９におけるインキ量を調整する（ｓ４１）。
【０１２５】
これにより、輪転機１９で印刷した連続紙Ｐをドライダウンした後の新聞紙３３の色調（インキ濃度値）と、新聞社で決めた適正印刷状態の色調と、が結果的にほぼ一致するようになる。
【０１２６】
なお、上記の説明では、画像処理部４３は、画像サーバ１３から２値網点画像データを読み出して、ＲＧＢＩｒ値の基準値を作成する構成としたが、これに限るものではなく、基準印刷物を光学センサ２５Ｆで撮像した基準画像データをそのままＲＧＢＩｒ値の基準値とするようにしても良い。
【０１２７】
＜補間付変化率ＬＵＴの構成方法＞
次に、色調監視装置２１の補間付変化率ＬＵＴ４８の構成方法について説明する。
【０１２８】
変化率ＬＵＴは、Ｎ×Ｎ×Ｎ×Ｎの格子点を持ち、網点面積率の４つの値Ac,Am,Ay,Akを入力とし、それに対応した以下の１６個の値を出力する。
【０１２９】
∂C/∂R    ∂C/∂G    ∂C/∂B    ∂C/∂Ir
∂M/∂R    ∂M/∂G    ∂M/∂B    ∂M/∂Ir
∂Y/∂R    ∂Y/∂G    ∂Y/∂B    ∂Y/∂Ir
∂K/∂R    ∂K/∂G    ∂K/∂B    ∂K/∂Ir
以下、便宜上、Ｎ＝５として説明する。ＬＵＴの各格子点は、
Ac ＝０，２５，５０，７５，１００
Am ＝０，２５，５０，７５，１００
Ay ＝０，２５，５０，７５，１００
Ak ＝０，２５，５０，７５，１００
の全ての組み合わせ５×５×５×５＝６２５個である。
【０１３０】
前記のCj( ),Mj( ),Yj( ),Kj()をj=0〜624として扱うことも可能であるが、ブロックのサイズが比較的大きいので、更に補間付きＬＵＴを用いることにより精度を上げることが可能となる。
【０１３１】
ＬＵＴへの入力が、Ｃ＝２５，Ｍ＝５０，Ｙ＝７５、Ｋ＝２５の場合、この値にちょうど対応するＬＵＴの格子点があるので、その格子点に登録されている、１６個の出力値を出力すれば済むが、Ｃ＝４０，Ｍ＝５０，Ｙ＝７５，Ｋ＝２５などの場合、格子点にのらないので、Ｃ＝２５，Ｍ＝５０，Ｙ＝７５、Ｋ＝２５の格子点と、Ｃ＝５０，Ｍ＝７５，Ｙ＝７５、Ｋ＝２５の格子点に登録されている１６個の出力値を補間して出力を決めることにする。
【０１３２】
補間方法はプリズム補間、２次補間、３次補間、スプライン補間などいろいろな方法が適用可能であるが、入力値の近傍１６点から線形補間するのが簡便である。
【０１３３】
ＬＵＴの格子点の設定方法としては以下に、代表的な２つの方法について記載する。
【０１３４】
（１）変化率ＬＵＴの格子点に設定する係数を、印刷機モデルにより求める方法
印刷機モデルは、面積率Ac,Am,Ay,Akとベタ濃度Dc,Dm,Dy,Dkを入力として、R,G,B,Irを返す関数とみなせる。
【０１３５】
ここで、印刷機モデルの詳細について説明する。図５は、印刷機モデルの概略の構成図である。印刷機モデル演算部６０は、ベタ混色モデル演算部６１、反射率演算部６３、ドットゲインモデル演算部６５、面積率演算部６７、及びノイゲバウアモデル演算部６９を備えている。
【０１３６】
ベタ混色モデル演算部６１は、各インキの単色ベタ濃度から、各インキのベタ重ね部の濃度値を求める混色モデルに基づいて混色ベタ部の濃度を計算する。また、本発明では一例として、ベタ混色モデル演算部６１に、単色ベタ濃度をインキ膜厚に換算する膜厚モデルを採用したものについて説明する。
【０１３７】
反射率演算部６３は、単色ベタ部及び混色ベタ部の反射率を計算する。
【０１３８】
ドットゲインモデル演算部６５は、各インキのベタ濃度値及び網点面積率値から、各重ね部の面積率を求めるドットゲインモデルに基づいて、紙面上での網点面積値を計算する。
【０１３９】
面積率演算部６７は、単色ベタ部及び混色ベタ部の面積率を計算する。
【０１４０】
ノイゲバウアモデル演算部６９は、周知のノイゲバウアの式を用いて、反射率演算部６３が出力した単色ベタ部及び混色ベタ部の反射率と、面積率演算部６７が出力した単色ベタ部及び混色ベタ部の面積率と、から、印刷物の画像の推定ＲＧＢＩｒ値を算出する。
【０１４１】
ここで、ドットゲインモデル、混色モデル、膜厚モデル、及びノイゲバウアモデルについて説明する。図６は、ドットゲインモデル及び混色モデルを説明するための模式図である。
【０１４２】
＜ドットゲインモデル＞
製版工程においては、図６（Ａ）に示すように、インキが乗る部分と乗らない部分を小さなドットの集合で表現し、これをマクロに見ることにより階調及び色の違いを表現することを目的としている。このインキが乗る小さなドットの占める面積の割合を網点面積率という。オフセット印刷では、版から紙に直接ではなく、一旦ゴムを貼ったローラ(ブランケット胴)に転写してから、紙に転写するが、図６（Ｂ）に示すように、実際にインキを乗せて印刷した場合は、インキが立体的に乗る。そのため、図６（Ｃ）に示すように、その量により版の上に形成されたドットからインキがはみ出る現象が生じる。この場合には、版の上に形成された網点面積率よりも、実際に印刷された紙面上の網点面積率の方が大きくなるドットゲイン（網点の太り）と呼ばれる現象が生じる。
【０１４３】
このドットゲインの割合は、図６（Ｄ）に示すように、同じインキを使った輪転機であっても元の網点面積率とインキの量によって変わる。また、紙面に直接インキが乗る場合と、多色インキで印刷済みの部分の上にインキ乗る場合とではドットゲインが異なる現象がある。このような特性を、単色ベタ濃度と網点面積率から計算上で再現して、各インキの重ね印刷部の面積率を求めるものをドットゲインモデルと称する。なお、ドットゲインの発生原因は複数あり、その主な要因の一つが上記の版の上に形成されたドットからインキがはみ出る現象によるものである。
【０１４４】
＜混色モデル＞
例えばＣＭＹＫの４色で印刷した場合、それぞれの版はインキが乗るか乗らないかの２値で表現されたドットの集合であるから、印刷されたものもミクロで見ると、インキが乗っている部分は１００％のベタで印刷されていると見ることができる。
【０１４５】
したがって、ＣＭＹＫ４色の印刷の場合は、それぞれが１００％のベタで印刷された部分か、白紙の部分かの２通りに分けられるため、２通りの４乗＝１６通りに分けることができる。つまり、全ての箇所をミクロで見れば、図６（Ｅ）に示すように、白紙部、Ｃ単色ベタ部、Ｍ単色ベタ部、Ｙ単色ベタ部、Ｋ単色ベタ部、ＣＭ混色部、ＣＹ混色部、ＣＫ混色部、ＭＹ混色部、ＭＫ混色部、ＹＫ混色部、ＣＭＹ混色部、ＣＭＫ混色部、ＣＹＫ混色部、ＭＹＫ混色部、ＣＭＹＫ混色部の１６通りのベタ混色(白紙を含む)で表現できる。一方、印刷機モデルへの入力はＣＭＹＫの各面積率およびＣＭＹＫの各単色濃度であるので、Ｃ単色ベタ部乃至ＣＭＹＫ混色部の１５通りのベタ混色部のＲＧＢＩｒ各反射率乃至反射濃度を入力のＣＭＹＫ各単色濃度から求める必要があり、これを計算上で再現するものがベタ混色モデルである。
【０１４６】
次に、このベタ混色部の反射率を表現する方法の一例として、膜厚モデルを用いた場合について説明する。
【０１４７】
＜膜厚モデル＞
図７〜図１２は、膜厚モデルを説明するための図である。膜厚モデルは、混色モデルの一例である。膜厚モデルは、インキの膜厚がいくら増えても、ある一定濃度以上には上がらず、紙及びインキの表面での反射を考慮した場合も、この表面の反射が最大濃度を決定するという特性を計算上で再現するものである。
【０１４８】
膜厚モデルは、以下の式を基本とし、これに以下の工夫を加えている。
【０１４９】
Ｄ＝Ｄ∞・（１−ｅ−ｍｔ）
ｔ：膜厚
ｍ：紙、インキ等に依存する係数
Ｄ∞：最大濃度
［工夫１］
インキの刷り順により、単色ベタ部、混色部の膜厚の関係が変わる。そのため、単色ベタ部および混色部の膜厚を図７のように考える。
【０１５０】
［工夫２］
単色の光源で、モノクロカメラで見た時に見える濃淡画像は、単色インキでその画像と等価なものを再現できることを意味しているため、各光源(ここではR,G,B,Irとするが、この限りではない)毎に代表インキ(ここでは、Ｒ光源に対してはCインキ、G光源に対してはMインキ、B光源に対してはYインキ、Ir光源に対してはKインキとするが、この限りではない)を決め、その代表インキでの膜厚に換算する（図７〜１１参照）。
【０１５１】
［工夫３］
前記の式および、膜厚の関係から、図１２の関係をつなげることにより、RGBIr画像からYMCK単色ベタ部の濃度(インキ量)を求める。
【０１５２】
＜ノイゲバウアモデル＞
周知のノイゲバウアの式を用いて網点印刷の色再現を行うモデルであり、白地、単色インキ部、２次、３次、４次色インキベタ部分を対象色として計測し、混色モデルでのＲＧＢＩｒ各反射率に各色部分の出現確率を乗じて任意のＣＭＹＫ組み合わせに対する再現色を求める。
【０１５３】
印刷機モデル演算部６０は、単色ベタ部のセンサ濃度と、上流システム（画像サーバ１３）から取得する刷版作成用デジタルデータの網点面積率との２つを入力とし、予想センサ値（推定ＲＧＢＩｒ値）を計算して出力する。
【０１５４】
単色ベタ部濃度は、単色１００％のベタ部が印刷物中にあった場合にその部分の濃度のことで、実際の印刷物中には単色ベタ部が存在していない場合もある。単色ベタ部濃度は、インキの膜厚に相関する量である。
【０１５５】
なお、上記のドットゲインモデルやノイゲバウアモデルとしては、より精度の高い様々な修正式が提案されており、これらを適用することも可能であり、実施例に限定されるものではない。
【０１５６】
また、混色モデルとして膜厚モデル以外に、近似多項式を用いた方法等も適用可能である。
【０１５７】
近似多項式の例：  単色ベタ濃度  Dc_r,Dm_g, Dy_b, Dk_irが与えられたとき
Dc_g = C2c_r・Dc_r2+ C1c_r・Dc_r + C0c_r
Dc_b = C2c_b・Dc_r2+ C1c_b・Dc_r + C0c_b
Dc_ir= C2c_ir・Dc_r2+ C1c_ir・Dc_r + C0c_ ir
Dm_r= C2m_r・Dm_g2+ C1m_r・Dm_g + C0c_r
Dm_b= C2m_b・Dm_g2+ C1m_b・Dm_g + C0c_b
Dm_ir= C2m_ir・Dm_g2+ C1m_ir・Dm_g + C0c_ir
Dy_r= C2y_r・Dy_b2+ C1y_r・Dy_b + C0y_r
Dy_g= C2y_g・Dy_b2+ C1y_g・Dy_b + C0y_g
Dy_ir= C2y_ir・Dy_b2+ C1y_ir・Dy_b + C0y_ir
Dk_r= C2k_r・Dk_ir2+ C1k_r・Dk_ir + C0k_r
Dk_g= C2k_g・Dk_ir2+ C1k_g・Dk_ir + C0k_r
Dk_b= C2k_b・Dk_ir2+ C1k_b・Dk_ir + C0k_r
（以下、xにはr,g,b,irが入るものとする。）
Dcm_x=C2cm_x・Dc_x+ C1cm_x・Dm_x + C0cm_x・Dc_x・Dm_x
Dcy_x=C2cy_x・Dc_x+ C1cy_x・Dy_x + C0cy_x・Dc_x・Dy_x
Dck_x= C2ck_x・Dc_x+ C1ck_x・Dk_x + C0ck_x・Dc_x・Dk_x
Dmy_x=C2my_x・Dm_x+ C1my _x・Dy_x + C0my _x・Dm_x・Dy_x
Dmk_x=C2mk_x・Dm_x+ C1mk_x・Dk_x + C0mk_x・Dm_x・Dk_x
Dyk_x=C2yk_x・Dy_x+ C1yk_x・Dk_x + C0yk_x・Dy_x・Dk_x
Dcmy_x=C7cmy_x・Dc_x+ C6cmy_x・Dm_x + C5cmy_x・Dy_x +
C4cmy_x・Dcm_x+ C3cmy_x・Dcy_x + C2cmy_x・Dmy_x+
C1cmy_x・Dc_x・D m_x・D y_x + C0cmy_x
Dcmk_x=C7cmk_x・Dc_x+ C6cmk_x・Dm_x + C5cmk_x・Dk_x +
C4cmk_x・Dcm_x+ C3cmk_x・Dck_x + C2cmk_x・Dmk_x+
C1cmk_x・Dc_x・Dm_x・Dk_x + C0cmk_x
Dcyk_x=C7cyk_x・Dc_x+ C6cyk_x・Dy_x + C5cyk_x・Dk_x +
C4cyk_x・Dcy_x+ C3cyk_x・Dck_x + C2cyk_x・Dyk_x+
C1cyk_x・Dc_x・Dy_x・Dk_x + C0cyk_x
Dmyk_x=C7myk_x・Dm_x+ C6myk_x・Dy_x + C5myk_x・Dk_x +
C4myk_x・Dmy_x+ C3myk_x・Dmk_x + C2myk_x・Dk_x+
C1myk_x・Dm_x・Dy_x・Dk_x + C0myk_x
Dcmyk_x=C15_x・Dc_x+ C14_x・Dm_x + C13_x・Dy_x +C12_x・Dk_x +
C11_x・Dcm_x+ C10_x・Dcy_x + C9_x・Dck_x +
C8_x・Dmy_x+ C7_x・Dmk_x + C6_x・Dyk_x +
C5_x・Dcmy_x+ C4_x・Dcmk_x + C3_x・Dcyk_x +C2_x・Dmyk_x +
C1_x・Dc_x・Dm_x・Dy_x・Dk_x +C0_x
近似多項式の係数（C1,C2,C3などで始まる係数）は、回帰分析手法で実際の印刷物にフィッティングさせて決める。近似多項式を用いた方法では、多項式の次数を変えることにより任意の精度を得ることが可能であり、モデルを仮定しないので特定の印刷方式に依存することがなく、あらゆる印刷方式に適用可能である。
【０１５８】
このように構成した印刷機モデルを用いることで、以下のように、補間付変化率ＬＵＴ４８の格子点に設定する変化率係数求めることができる。式にすると、
R  = PM_R(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk)
G  = PM_G(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk)
B  = PM_B(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk)
Ir =PM_Ir(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk)
となるので、数値微分法を用いて次のようにして求める。なお、下例では、1次の中心差分公式であるが、数値微分の計算方法については特に限定しない。
【０１５９】
面積率Ac,Am,Ay,Ak、標準印刷濃度をDc,Dm,Dy,Dk、微小な濃度差をdとして、
∂R/∂C=(PM_R(Ac,Am,Ay,Ak,Dc+d,Dm,Dy,Dk)-PM_R(Ac,Am,Ay,Ak,Dc-d,Dm,Dy,Dk))/ (2d)
∂G/∂C=(PM_G(Ac,Am,Ay,Ak,Dc+d,Dm,Dy,Dk)-PM_G(Ac,Am,Ay,Ak,Dc-d,Dm,Dy,Dk))/ (2d)
∂B/∂C=(PM_B(Ac,Am,Ay,Ak,Dc+d,Dm,Dy,Dk)-PM_B(Ac,Am,Ay,Ak,Dc-d,Dm,Dy,Dk))/ (2d)
∂Ir/∂C=(PM_Ir(Ac,Am,Ay,Ak,Dc+d,Dm,Dy,Dk) -PM_Ir(Ac,Am,Ay,Ak,Dc-d,Dm,Dy,Dk))/ (2d)
∂R/∂M=(PM_R(Ac,Am,Ay,Ak,Dc,Dm+d,Dy,Dk)-PM_R(Ac,Am,Ay,Ak,Dc,Dm-d,Dy,Dk))/ (2d)
∂G/∂M=(PM_G(Ac,Am,Ay,Ak,Dc,Dm+d,Dy,Dk)-PM_G(Ac,Am,Ay,Ak,Dc,Dm-d,Dy,Dk))/ (2d)
∂B/∂M=(PM_B(Ac,Am,Ay,Ak,Dc,Dm+d,Dy,Dk)-PM_B(Ac,Am,Ay,Ak,Dc,Dm-d,Dy,Dk))/ (2d)
∂Ir/∂M=(PM_Ir(Ac,Am,Ay,Ak,Dc,Dm+d,Dy,Dk) -PM_Ir(Ac,Am,Ay,Ak,Dc,Dm-d,Dy,Dk)) / (2d)
∂R/∂Y=(PM_R(Ac,Am,Ay,Ak,Dc,Dm,Dy+d,Dk)-PM_R(Ac,Am,Ay,Ak,Dc,Dm,Dy-d,Dk))/ (2d)
∂G/∂Y=(PM_G(Ac,Am,Ay,Ak,Dc,Dm,Dy+d,Dk)-PM_G(Ac,Am,Ay,Ak,Dc,Dm,Dy-d,Dk))/ (2d)
∂B/∂Y=(PM_B(Ac,Am,Ay,Ak,Dc,Dm,Dy+d,Dk)-PM_B(Ac,Am,Ay,Ak,Dc,Dm,Dy-d,Dk))/ (2d)
∂Ir/∂Y=(PM_Ir(Ac,Am,Ay,Ak,Dc,Dm,Dy+d,Dk) -PM_Ir(Ac,Am,Ay,Ak,Dc,Dm,Dy-d,Dk)) / (2d)
∂R/∂K=(PM_R(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk+d)- PM_R(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk-d))/(2d)
∂G/∂K=(PM_G(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk+d)-PM_G(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk-d))/ (2d)
∂B/∂K=(PM_B(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk+d)-PM_B(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk-d))/ (2d)
∂Ir/∂K=(PM_Ir(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk+d) - PM_Ir(Ac,Am,Ay,Ak,Dc,Dm,Dy,Dk-d))/ (2d)
を計算しておき、
∂R /∂C    ∂R /∂M  ∂R /∂Y    ∂R /∂K
∂G /∂C    ∂G /∂M  ∂G /∂Y    ∂G /∂K
∂B /∂C    ∂B /∂M  ∂B /∂Y    ∂B /∂K
∂Ir/∂C    ∂Ir/∂M  ∂Ir/∂Y    ∂Ir/∂K
の逆行列を計算することで、
∂C/∂R    ∂C/∂G    ∂C/∂B    ∂C/∂Ir
∂M/∂R    ∂M/∂G    ∂M/∂B    ∂M/∂Ir
∂Y/∂R    ∂Y/∂G    ∂Y/∂B    ∂Y/∂Ir
∂K/∂R    ∂K/∂G    ∂K/∂B    ∂K/∂Ir
を求めることができる。この計算をＬＵＴの格子点分実施することにより、ＬＵＴの格子点の係数が設定できる。
【０１６０】
ここで、実際には、逆行列を求める前の行列は、４行４列の正方行列にならない場合があり、例えばAc=0のパッチでは、Ｃのインキがないので、
∂R /∂M    ∂R /∂Y    ∂R /∂K
∂G /∂M    ∂G /∂Y    ∂G /∂K
∂B /∂M    ∂B /∂Y    ∂B /∂K
∂Ir/∂M    ∂Ir/∂Y    ∂Ir/∂K
のように、４行３列になってしまう。このような場合、一般化逆行列で計算が可能である。
【０１６１】
また、４色を同時に含まない（１色、２色、３色だけで印刷されている）パッチでは、１行１列、２行２列、３行３列の正方行列として扱い、逆行列を求めることで、十分な精度が得られる場合がある。
【０１６２】
また、
∂C/∂R    ∂C/∂G    ∂C/∂B    ∂C/∂Ir
∂M/∂R    ∂M/∂G    ∂M/∂B    ∂M/∂Ir
∂Y/∂R    ∂Y/∂G    ∂Y/∂B    ∂Y/∂Ir
∂K/∂R    ∂K/∂G    ∂K/∂B    ∂K/∂Ir
についても、Ac,Am,Ay,Akの組み合わせによっては、どれかの列は無視した方が好ましい結果が得られる場合があり、そのような列は強制的に０にする場合がある。
【０１６３】
なお、光学センサ２５Ｆで撮像する画像の波長領域数がＭで、印刷使用するインキ数がＮの場合には、下記のＮ×Ｍ行列
∂ink1/∂color1,∂ink1/∂Color2,・・・,∂ink1/∂ColorM
∂ink2/∂Color1,∂ink2/∂Color2,・・・,∂ink2/∂ColorM
・・・・・
∂inkN/∂Color1,∂inkN/∂Color2,・・・,∂inkN/∂ColorM
を求めることができる。この計算をＬＵＴの格子点分実施することにより、ＬＵＴの格子点の係数が設定できる。
【０１６４】
そして、実測画像データと基準画像データの差分ΔＣｏｌｏｒ１，ΔＣｏｌｏｒ２，・・・，ΔＣｏｌｏｒＭに対して、各インキの変化率係数の行列である上記のＮ×Ｍ行列を乗算することで、各インキのベタ濃度偏差Δｉｎｋ１，Δｉｎｋ２，・・・，ΔｉｎｋＮをインキの補正量として算出することができる。
【０１６５】
また、実測画像データと基準画像データの差分ΔＣｏｌｏｒ１，ΔＣｏｌｏｒ２，・・・，ΔＣｏｌｏｒＭに対して、各インキの変化率係数の行列である上記のＮ×Ｍ行列を乗算して、各インキのベタ濃度偏差Δｉｎｋ１，Δｉｎｋ２，・・・，ΔｉｎｋＮを算出し、これに基準画像のベタ濃度値、基準ｉｎｋ１、基準ｉｎｋ２、・・・、基準ｉｎｋＮを加算したものから、別途設定された目標濃度値、目標ｉｎｋ１、目標ｉｎｋ２、・・・、目標ｉｎｋＮを減算することで、インキの補正量として算出することができる。監視画像の濃度と任意に指定された目標濃度との差分を、新たにΔｉｎｋとして出力することにより、基準画像の濃度とは異なる目標濃度に合わせた印刷が可能になる。
【０１６６】
（２）変化率ＬＵＴの格子点に設定する係数を、実印刷データに対して回帰分析手法を用いて求める方法
図１３は、カラーチャートの一例、及びカラーチャートの濃度を変化させる方法を説明するための図である。ＬＵＴの格子点の６２５色を含むカラーチャートの濃度を変化させて印刷する。カラーチャートとしては、例えば図１３（Ａ）に示すようなものを使用する。このカラーチャートには、ＬＵＴの格子点の６２５パッチとＣＭＹＫのベタパッチを設けている。
【０１６７】
カラーチャートの濃度を変化させる方法としては、例えば、図１３（Ｂ）に示すように、標準濃度の画像に対して、Ｃ・Ｍ・Ｙ・Ｋの色毎に濃度を±０．１変化させて印刷する。すなわち、ステップ０では、ＣＭＹＫの４色とも標準濃度に合わせる。ステップ１では、Ｃの濃度を標準濃度から０．１下げる。ステップ２では、Ｃの濃度を標準濃度に戻し、Ｍの濃度を標準濃度から０．１下げる。ステップ３では、Ｍの濃度を標準濃度に戻し、Ｙの濃度を標準濃度から０．１下げる。ステップ４では、Ｙの濃度を標準濃度に戻し、Ｋの濃度を標準濃度から０．１下げる。ステップ５では、Ｋの濃度を標準濃度に戻し、Ｃの濃度を標準濃度から０．１上げる。ステップ６では、Ｃの濃度を標準濃度に戻し、Ｍの濃度を標準濃度から０．１上げる。ステップ７では、Ｍの濃度を標準濃度に戻し、Ｙの濃度を標準濃度から０．１上げる。ステップ８では、Ｙの濃度を標準濃度に戻し、Ｋの濃度を標準濃度から０．１上げる。ステップ９では、Ｋの濃度を標準濃度に戻し、ＣＭＹＫの４色とも標準濃度に合わせる。そして、上記のように濃度を変化させた画像を光学センサ２５Ｆで撮像して、得られたＲＧＢＩｒ値を用いて、以下のように演算する。
【０１６８】
上記印刷テスト中に、一定時間間隔で光学センサ画像の取得を行う。以下の説明では、１００点のセンサ画像が得られたとして、説明を続ける。（各センサ画像には、画像を取得した時刻の印刷物が写真撮影したように写っている。）
光学センサ画像の番号をi  （i=0〜99）
ＬＵＴの格子点の６２５パッチの番号をj  （j=0〜624）
パッチ番号jのパッチと同じ列に属するＣＭＹＫベタパッチの番号をxj（xj=0〜19）
i番目の画像のｊ番目の格子点パッチの光学センサ２５Ｆの値をRij,Gij,Bij,Irij
i番目の画像のxj番目のＣＭＹＫベタパッチの濃度をCixj,Mixj,Yixj,Kixj
光学センサ２５Ｆの画像の０番を基準画像
ΔRij = Rij -R0j
ΔGij = Gij -G0j
ΔBij = Bij -B0j
ΔIrij = Irij - Ir0j
ΔCij = Cixj - C0xj
ΔMij = Mixj - M0xj
ΔYij = Yixj - Y0xj
ΔKij = Kixj - K0xj
としたとき、各パッチｊについて
ΔCij ≒  a00×ΔRij + a01×ΔGij+ a02×ΔBij + a03×ΔIrij + a04
ΔMij ≒  a10×ΔRij + a11×ΔGij+ a12×ΔBij + a13×ΔIrij + a14
ΔYij ≒  a20×ΔRij + a21×ΔGij+ a22×ΔBij + a23×ΔIrij + a24
ΔKij ≒  a30×ΔRij + a31×ΔGij+ a32×ΔBij + a33×ΔIrij + a34
が、全てのiについてなるべく良い精度で成立するa00,〜,a34を求めることが目標となる。
【０１６９】
その解法としては、最小二乗法、正規方程式、最小絶対値偏差法、最小距離法等に関する文献及びその改良方法が複数知られているので、解き方の詳細は割愛する。
【０１７０】
ここでは、４変数１次式にフィッティングさせる最小二乗法を回帰分析手法の例として取り上げて簡単に説明する。
【０１７１】
【数３】
0004275716.tif 000004
【０１７２】
としたときの、ECjを最小にするa00,a01,a02,a03,a04を求めるために、上式をa00,a01,a02,a03,a04でそれぞれ偏微分し、極値である条件より、
∂ECj/∂a00 = 2Σ(a00×ΔRij +a01×ΔGij + a02×ΔBij + a03×ΔIrij + a04 -ΔCij) ×ΔRij=0
∂ECj/∂a00 = 2Σ(a00×ΔRij +a01×ΔGij + a02×ΔBij + a03×ΔIrij + a04 -ΔCij) ×ΔRij= 0
∂ECj/∂a01 = 2Σ(a00×ΔRij +a01×ΔGij + a02×ΔBij + a03×ΔIrij + a04 -ΔCij) ×ΔGij= 0
∂ECj/∂a02 = 2Σ(a00×ΔRij +a01×ΔGij + a02×ΔBij + a03×ΔIrij + a04 -ΔCij) ×ΔBij= 0
∂ECj/∂a03 = 2Σ(a00×ΔRij +a01×ΔGij + a02×ΔBij + a03×ΔIrij + a04 -ΔCij) ×ΔIrij= 0
∂ECj/∂a04 = 2Σ(a00×ΔRij +a01×ΔGij + a02×ΔBij + a03×ΔIrij + a04 -ΔCij) = 0
が得られ、これより解くべき連立方程式は掛算記号×を省略して記載すると、
a00Σ(ΔRij x0_i)+a01Σ(ΔGijΔRij)+a02Σ(ΔBijΔRij)+a03Σ(ΔIrijΔRij)+a04ΣΔRij -Σ(ΔCijΔRij)=0
a00Σ(ΔRij x1_i)+a01Σ(ΔGijΔGij)+a02Σ(ΔBijΔGij)+a03Σ(ΔIrijΔGij)+a04ΣΔGij-Σ(ΔCijΔGij)=0
a00Σ(ΔRij ΔBij)+a01Σ(ΔGijΔBij)+a02Σ(ΔBijΔBij)+a03Σ(ΔIrijΔBij)+a04ΣΔBij-Σ(ΔCijΔBij)=0
a00Σ(ΔRij ΔIrij)+a01Σ(ΔGijΔIrij)+a02Σ(ΔBijΔIrij)+a03Σ(ΔIrijΔIrij)+a04ΣΔIrij-Σ(ΔCijΔIrij)=0
a00ΣΔRij + a01ΣΔGij + a02ΣΔBij + a03ΣΔIrij + 100a04 - ΣΔCij = 0
の５元連立一次方程式となり、a00,a01,a02,a03,a04を連立方程式の解として求めることができる。（この解き方は多くの書物に記されているので省略する。）
以下、同様に、
【０１７３】
【数４】
0004275716.tif 000005
【０１７４】
として、
Emを最小にするa10,a11,a12,a13,a14、Eyを最小にするa20,a21,a22,a23,a24、Ekを最小にするa30,a31,a02,a33,a34を求めることができる。
【０１７５】
こうして求めったa00,〜,a34のうち、a04,a14,a24,a34は定数項であり、今回の場合、原理的に０または０に非常に近いはずのものであるため使用しない。
【０１７６】
j番目のパッチの網点面積率をAcj,Amj,Ayj,Akjとすると、ＬＵＴのAcj,Amj,Ayj,Akjの格子点のデータを、
∂C/∂R = a00    ∂C/∂G = a01    ∂C/∂B = a02    ∂C/∂Ir = a03
∂M/∂R = a10    ∂M/∂G = a11    ∂M/∂B = a12    ∂M/∂Ir = a13
∂Y/∂R = a20    ∂Y/∂G = a21    ∂Y/∂B = a22    ∂Y/∂Ir = a23
∂K/∂R = a30    ∂K/∂G = a31    ∂K/∂B = a32    ∂K/∂Ir = a33
のように設定する。
【０１７７】
以上のプロセスを各j（0 〜624）について実施することにより、ＬＵＴの全ての格子点に係数を設定することができる。
【０１７８】
上記は、各色独立に濃度を上げ下げし、その時に取得した連続した画像を基にして回帰分析手法で係数を求めているが、最終的にインキの上げ下げについて独立な結果が得られるのであれば、インキの濃度上げ下げの順序および組み合わせについてはこれに限定するものではない、例えば、CMを上げ、これを戻しMYを上げ、これを戻し、YKを上げ、これを戻しKCを上げ．．．といった組み合わせも可能である。
【０１７９】
また、連続した画像の回帰分析手法でなく、各濃度の上げ下げ時に各1個のデータを取得し、これによりデータを取得するといった手法も可能であるが、より精度を上げるためには本実施例のような形態をとった方が望ましい、しかし実施形態についてはこれに限定するものではない。
【０１８０】
なお、上記の説明では、変化率ＬＵＴの格子点に設定する係数を、実印刷データに対して求める場合について説明したが、これに限るものではなく、画像処理部４３が画像サーバ１３から読み出した２値網点画像データ（基準画像データ）に対して、変化率ＬＵＴの格子点に設定する係数を求めるようにしても良い。
【０１８１】
なお、補間付変化率ＬＵＴは、通常は輪転機または工場単位で１つ作成しても良いし、複数作成するようにしても良い。
【０１８２】
また、上記の説明では、変化率係数決定部４６の入力となる変化率係数画像は、画像サーバ１３から取得した上流データに基づいて、変化率係数画像作成部４５で作成したものを使用する場合について説明したが、これに限るものではなく、例えば、光学センサ２５Ｆが撮像した画像から推定して作成したものでも構わない。
【０１８３】
また、上流システムから画像ファイルが取得できない場合などには、図５（Ｂ）に示すように、初期設定時に、印刷機モデル演算部に対して、単色ベタ部濃度を与えるとともに、仮の網点面積率を与えて、網点面積率が一定値以下に収束するように、反復計算することで、ＣＭＹＫ網点の面積率の値が得られる。
【０１８４】
なお、以上の説明では、オフセット輪転機によりＣＭＹＫの４色のインキを用いて印刷した印刷物の色調監視や色調管理について説明したが、前記のヘキサクローム印刷やハイファイ印刷等のように、５色以上の多色印刷においても、それに合った印刷機モデルを作成することにより、印刷物の色調監視や色調管理が同様に可能である。
【０１８５】
すなわち、上記の説明ではＩｎｋ１＝Ｃ、Ｉｎｋ２＝Ｍ、Ｉｎｋ３＝Ｙ、Ｉｎｋ４＝Ｋの４色プロセス印刷の場合について説明を行ったが、５色以上使用している印刷、一例として、Ｉｎｋ１，Ｉｎｋ２，Ｉｎｋ３，Ｉｎｋ４、Ｉｎｋ５，Ｉｎｋ６の６色印刷の場合においても、４色以下の重ねで印刷されている領域がほとんどであるので、そのような部分に監視点を設定することにより同様に色調を監視することが可能である。
【０１８６】
例えば、高々Ｉｎｋ１，Ｉｎｋ２，Ｉｎｋ４，Ｉｎｋ６の４色重ねで印刷されている監視点については、上記の説明において、ＣをＩｎｋ１に、ＭをＩｎｋ２に、ＹをＩｎｋ４に、ＫをＩｎｋ６に置き換えればよい。この場合、印刷機モデル及びＬＵＴはＩｎｋ１乃至Ｉｎｋ６の６色対応とするか、または、６色中から４色を選び出す組合せは６通りあるので、上記説明のような４色対応のプリンタモデルおよびＬＵＴを６組用意しておき、切り替えて使用する。また、ＣＭＹＫと１つまたは複数の特色（スポットカラー、特練色とも言う。）インキによる特色印刷の場合も同様である。
【０１８７】
また、印刷システム１において、印刷物の色調の調整をオペレータが行う場合には、色調監視装置２１がインキキー制御装置２３に対して制御信号を出力する構成を備えていなくてもよい。例えば、各インキの監視点毎のＲＧＢＩｒ値の実測値と基準値の差分を確認できればよい場合には、色調監視装置２１の変化率係数決定部４６及び積演算部５０の代わりに、不図示の報知部を設けて、差分演算部４９が出力した監視点毎のＲＧＢＩｒ値の実測値と基準値の差分をモニタ２２に表示させる信号を出力するように構成するとよい。
【０１８８】
また、色調監視装置２１の集計部５１の代わりに、不図示の報知部を設けて、積演算部５０が算出した監視点毎のベタ濃度偏差、または差分演算部４９が出力した監視点毎のＲＧＢＩｒ値の実測値と基準値の差分の少なくとも一方をモニタ２２に表示させる信号を出力するように構成してもよい。
【０１８９】
さらに、集計部５１がインキキー制御装置２３に対して制御信号を出力しない構成にするとともに、集計部５１が、コラム毎のベタ濃度偏差、積演算部５０が算出した監視点毎のベタ濃度偏差、または差分演算部４９が出力した監視点毎のＲＧＢＩｒ値の実測値と基準値の差分の少なくともいずれかをモニタ２２に表示させる信号を出力するように構成してもよい。
【０１９０】
また、印刷システム１において、オペレータによる監視が不要の場合には、モニタ２２を設けず、集計部５１がモニタ２２に対して信号を出力しないようにしてもよい。
【０１９１】
また、以上の説明では、印刷機の一例としてオフセット輪転機を例に挙げて説明したが、本願発明はこれに限るものではなく他の印刷機で印刷した印刷物の色調を監視（管理）する場合にも、当然適用できる。例えば、インクジェットプリンタやレーザプリンタで印刷する場合にも、その装置にあった印刷機モデルを用いることにより、同様に色調を監視することが可能である。
【図面の簡単な説明】
【０１９２】
【図１】本発明の実施形態に係る色調管理システムの概略構成を示したブロック図である。
【図２】色調監視装置の構成を示す図である。
【図３】前処理の手順を示すフローチャートである。
【図４】色調監視処理の手順を示すフローチャートである。
【図５】印刷機モデルの概略の構成図である。
【図６】ドットゲインモデル及び混色モデルを説明するための模式図である。
【図７】膜厚モデルを説明するための図である。
【図８】膜厚モデルを説明するための図である。
【図９】膜厚モデルを説明するための図である。
【図１０】膜厚モデルを説明するための図である。
【図１１】膜厚モデルを説明するための図である。
【図１２】膜厚モデルを説明するための図である。
【図１３】カラーチャートの一例、及びカラーチャートの濃度を変化させる方法を説明するための図である。
【符号の説明】
【０１９３】
１…印刷システム    ３…色調管理システム    １１…紙面編集装置    １２…ＲＩＰ    １３…画像サーバ    １４…画像処理装置    １５…校正用プリンタ    １９…輪転機    ２１…色調監視装置    ２２…モニタ    ２３…インキキー制御装置    ２５Ｆ，２５Ｒ…光学センサ    ３１…校正紙    ３３…新聞紙    ４１…多値化処理部    ４３…画像処理部    ４４…プロファイル記憶部    ４５…変化率係数画像作成部    ４６…変化率係数決定部    ４７…監視点座標出力部    ４８…補間付変化率ＬＵＴ    ４９…差分演算部    ５０…積演算部    ５１…集計部    ５２…センサ値変換部    ５３…画像メモリ    ６１…ベタ混色モデル演算部    ６３…反射率演算部    ６５…ドットゲインモデル演算部    ６７…面積率演算部    ６９…ノイゲバウアモデル演算部    ７２Ａ〜７２Ｈ…インキキー    ７３…印刷ユニット
【図１】
0004275716.tif 000006
【図２】
0004275716.tif 000007
【図３】
0004275716.tif 000008
【図４】
0004275716.tif 000009
【図５】
0004275716.tif 000010
【図６】
0004275716.tif 000011
【図７】
0004275716.tif 000012
【図８】
0004275716.tif 000013
【図９】
0004275716.tif 000014
【図１０】
0004275716.tif 000015
【図１１】
0004275716.tif 000016
【図１２】
0004275716.tif 000017
【図１３】
0004275716.tif 000018