(57)【要約】【課題】筋むらに対する頑健性が高い閾値マトリクス作成装置、方法、プログラム及び画像形成装置を提供する。【解決手段】複数のドットが形成される位置誤差Ｅｒｐ１、Ｅｒｐ２又はサイズ誤差Ｅｒｓ１、Ｅｒｓ２に基づいて、閾値マトリクスＭｔの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した第１評価値ＥＶ１を算出する。そして、第１評価値ＥＶ１に基づいて閾値マトリクスＭｔを決定する。【選択図】図７
2014096611.tif 000002
【特許請求の範囲】
【請求項１】
複数のドットを形成可能である画像形成装置に応じたハーフトーン処理のための閾値マトリクスを作成する装置であって、
前記閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した第１評価値に基づいて前記閾値マトリクスを決定する閾値マトリクス決定部を有し、
前記第１評価値は、前記複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて算出された評価値である
ことを特徴とする閾値マトリクス作成装置。
【請求項２】
請求項１記載の装置において、
前記閾値マトリクス決定部は、
前記画像形成装置に入力される画像信号の複数の階調レベルに応じた複数のドットパターンを前記第１評価値に基づいて作成するドットパターン作成部と、
前記ドットパターン作成部により作成された前記複数のドットパターンにおけるドットの配列に従って前記閾値マトリクスを決定する閾値決定部と
を備える
ことを特徴とする閾値マトリクス作成装置。
【請求項３】
請求項２記載の装置において、
前記閾値マトリクス決定部は、幾何学的計算により得られた前記重なり程度から前記第１評価値を算出する評価値算出部をさらに備えることを特徴とする閾値マトリクス作成装置。
【請求項４】
請求項２記載の装置において、
前記閾値マトリクス決定部は、
各前記ドットパターンよりも小さいサイズのサブパターンを複数決定するサブパターン決定部と、
各前記ドットパターンの部分パターンと、前記サブパターン決定部により決定された各前記サブパターンとの一致度に基づいて前記第１評価値を算出する評価値算出部と
をさらに備える
ことを特徴とする閾値マトリクス作成装置。
【請求項５】
請求項４記載の装置において、
前記閾値マトリクス決定部は、
前記位置誤差の統計量又は前記サイズ誤差の統計量に基づいて各前記サブパターンにおける前記重なり程度の頑健性を定量化し、各得点を決定する頑健性得点決定部をさらに備え、
前記評価値算出部は、前記頑健性得点決定部により決定された各前記サブパターンに応じた前記得点に基づいて前記第１評価値を算出する
ことを特徴とする閾値マトリクス作成装置。
【請求項６】
請求項３〜５のいずれか１項に記載の装置において、
前記評価値算出部は、各前記ドットパターンのノイズ粒状性を定量化した第２評価値をさらに算出し、
前記ドットパターン作成部は、前記評価値算出部により算出された前記第１評価値及び前記第２評価値に基づいて前記複数のドットパターンを作成する
ことを特徴とする閾値マトリクス作成装置。
【請求項７】
請求項６記載の装置において、
前記第２評価値は、人間の視覚応答特性関数により補正されたＲＭＳ粒状度であることを特徴とする閾値マトリクス作成装置。
【請求項８】
請求項１〜７のいずれか１項に記載の装置において、
前記統計量には、標準偏差、平均値、最頻値、中心値、最大値及び最小値のうち少なくとも１つが含まれることを特徴とする閾値マトリクス作成装置。
【請求項９】
請求項１〜８のいずれか１項に記載の装置において、
前記第１評価値は、前記行方向及び前記列方向のうち誤差が大きい方の前記位置誤差の統計量に基づいて算出された前記一方向における評価値であることを特徴とする閾値マトリクス作成装置。
【請求項１０】
複数のドットを形成可能である画像形成装置に応じたハーフトーン処理のための閾値マトリクスを作成する方法であって、
前記複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて、前記閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した評価値を算出する算出ステップと、
算出された前記評価値に基づいて前記閾値マトリクスを決定する決定ステップと
を備えることを特徴とする閾値マトリクス作成方法。
【請求項１１】
請求項１０記載の閾値マトリクス作成方法において、
前記画像形成装置が、配列方向に沿って配列された複数のドット形成素子を備え、該複数のドット形成素子を用いて複数のドットを記録媒体上に形成するドット形成部と、前記ドット形成部及び前記記録媒体のうち少なくとも一方を所定の搬送方向に搬送することで前記ドット形成部と前記記録媒体とを相対移動させる搬送部とを有する場合、
前記複数のドット形成素子のうち１つのドット形成素子を用いて、前記搬送方向に延在する線画像を形成する画像形成ステップと、
形成された前記線画像から、前記１つのドット形成素子における前記位置誤差の統計量又は前記サイズ誤差の統計量を計測する計測ステップと
をさらに備えることを特徴とする閾値マトリクス作成方法。
【請求項１２】
複数のドットを形成可能である画像形成装置に応じたハーフトーン処理のための閾値マトリクスを作成するプログラムであって、
コンピュータに、前記閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した評価値に基づいて前記閾値マトリクスを決定するマトリクス決定部として機能させ、
前記第１評価値は、前記複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて算出された評価値である
ことを特徴とするプログラム。
【請求項１３】
配列方向に沿って配列された複数のドット形成素子を備え、該複数のドット形成素子を用いて複数のドットを記録媒体上に形成するドット形成部と、
前記ドット形成部及び前記記録媒体のうち少なくとも一方を所定の搬送方向に搬送することで前記ドット形成部と前記記録媒体とを相対移動させる搬送部と、
前記搬送部による前記相対移動の下、ハーフトーン処理に供される閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した評価値に基づいて前記閾値マトリクスを決定するマトリクス決定部と、
前記マトリクス決定部により決定された前記閾値マトリクスを用いた前記ハーフトーン処理により、入力された画像信号を、前記ドット形成部を制御する制御信号に変換する信号変換部と
を有し、
前記評価値は、前記複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて算出された評価値である
ことを特徴とする画像形成装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、ハーフトーン処理に供される閾値マトリクスを作成する閾値マトリクス作成装置、方法、プログラム及び画像形成装置に関する。
【背景技術】
【０００２】
近時、インクジェット技術の飛躍的進歩に伴い、インクジェット方式の画像形成装置による高速・高画質を両立したカラー大判印刷が可能になりつつある。この装置は、特にサイン・ディスプレイ用途において幅広い分野で用いられ、例えば、店頭ＰＯＰ（Point Of Purchase）や壁面ポスター、屋外広告・看板等の印刷にも適用可能である。インクジェット方式では、記録媒体上に複数種のインク（例えばＣＭＹＫインク）の液滴を吐出して多数のドットを形成することで、印刷物を得ることができる。
【０００３】
そして、前記インクジェット方式のうち、記録媒体の幅方向（以下、主方向という。）に沿って配列された複数のノズルを備える記録ヘッド（以下、ラインヘッドという。）を使用するシングルパス方式が特に注目されている。なぜならば、所定の直線パスに沿って記録媒体を１回だけ相対的に移動させることで記録媒体上の画像を完成可能であり、サイン・ディスプレイ用途で要求される各種仕様（高速化・低電力化・高画質化）をすべて両立し得るからである。一方、シングルパス方式では、何らかの変動要因により、各ノズルからのインク滴の着弾位置が目標位置から主方向にずれる結果、記録媒体の搬送方向に延在する濃い又は淡い筋むら（以下、単に筋むらという。）が発生し易いという課題があった。
【０００４】
この課題に対して、シングルパス方式であっても、前記筋むらの発生程度を低減可能なラインヘッド等の改良技術が種々提案されている。
【０００５】
特許文献１には、斜めの列方向に沿って複数のノズルが規則的な配列間隔でマトリクス状に多数配置された構造を有するラインヘッドが開示されている。これにより、複数のタイミングで各ドットを順次形成させて画像を完成することで、主方向での繰り返し周期を短縮可能であり、筋むらの発生を抑制できる。
【０００６】
このほか、各ノズルからインク滴を吐出制御するための制御信号を補正することで、筋むらの発生程度を低減する方法も種々提案されている。
【０００７】
特許文献２には、各ノズルの飛行曲がり量から主方向に隣接するドットの隙間を算出し、該隙間の大小に応じてハーフトーン処理に供される閾値マトリクス（ディザマトリクス）を補正する装置及び方法が開示されている。
【０００８】
特許文献３には、インク滴の着弾位置のずれ量を検出可能なテストチャートを印刷した後、その画像を読取処理部で読み取り、読取結果に基づいて画像信号を補正する装置及び方法が開示されている。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特開２００７−４４９６７号公報
【特許文献２】特開２００６−１５９８１０号公報
【特許文献３】特開２００１−１１３８０５号公報
【発明の概要】
【発明が解決しようとする課題】
【００１０】
ところで、各ノズルからのインク滴の着弾位置がずれる変動要因には、各ノズルの成形加工精度による恒久的な要因もあれば、ノズル付近での塵埃の付着等による一時的な要因もある。後者は、ノズル状態を回復させるべく、吸引等によるヘッドクリーニング処理を施すことで、インク滴の飛行曲がり量を低減できる場合がある。
【００１１】
しかし、上記した一時的な変動要因により、飛行曲がりの方向が突発的に反転する場合、特許文献２及び３に開示された装置等では、期待された補正の効果が得られないどころか、筋むらの発生が却って増長される懸念もある。すなわち、筋むらに対する頑健性（ロバストネス）が低いと言える。
【００１２】
本発明は上記した問題を解決するためになされたもので、筋むらに対する頑健性が高い閾値マトリクスを作成する閾値マトリクス作成装置、方法、プログラム及び画像形成装置を提供することを目的とする。
【課題を解決するための手段】
【００１３】
本発明は、複数のドットを形成可能である画像形成装置に応じたハーフトーン処理のための閾値マトリクスを作成する閾値マトリクス作成装置に関する。
【００１４】
前記閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した第１評価値に基づいて前記閾値マトリクスを決定する閾値マトリクス決定部を有し、前記第１評価値は、前記複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて算出された評価値であることを特徴とする。
【００１５】
このように、複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて、閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した第１評価値を算出するようにしたので、ドットの形成ばらつきに起因する画像への影響度を統計的且つ定量的に予測可能である。そして、前記第１評価値に基づいて閾値マトリクスを適切に決定することで、前記一方向に延在する筋むらに対する頑健性が高くなる。
【００１６】
また、前記閾値マトリクス決定部は、前記画像形成装置に入力される画像信号の複数の階調レベルに応じた複数のドットパターンを前記第１評価値に基づいて作成するドットパターン作成部と、前記ドットパターン作成部により作成された前記複数のドットパターンにおけるドットの配列に従って前記閾値マトリクスを決定する閾値決定部とを備えることが好ましい。
【００１７】
さらに、前記閾値マトリクス決定部は、幾何学的計算により得られた前記重なり程度から前記第１評価値を算出する評価値算出部をさらに備えることが好ましい。
【００１８】
また、前記閾値マトリクス決定部は、各前記ドットパターンよりも小さいサイズのサブパターンを複数決定するサブパターン決定部と、各前記ドットパターンの部分パターンと、前記サブパターン決定部により決定された各前記サブパターンとの一致度に基づいて前記第１評価値を算出する評価値算出部とをさらに備えることが好ましい。
【００１９】
さらに、前記閾値マトリクス決定部は、前記位置誤差の統計量又は前記サイズ誤差の統計量に基づいて各前記サブパターンにおける前記重なり程度の頑健性を定量化し、各得点を決定する頑健性得点決定部をさらに備え、前記評価値算出部は、前記頑健性得点決定部により決定された各前記サブパターンに応じた前記得点に基づいて前記第１評価値を算出することが好ましい。
【００２０】
さらに、前記評価値算出部は、各前記ドットパターンのノイズ粒状性を定量化した第２評価値をさらに算出し、前記ドットパターン作成部は、前記評価値算出部により算出された前記第１評価値及び前記第２評価値に基づいて前記複数のドットパターンを作成することが好ましい。
【００２１】
さらに、前記第２評価値は、人間の視覚応答特性関数により補正されたＲＭＳ粒状度であることが好ましい。
【００２２】
さらに、前記統計量には、標準偏差、平均値、最頻値、中心値、最大値及び最小値のうち少なくとも１つが含まれることが好ましい。
【００２３】
さらに、前記第１評価値は、前記行方向及び前記列方向のうち誤差が大きい方の前記位置誤差の統計量に基づいて算出された前記一方向における評価値であることが好ましい。
【００２４】
本発明に係る閾値マトリクス作成方法は、複数のドットを形成可能である画像形成装置に応じたハーフトーン処理のための閾値マトリクスを作成する方法であって、前記複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて、前記閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した評価値を算出する算出ステップと、算出された前記評価値に基づいて前記閾値マトリクスを決定する決定ステップとを備えることを特徴とする。
【００２５】
また、前記画像形成装置が、配列方向に沿って配列された複数のドット形成素子を備え、該複数のドット形成素子を用いて複数のドットを記録媒体上に形成するドット形成部と、前記ドット形成部及び前記記録媒体のうち少なくとも一方を所定の搬送方向に搬送することで前記ドット形成部と前記記録媒体とを相対移動させる搬送部とを有する場合、前記複数のドット形成素子のうち１つのドット形成素子を用いて、前記搬送方向に延在する線画像を形成する画像形成ステップと、形成された前記線画像から、前記１つのドット形成素子における前記位置誤差の統計量又は前記サイズ誤差の統計量を計測する計測ステップとをさらに備えることが好ましい。
【００２６】
本発明に係るプログラムは、複数のドットを形成可能である画像形成装置に応じたハーフトーン処理のための閾値マトリクスを作成するプログラムであって、コンピュータに、前記閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した評価値に基づいて前記閾値マトリクスを決定するマトリクス決定部として機能させ、前記第１評価値は、前記複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて算出された評価値であることを特徴とする。
【００２７】
本発明に係る画像形成装置は、配列方向に沿って配列された複数のドット形成素子を備え、該複数のドット形成素子を用いて複数のドットを記録媒体上に形成するドット形成部と、前記ドット形成部及び前記記録媒体のうち少なくとも一方を所定の搬送方向に搬送することで前記ドット形成部と前記記録媒体とを相対移動させる搬送部と、前記搬送部による前記相対移動の下、ハーフトーン処理に供される閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した評価値に基づいて前記閾値マトリクスを決定するマトリクス決定部と、前記マトリクス決定部により決定された前記閾値マトリクスを用いた前記ハーフトーン処理により、入力された画像信号を、前記ドット形成部を制御する制御信号に変換する信号変換部とを有し、前記評価値は、前記複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて算出された評価値であることを特徴とする。
【発明の効果】
【００２８】
本発明に係る閾値マトリクス作成装置、方法、プログラム及び画像形成装置によれば、複数のドットが形成される位置誤差の統計量又はサイズ誤差の統計量に基づいて、閾値マトリクスの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した評価値を算出するようにしたので、ドットの形成ばらつきに起因する画像への影響度を統計的且つ定量的に予測可能である。そして、前記評価値に基づいて閾値マトリクスを適切に決定することで、前記一方向に延在する筋むらに対する頑健性が高くなる。
【図面の簡単な説明】
【００２９】
【図１】本実施の形態に係る閾値マトリクス作成装置としての画像処理装置が組み込まれた画像形成システムの概略説明図である。
【図２】図１に示す画像処理装置の電気的なブロック図である。
【図３】図１に示す信号変換部における信号処理流れを表す概略説明図である。
【図４】組織的ディザ法によるハーフトーン処理の概略説明図である。
【図５】図１に示す画像形成装置の構成を表す断面側面図である。
【図６】図５に示すラインヘッドの構成例を表す平面透視図である。
【図７】図７Ａは、ドットパターンの一部を表す概略説明図である。図７Ｂは、第１評価値の算出方法を表す概略説明図である。
【図８】図８Ａは、ドットパターンを表す画像データを可視化した概略説明図である。図８Ｂは、図８Ａに示す画像データに対してＦＦＴを施して得られる二次元パワースペクトルの、動径方向における平均値である。
【図９】Ｄｏｏｌｅｙ−Ｓｈａｗ関数（観察距離３００ｍｍ）のグラフである。
【図１０】図２に示す閾値マトリクス決定部の機能ブロック図である。
【図１１】図１０に示す閾値マトリクス決定部の動作説明に供されるフローチャートである。
【図１２】図１１のステップＳ３における詳細フローチャートである。
【図１３】図１３Ａは、ドットの初期位置を表す説明図である。図１３Ｂは、図１３Ａに示す状態に対して１個のドットを移動した結果を表す説明図である。
【図１４】図１２のステップＳ３７における詳細フローチャートである。
【図１５】図１５Ａは、ドットの形成を模擬した画像データを可視化した概略説明図である。図１５Ｂは、図１５Ａの破線矢印に沿った仮想的な明度プロファイルである。
【図１６】図１６Ａは、所定の作成条件下で最適化されたドットパターンを表す画像データを可視化した概略説明図である。図１６Ｂは、図１６Ａに示す画像データに対してＦＦＴを施して得られる二次元パワースペクトルの、動径方向における平均値である。
【図１７】第１変形例に係る画像形成装置の機能ブロック図である。
【図１８】第２変形例に係る閾値マトリクス作成部の機能ブロック図である。
【図１９】図１９Ａは、筋むらが発生し易いサブパターンの一例を表す概略説明図である。図１９Ｂは、第２変形例に係る第１評価値の算出方法を表す概略説明図である。
【図２０】印字パターンの画像形成に供される制御信号を可視化した概略説明図である。
【発明を実施するための形態】
【００３０】
以下、本実施の形態に係る閾値マトリクス作成方法について、それを実施する閾値マトリクス作成装置及び印刷システムとの関係において好適な実施形態を挙げ、添付の図面を参照しながら詳細に説明する。本明細書において、画像を形成することを「印刷」という場合がある。
【００３１】
図１は、本実施の形態に係る閾値マトリクス作成装置としての画像処理装置１２を組み込んだ印刷システム１０の概略説明図である。印刷システム１０は、画像処理装置１２と、画像形成装置１４と、ＤＴＰ（Desktop Publishing）装置１６と、データベースサーバ１８とを基本的に備える。画像処理装置１２、ＤＴＰ装置１６及びデータベースサーバ１８は、有線又は無線によって相互に電気的に接続されている。
【００３２】
画像処理装置１２は、外部装置からの入力画像データ（デバイス色信号又はページ記述データ）を画像形成装置１４側に出力する。ここで、デバイス色信号とは、デバイス依存データで定義された色信号を意味し、例えば、４色（ＣＭＹＫ）或いは３色（ＲＧＢ）のカラーチャンネルを有するラスタ形式データ（ＴＩＦＦ、ビットマップ、ＲＡＷ等）である。また、画像形成装置１４に供給されるデバイス色信号として、任意のヘッダを付加した独自のフォーマットデータを用いてもよい。また、画像処理装置１２は、画像形成装置１４に応じたハーフトーン処理に供される閾値マトリクスＭｔを作成可能である。
【００３３】
画像形成装置１４は、画像処理装置１２に電気的に接続されている。この接続には、例えば、ＵＳＢ（Universal Serial Bus）、ＩＥＥＥ１３９４、イーサネット（登録商標）、無線ネットワーク等のシリアルインターフェイスや、セントロニクス等のパラレルインターフェイスを適用することができる。
【００３４】
画像形成装置１４は、記録媒体としての枚葉紙（以下、用紙２０という。）を所定方向に搬送させながら、ヘッドユニット２２からインク２４の液滴を吐出させることで用紙２０上に画像を形成し、印刷物２６を取得する、いわゆるインクジェットプリンタである。用紙２０の基材には、合成紙・厚紙・アルミ蒸着紙等の紙類、塩化ビニル・ＰＥＴ等の樹脂やターポリン等を用いることができる。
【００３５】
画像形成装置１４は、ヘッドユニット２２の他、画像処理装置１２から入力された画像信号（以下、入力画像信号という。）を、ヘッドユニット２２を制御するための制御信号に変換する信号変換部２８と、該信号変換部２８により取得された制御信号に基づいてヘッドユニット２２を吐出制御するヘッドドライバ３０とを備えている。
【００３６】
なお、ヘッドユニット２２による液滴の吐出機構として、種々の方式を採り得る。例えば、ピエゾ素子（圧電素子）等で構成されるアクチュエータの変形によってインク滴を吐出する方式を適用してもよい。また、ヒータ等の発熱体を介してインク２４を加熱することで気泡を発生させ、その圧力でインク滴を吐出するサーマルジェット方式を適用してもよい。また、ヘッドユニット２２は、ラインヘッドに限定されることなく、用紙２０の幅方向に往復走査しながら画像を形成させるマルチパス方式であってもよい。
【００３７】
ＤＴＰ装置１６は、文字、図形、絵柄や写真等から構成される素材を編集可能であり、該素材をページ毎に配置することで、ページ記述言語（以下、ＰＤＬという。）による電子原稿を生成する。ここで、ＰＤＬとは、印刷や表示等の出力単位である「ページ」内で文字、図形等の書式情報、位置情報、色情報（濃度情報を含む）等の画像情報を記述する言語である。ＤＴＰ装置１６は、ＰＤＬ形式の電子原稿に対してラスタライズ処理を施す。このラスタライズ処理には、ＰＤＬ形式からラスタ形式に変換するデータ形式変換処理と、ＩＣＣ（International Color Consortium）プロファイルを用いた色変換処理とが含まれる。
【００３８】
データベースサーバ１８は、電子原稿のジョブチケット｛例えば、ＪＤＦ（Job Definition Format）ファイル｝、色見本データ、ターゲットプロファイル、又は画像形成装置１４及び用紙２０の組合せに適した印刷プロファイル等のデータ管理を行う装置である。
【００３９】
図２は、図１に示す画像処理装置１２の電気的なブロック図である。
【００４０】
画像処理装置１２は、入力インターフェース３２と、制御部３４と、メモリ３６と、出力インターフェース３８とを備えている。メモリ３６には、閾値マトリクス作成装置として機能させるためのプログラムが格納されている。
【００４１】
入力インターフェース３２は、外部装置からの電気信号を受信する。例えば、ＤＴＰ装置１６で編集・作成されたデバイス色信号やＰＤＬデータを取得する。また、データベースサーバ１８で記憶されたＩＣＣプロファイル等の各種情報を取得する。出力インターフェース３８は、外部装置に電気信号を送信する。例えば、画像処理装置１２は、出力インターフェース３８を介して、画像処理されたデバイス色信号又は後述する閾値マトリクスＭｔのデータを画像形成装置１４側に供給する。
【００４２】
ＣＰＵ等の情報処理装置で構成される制御部３４は、閾値マトリクス決定部４０と、画像処理部４２とを備えている。閾値マトリクス決定部４０は、画像形成装置１４でのハーフトーン処理に供される閾値マトリクスＭｔを決定する。画像処理部４２は、入力されたデバイス色信号に対して必要に応じた画像処理を施す。
【００４３】
図３は、図１に示す信号変換部２８における画像処理流れを表す概略説明図である。信号変換部２８は、解像度変換部４４と、ＣＭＹＫ色変換部４６と、ハーフトーン処理部４８とを基本的に備える。
【００４４】
信号変換部２８に入力される画像信号（入力画像信号）は、複数のカラーチャンネルからなる多階調データである。例えば、８ビット（１画素当り２５６階調）ＲＧＢのＴＩＦＦ形式データであってもよい。
【００４５】
解像度変換部４４は、画像サイズを拡大又は縮小する画像拡縮処理を用いて、入力画像信号の解像度を、画像形成装置１４に応じた出力解像度に変換する。ここで得られる第１中間画像信号は、入力画像信号とデータ定義は同一であるが、データサイズが異なっている。この画像拡縮処理には、補間演算を含む公知のアルゴリズムを種々適用してもよい。
【００４６】
ＣＭＹＫ色変換部４６は、公知のカラーマッチング手法を用いて、解像度変換部４４から取得した第１中間画像信号を、画像形成装置１４で取り扱うデバイス色信号（ＣＭＹＫ色信号）に変換する。ここで得られる第２中間画像信号は、多階調のＣＭＹＫ色信号に相当する。
【００４７】
ハーフトーン処理部４８は、ＣＭＹＫ色変換部４６から取得した第２中間画像信号を、インク２４を適切に吐出制御するための制御信号（ヘッドユニット２２の制御に供される信号）に変換する。ここで得られる制御信号は、各ラインヘッド８４（図５参照）に対してインク２４の吐出動作の有無（オン・オフ）を時系列的に制御するＣＭＹＫ毎の２値データ（あるいは多値データ）である。本実施の形態では、画像処理装置１２から供給された閾値マトリクスＭｔを用いて、組織的ディザ法によるハーフトーン処理方法について説明する。
【００４８】
図４は、組織的ディザ法によるハーフトーン処理の概略説明図である。一例として、ベイヤー型の閾値マトリクスＭｔを用いた２値化の概念を示す。先ず、多値ＣＭＹＫ色信号の各アドレスと、閾値マトリクスＭｔの各行列要素とを対応付ける。そして、着目する画素での画素値と、着目する行列要素での閾値との大小関係をそれぞれ比較し、画素値の方が大きい場合には「１（オン）」を割り当て、それ以外の場合には「０（オフ）」を割り当てる。このようにして、画像信号の階調数を多値から２値に変換することができる。
【００４９】
図５は、図１に示す画像形成装置１４の構成を表す断面側面図である。
【００５０】
画像形成装置１４には、用紙２０の搬送方向上流側に、用紙２０を給紙搬送する給紙搬送部５４が設けられている。この給紙搬送部５４の下流側には、用紙２０の搬送方向に沿って、用紙２０の記録面（以下、画像形成面という。）に処理液を塗布する処理液塗布部５６と、前記画像形成面にインクを付着することで画像を形成する画像形成部５８と、用紙２０上に形成された処理液層のインクを乾燥させるインク乾燥部６０と、処理液層の画像を用紙２０に定着させる画像定着部６２と、画像が定着した用紙２０を排出する排出部６４とが設けられている。
【００５１】
給紙搬送部５４は、用紙２０を積載可能に設けられた積載部６６と、該積載部６６に積載された用紙２０を一枚ずつ給紙する給紙部６８と、該給紙部６８により給紙された用紙２０を処理液塗布部５６に搬送する搬送部７０とを備えている。
【００５２】
処理液塗布部５６は、回転可能に設けられた処理液塗布ドラム７２と、用紙２０の画像形成面に処理液を塗布する処理液塗布装置７４と、前記処理液を乾燥する処理液乾燥装置７６とを備えている。これにより、用紙２０の画像形成面上に薄膜の処理液層が塗布される。
【００５３】
処理液塗布部５６と画像形成部５８との間には、回転可能に設けられた第１中間搬送ドラム７８が配置されている。第１中間搬送ドラム７８の表面に用紙２０を保持した状態で該第１中間搬送ドラム７８を回転させることにより、処理液塗布部５６側から供給された用紙２０は、画像形成部５８側へ搬送される。
【００５４】
画像形成部５８は、回転可能に設けられた画像形成ドラム８０（搬送部）と、該画像形成ドラム８０により搬送される用紙２０にインク２４の液滴を吐出するヘッドユニット２２とを備えている。ヘッドユニット２２は、少なくとも基本色であるＹ（イエロー）、Ｍ（マゼンタ）、Ｃ（シアン）、Ｋ（ブラック）のラインヘッド８４（ドット形成部）を備えている。そして、各ラインヘッド８４は、画像形成ドラム８０の周方向に沿って配列されている。これにより、用紙２０の画像形成面上に塗布された処理液層上に、各色の画像を順次形成させる。なお、この処理液には、インク２４の溶媒中に分散した色材（顔料）及びラテックス粒子を凝集させる効果を持たせているので、用紙２０上での色材流れ等を防止可能である。
【００５５】
画像形成部５８とインク乾燥部６０との間には、回転可能に設けられた第２中間搬送ドラム８６が配置されている。第２中間搬送ドラム８６の表面に用紙２０を保持した状態で該第２中間搬送ドラム８６を回転させることにより、画像形成部５８側から供給された用紙２０は、インク乾燥部６０側へ搬送される。
【００５６】
インク乾燥部６０は、回転可能に設けられたインク乾燥ドラム８８と、用紙２０の処理液層を乾燥する複数の熱風ノズル９０と、複数の赤外線ヒータ９２とを備えている。これにより、用紙２０の処理液層に滞留するインク２４の溶媒を乾燥させる。
【００５７】
インク乾燥部６０と画像定着部６２との間には、回転可能に設けられた第３中間搬送ドラム９４が配置されている。第３中間搬送ドラム９４の表面に用紙２０を保持した状態で該第３中間搬送ドラム９４を回転させることにより、インク乾燥部６０側から供給された用紙２０は、画像定着部６２側へ搬送される。
【００５８】
画像定着部６２には、回転可能に設けられた画像定着ドラム９６と、該画像定着ドラム９６の表面に近接して配置された加熱ローラ９８と、該画像定着ドラム９６の表面に圧接した状態で配置された定着ローラ１００とを備えている。これにより、処理液層で凝集するラテックス粒子が加熱・加圧されて溶融し、用紙２０上に画像として固定・定着される。
【００５９】
上記した各工程を経て、画像形成面の画像が定着した用紙２０は、画像定着ドラム９６の回転により、画像定着部６２の下流側に設けられた排出部６４側へ搬送される。
【００６０】
図６は、図５に示すラインヘッド８４の構造例を表す平面透視図である。ラインヘッド８４は、千鳥でマトリクス状に配列された複数のインク室ユニット１２０（ドット形成素子）を備えている。各インク室ユニット１２０は、ノズル１２２と、圧力室１２４と、供給口１２６とをそれぞれ備えている。平面形状が概略正方形である圧力室１２４には、その対角線上の両隅部の一方にノズル１２２側への流出口が設けられ、他方に図示しない共通流路からの流入口（供給口１２６）が設けられている。
【００６１】
以下、ノズル１２２の配置の特徴について説明する。本図において、ラインヘッド８４の長手方向及び短手方向をそれぞれ矢印Ｘ方向（主方向）、矢印Ｙ方向（副方向）と定義する。このとき、用紙２０の搬送方向（図５参照）は、矢印Ｘ方向に直交するとともに、矢印Ｙ方向に平行する。
【００６２】
第Ｌ１列における各ノズル１２２は、矢印Ｘ方向に沿って所定間隔（４単位長に相当する。）おきに等間隔に配置されている。第Ｌ２〜第Ｌ４列における各ノズルについても、第Ｌ１列と同様に配置されている。以下、矢印Ｘ方向を、ノズル１２２（インク室ユニット１２０）の「配列方向」ともいう場合がある。
【００６３】
第Ｌ２列の各ノズル１２２は、第Ｌ１列の各ノズル１２２の位置を基準として、矢印Ｘの左方向に１単位長だけシフトした位置に配置されている。第Ｌ３列の各ノズル１２２は、第Ｌ２列の各ノズル１２２の位置を基準として、矢印Ｘの左方向に１単位長だけシフトした位置に配置されている。第Ｌ４列の各ノズル１２２は、第Ｌ３列の各ノズル１２２の位置を基準として、矢印Ｘの左方向に１単位長だけシフトした位置に配置されている。これにより、ラインヘッド８４の長手方向に沿って並ぶように投影される実質的なノズル１２２の間隔（投影ノズルピッチ）の高密度化を達成している。
【００６４】
以下、閾値マトリクス決定部４０（図２参照）による閾値マトリクスＭｔの作成方法について、図７Ａ〜図１５Ｂを参照しながら詳細に説明する。閾値マトリクスＭｔは、色版（例えばＣＭＹＫ）によらず同一であってもよいし、色版毎に異ならせてもよい。
【００６５】
閾値マトリクスＭｔの作成方法の説明に先立って、ドットパターンＤＰＴの評価方法について説明する。本実施の形態では、筋むらに対する頑健性を定量化した第１評価値ＥＶ１と、ドットの分散に起因するノイズ粒状性を定量化した第２評価値ＥＶ２とに基づいて評価を行う。
【００６６】
先ず、第１評価値ＥＶ１の算出方法について説明する。この第１評価値ＥＶ１は、閾値マトリクスＭｔの行方向及び列方向のいずれか一方向におけるドットの重なり程度を評価するための評価値である。以下、一方向として、図６の矢印Ｘ方向（主方向）を選択する。
【００６７】
図７Ａは、ドットパターンＤＰＴの一部を表す概略説明図である。
【００６８】
本実施の形態では、閾値マトリクスＭｔの各行列要素の閾値を算出するために、ドットのオン・オフ状態（図４参照）を模式化した仮想的な矩形格子を用意する。以下、説明の便宜のため、前記矩形格子の各領域のことを「セル」という。本図では、行方向、列方向にそれぞれ５個の計２５個のセルが表記されている。ハッチングを付したセル（図７Ａ例では、３個のセル）は、ドットが配置されたオン状態のセルに相当する。なお、用紙２０上に形成されるドットとそれぞれ対応する、各セル上に配置された仮想的なドット（オン状態）についても、単に「ドット」という場合がある。
【００６９】
図７Ｂは、第１評価値ＥＶ１の算出方法を表す概略説明図である。より詳細には、本図は、図７ＡのドットパターンＤＰＴに基づいて各ドットの配置状態を模擬した説明図である。以下、着目セル１３０の位置に形成されるドットを第１ドットという。また、着目セル１３０に隣接する隣接セル１３２の位置に形成されるドットを第２ドットという。
【００７０】
ここで、第１ドット及び第２ドットが形成される中心位置（目標値）をそれぞれ点Ｐ、点Ｑとする。なお、線分ＰＱの距離をＤｄとする。
【００７１】
また、第１ドット及び第２ドットのサイズ（半径の目標値）をそれぞれＲ１、Ｒ２とする。さらに、第１ドット及び第２ドットのサイズ誤差（半径の変動量）をそれぞれＥｒｓ１、Ｅｒｓ２とする。さらに、第１ドット及び第２ドットの着弾位置誤差（以下、単に「位置誤差」という場合がある。）をそれぞれＥｒｐ１、Ｅｒｐ２とする。ここで、着弾位置誤差とは、点Ｐ（又は点Ｑ）に対する着弾位置のずれ量を意味し、Ｅｒｐ１及びＥｒｐ２は、いずれも矢印Ｘ方向の成分を表す。
【００７２】
なお、ラインヘッド８４は、複数のサイズのドットを形成可能であってもよく、Ｒ１及びＲ２の値が異なる場合がある。また、ラインヘッド８４は、３種類以上のサイズを形成可能であってもよい。
【００７３】
各ドットのサイズ誤差に着目した場合における、第１ドットと第２ドットとの重なり量ＯＬは、幾何学的考察により、次の（１）式で表される。
【数１】
2014096611.tif 000003
【００７４】
各ドットの着弾位置誤差に着目した場合における、矢印Ｘ方向に投影された重なり量ＯＬｘは、次の（２）式で表される。なお、θは、Ｘ軸と線分ＰＱとのなす角である。
【数２】
2014096611.tif 000004
【００７５】
このように、着目セル１３０に応じた評価領域１３４内の各セルに対して、（２）式で示した重なり量ＯＬｘをそれぞれ算出する。図７Ａ例では、矢印Ｘ方向に隣接する可能性が高い２つの領域（太枠で囲んだ領域）を評価領域１３４として定義する。なお、評価領域１３４のサイズ及び形状は任意に決定してもよい。
【００７６】
そして、着目セル１３０での個別評価値ＥＶ（ｉ）は、次の（３）式で表される。下記のΣは、添字ｊに対する総和を意味する。
【数３】
2014096611.tif 000005
【００７７】
ここで、ｉは、着目セル１３０のインデックスを表す。ｊは、着目セル１３０に応じた評価領域１３４内の各セルのインデックスを表す。ｍａｘ（ａ，０）は、ａ及び０のうち大きい値を返す関数である。
【００７８】
最後に、ドットパターンＤＰＴにおける第１評価値ＥＶ１は、次の（４）式で表される。下記のΣは、添字ｉに対する総和を意味する。
【数４】
2014096611.tif 000006
【００７９】
（３）式から諒解されるように、重なり量ＯＬｘが正値であるほど、第１評価値ＥＶ１が総じて大きくなる。すなわち、この第１評価値ＥＶ１が大きいほど、隣接するドット間のＸ方向離間量が低減されるので、ドットパターンＤＴＰ全体として、矢印Ｙ方向に延在する筋むらに対する頑健性が高くなる。一方、第１評価値ＥＶ１が小さいほど、ドットパターンＤＴＰ全体として、矢印Ｙ方向に延在する筋むらに対する頑健性が低くなる。
【００８０】
なお、ドットのサイズ誤差Ｅｒｓ１（Ｅｒｓ２）として、種々の統計量Ｅｒｓを用いることができる。また、液滴の着弾位置誤差Ｅｒｐ１（Ｅｒｐ２）として、種々の統計量Ｅｒｐを用いることができる。ここで、「統計量」とは、統計学的手法を用いて算出した物理量であり、例えば、平均値、最頻値、中心値、最大値、最小値等であってもよい。これらの統計量を得るため、前記統計量と対応付けられた各種特性値（例えばノズルの加工精度等）から推定してもよい。
【００８１】
また、第１評価値ＥＶ１は、主方向（矢印Ｘ方向）又は副方向（矢印Ｙ方向）のいずれか一方を考慮して算出してよいし、両方を考慮して算出してもよい。特に、各種誤差（Ｅｒｐ、Ｅｒｓ）が相対的に大きい方向に関する統計量であることが好ましい。
【００８２】
続いて、第２評価値ＥＶ２の算出方法について説明する。この第２評価値ＥＶ２は、各ドットの分散に起因するノイズ粒状性を評価するための評価値である。
【００８３】
図８Ａは、ドットパターンＤＰＴを表す画像データを可視化した概略説明図である。図８Ｂは、図８Ａに示す画像データに対してＦＦＴを施して得られる二次元パワースペクトルの、動径方向における平均値である。以下、平均化されたパワースペクトルを単にスペクトルＳＰＣという。
【００８４】
スペクトルＳＰＣは、低〜中空間周波数帯域にかけて強度レベルが小さくなり、高い空間周波数帯域のみ強度レベルが高くなるいわゆるハイパス型（いわゆるブルーノイズ型）の空間周波数特性を有している。
【００８５】
ノイズ粒状性を評価する例として、ドットパターンＤＰＴの各セルの画素値（オン・オフ状態値）に対してＲＭＳ（Root Mean Square）を求めてもよい。しかし、図８Ｂに示すように、ブルーノイズ型のスペクトルＳＰＣを有するドットパターンＤＰＴでは、観察者による視認性が良好であるにもかかわらず、ＲＭＳの値が大きくなる傾向がある。すなわち、評価値としてのＲＭＳと官能評価結果とが対応しない場合がある。そこで、本実施の形態では、人間の視覚応答特性を評価に採り入れ、さらに改良した第２評価値ＥＶ２を用いている。
【００８６】
図９は、人間の標準視覚応答特性の一例である、ドゥーリー・ショー（Dooley-Shaw）関数（観察距離３００ｍｍ）のグラフである。
【００８７】
本実施の形態では、人間の標準視覚応答特性として、観察距離３００ｍｍでのドゥーリー・ショー関数を用いている。ドゥーリー・ショー関数は、ＶＴＦ（Visual Transfer Function）の一種であり、人間の標準視覚応答特性を模した代表的な関数である。具体的には、輝度のコントラスト比特性の２乗値に相当する。なお、グラフの横軸は空間周波数（単位：ｃｙｃｌｅ／ｍｍ）であり、縦軸はＶＴＦの値（単位は無次元）である。
【００８８】
観察距離を３００ｍｍとすると、０〜１．０ｃｙｃｌｅ／ｍｍの範囲ではＶＴＦの値は一定（１に等しい。）であり、空間周波数が高くなるにつれて次第にＶＴＦの値が減少する傾向がある。すなわち、この関数は、中〜高空間周波数帯域を遮断するローパスフィルタとして機能する。
【００８９】
第２評価値ＥＶ２は、スペクトルＳＰＣの値をＦ（Ｕｘ，Ｕｙ）とするとき、次の（５）式で算出される。
【数５】
2014096611.tif 000007
【００９０】
ウィナー・ヒンチン（Wiener-Khintchen）の定理によれば、スペクトルＳＰＣを全空間周波数帯域で積分した値は、ＲＭＳの２乗値に一致する。このスペクトルＳＰＣに対してＶＴＦを乗算し、この新たなスペクトルＳＰＣを全空間周波数帯域で積分した値は、人間の視覚特性に略一致する評価指標となる。この第２評価値ＥＶ２は、人間の視覚応答特性で補正したＲＭＳということができる。通常のＲＭＳと同様に、第２評価値ＥＶ２は、常に０以上の値を取り、０に近づくほどノイズ粒状性が良好であるといえる。
【００９１】
また、図９に示すＶＴＦに対して逆フーリエ変換（例えば、ＩＦＦＴ）を施すことで、ＶＴＦに対応する実空間上のマスクを算出し、評価しようとするドットパターンＤＰＴに対して該マスクを作用して畳み込み演算を行い、新たなドットパターンＤＰＴに対してＲＭＳを求めてもよい。これにより、（５）式を用いた上記方法と同等の演算結果を得ることができる。
【００９２】
以上説明した第１評価値ＥＶ１と、第２評価値ＥＶ２とを組み合わせて、ドットパターンＤＰＴについての全体評価値ＥＶＡを算出することができる。全体評価値ＥＶＡの一例として、次の（６）式で与えられる。
【数６】
2014096611.tif 000008
【００９３】
この（６）式において、各評価項目の重要度に応じてα及びβの各値を定めることができる。例えば、非負値である第１評価値ＥＶ１は値が大きいほど望ましい特性であるとし、非負値である第２評価値ＥＶ２は０に近づくほど望ましい特性であるとする。このとき、α＝１、β＝１と定めることにより、筋むらの頑健性（第１評価値ＥＶ１）及びノイズ粒状性（第２評価値ＥＶ２）のバランスが十分取れた状態下において、全体評価値ＥＶＡの値は最大となる。ここで、α及びβは、理論的に定めた値であるか、経験的に見出した値であるかは問わない。
【００９４】
なお、ドットパターンＤＰＴを決定するための目標レベル（許容範囲）や評価関数に応じて、第１評価値ＥＶ１、第２評価値ＥＶ２及び全体評価値ＥＶＡの算出式を種々変更し得ることはいうまでもない。
【００９５】
図１０は、図２に示す閾値マトリクス決定部４０の機能ブロック図である。
【００９６】
閾値マトリクス決定部４０は、未だ設定していない階調レベル（後述する指定レベルｌｖ）を順次指定する階調レベル指定部１４０と、該階調レベル指定部１４０により指定された階調レベルでのドットパターンＤＰＴにおける制約条件を設定する制約条件設定部１４２と、指定された前記階調レベルに応じたドットパターンＤＰＴｆｉｘを決定するドットパターン決定部１４４と、該ドットパターン決定部１４４により決定されたドットパターンＤＰＴｆｉｘを格納するドットパターン格納部１４６と、該ドットパターン格納部１４６により格納されたドットパターンＤＰＴｆｉｘにおけるドットの配列に従って閾値マトリクスＭｔを決定する閾値決定部１４８と、前記ドットパターン決定部１４４での演算過程で生じる各種中間データを適宜格納する中間データ格納部１５０と、複数のドットが形成される位置誤差の統計量Ｅｒｐ、サイズ誤差の統計量Ｅｒｓ等のドット形成情報を格納するドット形成情報格納部１５１とを備える。
【００９７】
さらに、ドットパターン決定部１４４は、ドット数決定部１５２と、初期パターン作成部１５４と、ドット位置移動部１５６（ドットパターン作成部）と、評価値マップ算出部１５８と、全体評価値算出部１６０（評価値算出部）と、更新可否判別部１６２とを備えている。
【００９８】
以下、閾値マトリクス決定部４０の動作の概略について、図１１、図１２及び図１４のフローチャート、並びに図１０の機能ブロック図を参照しながら説明する。
【００９９】
本実施の形態では、閾値マトリクスＭｔの矩形格子上に配置するドット数を階調レベルに応じて決定し、ドットの配置とその評価を順次繰り返すことで、ドットの配置順を最適化する。そして、最適化された前記ドットの配置順に基づいて、閾値マトリクスＭｔの各行列要素での閾値を決定する。
【０１００】
図１１において、先ず、階調レベル指定部１４０は、未だ設定していない階調レベルを１つ指定する（ステップＳ１）。以下、この指定された階調レベルを「指定レベルｌｖ」という場合がある。この指定は、階調レベルの順番（昇順又は降順）に従ってもよいし、ランダムに順次指定してもよい。
【０１０１】
次いで、制約条件設定部１４２は、指定レベルｌｖでのドットパターンＤＰＴ（後述するＤＰＴｆｉｘ、ＤＰＴｉｎｉ及びＤＰＴｔｍｐも含まれる。）における制約条件を設定する（ステップＳ２）。ここで、制約条件とは、閾値マトリクスＭｔの各行列要素における閾値を矛盾なく決定するため、各階調レベルでのドットの配置を制約する条件である。具体的には、指定レベルｌｖよりも低い階調レベルで既に決定したドットの配置及び順番を、該指定レベルｌｖでも継承する必要がある。また、指定レベルｌｖよりも高い階調レベルで既に決定したドットの配置の中から、ドットの位置を選択する必要がある。このような制約条件を全て満たすように、既に確定した各階調レベルのドットパターンＤＰＴｆｉｘを読み出し、参照してもよい。
【０１０２】
次いで、ドットパターン決定部１４４は、指定レベルｌｖに応じたドットパターンＤＰＴｆｉｘ［ｌｖ］を決定する（ステップＳ３）。具体的な決定方法に関しては後述する。なお、ドットパターン格納部１４６は、ドットパターン決定部１４４で決定されたドットパターンＤＰＴｆｉｘ［ｌｖ］を取得し、一時的に格納する。
【０１０３】
次いで、階調レベル指定部１４０は、すべての階調レベルでドットパターンＤＰＴｆｉｘ［ｌｖ］が確定したか否かを判別する（ステップＳ４）。確定していないと判別された場合、ステップＳ１に戻って、以下、確定するまでステップＳ１〜Ｓ３を繰り返す。
【０１０４】
一方、確定したと判別された場合、閾値決定部１４８は、ドットパターン格納部１４６から取得したドットパターンＤＰＴｆｉｘを用いて、閾値マトリクスＭｔを決定する（ステップＳ５）。具体的には、指定レベルｌｖにおいて、ドットパターンＤＰＴｆｉｘ［ｌｖ＋１］と、ドットパターンＤＰＴｆｉｘ［ｌｖ］との差分を算出し、新たにドットが発生した１つの行列要素に対し、閾値（＝ｌｖ）を割り当てる。
【０１０５】
このようにして、閾値マトリクス決定部４０は、閾値マトリクスＭｔを決定する。
【０１０６】
続いて、図１１のステップＳ３における、指定レベルｌｖに応じたドットパターンＤＰＴｆｉｘ［ｌｖ］を決定する具体的方法について、図１２のフローチャート及び図１０の機能ブロック図を参照しながら詳細に説明する。
【０１０７】
本実施の形態では、ドットの配置を種々変化させたドットパターンＤＰＴの作成と、評価値（後述する全体評価値ＥＶＡ及び評価値マップＥＶ＿ＭＡＰ）による評価とを順次繰り返す方法を用いる。この場合、ドットパターンＤＰＴを決定する最適化問題として、構成的アルゴリズムや逐次改善アルゴリズム等の種々の探索アルゴリズムを用いることができる。
【０１０８】
ここでは、ボイドアンドクラスター法（Void-and-Cluster Method；以下、ＶＣ法という。）によるドットパターンＤＰＴの最適化方法について説明する。
【０１０９】
先ず、ドット数決定部１５２は、指定レベルｌｖに応じて、セルに配置する予定であるドットの個数を決定する（ステップＳ３１）。例えば、指定レベルｌｖの値又は閾値マトリクスＭｔのサイズ（セルの総数）と予め対応付けておく。
【０１１０】
次いで、初期パターン作成部１５４は、初期データとしてのドットパターンＤＰＴｉｎｉを作成する（ステップＳ３２）。初期パターン作成部１５４は、決定した個数のドットをセル上に配置することで、図１３Ａに例示するドットパターンＤＰＴｉｎｉを得る。
【０１１１】
ドットの配置は、例えば、擬似乱数の発生アルゴリズムを用いて乱数値を発生させ、該乱数値に基づいてランダムに決定してもよい。擬似乱数の発生アルゴリズムとして、メルセンヌ・ツイスタ（Mersenne Twister）、ＳＦＭＴ（SIMD-oriented Fast Mersenne Twister）やＸｏｒｓｈｉｆｔ法等の種々のアルゴリズムを用いてもよい。
【０１１２】
次いで、評価値マップ算出部１５８は、ステップＳ３２で作成したドットパターンＤＰＴｉｎｉから、初期データとしての評価値マップＥＶ＿ＭＡＰｉｎｉを算出する（ステップＳ３３）。評価値マップＥＶ＿ＭＡＰｉｎｉは、所定の評価関数に基づいてセル毎に算出された各評価値で構成されている。この評価値が大きいほど、良好な各種性能が得られるように設計されている。例えば、評価値マップＥＶ＿ＭＡＰｉｎｉは、セル毎の評価値ＥＶ（ｉ）｛（３）式及び図７Ｂを参照｝の集合で与えられる。評価値ＥＶ（ｉ）を算出する際には、ドット形成情報格納部１５１からドット形成情報を予め読み出しておく。ここで、「ドット形成情報」とは、ドットの形成に関連する情報であり、位置誤差の統計量Ｅｒｐ、サイズ誤差の統計量Ｅｒｓ、出力解像度（矢印Ｘ方向・矢印Ｙ方向）が含まれてもよい。また、ドット形成情報格納部１５１は、記録媒体の種類（又はインクの種類）に応じたドット形成情報をそれぞれ格納してもよい。
【０１１３】
次いで、全体評価値算出部１６０は、ステップＳ３３で作成した評価値マップＥＶ＿ＭＡＰｉｎｉから、初期データとしての全体評価値ＥＶＡを算出する（ステップＳ３４）。評価値マップＥＶ＿ＭＡＰｉｎｉは、各セルでの個別の評価値であるのに対し、全体評価値ＥＶＡは、ドットの全体配列に対する評価値である。全体評価値ＥＶＡの算出方法は、種々の方法を取り得る。例えば、全体評価値ＥＶＡは、上記した（６）式で与えられる。
【０１１４】
次いで、更新可否判別部１６２は、整数である変数Ｋに０を代入し、Ｋを初期化する（ステップＳ３６）。ここで、Ｋは、ドットパターンＤＰＴの更新要否の判別回数を表すカウンタである。
【０１１５】
次いで、ドット位置移動部１５６は、第１評価値ＥＶ１及び第２評価値ＥＶ２に基づいた所定の評価基準に従いつつ、制約条件を満たすようにドットの位置を移動することで、ドットパターンＤＰＴｔｍｐを取得する（ステップＳ３７）。この取得方法について、図１４のフローチャートを参照しながらさらに詳細に説明する。
【０１１６】
本実施の形態では、ドットの移動（セル対の交換）を決定するための評価値マップＥＶ＿ＭＡＰに対し、人間の視覚応答特性に応じた補正を行う。具体的には、視認され易いセル（密集したドット）における評価値が最大となり、ドットがオフ状態であるセルにおける評価値が最小となる。これにより、ドットが移動する成功確率が向上し、処理演算を高速化できる。
【０１１７】
図１４において、先ず、ドット位置移動部１５６は、ドットの形成を模擬した画像データを作成する（ステップＳ３７１）。以下、この画像データを模擬画像ＳＩＭ＿ＩＭＧという。
【０１１８】
図１５Ａは、図７ＡにおけるドットパターンＤＰＴに基づいて各ドットの配置状態を模擬した図である。例えば、第１ドットの形成領域２００及び第２ドットの形成領域２０２は、それぞれ均一な色を有すると仮定する。このとき、図１５Ｂに示すように、破線矢印に沿った明度（Ｌ*）プロファイルが推定される。このグラフは、信号値が明度で定義された模擬画像ＳＩＭ＿ＩＭＧの断面プロファイルである。
【０１１９】
このように、模擬画像ＳＩＭ＿ＩＭＧのサンプル間隔を、閾値マトリクスＭｔのサンプル間隔よりも小さく設定することで、例えば、形成領域２００と形成領域２０２との重複領域２０４での色変化を精度良く模擬できる。
【０１２０】
次いで、ドット位置移動部１５６は、模擬画像ＳＩＭ＿ＩＭＧに対して、視覚フィルタを作用する（ステップＳ３７２）。視覚フィルタとして、例えば、図９に示すＶＴＦの逆変換特性を用いてもよい。各セルに対応する位置での算出値は、視覚補正マップＶＦ＿ＭＡＰとして中間データ格納部１５０に一時的に格納される。なお、人間の視覚応答特性と精度良く対応させるために、模擬画像ＳＩＭ＿ＩＭＧのサンプル間隔に応じた視覚フィルタを用いることが好ましい。
【０１２１】
次いで、ドット位置移動部１５６は、評価値マップＥＶ＿ＭＡＰを補正する（ステップＳ３７３）。具体的には、現時点での評価値マップＥＶ＿ＭＡＰに対して、ステップＳ３７２で算出された視覚補正マップＶＦ＿ＭＡＰを乗算する演算をセル毎に実行する。以下、得られたマップを、補正済みマップＥＶＣ＿ＭＡＰという。
【０１２２】
次いで、ドット位置移動部１５６は、すべての交換組合せを有効にする（ステップＳ３７４）。交換組合せとは、２個のセル間におけるドットのオン・オフ状態を交換する組合せである。図１３Ａ例では、２５６個の行列要素のうちから２個の組合せを選択することから、３２６４０（＝２５６・２５５／２）通りの組合せを有効にする。
【０１２３】
次いで、ドット位置移動部１５６は、補正済みマップＥＶＣ＿ＭＡＰの差分値が最大となる一対のセル（セル対）を決定する（ステップＳ３７５）。すなわち、補正済みマップＥＶＣ＿ＭＡＰの値が最大であるセルと、最小であるセルとを１個ずつ抽出する。
【０１２４】
次いで、ドット位置移動部１５６は、ステップＳ３７５で選択したセル対を交換することで、現時点でのドットパターンＤＰＴが変化するか否かを判別する（ステップＳ３７６）。具体的には、選択されたセル対のうちドット（オン状態）が０個又は２個である場合、ドットパターンＤＰＴは変化しない。また、選択されたセル対のうちドット（オン状態）が１個である場合、ドットパターンＤＰＴは変化する。
【０１２５】
変化しないと判別された場合、ステップＳ３７５で選択されたセル対による交換組合せを除外し（ステップＳ３７８）、補正済みマップＥＶＣ＿ＭＡＰの差分値が次に最大となる他のセル対を順次選択する（ステップＳ３７５）。
【０１２６】
一方、変換すると判別された場合、ドット位置移動部１５６は、指定レベルｌｖでの制約条件を満たすか否かをさらに判別する（ステップＳ３７７）。制約条件を満たさないと判別された場合、ステップＳ３７５で決定されたセル対による交換組合せを除外し（ステップＳ３７８）、補正済みマップＥＶＣ＿ＭＡＰの差分値が次に最大となる他のセル対を順次選択する（ステップＳ３７５）。
【０１２７】
制約条件を満たすと判別された場合、ドット位置移動部１５６は、ドットパターンＤＰＴに対してセル対のオン・オフ状態を交換し、新たなドットパターンＤＰＴｔｍｐを取得する（ステップＳ３７９）。
【０１２８】
図１３Ａにおいて、第１セル２０６における補正済みマップＥＶＣ＿ＭＡＰが最大であり、第２セル２０８における補正済みマップＥＶＣ＿ＭＡＰが最小であったとする。このとき、図１３Ｂに示すように、第１セル２０６及び第２セル２０８のオン・オフ状態が反転することで、第１セル２０６の位置に存在するドットが第２セル２０８の位置に移動する。
【０１２９】
このようにして、ドット位置移動部１５６は、制約条件を満たすようにドットの位置を移動することで、ドットパターンＤＰＴｔｍｐを取得する（ステップＳ３７）。
【０１３０】
次いで、図１２に戻って、評価値マップ算出部１５８は、ドットパターンＤＰＴｔｍｐから評価値マップＥＶ＿ＭＡＰｔｍｐを算出する（ステップＳ３８）。そして、全体評価値算出部１６０は、評価値マップＥＶ＿ＭＡＰｔｍｐから全体評価値ＥＶＡを算出する（ステップＳ３８）。これらの算出方法は、ステップＳ３２及びＳ３３での算出方法と同一である。
【０１３１】
次いで、更新可否判別部１６２は、ドットパターンＤＰＴｔｍｐの全体評価値ＥＶＡｔｍｐと、現時点での最大値である全体評価値ＥＶＡとの大小関係を比較する（ステップＳ３９）。ＥＶＡｔｍｐ＞ＥＶＡを満たさない（すなわち、ＥＶＡｔｍｐ≦ＥＶＡである）場合は、ステップＳ４０を行うことなく、次のステップＳ４１に進む。
【０１３２】
一方、ＥＶＡｔｍｐ＞ＥＶＡを満たす場合、ステップＳ３７で求めたドットパターンＤＰＴｔｍｐを、現時点での最適なドットパターンであると判別する。このとき、中間データ格納部１５０は、現在記憶しているドットパターンＤＰＴに対し、ステップＳ３７で求めたドットパターンＤＰＴｔｍｐを上書き更新する（ステップＳ４０）。また、中間データ格納部１５０は、現在記憶している全体評価値ＥＶＡに対し、ステップＳ３８で求めた評価値ＥＶＡｔｍｐを上書き更新する（ステップＳ４０）。さらに、中間データ格納部１５０は、現在記憶している評価値マップＥＶ＿ＭＡＰに対し、ステップＳ３８で求めた評価値マップＥＶ＿ＭＡＰｔｍｐを上書き更新する（ステップＳ４０）。その後、次のステップＳ４１に進む。
【０１３３】
次いで、更新可否判別部１６２は、現時点でのＫの値を１だけ加算する（ステップＳ４１）。
【０１３４】
次いで、更新可否判別部１６２は、現時点でのＫの値と予め定められたＫｍａｘの値との大小関係を比較する（ステップＳ４２）。Ｋの値の方が小さい場合は、更新可否判別部１６２は、現時点での全体評価値ＥＶＡと、予め定められた許容値ＥＶＡ＿ＯＫとの大小関係をさらに比較する（ステップＳ４３）。ＥＶＡ＞ＥＶＡ＿ＯＫを満たさない（すなわち、ＥＶＡ≦ＥＶＡ＿ＯＫである）場合、ステップＳ３７まで戻り、以下ステップＳ３７〜Ｓ４１を順次繰り返す。なお、この最適化演算における収束性を十分確保するため、例えば、Ｋｍａｘ＝１００００と定めることができる。
【０１３５】
第１の終了条件（Ｋ＞Ｋｍａｘ）又は第２の終了条件（ＥＶＡ＞ＥＶＡ＿ＯＫ）の少なくとも一方を満たす場合、更新可否判別部１６２は、現時点でのドットパターンＤＰＴを、暫定的に確定したドットパターンＤＰＴとして決定する（ステップＳ４４）。
【０１３６】
このようにして、ドットパターン決定部１４４は、最後に更新されたドットパターンＤＰＴを、ドットパターンＤＰＴｆｉｘ［ｌｖ］として確定する（図１１のステップＳ３参照）。
【０１３７】
図１６Ａは、所定の作成条件下で最適化されたドットパターンＤＰＴｆｉｘを表す画像データを可視化した概略説明図である。図１６Ｂは、図１６Ａに示す画像データに対してＦＦＴを施して得られる二次元パワースペクトルの、動径方向における平均値である。
【０１３８】
図１６Ａに示すドットパターンＤＰＴｆｉｘは、図７ＡのドットパターンＤＰＴと比べて、一方向（行方向）に連結するドットの割合が多い。これにより、搬送方向（矢印Ｙ方向）に延在する筋むらに対する頑健性が向上する。
【０１３９】
このように、複数のドットが形成される位置誤差（Ｅｒｐ１、Ｅｒｐ２）の統計量Ｅｒｐ又はサイズ誤差（Ｅｒｓ１、Ｅｒｓ２）の統計量Ｅｒｓに基づいて、閾値マトリクスＭｔの行方向及び列方向のいずれか一方向におけるドットの重なり程度を定量化した第１評価値ＥＶ１を算出するようにしたので、ドットの形成ばらつきに起因する画像への影響度を統計的且つ定量的に予測可能である。そして、第１評価値ＥＶ１に基づいて閾値マトリクスＭｔを適切に決定することで、前記一方向に延在する筋むらに対する頑健性が高くなる。
【０１４０】
本実施の形態は、以下の変形例１〜３のように構成してもよい。なお、本実施の形態に係る画像処理装置１２及び画像形成装置１４（図１〜図１６Ｂ参照）と同じ構成要素については、同じ参照符号を付して、その詳細な説明を省略する。
【０１４１】
＜第１変形例＞
本実施の形態では、閾値マトリクス決定部４０は、画像処理装置１２に設けてられているが（図２参照）、画像形成装置２２０に一体的に組み込んでもよい。
【０１４２】
図１７は、第１変形例に係る画像形成装置２２０の機能ブロック図である。
【０１４３】
画像形成装置２２０は、信号変換部２８、ヘッドドライバ３０及びヘッドユニット２２の他、閾値マトリクス決定部４０（図２及び図１０参照）と同等の機能を有する閾値マトリクス決定部４０Ａをさらに備えている。これにより、画像形成装置２２０単体でハーフトーン処理を施すことが可能であり、別装置による演算部（図２の閾値マトリクス決定部４０）が不要となる。
【０１４４】
＜第２変形例＞
本実施の形態では、全体評価値算出部１６０は、幾何学的計算により第１評価値ＥＶ１を算出しているが、他の手法（例えばパターンマッチング法）を用いて算出してもよい。
【０１４５】
図１８は、第２変形例に係る閾値マトリクス決定部４０Ｂの機能ブロック図である。閾値マトリクス決定部４０Ｂは、図１０に示す各部の他、閾値マトリクスＭｔよりも小さいサイズの部分的なドットパターン（以下、サブパターンＳＰＴという。）を複数決定するサブパターン決定部２３０と、該サブパターン決定部２３０により決定された各サブパターンＳＰＴにおけるドットの重なり程度の頑健性を定量化し、各得点を決定する頑健性得点決定部２３２とをさらに備えている。そして、評価値マップ算出部１５８Ｂ（又は全体評価値算出部１６０Ｂ）は、第２変形例における算出方法を用いて、評価値マップＥＶ＿ＭＡＰ等（又は全体評価値ＥＶＡ等）を算出する。
【０１４６】
以下、閾値マトリクス決定部４０Ｂの動作について説明する。この動作は、図１１、図１２及び図１４で示したフローチャートと基本的に同様であるので、全体の説明を省略する。本実施の形態と異なるステップＳ３３、Ｓ３４及びＳ３８（図１２参照）について、図１９Ａ及び図１９Ｂを参照しながら説明する。
【０１４７】
ステップＳ３３（又はＳ３８）において、評価値マップＥＶ＿ＭＡＰｉｎｉ（ＥＶ＿ＭＡＰｔｍｐ）の算出に先立ち、サブパターン決定部２３０は、後述するパターンマッチングに用いるサブパターンＳＰＴを決定する。そして、頑健性得点決定部２３２は、サブパターン決定部２３０により決定された各サブパターンＳＰＴにおけるドットの重なり程度の頑健性を定量化する。
【０１４８】
図１９Ａは、筋むらが発生し易いサブパターンＳＰＴの一例を表す概略説明図である。本図では、行方向、列方向ともに３セルであるサブパターンＳＰＴを表記している。各サブパターンＳＰＴは、各ドットは上下方向に隣接せずに、斜め方向に隣接する特徴がある。これらの場合、着弾位置誤差Ｅｒｐ１、Ｅｒｐ２が次第に大きくなるにつれ、一方向（行方向）に延存する筋むらが発生する可能性が急激に高くなる。このように、筋むらが発生し易いサブパターンＳＰＴを選択的に抽出する。
【０１４９】
なお、サブパターンの決定又は定量化に際し、画像形成装置１４の出力特性（出力解像度、ドットのサイズＲ１、Ｒ２、位置誤差の統計量Ｅｒｐ、サイズ誤差の統計量Ｅｒｓに応じて種々変更してもよい。
【０１５０】
そして、評価値マップ算出部１５８Ｂは、パターンマッチング手法を用いて、各セルでの個別評価値ＥＶ（ｉ）を算出する。図１９Ｂは、第２変形例に係る第１評価値ＥＶ１の算出方法を表す概略説明図である。
【０１５１】
図１９Ｂに示すように、ドットパターンＤＰＴｉｎｉ（又はＤＰＴｔｍｐ）の部分パターンと、サブパターン決定部２３０で予め決定されたサブパターンＳＰＴとを照合する。ここで、部分パターンとは、ドットパターンＤＰＴｉｎｉ（又はＤＰＴｔｍｐ）のうち一部のドットパターンをいう。本変形例では、着目セル（部分パターンの中心セル）の位置を順次変更しながらパターンマッチングを行う。本図例のように、横方向にラスタ走査してもよい。
【０１５２】
部分パターンとサブパターンＳＰＴとのドット状態が完全に一致する場合、評価値マップ算出部１５８Ｂは、該サブパターンＳＰＴに応じた得点を加算する。すべてのサブパターンＳＰＴでの照合が完了した後、評価値マップ算出部１５８Ｂは、その得点の合計を、前記着目セルでの個別評価値ＥＶ（ｉ）として算出する。そして、着目セルを順次変更しながら、個別評価値ＥＶ１（ｉ）を算出することで、評価値マップＥＶ＿ＭＡＰを得る。
【０１５３】
同様に、ステップＳ３４（又はＳ３８）において、全体評価値算出部１６０Ｂは、評価値マップ算出部１５８Ｂにより算出された評価値マップＥＶ＿ＭＡＰに基づいて全体評価値ＥＶＡを算出する。
【０１５４】
このように構成することで、演算量が少ない比較的簡便な方法で第１評価値ＥＶ１を算出できる。
【０１５５】
＜第３変形例＞
本実施の形態では、各種特性値（ノズルの成形加工精度等）から統計量Ｅｒｐ、Ｅｒｓを推定しているが、実際の印刷結果に基づいて計測してもよい。例えば、特定の印字パターン２４０を印刷した後に画像を解析することで、各種統計量を簡便に取得可能である。
【０１５６】
図２０は、印字パターン２４０の吐出形成に供される制御信号を可視化した概略説明図である。信号領域２４２内に表記された複数のセルは、各ノズル１２２（図６参照）における制御信号を時系列的に表している。具体的には、黒色のセルはオン状態を表し、白色のセルはオフ状態を表す。すなわち、用紙２０上において、黒色のセルに対応する二次元位置に各ドットが形成されることとなる。
【０１５７】
信号領域２４２には、矢印Ｙ方向（搬送方向）に延在する複数の線画像２４４が規則的に形成されている。複数のノズル１２２のうち１つのノズルを選択し、そのノズル単体で、矢印Ｙ方向（搬送方向）に延在する線画像２４４を１本形成する。ノズル１２２毎に線画像２４４を形成しているので、線画像２４４の総数は、ラインヘッド８４（図６参照）に設けられたノズル１２２の総数と一致する。なお、図２０例では、線画像２４４の幅は１セル（ドット）に相当し、長さは１０セル（ドット）に相当する。
【０１５８】
以下、図２０に示す印字パターン２４０を用いて各種統計量を計測する方法について具体的に説明する。印字パターン２４０を表す制御信号をヘッドドライバ３０（図１参照）に供給し、適切な印刷制御をすることで、用紙２０上に印字パターン２４０が形成される。その後、図示しない画像読取装置（いわゆるスキャナ）で印字パターン２４０を読み取った後、読取画像データに対して各種解析手法を施して各種統計量を計測する。以下、具体例を示す。
【０１５９】
先ず、公知のパターンマッチング手法を用いて、読取画像データと、図２０に示す制御信号データとの位置合わせ処理を行う。なお、制御信号データは、画像読取装置の入力解像度と一致させるように拡縮する必要がある。
【０１６０】
次いで、読取画像データから線画像２４４をそれぞれ検出し、制御信号データでの位置（理想位置）に対する読取画像データでの位置ずれ量を計測する。実際の位置ずれ量は、読取画像データでの画素数と、単位画素サイズとの積で計算可能である。このように、各ノズル１２２でのドットの位置誤差（統計量Ｅｒｐ）を計測することができる。
【０１６１】
次いで、検出された各線画像２４４の幅及び長さを計測する。実際の各幅及び長さは、読取画像データでの画素数と、単位画素サイズとの積で計算可能である。このように、各ノズル１２２でのドットのサイズ誤差（統計量Ｅｒｓ）を計測することができる。
【０１６２】
最後に、位置誤差の統計量Ｅｒｐ及びサイズ誤差の統計量Ｅｒｓをノズル１２２毎に計測した後、画像形成装置１４の図示しないメモリに格納する。あるいは、閾値マトリクスＭｔを作成するため、画像処理装置１２（例えば、図１０のドット形成情報格納部１５１）に格納する。これにより、閾値マトリクスＭｔの作成の際、ドット形成情報を参照できるようになる。
【０１６３】
なお、印字パターン２４０の計測には、上記した画像読取装置の他、種々の計測装置を用いてもよい。画像読取装置を用いる場合、その種類に応じて線画像２４４の幅又は長さを適宜変更してもよい。例えば、線画像２４４を充分に長くすることで、画像を読み取る際に、塵埃の付着等による誤差の影響を受け難くなる。ラインヘッド８４により複数のサイズのドットを形成可能である場合、ドットサイズに応じて統計量Ｅｒｐ、Ｅｒｓを計測してもよい。
【０１６４】
また、印字パターン２４０の形成の際、インク滴同士の着弾干渉が発生しないように、線画像２４４の配置を適宜変更してもよい。例えば、サイズが大きいドットを形成する場合、線画像２４４の間隔を通常よりも大きくしてもよい。互いに隣接するドットにおいて、インク滴の着弾時点に十分な時間差を設けても同様の作用効果が得られる。
【０１６５】
このように、図２０に示す印字パターン２４０を用いることで、位置誤差の統計量Ｅｒｐ及びサイズ誤差の統計量Ｅｒｓを簡便に計測できる。
【０１６６】
なお、この発明は、上述した実施形態に限定されるものではなく、この発明の主旨を逸脱しない範囲で自由に変更できることは勿論である。
【０１６７】
例えば、本実施の形態では主にＣＭＹＫ（４つの色版）を中心に説明したが、これに限定されることなく、任意の色版の種類及び版数に設計変更できる。例えば、ＣＭＹＫの標準インクと、ＬＣ、ＬＭ等の淡色やＷ（白色）等のオプションインクとを組み合わせてもよい。
【０１６８】
また、本実施の形態では、画像形成ドラム８０の回転により用紙２０のみを搬送させているが、ヘッドユニット２２及び用紙２０のうち少なくとも一方を搬送させればよい。両者を相対移動させる構成であれば、本発明を適用できるからである。
【０１６９】
さらに、閾値マトリクスＭｔの行方向若しくは列方向のサイズ、又は閾値のレベル数は、本実施の形態に限定されず、適宜変更してもよいことは言うまでもない。
【符号の説明】
【０１７０】
１０…印刷システム                        １２…画像処理装置
１４、２２０…画像形成装置                ２２…ヘッドユニット
４０（Ａ、Ｂ）…閾値マトリクス決定部      １２２…ノズル
１４４…ドットパターン決定部              １４６…ドットパターン格納部
１４８…閾値決定部                        １５１…ドット形成情報格納部
１５２…ドット数決定部                    １５４…初期パターン作成部
１５６…ドット位置移動部                  １５８（Ｂ）…評価値マップ算出部
１６０（Ｂ）…全体評価値算出部            １６２…更新可否判別部
２３０…サブパターン決定部                ２３２…頑健性得点決定部
ＤＰＴ、ＤＰＴｆｉｘ、ＤＰＴｉｎｉ、ＤＰＴｔｍｐ…ドットパターン
Ｅｒｐ、Ｅｒｓ…統計量                    Ｍｔ…閾値マトリクス
【図１】
2014096611.tif 000009
【図２】
2014096611.tif 000010
【図３】
2014096611.tif 000011
【図４】
2014096611.tif 000012
【図５】
2014096611.tif 000013
【図６】
2014096611.tif 000014
【図７】
2014096611.tif 000015
【図８】
2014096611.tif 000016
【図９】
2014096611.tif 000017
【図１０】
2014096611.tif 000018
【図１１】
2014096611.tif 000019
【図１２】
2014096611.tif 000020
【図１３】
2014096611.tif 000021
【図１４】
2014096611.tif 000022
【図１５】
2014096611.tif 000023
【図１６】
2014096611.tif 000024
【図１７】
2014096611.tif 000025
【図１８】
2014096611.tif 000026
【図１９】
2014096611.tif 000027
【図２０】
2014096611.tif 000028