(57)【要約】      （修正有）【課題】誤差拡散法のドット発生の遅延や中脳度領域でのテクスチャ発生を抑制して高い画質の二値化データに二値化することができる画像処理装置及び画像処理プログラムを提供する。【解決手段】入力画像データを画素毎の二値化データ１９に変換する画像処理装置において，（１）多階調画像データ１３の最大階調値の１／２未満の基準値より処理対象画素の多階調画像データの階調値が低い場合は，複数の低い閾値を集中させることなく複数の閾値が分散配置された閾値マトリクスに基づいて当該処理対象画素をドット有りとドットなしの二値化データに変換する第１の二値化処理を行い，（２）高い場合は，当該処理対象画素をドット有りの二値化データに変換する第２の二値化処理と，処理対象画素の入力階調値と最大階調値との差からなる誤差を第１の誤差伝搬マトリクスに基づいて周囲の画素の多階調画像データに加算する誤差伝搬処理とを行う。【選択図】図１
2007006271.tif 000002
【特許請求の範囲】
【請求項１】
所定の解像度の画素毎に多階調画像データを有する入力画像データを前記画素毎の二値化データに変換する画像処理装置において，
前記多階調画像データの最大階調値の１／２未満の基準値より処理対象画素の多階調画像データの階調値が低い場合は，複数の低い閾値を集中させることなく複数の閾値が分散配置された閾値マトリクスに基づいて当該処理対象画素をドット有りとドットなしの二値化データに変換する第１の二値化処理を行い，
前記基準値より処理対象画素の多階調画像データの階調値が高い場合は，当該処理対象画素をドット有りの二値化データに変換する第２の二値化処理と，前記処理対象画素の入力階調値と最大階調値との差からなる誤差を第１の誤差伝搬マトリクスに基づいて周囲の画素の前記多階調画像データに加算する誤差伝搬処理とを行う二値化ユニットを有する画像処理装置。
【請求項２】
請求項１において，
前記二値化ユニットは，前記第１の二値化処理と共に，当該第１の二値化処理が行われた処理対象画素の誤差を周囲の画素の前記多階調画像データに加算する誤差伝搬処理を行うことを特徴とする画像処理装置。
【請求項３】
請求項１において，
前記二値化ユニットは，前記第１の二値化処理においてドット有りに二値化された処理対象画素の誤差を，前記第１の誤差伝搬マトリクスに基づいて周囲の画素の前記多階調画像データに加算する誤差伝搬処理を行うことを特徴とする画像処理装置。
【請求項４】
請求項３において，
前記第１の誤差伝搬マトリクスは，処理対象画素から第１の方向に位置する周囲の画素により高い重み付けで前記誤差の全てを伝搬することを特徴とする画像処理装置。
【請求項５】
請求項１において，
前記二値化ユニットは，前記第１の二値化処理においてドットなしに二値化された処理対象画素の誤差を，第２の誤差伝搬マトリクスに基づいて周囲の画素の前記多階調画像データに加算する誤差伝搬処理を行うことを特徴とする画像処理装置。
【請求項６】
請求項５において，
前記第２の誤差伝搬マトリクスは，処理対象画素から第２の方向に位置する周囲の画素により高い重み付けで前記誤差の一部を伝搬することを特徴とする画像処理装置。
【請求項７】
請求項１において，
更に，前記画素毎に多階調画像データを有する入力画像データを，前記二値化ユニットによる二値化処理に伴って発生する，入力画像データの階調値と二値化データにより再現される画像の階調値とのノンリニア特性を補正するガンマ補正処理ユニットを有し，
前記ガンマ補正処理ユニットにより補正された補正入力画像データが，前記二値化ユニットにより二値化されることを特徴とする画像処理装置。
【請求項８】
請求項７において，
前記ガンマ補正処理ユニットは，更に，前記入力画像データに対して低濃度階調領域でより低濃度階調に補正するドット除去補正と，高濃度階調領域でより高濃度階調に補正する黒強調補正を行うことを特徴とする画像処理装置。
【請求項９】
請求項８において，
前記ガンマ補正処理ユニットは，前記ノンリニア特性の補正と前記ドット除去補正と前記黒強調補正とが合成された合成補正テーブルに基づいて，前記入力画像データを補正することを特徴とする画像処理装置。
【請求項１０】
請求項１，２，３，５のいずれかにおいて，処理対象画素の多階調データが複数ビットからなるバイナリデータ形式で処理され，前記誤差伝搬処理において，当該多階調データを所定数ビット右シフト処理することで，伝搬すべき誤差を算出することを特徴とする画像処理装置。
【請求項１１】
所定の解像度の画素毎に多階調画像データを有する入力画像データを前記画素毎の二値化データに変換する画像処理プログラムにおいて，
前記多階調画像データの最大階調値の１／２未満の基準値より処理対象画素の多階調画像データの階調値が低い場合は，複数の低い閾値を集中させることなく複数の閾値が分散配置された閾値マトリクスに基づいて当該処理対象画素をドット有りとドットなしの二値化データに変換する第１の二値化処理工程をコンピュータに実行させ，
前記基準値より処理対象画素の多階調画像データの階調値が高い場合は，当該処理対象画素をドット有りの二値化データに変換する第２の二値化処理工程と，前記処理対象画素の入力階調値と最大階調値との差からなる誤差を第１の誤差伝搬マトリクスに基づいて周囲の画素の前記多階調画像データに加算する誤差伝搬処理工程とをコンピュータに実行させる画像処理プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は，多階調画像データを中間階調表示可能な二値化データに変換する二値化処理（ハーフトーニング処理）を行う画像処理装置及び画像処理プログラムに関し，特に，低解像度画像データを高画質の二値化データに変換することができる画像処理装置及び画像処理プログラムに関する。
【背景技術】
【０００２】
ファクシミリや複写機は，紙などの印刷媒体に形成された画像をスキャナユニットにより読み取り多階調画像データを取得し，それを中間階調表示可能な二値化データに変換する二値化処理を行う。ファクシミリは，その二値化データを通信回線を経由して送信する。また，複写機は，二値化データを印刷ユニットにより画像再生する。また，ファクシミリや複写機以外にも，パーソナルコンピュータシステムが，スキャナ機能を利用して印刷媒体の画像を読み取り，読み取った多階調画像データを二値化データに変換して，例えば表示装置に画像再生して表示する。
【０００３】
特に，ファクシミリの場合，標準化規格により200dpi程度の低解像度の二値化データを生成することが要求されており，そのスキャナユニットはそれに対応して200dpi程度の低解像度の多階調画像データを出力する。このようにスキャナユニットにより取得された低解像度の多階調画像データを二値化データに変換する場合，細線や文字の視認性を向上しつつ写真の階調性を維持することが必要になる。
【０００４】
一般に，多階調画像データの二値化処理として，固定された閾値マトリクスを使用する固定ディザ法や誤差拡散法などがある。さらに，固定ディザ法には，閾値マトリクスの低閾値を中心部分に配置することで入力階調が高くなるにしたがって網点または線が成長するドット集中型固定ディザ法と，閾値マトリクスの閾値を分散させたドット分散型固定ディザ法とがある。ドット集中型は，ドットが集中した網点の大きさで階調を表現することからＡＭ（振幅変調）スクリーンと呼ばれ，ドット分散型は，孤立ドットの間隔の粗密により画像を再現することからＦＭ（周期変調）スクリーンと呼ばれている。
【０００５】
ドット集中型固定ディザ法では，階調性を重視してスクリーン線数を増加させれば解像性が低下し，解像性を重視してスクリーン線数を低下させれば階調性が低下するので，低解像度画像データを二値化する場合，細線や文字の視認性を向上しつつ写真の階調性も維持することは困難である。
【０００６】
また，インクジェットプリンタなどで採用されているドット分散型固定ディザ法は，ドットを分散させてドットの粗密により階調を表現するので階調性は向上するものの，低解像度画像データの場合ドット抜けが目立って文字や細線の視認性を確保することができない。
【０００７】
一方，誤差拡散法は，二値化により発生した入力画像データとの誤差を周囲の画素に伝搬しながら二値化を繰り返す処理であり，階調性と解像性を両立させることができる二値化処理として知られている。かかる誤差拡散法については，例えば特許文献１に記載されている。
【特許文献１】特開平８−３０７６８０号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
しかしながら，誤差拡散法は，二値化により発生した誤差を周囲の画素に伝搬するので，低濃度領域においてドットの生成が遅れるドット生成の遅延や，中濃度領域において誤差拡散法特有のテクスチャの発生があり，画質の低下を招いている。誤差拡散法では固定された誤差伝搬マトリクスで誤差を周囲の画素に伝搬するので，中濃度領域ではドットの発生に規則性が生じ，ワームノイズと呼ばれるテクスチャが生成されやすい。したがって，誤差拡散法は，文字や細線の解像性は優れてはいるものの，低濃度領域でのドット生成の遅延や中濃度領域でのテクスチャ発生に問題がある。
【０００９】
そこで，本発明の目的は，誤差拡散法のドット発生の遅延や中濃度領域でのテクスチャ発生を抑制して高い画質の二値化データに二値化することができる画像処理装置及び画像処理プログラムを提供することにある。
【課題を解決するための手段】
【００１０】
上記の目的を達成するために，本発明の第１の側面によれば，所定の解像度の画素毎に多階調画像データを有する入力画像データを前記画素毎の二値化データに変換する画像処理装置において，
前記多階調画像データの最大階調値の１／２未満の基準値より処理対象画素の多階調画像データの階調値が低い場合は，複数の低い閾値を集中させることなく複数の閾値が分散配置された閾値マトリクスに基づいて当該処理対象画素をドット有りとドットなしの二値化データに変換する第１の二値化処理を行い，
前記基準値より処理対象画素の多階調画像データの階調値が高い場合は，当該処理対象画素をドット有りの二値化データに変換する第２の二値化処理と，前記処理対象画素の入力階調値と最大階調値との差からなる誤差を第１の誤差伝搬マトリクスに基づいて周囲の画素の前記多階調画像データに加算する誤差伝搬処理とを行う二値化ユニットを有する画像処理装置である。
【００１１】
上記の第１の側面によれば，基準値より低い階調値の場合は，ドット分散型固定ディザ法により二値化するので，誤差拡散法に特有の低階調濃度領域におけるドット再現の遅延をなくし中階調濃度領域におけるテクスチャの発生を抑制することができる。また，基準値より高い階調値の場合は，常にドット有りの二値化データに変換して誤差を周囲の画素に伝搬する誤差拡散法による二値化処理をするので，低解像度の入力画像データであっても，文字や細線のドット抜けによる視認性の低下を抑制することができる。
【００１２】
第１の側面において，好ましい態様によれば，前記二値化ユニットは，前記第１の二値化処理と共に，当該第１の二値化処理が行われた処理対象画素の誤差を周囲の画素の前記多階調画像データに加算する誤差伝搬処理を行う。かかる態様によれば，低濃度領域の画素についても二値化に伴う誤差を周囲の画素に伝搬するので，高濃度領域の画素での誤差伝搬の連続性を保つことができ，ドットの粗密な部分で生じるドットむらを抑制することができる。文字や細線の視認性や低濃度領域でのドット発生の遅延については，必ずしもこの誤差伝搬処理には依存しない。
【発明を実施するための最良の形態】
【００１３】
以下，図面にしたがって本発明の実施の形態について説明する。但し，本発明の技術的範囲はこれらの実施の形態に限定されず，特許請求の範囲に記載された事項とその均等物まで及ぶものである。
【００１４】
以下の実施の形態では，低解像度画像を二値化する例として，ファクシミリの画像処理を一例にして説明する。ただし，ファクシミリ以外でも低解像度画像を二値化する場合に本実施の形態を適用することができる。
【００１５】
図１は，本実施の形態における画像処理ユニットを有するファクシミリ装置の構成図である。図中，送信側１００の構成と受信側２００の構成とが示されている。送信側の構成１００は，オリジナル画像１０をスキャンして画素毎の多階調画像データ１３を生成するスキャナ１２と，スキャナ１２が生成した多階調画像データ１３を画像処理して画素毎の二値化データ（ドット有り，ドットなしのデータ）１９を生成する画像処理ユニット１４と，二値化データ１９を送出する送信ユニット２０とを有する。多階調画像データ１３は，例えば２００ｄｐｉの低解像度の画素毎に８ビット，２５６階調の階調表現を行うことができるデータである。ファクシミリ装置の場合，標準規格により送信データの解像度が２００ｄｐｉと定められている関係上，スキャナ１２はそれと同程度の解像度を有するのが一般的である。それより高い解像度で多階調画像データを生成しても二値化データの解像度がそれより低い場合は，スキャナを高解像度にする意味がないからである。
【００１６】
本実施の形態における画像処理ユニット１４は，スキャナ１２が生成した多階調画像データ１３を補正するガンマ補正ユニット１６と，補正された多階調画像データを二値化する二値化ユニット１８とで構成される。画像処理ユニット１４は，それぞれ専用のハードウエアによるガンマ補正ユニットと二値化ユニットで構成される場合もあれば，内蔵するＣＰＵユニットがガンマ補正プログラムと二値化プログラムを実行することにより，ガンマ補正処理と二値化処理を行う場合もある。
【００１７】
受信側の構成２００は，送信された二値化データを受信する通信ユニット２２と，受信した二値化データ２３にしたがって画素にドット有り，ドットなしで画像再生を行う印刷ユニット２４とを有する。印刷ユニット２４が二値化データに基づいて画像再生を行うことで，受信画像２６が形成される。
【００１８】
［ハイブリッド誤差拡散による二値化処理］
本実施の形態における二値化処理は，誤差拡散法とドット分散型固定ディザ法とを併用することで，誤差拡散法による文字や細線の高解像度性を維持したまま，低濃度領域におけるドット生成の遅延を抑制し，中濃度領域でのテクスチャ（ワームノイズ）の発生を抑制する。そして，二値化された画素の誤差を周囲の画素の階調値に加算する誤差拡散処理を行う。
【００１９】
図２は，本実施の形態における二値化ユニットの二値化処理を示すフローチャート図である。この二値化処理によれば，処理済みの画素から伝搬された誤差を加えた注目画素の階調値（濃度）が，最大階調値２５５の半分より小さい基準値，例えば６４，より小さい低濃度領域か，基準値より大きい高濃度領域かに応じて，固定ディザ法による二値化処理と，常にドット有りとして誤差を拡散する誤差拡散処理とを区別して行う。つまり，注目されている処理対象画素の濃度が基準値（６４）より小さい場合は（Ｓ１０のＹＥＳ），固定ディザ法による二値化処理を行う（Ｓ１２）。この固定ディザ法による二値化処理は，例えばドット集中型固定ディザ法よりもドット分散型固定ディザ法のほうが好ましい。そして，固定ディザ法による二値化処理した場合に，その二値化により生成される誤差を後述する誤差伝搬マトリクスによって周囲の画素の多階調画像データの階調値に加算（伝搬）する。
【００２０】
また，処理対象画素の濃度が基準値以上の場合は（Ｓ１０のＮＯ），基準値を閾値とする誤差拡散法により二値化処理する（Ｓ１６）。つまり，処理対象画素をドット有りの二値化データに変換し，それに伴う誤差を第１の誤差伝搬マトリクスに基づいて周囲の画素の多階調画像データの階調値に加算（伝搬）する。即ち，この誤差拡散法により二値化処理Ｓ１６は，基準値を閾値とする誤差拡散法と同等である。
【００２１】
図３は，ドット分散型固定ディザ法における閾値マトリクスの例を示す図である。この例では，多階調画像データの階調値が０〜１５を前提にしており，閾値マトリクスには４ｘ４の画素マトリクスにそれぞれ閾値（０〜１５）が分散して配置されている。多階調画像データの階調値が０〜２５５の場合は，閾値マトリクスのサイズに応じた離散的値が分散して配置される。例えば図３の４ｘ４のマトリクス例では０〜１５の閾値が、０〜２５５に正規化した１７刻みの離散的閾値となる。このように，注目画素の濃度が低濃度領域の場合にドット分散型固定ディザ法による二値化処理を行うことで，低解像度の画像データであっても文字や細線のドット欠けを減らして視認性を高めることができる。
【００２２】
そして，固定ディザ法による二値化処理をした場合は，注目されている処理対象画素がドット有りの二値化データに変換されると，それによる誤差を第１の誤差伝搬マトリクスで周囲の画素の多階調画像データの階調値に加算（伝搬）する。一方，処理対象画素がドットなしの二値化データに変換されると，それによる誤差を第１の誤差伝搬マトリクスとは異なる第２の誤差伝搬マトリクスで伝搬する（Ｓ１４）。つまり，処理対象画素の多階調画像データの階調値が図３に示した閾値マトリクス（但し閾値は０〜２５５）の閾値以上の場合はドット有り（階調値＝２５５）の二値化データに変換されると，処理対象画素の入力階調値とドット有りの階調値（最大階調値２５５）との差であるマイナス誤差は，高濃度画素の誤差伝搬に使用する誤差伝搬マトリクスと同じ第１の誤差伝搬マトリクスにしたがって周囲の画素に伝搬する。一方，処理対象画素の多階調画像データの階調値が閾値マトリクスの閾値より小さい場合はドットなし（階調値＝０）の二値化データに変換されると，処理対象画素の入力階調値とドットなしの階調値（最小階調値０）との差であるプラス誤差は，第１の誤差伝搬マトリクスとは誤差伝搬方向が異なり，伝搬する誤差の量も異なる第２の誤差伝搬マトリクスにしたがって周囲の画素に伝搬する。
【００２３】
図４は，第１及び第２の誤差伝搬マトリクスの一例を示す図である。第１の閾値伝搬マトリクスＥＴＭ１は，処理対象画素である注目画素（図中※）で二値化により発生した誤差（注目画素の入力階調値と二値化後の階調値（０または２５５）の差分）を周囲の画素に伝搬する量の割合を示している。この例は，Javis-Judice-Ninke型の誤差伝搬マトリクスであり，誤差を３行５列の周囲の画素に伝搬する割合を有し，その割合は，右横と下の画素に７／４８，右横２番目と右下と左下の画素に５／４８，右下２番目と左下２番目とに３／４８，そして，右下３番目と左下３番目とに１／４８が配置される。したがって，第１の誤差伝搬マトリクスＥＴＭ１は，誤差の全てを周囲の画素に伝搬し，伝搬方向が注目画素を中心に且つやや右下方向に広がっているところに特徴がある。
【００２４】
一方，第２の誤差伝搬マトリクスＥＴＭ２は，注目画素で二値化により発生した誤差の一部を周囲の画素に伝搬している。つまり，誤差は第１の誤差伝搬マトリクスＥＴＭ１より小さい２行３列のマトリクスサイズであり，その割合は，注目画素の左下に５／１６，下に３／１６，左下２番目に１／１６が配置される。したがって，第２の誤差伝搬マトリクスＥＴＭ２は，誤差の一部９／１６を，周囲の３つの画素に図示された位置に伝搬する。この伝搬方向は，第１の誤差伝搬マトリクスのやや右下とは異なるやや左下方向となっている。
【００２５】
本実施の形態における二値化処理の意味するところは以下の通りである。まず，注目されている処理対象画素の濃度が中低濃度の場合（階調値が６４未満），ドット分散型固定ディザ法により，閾値マトリクスの閾値より高い階調値ならドット有りに，低い階調値ならドットなしに二値化する。このようにドット分散型固定ディザ法による二値化処理を行うことで，文字や細線のドット欠けを抑制し視認性を向上させることができる。また，誤差拡散法に特有の低濃度領域でのドット発生の後れをなくすことができる。
【００２６】
そして，ドット分散型固定ディザ法で二値化した画素について，固定ディザ法では誤差の拡散は通常行われないにもかかわらず，本実施の形態では誤差を周囲の画素に伝搬する。しかも，誤差の伝搬は，ドット有りに二値化された画素のマイナス誤差は，全誤差を所定の割合で配置して且つやや右下方向に伝搬する第１の誤差伝搬マトリクスＥＴＭ１に基づいて周囲の画素に伝搬する。一方，ドットなしに二値化された画素のプラス誤差は，一部の誤差を所定の割合で配置して且つやや左下方向に伝搬する第２の誤差伝搬マトリクスＥＴＭ２に基づいて周囲の画素に伝搬する。つまり，第１に，マイナス誤差とプラス誤差の誤差拡散の方向を異ならせることで，誤差拡散法に特有の固定方向に誤差伝搬することによるテクスチャ（ワームノイズ）を抑制することができるものと考えられる。第２に，最大階調値の１／２より低い基準値より更に低い濃度の画素に対して閾値マトリクスにより二値化しているので，マイナス誤差は全誤差を伝搬して画像がドット過剰になるのを防止するとともに，プラス誤差は一部の誤差だけを伝搬して低濃度領域でドット不足になるのを防止する。
【００２７】
次に，注目されている処理対象画素の濃度が中高濃度の場合（階調値が６４以上），画素を全てドット有りの二値化データにし，そのマイナス誤差を第１の誤差伝搬マトリクスＥＴＭ１に基づいて周囲の画素に伝搬する。つまり，基準値６４を閾値とする誤差拡散法による二値化処理である。これにより，最大階調値の１／２を閾値とする誤差拡散法に特有の中濃度領域でのテクスチャ（ワームノイズ）の発生を抑制することができる。つまり，最大階調値の１／２より低い濃度であっても６４以上であればドット有りに二値化されるので，テクスチャの発生は抑制される。ただし，全誤差を第１の誤差伝搬マトリクスで伝搬させることで，中濃度領域での適切な中間階調を二値化データで表現することができる。また，基準値を６４と低めに設定することで，文字や細線に対してドットの発生が起きやすくなり視認性を向上できる。
【００２８】
さらに，画素の濃度が低濃度領域，つまり階調値が６４以下で且つ固定ディザ法でドット有りに二値化された画素については，中高濃度領域の誤差拡散法で利用した第１の誤差伝搬マトリクスを使って誤差を周囲の画素に伝搬しているので，中濃度から高濃度領域での誤差伝搬の連続性を保って，異なる二値化処理による不自然さを抑制している。
【００２９】
また，画素の濃度が中低濃度領域では，誤差拡散法を採用せずに固定ディザ法を採用するので，誤差拡散法に特有の低濃度領域でのドット発生の遅延を抑制することができ，一律に誤差拡散法で処理するよりも中低濃度領域の中間階調の画質を向上させることができる。つまり，低濃度であっても閾値マトリクスの閾値よりも大きければドットが確実に発生するので，ドット発生の遅延を抑制できるのである。
【００３０】
［ガンマ補正］
図１に示したとおり，本実施の形態の画像処理ユニット１４では，二値化ユニット１８による二値化処理に先立って，ガンマ補正ユニット１６が入力される多階調画像データ１３をガンマ補正している。このガンマ補正により，本実施の形態において採用したハイブリッド誤差拡散による二値化処理で生じるノンリニア特性を補正する。つまり，ハイブリッド誤差拡散法では，基準値より高濃度は常にドット有りに二値化しその誤差を伝搬し，基準値より低濃度は閾値マトリクスにより二値化して誤差を伝搬させている。このハイブリッド誤差拡散法により二値化前の多階調画像データの入力階調値と二値化データの出力階調値との間にノンリニアな特性が含まれる。そこで，このノンリニア特性を是正するために，ガンマ補正が行われる。
【００３１】
図５は，本実施の形態におけるガンマ補正処理を示す図である。第１象限に示した濃度特性曲線３０が，ハイブリッド誤差拡散による二値化処理の入出力特性である。二値化処理の入力階調値に対して二値データの出力階調値とは，例えば特性曲線３０に示されるような下に凸のノンリニア特性を有する。つまり，理想的なリニア特性３２に比較すると，入力階調ＩＮに対する二値化後の出力階調ＯＵＴは，理想的なリニア特性３２による二値化された場合の出力階調ＣＯＵＴよりも小さくなっている。
【００３２】
そこで，本実施の形態では，第４象限に示される補正ガンマテーブル３４によって入力階調ＩＮを補正し，その補正入力ＣＩＮを二値化ユニットにより二値化処理する。この補正ガンマテーブル３４は，特性曲線３０とは逆に上に凸の特性を有する一次元テーブルであり，入力ＩＮがこの補正ガンマテーブル３４で補正されると補正入力ＣＩＮは入力ＩＮよりも大きい値になる。
【００３３】
第３象限は補正入力ＣＩＮを第２象限の横軸に移動するための折り返し線３６を有する。つまり，第３，４象限の縦軸の補正入力ＣＩＮと第２象限の横軸の補正入力ＣＩＮとは同じ値である。
【００３４】
最後に，第２象限には特性曲線３０と左右対称の特性曲線３０Ｒが示され，補正入力ＣＩＮをこの特性曲線３０Ｒで二値化処理された補正出力ＣＯＵＴは，入力ＩＮを理想的なリニア特性３２で二値化した場合の出力値に一致する。
【００３５】
上記の補正ガンマテーブル３４は，次のようにして求められる。まず，入力画像データとして所定の階調値のパッチ画像データを生成する。このパッチ画像データは，白（階調値０）〜黒ベタ（階調値２５５）までを所定の階調刻みの階調値を有する複数のパッチパターンデータからなる。そして，このパッチ画像データを二値化ユニット１８によりハイブリッド誤差拡散法で二値化データに変換し，その二値化データで画像を形成し，形成された画像の各階調のパッチパターン濃度を測定することで，二値化ユニットの特性曲線３０を検出する。この二値化データの画像形成には，ファクシミリ装置の印刷ユニットを利用して紙媒体に画像を印刷することで行われる。または，二値化データによる画像を表示装置に表示しそのパッチパターン濃度を測定することでも良い。印刷ユニットを利用した場合は，特性曲線３０に印刷ユニットの特性も含まれてしまう。また，表示装置に表示した場合も，同様に特性曲線３０にその表示装置の特性が含まれる。このようにして検出した特性曲線３０に対して，理想的なリニア特性が得られるように，補正ガンマテーブル３４が求められる。これが濃度補正用ガンマテーブルである。
【００３６】
図６は，本実施の形態における濃度補正用ガンマテーブルの具体例を示す図である。この具体例では，入力画像データとして，階調値３３刻みで白（０）から黒（２５５）までの８ビットのパッチ画像データを生成し，ハイブリッド誤差拡散法で二値化データを生成し，その二値化データを印刷して，各階調のパッチ画像の出力濃度を測定した。この測定値は６４ビットの浮動小数点による実数精度の出力濃度データである。この実数精度の出力濃度データから濃度補正用ガンマテーブル３４が求められる。
【００３７】
この濃度補正用ガンマテーブル３４に基づいて，ガンマ補正ユニット１６が入力される多階調画像データの階調値を補正入力階調値に変換する。この時，補正入力階調値は浮動小数点のビット形式による実数精度である。したがって，より高精度の補正が可能である。補正された入力階調値が二値化ユニット１８により二値化されるので，二値化ユニット１８は，実数精度の補正入力階調値に対応して，各画素の多階調画像データを６４ビットで処理する。そのため，二値化処理には多くの処理負荷が発生する。そこで，より好ましい実施の形態の形態では，後述するとおり，６４ビットの実数精度を１６ビットの固定小数点のビット形式の整数精度にし，拡散する誤差の計算を除算ではなく右ビットシフトを利用する。
【００３８】
図７は，本実施の形態における背景除去と黒強調の画像処理用ガンマテーブルを示す図である。モノクロファクシミリの場合，オリジナル画像の印刷物の紙白部分に濃度情報が発生する。この濃度情報が，二値化データでドットを生成させ，画質の低下と二値化データの圧縮率（例えばランレングス圧縮の場合白（階調値０）の繰り返し頻度）が低下する。そこで，紙白領域にドットを発生させないための背景除去処理を行う。さらに，黒文字の品質を向上するために黒強調処理を行う。これらの処理は，いずれもスキャナ１２の特性も考慮して行う必要がある。
【００３９】
図７に示される変換用のガンマテーブル４０は，入力が低濃度領域４０Ａでは補正出力がより低く，具体的には階調値０になるようにし，入力が高濃度領域４０Ｂでは補正出力がより高く，具体的には階調値２５５になるようにする。このようなガンマテーブル４０を利用して入力の多階調画像データ１３を補正出力に変換することで，低濃度領域での背景除去処理と高濃度領域での黒強調処理を行うことができる。
【００４０】
上記の背景除去と黒強調の画像処理は，ガンマ補正ユニット１６により行われる。この画像処理用ガンマテーブル４０は固定テーブルとし，濃度補正用の補正ガンマテーブル３４は再調整可能なような変動テーブルとすることが望ましい。補正ガンマテーブル３４にはハイブリッド誤差拡散によるノンリニア特性と印刷ユニットのノンリニア特性とが含まれ，印刷ユニットのノンリニア特性は経年変化により変動するからである。
【００４１】
また，画像処理用ガンマテーブル４０と補正ガンマテーブル３６とを合成して合成ガンマテーブルを生成しておき，その合成ガンマテーブルを参照して入力の多階調画像データの階調値を補正するようにしても良い。合成ガンマテーブルを利用することでガンマ補正ユニットでの処理工数を減らすことができる。
【００４２】
図８は，合成ガンマテーブルの具体例を示す図である。図中，濃度補正用の補正ガンマテーブル３６と画像処理用ガンマテーブル４０の例が示され，これらを合成した合成ガンマテーブル４４が示されている。合成するためには，画像処理用ガンマテーブル４０も出力を６４ビットの実数精度にしておく必要がある。
【００４３】
図９は，本実施の形態における画像処理ユニットの構成図である。この画像処理ユニットは，入力多階調画像データ１３を受信するインターフェースＩＦと，中央処理ユニットＣＰＵと，合成ガンマテーブルが格納されているメモリ４４と，二値化プログラムが格納されているメモリ５０と，誤差伝搬マトリクスが格納されているメモリＥＴＭと，二値化処理に利用される画像データバッファ５４及び二値化データが格納される二値化データバッファ５６を有するワークメモリ５２と，それらを接続する内部バスＢＵＳとを有する。
【００４４】
画像データバッファ５４は，図４に示した誤差伝搬マトリクスＥＴＭ１にしたがって二値化処理に伴う誤差を３行５列の周囲の画素に伝搬するので，画像データの水平方向の画素数に両側に２画素のダミー画素を加えた３行の容量を有する。合成ガンマテーブル４４により８ビットの入力多階調画像データ１３が６４ビットの補正された多階調画像データに変換されるので，各画素は６４ビットのデータを格納する容量が必要になる。
【００４５】
画像処理ユニット１４は，合成ガンマテーブル４４を参照して入力多階調画像データ１３を補正された多階調画像データに変換する。この多階調画像データは，前述したとおり浮動小数点の実数精度であり１画素６４ビットのデータ長を有する。補正された多階調画像データは，画像データバッファ５４に一時的に格納される。次に，画像処理ユニット１４の中央処理ユニットＣＰＵは，二値化プログラム５０を実行して，補正された多階調画像データをハイブリッド誤差拡散による二値化処理を行って，画素毎にドット有り，なしの二値化データに変換し，二値化データバッファ５６に格納する。このハイブリッド誤差拡散による二値化処理中，画像データバッファ５４内の各画素の多階調画像データには，二値化処理で発生した誤差が誤差伝搬マトリクスＥＴＭに基づいて加算され，その誤差が加算された画像データに対して二値化処理が行われる。
【００４６】
図１０，図１１は，本実施の形態における二値化プログラムのフローチャート図である。この二値化プログラムは，ガンマ補正処理と二値化処理の両方を含んでいるが，ガンマ補正された画像データに対して二値化処理するようにしても良い。
【００４７】
二値化プログラムは，まず入力画像幅Ｗｄと高さＨｇを取得し（Ｓ２０），幅Ｗｄ＋４，高さ３行の画像データバッファ５４を確保し，初期化する（Ｓ２２）。そして，変数の処理画素ｗ，ｈをｗ＝１，ｈ＝１に初期化する。この変数ｗ，ｈは画像データバッファ５４内での処理対象の画素の位置を示すものである。次に，１画素分の入力画像データ１３を入力し，変数ｗをｗ＝ｗ＋１に更新する（Ｓ２４）。入力画像データ１３は，合成ガンマテーブルを参照して補正された画像データに変換される（Ｓ２６）。これにより入力画像データＸが補正された画像データＸ‘に変換される。
【００４８】
そして，変換された画像データＸ‘は，画像データバッファ内のｈライン上に格納される（Ｓ３０）。上記の補正と画像データバッファ内への保存は，変数ｈが３に達するまで繰り返される（Ｓ２８，Ｓ３２，Ｓ３４）。変数ｈ＝３になると，変換画像データＸ’が画像データバッファの３ライン上に格納される（Ｓ３６，Ｓ３７，Ｓ３８）。変数ｈ＝３になると，それ以降は常に変換画像データＸ‘は画像データバッファの３ライン上に格納される（Ｓ３６）。
【００４９】
画像データバッファ内に３行分の補正画像データが格納されると，バッファ内の１ライン目の画素の補正画像データについてハイブリッド誤差拡散による二値化処理が行われる（Ｓ４０）。この二値化処理は，図１１に詳述されている。１ライン目の画素の二値化処理が完了すると，画像データバッファの更新が行われる（Ｓ４２）。具体的には，画像データバッファの更新処理では，２ライン，３ラインの画像データを１ライン，２ラインに移動し，３ライン目を初期化する。
【００５０】
最後に，変数ｈ＝Ｈｇに達すると（Ｓ４４），画像データバッファ内の３ライン分の画素の画像データがハイブリッド誤差拡散による二値化処理され（Ｓ４６），二値化プログラムを終了する。
【００５１】
上記のように，図１０の処理では，３ライン分の補正画像データを画像データバッファに格納し，二値化処理を行う。また，二値化処理中に二値化に伴う誤差が周囲の画素の階調値に加算され画像データバッファ内に再度格納される。ただし，誤差バッファを別に設けることで，画像データバッファ内には補正画像データを格納し，二値化処理で発生する誤差を誤差バッファに格納し，二値化処理は補正画像データと誤差とを加算した階調値に対して行うようにしても良い。
【００５２】
図１１のハイブリッド誤差拡散による二値化処理を説明する。まず，処理画素の位置を示す変数ｗをｗ＝０に初期化する（Ｓ５０）。そして，注目画素の多階調画像データＸ‘を画像データバッファから読み出して，処理画素の階調値を示す変数ＴａｇをＴａｇ＝Ｘ’とする（Ｓ５２）。この階調値Ｔａｇには，補正された画像データに二値化処理で伝搬された誤差が加算されている。
【００５３】
そして，注目画素の階調値Ｔａｇが，最大階調値２５５の１／２より小さい６４未満か否かの判定が行われる（Ｓ５４）。ここで，注目画素の階調値Ｔａｇが６４以上の場合は（Ｓ５４のＮＯ），その画素の二値化データは常にドット有り（黒）に分類される（Ｓ５６）。つまり，二値化後の階調値は最大階調値のｘ＝２５５となる。そして，その二値化により発生したマイナス誤差Ｅ＝Ｔａｇ−ｘが算出され（Ｓ５８），第１の誤差伝搬マトリクスＥＴＭ１に基づいて周辺画素に伝搬され，誤差Ｅが周辺画素の階調値に加算される（Ｓ６０）。第１の誤差伝搬マトリクスＥＴＭ１の例は，図４で示したとおりであり，全誤差をマトリクスに示した割合で分配する。この誤差伝搬マトリクスにより，誤差はやや右下方向に伝搬される。
【００５４】
このように階調値Ｔａｇが基準値６４以上の場合に，その画素をドット有りに二値化することで，閾値を最大階調値の１／２である１２７にした誤差拡散の場合に発生する中濃度領域でのテクスチャ（ワームノイズ）の発生を抑制することができる。
【００５５】
一方，注目画素の階調値Ｔａｇが６４未満の場合は（Ｓ５のＹＥＳ），図３に示したドット分散型ディザマトリクスから注目画素に対応する閾値Ｄｔｈを取得し（Ｓ６２），注目画素の階調値Ｔａｇが閾値Ｄｔｈ以下か否かが判定される（Ｓ６４）。注目画素の階調値Ｔａｇが閾値Ｄｔｈを越えている場合は（Ｓ６４のＮＯ），その画素の二値化データはドット有り（黒）に分類され（Ｓ５４），そのマイナス誤差Ｅが算出され（Ｓ５６），誤差Ｅが第１の誤差伝搬マトリクスＥＴＭ１にしたがって周囲の画素に伝搬される（Ｓ６０）。一方，注目画素の階調値Ｔａｇが閾値Ｄｔｈ以下の場合は（Ｓ６４のＹＥＳ），その画素の二値化データはドットなし（白）に分類される（Ｓ６６）。つまり二値化後の階調値は最小階調値のｘ＝０になる。そして，その二値化処理により発生したプラスの誤差Ｅ＝Ｔａｇ−ｘが算出され（Ｓ６８），そのプラス誤差Ｅが，第２の誤差伝搬マトリクスにより周辺の画素に伝搬される（Ｓ７０）。第２の誤差伝搬マトリクスＥＴＭ２の例は，図４で示したとおりであり，誤差の一部をマトリクスに示した割合で分配する。この誤差伝搬マトリクスにより，誤差はやや左下方向に伝搬される。
【００５６】
上記のように，階調値Ｔａｇが基準値未満になる中低濃度領域では，ドット分散型固定ディザ法により二値化処理が行われる。固定ディザ法の場合ではあるが，二値化により発生したマイナス誤差は，階調値Ｔａｇが基準値以上の場合との連続性を考慮して，第１の誤差伝搬マトリクスによって周囲の画素に伝搬される。これにより，基準値を境にして異なる二値化処理をしてもそれによる不連続な階調表現にはされない。また，マイナス誤差については，プラス誤差の第１の誤差伝搬マトリクスとは誤差伝搬方向と伝搬量が異なる第２の誤差伝搬マトリクスにより周囲の画素に伝搬する。このように伝搬することで，本発明者は経験的に中低濃度領域でのテクスチャが抑制されたことを確認した。これは，プラスとマイナスの誤差伝搬が互いに打ち消し合うことによる効果と考えられる。
【００５７】
上記のハイブリッド誤差拡散による二値化処理が，処理画素の位置を示す変数ｗがｗ＝Ｗｄに達するまで繰り返される（Ｓ７２，Ｓ７４）。以上で，１ライン分の画素のハイブリッド誤差拡散による二値化処理が完了する。
【００５８】
［階調値を固定小数点のビット構成にした高速化処理］
上記の実施の形態では，合成ガンマテーブルによる補正精度を上げるために補正後の階調値を浮動小数点によるビット形式（倍精度の６４ビット）にした。しかし，浮動小数点のビット形式の実数を誤差伝搬などのために４８で除算する処理は，コンピュータには多くの処理工数を要する。そこで，高速処理を可能にするために，階調値を整数の固定小数点によるビット形式（１６ビット，−３２７６７〜＋３２７６７）にし，誤差伝搬マトリクスの割合を２のべき乗で除算した割合にする。
【００５９】
すなわち，まず，ドット分散型閾値マトリクス内の閾値を８ビットから１６ビット（０〜３２７６７）に，実数精度であった合成ガンマテーブルも１６ビット（０〜３２７６７）に正規化しておく。これらは正の整数であるため、１６ビット中の負の整数領域（−３２７６７〜−１）は使用されない。さらに，画像データバッファの各画素を１６ビットの容量にする。また，低濃度領域判定のための基準値６４を，０〜３２７６７内の基準値に正規化した８２２４（６４×３２７６７／２５５）にするし、伝播誤差が１６ビット（最大で２16／２−１＝−３２７６７〜＋３２７６７）で伝搬されるようにする。
【００６０】
そして，誤差を第１の誤差伝搬マトリクスで伝搬するとき，誤差を２のべき乗分のｎ（ｎは整数）の割合で配置した誤差伝搬マトリクスを利用する。図１２は，２の５乗（６４）分のｎの割合で分配する第１の誤差伝搬マトリクスの例を示す図である。この例は，図４の第１の誤差伝搬マトリクスＥＴＭ１と同様に，全誤差をほぼ右下方向に伝搬するものである。これによれば，伝搬すべき誤差は，二値化により発生した誤差（Ｔａｇ−ｘ）を６４で除算し整数倍（９，７，４倍）した値になる。したがって，誤差伝搬処理において発生した誤差のバイナリデータを６ビット右シフトする処理により６４で除算し，それにそれぞれの整数倍することで，伝搬すべき誤差を簡単に求めることができ，プロセッサによる処理工数を削減することができる。
【００６１】
また，プラス誤差については，図４の第２の誤差伝搬マトリクスＷＥＴＭ２をそのまま利用することができる。この例では，発生した誤差のバイナリデータを４ビット右シフトして１６で除算することができる。このように，固定小数点によるビット形式にし，分配誤差の計算に右シフトを利用することで，演算処理工数を大きく削減することができる。
【００６２】
[二値化処理結果の比較]
本実施の形態のハイブリッド誤差拡散による二値化処理した結果と，従来の固定ディザ法，誤差拡散法による二値化結果とを比較する。第１の比較例は，中高濃度領域はドット有りに二値化して誤差を拡散し，中低濃度領域はディザ法による閾値マトリクスによる二値化を行う処理をした場合である。この第１の比較例では，中低濃度領域でディザ法で二値化した時に発生する誤差は伝搬させていない。第２の比較例は，同様に，中高濃度領域はドット有りに二値化して誤差を伝搬し，中低濃度領域はディザ法による閾値マトリクスによる二値化してその誤差を伝搬する処理をした場合である。
【００６３】
［第１の比較例］
図１３，１４は，中濃度領域のパッチ画像についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す。オリジナル画像（Ａ）には，階調値６８，６９，７０のパッチパターンが上段に，階調値８４，８５，８６のパッチパターンが下段に配置されている。これらは，低い基準値による二値化に伴うマイナス誤差の拡散で，二値化処理時には多くの画素の階調値が基準値未満になる。したがって，低中濃度領域での分散ディザ法とその誤差拡散の処理結果と，中濃度領域での分散ディザ法と誤差拡散法の境界部分での連続性とを比較することができる。ドット分散型固定ディザ法（Ｂ）では，やや大きなテクスチャ（ワームノイズ）が見られる。また，誤差拡散法（Ｃ）では，処理開始時のドット発生後れと大きなテクスチャ（ワームノイズ）が見られる。それに対して，本実施の形態（Ｄ）では，誤差拡散特有のワームノイズは抑制されているものの，誤差伝搬の連続性がないのでドットの粗密な部分で生じるムラが目立っている。
【００６４】
図１５は，文字についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す。ドット分散型固定ディザ法（Ｂ）では，予想とおり文字にドット欠けが目立ち視認性が低下している。誤差拡散法（Ｃ）では，ドット欠けが少なく視認性は高い。そして，本実施の形態（Ｄ）でも誤差拡散法と同様にドット欠けが少なく視認性が改善されている。
【００６５】
図１６は，低濃度領域を有するパッチパターンについて，オリジナル画像（Ａ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す。誤差拡散法（Ｃ）では，予想とおり低濃度領域でドット発生の遅延（領域６０）が生じている。それに対して，本実施の形態（Ｄ）では，低濃度領域をドット分散型固定ディザ法で二値化しているので，ドット発生の遅延は生じていない。
【００６６】
［第２の比較例］
図１７，１８は，中濃度領域のパッチ画像についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す。オリジナル画像（Ａ）には，階調値６８，６９，７０のパッチパターンが上段に，階調値８４，８５，８６のパッチパターンが下段に配置されている。これらは，低い基準値による二値化に伴うマイナス誤差の拡散で，二値化処理時には多くの画素の階調値が基準値未満になる。したがって，低中濃度領域での分散ディザ法とその誤差拡散の処理結果と，中濃度領域での分散ディザ法と誤差拡散法の境界部分での連続性とを比較することができる。ドット分散型固定ディザ法（Ｂ）では，やや大きなテクスチャ（ワームノイズ）が見られる。また，誤差拡散法（Ｃ）では，処理開始時のドット発生後れと大きなテクスチャ（ワームノイズ）が見られる。それに対して，本実施の形態（Ｄ）では，誤差伝搬の連続性があるのでテクスチャが抑制され画質が向上しているのがわかる。
【００６７】
図１９は，文字についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す。ドット分散型固定ディザ法（Ｂ）では，予想とおり文字にドット欠けが目立ち視認性が低下している。誤差拡散法（Ｃ）では，ドット欠けが少なく視認性は高い。そして，本実施の形態（Ｄ）でも誤差拡散法と同様にドット欠けが少なく視認性が改善されている。
【００６８】
図２０は，低濃度領域を有するパッチパターンについて，オリジナル画像（Ａ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す。誤差拡散法（Ｃ）では，予想とおり低濃度領域でドット発生の遅延（領域６０）が生じている。それに対して，本実施の形態（Ｄ）では，低濃度領域をドット分散型固定ディザ法で二値化しているので，ドット発生の遅延は生じていない。
【００６９】
以上説明したとおり，本実施の形態によれば，中低濃度領域ではドット分散型固定ディザ法で二値化し，中高濃度領域ではドット有りに二値化しそこで発生したマイナス誤差を第１の誤差伝搬マトリクスで伝搬しているので，誤差拡散法で低濃度領域で発生するドット発生の遅延や中濃度領域で発生するテクスチャ（ワームノイズ）の発生を抑えることができる。
【図面の簡単な説明】
【００７０】
【図１】本実施の形態における画像処理ユニットを有するファクシミリ装置の構成図である。
【図２】本実施の形態における二値化ユニットの二値化処理を示すフローチャート図である。
【図３】ドット分散型固定ディザ法における閾値マトリクスの例を示す図である。
【図４】第１及び第２の誤差伝搬マトリクスの一例を示す図である。
【図５】本実施の形態におけるガンマ補正処理を示す図である。
【図６】本実施の形態における濃度補正用ガンマテーブルの具体例を示す図である。
【図７】本実施の形態における背景除去と黒強調の画像処理用ガンマテーブルを示す図である。
【図８】合成ガンマテーブルの具体例を示す図である。
【図９】本実施の形態における画像処理ユニットの構成図である。
【図１０】本実施の形態における二値化プログラムのフローチャート図である。
【図１１】本実施の形態における二値化プログラムのフローチャート図である。
【図１２】高速化処理用の第１の誤差伝搬マトリクスの例を示す図である。
【図１３】中濃度領域のパッチ画像についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す図である。
【図１４】中濃度領域のパッチ画像についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す図である。
【図１５】文字についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す図である。
【図１６】低濃度領域を有するパッチパターンについて，オリジナル画像（Ａ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す図である。
【図１７】中濃度領域のパッチ画像についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す図である。
【図１８】中濃度領域のパッチ画像についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す図である。
【図１９】文字についてオリジナル画像（Ａ），ドット分散型固定ディザ（Ｂ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す図である。
【図２０】低濃度領域を有するパッチパターンについて，オリジナル画像（Ａ），誤差拡散（Ｃ），ハイブリッド誤差拡散（Ｄ）の二値化結果を示す図である。
【符号の説明】
【００７１】
１２：スキャナ            １３：入力多階調画像データ
１４：画像処理ユニット（画像処理装置）
１６：ガンマ補正ユニット  １８：二値化ユニット
１９：二値化データ
【図１】
2007006271.tif 000003
【図２】
2007006271.tif 000004
【図３】
2007006271.tif 000005
【図４】
2007006271.tif 000006
【図５】
2007006271.tif 000007
【図６】
2007006271.tif 000008
【図７】
2007006271.tif 000009
【図８】
2007006271.tif 000010
【図９】
2007006271.tif 000011
【図１０】
2007006271.tif 000012
【図１１】
2007006271.tif 000013
【図１２】
2007006271.tif 000014
【図１３】
2007006271.tif 000015
【図１４】
2007006271.tif 000016
【図１５】
2007006271.tif 000017
【図１６】
2007006271.tif 000018
【図１７】
2007006271.tif 000019
【図１８】
2007006271.tif 000020
【図１９】
2007006271.tif 000021
【図２０】
2007006271.tif 000022