(57)【要約】【課題】ブロックノイズに対する抑圧度を向上させながら、信号成分の損失を少なくすることができる画像信号処理装置及び画像信号処理方法を提供する。【解決手段】制御信号生成手段３０は、輪郭検出フィルタ１０からの輪郭検出信号と、ばらつき量検出手段２０からの標準偏差とを用いて画質制御値Ｇを算出し、ブロックノイズ調整手段５０に供給する。ブロックノイズ調整手段５０は、画質制御値Ｇを所定倍したシフト値が所定範囲内の場合にはパターン値を「１」とし、その他の場合には「０」とする。ブロックノイズ調整手段５０は、マッチング数を合計してブロック境界を特定し、ブロック境界で隣接する画素のエッジ変化量に応じた基準係数を特定し、これからフィルタ係数ｋを特定する。第１制御部６０は、平滑化フィルタ４０からの出力信号に対してフィルタ係数ｋを乗算した信号と、入力画像信号に対して「１−ｋ」を乗算した信号とを加算した信号を出力する。【選択図】図１
2009118081.tif 000002
【特許請求の範囲】
【請求項１】
要素画素が２次元に並べられたブロックから構成された入力画像信号を取得し、この入力画像信号の高周波数成分を平滑化した平滑化信号を生成する平滑化手段と、
前記入力画像信号における輪郭成分を検出する輪郭検出手段と、
前記ブロック内の各要素画素のばらつき量を算出する分散量検出手段と、
前記ブロック内の中心画素について、前記輪郭成分と前記ばらつき量との比較に基づいて画質制御値を算出する制御信号生成手段と、
前記画質制御値を、画質制御値範囲の中間値を含む所定の基準範囲と比較することにより２値化した画素のパターン値を算出する手段と、
前記パターン値のマップとブロック規定パターンとの比較により、前記ブロック規定パターンの周期的配置を検知してブロック境界を特定するブロック境界検出手段と、
前記ブロック境界において隣接する画素の信号の差分値を算出する手段と、
算出された差分値に応じて、前記平滑化信号と前記入力画像信号とを用いて画像出力信号を生成する画像調整手段と
を備えたことを特徴とする画像信号処理装置。
【請求項２】
前記ブロック境界において隣接する画素の信号の差分値に基づいて基準係数を算出する係数算出手段を更に備え、
前記画像調整手段は、
前記入力画像信号に対して、前記基準係数を用いて重み付けした値を乗算して入力画像調整信号を生成する第１乗算器と、
前記平滑化信号に基づいて生成された信号に対して、前記基準係数を乗算して平滑化調整信号を生成する第２乗算器と、
前記入力画像調整信号と前記平滑化調整信号とを加算して、前記画像出力信号を生成する第１加算手段と
を備えたことを特徴とする請求項１に記載の画像信号処理装置。
【請求項３】
前記係数算出手段は、算出した基準係数を用いて、前記ブロック境界からの画素数に応じて小さくなるフィルタ係数を算出し、
前記第１乗算器は、算出した前記フィルタ係数を前記基準係数として用いて入力画像調整信号を生成し、
前記第２乗算器は、算出した前記フィルタ係数を前記基準係数として用いて前記平滑化調整信号を生成することを特徴とする請求項２に記載の画像信号処理装置。
【請求項４】
前記入力画像信号から高域周波数成分を抽出する高域周波数抽出手段と、
高域周波数成分と画質制御値とを乗算して輪郭成分補正信号を生成する輪郭補正信号生成手段と、
前記平滑化信号と前記輪郭補正信号とを加算する第２加算手段と、
入力画像信号に対して、前記画質制御値を乗算する第３乗算器と、
前記第２加算手段の出力信号に対して、前記画質制御値を用いて重み付けした値を乗算する第４乗算器と、
前記第３乗算器の出力信号と前記第４乗算器の出力信号とを加算して、前記平滑化信号に基づいて生成された信号として、前記第２乗算器に供給する第３加算手段とを更に備えたことを特徴とする請求項２又は３に記載の画像信号処理装置。
【請求項５】
要素画素が２次元に並べられたブロックから構成された入力画像信号を取得し、この入力画像信号の高周波数成分を平滑化した平滑化信号を生成する平滑化段階と、
前記入力画像信号における輪郭成分を検出する輪郭検出段階と、
前記ブロック内の各要素画素のばらつき量を算出する分散量検出段階と、
前記ブロック内の中心画素について、前記輪郭成分と前記ばらつき量との比較に基づいて画質制御値を算出する制御信号生成段階と、
前記画質制御値を、画質制御値範囲の中間値を含む所定の基準範囲と比較することにより２値化した画素のパターン値を算出する段階と、
前記パターン値のマップとブロック規定パターンとの比較により、前記ブロック規定パターンの周期的配置を検知してブロック境界を特定するブロック境界検出段階と、
前記ブロック境界において隣接する画素の信号の差分値を算出する段階と、
算出された差分値に応じて、前記平滑化信号と前記入力画像信号とを用いて画像出力信号を生成する画像調整段階と
を備えたことを特徴とする画像信号処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像の信号圧縮符号化・復号化に起因して発生する画像のノイズを低減するための画像信号処理装置及び画像信号処理方法に関する。
【背景技術】
【０００２】
従来から、ノイズ除去回路として、フレームメモリを用いた巡回型のノイズリダクション装置が知られている（例えば、特許文献１参照。）。このノイズリダクション装置の構成の一例を図１５に示す。図１５に示すノイズリダクション装置は、ゲインＫを乗算する乗算器９１と、１つ前のフレームの画像を取得する遅延回路９２と、ゲイン（１−Ｋ）を乗算する乗算器９３と、乗算器９１，９３の出力信号を加算する加算器９４とから構成されている。そして、入力信号に対してゲインＫを乗算した値と、遅延回路９２において取得した前のフレームの信号に対してゲイン（１−Ｋ）を乗算した値とを加算した信号を出力信号として出力する。
【０００３】
また、従来のノイズ除去回路として、コアリングをシャープネス回路に応用したノイズリダクション装置がある（例えば、特許文献２参照。）。このコアリング処理部においては、周波数の広帯域に拡散した信号をホワイトノイズとして固定レベル（コアリングレベル）でカットする。このノイズリダクション装置の構成の一例を図１６に示す。図１６に示すノイズリダクション装置は、高周波フィルタ９５及び低周波フィルタ９６を備える。入力信号のうち高周波フィルタ９５の成分のうち、コアリングレベル以下の成分が除去され、残りの成分に対してゲイン調整される。そして、入力信号の低周波フィルタ９６を通過した信号と加算されて出力信号が生成される。
【０００４】
ところで、近年、従来のアナログ放送に代わるデジタル放送が開始されている。デジタル放送は、デジタル処理により画像を表示するため、従来のアナログ放送では発生しないノイズが発生することが知られている。例えば、ブロックサイズの境界が検知されるノイズ（以下、ブロックノイズという。）がある。このノイズは、画像をＭＰＥＧ２方式による圧縮を行なうため、圧縮率との関係で発生する平均輝度成分の量子化エラーによるものである。
【０００５】
このブロックノイズは、デジタル圧縮方式で記録媒体に記録された信号が、アナログ出力として変換されてテレビ受信機に供給される場合でも同様に発生される。言い換えれば、ブロックノイズは、デジタル又はアナログのインターフェースの形態に係わらず、圧縮復号された信号に内在する。更に、ブロックノイズには、以下の２つの特徴がある。
【０００６】
（ａ）固定パターンで発生する。
（ｂ）ブロックの境界における変化量が微小である。
ここで、（ａ）の特徴は、例えば８×８のブロックサイズで周期的に画面上にブロックノイズが現われることを意味する。（ｂ）の特徴は、そのブロックサイズの境界前後での変化量が微小であり、通常はあまり検知されないが、シーンチェンジの変化等に応じて圧縮率が変化する場合にその境界レベルが変動しノイズとして検知されることがある。特に、アナログ信号接続の場合は、アナログ信号入力のばらつき等の影響を受け、ブロック境界の変化量もそのばらつき影響を受け変動する。従って、従来のノイズ除去回路では、微小レベルが必ずしも一定ではないためにアナログ系ノイズとの識別が困難であり、更に固定パターンであるために信号との識別が困難である。
【特許文献１】特開平７−１６２８０７号公報（図３）
【特許文献２】特開平１０−２００７８７号公報（第１頁）
【発明の開示】
【発明が解決しようとする課題】
【０００７】
ところで、特許文献１に記載のノイズリダクション装置を用いた場合、静止画におけるランダムノイズの除去性能が高い。しかし、固定パターンを有するブロックノイズには効果がない。
【０００８】
更に、上述のコアリングの技術を用いる場合、抑圧度を向上させるためには、これらノイズの特性からコアリングレベルを大幅に上げる必要がある。
特許文献２に記載のノイズリダクション装置では、コアリングレベルは絶対値であるため、コアリングレベル以下の信号成分も失われる。更に、同一入力信号であっても、系によりコントラスト（ゲイン）が変化している場合、コアリングにより失われる信号成分も変化するという問題が生じる。
【０００９】
本発明は、上述の問題に鑑みてなされ、ブロックノイズと信号とを識別してブロックノイズに対する抑圧度を向上させながら、信号成分の損失を少なくすることができる画像信号処理装置及び画像信号処理方法を提供することにある。
【課題を解決するための手段】
【００１０】
上記問題点を解決するために、請求項１に記載の発明は、要素画素が２次元に並べられたブロックから構成された入力画像信号を取得し、この入力画像信号の高周波数成分を平滑化した平滑化信号を生成する平滑化手段と、前記入力画像信号における輪郭成分を検出する輪郭検出手段と、前記ブロック内の各要素画素のばらつき量を算出する分散量検出手段と、前記ブロック内の中心画素について、前記輪郭成分と前記ばらつき量との比較に基づいて画質制御値を算出する制御信号生成手段と、前記画質制御値を、画質制御値範囲の中間値を含む所定の基準範囲と比較することにより２値化した画素のパターン値を算出する手段と、前記パターン値のマップとブロック規定パターンとの比較により、前記ブロック規定パターンの周期的配置を検知してブロック境界を特定するブロック境界検出手段と、前記ブロック境界において隣接する画素の信号の差分値を算出する手段と、算出された差分値に応じて、前記平滑化信号と前記入力画像信号とを用いて画像出力信号を生成する画像調整手段とを備えたことを要旨とする。
【００１１】
請求項２に記載の発明は、請求項１に記載の画像信号処理装置において、前記ブロック境界において隣接する画素の信号の差分値に基づいて基準係数を算出する係数算出手段を更に備え、前記画像調整手段は、前記入力画像信号に対して、前記基準係数を用いて重み付けした値を乗算して入力画像調整信号を生成する第１乗算器と、前記平滑化信号に基づいて生成された信号に対して、前記基準係数を乗算して平滑化調整信号を生成する第２乗算器と、前記入力画像調整信号と前記平滑化調整信号とを加算して、前記画像出力信号を生成する第１加算手段とを備えたことを要旨とする。
【００１２】
請求項３に記載の発明は、請求項２に記載の画像信号処理装置において、前記係数算出手段は、算出した基準係数を用いて、前記ブロック境界からの画素数に応じて小さくなるフィルタ係数を算出し、前記第１乗算器は、算出した前記フィルタ係数を前記基準係数として用いて入力画像調整信号を生成し、前記第２乗算器は、算出した前記フィルタ係数を前記基準係数として用いて前記平滑化調整信号を生成することを要旨とする。
【００１３】
請求項４に記載の発明は、請求項２又は３に記載の画像信号処理装置において、前記入力画像信号から高域周波数成分を抽出する高域周波数抽出手段と、高域周波数成分と画質制御値とを乗算して輪郭成分補正信号を生成する輪郭補正信号生成手段と、前記平滑化信号と前記輪郭補正信号とを加算する第２加算手段と、入力画像信号に対して、前記画質制
御値を乗算する第３乗算器と、前記第２加算手段の出力信号に対して、前記画質制御値を用いて重み付けした値を乗算する第４乗算器と、前記第３乗算器の出力信号と前記第４乗算器の出力信号とを加算して、前記平滑化信号に基づいて生成された信号として、前記第２乗算器に供給する第３加算手段とを更に備えたことを要旨とする。
【００１４】
請求項５に記載の発明は、要素画素が２次元に並べられたブロックから構成された入力画像信号を取得し、この入力画像信号の高周波数成分を平滑化した平滑化信号を生成する平滑化段階と、前記入力画像信号における輪郭成分を検出する輪郭検出段階と、前記ブロック内の各要素画素のばらつき量を算出する分散量検出段階と、前記ブロック内の中心画素について、前記輪郭成分と前記ばらつき量との比較に基づいて画質制御値を算出する制御信号生成段階と、前記画質制御値を、画質制御値範囲の中間値を含む所定の基準範囲と比較することにより２値化した画素のパターン値を算出する段階と、前記パターン値のマップとブロック規定パターンとの比較により、前記ブロック規定パターンの周期的配置を検知してブロック境界を特定するブロック境界検出段階と、前記ブロック境界において隣接する画素の信号の差分値を算出する段階と、算出された差分値に応じて、前記平滑化信号と前記入力画像信号とを用いて画像出力信号を生成する画像調整段階とを備えたことを要旨とする。
【００１５】
（作用）
請求項１又は５に記載の発明によれば、入力画像信号における輪郭成分を検出し、ブロック内の要素画素のばらつき量を算出し、ブロック内の中心画素について、輪郭成分とばらつき量との比較に基づいて画質制御値を算出する。この画質制御値を、画質制御値範囲の中間値を含む所定の基準範囲と比較することにより２値化した画素のパターン値を算出し、パターン値のマップとブロック規定パターンとの比較により、ブロック境界を特定する。ブロック境界において隣接する画素の信号の差分値に応じて、平滑化信号と前記入力画像信号とを用いて画像出力信号を生成する。画質制御値は、ブロックに輪郭がない場合にはほぼ「０」になり、ブロックに輪郭がある場合にはほぼ「１」になるため、画質制御値が「０」と「１」との間にある場合にはノイズであることを示すパターン値を与える。更に、このパターン値のマップとブロック規定パターンとを比較してブロック境界にあるノイズ（ブロックノイズ）を特定することができる。特定したブロック境界で隣接する画素の信号の差分値に基づいて、平滑化信号の割合を自動的に変更することができる。従って、ブロックノイズと信号とを識別してブロックノイズに対応する抑圧度を向上させながら、入力画像信号に含まれている信号成分の損失を少なくすることができる。
【００１６】
請求項２に記載の発明によれば、ブロック境界において隣接する画素の信号の差分値に基づいて基準係数を算出する。画像調整手段は、入力画像信号に対して基準係数を用いて重み付けした値を乗算して入力画像調整信号を生成し、平滑化信号に基づいて生成された信号に対して基準係数を乗算して平滑化調整信号を生成し、これらを加算して画像出力信号を生成する。これにより、基準係数に基づいて平滑化信号と入力画像信号との調整を実現することができる。
【００１７】
請求項３に記載の発明によれば、係数算出手段は、算出した基準係数を用いて、前記ブロック境界からの画素数に応じて小さくなるフィルタ係数を算出する。第１乗算器及び第２乗算器は、算出したフィルタ係数を基準係数として用いて、平滑化調整信号と入力画像信号を生成する。このため、ブロックノイズの大きさに応じて、かつブロック境界から遠くなるに従ってフィルタ係数が小さくなり、これにより入力画像信号の割合が大きくなるように調整が行なわれる。従って、ブロックノイズを段階的に拡散するように調整することができるので、ブロックノイズの抑圧度の向上のためにブロック境界の画素の信号を調整して、この画素に隣接する画素との間において新たなノイズの発生を抑制することができる。
【００１８】
請求項４に記載の発明によれば、入力画像信号から高域周波数成分を抽出し、この高域周波数成分と画質制御値とを乗算して輪郭成分補正信号を生成し、これと平滑化信号とを加算する。入力画像信号に対して画質制御値を乗算し、第２加算手段の出力信号に対して画質制御値を用いて重み付けした値を乗算し、これらを加算して第２乗算器に供給する。平滑化信号は、高周波数成分を平滑化するため、ノイズが除去されているが、輪郭成分も失われてしまう。このため、抽出した高域周波数成分に基づく輪郭成分補正信号を平滑化信号に加算することにより、ノイズの状況に応じて加算される高域周波数成分を変更しながら平滑化調整信号を生成することができる。この平滑化調整信号を用いて画像出力信号が生成されるので、ノイズが除去され、高周波周波数成分がある平滑化信号を用いることができる。
【発明の効果】
【００１９】
本発明によれば、ブロックノイズと信号とを識別してブロックノイズに対する抑圧度を向上させながら、信号成分の損失を少なくすることができる。
【発明を実施するための最良の形態】
【００２０】
（第１実施形態）
以下、本発明を具体化したノイズリダクション装置の一実施形態を図１〜図１３に基づいて説明する。
【００２１】
図１に示すように、本実施形態のノイズリダクション装置１００は、輪郭検出フィルタ１０、ばらつき量検出手段２０、制御信号生成手段３０、平滑化フィルタ４０、ブロックノイズ調整手段５０及び第１制御部６０を備えている。本実施形態では、輪郭検出フィルタ１０が輪郭検出手段として機能し、ばらつき量検出手段２０が分散量検出手段として機能する。更に、平滑化フィルタ４０が平滑化手段として機能する。
【００２２】
このノイズリダクション装置１００には、コンポーネント映像信号が入力される。この映像信号のうち、ブロック構成部によって、ｎ行ｍ列（ｎ，ｍは正の整数）のブロックが、輪郭検出フィルタ１０、ばらつき量検出手段２０、平滑化フィルタ４０及び第１制御部６０に供給される。本実施形態では、ｎ＝ｍ＝３として、３×３のブロックで構成される要素画素の入力画像信号が各機能部に供給される。
【００２３】
〔ノイズリダクション装置１００の構成〕
図２に示すように、輪郭検出フィルタ１０は、隣接するピクセルの輝度の差分から、ブロック内のエッジ（輪郭）を検出する。本実施形態では、この輪郭検出フィルタ１０は、公知のプレヴィット（Ｐｒｅｗｉｔｔ）方式のフィルタを用いる。具体的には、輪郭検出フィルタ１０は、第１列、第２列及び第３列の各要素が「−１」、「０」、「１」の水平方向フィルタ１１と、第１行、第２行及び第３行の各要素が「−１」、「０」、「１」の垂直方向フィルタ１２とを備える。
【００２４】
このような輪郭検出フィルタ１０において、ブロックを構成する各画素の信号が水平方向フィルタ１１を通過すると、水平方向で隣接しているピクセルの輝度の差分を算出する。また、ブロックを構成する各画素の信号が垂直方向フィルタ１２を通過すると、垂直方向で隣接しているピクセルの輝度の差分を算出する。
【００２５】
更に、輪郭検出フィルタ１０は、輪郭検出信号ＹＥの算出部１３を備える。この算出部１３は、水平方向フィルタ１１及び垂直方向フィルタ１２を通過させて得た値と、図２の（１）式とを用いて輪郭検出信号ＹＥを生成する。輪郭検出フィルタ１０は、生成した輪郭検出信号ＹＥを制御信号生成手段３０に供給する。
【００２６】
ばらつき量検出手段２０は、ブロックを構成する画素のばらつき量（標準偏差σ）を算出する。ブロック内に輪郭がない場合には、通常はばらつき量は小さくなるため、ばらつき量が大きい場合にはノイズが発生していることになる。また、ブロック内に輪郭がある場合には、輪郭の量に応じてばらつき量が大きくなる。
【００２７】
本実施形態では、標準偏差σを算出する際には、平滑化フィルタ４０の統計平均値Ｓを用いる。ばらつき量検出手段２０は、図３の（２）式を用いて、ばらつき量（標準偏差σ）を算出する。ばらつき量検出手段２０は、取得した入力画像信号のブロックを構成する画素の信号からばらつき量（標準偏差σ）を算出し、算出した標準偏差σを制御信号生成手段３０に供給する。
【００２８】
制御信号生成手段３０は制御信号contを出力する。この制御信号contには、輪郭検出信号ＹＥを標準偏差σで正規化した画質制御値Ｇを含める。具体的には、図４に示すように、制御信号生成手段３０は、比較手段３１、上限設定手段３２及び制御値算出手段３３を備える。比較手段３１は、輪郭検出フィルタ１０から取得した輪郭検出信号ＹＥと、ばらつき量検出手段２０から取得した標準偏差σとを比較する。比較手段３１が、輪郭検出信号ＹＥが標準偏差σ以上であると判断した場合（「ＹＥＳ」の場合）、上限設定手段３２が、輪郭検出信号ＹＥの値を標準偏差σと同じ値に置換する。
【００２９】
制御値算出手段３３は、図４の（３）式に輪郭検出信号ＹＥと標準偏差σを代入して制御信号contの画質制御値Ｇを算出する。この画質制御値Ｇは、輪郭検出信号ＹＥが標準偏差σ以上のときには「１」になるので、「０」〜「１」の範囲の大きさになる。更に、この画質制御値Ｇは、信号のレベルが変動しても、輪郭検出信号ＹＥが標準偏差σより小さい場合には、「０」に近づき、輪郭検出信号ＹＥが標準偏差σより大きい場合には、「１」に近づく。この制御信号生成手段３０は、算出した画質制御値Ｇを、ブロックノイズ調整手段５０及び第１制御部６０に供給する。
【００３０】
一方、平滑化フィルタ４０は、入力画像信号から平坦領域を取り出すフィルタである。この平滑化フィルタ４０は、ブロック内の統計平均値Ｓを算出して、この統計平均値Ｓに応じた信号を出力信号とする。本実施形態では、９画素のメディアン値をばらつき量（標準偏差σ）を算出するための統計平均値Ｓとして用いる。具体的には、平滑化フィルタ４０は、図３に示すように、入力画像信号から取得した３×３のブロックを構成する画素（ｓ１〜ｓ９）のメディアン値を算出し、このメディアン値をブロックの統計平均値Ｓとする。メディアン値を特定することにより、比較的大きなエッジ成分（すなわち明らかな信号成分）に対しては可能な限り高次スペクトル成分を保持することができる。平滑化フィルタ４０は、算出したメディアン値の平滑化信号を、ばらつき量検出手段２０及び第１制御部６０に供給する。
【００３１】
図１に示すブロックノイズ調整手段５０は、パターン値決定手段５１、ブロック境界検出手段としてのパターンマッチング手段５２、前後画素差分算出手段５３及び係数算出手段５４を備えている。パターン値決定手段５１は、パターンマッチングに用いるパターン値を特定する。具体的には、パターン値決定手段５１は、図５に示すように、レベルシフタ５１１、パターン値比較手段５１２、第１パターン値特定手段５１３及び第２パターン値特定手段５１４を備えている。
【００３２】
レベルシフタ５１１は、制御信号contの画質制御値Ｇを、画像信号のレンジに応じた量子化数を乗算したシフト値Ｄを算出する。例えば、８ビットで量子化される画像信号の場合には、画質制御値Ｇを、２５５倍（「２」の８乗倍）にする。これにより、「０」〜「１」の範囲を取る画質制御値Ｇが、「０」〜「２５５」の範囲にシフトする。
【００３３】
パターン値比較手段５１２は、算出したシフト値Ｄが、所定の基準範囲Ｌ内であるか否かを判断する。本実施形態の基準範囲Ｌは、エッジ成分がないと認識される「０」の近傍の値（輪郭判断値ｄ１より大きい値）を除き、かつエッジ成分があると認識される「２５５」の近傍の値（輪郭判断値ｄ２より小さい値）を除く中間値の範囲（ｄ１〜ｄ２の範囲）である。
【００３４】
第１パターン値特定手段５１３は、基準範囲Ｌにある（パターン値比較手段５１２が「ＹＥＳ」と判断した）シフト値Ｄのパターン値を「１」と特定する。第２パターン値特定手段５１４は、基準範囲Ｌにない（パターン値比較手段５１２が「ＮＯ」と判断した）シフト値Ｄのパターン値を「０」と特定する。そして、パターン値決定手段５１は、特定したパターン値（「１」又は「０」）を画素の位置（座標）とともに一時記憶する。パターン値決定手段５１は、入力画像信号のうち最初の８×８のブロックを構成する画素のパターン値をすべて特定すると、ブロックを構成する画素のパターン値のマップをパターンマッチング手段５２に供給する。
【００３５】
パターンマッチング手段５２は、パターン値決定手段５１から供給されたパターン値のマップと、予め記憶されている配列パターン（ブロック規定パターン）とのマッチングを行なって、ブロック境界位置を検出する。本実施形態では、配列パターンとして、図７に示す８×８のブロックパターンを用いる。この配列パターンは、ブロック境界に隣接する外周の画素Ｐ１（ハッチングで示す画素）のパターン値が「１」で、ブロックの中央の画素Ｐ０（無地の画素）のパターン値が「０」になっている。なお、画素Ｐ１と中央の画素Ｐ０との間に位置する画素（ドット模様で示す画素）のパターン値は、「１」又は「０」の両方を取り得るので、この画素のパターン値はマッチングに用いない。
【００３６】
パターンマッチング手段５２は、パターン値決定手段５１から供給されたパターン値のマップに対して、配列パターンを１画素ずつずらして当てはめる。具体的には、パターンマッチング手段５２は、開始点のｉ行、ｊ列の座標（ｉ，ｊ）を、「８」の整数倍＋Ｚ（Ｚ＝０，１，・・・，７）という座標に順次変換する。
【００３７】
そして、パターンマッチング手段５２は、変換した各座標の画素のパターン値と、配列パターンのパターン値とが一致した場合には、マッチング数をカウントアップする。ここで、マッチング数とは、変換した各座標の画素のパターン値と、配列パターンのパターン値とが一致した画素の数である。そして、パターンマッチング手段５２は、配置パターンをそれぞれ当てはめたときのマッチング数の合計を算出する。この結果、マッチング数の合計が１箇所に集中する行と列が出てくる。例えば、図８のように、マッチング数の合計は、行Ｑａで、列Ｑｂに集中し、その前後の行［ｑａ±］と列［ｑｂ±］のマッチング数の合計より大きく上回る。この集中した一箇所（図８の点ＳＴ）の要素行Ｒａ及び要素列Ｒｂに対応する位置が境界位置に相当する。
【００３８】
この場合、入力画像信号の有効画面の開始点を（１，１）とし、マッチング数が集中した要素行の番号をｒａ、要素列の番号をｒｂとすると、境界位置は、（ｒａ＋８＊ｎ’，ｒｂ＋８＊ｍ’）で示される。ここで、ｎ’，ｍ’は「０」以上の整数である。図８の場合では、ｒａ＝３、ｒｂ＝５になっている。
【００３９】
そこで、パターンマッチング手段５２は、マッチング数の合計の最大値を特定し、この最大値に対応する要素列の番号ｒａと要素行の番号ｒｂを特定する。そして、パターンマッチング手段５２は、要素列の番号ｒａと要素行の番号ｒｂを用いて、境界位置（ｒａ＋８＊ｎ’，ｒｂ＋８＊ｍ’）を算出し、ブロック境界の位置を特定する。
【００４０】
次に、パターンマッチング手段５２は、特定したブロック境界の位置を示す検出信号を、検出ウィンドウの形で出力し、前後画素差分算出手段５３に供給する。ここで、出力される検出信号は、図９に示すように、境界位置に対して前後１画素分の幅を有するパルス波形の検出信号Ｐである。
【００４１】
前後画素差分算出手段５３は、ブロック境界において隣接する画素のエッジ変化量（差分値）を算出する。具体的には、前後画素差分算出手段５３は、供給された検出信号Ｐが出力されたときの画素（ブロック境界において前後１画素）の信号を取得し、隣接する画素の差分の絶対量を算出する。前後画素差分算出手段５３は、算出したエッジ変化量を係数算出手段５４に供給する。
【００４２】
係数算出手段５４は、算出したエッジ変化量と各閾値とを比較する。本実施形態では、水平方向の閾値として３つの閾値ｈ１，ｈ２，ｈ３を用いる。更に、本実施形態では、垂直方向の閾値として３つの閾値ｖ１，ｖ２，ｖ３を用いる。
【００４３】
更に、係数算出手段５４は、エッジ変化量の大きさに対応する基準係数ｄｈ（ｄｖ）のテーブルを備えている。本実施形態では、エッジ変化量が大きければ基準係数ｄｈ（ｄｖ）が大きくなるように設定されている。具体的には、図１０に示すように、エッジ変化量が閾値ｈ１（ｖ１）より小さい場合には、「０」の基準係数ｄｈ（ｄｖ）が対応付けられている（領域（１））。エッジ変化量が閾値ｈ１（ｖ１）以上で閾値ｈ２（ｖ２）より小さい場合には、「１／４」の基準係数ｄｈ（ｄｖ）が対応付けられている（領域（２））。エッジ変化量が閾値ｈ２（ｖ２）以上で閾値ｈ３（ｖ３）より小さい場合には、「１／２」の基準係数ｄｈ（ｄｖ）が対応付けられている（領域（３））。エッジ変化量が閾値ｈ３（ｖ３）以上の場合には、「１」の基準係数ｄｈ（ｄｖ）が対応付けられている（領域（４））。
【００４４】
係数算出手段５４は、エッジ変化量と各閾値とを比較した結果、算出されたエッジ変化量に対応する基準係数ｄｈ（ｄｖ）をテーブルから特定する。
更に、係数算出手段５４は、基準係数を用いてフィルタ係数ｋを算出する。ここでは、ブロック境界にある第１画素に隣接する第２画素、これに隣接する第３画素というように、水平方向及び垂直方向へ広がる隣接画素に対してフィルタ係数ｋが小さくなるように設定する。本実施形態では、係数算出手段５４は、ブロック境界からの画素の数を「ｎ”」とした場合、「１／２」の（ｎ”−１）乗のフィルタ係数ｋを算出する。具体的には、図１１に示すように、係数算出手段５４は、境界に位置する第１画素（「±１」の画素）に対しては、基準係数ｄｈ（ｄｖ）と同じ値をフィルタ係数ｋとして算出する。第１画素に隣接し第１画素よりも境界から離れた第２画素（図１１の「±２」の画素）に対しては、基準係数ｄｈ（ｄｖ）の「１／２」の値をフィルタ係数ｋとして算出する。更に、第２画素に隣接し第２画素よりも境界から離れた第３画素（図１１の「±３」の画素）に対しては、基準係数ｄｈ（ｄｖ）の「１／４」の値をフィルタ係数ｋとして算出する。更に、本実施形態では、それ以外の画素に対しては、「０」をフィルタ係数ｋとして算出する。そして、係数算出手段５４は、算出した各画素のフィルタ係数ｋをブロックノイズ調整手段５０の出力として第１制御部６０に供給する。
【００４５】
第１制御部６０は、画像調整手段として機能し、第１乗算器６１、第２乗算器６２及び加算器６３を備えている。第１乗算器６１及び第２乗算器６２は、フィルタ係数ｋを取得する。第１乗算器６１は、入力画像信号に対して「１−ｋ」（フィルタ係数ｋを用いて重み付けした値）を乗算した入力画像調整信号を出力する。第２乗算器６２は、平滑化フィルタ４０からの出力信号に対して、フィルタ係数ｋを乗算した平滑化調整信号を出力する。加算器６３は、第１乗算器６１からの入力画像調整信号と第２乗算器６２から平滑化調整信号とが供給され、これらを加算した信号を、第１制御部６０の出力信号として出力す
る。この第１制御部６０の出力信号が、ノイズリダクション装置１００の出力信号（画像出力信号）になる。
【００４６】
〔ノイズリダクション装置１００の動作〕
次に、このノイズリダクション装置１００の処理動作について説明する。
映像信号がノイズリダクション装置１００に入力された場合、輪郭検出フィルタ１０、ばらつき量検出手段２０、平滑化フィルタ４０及び第１制御部６０には、映像信号のうち３行３列のブロックが順次、供給される。
【００４７】
輪郭検出フィルタ１０は、輪郭検出信号ＹＥを算出して制御信号生成手段３０に供給する。ばらつき量検出手段２０は、標準偏差σを算出して制御信号生成手段３０に供給する。制御信号生成手段３０は、輪郭検出信号ＹＥと標準偏差σとを用いて画質制御値Ｇを算出する。そして、この画質制御値Ｇを含めた制御信号contを、ブロックノイズ調整手段５０に供給する。
【００４８】
平滑化フィルタ４０は、ブロックを構成する要素画素のメディアン値を算出し、ばらつき量検出手段２０と第１制御部６０の第２乗算器６２に供給する。この場合、平滑化フィルタ４０の出力信号は、モスキートノイズやブロックノイズが抑圧されているが、輪郭成分も失われている。
【００４９】
ブロックノイズ調整手段５０のパターン値決定手段５１は、制御信号contの画質制御値Ｇから各画素のパターン値を特定して、パターンマッチング手段５２に供給する。パターンマッチング手段５２は、マッチング数を合計してブロック境界の位置を示す検出信号を出力する。前後画素差分算出手段５３は、検出信号に基づいて境界位置の前後の画素の差分を算出する。係数算出手段５４は、前後の画素の差分から境界に位置する各画素の基準係数ｄｈ（ｄｖ）を算出し、各画素のフィルタ係数ｋを算出して第１制御部６０に供給する。
【００５０】
第１制御部６０の第１乗算器６１は、入力画像信号に対して「１−ｋ」を乗算して、加算器６３に供給する。第１制御部６０の第２乗算器６２は、平滑化フィルタ４０からの出力信号に対して「ｋ」を乗算した信号を算出して、加算器６３に供給する。加算器６３は、第１乗算器６１からの出力と第２乗算器６２からの出力とを加算して、ノイズリダクション装置１００の出力信号として出力する。
【００５１】
発生したブロックノイズにより、ブロック境界の前後でエッジ変化量が大きい場合には、ブロックノイズ調整手段５０が出力するフィルタ係数ｋが大きくなる。この場合、ノイズリダクション装置１００は、平滑化フィルタ４０からの出力信号の割合が大きく、入力画像信号の出力信号の割合が小さい合成信号を生成して出力する。
【００５２】
一方、ブロック境界の前後でエッジ変化量が小さい場合には、ブロックノイズ調整手段５０が出力するフィルタ係数ｋが小さくなる。この場合、ノイズリダクション装置１００は、平滑化フィルタ４０からの出力信号の割合が小さく、入力画像信号の出力信号の割合が大きい合成信号を生成して出力する。
【００５３】
図１２には、ノイズリダクション装置１００に入力される映像信号をそのまま用いて表示した参考画像を示す。図１３には、この映像信号が入力されたノイズリダクション装置１００から出力される画像出力信号を用いて表示した参考画像を示す。図１２の画像では、図１３の画像に比べて、ブロックノイズが低減していることがわかる。
【００５４】
本実施形態によれば、以下のような効果を得ることができる。
（１）  本実施形態では、制御信号生成手段３０は、輪郭検出信号ＹＥと標準偏差σとの比に応じた画質制御値Ｇの制御信号contを生成し、ブロックノイズ調整手段５０に供給する。ブロックノイズの場合は、図６に示すテーブル６００の中で、水平方向又は垂直方向に連続する数画素とその左右上下画素との輝度変化パターンが（σ小、ＹＥ大）の領域に含まれる。σ＝ＹＥの条件で、輪郭検出信号ＹＥ及び標準偏差σを徐々に減少させるように入力画像信号のゲインを変化させていくと、信号コントラストは「０」に近づくが、信号画素間の相対的関係は変化していないため、「Ｇ＝１」が保たれる。これは信号の画素間の相対的レベル比が変化しなければ、絶対値のコントラストには依存せずに目的とする輪郭成分（信号の差分）の大きさを抽出し、係数化できることを意味する。すなわち、ノイズがなくブロック内に輪郭がある場合には、「Ｇ≒１」になり、ノイズがなくブロック内に輪郭がない場合には、「Ｇ≒０」になる。そして、ノイズがある場合には、画質制御値Ｇは「０」と「１」の中間値になる。
【００５５】
ブロックノイズ調整手段５０のパターン値決定手段５１は、制御信号contの画質制御値Ｇが、「０」と「１」の中間値である基準範囲Ｌにある場合にはパターン値を「１」とする。パターンマッチング手段５２は、配列パターンとのマッチングを行なって、ブロック境界の位置を検出する。前後画素差分算出手段５３は、ブロック境界において隣接する画素のエッジ変化量を算出する。係数算出手段５４は、算出したエッジ変化量に対応する基準係数ｄｈ（ｄｖ）を算出しフィルタ係数ｋを算出する。第１制御部６０は、入力画像信号に対して「１−ｋ」を乗算した入力画像調整信号と、平滑化フィルタ４０からの出力信号に対して、「ｋ」を乗算した平滑化調整信号とを加算した信号を、ノイズリダクション装置１００の出力信号として出力する。このため、ブロック境界を特定し、ブロック境界において隣接する画素のエッジ変化量からブロックノイズの大きさを特定し、この大きさに応じて入力画像信号と平滑化調整信号との割合が自動的に変更される。従って、ノイズが大きいときには平滑化調整信号の割合が自動的に多くなるので、ブロックノイズの抑圧度を向上することができる。また、ノイズが小さいときには、入力画像信号の割合が自動的に多くなるので、入力画像信号に含まれている信号成分の損失を少なくすることができる。
【００５６】
（２）  本実施形態では、ブロックノイズ調整手段５０のパターン値決定手段５１は、制御信号contの画質制御値Ｇを、画像信号のレンジに応じた量子化数を乗算したシフト値Ｄを算出した後、算出したシフト値Ｄが基準範囲Ｌ内であるか否かを判断する。このため、微小レベルの画質制御値Ｇが、「０」〜「１」の範囲から「０」〜「２５５」の範囲に変換されることになる。画質制御値Ｇが、「０」〜「１」の範囲では、低輝度（黒）であるため、信号の判別が難しく、境界位置の検出が困難であった。「０」〜「２５５」の範囲に変換されることにより、ブロック境界を検出し易くできる。
【００５７】
（３）  本実施形態では、ブロックノイズ調整手段５０の係数算出手段５４は、ブロック境界からの画素の数を「ｎ”」とした場合、「１／２」の（ｎ”−１）乗のフィルタ係数ｋを算出する。このため、ブロック境界に位置する第１画素以外の画素についても調整を行なうので、ブロックノイズが大きい場合には、そのノイズをブロックの内側に拡散していくように信号を調整することができる。従って、ブロックノイズの抑圧度を向上させるためにブロック境界の画素の信号を調整したときに、調整した画素とこれに隣接する画素との間において新たなノイズの発生を抑制することができる。
【００５８】
（４）  本実施形態では、ばらつき量検出手段２０は、統計平均値Ｓを算出し、これに応じた信号を出力信号とする。本実施形態では、第１制御部６０は、平滑化フィルタ４０からの出力信号に対して、「ｋ」を乗算した平滑化調整信号と、入力画像調整信号とを加算した信号を出力信号として出力する。このため、平滑化フィルタ４０とばらつき量検出手段２０とが、同じ統計平均値Ｓを用いた演算を行なうことにより、整合性を良好にする
ことができる。
【００５９】
（５）  本実施形態では、制御信号生成手段３０において、比較手段３１が、輪郭検出信号ＹＥが標準偏差σ以上であると判断した場合（「ＹＥＳ」の場合）、上限設定手段３２が、輪郭検出信号ＹＥの値を標準偏差σと同じ値に置換する。制御信号生成手段３０の制御値算出手段３３は、図４の（３）式に輪郭検出信号ＹＥと標準偏差σを代入して制御信号contの画質制御値Ｇを算出する。このため、輪郭検出信号ＹＥが標準偏差σより大きな範囲を取る場合であっても、制御信号contの画質制御値Ｇを「０」〜「１」の間に設定することができるので、ノイズリダクション装置１００からの出力信号を、入力画像信号に応じた大きさで調整することができる。
【００６０】
（６）  本実施形態では、輪郭検出フィルタ１０は、センタ行又はセンタ列のみ変化している場合には反応しない公知のプレヴィット方式のフィルタを用いる。ブロックノイズは、水平方向又は垂直方向に延びるエッジであるため、プレヴィット方式を用いることにより、ブロックノイズの特徴と重なる特定のパターンを輪郭として検出しにくい。従って、隣接する画素の輝度の差分が高いノイズを輪郭として検出し難いので、ブロックノイズと輪郭とをより確実に識別することができる。
【００６１】
（第２実施形態）
次に、本発明を具体化した第２実施形態を図１４に基づいて説明する。本実施形態は、上記実施形態と、フィルタ係数ｋを乗算する調整信号の生成の方法が異なる。また、本実施形態において、上記実施形態と同様な部分については、同一の符号を用いて、その詳細な説明は省略する。
【００６２】
図１４に示すように、本実施形態のノイズリダクション装置２００は、高域周波数抽出手段としての高域周波数抽出フィルタ７０、輪郭補正信号生成手段としての乗算器７１、加算器７２及び第２制御部８０を更に備えている。高域周波数抽出フィルタ７０には、入力画像信号に基づくｎ行ｍ列のブロックが供給される。本実施形態では、上記実施形態と同様に、３行３列のブロックが供給される。乗算器７１及び第２制御部８０には、制御信号生成手段３０が生成した制御信号contが供給され、乗算器７１及び第２制御部８０は制御信号contの画質制御値Ｇを取得する。
【００６３】
高域周波数抽出フィルタ７０は、高域周波数に相当する成分のうち、輪郭成分の中核的な周波数のみを抽出する。高域周波数抽出フィルタ７０は、入力画像信号から所定の範囲の周波数成分を乗算器７１に供給する。乗算器７１は、高域周波数抽出フィルタ７０からの出力信号に対して、制御信号contの画質制御値Ｇを乗算した輪郭成分補正信号を出力する。加算器７２は、第２加算手段として機能し、平滑化フィルタ４０からの平滑化信号と、乗算器７１からの輪郭成分補正信号とを加算する。この加算器７２の出力信号が第２制御部８０に供給される。従って、この加算器７２からの出力信号は、画質制御値Ｇに応じて輪郭成分の信号強度を制御し、これを平滑化信号に重畳した信号になる。
【００６４】
第２制御部８０は、第３乗算器８１、第４乗算器８２及び、第３加算手段としての加算器８３を備えている。第３乗算器８１及び第４乗算器８２は、制御信号生成手段３０が生成した制御信号contが供給されて、画質制御値Ｇを取得する。第３乗算器８１は、入力画像信号に対して画質制御値Ｇを乗算した信号を出力する。第４乗算器８２は、加算器７２からの出力信号に対して、「１−Ｇ」（画質制御値を用いて重み付けした値）を乗算した信号を出力する。加算器８３には、第３乗算器８１からの信号と第４乗算器８２から信号とが供給され、これら信号を加算した信号を、第２制御部８０の出力信号として出力する。この第２制御部８０の出力信号が、第１制御部６０の第２乗算器６２に供給される。
【００６５】
本実施形態においては、上記第１実施形態における（１）〜（６）と同様な効果に加えて、以下の効果を得ることができる。
（７）  本実施形態では、制御信号生成手段３０は、生成した制御信号contの画質制御値Ｇは、乗算器７１及び第２制御部８０に供給する。第２制御部８０は、ノイズの状況に応じて変化する画質制御値Ｇに基づいて、第３乗算器８１からの信号と第４乗算器８２から信号とを加算する割合を変化させて、第１制御部６０に供給する。第１制御部６０は、フィルタ係数ｋに基づいて、第２制御部８０からの出力信号に基づく平滑化調整信号と入力画像信号に基づく入力画像調整信号とを加算する割合を変化させた信号を、ノイズリダクション装置２００の出力信号（画像出力信号）として出力する。このため、ブロック境界におけるノイズは、ブロック境界に位置する画素のエッジ変化量だけでなく、ノイズの状況に応じて変化する画質制御値Ｇからの調整が行なわれる。従って、より的確な信号を用いて画像を表示することができる。
【００６６】
（８）  本実施形態では、加算器７２は、高域周波数抽出フィルタ７０からの出力信号に対して画質制御値Ｇを乗算し、これを平滑化フィルタ４０に加算して平滑化調整信号を生成する。平滑化フィルタ４０は、ノイズはないが輪郭も損失している出力信号を出力する。このため、高域周波数成分に画質制御値Ｇを乗算した値を平滑化信号に加算することにより、ノイズの状況に応じて高域周波数成分を加えて生成される信号を第２制御部８０に供給し、この信号に応じた出力を行なうことができる。
【００６７】
また、上記実施形態は、以下のように変更してもよい。
○  上記実施形態においては、３つの閾値とエッジ変化量とを比較して基準係数ｄｈ（ｄｖ）を決定した。これに限らず、基準係数の算出方法は、他の方法であってもよい。例えば、エッジ変化量を変数とする基準係数の算出式に、算出したエッジ変化量を代入して基準係数を算出してもよい。
【００６８】
○  上記実施形態においては、係数算出手段５４は、ブロック境界からの画素の数を「ｎ”」とした場合、「１／２」の（ｎ”−１）乗のフィルタ係数ｋを算出する。これに限らず、ブロック境界からの画素数が大きくなるに従って、フィルタ係数ｋが小さくなるように設定されれば、他の算出方法を用いてもよい。更に、上記実施形態では、ブロック境界から画素が４つ以上離散した場合にはフィルタ係数ｋを「０」とした。フィルタ係数ｋを「０」とする画素数もこれに限られない。例えば、基準係数の大きさによってフィルタ係数ｋを「０」とする画素数を変更してもよい。
【００６９】
○  上記実施形態においては、平滑化フィルタ４０を用いたが、モスキートノイズやブロックノイズのない信号の平坦領域だけを取り出すことができれば、これに限られない。例えば、カットオフの低い単純なローパスフィルタを採用することも可能である。また、本実施形態においては、９画素のメディアン値を、標準偏差を算出するための統計平均値Ｓとして用いたが、簡易手法によって統計平均値Ｓを求めてもよい。例えば、図３における３×３のブロックを構成する画素のうち、５つの画素（図３の入力画像信号の画素ｓ１，ｓ３，ｓ５，ｓ７，ｓ９）のメディアン値と、残りの４つの画素（ｓ２，ｓ４，ｓ６，ｓ８）の平均値とを加算した値の半分の値を統計平均値Ｓとして用いてもよい。
【００７０】
○  上記実施形態において、ばらつき量検出手段２０は、平滑化フィルタ４０の出力を用いて算出した。これに限らず、ばらつき量検出手段２０は、ブロックの全画素の平均値として、単純な９画素の加算平均を用いてもよい。
【００７１】
○  上記実施形態においては、制御信号生成手段３０は、ばらつき量検出手段２０が算出した標準偏差σを用いて、制御信号contの画質制御値Ｇを算出した。これに限らず、制御信号生成手段３０は、標準偏差σに限られず、ばらつき量を示す値であれば、他のもの
であってもよい。例えば、標準偏差σの２乗の値である分散を用いて、制御信号contの値Ｇを算出してもよい。
【図面の簡単な説明】
【００７２】
【図１】第１実施形態のノイズリダクション装置の構成を説明するブロック図。
【図２】輪郭検出フィルタの構成を説明するための処理ブロック図。
【図３】ばらつき量検出手段の構成を説明する処理ブロック図。
【図４】制御信号生成手段の構成を説明する処理ブロック図。
【図５】パターン値決定手段の構成を説明する処理ブロック図。
【図６】輪郭検出信号及び標準偏差の大きさに応じた信号を説明するための説明図。
【図７】ブロック内のパターン値の分布を説明する説明図。
【図８】ブロック境界を特定する方法を説明するための説明図。
【図９】ブロック境界で隣接する画素のエッジ変化量を説明するための説明図。
【図１０】エッジ変化量と基準係数との関係を説明する説明図。
【図１１】画素とフィルタ係数との関係を説明するための説明図。
【図１２】入力画像信号に基づいて表示を行なった画像を説明する説明図。
【図１３】画像出力信号に基づいて表示を行なった画像を説明する説明図。
【図１４】第２実施形態のノイズリダクション装置の構成を説明するブロック図。
【図１５】従来技術におけるノイズリダクション装置の構成を説明する説明図。
【図１６】従来技術におけるノイズリダクション装置の構成を説明する説明図。
【符号の説明】
【００７３】
ｄｈ，ｄｖ…基準係数、Ｇ…画質制御値、ｋ…フィルタ係数、Ｌ…基準範囲、Ｐ０，Ｐ１…画素、１０…輪郭検出手段としての輪郭検出フィルタ、２０…分散量検出手段としてのばらつき量検出手段、３０…制御信号生成手段、４０…平滑化手段としての平滑化フィルタ、５０…ブロックノイズ調整手段、５１…パターン決定手段、５２…ブロック境界検出手段としてのパターンマッチング手段、５３…前後画素差分算出手段、５４…係数算出手段、６０…画像調整手段としての第１制御部、６１…第１乗算器、６２…第２乗算器、６３…第１加算手段としての加算器、７０…高域周波数抽出手段としての高域周波数抽出フィルタ、７１…輪郭補正信号生成手段としての乗算器、７２…第２加算手段としての加算器、８０…第２制御部、８１…第３乗算器、８２…第４乗算器、８３…第３加算手段としての加算器、１００，２００…画像信号処理装置としてのノイズリダクション装置。
【図１】
2009118081.tif 000003
【図２】
2009118081.tif 000004
【図３】
2009118081.tif 000005
【図４】
2009118081.tif 000006
【図５】
2009118081.tif 000007
【図６】
2009118081.tif 000008
【図７】
2009118081.tif 000009
【図８】
2009118081.tif 000010
【図９】
2009118081.tif 000011
【図１０】
2009118081.tif 000012
【図１１】
2009118081.tif 000013
【図１４】
2009118081.tif 000014
【図１５】
2009118081.tif 000015
【図１６】
2009118081.tif 000016
【図１２】
2009118081.tif 000017
【図１３】
2009118081.tif 000018