(57)【特許請求の範囲】
【請求項１】
現フィールド画像及び前フィールド画像を、偽輪郭の発生する信号レベル近傍の区間をしきい値により２値化することによってそれぞれ現フィールド２値化画像及び前フィールド２値化画像とし、前記現フィールド２値化画像を複数の検出ブロックに分割し、前記検出ブロック毎に前記前フィールド２値化画像内に参照領域を設定し、前記参照領域内に参照ブロックを複数設定し、前記検出ブロックと前記参照ブロックとの一致度を評価し、評価値の最も高い参照ブロックと前記検出ブロックとの位置関係から動き画素の移動画素数及び移動方向を検出し、前記現フィールド画像で偽輪郭の発生する信号レベルを持った動き画素をその移動画素数および移動方向に応じて補正する動画像表示方法であって、
前記検出ブロックと前記参照ブロックとの一致度の評価において動きが検出されなかった前記検出ブロックを不定ブロックとして仮登録し、動きが検出された前記検出ブロックを既知ブロックとし、前記既知ブロックに囲まれた前記不定ブロックを前記既知ブロックと同じ移動画素数及び移動方向にすることを特徴とする動画像表示方法。
【請求項２】
前記不定ブロックの移動画素数及び移動方向は、前記不定ブロックを囲む前記既知ブロックの移動画素数及び移動方向により線形補間により求めることを特徴とする請求項１記載の動画像表示方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、プラズマディスプレイパネル（以下、単に「ＰＤＰ」と記述する）等において画像が移動する場合に発生する偽輪郭の発生を抑止する動画像表示方法及びその装置に関する。
【０００２】
【従来の技術】
近年の表示装置の大型化要望に応えるものとして、ＰＤＰ，ＥＬ表示素子、蛍光表示管、液晶表示素子等の薄型のマトリックスパネルが提供され始めた。かかる薄型の表示装置の中で、特にＰＤＰは大画面で直視型の表示デバイスとしての期待が非常に大きい。
【０００３】
ところで、ＰＤＰの中間調表示方法の一つとして、フィールド内時間分割法がある。この中間調表示方法は１フィールドを輝度の重みの異なるＮ枚の画面（これらを以下サブフィールドと呼ぶ）で構成する。輝度の重みの小さい側からＳＦ０，ＳＦ１，ＳＦ２，…・・，ＳＦ（Ｎ−１）と呼ばれ、それらサブフィールドの輝度の重みの比はそれぞれ、２0，２1，２2，…・・，２N-1である。１フィールド内の中間輝度は、これらのサブフィールドの発光の有無を選択する事により行われ、人間の視覚特性（残光特性）により、人間の目に感じる輝度は発光サブフィールドの各々の輝度の和で表せる。この中間調表示方法で表現出来る中間調数は１フィールド内のサブフィールド数、即ち２N通りである。
【０００４】
図２１に上記中間調表示方法を用いた１フィールド内の表示シーケンスを示す。１フィールドは８枚（Ｎ＝８）の輝度の重みの異なるサブフィールドで構成され、輝度の重みの大きいほうからＳＦ７，ＳＦ６，…・・，ＳＦ０と呼ばれている。ここで、ＳＦ７を最上位ビット（ＭＳＢ）側、ＳＦ０を最下位ビット（ＬＳＢ）側と呼んでいる。各々のサブフィールドは１フィールドの中に、ＳＦ０，ＳＦ１，…・・，ＳＦ７と輝度の重みの小さいものから順に並んで発光を制御する場合が多く用いられている。即ち、各サブフィールドの発光回数の比は、ＳＦ０を“１”とすると、ＳＦ１は“２”，ＳＦ２は“４”，…・・ＳＦ６は“６４”，ＳＦ７は“１２８”である。このサブフィールド数が８個のときは２５６階調まで表現できる。
【０００５】
ところで、上述したサブフィールド法による中間調表示方法は、１と０の２つの階調しか表現出来ないＰＤＰのような２値表示デバイスでも多階調表現が可能な技術として優れた方法であり、このサブフィールドを用いた表示方法により、テレビ画像とほぼ同様な画質がＰＤＰにおいても得られるようになった。
【０００６】
しかしながら、例えば、濃淡が緩やかに変化している被写体で動きのある映像が表示された場合、テレビでは見られないＰＤＰ画像に特有のいわゆる偽輪郭が発生する問題がある。この偽輪郭発生現象は視覚の特性からくるもので、その映像信号レベルが２５６階調表示のとき、上記の１２８，６４，３２，１６などといった２のＮ乗の境界付近に沿って、あたかも階調が失われたような状態で、更には本来表示すべき色と違った色が縞状となって見られる現象である。しかし、静止画像を表示した場合には偽輪郭は感じられない。動きのある部分でかつ上記信号レベルの周辺でのみ認知されるのが偽輪郭の特徴である。
【０００７】
図２２を用いてサブフィールド法による階調表示方法で偽輪郭が発生する原理について説明する。
図２２（ａ）では１フィールド内のサブフィールド数が８個でその配列が輝度の重みの小さい方、即ちＳＦ０，ＳＦ１，ＳＦ２，…，ＳＦ７の順に並ぶ場合を示している。ある画素位置の信号レベルが１２７から１２８に変化しているときに、この動画像が１フィールドで３画素移動している。図２２（ｂ）は、観測者が画面上でこの動画像を観測した結果を示す。
【０００８】
このように、信号レベル１２７（ＳＦ０からＳＦ６までの発光）と信号レベル１２８（ＳＦ７のみが発光）が隣り合っている場合、その階調差は１ＬＳＢ（１／２５６）であるが、人間の網膜上で感じる発光値はこの発光時間の不均一性により画像の移動した画素分だけ各々の信号レベルの発光が重なり合う結果、空間的拡がりとなって網膜上には大きな値(積分値)として感じる。即ち、本来同じ画素で発光しているべき各々のサブフィールドの発光が動画像部では異なった画素位置で発光していることになり、画素の中間調輝度が単に各サブフィールドの和で表現出来なくなる。これが偽輪郭として感じられる理由である。
【０００９】
図２２に示した様に、動画像が表示画面の左側から右側へスクロールすると、上述の信号レベルの境界部は明るい線として感じられ、反対に動画像が表示画面の右側から左側へスクロールすると、上述の信号レベルの境界部はサブフィールドの空間的分離として感じられ、暗い線として感じられることになる。一方、サブフィールドの配列が輝度の重みの大きい方、即ちＳＦ７、ＳＦ６，ＳＦ５，…，ＳＦ０と順に並んでいる表示方法においては、動画像が表示画面の左側から右側へスクロールすると、信号レベルの境界部は暗い線として感じられ、反対に動画像が表示画面の右側から左側へスクロールすると、信号レベルの境界部は明るい線として感じられることになる。つまり、表示画面の動画像の移動方向によって、偽輪郭の見え方が異なることになる。
【００１０】
更に、この偽輪郭の発生は動画像の動き速度にも依存し、動き速度が速い程、偽輪郭の及ぶ範囲は大きい。例えば、１フィールド中に１０画素移動する動画像では偽輪郭の及ぶ画素幅は１０画素にも及ぶ。
【００１１】
従来より、この偽輪郭に対する対策として各種の提案がされており、特開平７−２７１３２５号公報では、サブフィールドの表示順を、パルス数比率が１，２，４，８，１６，３２，６４，１２８のような単純増加でなく、偽輪郭が目立たないような順序に並び替える技術を開示している。例えば、サブフィールド中で一番表示期間の長いサブフィールドをフィールドの中心に配置するような順番に表示したり、１フィールド毎に表示順番を変化させたりする方法などである。
【００１２】
しかし、サブフィールドの並び替えやフィールド毎にサブフィールドの発光シーケンスを変えたのでは、動きの速い動画像に対する偽輪郭には対応出来ないなど効果が極めて限定される。
【００１３】
また、特開平８−１２３３５５号公報には動き検出を利用した偽輪郭の抑止技術が開示されている。これは連続した２枚のフィールド画面の動画と背景画に対応した画面から動き量と方向を求め、この検出値と各サブフィールド画面の単位時間における分割時間割合に基づいて動き補正量を求め、この補正量だけ対応するサブフィールド画面の発光パターンを移動させるといった技術である。
【００１４】
また、特開平８−２１１８４８号公報では、動きベクトルをフィールド間の表示データにより画素ブロック毎に検出し、フィールド内の先頭サブフィールドは入力データのそれに対応するデータを表示し、それに続く各サブフィールドは各々の先頭サブフィールドからの遅れ時間をフィールド周期で割った値を動きベクトルに掛け算した値を用いて表示データを移動させ画像を表示する技術内容が開示されている。
【００１５】
上記したような動き量に応じてサブフィールドの発光パターンを移動させたり表示データを変えるだけでは、動画像のブロック境界で不自然な階調変化や髭が発生し、かえって副作用が目立つ場合がある。また、視覚光量とのマッチングも難しく、動き量により単にデータの移動だけでは対応出来ないことが視覚実験により判った。また、動き検出を利用した偽輪郭抑制では、偽輪郭を防止するための決め手は如何に精度良く動き量を検出するかにかかっているにも拘わらず、実用性のある動き検出の具体的構成が十分に開示されていない。
【００１６】
また、特開平８−２３４６９４号公報に開示された偽輪郭補正方法では、同一画素に対応する単位画素データを少なくとも１フレーム期間だけ離れた前回値と今回値とを比較して、双方の発光論理値の最上位ビットの桁位置がお互いに異なる場合、今回値について補正データを加算あるいは減算している。
【００１７】
しかし、この偽輪郭補正方法では動画像の動き方向が特定出来なければ効果が反対になる可能性がある。例えば、上位方向にビットの桁位置を検出した時補正データを減算するが、前述したように画像が左方向に動いているときに上記演算を行うとかえって偽輪郭が強調され効果が反対になる場合が生じてくる。同様に下位方向にビットを検出した時補正データを加算するが、画像が反対方向に移動している場合は効果が逆となる。また、速度の速い動画像にも対応できない欠点を有している。
【００１８】
【発明が解決しようとする課題】
このように、従来からある偽輪郭の抑止に関する技術では、動きの速い動画像に発生する偽輪郭の防止には十分に対応することができなかった。
【００１９】
本発明は、以上のような実情に鑑みてなされたもので、サブフィールド法で階調表示を行う表示装置において、動画を目で追従したときの偽輪郭の発生を大幅に抑制できる動画像表示方法及び動画像表示装置を提供することを目的とする。
【００２０】
【課題を解決するための手段】
上記課題を達成するために、本発明は以下のような手段を講じた。
【００２１】
請求項１記載の発明は、偽輪郭の発生する信号レベル近傍の画素の動きを捉え、画素の動きに応じて現画像データを補正する構成を採る。
この構成によれば、偽輪郭の発生する信号レベル近傍で画素の動きを捉えることで、画素値の分布の偏りに影響されることなく、画像の局所的な変化を高速に検出することができる。偽輪郭部の発生する画素の動きは極めて局所的な検出で十分であり、計算時間、回路構成等を簡単にできる。
【００２２】
請求項２記載の発明は、偽輪郭の発生する信号レベル近傍の区間をしきい値により現フィールド画像及び前フィールド画像を２値化し、２値化画像を比較して動いた画素の移動画素数及び移動方向を検出し、現フィールド画像で偽輪郭の発生する信号レベルを持った動き画素をその移動画素数及び移動方向に応じて補正する構成を採る。
この構成によれば、偽輪郭の発生する信号レベル近傍のしきい値にて現フィールド画像及び前フィールド画像を２値化するので、少領域でも原画像の特徴を反映したものとなり、偽輪郭部の発生する画素の動きを検出できる。また、動きの検出された画素の移動画素数及び移動方向に応じて補正するので偽輪郭の発生を防止できる。
【００２３】
請求項３記載の発明は、現フィールド２値化画像を複数の検出ブロックに分割し、検出ブロック毎に前フィールド２値化画像内に参照領域を設定し、参照領域内に複数設定した参照ブロックと検出ブロックとの一致度を評価し、評価値の最も高い参照ブロックと前記検出ブロックとの位置関係から動き画素の移動画素数及び移動方向を検出する構成を採る。
この構成によれば、偽輪郭が発生しない画像の動きベクトルは無視することができ、偽輪郭の発生する近傍の２値画像で動きを検出することから、例えば８ビットで表された画素値（多値画像）を用いる場合に比べて減算、比較等の演算時間を大幅に減らすことができ、動きの速い動画像に対しても十分に対応できる。
【００２４】
請求項４記載の発明は、検出ブロックと参照ブロックとの比較で動きの検出されなかった検出ブロックを不定ブロックとして仮登録し、動きが検出された既知ブロックに囲まれた不定ブロックを既知ブロックと同じ移動画素数及び移動方向にする構成を採る。
この構成によれば、動きが検出された既知ブロックに囲まれた不定ブロックを既知ブロックと同じ動きがあったものとして扱うので、動画像のエッジ部を含む複数のブロックが１つ又は複数のブロックを囲む場合であっても、囲まれたブロックの動きを検出できる。
【００２５】
請求項５記載の発明は、不定ブロックの処理は、それを囲む既知ブロックの移動画素数及び移動方向により線形補間により、当該検出ブロックの移動画素数及び移動方向を求める構成を採る。
この構成によれば、動きが検出された既知ブロックに囲まれた不定ブロックをその近傍既知ブロックから線形補間により求めるので、動画像のエッジ部を含む複数のブロックが１つ又は複数の不定ブロックを囲む場合であっても、囲まれたブロックの動きをさらに高精度に検出できる。
【００２６】
請求項６記載の発明は、検出ブロック毎にＲ，Ｇ，Ｂの各色成分について移動画素数及び移動方向を検出し、多数決判定により当該検出ブロックの移動画素数及び移動方向を確定する構成を採る。
この構成によれば、多数決判定により当該検出ブロックの移動画素数及び移動方向を確定するので、そのブロックの動き判定が実現できる。
【００２７】
請求項７記載の発明は、各色成分の当該検出ブロックとその近傍ブロックの移動画素数及び移動方向を取り出す際に、各色成分毎に数値の大きい順に並び替えを行い、その中央値を代表値として求める構成を採る。
この構成によれば、当該検出ブロックの移動画素数及び移動方向を、各色成分毎に求められた移動画素数及び移動方向とそれを囲む近傍ブロックの移動画素数及び移動方向の複数の値を使用して求めるため、さらに高精度の動き判定が実現できる。
【００２８】
請求項８記載の発明は、現フィールド画像及び前フィールド画像を、２のＮ乗の信号レベル又はこれら信号レベルを組み合わせた信号レベルを中心とした多階層のしきい値で２値化し、多階層画像データ毎に動きを検出する構成を採る。
この構成によれば、２のＮ乗の信号レベル境界付近で偽輪郭が発生するので、偽輪郭が発生する各信号レベル近傍の２値画像を抽出することができ、多階層画像データ毎に動きを検出するので各信号レベルで発生する偽輪郭をそれぞれ抑止できる。
【００２９】
請求項９記載の発明は、現フィールド画像及び前フィールド画像を２値化前に平滑化処理する構成を採る。
この構成によれば、現フィールド画像及び前フィールド画像のノイズを除去できるので、正確に動きを検出できる。
【００３０】
請求項１０記載の発明は、偽輪郭の発生する信号レベルに対応して移動画素数及び移動方向に応じた補正データを補正テーブルに登録しておき、検出した動き画素の移動画素数及び移動方向に基づいて前記補正テーブルから補正データを取り出す構成を採る。
この構成によれば、動き画素の移動画素数及び移動方向を検出でき、移動画素数及び移動方向に基づいて細かく補正データを設定した補正テーブルを使って補正することができる。
【００３１】
請求項１１記載の発明は、検出ブロック毎に検出した移動画素数及び移動方向を動きベクトルテーブルに登録する一方、現フィールド画像から偽輪郭の発生する信号レベルの座標データを検出して補正候補画素を抽出し、この補正候補画素が属する検出ブロックの移動画素数及び移動方向を前記動きベクトルテーブルから取りだし、取り出された移動画素数及び移動方向と前記偽輪郭発生信号レベルとに対応した補正データを前記補正テーブルから取り出す構成を採る。
この構成によれば、偽輪郭の発生する信号レベルを持ち、かつ動きの検出された画素について、信号レベル、動き画素数及び方向に応じた最適な補正データを取得することができる。
【００３２】
請求項１２記載の発明は、偽輪郭の発生する信号レベルの検出密度に基づいて補正候補画素を決定する構成を採る。
この構成によれば、視覚的に偽輪郭が確認される画素密度の領域にある画素だけ補正することができる。
【００３３】
請求項１３記載の発明は、偽輪郭の発生する信号レベル近傍のしきい値で現フィールド画像及び前フィールド画像を２値化する２値化処理手段と、フィールド間で２値化画像を比較して画素の動きを検出する動き検出手段と、偽輪郭の発生する信号レベル毎に移動画素数及び移動方向に対応した補正データが登録された補正テーブルと、現フィールド画像から偽輪郭の発生する信号レベルの画素を抽出する抽出手段と、この抽出画素の動き検出結果及び信号レベルに対応した補正データを前記補正テーブルから取り出して当該画素を補正する補正手段とを具備する構成を採る。
この構成によれば、ＰＤＰ等において偽輪郭が見えない動画像表示を確実に得ることができる。
【００３４】
【発明の実施の形態】
次に、本発明の実施の形態を図面を参照して詳細に説明する。
図１は、本発明の動画像表示方法を適用した階調画像表示装置の全体構成を示している。この階調画像表示装置は、ビデオ信号処理部１でビデオ信号をＲ，Ｇ，Ｂの各色成分に分離し、Ａ／Ｄ変換部２でＲ，Ｇ，Ｂの画像データに変換してから動き検出処理部３に入力する。動き検出処理部３で、平滑化処理、多階層２値化処理、ブロックマッチング処理、多数決・統合判定処理を実施して入力画像の動き画素数及び移動方向を検出する。動き検出処理部３で検出した動き画素数、移動方向及び検出ブロック情報を、動き画素数及び移動方向に応じた補正量が設定された補正テーブルを備えたデータ補正処理部４へ入力して偽輪郭の発生する画素の階調データを補正して出力処理部５へ出力する。出力処理部５は、各画素の階調データを電圧印可時間幅に対応したパルス数に変換してＸスキャンドライバ７及びＹスキャンドライバ８に与えることでＰＤＰで構成された画像表示部６に中間調表示を行っている。
【００３５】
なお、同期分離処理部９においてビデオ信号から同期信号を分離して入力ビデオ信号に同期したタインミング信号をタイミング発生部１０で生成して各部に与えている。
動き検出処理部３の具体的な構成を図２に示す。Ａ／Ｄ変換部２から入力した原画像データを、現フィールドブロックＢ１に入力するとともに、遅延回路２０で１フィールド分遅延させてから前フィールドブロックＢ２へ入力する。現フィールドブロックＢ１では、現フィールドの原画像データをＲ，Ｇ，Ｂの３つの平滑化フィルタからなる平滑化フィルタブロック２１−１に入力し、各色毎に平滑化処理して原画像に含まれているノイズ成分を除去する。平滑化処理した原画像データをＲ，Ｇ，Ｂ毎に設けられた多階層２値化処理部からなる多階層２値化ブロック２２−１に入力して各色毎に複数の２値化しきい値で２値化処理する。ここで、２値化しきい値は偽輪郭の発生する信号レベル近傍に設定する。例えば、８ビットで２５６階調の表示を行う場合であれば、３１から３２、６３から６４、１２７から１２８及びその逆の信号レベルを含んだ値に設定する。各しきい値毎に得られた２値化画像データ（多階層画像データ）はそれぞれ多階層メモリブロック２３−１の各色の画像メモリに格納する。
【００３６】
１フィールド前の原画像を処理する前フィールドブロックＢ２は上記現フィールドブロックＢ１と同様に構成されている。すなわち、平滑化フィルタブロック２１−２で１フィールド前の原画像データを各色毎に平滑化処理し、Ｒ，Ｇ，Ｂ毎に設けられた多階層２値化処理部からなる多階層２値化ブロック２２−２で各色毎に複数の２値化しきい値で２値化処理し、各しきい値毎に得られた多階層画像データを多階層メモリブロック２３−２の各色の画像メモリに格納する。
【００３７】
動き検出のためのブロックマッチング処理において、現フィールド画像から切り出す検出ブロックＫＢのアドレスマッピングをアドレスマッピング部２４−１が行い、前フィールド画像から切り出す参照ブロックＲＢのブロックのアドレスマッピングをアドレスマッピング部２４−２が行う。検出ブロックＫＢ及び参照ブロックＲＢの各画像データは動きベクトル検出部２５へ入力する。
【００３８】
動きベクトル検出部２５は、２値演算部が各色毎に設けられた２値演算ブロック２６と、比較部が各色毎に設けられた比較検出ブロック２７とで構成され、ブロックマッチングによってフィールド間の画像の動きを求める。
【００３９】
多数決・統合判定部２８は、各検出ブロックＫＢの動き画素数及び動き方向を判定して、判定結果を動きベクトルテーブル２９に登録する。
【００４０】
データ補正処理部４は、動きベクトルテーブル２９及び予め視覚実験に基づいて作り込んだ補正テーブルにより偽輪郭部の画素データを補正する。
【００４１】
以上のように構成された階調画像表示装置の動作の詳細について説明する。
図１０は、図２に示す動き検出処理部３の各ブロックでの処理内容を概念的に抜き出した図である。以下、図１０に示す動き検出処理部３の各ブロックでの処理内容について説明する。
【００４２】
現フィールド画像及び前フィールド画像のそれぞれにつき、偽輪郭の発生する信号レベルに設定した２値化しきい値により２値化処理がなされる。
【００４３】
ここで、多階層２値化画像の概念について図１１を参照して説明する。
図１１は横軸が画素位置、縦軸が画素値を示す。画素値の変域はｎ個のしきい値Ｔｈ１からＴｈｎ（同図ではｎ＝５）により、（ｎ＋１）個の区間に分割されており、隣接区間では値が異なるよう２値化画素値が割り振られ、各画素は画素値がどの区間に属するかに応じて２値化される。
【００４４】
この２値化方法は画素値の分布の偏りに影響されることなく、画像の局所的な変化を表すことができるので、２値化された画像は動きベクトル検出のブロックのような少量域内でも原画像の特徴を反映したものとなる。偽輪郭の発生する画素の動きは極めて局所的な検出で十分であり、計算時間、回路構成等を簡単にできる。
【００４５】
次に動きベクトル検出部２５でのブロックマッチング処理について説明する。フィールド間の画像の動きを求める方法としてブロックマッチング方法を用いている。この方法は、前述したように検出ブロックＫＢを水平１６画素、垂直１６画素に、参照ブロックＲＢも検出ブロックＫＢと同じサイズに設定したとき、参照ブロックＲＢを含む参照領域Ｒを水平４８画素、垂直４８画素とすれば、水平方向（ｘ方向）の動き画素数は−１６画素から＋１６画素まで検出でき、垂直方向（ｙ方向）の動き画素数も同様に−１６画素から＋１６画素まで検出できる。動画像が表示画面で動くとき、その偽輪郭が目立つのは、そのフィールド間の動きが６画素から１５画素近傍であるので、動き画素数の検出能力は上記の程度必要である。図１２に示すように、現フィールドと前フィールド間の画像の動きを矩形のブロック単位で求め、現フィールドを複数の参照領域Ｒに分割した各検出ブロックＫＢについて最もよく一致する部分を前フィールドの参照領域Ｒの中に含まれる参照ブロックＲＢの中から見つけ、この間の量を動きベクトルとして求める。
【００４６】
この際、ブロックの一致の度合いは定められた評価関数値の大小により判定し動きベクトルの検出は参照領域Ｒに含まれる多数の参照ブロックＲＢの中から評価関数最小値を与えるものを探し出すことにより行う。
【００４７】
上記方法に基づいた構成を濃淡画像で行う場合、評価関数値の計算に多値（例えば８ビット）で表された画素値を用いて減算、比較等の多量の演算を行うことにより実時間処理で評価出来る参照ブロックＲＢの数が限られているため、正しい動きベクトルを与える参照ブロックＲＢが評価されない場合が生じ、結果として動きベクトルの検出精度が低くなるという問題を有していた。
【００４８】
ＰＤＰにおける偽輪郭部抑制のためには、偽輪郭の発生する信号レベルが判っているため、偽輪郭の発生する近傍の画像の動きを求めるためには、この近傍の画像を２値化し、その２値画像の動きベクトルを求めることにより、計算時間及び精度が向上する。また、偽輪郭が発生しない画像の動きベクトルは無視することができる。
【００４９】
２値演算ブロック２６において、一致度を示す評価関数を計算する。図３に、検出ブロックＫＢサイズが４＊４画素の場合の２値演算部の構成を示す。なお、検出ブロックＫＢサイズは１６＊１６、参照領域Ｒは４８＊４８で説明するが、ここでは説明を簡単にするために４＊４で説明する。
【００５０】
２値演算部は、ブロック内の４ラインについてフィールド間のマッチングを検出する２値ラインマッチング部３０−１、３０−２、３０−３、３０−４を備えている。２値ラインマッチング部３０−１は、現フィールドと前フィールドとの間で第１ラインのラインマッチングをみる。検出ブロックＫＢの第１ラインのラインデータをシフトレジスタ３１−１に入力し、参照ブロックＲＢの第１ラインのラインデータをシフトレジスタ３１−２に入力する。シフトレジスタ３１−１、３１−２は、４つの画素データを保持することができ、２つのシフトレジスタから同一画素位置の画素データが同一の排他的論理和回路３２−１、３２−２、３２−３、３２−４へ入力するようにしている。排他的論理和回路の出力をビット加算回路３３で加算してマルチプレクサ３４へ入力する。
【００５１】
２値ラインマッチング部３０−２、３０−３、３０−４も、上記２値ラインマッチング部３０−１と同じ構成を有しており、それぞれ第２から第４ラインのラインマッチングを評価してライン評価値となるビット加算値をマルチプレクサ３４へ入力する。マルチプレクサ３４を介して２値ラインマッチング部３０−１、３０−２、３０−３、３０−４の出力を選択的に加算器３５へ入力して４ラインのライン評価値を加算した加算値を１ブロック分のブロック評価値としてレジスタ３６に保存する。レジスタ３７は出力制御のためのものである。
【００５２】
上記２値演算部において、多階層画像メモリから信号線を介して供給された検出ブロックＫＢ、参照ブロックＲＢのデータは走査線単位で２値ラインマッチング部３０に送られる。各ラインマッチング部３０−１から３０−４では、シフトレジスタ３１−１、３１−２を用いて画素毎の２値データが取り出され、ブロック内で同位置の画素同士で排他的論理輪回路３２−１から３２−４により一致、不一致が評価され、不一致のとき値１が、一致のとき値０がビット加算回路３３に供給される。
【００５３】
ビット加算回路３３ではこれらの和が求められる。このビット加算値は走査線毎の不一致画素数を示す。この和が信号線を介して２値ラインマッチング部３０より出力され、マルチプレクサ３４に供給される。マルチプレクサ３４では、選択制御線より送られた選択制御信号により各ラインのビット加算回路３３の出力を順次選択していき、信号線を介して加算器３５に供給する。加算器３５及びレジスタ３６によりこの入力値の和が求められ、ブロック間の不一致画素数が評価関数値として求められる。
【００５４】
比較検出ブロック２７の動作について説明する。
図４に比較検出ブロック２７の比較部の構成を示している。比較部は、最小評価関数値を保持するレジスタ４１、最小評価関数値を示した参照ブロックＲＢのシフト量を保持するレジスタ４２、現在の最小評価関数値と今回比較対象の２値演算部出力（あるシフト量のブロック評価関数値）とを比較する比較器４３、マルチプレクサ４４、４５で構成されている。
【００５５】
比較部において、レジスタ４１に各時点における評価関数値の最小値が保持され、レジスタ４２にその最小値に対応する参照ブロックＲＢのシフト量が保持される。比較器４３では信号線を介して供給される評価関数値とレジスタ４１から供給される評価関数最小値とが比較され、その比較結果が信号線を介して２つのマルチプレクサ４４、４５へ選択制御信号として送出される。２値演算部からの入力が小のとき、マルチプレクサ４４は２値演算部からの入力評価関数値で、レジスタ４１の内容を更新し、かつマルチプレクサ４５は２値演算部からの入力シフト量でレジスタ４２の内容を更新する。最終的にレジスタ４２に保持されているシフト量がこの検出ブロックＫＢの動きベクトルとして多数決・統合判定部２８へ送出される。
【００５６】
多数決・統合判定部２８の動作について説明する。
ここでは各Ｒ，Ｇ，Ｂ成分の検出ブロックと参照ブロックとの比較で動きの検出されなかった検出ブロック、即ち、不定ブロックを周囲の既知ブロックの情報から演算する処理と各Ｒ，Ｇ，Ｂ成分の検出ブロックの動きベクトル情報から一つの動きベクトルに統合演算する処理を行う。
【００５７】
多数決・統合判定部２８に動きベクトル検出部２５が現フィールドを複数分割した各検出ブロックＫＢについて検出した動きベクトル情報が入力する。多数決・統合判定部２８は、不定ブロック処理及び多数決処理により各検出ブロックＫＢの動き画素数及び動き方向を判定して、判定結果を動きベクトルテーブル２９に登録する。
【００５８】
多数決・統合判定部２８で実行される不定ブロック処理について図１３から図１５を参照して説明する。
【００５９】
多階層２値画像をブロック毎に区切った場合、そのエッジ部がブロック内にある場合は動き画素数を求めることができるが、ブロック内の画像データが全て“１”の領域およびその領域が隣接している場合（図１３のＦの符号を付した領域）は、そのブロックの領域の動き画素数が見つけられない。
【００６０】
このような場合は、動きベクトルの検出時にこのブロック域での動き画素数を不定（例えば、フラグを１）とし、他のブロックと区別して仮登録しておく。そして、既知ブロックＫＴ領域で挟まれた不定ブロックＦＴ領域の動き画素数及び方向を既知ブロックＫＴから求める。不定ブロックＦＴ領域は既知ブロックＫＴ領域と同じ動きをするから、これらを囲む既知ブロックＫＴの動き画素数と同一値を採用する。
【００６１】
この方法として、イメージマップ法が考えられる。例えば、２値画像サイズが水平６４０画素、垂直４８０画素のいわゆるＶＧＡタイプの場合、動きベクトルの検出ブロックＫＢが１６＊１６画素であると、図１４に示すようにブロック数は水平４０、垂直３０の計１２００個のブロック域に分割できる。このため、４０＊３０画素のイメージマップとして不定ブロックＦＴを囲む既知ブロックＫＴの形状からその不定ブロックＦＴを求めることができる。ここで、不定ブロックＦＴを２値画像でデータ“１”の領域とし、既知ブロックＫＴをデータ“０”の領域とする。
【００６２】
この方法では、図１５に示す３＊３のエッジ検出ウィンドウオペレータを用いる。図１５には３＊３のエッジ検出パターンの組み合わせのうち９例を示している。このウィンドウで４０＊３０のブロックイメージを走査し、注目点を基準にしてエッジパターンと一致したところがあれば、それは不定ブロックＦＴ域と既知ブロックＫＴ域の境界であるから、その不定ブロックＦＴを既知ブロックＫＴに置き換えていく。その不定ブロックＦＴの動き画素数は注目点を含む既知ブロックＫＴの動き画素数と等しくなる。
【００６３】
次に、さらに、高精度な方法である線形補間法について説明する。
図１６はこの方法による不定ブロックを含んだブロック関係を示す。この場合の線形補間の手順は次のように行われる。
【００６４】
図１６において、まず注目ブロック（図中の＊印）を基準として右方向へフラグ’−１’の既知ブロックの探索を行い、このブロックが存在すればその動きベクトルを取り出し、参照ブロック１とする。この時、参照ブロック１の位置の動きベクトルを（ｘ１，ｙ１）で示し、注目ブロックからの距離をｄ１とする。なお、動きベクトルは、動き画素数と動き方向（＋／−）で示すものとする。
【００６５】
次に、注目ブロックを基準として左方向へフラグ’−１’の既知ブロックの探索を行い、このブロックが存在すればその動きベクトルを取り出し、参照ブロック２とする。この時、参照ブロック２の位置の動きベクトルを（ｘ２，ｙ２）で示し、注目ブロックからの距離をｄ２とする。
さらに、注目ブロックを基準として上方向へフラグ’−１’の既知ブロックの探索を行い、このブロックが存在すればその動きベクトルを取り出し、参照ブロック３とする。この時、参照ブロック３の位置の動きベクトルを（ｘ３，ｙ３）で示し、注目ブロックからの距離をｄ３とする。
最後に、注目ブロックを基準として下方向へフラグ’−１’の既知ブロックの探索を行い、このブロックが存在すればその動きベクトルを取り出し、参照ブロック４とする。この時、参照ブロック４の位置の動きベクトルを（ｘ４，ｙ４）で示し、注目ブロックからの距離をｄ４とする。
【００６６】
このように、左右上下方向の探索により取り出した参照ブロックと注目ブロックとの距離から、各参照ブロックへの重みを算出し、その重みと参照ブロックの動きベクトルを用いて、以下に示す線形補間による演算式に従って注目ブロックの動きベクトルを求める。
【００６７】
図１６から、各参照ブロックへの重みｗは、次式で与えられる。
【００６８】
【数１】
0004203136.tif 000002
そして、線形補間によって求める注目ブロックの動きベクトル（ｍｘ，ｍｙ）は
【００６９】
【数２】
0004203136.tif 000003
となる。
【００７０】
ここで、
【００７１】
【数３】
0004203136.tif 000004
である。これをブロック個数分繰り返し行い、不定ブロックの動きベクトルが線形補間によって求められる。
【００７２】
図１７は、線形補間の図形的位置関係を示す。
ここで、ｉ，ｊ，ｋ，ｌは注目点Ｐからの距離で、Ｐｉ，Ｐｊ，Ｐｋ，ＰｌはＰからそれぞれｉ，ｊ，ｋ，ｌだけ離れた点の値とすると、注目ブロック位置Ｐの補間式は以下ようになる。
【００７３】
各点への重みをＩとすると、
【００７４】
【数４】
0004203136.tif 000005
で与えられるから、求める点の値（Ｐｘ、Ｐｙ）は
【００７５】
【数５】
0004203136.tif 000006
となる。ここで、
【００７６】
【数６】
0004203136.tif 000007
である。
【００７７】
以上のように、Ｒ，Ｇ，Ｂ画像毎に不定ブロックＦＴ領域の動きベクトルを周辺の既知ブロックＫＴの動き画素数から演算式により求める。そして、次にＲ，Ｇ，Ｂそれぞれの動きベクトルから多数決判定などの統合処理により各ブロックの動きベクトル値を求める。
【００７８】
ここで、図１８に示すフローチャートを参照して多数決処理による統合処理の具体的な内容を説明する。
【００７９】
各ブロックはＲ，Ｇ，Ｂ画像について同じブロックであるため、当然動き方向及び動き画素数は同じはずである。しかしながら、対象画像のブロックにおける２値演算及び比較部の演算誤差などにより差違が発生する場合がある。
【００８０】
そこで、ブロック毎のｘ、ｙ方向の動き画素数を取り込み（Ｓ1）、Ｒ、Ｇ、Ｂ毎にｘ、ｙ移動方向が一致しているかどうか判定する（Ｓ２）。この場合、移動方向の補正・選定では（Ｓ３）、Ｒ，Ｇ，Ｂ各ブロックの２値演算、比較検出ブロックから出力されたｘ、ｙの動き方向が第１象限から第４象限までの座標軸を基準としてその符号が付与されているため、これがＲ，Ｇ，Ｂで同じかどうか判断する。これは、Ｒ，Ｇ，Ｂの動き方向のｘ方向とｙ方向のそれぞれについて、該当のブロックとそれに隣接するブロックの動きの符号を含めて、多数決判定し動き方向として動きベクトルテーブル２９に登録する（Ｓ６）。このような手段により動き方向の精度を向上させる。
【００８１】
さらに、画素数の補正・選定（Ｓ５）も同様にＲ，Ｇ，Ｂ各ブロックの２値演算、比較検出ブロック２６、２７から出力されたｘ、ｙの動き画素数もＲ，Ｇ，Ｂで同じかどうか判定する（Ｓ４）。この場合、Ｒ，Ｇ，Ｂの動き画素数が全く違っていれば、お互いが近い画素数で、かつその値の中央値から順に重み付けした値の平均値を採用し、動き画素数として動きベクトルテーブル２９に登録する（Ｓ６）。この平均値は通常四捨五入した値をとる。
【００８２】
例えば、ｘ方向の動き画素数が、Ｒ，Ｇ，Ｂそれぞれ６、５、２であるとし、その重みを２、３、１とする。ｘ方向の動き画素数は、（２＊６＋３＊５＋１＊２）を（２＋３＋１）で割った値、すなわち５画素とする。
【００８３】
また、Ｒ，Ｇ，Ｂについて同じ動き画素数値が２つ存在した場合も同様である。例えば、ｘ方向の動き画素数が、Ｒ，Ｇ，Ｂそれぞれ２、５、２であるとし、その重みは３、１、２であるから、ｘ方向の動き画素数は（３＊２＋１＊５＋２＊２）を（３＋１＋２）で割った値、すなわち３画素とする。ｙ方向についても同様に行えばよい。
【００８４】
さらに、動きベクトルの精度を向上させるための統合処理方法を以下に説明する。
【００８５】
図１９はＲ，Ｇ，Ｂの各動きベクトル値から、一つの動きベクトルを求めるための動作フローチャートである。
【００８６】
前述の図１８に示した統合処理では、Ｒ，Ｇ，Ｂ画像の各注目ブロックのみによる演算処理のため、得られた動きベクトルの値が精度の点で実際の値と一致しないブロックが発生する場合があり、この検証結果では、精度向上のために注目ブロックが不定ブロックなのか、更には既知ブロックなのかということと、注目ブロックの周辺ブロックの動きベクトルをも考慮に入れる必要があるとういうことが示唆された。
【００８７】
そこで、図１９の動作フローチャートで示すように、この統合処理では注目ブロックのフラグをまずチェックし（Ｓ１０）、フラグが’１’（Ｓ１１）と’０’（Ｓ１２）の場合は注目ブロックの動きベクトルを取り出さず、フラグが’−１’（Ｓ１３）の場合にのみ注目ブロックの動きベクトルを取り出す（Ｓ１４）。しかる後に、注目ブロックを囲む周辺８近傍の動きベクトルを取り出し（Ｓ１５）、取り出した複数の動きベクトルをｘ，ｙ成分毎にそれぞれ並び替えを行い、その中央値を代表値とする（Ｓ１６）。
【００８８】
これを統合すべきＲ，Ｇ，Ｂ数まで行い（Ｓ１７）、その後代表値の個数に応じて処理を分ける。代表値は各Ｒ，Ｇ，Ｂ成分毎に示されるものであり、この時、代表値が１個の場合には（Ｓ１８）、統合結果として注目ブロックの動きベクトル値はフラグ＝−１として、代表値そのままとする（Ｓ２３）。又、対象ブロックが偶数の場合など代表値が２個の場合は（Ｓ１９）、２つの代表値の平均値を求め（Ｓ２１）る。統合結果として注目ブロックの動きベクトル値はフラグ＝−１として、その平均値演算結果を採る（Ｓ２４）。さらに、代表値が３個の場合は（Ｓ２０）、３つの代表値の平均値を求め（Ｓ２２）る。統合結果は注目ブロックの動きベクトル値をフラグ＝−１として、その平均値演算結果を採る（Ｓ２４）。この処理を全ブロック個数まで求め（Ｓ２６）、注目ブロックの動きベクトル値とする。
【００８９】
さらに、図２０には、一例として注目ブロックとその８近傍のブロックの動きベクトルの値から一つの動きベクトルを求める統合処理の例を示す。
【００９０】
図２０（ａ）はＲ，Ｇ，Ｂそれぞれの動きベクトルについて、注目ブロックとその８近傍ブロック動きベクトル値を取り出した時の例を示す。この時、網線で記したブロックが注目ブロックである。又、図中（−、−）で記したブロックは動きベクトルが検出されない領域（フラグ＝０）のブロックを示す。
図２０（ｂ）はＲ，Ｇ，Ｂベクトルから動きベクトルの有効なブロックを取り出した時を示す。
【００９１】
図２０（ｃ）はＲ，Ｇ，Ｂそれぞれについて数値の大きい順に並び替える手順を示す。
図２０（ｄ）はＲ，Ｇ，Ｂそれぞれの代表値を算出する例を示す。その際、動きベクトルのデータ数が奇数個の場合は配列の中央値を選択し、偶数個の場合は中央２個の平均値を算出する手順を示す。このとき、動きベクトル値は小数点以下は切り捨てとする。このようにして、各Ｒ，Ｇ，Ｂの代表値が求められる。
【００９２】
この結果を用いて、図１９に示した様に統合処理を行う。つまり、代表値が１つの場合は、そのままの値を結果とし、代表値が２つの場合は、２つの値で平均値を採る。代表値が３つの場合は、その３つの値で平均値をとる。上記の平均値演算はいずれも小数点以下は四捨五入を行う。
【００９３】
この場合の動きベクトル結果は次式で与えられる。
【００９４】
【数７】
0004203136.tif 000008
従って、Ｘ方向の動きベクトルは’７’、Ｙ方向の動きベクトルは’１’と求まる。
【００９５】
上記統合処理の結果では、６４０画素ｘ４８０画素のＶＧＡタイプでは、検出ブロックが１６画素ｘ１６画素であるため、総検出ブロック数は１２００個である。この統合処理によって得られた平均誤差は先述の方式に比べ約２０％改善され、そのバラツキも極端な動きベクトルを示すことはなく、全体的な動きベクトルが正しく反映した結果が得られた。
【００９６】
以上のようにＲ，Ｇ，Ｂ画像毎に求められた動き画素数を各ブロック毎に１つの値に統合処理し、動き画素数として動きベクトルテーブルに登録する。
【００９７】
次に、補正処理ブロック４の処理内容について説明する。
【００９８】
図５にデータ補正処理部４の構成を示している。動きベクトルテーブル２９には現フィールドの各検出ブロックＫＢの動きベクトル情報が登録されている。
【００９９】
図６に動きベクトルテーブル２９の構成を示している。図６は、例えば水平６４０画素、垂直４８０画素の表示ができるＶＧＡタイプの場合の動きテーブルの例を示している。この場合、検出ブロックが水平１６画素、垂直１６画素であれば、合計１２００個のブロックに分割され、その各々のブロックについて動き画素数と方向がテーブル化されている。ここで、各ブロックの左上を座標の原点としている。現フィールド画像について動き検出の済んだ検出ブロック順に、検出ブロックＫＢのブロック番号と、その検出ブロックＫＢの原点からのオフセット値と、その検出ブロックＫＢの前フィールドからの動き画素数（動き方向を含む）とを登録している。
【０１００】
図７に補正テーブル５３の部分的な構成を示しており、図８及び図９に具体的なテーブル構成を示している。補正テーブル５３は、検出ブロックＫＢの動き方向を第１象限から第４象限に分け補正データを登録している。第１象限には、水平方向となるｘ方向（＋方向）の動き画素数に応じて各補正データが登録され、第２象限には、垂直方向となるｙ方向（＋方向）の動き画素数に応じて各補正データが登録されている。また、第３象限には、水平方向となるｘ方向（−方向）の動き画素数に応じて各補正データが登録され、第４象限には、垂直方向となるｙ方向（−方向）の動き画素数に応じて各補正データが登録されている。
【０１０１】
図８はｘ方向への動きが＋及び−方向、即ち、第１象限と第３象限への動きに対する補正データが、その信号レベルの変化とその位置に対応してマトリックス表になっている。
同様に図９はｙ方向への動きが＋及び−方向、即ち、第２象限と第４象限への動きに対する補正データが、その信号レベルの変化とその位置に対応してマトリックス表になっている。本説明では、補正データは記入していない。
【０１０２】
補正処理ブロック４では、現フィールドの原画像データを偽輪郭発生レベル検出部５１へ入力して、現フィールド画像の中から各画素の信号レベルにより偽輪郭が発生すると予測される偽輪郭発生予測点のｘ、ｙ座標を検出する。上述したように、偽輪郭の発生する信号レベルは、２のＮ乗の信号レベル境界付近であるので、２のＮ乗の信号レベル境界付近の画素を抽出する。
【０１０３】
偽輪郭発生レベル検出部５１で検出した全ての偽輪郭発生予測点を偽輪郭発生画素候補検出部５２に供給し、視覚特性にあわせて実際の表示で偽輪郭として認識される画素位置を偽輪郭発生予測点の中から特定して偽輪郭発生画素候補とする。実際の表示で偽輪郭として認識される画素位置を選定するに当たり、偽輪郭発生予測点の発生密度を反映させる。
【０１０４】
この場合の処理として、偽輪郭発生予測点のｘ，ｙ座標を中心に３＊３又は５＊５のウインドウオペレータによりその候補点の密度を検出することが可能である。この検出点の数により候補点が１点で孤立している場合に補正処理を行うかどうかを判定することができる。
【０１０５】
偽輪郭発生画素候補に現実に偽輪郭が発生するのは、当該偽輪郭発生画素候補の属する検出ブロックＫＢが実際に移動している場合である。
【０１０６】
偽輪郭発生画素候補の属する検出ブロックＫＢが移動ブロックの場合は、動きベクトルテーブル２９から検出ブロック番号、ｘ、ｙの動き画素数を取り出して補正テーブル５３へ供給する。補正テーブル５３には動き画素数に応じた補正量が登録されている。補正テーブル５３は予めＰＤＰ装置を用いて視覚実験によりその補正量を決めて、前述した図８、９の構造に従ってテーブル化しておく。
【０１０７】
偽輪郭発生画素候補検出部５２から出力される偽輪郭発生レベルと動きベクトルテーブル２９から取り出されるｘ、ｙの動き画素数とから特定される補正データを補正テーブル５３から取り出して補正演算部５４へ与える。補正演算部５４は、原画像データを補正するに当たり、偽輪郭発生レベルで且つ動きのあった画素は補正テーブル５３から動き画素数及び動き方向に応じた補正データが与えられるので視覚光量に合わせた信号レベルの補正する。
【０１０８】
以上のように、本発明の構成では動き検出処理部とデータ補正処理部と分割されており、各処理部は従属関係を持っているので、柔軟な補正処理が可能である。
【０１０９】
【発明の効果】
以上の説明から明らかなように、本発明によれば、サブフィールド法で階調表示を行う表示装置において、動画を目で追従したときの偽輪郭の発生を大幅に抑制できる動画像表示方法及び動画像表示装置を提供できる。
【図面の簡単な説明】
【図１】発明の一実施の形態に係る階調画像表示装置の全体構成図
【図２】上記一実施の形態に係る階調画像表示装置における動き検出処理部の機能ブロック図
【図３】動き検出処理部における２値演算部の回路構成図
【図４】動き検出処理部における比較部の構成図
【図５】データ補正処理部の構成図
【図６】動きベクトルテーブルの構成図
【図７】補正テーブルの構成図
【図８】補正テーブルにおける第１象限、第２象限のテーブル構成図
【図９】補正テーブルにおける第３象限、第４象限のテーブル構成図
【図１０】上記一実施の形態に係る階調画像表示装置の全体的な処理の流れを示す図
【図１１】サンプル画像での多階層画像としきい値との関係を示す図
【図１２】ブロックマッチングの概念図
【図１３】既知ブロックＫＴ処理の概念図
【図１４】不定ブロック処理の既知ブロックと不定ブロックの具体例を示す図
【図１５】エッジ検出ウィンドウの具体例を示す図
【図１６】不定ブロック処理のブロック関係を示す図
【図１７】不定ブロック処理における線形補間の位置関係を示す図
【図１８】多数決判定処理のフローチャート
【図１９】統合処理の動作フローチャート
【図２０】統合処理の演算を示す図
【図２１】サブフィールドの輝度比を示す図
【図２２】サブフィールド法での偽輪郭発生原理を示す図
【符号の説明】
１  ビデオ信号処理部
２  Ａ／Ｄ変換部
３  動き検出処理部
４  データ補正処理部
５  出力処理部
６  画像表示部
７  Ｘスキャンドライバ
８  Ｙスキャンドライバ
２０  遅延回路
２１−１、２１−２  平滑化フィルタブロック
２２−１、２２−２  多階層２値化ブロック
２３−１、２３−２  多階層メモリブロック
２４−１、２４−２  アドレスマッピング部
２５  動き検出ブロック
２６  ２値演算ブロック
２７  比較検出ブロック
２８  多数決・統合判定部
２９  動きベクトルテーブル
３０−１、３０−２  ２値ラインマッチング部
３０−３、３０−４  ２値ラインマッチング部
３４  マルチプレクサ部
３５  加算器
３６、３７  レジスタ
５１  偽輪郭発生レベル検出部
５２  偽輪郭発生画素候補検出部
５３  補正テーブル
５４  補正演算部
【図１】
0004203136.tif 000009
【図２】
0004203136.tif 000010
【図３】
0004203136.tif 000011
【図４】
0004203136.tif 000012
【図５】
0004203136.tif 000013
【図６】
0004203136.tif 000014
【図７】
0004203136.tif 000015
【図８】
0004203136.tif 000016
【図９】
0004203136.tif 000017
【図１０】
0004203136.tif 000018
【図１１】
0004203136.tif 000019
【図１２】
0004203136.tif 000020
【図１３】
0004203136.tif 000021
【図１４】
0004203136.tif 000022
【図１５】
0004203136.tif 000023
【図１６】
0004203136.tif 000024
【図１７】
0004203136.tif 000025
【図１８】
0004203136.tif 000026
【図１９】
0004203136.tif 000027
【図２０】
0004203136.tif 000028
【図２１】
0004203136.tif 000029
【図２２】
0004203136.tif 000030