(57)【要約】【課題】  特殊な計測装置を使うことがなく、ドットの位置及び位置ずれを短時間で、しかも高精度に測定できるドット位置検出装置を提供する。【解決手段】  インクジェットプリンタで形成された画像におけるドット位置を予測するドット位置予測部３０５、予測されたドット位置に基づいて画像におけるドットの位置を重心として算出するドット重心位置計算部３０７、算出されたドットの重心に基づいて、ドットが形成されるべき理想位置を算出する理想位置計算部３０９、算出されたドットの理想位置とのずれ量を算出する位置ずれ計算部３１１とを備える。【選択図】  図３
2006142546.tif 000002
【特許請求の範囲】
【請求項１】
インクジェットプリンタのドット位置を検出するドット位置検出装置であって、
インクジェットプリンタで形成された画像におけるドット位置を予測するドット位置予測手段と、
前記ドット位置予測手段によって予測されたドット位置に基づいて、前記画像におけるドットの位置を算出するドット位置算出手段と、
前記ドット位置算出手段によって算出されたドットの位置に基づいて、ドットが形成されるべき理想位置を算出する理想位置算出手段と、
前記ドット位置算出手段によって算出されたドットの位置と、前記理想位置算出手段によって算出されたドットの理想位置とのずれ量を算出するドット位置ずれ量算出手段と、
を備えることを特徴とするドット位置検出装置。
【請求項２】
前記ドット位置予測手段が画像中のドットとドットが無い領域とを輝度にかかる情報に基づいて判断すると共に、ドットとドットが無い領域とを判断するための輝度にかかるしきい値を設定するしきい値設定手段を備えることを特徴とする請求項１に記載のドット位置検出装置。
【請求項３】
前記ドット位置予測手段は、前記画像における第１の方向について輝度の値のヒストグラムを作成すると共に該ヒストグラムの極小点を通る第１の直線を求め、
前記第１の方向と直交する第２の方向について輝度の値のヒストグラムを作成すると共に該ヒストグラムの極小点を通る第２の直線を求め、
前記第１の直線と前記第２の直線との交点をドットの位置であると予測することを特徴とする請求項１または２に記載のドット位置検出装置。
【請求項４】
前記ドット位置算出手段は、前記ドット位置予測手段によって予測されたドット位置に基づいて、前記画像におけるドットの半径をさらに算出することを特徴とする請求項１ないし３のいずれか１項に記載のドット位置検出装置。
【請求項５】
前記ドット位置算出手段は、前記ドット位置予測手段によって予測されたドット位置を含む領域の輝度にかかるデータを互いに直交する二方向について取得することによってドットの範囲を計測し、計測された範囲からドットの重心点を算出することを特徴とする請求項１から４のいずれか１項に記載のドット位置検出装置。
【請求項６】
前記ドット位置算出手段は、予測されたドット位置を含む領域の輝度にかかるデータを画像に対する水平方向と垂直方向、水平方向と垂直方向とに対して傾きを持った二方向の少なくとも一方について取得することを特徴とする請求項５に記載のドット位置検出装置。
【請求項７】
前記ドット位置算出手段は、前記ドット位置予測手段によって予測されたドット位置に基づいて算出されたドットの位置に基づいて、さらに前記画像におけるドット位置を算出することを特徴とする請求項１から６のいずれか１項に記載のドット位置検出装置。
【請求項８】
前記理想位置算出手段は、前記ドット位置算出手段によって算出されたドットの位置に基づいてドットを通る直線を近似して算出し、
算出された直線の傾きの平均を傾きとする理想直線を算出し、算出された理想直線上をドットの理想位置とすることを特徴とする請求項１から７のいずれか１項に記載のドット位置検出装置。
【請求項９】
前記理想位置算出手段は、前記理想直線とドットの形成領域上で設定された軸との交点から理想直線の設定開始位置を求めることを特徴とする請求項８に記載のドット位置検出装置。
【請求項１０】
インクジェットプリンタのドット位置を検出するドット位置検出方法をコンピュータに実行させるためのプログラムであって、
インクジェットプリンタで形成された画像におけるドット位置を予測するドット位置予測ステップと、
前記ドット位置予測ステップにおいて予測されたドット位置に基づいて、前記画像におけるドットの位置を算出するドット位置算出ステップと、
前記ドット位置算出ステップにおいて算出されたドットの位置に基づいて、ドットが形成されるべき理想位置を算出する理想位置算出ステップと、
前記ドット位置算出ステップにおいて算出されたドットの位置と、前記理想位置算出ステップにおいて算出されたドットの理想位置とのずれ量を算出するドット位置ずれ量算出ステップと、
を含むことを特徴とするドット位置検出方法をコンピュータに実行させるためのプログラム。
【請求項１１】
前記ドット位置予測ステップは、前記画像における第１の方向について輝度の値のヒストグラムを作成すると共に該ヒストグラムの極小点を通る第１の直線を求めるステップと、
前記第１の方向と直交する第２の方向について輝度の値のヒストグラムを作成すると共に該ヒストグラムの極小点を通る第２の直線を求めるステップと、を含み、
前記第１の直線と前記第２の直線との交点をドットの位置であると予測することを特徴とする請求項１０に記載のドット位置検出方法をコンピュータに実行させるためのプログラム。
【請求項１２】
前記理想位置算出ステップは、前記ドット位置算出ステップにおいて算出されたドットの位置に基づいてドットを通る直線を近似して算出するステップと、
算出された直線の傾きの平均を傾きとする理想直線を算出するステップと、を含み、
算出された理想直線上をドットの理想位置とすることを特徴とする請求項１０または１１に記載のドット位置検出方法をコンピュータに実行させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、インクジェットプリンタで印刷された画像を構成するドットの位置を検出するドット位置検出装置及びドット位置検出方法及をコンピュータに実行させるためのプログラムにかかり、特に印刷された画像からドットの位置を検出するドット位置検出装置及びドット位置検出方法及をコンピュータに実行させるためのプログラムに関する。
【背景技術】
【０００２】
インクジェット式のプリンタは、ヘッドに多数の微細なノズルを設け、ノズルによって噴出されたインクの液滴を用紙に付着させ、用紙上に画像を形成する（用紙上への画像の形成を本明細書では印刷とも記す。）。噴出された液滴は、用紙に到達し、到達点を中心にして略円形に広がってドットを形成する。インクジェットプリンタで印刷される画像の品質は、ドットの形状や大きさが均一であるほど高くなる。
【０００３】
ところで、近年、インクジェットプリンタには益々印刷の高速化が求められている。高速化の要請に応えるため、ノズルが設けられているヘッドを複数並べたラインヘッドプリントが普及してきている。ラインヘッドプリンタは、特にノズル特性の影響がスジやシェーディングといった不具合として画像に表れやすい。このため、ラインヘッドプリンタに使用されるヘッドには、特に形成されるドットの均一性が求められる。
【０００４】
ドット形状や大きさを均一化する方法としては、ノズルのインク噴出孔の位置や角度、大きさ（ノズル特性）をそろえることによってより均一化する方法と、プリンタ等、画像制御装置の画像形成条件を制御してドットの大きさ等を変化させる方法とがある。なお、画像形成条件の制御を、本明細書では、以降画像処理と記す。
ノズル特性の均一化は、加工精度やノズルの生産性との関係により制限される。このため、特にノズル特性に高い均一性が求められる場合、実際のノズル特性のばらつきを測定し、測定の結果得られたばらつきを画像処理によって補正することが考えられる。このような技術においては、ノズル特性のばらつきを正確に測定することが重要である。
【０００５】
ノズル特性を測定する従来技術としては、例えば、特許文献１、特許文献２が挙げられる。特許文献１に記載された発明は、ノズルからインクの液滴を噴出させ、この液滴が用紙等の記録媒体に付着する様子をカメラで記録する。そして、記録された映像を画像処理し、液滴が記録媒体に到達した位置とノズルの噴出口との位置関係を計測によって得ている。
【０００６】
また、特許文献２に記載された発明は、評価される画像のドットの重心位置を算出する。また、重心位置を中心にドットの配列周期よりも小さいマスク領域を設定し、マスク領域として設定された範囲の物理量の分布からドットごとの特徴量を算出している。
【特許文献１】特開２０００−６２１５８号公報
【特許文献２】特開２００１−２９１１０７号公報
【発明の開示】
【発明が解決しようとする課題】
【０００７】
しかしながら、特許文献１の発明は、液滴の撮影をするために特殊な撮影装置が必要になる。このため、ドットの位置検出にかかるコストが増加し、ひいてはプリンタヘッドの製造コストを高めることになる。
また、特許文献２の発明は、重心位置の算出にかかる計算量が多いという課題がある。すなわち、重心位置は、一般的に円形のテンプレートとのマッチングによって行われる。テンプレートマッチングは、ドットの大きさが分からないとき、さまざまな径及び大きさの円形テンプレートとドットとをマッチングする。このため、１つのドットの重心の算出に比較的時間がかかることになる。
【０００８】
さらに、近年ではヘッドが長大化する傾向があり、１インチ幅のヘッドに３６０個のノズルが設けられているものもある。このようなヘッドのノズル特性を測定する場合、各ノズルによって形成されたドットを１００個程度測定するとしても、３６０００個のドットを測定することになる。
このように多くのドットを測定する必要があるノズル特性の測定にあっては、ヘッド全体のノズル特性を測定するためには、１つのドットについての位置や位置ずれの検出にかかる時間は短いほど望ましい。
【０００９】
本発明は、このような点に鑑みてなされたものであり、特殊な計測装置を使うことがなく、ドットの位置及び位置ずれを短時間で、しかも高精度に測定できるドット位置検出装置及びドット位置検出方法をコンピュータに実行させるためのプログラム提供することを目的とする。
【課題を解決するための手段】
【００１０】
以上の課題を解決するため、本発明のドット位置検出装置は、インクジェットプリンタのドット位置を検出するドット位置検出装置であって、インクジェットプリンタで形成された画像におけるドット位置を予測するドット位置予測手段と、前記ドット位置予測手段によって予測されたドット位置に基づいて、前記画像におけるドットの位置を算出するドット位置算出手段と、前記ドット位置算出手段によって算出されたドットの位置に基づいて、ドットが形成されるべき理想位置を算出する理想位置算出手段と、前記ドット位置算出手段によって算出されたドットの位置と、前記理想位置算出手段によって算出されたドットの理想位置とのずれ量を算出するドット位置ずれ量算出手段と、を備えることを特徴とする。
【００１１】
このような発明によれば、演算処理を実行するプログラムのみでドットの位置及び位置ずれを検出することができる。また、この処理においてドットのパターンマッチング等の処理を用いることなくドットの位置及び位置ずれを検出することができる。したがって、本発明は、特殊な計測装置を使うことがなく、ドットの位置及び位置ずれを短時間で、しかも高精度に測定することができるドット位置検出装置を提供することができる。
【００１２】
また、本発明のドット位置検出装置は、前記ドット位置予測手段が画像中のドットとドットが無い領域とを輝度にかかる情報に基づいて判断すると共に、ドットとドットが無い領域とを判断するための輝度にかかるしきい値を設定するしきい値設定手段を備えることを特徴とする。
このような発明によれば、画像の濃度や特性に対応したしきい値を設定することができ、ドットの位置を正確に判断することができる。
【００１３】
また、本発明のドット位置検出装置は、前記ドット位置予測手段が、前記画像における第１の方向について輝度の値のヒストグラムを作成すると共に該ヒストグラムの極小点を通る第１の直線を求め、前記第１の方向と直交する第２の方向について輝度の値のヒストグラムを作成すると共に該ヒストグラムの極小点を通る第２の直線を求め、前記第１の直線と前記第２の直線との交点をドットの位置であると予測することを特徴とする。
【００１４】
このような発明によれば、比較的簡単にドットの位置を予測することができる。
また、本発明のドット位置検出装置は、前記ドット位置算出手段が、前記ドット位置予測手段によって予測されたドット位置に基づいて、前記画像におけるドットの半径をさらに算出することを特徴とする。
このような発明によれば、ドットにかかるより詳細なデータを取得することができ、ドットの位置をより正確に検出することができる。
【００１５】
また、本発明のドット位置検出装置は、前記ドット位置算出手段が、前記ドット位置予測手段によって予測されたドット位置を含む領域の輝度にかかるデータを互いに直交する二方向について取得することによってドットの範囲を計測し、計測された範囲からドットの重心点を算出することを特徴とする。
このような発明によれば、二方向からドットの位置を計測し、この結果を合わせてドットの重心を比較的簡易に算出することができる。
【００１６】
また、本発明のドット位置検出装置は、前記ドット位置算出手段が、予測されたドット位置を含む領域の輝度にかかるデータを画像に対する水平方向と垂直方向、水平方向と垂直方向とに対して傾きを持った二方向の少なくとも一方について取得することを特徴とする。
このような発明によれば、ドットの粗密等の状態や形状に応じてドットの位置をより正確に算出可能な方向を選択することができる。
【００１７】
また、本発明のドット位置検出装置は、前記ドット位置算出手段が、前記ドット位置予測手段によって予測されたドット位置に基づいて算出されたドットの位置に基づいて、さらに前記画像におけるドット位置を算出することを特徴とする。
このような発明によれば、算出されたドット位置の精度をより高め、ドット位置をより正確に算出することができる。
【００１８】
また、本発明のドット位置検出装置は、前記理想位置算出手段が、前記ドット位置算出手段によって算出されたドットの位置に基づいてドットを通る直線を近似して算出し、算出された直線の傾きの平均を傾きとする理想直線を算出し、算出された理想直線上をドットの理想位置とすることを特徴とする。
このような発明によれば、ドットの理想位置を比較的簡単、かつ短時間で算出することができる。
【００１９】
また、本発明のドット位置検出装置は、前記理想位置算出手段が、前記理想直線とドットの形成領域上で設定された軸との交点から理想直線の設定開始位置を求めることを特徴とする。
このような発明によれば、理想直線を自動的、かつ簡易に設定することができる。
また、本発明のドット位置検出方法をコンピュータに実行させるためのプログラムは、インクジェットプリンタのドット位置を検出するドット位置検出方法をコンピュータに実行させるためのプログラムであって、インクジェットプリンタで形成された画像におけるドット位置を予測するドット位置予測ステップと、前記ドット位置予測ステップにおいて予測されたドット位置に基づいて、前記画像におけるドットの位置を算出するドット位置算出ステップと、前記ドット位置算出ステップにおいて算出されたドットの位置に基づいて、ドットが形成されるべき理想位置を算出する理想位置算出ステップと、
前記ドット位置算出ステップにおいて算出されたドットの位置と、前記理想位置算出ステップにおいて算出されたドットの理想位置とのずれ量を算出するドット位置ずれ量算出ステップと、を含むことを特徴とする。
【００２０】
このような発明によれば、演算処理を実行するプログラムのみでドットの位置及び位置ずれを検出することができる。また、この処理においてドットのパターンマッチング等の処理を用いることなくドットの位置及び位置ずれを検出することができる。したがって、本発明は、特殊な計測装置を使うことがなく、ドットの位置及び位置ずれを短時間で、しかも高精度に測定することができるドット位置検出方法をコンピュータに実行させるためのプログラムを提供することができる。
【００２１】
また、本発明のドット位置検出方法をコンピュータに実行させるためのプログラムは、前記ドット位置予測ステップは、前記画像における第１の方向について輝度の値のヒストグラムを作成すると共に該ヒストグラムの極小点を通る第１の直線を求めるステップと、前記第１の方向と直交する第２の方向について輝度の値のヒストグラムを作成すると共に該ヒストグラムの極小点を通る第２の直線を求めるステップと、を含み、前記第１の直線と前記第２の直線との交点をドットの位置であると予測することを特徴とする。
【００２２】
このような発明によれば、比較的簡単にドットの位置を予測することができる。
また、本発明のドット位置検出方法をコンピュータに実行させるためのプログラムは、前記理想位置算出ステップが、前記ドット位置算出ステップにおいて算出されたドットの位置に基づいてドットを通る直線を近似して算出するステップと、算出された直線の傾きの平均を傾きとする理想直線を算出するステップと、を含み、算出された理想直線上をドットの理想位置とすることを特徴とする。
【００２３】
このような発明によれば、ドットの理想位置を比較的簡単、かつ短時間で算出することができる。
【発明を実施するための最良の形態】
【００２４】
以下、図を参照して本発明に係るドット位置検出装置及びドット位置検出方法及をコンピュータに実行させるためのプログラムの実施形態を説明する。図１は、本発明の一実施形態のドット位置検出装置及びドット位置検出方法をコンピュータに実行させるためのプログラムを説明するための図である。ドット位置検出方法をコンピュータに実行させるためのプログラムは、画像の解析アプリケーション２００として構成される。このため、本実施形態では、汎用的なＰＣ（Personal Computer）に解析アプリケーション２００を格納し、ドット位置検出装置として機能させるものとする。
【００２５】
解析アプリケーション２００は、解析すべき画像をスキャナ１０１で取り込む解析画像入力部２０１、ドット領域と背景領域とを判定するための輝度のしきい値を設定するしきい値設定部２０３、ドットの位置を検出するドット解析部２０５、ドット解析部２０５で解析された結果を出力する結果出力部２０７を備えている。本実施形態のドット位置検出方法は、以上の機能のうち、特にドット解析部２０５におけるドットの解析方法に特徴を有するものである。
【００２６】
しきい値設定部２０３におけるしきい値は、解析アプリケーション２００が平均分析２値化値法などの一般的なアルゴリズムを用いて設定するものであってもよい。また、オペレータが２値化しきい値を設定することも可能である。このとき、しきい値の設定は、ＰＣが一般的に備えているキーボードなどの入力装置２０９によって行われる。また、ドット解析の結果出力部２０７は、ディスプレイ等の出力装置２１１へ結果を表示データやファイルとして出力する。
【００２７】
図２は、図１に示した解析アプリケーション２００を動作させるハードウェア構成を説明するための図である。図２（ａ）は、汎用ＰＣのハードウェア構成を説明するための図である。図示した汎用ＰＣは、スキャナ１０１で読み取られた画像データ等を外部から入力するＵＳＢ１０３、入力された画像データを処理するためのプログラムやデータを記憶するメモリ１０７、メモリ１０７に記憶されているプログラムやデータを使って動作するハードディスクと読取装置とが一体化して構成されるハードディスク・ドライバ（ＨＤＤ）１０９、以上の構成を統括的に制御するＣＰＵ１０５を備えている。
【００２８】
また、図２（ａ）に示したＰＣは、オペレータの操作画面や処理結果を表示する表示画面として機能するディスプレイ１１５、キーボードやタッチパネル等の入力装置２０９を備えている。さらに、ソフトドライバをインストールするためのＣＤを駆動するＣＤドライバやＤＶＤを駆動するためのＤＶＤドライバ１１３を備える場合もある。
なお、本実施形態のドット位置検出装置は、汎用ＰＣとして構成されるものに限定されるものではない。例えば、図２（ｂ）のように、プリンタとして構成してもよい。図２（ｂ）に示したプリンタは、ハードウェア構成として、ＰＣと同様に、ＵＳＢ１０３、ＣＰＵ１０５、メモリ１０７を備えている。また、プリンタは、インクジェット方式であればインクを噴出するノズルが多数設けられたノズルヘッド１１９、用紙を搬送する機械的構成である搬送メカ部１２１、プリンタ動作に必要なプログラム等を記憶するためのＲＯＭ１２３を備えている。
【００２９】
解析アプリケーション２００は、図２（ａ）に示したＰＣにあってはＨＤＤ１０９のハードディスクにソフトウェアとして格納されるようにしてもよい。また、図２（ｂ）に示したプリンタにあっては、ＲＯＭ１２３にソフトウェアとして格納されるようにしてもよい。いずれの場合にあっても、解析アプリケーション２００は、ＣＰＵ１０５にロードされて実行される。また、ドット位置検出装置をＰＣとして構成した場合、解析される画像データは、スキャナ１０１によって読み込まれるものの他、ＨＤＤ１０９、ＣＤ１１１、ＤＶＤ１１３などの記録メディアに保存されたものを使うことも可能である。
【００３０】
図３は、図２に示したドット解析部２０５の機能ブロック図である。また、図４は、図３に示したドット解析部２０５でなされる処理の概略を示すフローチャートである。図３に示すように、ドット解析部２０５は、ドット位置予測部３０５、ドット重心位置計算部３０７、計測したドット位置から理想着弾位置を推測する理想位置計算部３０９を備えている。なお、本実施形態では、ドットが円形であり、ドットの重心は円の中心に一致するが、文言を重心に統一するものとする。また、本明細書では、ノズルから噴出されたインクの液滴が用紙等の媒体に達することを着弾という。また、液滴が用紙に最初に着弾した位置を着弾位置というものとする。理想着弾位置とは、画像品質の観点から最も適切であると考えられる着弾位置をいう。
【００３１】
さらに、ドット解析部２０５は、位置ズレ計算部３１１を備える。位置ズレ計算部３１１の計算の結果得られた位置ずれのデータは、出力部３１３を介して結果出力部２０７に出力される。
ドット解析部２０５には、解析画像入力部２０１から解析画像を取り込むと共に、しきい値設定部２０３が設定したしきい値が取り込まれる。そして、図４に示すように、ドット位置予測部３０５が画像データからドットのおおよその位置を予測する（Ｓ４０１）。また、ドット重心位置計算部３０７は、しきい値を使ってドットの重心位置とドットの径を計測する（Ｓ４０２）。理想位置計算部３０９は、ドット重心位置計算部３０７によって計測されたドット位置から理想着弾位置を推測する（Ｓ４０３）。
【００３２】
さらに、位置ズレ計算部３１１は、実際のドットの位置と理想着弾位置とのずれを計算する（Ｓ４０４）。なお、本実施形態では、ノズルから噴出されたインクの液滴が用紙等の媒体に達することを着弾という。また、液滴が用紙に最初に着弾した位置を着弾位置というものとする。理想着弾位置とは、画像品質の観点から最も適切であると考えられる着弾位置をいう。
【００３３】
位置ズレ計算部３１１の計算の結果得られた位置ずれのデータは、出力部３１３を介して結果出力部２０７に出力される。
次に、以上述べたドット解析部２０５によってなされる処理を、ドット中心位置予測、ドット重心位置計算、理想位置計算、位置ずれ計算の順により詳細に説明する。
（１）ドット中心位置予測
図５は、本実施形態のドット中心予測の処理を説明するための図である。図５に示すように、ドット位置予測部３０５は、入力された画像データの輝度を、列方向（図中にｘで示す）に積算する。また、同様に、輝度を行方向（図中にｙ方向で示す）に積算する。積算された輝度を、図中に積算輝度として示す。
【００３４】
本実施形態では、画像データが、ドットの輝度が背景の輝度より低い値を持つものとする。このため、積算輝度は、ｘ方向、ｙ方向共にドットの粗密の程度によって極大点、極小点を持つ二次曲線を示す。二次曲線の極大点は輝度が高い、つまりドットのない背景の領域を示し、極小点はドットが密集していることを示す。
ドット位置予測部３０５は、ｘ方向及びｙ方向の積算輝度から極小点の位置を検出する。そして、この位置を通る直線（図中にｌｘ、ｌｙとして示す）を列方向、行方向の各々について仮想的に設定し、直線ｌｘとｌｙとの交点をドット位置であると予測する。本実施形態では、以降、予測されたドット位置の座標を（ｘ1，ｙ1）、（ｘ2，ｙ2）、…（ｘn，ｙm）と記す。
【００３５】
なお、しきい値設定部２０３で設定されたしきい値は、極小値の検出に使用することもできる。つまり、本実施形態では、二次曲線上で下に凸の値を有し、かつ、図中に示すしきい値ｔｈよりも小さい値を持つ点を極小点とする。
しきい値ｔｈを設定したことにより、本実施形態のドット位置検出装置は、極小点の値の検出に対する外乱等の影響を抑え、極小点の位置、つまりドット位置をより正確に検出することができる。
図６は、以上述べた処理の手順をより具体的に説明するためのフローチャートである。図示するように、ドット位置予測部３０５は、積算輝度値をｘ方向、ｙ方向の各々について求める。具体的には、画像データの列方向（ｘ方向）、行方向（ｙ方向）の各々について、輝度のヒストグラムを作成する（Ｓ６０１）。続いて、列方向のヒストグラムから検出された最小点の座標（ｘ方向のみの座標）を全てファイルＰｘ[Ｗ]に保存する。また、同様に、行方向のヒストグラムから検出された最小点の座標（ｙ方向のみの座標）を全てファイルＰｙ[Ｈ]に保存する（Ｓ６０３）。
【００３６】
また、ドット位置予測部３０５は、予測された中心位置からドット数をも予測することができる。列方向、行方向それぞれの極小点の数を求め、列方向の極小点の個数をｎ個、行方向の極小点の個数をｍ個として保存する（Ｓ６０４）。そして、列方向、行方向の各々について、極小点間の距離の平均値ｐ（平均距離ｐ）を計算し、保存する。このとき、列方向の平均距離と行方向の平均距離が異なる場合、より小さい方の距離を平均距離ｐとして採用する。
【００３７】
平均距離ｐの算出は、後のドット重心位置計算で取得する輝度データ数を求めるために行うものである。すなわち、ドット重心位置計算では、行、列方向に配置されている複数のドットの個々の輝度を検出して各ドットの重心を計算している。このため、隣接する２つのドットに渡って検出された輝度データを１つのドットのものとした場合、重心の正確な位置を検出することができなくなる。このため、本実施形態では、隣接するドット間の平均距離を予め算出しておき、あるドットの輝度データの取得にあたって隣接するドットの輝度データを取得することを防いでいる。
【００３８】
なお、ドットの大きさや間隔がある程度分かっている場合、重心位置判定での処理量を考慮しない場合に、平均距離ｐを固定値とすることも可能である。また、平均距離に基づいて決定される輝度データの取得数は、奇数であることが後の重心計算処理の計算上望ましい。そこで、本実施形態では、データ取得数が偶数の場合、１を加算する等の処理を行っている。
【００３９】
以上の処理により、入力された画像データによって形成される画像のドット位置を予測する予測座標（ｘ1，ｙ1）、（ｘ2，ｙ2）、…（ｘn，ｙm）の予測が完了する。
（２）ドット重心位置計算
次に、本実施形態のドット位置検出装置のドット重心位置計算部３０７で行われる処理を説明する。ドット重心位置計算部３０７は、予想された座標を通り、かつ直交する２直線に沿ってドットの濃度を検出することによってドットの外縁を判定する。この際、ドット重心位置計算部３０７は、画像の特性等によって直交する２直線の傾きを任意に設定することができる。
【００４０】
図７（ａ）、（ｂ）は、いずれもドットの重心の位置を計算する処理を説明するための図である。図７（ａ）は、直交する２直線が各々垂直、水平方向に延びる例を示し、（ｂ）は、直交する２直線が、（ａ）に示した２直線に対して４５度の傾きを持つ例を示している。図７（ａ）、（ｂ）に示した例は、いずれの場合にも同様にして重心を算出することが可能であるが、本実施形態では（ｂ）に示した例について重心の位置を計算する処理を説明するものとする。
【００４１】
図７（ｂ）は、ドットｄ上に予測された予測座標（Ｘｉ，Ｙｉ）があり、予測座標（Ｘｉ，Ｙｉ）とずれた位置にドットｄの重心（Ｘｏ，Ｙｏ）があることを示している。ドット重心位置計算部３０７は、直線ｌｙに沿ってデータを取得する。データの取得は、画素単位で行われる。図中にデータの取得点をｑ、データ取得点ｑに対応する測定点ｑ’を示す。
【００４２】
ドット重心位置計算部３０７は、直線ｌｙに沿って画像の輝度データをｒ個取得する。そして、取得されたｒ個のデータのうち、輝度のしきい値ｔｈ以下の輝度を示す輝度データを抽出し、抽出された輝度データがドット上で取得されたものとみなす。つまり、直線ｌｙに沿う輝度データの取得において、しきい値ｔｈ以下の輝度データの取得が開始される点（ｘ０，ｙ０）の測定後、しきい値ｔｈ以下の輝度データが取得されなくなる点（ｘ１，ｙ１）までの測定距離がドットの外縁から外縁に相当する。
【００４３】
本実施形態では、ドット重心位置計算部３０７が、以下の式を使ってドットｄ上の直線ｌｙに沿う弦の距離を測定し、この弦の中心座標（ｄｘ１，ｄｙ１）と弦の半分の長さＲ１とを求める。
ｄｘ１＝（ｘ１−ｘ０）／２                …（式１）
ｄｙ１＝（ｙ１−ｙ０）／２                …（式２）
Ｒ１＝√２（（ｘ１−ｘ０）＋（ｙ１−ｙ０））／４      …（式３）
また、ドット重心位置計算部３０７は、直線ｌｘに沿って画像の輝度データをｒ個取得する。そして、取得されたｒ個のデータのうち、輝度のしきい値ｔｈ以下の輝度を示す輝度データを抽出し、抽出された輝度データがドット上で取得されたものとみなす。つまり、直線ｌｙに沿う輝度データの取得において、しきい値ｔｈ以下の輝度データの取得が開始される点（ｘ０，ｙ０）の測定後、しきい値ｔｈ以下の輝度データが取得されなくなる点（ｘ１，ｙ１）までの測定距離がドットの外縁から外縁に相当する。
【００４４】
本実施形態では、ドット重心位置計算部３０７が、以下の式を使ってドットｄ上の直線ｌｙに沿う弦の距離を測定し、この弦の中心座標（ｄｘ２，ｄｙ２）と弦の半分の長さＲ２とを求める。
ｄｘ２＝（ｘ１−ｘ０）／２                …（式４）
ｄｙ２＝（ｙ０−ｙ１）／２                …（式５）
Ｒ２＝√２（（ｘ１−ｘ０）＋（ｙ０−ｙ１））／４      …（式６）
以上の計算によって得られた値を用い、ドットｄの重心（Ｘｏ，Ｙｏ）及びドットｄの半径（Ｒｏ）は、以下の式により求められる。
【００４５】
Ｘｏ＝Ｘｉ＋ｄｘ１＋ｄｘ２
Ｙｏ＝Ｙｉ＋ｄｙ１＋ｄｙ２
Ｒｏ＝（Ｒ１＋Ｒ２）／２
以上述べた演算は、複数回行うことによってさらにドットの重心の精度を高めることができる。すなわち、ドット重心位置計算部３０７は、ドット位置予測部３０５によって予測されたドット位置に基づいてドットの位置を算出する。そして、算出されたドットの位置を予測されたドットの位置に置き換えてドット位置を算出してもよい。また、以上の処理においては、ドット間のデータを線形補完などの手法で補完すると重心の計算精度を高めることができる。
【００４６】
なお、上記したｄｘ１，ｄｙ１及びｄｘ２，ｄｙ２と図７（ａ）、（ｂ）中に示したＸｏ，Ｙｏ及びＸｉ，Ｙｉの関係を図７（ｃ）として示す。
ドット重心位置計算部３０７は、以上述べた処理で得られた得られるｄｘ、ｄｙ、ドットｄの重心座標、ドット半径Ｒをｎ×ｍ個のドットの全てについて算出し、結果を保存する。結果の保存のため、ドット重心位置計算部３０７は、情報バッファＰｔ［ｎ］［ｍ］を作成する。情報バッファＰｔ［ｎ］［ｍ］とは、各ドットについて、少なくとも以下の情報を格納するバッファである。
【００４７】
重心の座標（ｘ，ｙ）    Ｐｔ［Ｉ］［ｊ］．ｘ、Ｐｔ［Ｉ］［ｊ］．ｙ
ドット半径  Ｒ          Ｐｔ［Ｉ］［ｊ］．Ｒ
位置ズレ量（ｄｘ、ｄｙ）    Ｐｔ［Ｉ］［ｊ］．ｄｘ、Ｐｔ［Ｉ］［ｊ］．ｄｙ
図８は、以上述べたドット重心位置計算部３０７でなされる処理を示したフローチャートである。図８のフローチャートによれば、以上の処理は、具体的には以下のように行われる。すなわち、ドット重心位置計算部３０７は、先ず、情報バッファＰｔ［ｎ］［ｍ］を作成する（Ｓ８０１）。そして、座標（ｘ，ｙ）を設定し、カウントｘ及びカウントｙの値を０に初期化する（Ｓ８０２）。
【００４８】
そして、ドット重心位置計算部３０７は、この座標（ｘ，ｙ）にｍ×ｎ個のドットの予測位置（Ｘｉ，Ｙｉ）を順次代入し、代入された予測位置で各ドットの重心位置（Ｘｏ，Ｙｏ）及び半径Ｒｏを算出する（Ｓ８０３）。本実施形態では、ステップＳ８０３の処理をサブルーチンとして図８に示し、サブルーチンの内容を、図９を使って説明する。
次に、ドット重心位置計算部３０７は、情報バッファＰｔ［ｎ］［ｍ］に得られた重心の座標と半径とを保存する（Ｓ８０４）。さらに、ステップ８０２で０に設定されたカウントｘをカウントアップし（Ｓ８０５）、カウントｘがｎに達したか否か判断する（Ｓ８０６）。この判断の結果、カウントｘがｎに達していない場合（Ｓ８０６：Ｎｏ）、ドット重心位置計算部３０７は、カウントアップされたｘの位置にあるドットの重心の計算をする。
【００４９】
また、ドット重心位置計算部３０７は、カウントｘのカウント値がｘに達した場合（Ｓ８０６：Ｙｅｓ）、カウントｘを０にすると共にカウントｙの値をカウントアップする（Ｓ８０７）。そして、カウントアップされたカウントｙの値がｍに達したか否か判断し（Ｓ８０８）、カウントｙがｍに達していない場合には引き続きドットの重心の計算をする。また、カウントｙがｍに達した場合（Ｓ８０８：Ｙｅｓ）、処理を終了する。
【００５０】
図９は、図８に示したステップＳ８０３の処理を説明するためのフローチャートである。図９に示した処理は、予測座標（Ｘｉ，Ｙｉ）を入力し、ドットの重心座標（Ｘｏ，Ｙｏ）とドット半径Ｒｏとを出力する。なお、図９に示したフローチャートは、ドット１つにつき行われる処理であって、１つのドットの重心座標（Ｘｏ，Ｙｏ）及びドット半径Ｒｏが算出されるたびに終了し、次のドットの予測座標（Ｘｉ，Ｙｉ）が入力されるたびに開始される。
【００５１】
ステップＳ８０３では、ドット重心位置計算部３０７が、先ず、カウンタ（Ｃｎｔ）を０にリセットする（Ｓ９０１）。次に、カウンタの値が０、１、２のいずれであるか判断し（Ｓ９０２）、この値が０である場合には、推測されたドットの位置を示す２直線のうちの一方に沿ってｒ個の輝度のデータを取得する（Ｓ９０３）。
次に、ドット重心位置計算部３０７は、取得されたデータのうち、輝度がしきい値ｔｈ以下のデータをドット上で取得されたデータとみなし、しきい値ｔｈ以下の値が取得される距離を測定する。そして、測定された距離を持つ弦の中心位置（Ｘ１，Ｙ１）及び弦の半分の長さＲ１を算出する（Ｓ９０４）。そして、後にドットの重心位置の算出に使用されるｄｘ１、ｄｙ１の値を算出する（Ｓ９０５）。
【００５２】
以上の処理の後、ドット重心位置計算部３０７は、カウンタをカウントアップし（Ｓ９０６）、ステップＳ９０２に戻ってカウンタのカウント値を判断する。そして、カウンタ値が１であった場合、ステップＳ９０３でデータを取得した直線と異なる他の直線に沿ってｒ個の輝度のデータを取得する（Ｓ９０７）。
次に、ドット重心位置計算部３０７は、取得されたデータのうち、輝度がしきい値ｔｈ以下のデータをドット上で取得されたデータとみなし、しきい値ｔｈ以下の値が取得される距離を測定する。そして、測定された距離を持つ弦の中心位置（Ｘ２，Ｙ２）及び弦の半分の長さＲ２を算出する（Ｓ９０８）。そして、後にドットの重心位置の算出に使用されるｄｘ２、ｄｙ２の値を算出する（Ｓ９０９）。
【００５３】
以上の処理の後、ドット重心位置計算部３０７は、カウンタをカウントアップし（Ｓ９０６）、ステップＳ９０２に戻ってカウンタのカウント値を判断する。そして、カウンタ値が２であった場合、ステップＳ９０６で算出されたｄｘ１、ｄｙ１、ステップＳ９０９で算出されたｄｘ２、ｄｙ２を使ってドットの重心座標（Ｘｏ，Ｙｏ）とドット半径Ｒｏとを算出する（Ｓ９１０）。
【００５４】
以上の処理により、本実施形態のドット位置検出装置は、ドット重心位置計算部３０７によって実際のドットの重心の位置を示す座標（Ｘｏ，Ｙｏ）と、半径Ｒｏとを算出することができる。
（３）理想位置計算
次に、本実施形態のドット位置検出装置の理想位置計算部３０９で行われる処理を説明する。図１０は、本実施形態の理想ドット位置の表し方を説明するための図である。本実施形態では、図示するように、ドットの配置を、ドットを通って、かつ直交する２方向の直線群で示す。本実施形態では、この直線群を各々理想直線と記す。図１０に示した理想直線同士の交点は、噴出されたインクの理想的な着弾位置を示す。ドットの理想的な配置とは、ドットが極力等間隔で、かつ理想直線上に配列された状態をいう。また、理想直線は、各々一次式の傾きや切片を使って以下のように表される。
【００５５】
ｋ行目の理想直線
Ｙｉ＝Ａｘ・Ｘ+Ｂｘ（ｋ）
ただし、Ｂｘ（ｋ）＝Ｙ０+Ｐｘ・ｋ
ｋ列目の理想直線
Ｘｉ＝Ａｙ・Ｙ+Ｂｙ（ｋ）
ただし、Ｂｙ（ｋ）＝Ｘ０+Ｐｙ・ｋ
図１１は、ドットを通る理想直線の傾きと切片との求め方を説明するための図である。本実施形態では、実際に配置されたドットの各行、各列に基づき、各ドット中心を通る直線を最少二乗法等を使った近似式によって求める。求められた直線を、図中に直線Ｘ、直線Ｙとして示す。なお、直線Ｘは行方向に延びる直線を示し、直線Ｙは、列方向に延びる直線を示す。理想位置計算部３０９は、得られた直線Ｘ、直線Ｙの各々について、傾きと切片とを算出する。そして、算出された複数の直線の傾きの平均値（Ａｘ、Ａｙ）と、切片間の距離についての平均値（Ｐｘ、Ｐｙ）を算出する。
【００５６】
理想直線は、算出された傾きの平均値を傾きとし、理想直線の番号ｋと切片間の距離の平均値の積と切片初期位置Y0の和（Ｂｘ（ｋ）、Ｂｙ（ｋ））を切片とする直線であって、直線Ｘ、直線Ｙの各々について定められる。直線Ｘの理想直線を理想直線Ｘｉ、直線Ｙの理想直線を理想直線Ｙｉと以降記すものとする。理想直線Ｘｉ、理想直線Ｙｉは、ドットの各行、各列に設定され、後でドットとのずれ量の算出に使用される。また、ｘ方向に配置された直線間の距離の平均を水平方向ドットピッチとし、ｙ方向に配置された直線間の距離の平均を水直方向ドットピッチとする。
【００５７】
次に、理想位置計算部３０９は、理想直線のうち、垂直方向に配列される直線Ｙ、垂平方向に配列される直線Ｘの基準点を決定する。すなわち、理想位置計算部３０９は、ドットが配置される配置領域１２００を規定する。そして、配置領域１２００における左端の縦軸と直線Ｙとが交差する点である初期切片Ｙ０を算出する。また、理想位置計算部３０９は、配置領域１２００における左端の縦軸と直線Ｘとが交差する点である初期切片Ｘ０を算出する。
【００５８】
図１２は、初期切片Ｙ０の算出方法を示す図であって、図１３は初期切片Ｘ０の算出方法を示す図である。初期切片Ｙ０は、理想直線Ｙｉ交差するドットの重心の座標と理想直線Ｙｉとの線形最少二乗誤差が最少の値をとるように算出される。また、初期切片Ｘ０は、理想直線Ｘｉと交差するドットの重心の座標と理想直線Ｘｉとの線形最少二乗誤差が最少の値をとるように算出される。
【００５９】
理想直線Ｙｉとの線形最少二乗誤差は、以下の式によって算出できる。
Ｌ行目の理想直線Ｙｉ
Ｙ＝Ａｘ_ave・Ｘ＋Ｂｘ（Ｌ）
（但し、Ｂｘ（Ｌ）＝Ｙ０＋Ｂｘ_ave・Ｌ）
この直線と、ドット中心座標（ｘ，ｙ）の誤差：Ｅｒｒは、
Ｅｒｒ＝Ａｘ_ave・ｘ＋Ｂｘ（Ｌ）−ｙ
また、理想直線Ｘｉとの線形最少二乗誤差は、以下の式によって算出できる。
【００６０】
Ｌ列目の理想直線Ｙｉ
Ｘ＝Ａｙ_ave・Ｙ＋Ｂｙ（Ｌ）
（但し、Ｂｙ（Ｌ）＝Ｘ０＋Ｂｙ_ave・Ｌ）
この直線と、ドット中心座標（ｘ，ｙ）の誤差：Ｅｒｒは、
Ｅｒｒ＝Ａｙ_ave・ｙ＋Ｂｙ（Ｌ）−ｘ
本実施形態では、以上の式によってｎ×ｍ個の全ドットの２乗誤差の総和（Ｓｕｍ）を最小にする初期切片Ｘ０、Ｙ０を求めている。
【００６１】
図１４は、以上述べた理想位置計算部３０９でなされる処理全般を説明するためのフローチャートである。図１４に示したように、理想位置計算部３０９は、先ず、情報バッファＰｔ［ｎ］［ｍ］に保存されているドットの重心位置や半径等の座標からドットの配置を検出し、このドットの配置に基づいて、ｎ×ｍ個のドットを通る複数の直線の傾き及び切片を算出する（Ｓ１４０１）。
【００６２】
次に、理想位置計算部３０９は、複数の直線の傾き及び切片から、行方向に配置される直線Ｙ、列方向に配置される直線Ｘについて平均値を算出する。つまり、直線Ｙの傾きの平均値（Ａｘ_ave）と、各切片間の距離の平均値（Ｂｘ_ave）を求める。また、直線Ｘの傾きの平均値（Ａｙ_ave）と、各切片間の距離の平均値（Ｂｙ_ave）を求める。（Ｓ１４０２）。
【００６３】
次に、理想位置計算部３０９は、ドットが配置される領域基準となる縦軸と理想直線Ｘｉとが交差する点を示す初期切片を算出する（Ｓ１４０３）。また、理想位置計算部３０９は、ドットが配置される領域基準となる縦軸と理想直線Ｙｉとが交差する点を示す初期切片を算出する（Ｓ１４０４）。本実施形態のドット位置検出装置は、初期切片を算出することによってドットの形成開始位置を定めることができる。なお、初期切片の算出は、線形最少誤差を算出し、この値が最小値をとるように定められる。
【００６４】
図１５は、図１４に示したステップＳ１４０１の処理を説明するためのフローチャートである。先ず、理想位置計算部３０９は、列をカウントするカウンタｘを０にセットする（Ｓ１５０５）。そして、情報バッファＰｔ［ｎ］［ｍ］から１列分のｍ個のドットのデータを取得する（Ｓ１５０６）。さらに、各ドットの重心の座標から線形最少二乗法を使って各ドットの重心を通るものと近似される直線の傾き及び切片を算出する（Ｓ１５０７）。
【００６５】
次に、理想位置計算部３０９は、カウンタｘを１つカウントアップし（Ｓ１５０８）、カウンタｘがｎに達したか否か判断する（Ｓ１５０９）。この判断の結果、カウンタｘがｎに達していない場合には（Ｓ１５０９：Ｎｏ）、次の列のドットをｍ個取得する。
また、理想位置計算部３０９は、カウンタｘがｎに達している場合（Ｓ１５０９：Ｙｅｓ）、行方向に配置されたドットの処理に移行する。そして、行をカウントするカウンタｙを０にセットする（Ｓ１５１０）。さらに、情報バッファＰｔ［ｎ］［ｍ］から１行分のｎ個のドットのデータを取得する（Ｓ１５１１）。そして、各ドットの重心の座標から線形最少二乗法を使って各ドットの重心を通るものと近似される直線の傾き及び切片を算出する（Ｓ１５１２）。
【００６６】
次に、理想位置計算部３０９は、カウンタｙを１つカウントアップし（Ｓ１５１３）、カウンタｙがｍに達したか否か判断する（Ｓ１５１４）。この判断の結果、カウンタｘがｎに達していない場合には（Ｓ１５１４：Ｎｏ）、次の列のドットをｎ個取得する。また、カウンタｙがｍに達している場合（Ｓ１５１４：Ｙｅｓ）、処理を終了する。
図１６及び１７は、図１４に示したステップＳ１４０３の処理を説明するためのフローチャートであって、ドットの水平方向について行われる処理を示している。本実施形態では、図示を省くが、ドットの垂直方向について行われるＳ１４０４の処理のフローチャートは、図１６及び１７のフローチャート中に記したＹ０→Ｘ０、ｄｙ→ｄｘに変更することによって得ることができる。
【００６７】
理想位置計算部３０９は、先ず、座標Ｙ０、線形二乗誤差の総和Ｓｕｍ、前回算出された線形二乗誤差の総和ＰｒｅＳｕｍの値をいずれも０にセットする（Ｓ１６０１）。なお、Ｙ０は、理想直線Ｙｉの初期切片である。
次に、理想位置計算部３０９は、線形二乗誤差の総和Ｓｕｍを算出する（Ｓ１６０２）。この算出については、図１７を用いて後に説明する。総和Ｓｕｍの算出後、理想位置計算部３０９は、座標Ｙ０を所定の増分値ｄｙだけ増加させる（Ｓ１６０３）。また、算出された総和Ｓｕｍを総和ＰｒｅＳｕｍとする（Ｓ１６０６）。そして、増加された新たなＹ０を使って線形二乗誤差の総和Ｓｕｍを算出する（Ｓ１６０５）。さらに、ステップＳ１６０５で算出された総和Ｓｕｍと先に算出された総和Ｓｕｍとの差分ｄｓを求め（Ｓ１６０６）、差分ｄｓの絶対値を予め設定されたしきい値ＴＨと比較する（Ｓ１６０７）。
【００６８】
なお、差分ｄｓの絶対値と比較されるしきい値ＴＨは、二回の計算で算出された総和Ｓｕｍの違いとして許容できる値として設定されたしきい値である。しきい値ＴＨとしては、例えば、０．０１が設定される。また、増分値ｄｙは、１００程度に設定される。
差分ｄｓがしきい値ＴＨよりも小さい場合（Ｓ１６０７：Ｙｅｓ）、理想直線Ｙｉの初期切片が適正な値であるとして処理を終了する。また、差分ｄｓがしきい値ＴＨよりも大きい場合（Ｓ１６０７：Ｎｏ）、差分ｄｓが負であるか否か判断する（Ｓ１６０８）。差分ｄｓが負である場合（Ｓ１６０８：Ｙｅｓ）、再度初期切片Ｙ０を増加させて総和Ｓｕｍを算出する。また、差分ｄｓが正である場合（Ｓ１６０８：Ｎｏ）、増加分ｄｙに−０．８程度を乗じることによってｄｙの増加方向を変更し（Ｓ１６０９）、再度初期切片Ｙ０を増加させて総和Ｓｕｍを算出する。
【００６９】
図１７は、図１６に示したステップＳ１６０２の処理を説明するためのフローチャートである。図１６に示したフローチャートの処理では、理想位置計算部３０９が、線形二乗誤差の総和Ｓｕｍ、カウンタｘ、カウンタｙを各々０に設定する（Ｓ１７０１）。そして、情報バッファＰｔ［ｎ］［ｍ］に保存されているドットの重心の座標（Ｘｏ，Ｙｏ）と、ｘ列目の理想直線（ｘ列目のドットによって形成され得る直線の平均の傾きと切片とを持つ直線）との誤差Ｅｒｒを算出する（Ｓ１７０２）。
【００７０】
なお、ｘ列目の理想直線Ｙxは、
ｙ＝Ａｘ_ave・x＋Ｂｘ
ただし、Ｂｘ＝Ｙｏ＋Ｂx_ave・ｘ
ｘ列目の理想直線とドットの重心の座標（Ｘｏ，Ｙｏ）との誤差Ｅｒｒは、以下のように表される。
【００７１】
Ｅｒｒ＝Ａｘ_ave・Ｘｏ＋Ｂｘ−Ｙｏ
次に、理想位置計算部３０９は、カウンタｘを１つカウントアップする（Ｓ１７０４）。そして、カウントアップ後のカウンタｘがｎに達したか否か判断し（Ｓ１７０５）、達していない場合には次のドットについての処理を実行する（Ｓ１７０５：Ｎｏ）。また、カウンタｘがｎに達した場合には（Ｓ１７０５：Ｙｅｓ）、カウンタｙをカウントアップする（Ｓ１７０６）。そして、カウンタｙがｍに達したか否か判断する（Ｓ１７０７）。
【００７２】
（４）位置ずれ計算
次に、本実施形態のドット位置検出装置の位置ずれ計算部３１１で行われる処理を説明する。図１８（ａ）、（ｂ）、（ｃ）は、位置ずれ計算の手法を説明するための図であって、（ａ）は実際のドットの配置を示す。また、図１８（ｂ）は、（ａ）に示したドット上に破線で示す理想直線を重ねて示した図であって、理想直線の交点が理想的な着弾位置を示す。また、図１８（ｃ）は、本実施形態でいうドットの位置ずれの定義を説明するための図である。
【００７３】
図示したように、ドットは、Ｘ方向に６個、Ｙ方向に５個配置された６×５のマトリックスを構成して配置されている。また、ドットの配置領域は、縦に長さＹを有し、横に長さＨを有するものとする。図１８（ａ）のように配置されたドットに対し、理想直線Ｘｉ、理想直線Ｙｉを重ねると、図１８に示すドットの位置ずれ量が算出できる。
図１８（ｃ）に示したように、ドットｄは、重心点Ｐｒを中心にした半径Ｒの円として形成される。理想直線Ｘｉ、理想直線Ｙｉは、ドットｄ上の点Ｐｉで交差しており、重心点Ｐｒと点Ｐｉとの相違のｘ成分をｄｘ、ｙ成分をｄｙと記す。位置ずれ計算部３１１は、理想直線Ｘｉ、理想直線Ｙｉの交点と重心点とのズレ量を全ドットについて計算する。この計算は、以下のようにして行われる。
【００７４】
ｉ列ｊ行目の理想着弾位置（理想直線Ｘｉｉ、理想直線Ｙｉｊの交点（ｘ，ｙ））
ｙ＝Ａｘ_ave・x＋Ｂｘ（ｉ）
（Ｂｘ（ｉ）＝Ｙ０＋Ｂx_ave・i）
ｘ＝Ａｙ_ave・y＋Ｂｙｊ
（Ｂｙ（ｊ）＝Ｘ０＋Ｂｙ_ave・j）より、
ｘ＝（Ａｙ_ave・Ｂｘ（ｉ）＋Ｂｙ（ｊ））／（１−Ａｙ_ave・Ａｘ_ave）
したがって、位置ズレ量ｄｘ，ｄｙは、以下のように求まる。
【００７５】
Ｐｔ［ｉ］［ｊ］．ｄｘ＝ｘ−Ｐｔ［ｉ］［ｊ］．ｘ
Ｐｔ［ｉ］［ｊ］．ｄｙ＝ｙ−Ｐｔ［ｉ］［ｊ］．ｙ
以上の計算式によって算出されたずれ量は、例えば、以下の形式の出力データとして出力される。
ドット位置（ｍ，ｎ）：（ずれ量（ｄｘ，ｄｙ），ドット径
（０，０）：（０，１），１２
（０，１）：（０，２），１２
（０，２）：（０，１），１２
（０，３）：（０，１），１２
・
・
・
（ｎ，ｍ）：（−１，０），１２
図１９は、位置ずれ計算部３１１で行われる処理を説明するためのフローチャートである。位置ずれ計算部３１１は、行をカウントするカウンタｘ、列をカウントするカウンタｙをそれぞれ０に設定する（Ｓ１９０１）。そして、情報バッファＰｔ［ｎ］［ｍ］に保存されている重心の座標を使って理想直線の交点との相違、すなわち位置ずれ量を算出する（Ｓ１９０２）。そして、カウンタｘをカウントアップし（Ｓ１９０３）、カウンタｘがｎに達したか否か判断する（Ｓ１９０４）。
【００７６】
ステップＳ１９０４の判断の結果、カウンタｘがｎに達していない場合（Ｓ１９０４：Ｎｏ）、次の行について位置ずれ量の算出を実行する。一方、カウンタｘがｎに達している場合（Ｓ１９０４：Ｙｅｓ）、カウンタｘを０にリセットすると共にカウンタｙをカウントアップする（Ｓ１９０５）。そして、カウンタｙの値がｍに達したか否か判断し（Ｓ１９０６）、カウンタｙがｍに達していない場合（Ｓ１９０６：Ｎｏ）、次の列について位置ずれ量の算出を実行する。また、カウンタｙがｍに達している場合（Ｓ１９０６：Ｙｅｓ）、処理を終了する。
【００７７】
なお、以上述べた本実施形態のドット位置検出方法をコンピュータに実行させるプログラムは、インストール可能な形式又は実行可能な形式のファイルでＣＤ−ＲＯＭ、フロッピー（登録商標）ディスク（ＦＤ）、ＤＶＤ等のコンピュータで読み取り可能な記録媒体に記録されて提供される。また、本実施形態のドット位置検出方法をコンピュータに実行させるプログラムを、インターネット等のネットワークに接続されたコンピュータ上に格納し、ネットワーク経由でダウンロードさせることにより提供するように構成しても良い。
【００７８】
また、本発明は、以上述べた構成に限定されるものではない。すなわち、以上述べた実施形態で想定される画像よりも大きな画像に本発明を適用した場合、ヒストグラムによるドット位置予測ができなくなる。このような場合、本発明のドット位置検出装置及びドット位置検出方法をコンピュータに実行させるプログラムでは、画像の一部を用いて理想着弾位置を計算し、その理想着弾位置を全画像への推測ドット位置として適用する。次いで、全画像に対して重心位置から位置ズレまでを計算してもよい。このような実施形態を図２０に模式的に示す。
【図面の簡単な説明】
【００７９】
【図１】本発明の一実施形態のドット位置検出装置を説明するための図である。
【図２】図１に示した解析アプリケーションを動作させるハードウェア構成を説明するための図である。
【図３】図２に示したドット解析部の機能ブロック図である。
【図４】図３に示したドット解析部でなされる処理の概略を示すフローチャートである。
【図５】本発明の一実施形態のドット中心予測の処理を説明するための図である。
【図６】本発明の一実施形態のドット中心予測の処理のフローチャートである。
【図７】本発明の一実施形態のドットの重心の位置を計算する処理を説明するための図である。
【図８】本発明の一実施形態のドットの重心の位置を計算する処理のフローチャートである。
【図９】図８に示したステップＳ８０３の処理を説明するためのフローチャートである。
【図１０】本発明の一実施形態の理想ドット位置の表し方を説明するための図である。
【図１１】本発明の一実施形態のドットを通る理想直線の傾きと切片との求め方を説明するための図である。
【図１２】本発明の一実施形態の初期切片Ｙ０を示す図である。
【図１３】本発明の一実施形態の初期切片Ｘ０を示す図である。
【図１４】本発明の一実施形態の理想位置計算部でなされる処理全般を説明するためのフローチャートである。
【図１５】図１４に示したステップＳ１４０１の処理を説明するためのフローチャートである。
【図１６】図１４に示したステップＳ１４０３の処理を説明するためのフローチャートである。
【図１７】図１６に示したステップＳ１６０２の処理を説明するためのフローチャートである。
【図１８】本発明の一実施形態の位置ずれ計算の手法を説明するための図である。
【図１９】本発明の一実施の位置ずれ計算部で行われる処理を説明するためのフローチャートである。
【図２０】本発明の一実施形態の応用例を示した図である。
【符号の説明】
【００８０】
２００  解析アプリケーション、２０１  解析画像入力部、２０３  しきい値設定部
２０５  ドット解析部
２０７  結果出力部、２０９  入力装置、２１１  出力装置、
３０５  ドット位置予測部、３０７  ドット重心位置計算部、３０９  理想位置計算部
３１１  位置ズレ計算部、３１３  出力部
【図１】
2006142546.tif 000003
【図２】
2006142546.tif 000004
【図３】
2006142546.tif 000005
【図４】
2006142546.tif 000006
【図５】
2006142546.tif 000007
【図６】
2006142546.tif 000008
【図７】
2006142546.tif 000009
【図８】
2006142546.tif 000010
【図９】
2006142546.tif 000011
【図１０】
2006142546.tif 000012
【図１１】
2006142546.tif 000013
【図１２】
2006142546.tif 000014
【図１３】
2006142546.tif 000015
【図１４】
2006142546.tif 000016
【図１５】
2006142546.tif 000017
【図１６】
2006142546.tif 000018
【図１７】
2006142546.tif 000019
【図１８】
2006142546.tif 000020
【図１９】
2006142546.tif 000021
【図２０】
2006142546.tif 000022