(57)【特許請求の範囲】
【請求項１】
画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する装置であって、
前記連続調画像信号に紐付けられた前記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得する閾値候補取得部と、
前記２値画像信号の作成後から前記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得する応答特性取得部と、
前記応答特性取得部により取得された前記応答特性を作用することで、前記網点形状の変化が模擬的に再現された模擬画像信号を作成する模擬画像作成部と、
前記模擬画像作成部により作成された前記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出する定量値算出部と、
前記閾値候補取得部により取得された前記各閾値候補を適用し、前記定量値算出部にてそれぞれ算出された前記定量値に基づく比較及び／又は評価によって、前記各閾値候補の中から設定対象である閾値データを決定する閾値データ決定部と
を備えることを特徴とする閾値データ設定装置。
【請求項２】
請求項１記載の閾値データ設定装置において、
前記模擬画像作成部は、平網領域を少なくとも１つ含む前記模擬画像信号を作成し、
前記定量値算出部は、前記模擬画像信号の前記平網領域でのパワースペクトルに関し、前記各網点に対応する成分を除いた主成分に基づいて前記定量値を算出する
ことを特徴とする閾値データ設定装置。
【請求項３】
請求項２記載の閾値データ設定装置において、
前記模擬画像作成部は、階調レベルが異なる平網領域を２つ以上含む前記模擬画像信号を作成し、
前記定量値算出部は、各前記平網領域における前記主成分についての統計値に基づいて前記定量値を算出する
ことを特徴とする閾値データ設定装置。
【請求項４】
請求項１〜３のいずれか１項に記載の閾値データ設定装置において、
前記画像形成過程にて前記網点形状に変化を及ぼす事象に応じた前記応答特性を推定する応答特性推定部を更に備え、
前記応答特性取得部は、前記応答特性推定部により推定された前記応答特性を取得する
ことを特徴とする閾値データ設定装置。
【請求項５】
請求項４記載の閾値データ設定装置において、
前記応答特性推定部は、刷版を介して被印刷材に色材を転移させる際の前記応答特性を推定することを特徴とする閾値データ設定装置。
【請求項６】
請求項５記載の閾値データ設定装置において、
前記応答特性推定部は、版材の主面を彫刻することで前記刷版を形成する際の前記応答特性を推定することを特徴とする閾値データ設定装置。
【請求項７】
請求項１〜６のいずれか１項に記載の閾値データ設定装置において、
前記閾値データ決定部により決定された前記閾値データを設定することで、前記連続調画像信号を前記２値画像信号に変換する２値化処理部を更に備えることを特徴とする閾値データ設定装置。
【請求項８】
請求項７記載の閾値データ設定装置と、
前記閾値データ設定装置により変換された前記２値画像信号に基づいて前記画像を形成する画像形成装置と
を備えることを特徴とする画像形成システム。
【請求項９】
画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する方法であって、
前記連続調画像信号に紐付けられた前記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得するステップと、
前記２値画像信号の作成後から前記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得するステップと、
取得された前記応答特性を作用することで、前記網点形状の変化が模擬的に再現された模擬画像信号を作成するステップと、
作成された前記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出するステップと、
取得された前記各閾値候補を適用し、それぞれ算出された前記定量値に基づく比較及び／又は評価によって、前記各閾値候補の中から設定対象である閾値データを決定するステップと
をコンピュータに実行させることを特徴とする閾値データ設定方法。
【請求項１０】
画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定するためのプログラムであって、
前記連続調画像信号に紐付けられた前記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得するステップと、
前記２値画像信号の作成後から前記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得するステップと、
取得された前記応答特性を作用することで、前記網点形状の変化が模擬的に再現された模擬画像信号を作成するステップと、
作成された前記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出するステップと、
取得された前記各閾値候補を適用し、それぞれ算出された前記定量値に基づく比較及び／又は評価によって、前記各閾値候補の中から設定対象である閾値データを決定するステップと
をコンピュータに実行させることを特徴とする閾値データ設定プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、画像を構成する各網点の形状・配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する閾値データ設定装置、方法及びプログラム、並びに画像形成システムに関する。
【背景技術】
【０００２】
印刷分野において、中間成果物を生成せずに電子データにより直接刷版を形成するＣＴＰ（Computer to Plate）が普及している。通常、刷版を形成する過程にて、連続調である画像信号に対して、組織的ディザ法を含むスクリーニング処理が実行される。従来から、高画質の印刷物を得るために、印刷システムの全体又は一部の特性を考慮した閾値データを決定する手法が種々提案されている。
【０００３】
特許文献１では、（１）ドットの擬似発生、及び（２）該ドットを含む濃度パターンの特定周波数成分の評価、を順次行うことで、モアレ等の発生が起こり難い閾値配列を決定する方法が提案されている。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００２−３６８９９５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
近時、使用可能な印刷装置又は印刷材（色材、印刷媒体、刷版等）の種類が多種多様になっている。例えば、印刷装置の種類として、オフセット印刷機のみならず、グラビア印刷機、フレキソ印刷機、デジタル印刷機（オンデマンド印刷機）等が存在する。また、印刷媒体の種類として、合成紙・厚紙・アルミ蒸着紙等の紙類のみならず、塩化ビニル・ＰＥＴ（ポリエチレンテレフタラート）等の樹脂、ターポリン、金属シート等が存在する。
【０００６】
ところで、一般的には、使用される印刷装置・印刷材の組み合わせに応じて、印刷特性が変化することが知られている。例えば、オフセット印刷・フレキソ印刷では、版材の種類（特に、弾性特性）が異なるため、ドットゲインが変化する程度が異なる傾向がみられる。また、刷版が不要なデジタル印刷機の場合、当該刷版の形成に起因する画質変化を何ら考慮する必要はない。
【０００７】
しかし、特許文献１に記載の方法によれば、印刷装置・印刷材の組み合わせ毎に適した閾値データを作成することを前提としているため、その準備及び作成処理に多くの作業工数を必要とし、改良の余地が十分にあった。
【０００８】
本発明は、上記した課題を解決するためになされたものであり、印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データを簡便且つ容易に決定可能な閾値データ設定装置、方法及びプログラム、並びに画像形成システムを提供することを目的とする。
【課題を解決するための手段】
【０００９】
本発明に係る閾値データ設定装置は、画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する装置であって、上記連続調画像信号に紐付けられた上記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得する閾値候補取得部と、上記２値画像信号の作成後から上記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得する応答特性取得部と、上記応答特性取得部により取得された上記応答特性を作用することで、上記網点形状の変化が模擬的に再現された模擬画像信号を作成する模擬画像作成部と、上記模擬画像作成部により作成された上記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出する定量値算出部と、上記閾値候補取得部により取得された上記各閾値候補を適用し、上記定量値算出部によりそれぞれ算出された上記定量値に基づく比較及び／又は評価によって、上記各閾値候補の中から設定対象である閾値データを決定する閾値データ決定部とを備える。
【００１０】
このように、２値画像信号の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得し、当該応答特性を作用することで、網点形状の変化が模擬的に再現された模擬画像信号を作成するので、網点形状を模擬画像上に簡便に且つ精度よく再現可能である。そして、各閾値候補として取得された複数の閾値データを、画質の程度を示す定量値に基づき比較及び／又は評価するので、印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データを簡便且つ容易に決定できる。
【００１１】
また、上記模擬画像作成部は、平網領域を少なくとも１つ含む上記模擬画像信号を作成し、上記定量値算出部は、上記模擬画像信号の上記平網領域でのパワースペクトルに関し、上記各網点に対応する成分を除いた主成分に基づいて上記定量値を算出することが好ましい。
【００１２】
また、上記模擬画像作成部は、階調レベルが異なる平網領域を２つ以上含む上記模擬画像信号を作成し、上記定量値算出部は、各上記平網領域における上記主成分についての統計値に基づいて上記定量値を算出することが好ましい。
【００１３】
また、上記画像形成過程にて上記網点形状に変化を及ぼす事象に応じた上記応答特性を推定する応答特性推定部を更に備え、上記応答特性取得部は、上記応答特性推定部により推定された上記応答特性を取得することが好ましい。
【００１４】
また、上記応答特性推定部は、刷版を介して被印刷材に色材を転移させる際の上記応答特性を推定することが好ましい。
【００１５】
また、上記応答特性推定部は、版材の主面を彫刻することで上記刷版を形成する際の上記応答特性を推定することが好ましい。
【００１６】
また、上記閾値データ決定部により決定された上記閾値データを設定することで、上記連続調画像信号を上記２値画像信号に変換する２値化処理部を更に備えることが好ましい。
【００１７】
本発明に係る画像形成システムは、上記の閾値データ設定装置と、上記閾値データ設定装置により変換された上記２値画像信号に基づいて上記画像を形成する画像形成装置とを備える。
【００１８】
本発明に係る閾値データ設定方法は、画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定する方法であって、上記連続調画像信号に紐付けられた上記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得するステップと、上記２値画像信号の作成後から上記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得するステップと、取得された上記応答特性を作用することで、上記網点形状の変化が模擬的に再現された模擬画像信号を作成するステップと、作成された上記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出するステップと、取得された上記各閾値候補を適用し、それぞれ算出された上記定量値に基づく比較及び／又は評価によって、上記各閾値候補の中から設定対象である閾値データを決定するステップとをコンピュータに実行させる。
【００１９】
本発明に係る閾値データ設定プログラムは、画像を構成する各網点の形状及び配置を表す２値画像信号を、連続調画像信号に基づき作成するための閾値データを設定するためのプログラムであって、上記連続調画像信号に紐付けられた上記各網点の形状及び配置に関する網点情報を取得し、該網点情報に合致する複数の閾値データを各閾値候補として取得するステップと、上記２値画像信号の作成後から上記画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得するステップと、取得された上記応答特性を作用することで、上記網点形状の変化が模擬的に再現された模擬画像信号を作成するステップと、作成された上記模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出するステップと、取得された上記各閾値候補を適用し、それぞれ算出された上記定量値に基づく比較及び／又は評価によって、上記各閾値候補の中から設定対象である閾値データを決定するステップとをコンピュータに実行させる。
【発明の効果】
【００２０】
本発明に係る閾値データ設定装置、方法及びプログラム、並びに画像形成システムによれば、２値画像信号の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性として取得し、当該応答特性を作用することで、網点形状の変化が模擬的に再現された模擬画像信号を作成するので、網点形状を模擬画像上に簡便に且つ精度よく再現可能である。そして、各閾値候補として取得された複数の閾値データを、画質の程度を示す定量値に基づき比較・評価するので、印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データを簡便且つ容易に決定できる。
【図面の簡単な説明】
【００２１】
【図１】この実施形態に係る画像形成システムの全体構成図である。
【図２】図１に示す２値化処理部の動作に関する概略説明図である。
【図３】網点情報に関する各パラメータの説明図である。
【図４】図１に示す凸版形成装置を構成するレーザ彫刻機の概略構成図である。
【図５】図１に示すフレキソ印刷機の概略構成図である。
【図６】図１に示す画像処理装置の電気的なブロック図である。
【図７】図１及び図６に示す画像処理装置の動作説明に供されるフローチャートである。
【図８】図８Ａ及び図８Ｂは、彫刻応答特性の一例を示すグラフである。
【図９】各テスト画像信号に対して２値化処理を施した結果を示す説明図である。
【図１０】図６に示す模擬画像作成部の動作に関する概略説明図である。
【図１１】観察距離が３００ｍｍであるドゥーリー・ショー関数のグラフである。
【図１２】図１２Ａは、模擬画像信号におけるパワースペクトルの一例を示すグラフである。図１２Ｂは、各閾値候補におけるパワースペクトルのピーク特性を示すグラフである。
【図１３】各閾値候補に対する定量結果を示すテーブルである。
【図１４】応答特性推定システムの全体構成図である。
【図１５】図１５Ａは、標本画像網点を可視化した概略正面図である。図１５Ｂは、図１５Ａの標本画像網点に基づき形成された標本凸版の正面図である。図１５Ｃは、図１５Ｂの標本凸版に基づき形成された標本印刷物の正面図である。
【図１６】図１６Ａは、網点突起部の高さ形状の推定モデルに関する第１の説明図である。図１６Ｂは、図１６Ａに示す高さ形状のフーリエ変換特性を示すグラフである。
【図１７】図１７Ａは、網点突起部の高さ形状の推定モデルに関する第２の説明図である。図１７Ｂは、図１７Ａに示す高さ形状のフーリエ変換特性を示すグラフである。
【発明を実施するための形態】
【００２２】
以下、本発明に係る閾値データ設定方法について、それを実施する閾値データ設定装置、プログラム及び画像形成システムとの関係において好適な実施形態を挙げ、添付の図面を参照しながら詳細に説明する。本明細書において、画像を形成することを「印刷」又は「印字」という場合がある。
【００２３】
［画像形成システム１０の全体構成］
図１は、この実施形態に係る画像形成システム１０の全体構成図である。この画像形成システム１０は、平板状の版材１２を彫刻することで凸版１４（刷版）を形成し、該凸版１４を介して被印刷体１６にインキ（色材）を転移させることで印刷物１８を生産するシステムである。この印刷物１８は、被印刷体１６上に少なくとも１色の網点が形成されてなる。
【００２４】
画像形成システム１０は、閾値データ設定装置としての画像処理装置２０と、凸版形成装置２２と、フレキソ印刷機２４と、閾値データベース（以下、閾値ＤＢ２６）と、応答特性データベース（以下、応答特性ＤＢ２８）を基本的に備えている。この実施形態にて、凸版形成装置２２及びフレキソ印刷機２４は、被印刷体１６上に画像を形成する画像形成装置を構成する。
【００２５】
画像処理装置２０は、入力された印刷用データ３０に基づいて、凸版１４の主面上の三次元形状を表す深さデータを作成し、当該データを凸版形成装置２２側に送出する。画像処理装置２０は、より詳細には、ＲＩＰ（Raster Image Processor）処理部３２、２値化処理部３４、閾値データ設定部３６及び深さデータ作成部３８を備えている。
【００２６】
ＲＩＰ処理部３２は、ページ記述言語（ＰＤＬ）で表現された連続調画像信号４０に対してラスタライズ処理を施すことで、ラスタ形式の画像信号（以下、ラスタ画像信号４２ともいう）に変換する。
【００２７】
２値化処理部３４は、ラスタ画像信号４２に対してスクリーニング処理を施すことで、２値の画像信号（以下、２値画像信号４６ともいう）に変換する。この２値画像信号４６は、画像を構成する各網点の形状及び配置を表している。この実施形態では、スクリーニング処理手法として、閾値データ４４による組織的ディザ法を採用する。
【００２８】
図２は、図１に示す２値化処理部３４の動作に関する概略説明図である。本図は、閾値データ４４（閾値マトリクス）を用いた２値化の概念を示す。先ず、ラスタ画像信号４２の各アドレスに閾値データ４４の各行列要素を対応付ける。そして、２値化処理部３４は、着目する画素での画素値Ｇと、対応する行列要素での閾値Ｔの大小関係をそれぞれ比較する。Ｇ＞Ｔを満たす場合には「１（オン）」を割り当て、それ以外（Ｇ≦Ｔ）の場合には「０（オフ）」を割り当てる。このようにして、階調レベルの数が「２６」であるラスタ画像信号４２は、階調レベルの数が「２」である２値画像信号４６に変換される。本図例では、２５個の画素を単位とする正方領域毎に、ハッチングが付された１３個の「オン」画素のクラスタ（ＡＭ網点）が形成されている。
【００２９】
図１に戻って、閾値データ設定部３６は、上記したスクリーニング処理に先立ち、閾値ＤＢ２６から印刷用データ３０に適した閾値データ４４を取得し、２値化処理部３４側に供給する。この印刷用データ３０には、上記した連続調画像信号４０の他、印刷物１８の形成に関する各種条件（以下、印刷条件４８）が含まれている。そして、印刷条件４８は、網点の形状・配置に関する情報（以下、網点情報５０）と、画像形成システム１０に関する装置又は印刷材に関する情報（以下、装置・印刷材情報５２）で構成される。
【００３０】
図３は、網点情報５０に関する各パラメータの説明図である。「スクリーン線数」は単位長さ当たりの網点の密度であり、単位としてＬＰＩ（Line Per Inch）が一般に用いられる。「スクリーン角度」は網点の配列方向と垂直軸（又は水平軸）とのなす角であり、単位として度（degree）が一般に用いられる。
【００３１】
図１に戻って、深さデータ作成部３８は、２値化処理部３４から取得した２値画像信号４６に基づいて、版材１２の彫刻量の２次元分布に相当する深さデータを作成する。深さデータの値又は範囲は、凸版形成装置２２に応じて任意に定義することができる。また、彫刻量に代わって、版材１２の基準面からの高さとして定義してもよい。
【００３２】
［凸版形成装置２２の概略構成］
凸版形成装置２２は、画像処理装置２０から供給された深さデータに基づき、版材１２の主面を彫刻することで、複数の網点突起部（図示しない）を備える凸版１４を形成可能なＣＴＰ（Computer To Plate）描画装置である。フレキソ印刷の場合は、版材１２として、ゴムシート、光硬化性樹脂シート等を含む高弾性材料が用いられる。
【００３３】
図４は、図１の凸版形成装置２２を構成するレーザ彫刻機６０の概略構成図である。レーザ彫刻機６０は、露光ヘッド６２と、ピント位置変更機構６４と、矢印Ｙ方向（副走査方向）への間欠送り機構６６とを備える。
【００３４】
ピント位置変更機構６４は、版材１２が取り付けられたドラム６８に対して、露光ヘッド６２を前後移動させるモータ７０及びボールネジ７２を有する。ピント位置変更機構６４は、モータ７０の制御によりピント位置を移動させる。
【００３５】
間欠送り機構６６は、ボールネジ７６及び該ボールネジ７６を回転させる副走査モータ７８を有する。間欠送り機構６６は、副走査モータ７８の制御により、露光ヘッド６２が搭載されたステージ７４を矢印Ｙ方向に移動させる。これにより、露光ヘッド６２は、ドラム６８の軸線８０の方向に間欠送りされる。
【００３６】
版材１２は、ドラム６８上でチャック部材８２によりチャックされる。チャック部材８２の位置は、露光ヘッド６２による露光が行われない領域に設定されている。この場合、軸線８０を中心としてドラム６８を回転させながら、該ドラム６８上の版材１２に対して、露光ヘッド６２からレーザビーム８４を照射することにより、版材１２の主面がレーザ彫刻される。そして、ドラム６８の回転により、露光ヘッド６２の前方をチャック部材８２が通過するときに、矢印Ｙ方向に間欠送りを行うことで、次のライン分のレーザ彫刻が行われる。
【００３７】
ドラム６８の回転による版材１２の矢印Ｘ方向（主走査方向）への送りと、露光ヘッド６２の矢印Ｙ方向（副走査方向）の間欠送りを逐次繰り返すことにより、露光走査の位置が制御される。この制御と同期して、画像処理装置２０（図１）から供給された位置毎の深さデータ（すなわち、彫刻量）に基づいて、レーザビーム８４のオン・オフ又は発光強度が制御される。この結果、版材１２の主面（版面）に所望の形状を有するレリーフ（網点突起部を含む）が形成される。
【００３８】
このようにして、凸版形成装置２２は、印刷に供される凸版１４を色版毎に形成する。この凸版１４は、フレキソ印刷機２４に装着される。
【００３９】
［フレキソ印刷機２４の概略構成］
フレキソ印刷機２４は、凸版形成装置２２により形成された凸版１４を介し、被印刷体１６の版面にインキを転移することで、印刷物１８を形成可能な印刷装置である。フレキソ印刷の場合は、被印刷体１６として、印刷用途の枚葉紙、ロール紙のみならず、表面の凹凸形状が大きい印刷媒体、例えば、段ボール紙、フイルム、布等を用いることができる。
【００４０】
図５は、図１に示すフレキソ印刷機２４の概略構成図である。フレキソ印刷機２４は、上記した凸版１４と、クッションテープ８６を介して凸版１４が取り付けられる版胴８８と、ドクターチャンバ９０によりインキが供給されるアニロックスローラ９２と、被印刷体１６を加圧するための圧胴９４とから構成される。
【００４１】
ドクターチャンバ９０から供給されたインキは、アニロックスローラ９２、凸版１４（網点突起部の頂面）の順に転移される。これに併せて、被印刷体１６は、版胴８８と圧胴９４との間に挟持された状態下に搬送される。すなわち、凸版１４に転移されたインキは被印刷体１６に更に転移されることで、多数の網点からなる画像が被印刷体１６上に形成される。
【００４２】
［画像処理装置２０の電気的なブロック図］
図６は、図１に示す画像処理装置２０の電気的なブロック図である。画像処理装置２０は、制御部１００と、通信Ｉ／Ｆ１０２と、入力部１０４と、表示部１０６と、出力Ｉ／Ｆ１０８と、メモリ１１０（記憶媒体）とを備えるコンピュータである。
【００４３】
通信Ｉ／Ｆ１０２は、外部装置からの電気信号を送受信するインターフェース（Ｉ／Ｆ）である。これにより、画像処理装置２０は、印刷用データ３０を含む各種データを外部装置から取得可能である。
【００４４】
入力部１０４は、マウス、トラックボール、キーボード、タッチパネル等の種々の入力デバイスで構成される。表示部１０６は、図示しない制御回路を介して、表示制御信号に応じた画像を表示する。入力部１０４による入力機能及び表示部１０６による表示機能を組み合わせることで、ユーザ・インターフェースを実現する。
【００４５】
出力Ｉ／Ｆ１０８は、電気信号としての深さデータを凸版形成装置２２に向けて送信するインターフェース（Ｉ／Ｆ）である。
【００４６】
メモリ１１０は、制御部１００が各構成要素を制御するのに必要なプログラム及びデータ等を記憶している。本図例では、印刷用データ３０、複数種類の閾値データ４４及び後述する応答特性データ１１２がそれぞれ格納されている。
【００４７】
メモリ１１０は、非一過性であり、且つ、コンピュータ読み取り可能な記憶媒体で構成されてもよい。ここで、コンピュータ読み取り可能な記憶媒体は、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ、フラッシュメモリ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置である。また、この記憶媒体は、短時間に且つ動的にプログラムを保持するものであっても、一定時間プログラムを保持するものであってもよい。
【００４８】
制御部１００は、ＣＰＵ（Central Processing Unit）等のプロセッサによって構成されている。制御部１００は、メモリ１１０に格納されたプログラムを読み出し実行することで、ＲＩＰ処理部３２、２値化処理部３４、閾値データ設定部３６、深さデータ作成部３８（図１）及び応答特性推定部１１４（図１４）の各機能を実現可能である。
【００４９】
閾値データ設定部３６は、設定対象の候補である閾値データ４４（以下、閾値候補）を取得する閾値候補取得部１１６、網点の形状変化を表現する空間周波数応答特性（以下、単に「応答特性」ともいう）を取得する応答特性取得部１１８、及び、各閾値候補における網点形状の再現性を定量評価する閾値候補評価部１２０を備える。
【００５０】
閾値候補評価部１２０は、２値網点作成部１２２、模擬画像作成部１２４、定量値算出部１２６及び閾値データ決定部１２８を備える。なお、各部の機能については後述する。
【００５１】
［画像処理装置２０の動作］
この実施形態に係る網点データ設定装置としての画像処理装置２０は、以上のように構成される。続いて、図１及び図６に示す画像処理装置２０（主に、閾値データ設定部３６）の動作について、図７のフローチャートを主に参照しながら詳細に説明する。
【００５２】
ステップＳ１において、閾値データ設定部３６は、連続調画像信号４０に紐付けられた印刷条件４８を取得する。図１に示すように、印刷条件４８には、網点情報５０と、装置・印刷材情報５２が含まれる。
【００５３】
網点情報５０として、例えば、スクリーン線数、スクリーン角度、網点の目標形状（円形、楕円形、矩形、これらの組合せ）、スクリーン特性（ＡＭ又はＦＭ）が挙げられる。装置・印刷材情報５２として、例えば、画像形成システム１０の装置構成、被印刷体１６、版材１２、色材、又は装置の構成部品（例えば、図４の露光ヘッド６２）の種類が挙げられる。
【００５４】
ステップＳ２において、閾値候補取得部１１６は、ステップＳ１にて得た印刷条件４８と部分的に一致する閾値データ４４（以下、閾値候補ともいう）を取得する。この取得に先立ち、閾値ＤＢ２６は、画像処理装置２０からの要求に応じてデータベースの検索処理を実行し、網点情報５０に合致する閾値データ４４を画像処理装置２０側に供給する。その後、閾値候補取得部１１６は、メモリ１１０に一時的に記憶された閾値候補をすべて読み出し取得する。ここでは、３つの閾値候補Ａ、Ｂ、Ｃをそれぞれ取得したとする。
【００５５】
ステップＳ３において、応答特性取得部１１８は、ステップＳ１にて得た印刷条件４８と部分的に一致する応答特性データ１１２を取得する。この取得に先立ち、応答特性ＤＢ２８は、画像処理装置２０からの要求に応じてデータベースの検索処理を実行し、装置・印刷材情報５２に合致する応答特性データ１１２を画像処理装置２０側に供給する。その後、応答特性取得部１１８は、メモリ１１０に一時的に記憶された応答特性データ１１２を読み出し取得する。
【００５６】
なお、応答特性データ１１２は、２値画像信号４６（図２）の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を示すデータであり、画像形成システム１０の構成に依存する点に留意する。図１例では、応答特性取得部１１８は、応答特性データ１１２として、レーザ彫刻機６０（図４）の光学的特性に相当する「彫刻応答特性」、及びフレキソ印刷機２４（図５）の転移特性に相当する「転移応答特性」をそれぞれ取得する。
【００５７】
図８Ａ及び図８Ｂは、彫刻応答特性の一例を示すグラフである。具体的には、図８Ａは点像分布関数｛以下、ＰＳＦ（Point Spread Function）｝であり、図８Ｂは図８Ａに示すＰＳＦのフーリエ変換関数の断面図に相当する。また、本グラフにおけるＸ軸、Ｙ軸は、図４に示す矢印Ｘ方向、矢印Ｙ方向にそれぞれ対応する。
【００５８】
図８Ａに示すように、ＰＳＦは、原点Ｏをピークとする概略円錐形状の関数である。ＰＳＦは、理想系ではデルタ関数であるが、実際系では種々の光学的因子により点像が拡がる。Ｘ軸方向の拡がり程度は、レーザビーム８４（図４）の矢印Ｘ方向の有効ビーム径、ドラム６８（同図）の回転速度、レーザ光源（不図示）の立ち上がり／立ち下がり特性等に依存する。また、Ｙ軸方向の拡がり程度は、レーザビーム８４の矢印Ｙ方向の有効ビーム径、レーザビーム８４の露光多重度等に依存する。
【００５９】
図８Ｂに示すＭＴＦ（Modulation Transfer Function）は、空間周波数毎の正弦波応答特性に相当する。本図から理解されるように、空間周波数が高くなるにつれてＭＴＦが次第に小さくなっている。そして、実線で示す「副走査」のＭＴＦは、破線で示す「主走査」よりも高くなる傾向を有する。
【００６０】
ステップＳ４において、閾値候補評価部１２０は、各閾値候補の定量評価に供されるテスト画像信号を作成する。テスト画像信号として任意のラスタ画像信号４２を用いてもよいが、ここでは、少なくとも１つの平網領域（ハーフトーンべた領域）を含む画像信号を作成しておく。
【００６１】
ステップＳ５において、閾値候補評価部１２０は、ステップＳ２により取得された複数の閾値候補（３つの閾値候補Ａ〜Ｃ）のうち、未だ評価されていない１つの閾値候補を指定する。以下、指定された閾値候補のことを「指定閾値候補」と称する。
【００６２】
ステップＳ６において、２値網点作成部１２２は、ステップＳ４で生成された各テスト画像信号に対して、指定閾値候補を用いたスクリーニング処理を施すことで、階調レベル毎の２値網点信号を作成する。
【００６３】
図９は、テスト画像信号に対して２値化処理を施した結果を示す説明図である。本図例では、階調レベルが１０％、３０％、５０％、７０％、９０％である５種類の平網領域を含むテスト画像信号から、濃淡が異なる５種類の平網画像を表す２値網点信号を作成する。
【００６４】
ステップＳ７において、模擬画像作成部１２４は、ステップＳ６で作成された２値網点信号から、印刷物１８上での網点形状が模擬的に再現された模擬画像信号を作成する。
【００６５】
図１０は、図６に示す模擬画像作成部１２４の動作に関する概略説明図である。この図は、２値網点信号から模擬画像信号を作成するまでの処理過程を模式的に示している。
【００６６】
先ず、模擬画像作成部１２４は、各画素のオン・オフ状態を示す２値網点信号を、凸版１４の高さ分布に相当する多値網点信号に変換する。その後、模擬画像作成部１２４は、得られた多値網点信号に対して空間周波数領域での応答特性を作用する処理を実行する。この作用処理には、マスク処理、フィルタ処理、又はこれらの組合せを適用できる。なぜならば、フーリエ変換の写像の性質上、ユークリッド空間上又はフーリエ空間上での演算結果は実質的に等価だからである。
【００６７】
ここで、「マスク処理」は、ユークリッド空間上においてマスク演算（畳み込み）を実行することで、画像の空間周波数成分（スペクトル強度）を変調する画像処理を意味する。また、「フィルタ処理」は、フーリエ空間上においてフィルタ演算（乗算）を実行することで、画像信号の空間周波数成分（スペクトル強度）を変調する画像処理を意味する。
【００６８】
本図上段に示す処理流れによれば、模擬画像作成部１２４は、彫刻応答特性を示すマスクを用いた第１マスク処理、転移応答特性を示すマスクを用いた第２マスク処理を順次施すことで模擬画像信号を作成する。また、本図下段の処理流れによれば、模擬画像作成部１２４は、彫刻応答特性を示すフィルタを用いた第１フィルタ処理、転移応答特性を示すフィルタを用いた第２フィルタ処理を順次施すことで模擬画像信号を作成する。ここで、フィルタ処理の実行前にＦＦＴ（Fast Fourier Transform）を、フィルタ処理の実行後にＩＦＦＴ（Inverse Fast Fourier Transform）をそれぞれ施す点に留意する。
【００６９】
ステップＳ８において、定量値算出部１２６は、ステップＳ７で作成された模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出する。この実施形態では、定量値を算出するために、模擬画像信号のパワースペクトルを利用する。例えば、視認性との相関度を高めるため、このパワースペクトルに対して人間の標準視覚応答特性を乗算する前処理を行ってもよい。
【００７０】
図１１は、観察距離が３００ｍｍであるドゥーリー・ショー（Dooley-Shaw）関数のグラフである。この関数は、ＶＴＦ（Visual Transfer Function）の一種であり、人間の標準視覚応答特性を模した代表的な関数である。具体的には、輝度のコントラスト比特性の２乗値に相当する。グラフの横軸は空間周波数（単位：Ｃｙｃｌｅ／ｍｍ）であり、縦軸はＶＴＦの値（単位は無次元）である。
【００７１】
なお、ＶＴＦの関数形状はこれに限られず、数理モデルや実験データ等から導出された種々の特性を適用してもよい。また、観察距離は３００ｍｍ相当のみならず、画像の観察態様又は評価基準等に応じて種々変更してもよい。
【００７２】
また、上記の前処理による視覚的な補正効果を適切に反映させるため、模擬画像信号の画素値を、画像を反射又は透過する光量に対して相関性が高い（好ましくは線形的な）量に変換しておいてもよい。その一例として、ＲＧＢ値、三刺激値（ＸＹＺ）、反射原稿の場合は光反射率、透過原稿の場合は光透過率等を用いてもよい。
【００７３】
図１２Ａは、模擬画像信号におけるパワースペクトルの一例を示すグラフである。グラフの横軸は空間周波数（単位：Ｃｙｃｌｅ／ｍｍ）であり、グラフの縦軸はスペクトル強度（単位は任意；以下、単に「成分」ともいう）である。本図に示すパワースペクトルは、空間周波数が約１．９、約２．４、約３．６、約７．４Ｃｙｃｌｅ／ｍｍの位置に４つのピークを有する。定量値算出部１２６は、各網点に対応する成分（ここでは、約７．４Ｃｙｃｌｅ／ｍｍ周辺の成分）を除いた上で、スペクトル強度の最大値（主成分に相当；以下、最大強度という）を抽出する。この場合、定量値算出部１２６は、約２．４Ｃｙｃｌｅ／ｍｍにおけるピーク最大値である「６０」を、最大強度として抽出しておく。この最大強度の値が大きいほど画質の程度が低く、値が小さいほど画質の程度が高いと評価する。この定量値は、印刷物１８の画質に対して定量的又は定性的に相関する値である。なぜならば、各網点に対応する成分以外の成分は、色ムラ、モアレ、粒状感等を含む画像ノイズとして寄与するからである。
【００７４】
なお、本図例のグラフは、閾値候補Ｂを用いて、階調レベルが７０％であるテスト画像信号から作成された模擬画像信号のパワースペクトルに相当する。定量値算出部１２６は、７０％の場合も含め、階調レベル毎の模擬画像信号に対して最大強度をそれぞれ抽出する。
【００７５】
その後、定量値算出部１２６は、上記のように得られた階調レベル毎の最大強度に基づいて、指定閾値候補の定量値を算出する。ここでは、定量値算出部１２６は、全階調レベルにわたる平均値及び／又は最大値を「定量値」として算出する。なお、定量値の算出式はこれに限られず、各平網領域における最大強度についての統計値、又はこれらの組み合わせ（例えば、平均値と最大値の重み付け加算値）に基づいて種々算出してもよい。また、模擬画像信号に対する評価処理は上記した手法に限られず、粒状性｛ＲＭＳ（Root Mean Square）等｝、鮮鋭性（ＭＴＦ等）、階調性（階調滑らかさ等）を含む公知の画像評価手法を種々適用できることは言うまでもない。
【００７６】
ステップＳ９において、閾値候補評価部１２０は、すべての閾値候補に対する指定が完了したか否かを判定する。まだ完了していないと判定された場合（ステップＳ９：ＮＯ）、ステップＳ５に戻って、以下ステップＳ５〜Ｓ８を順次繰り返す。一方、すべて完了したと判定された場合（ステップＳ９：ＹＥＳ）、次のステップ（Ｓ１０）に進む。
【００７７】
ステップＳ１０において、閾値データ決定部１２８は、閾値候補の比較・評価を行うことで、１つの閾値データ４４を決定する。
【００７８】
図１２Ｂは、各閾値候補Ａ〜Ｃにおける階調レベル毎の最大強度を示すグラフである。図１３は、各閾値候補Ａ〜Ｃに対する定量値の算出結果を示すテーブルである。閾値候補評価部１２０は、図１３に示す評価結果を踏まえ、「平均値」或いは「最大値」が最も小さい閾値候補Ａを、閾値データ４４として決定する。
【００７９】
ステップＳ１１において、閾値データ設定部３６は、ステップＳ１０で決定された閾値データ４４を印刷用データ３０に対応付ける設定を行う。これにより、印刷用データ３０から深さデータを作成する際に、画像形成システム１０に適した閾値データ４４を自動的に設定することができる。このようにして、画像処理装置２０（主に閾値データ設定部３６）の動作は完了する。
【００８０】
［応答特性推定システム１５０の構成及び動作］
続いて、応答特性推定部１１４（図６）を組み込んだ応答特性推定システム１５０の構成及び動作について、図１４〜図１７Ｂを参照しながら説明する。
【００８１】
図１４は、応答特性推定システム１５０の全体構成図である。この応答特性推定システム１５０は、標本網点信号１５２に基づき標本凸版１５４及び標本印刷物１５６を順次形成し、第１測定器１５８、第２測定器１６０による測定結果を用いて、画像形成システム１０の全体又は一部における応答特性を推定するシステムである。
【００８２】
標本網点信号１５２は、任意の画像を構成する各網点の形状・配置を表す２値の画像信号である。ここでは、後述する推定処理を容易にするため、標本網点信号１５２は、２値網点信号（図９）と同一の信号であってもよい。
【００８３】
深さデータ作成部３８は、入力された標本網点信号１５２に基づいて、版材１２の彫刻量の２次元分布に相当する深さデータを作成する。そして、凸版形成装置２２は、深さデータ作成部３８から供給された深さデータに基づき、版材１２の主面を彫刻することで、複数の網点突起部（図示しない）を備える標本凸版１５４を形成する。そして、フレキソ印刷機２４は、凸版形成装置２２により形成された標本凸版１５４を介し、被印刷体１６の版面にインキを転移することで、標本印刷物１５６を形成する。
【００８４】
第１測定器１５８は、標本凸版１５４の凹凸形状を光学的に読み取る装置であり、標本凸版１５４上での網点面積率を測定可能である。また、第２測定器１６０は、標本印刷物１５６の色特性を光学的に読み取る装置であり、標本印刷物１５６上での網点面積率を測定可能である。
【００８５】
応答特性推定部１１４は、画像形成過程にて網点形状に変化を及ぼす事象に応じた応答特性を推定する。網点形状に変化を及ぼす事象として、画像形成装置の記録方式に応じた１つ又は２つ以上の事象を考慮してもよい。ここでは、応答特性推定部１１４は、版材１２の主面を彫刻することで凸版１４を形成する際の「彫刻応答特性」、凸版１４を介して被印刷体１６にインキを転移させる際の「転移応答特性」をそれぞれ推定する。
【００８６】
応答特性推定部１１４は、第１測定器１５８により測定された標本凸版１５４上での網点面積率に基づき彫刻応答特性を推定する彫刻特性推定部１６２と、第２測定器１６０により測定された標本印刷物１５６上での網点面積率に基づき転移応答特性を推定する転移特性推定部１６４とを備える。
【００８７】
以下、彫刻特性推定部１６２による推定方法について、図１５Ａ〜図１６Ｂを参照しながら詳細に説明する。
【００８８】
図１５Ａは、図１４に示す標本網点信号１５２を可視化した概略正面図である。ここで、図１５Ａ及び後述する図１５Ｂ、図１５Ｃにおける塗り潰し領域は、各網点形状の決定に寄与する部位を示す。本図例では、各網点が円形状であり、直径がＤ（既知の値）である場合を想定する。
【００８９】
図１５Ｂは、図１５Ａの標本網点信号１５２に基づき形成された標本凸版１５４の概略正面図である。図８Ａ及び図８Ｂに示す応答特性の作用下、標本凸版１５４上での網点面積率は、網点面積率の設計値よりも小さくなる。その結果、測定された網点面積率に相当する平均直径が（Ｄ−ｄ１）であったとする。
【００９０】
図１５Ｃは、図１５Ｂの標本凸版１５４に基づき形成された標本印刷物１５６の概略正面図である。版材１２として高弾性材料が用いられる場合、インキを転移する際にドットゲインが増加する傾向がある。このとき、標本印刷物１５６上での網点面積率は、標本凸版１５４上での網点面積率よりも大きくなる。その結果、測定された網点面積率に相当する平均直径が（Ｄ−ｄ２）（ただし、ｄ２＜ｄ１）であったとする。
【００９１】
ここで、測定された網点面積率から網点突起部の高さ形状を推定し、該高さ形状の変化に基づいて彫刻応答特性を推定する手法を適用できる。以下、具体的な方法について、図１６Ａ〜図１７Ｂを参照しながら説明する。
【００９２】
図１６Ａは、網点突起部の高さ形状の推定モデルに関する第１の説明図である。グラフの横軸はＸ軸（又はＹ軸）方向の位置であり、グラフの縦軸はＺ軸（高さ）方向の位置である。本図の実線は、標本網点信号１５２における各網点突起部の理想的な高さ形状（以下、理想網点形状１７０）を示している。本図の破線は、標本凸版１５４における各網点突起部の仮想的な高さ形状（以下、仮想網点形状１７２）を示している。この推定モデルは、網点突起部の高さが常に一定である前提の下に構築されており、主にドットゲインの現象を模擬している。
【００９３】
図１６Ｂは、図１６Ａに示す高さ形状のフーリエ変換特性を示すグラフである。グラフの横軸は空間周波数であり、グラフの縦軸はフーリエ変換値である。本図の実線は、理想網点形状１７０に対するフーリエ変換値（以下、理想周波数成分１７０ｆ）を示している。本図の破線は、仮想網点形状１７２に対するフーリエ変換値（以下、仮想周波数成分１７２ｆ）を示している。
【００９４】
この場合、彫刻特性推定部１６２は、理想周波数成分１７０ｆに対する仮想周波数成分１７２ｆの比、すなわち、空間周波数毎の（仮想周波数成分１７２ｆ／理想周波数成分１７０ｆ）の値を、彫刻応答特性として推定する。
【００９５】
図１７Ａは、網点突起部の高さ形状の推定モデルに関する第２の説明図である。グラフの定義は図１６Ａと同様であるため、その説明を省略する。本図の実線は理想網点形状１７０を示すと共に、本図の破線は仮想網点形状１７４を示している。
【００９６】
この推定モデルは、頂部の直径にかかわらず、網点突起部の体積（面積）が常に保存される前提の下に構築されている。本図から理解されるように、頂部の直径が小さくなるにつれて側面の傾斜が小さくなっている。
【００９７】
図１７Ｂは、図１７Ａに示す高さ形状のフーリエ変換特性を示すグラフである。グラフの定義は図１６Ｂと同様であるため、その説明を省略する。本図の実線は理想周波数成分１７０ｆを示すと共に、本図の破線は仮想周波数成分１７４ｆを示している。
【００９８】
この場合、彫刻特性推定部１６２は、理想周波数成分１７０ｆに対する仮想周波数成分１７４ｆの比、すなわち、空間周波数毎の（仮想周波数成分１７４ｆ／理想周波数成分１７０ｆ）の値を、彫刻応答特性として推定する。
【００９９】
図１４に戻って、転移特性推定部１６４は、応答特性の乗算則（カスケードモデル）に基づいて転移応答特性を推定する。先ず、転移特性推定部１６４は、図１６Ａ〜図１７Ｂと同一の又は異なる推定モデルを用いて、画像形成システム１０全体における応答特性（以下、全体応答特性）を推定する。そして、転移特性推定部１６４は、全体応答特性から彫刻応答特性を除算することで、転移応答特性を推定する。
【０１００】
最後に、応答特性推定部１１４は、それぞれ推定された彫刻応答特性及び転移応答特性を応答特性ＤＢ２８側に供給する。そして、応答特性ＤＢ２８は、取得した彫刻応答特性を、版材１２及び凸版形成装置２２の種類に対応付ける。応答特性ＤＢ２８は、取得した転移応答特性を、版材１２及びフレキソ印刷機２４の種類に対応付ける。このようにして、応答特性データ１１２は、装置・印刷材情報５２に紐付けられた状態下、応答特性ＤＢ２８内にて保存・管理される。
【０１０１】
［本発明の効果］
閾値データ設定装置としての画像処理装置２０は、画像を構成する各網点の形状・配置を表す２値画像信号４６を、連続調画像信号４０に基づき作成するための閾値データ４４を設定する装置である。
【０１０２】
そして、連続調画像信号４０に紐付けられた各網点の形状・配置に関する網点情報５０を取得し、網点情報５０に合致する複数の閾値データ４４を各閾値候補として取得する閾値候補取得部１１６と、２値画像信号４６の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性データ１１２として取得する応答特性取得部１１８と、応答特性データ１１２を作用することで、網点形状の変化が模擬的に再現された模擬画像信号を作成する模擬画像作成部１２４と、模擬画像信号に対して所定の評価処理を実行して得られる評価値を、画質の程度を示す定量値として算出する定量値算出部１２６と、各閾値候補を適用しそれぞれ算出された定量値に基づく比較・評価によって、各閾値候補の中から設定対象である閾値データ４４を決定する閾値データ決定部１２８とを備える。
【０１０３】
このように、２値画像信号４６の作成後から画像を形成するまでの画像形成過程にて生じる網点形状の変化を、空間周波数領域での応答特性データ１１２として取得し、当該応答特性データ１１２を作用することで、網点形状の変化が模擬的に再現された模擬画像信号を作成するので、網点形状を模擬画像上に簡便に且つ精度よく再現可能である。そして、各閾値候補として取得された複数の閾値データ４４を、画質の程度を示す定量値に基づき比較・評価するので、印刷装置及び印刷材の種類数が膨大な場合であっても、これらの組み合わせに適した閾値データ４４を簡便且つ容易に決定できる。
【０１０４】
［補足］
なお、この発明は、上述した実施形態に限定されるものではなく、この発明の主旨を逸脱しない範囲で自由に変更できることは勿論である。
【０１０５】
この実施形態ではフレキソ印刷を例に説明したが、画像形成装置はこの印刷方式に限られず、オフセット印刷、グラビア印刷、デジタル印刷機（オンデマンド印刷機）等にも適用できることは言うまでもない。
【符号の説明】
【０１０６】
１０…画像形成システム                    １４…凸版
１６…被印刷体                            １８…印刷物
２０…画像処理装置                        ２２…凸版形成装置
２４…フレキソ印刷機                      ３０…印刷用データ
３４…２値化処理部                        ３６…閾値データ設定部
４０…連続調画像信号                      ４４…閾値データ
５０…網点情報                            ５２…装置・印刷材情報
６０…レーザ彫刻機                        １００…制御部
１１２…応答特性データ                    １１４…応答特性推定部
１１６…閾値候補取得部                    １１８…応答特性取得部
１２０…閾値候補評価部                    １２２…２値網点作成部
１２４…模擬画像作成部                    １２６…定量値算出部
１２８…閾値データ決定部                  １５０…応答特性推定システム
１５２…標本網点信号                      １５４…標本凸版
１５６…標本印刷物                        １５８…第１測定器
１６０…第２測定器                        １７０…理想網点形状
１７２、１７４…仮想網点形状
【図１】
0005977716.tif 000002
【図２】
0005977716.tif 000003
【図３】
0005977716.tif 000004
【図４】
0005977716.tif 000005
【図５】
0005977716.tif 000006
【図６】
0005977716.tif 000007
【図７】
0005977716.tif 000008
【図８】
0005977716.tif 000009
【図９】
0005977716.tif 000010
【図１０】
0005977716.tif 000011
【図１１】
0005977716.tif 000012
【図１２】
0005977716.tif 000013
【図１３】
0005977716.tif 000014
【図１４】
0005977716.tif 000015
【図１５】
0005977716.tif 000016
【図１６】
0005977716.tif 000017
【図１７】
0005977716.tif 000018