(57)【要約】【課題】効率的に量子化を行う。【解決手段】ＭＰＥＧ４の量子化装置１０は、量子化スケール値（ｑ＿ｓｃａｌｅ）を発生する量子化スケール発生部１１と、１／２丸め演算を行う第１の演算部１２と、２／３丸め演算を行う第２の演算部１３と、丸め演算式を選択するセレクタ１４と、コントローラ１５とを備えている。量子化装置１０は、整数値である元値（ｃｏｅｆ）を０又は自然数である所定の量子化スケール値（ｑ＿ｓｃａｌｅ）に基づき量子化し、出力値（Ｃ）を出力する。この量子化装置では、３／４丸め演算を行う場合に、下式の演算を用いる。Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　｛｜Ｘ｜　−　（ｑ＿ｓｃａｌｅ　＋　１）　／　２｝　／　（２　＊　ｑ＿ｓｃａｌｅ）【選択図】　図１
2004048191.tif 000002
【特許請求の範囲】
【請求項１】
整数値である元値（Ｘ）を０又は自然数である所定の量子化スケール値（ｑ＿ｓｃａｌｅ）に基づき量子化し、出力値（Ｙ）を出力する量子化装置において、
下記式（１）に基づき出力値（Ｙ）を算出すること
を特徴とする量子化装置。
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　｛｜Ｘ｜　−　（ｑ＿ｓｃａｌｅ　＋　１）　／　２｝　／　（２　＊　ｑ＿ｓｃａｌｅ）　…（１）
なお、上記式（１）では、Ａ，Ｂが任意の整数値であるとした場合、“Ａ　＋　Ｂ”はＡとＢとを加算した値、“Ａ　−　Ｂ”はＡからＢを減算した値、”Ａ　＊　Ｂ”はＡとＢとを乗算した値、“Ａ　／　Ｂ”はＡをＢで除算して小数点以下を切り捨てた値、｜Ａ｜はＡの絶対値、ｓｉｇｎ（Ａ）は絶対値が１で正負符号がＡと同一の値、を意味している。
【請求項２】
下記式（２）に基づき出力値（Ｙ）を算出する第１の丸め演算部と、
上記式（１）に基づき出力値（Ｙ）を算出する第２の丸め演算部と、
第１の丸め演算部からの出力値（Ｙ）と第２の丸め演算部からの出力値（Ｙ）のうち一方の出力値を選択して出力する選択部とを備えること
を特徴とする請求項１記載の量子化装置。
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　（｜Ｘ｜）／　（２　＊　ｑ＿ｓｃａｌｅ）　…（２）
なお、上記式（２）では、Ａ，Ｂが任意の整数値であるとした場合、”Ａ　＊　Ｂ”はＡとＢとを乗算した値、“Ａ　／　Ｂ”はＡをＢで除算して小数点以下を切り捨てた値、｜Ａ｜はＡの絶対値、ｓｉｇｎ（Ａ）は絶対値が１で正負符号がＡと同一の値、を意味している。
【請求項３】
上記選択部は、元値（Ｘ）がイントラマクロブロックに含まれている画像データの場合には上記第１の丸め演算部の出力値（Ｙ）を選択し、元値（Ｘ）がインターマクロブロックに含まれている画像データの場合には上記第２の丸め演算部の出力値（Ｙ）を選択すること
を特徴とする請求項２記載の量子化装置。
【請求項４】
上記選択部は、当該量子化装置が符号化装置に用いられている際に、その符号化装置によるデータの圧縮率が所定の値よりも小さい場合には上記第１の丸め演算部の出力値（Ｙ）を選択し、圧縮率が所定の値以上の場合には上記第２の丸め演算部の出力値（Ｙ）を選択すること
を特徴とする請求項２記載の量子化装置。
【請求項５】
上記選択部は、出力値（Ｙ）が所定の閾値より大きい場合には上記第１の丸め演算部の出力値（Ｙ）を選択し、出力値（Ｙ）が所定の閾値以下の場合には上記第２の丸め演算部の出力値（Ｙ）を選択すること
を特徴とする請求項２記載の量子化装置。
【請求項６】
元値（Ｘ）は、離散コサイン変換された後の係数であり、
上記選択部は、元値（Ｘ）の周波数に応じて、第１の丸め演算部からの出力値（Ｙ）又は第２の丸め演算部からの出力値（Ｙ）のうちいずれか一つの出力値を選択すること
を特徴とする請求項２記載の量子化装置。
【請求項７】
下記式（２）に基づき出力値（Ｙ）を算出する第１の丸め演算部と、
上記式（１）に基づき出力値（Ｙ）を算出する第２の丸め演算部と、
下記式（３）に基づき出力値（Ｙ）を算出する第３の丸め演算部と、
第１の丸め演算部からの出力値（Ｙ）、第２の丸め演算部からの出力値（Ｙ）又は第３の丸め演算部からの出力値（Ｙ）のうちいずれか一つの出力値を選択して出力する選択部とを備えること
を特徴とする請求項１記載の量子化装置。
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　（｜Ｘ｜）／　（２　＊　ｑ＿ｓｃａｌｅ）　…（２）
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　（｜Ｘ｜　−　ｑ＿ｓｃａｌｅ）／　（２　＊　ｑ＿ｓｃａｌｅ）　…（３）
なお、上記式（２），（３）では、Ａ，Ｂが任意の整数値であるとした場合、“Ａ−　Ｂ”はＡからＢを減算した値、”Ａ　＊　Ｂ”はＡとＢとを乗算した値、“Ａ　／　Ｂ”はＡをＢで除算して小数点以下を切り捨てた値、｜Ａ｜はＡの絶対値、ｓｉｇｎ（Ａ）は絶対値が１で正負符号がＡと同一の値、を意味している。
【請求項８】
元値（Ｘ）は、離散コサイン変換された後の係数であり、
上記選択部は、元値（Ｘ）の周波数に応じて、第１の丸め演算部からの出力値（Ｙ）、第２の丸め演算部からの出力値（Ｙ）又は第３の丸め演算部からの出力値（Ｙ）のうちいずれか一つの出力値を選択すること
を特徴とする請求項７記載の量子化装置。
【請求項９】
整数値である元値（Ｘ）を０又は自然数である所定の量子化スケール値（ｑ＿ｓｃａｌｅ）に基づき量子化し、出力値（Ｙ）を出力する量子化方法において、
下記式（１）に基づき出力値（Ｙ）を算出すること
を特徴とする量子化方法。
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　｛｜Ｘ｜　−　（ｑ＿ｓｃａｌｅ　＋　１）　／　２｝　／　（２　＊　ｑ＿ｓｃａｌｅ）　…（１）
なお、上記式（１）では、Ａ，Ｂが任意の整数値であるとした場合、“Ａ　＋　Ｂ”はＡとＢとを加算した値、“Ａ　−　Ｂ”はＡからＢを減算した値、”Ａ　＊　Ｂ”はＡとＢとを乗算した値、“Ａ　／　Ｂ”はＡをＢで除算して小数点以下を切り捨てた値、｜Ａ｜はＡの絶対値、ｓｉｇｎ（Ａ）は絶対値が１で正負符号がＡと同一の値、を意味している。
【請求項１０】
下記式（２）に基づき出力値（Ｙ）を算出する第１の丸め演算、及び、上記式（１）に基づき出力値（Ｙ）を算出する第２の丸め演算のうち一方の出力値を選択して出力すること
を特徴とする請求項９記載の量子化方法。
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　（｜Ｘ｜）／　（２　＊　ｑ＿ｓｃａｌｅ）　…（２）
なお、上記式（２）では、Ａ，Ｂが任意の整数値であるとした場合、”Ａ　＊　Ｂ”はＡとＢとを乗算した値、“Ａ　／　Ｂ”はＡをＢで除算して小数点以下を切り捨てた値、｜Ａ｜はＡの絶対値、ｓｉｇｎ（Ａ）は絶対値が１で正負符号がＡと同一の値、を意味している。
【請求項１１】
下記式（２）に基づき出力値（Ｙ）を算出する第１の丸め演算、上記式（１）に基づき出力値（Ｙ）を算出する第２の丸め演算、及び、下記式（３）に基づき出力値（Ｙ）を算出する第３の丸め演算のうちいずれか一つの出力値を選択して出力すること
を特徴とする請求項９記載の量子化方法。
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　（｜Ｘ｜）／　（２　＊　ｑ＿ｓｃａｌｅ）　…（２）
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　（｜Ｘ｜　−　ｑ＿ｓｃａｌｅ）／　（２　＊　ｑ＿ｓｃａｌｅ）　…（３）
なお、上記式（２），（３）では、Ａ，Ｂが任意の整数値であるとした場合、“Ａ−　Ｂ”はＡからＢを減算した値、”Ａ　＊　Ｂ”はＡとＢとを乗算した値、“Ａ　／　Ｂ”はＡをＢで除算して小数点以下を切り捨てた値、｜Ａ｜はＡの絶対値、ｓｉｇｎ（Ａ）は絶対値が１で正負符号がＡと同一の値、を意味している。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、例えば画像や音声等の符号化データの符号量の削減をする量子化装置及び方法に関するものである。
【０００２】
【従来の技術】
ＭＰＥＧ２（ＩＳＯ／ＩＥＣ　１３８１８−２）及びＭＰＥＧ４（ＩＳＯ／ＩＥＣ　１４４９６−２）の画像符号化規格に対応した画像符号化装置を図２３に示し、画像復号装置のブロック構成図を図２４に示す。
【０００３】
画像符号化装置１０１は、図２３に示すように、第１の加算回路１０２と、離散コサイン変換（ＤＣＴ）回路１０３と、量子化回路１０４と、逆量子化回路１０５と、逆離散コサイン変換（ＩＤＣＴ）回路１０６と、第２の加算回路１０７と、フレームメモリ１０８と、動き補償回路１０９とを備えて構成される。
【０００４】
第１の加算回路１０２には、符号化対象となる空間領域の画像データ（ｐｉｘｅｌ）がフレーム単位で入力される。また、第１の加算回路１０２には、符号化対象となる画像データに対して、フレーム間の相関を利用した符号化処理を行う場合には、つまり、Ｐピクチャ又はＢピクチャの符号化をする場合には、動き補償回路１０９から予測画像データも入力される。第１の加算回路１０２は、インターマクロブロックに対して処理を行う場合には、入力された画像データから予測画像データを減算する。第１の加算回路１０２は、イントラマクロブロックに対して処理を行う場合には、入力された画像データをそのまま出力する。
【０００５】
ＤＣＴ回路１０３は、第１の加算回路１０２から出力された画像データに対して離散コサイン変換を行い、周波数領域の画像データであるＤＣＴ係数データ（ｃｏｅｆ）を生成する。ＤＣＴ回路１０３は、生成したＤＣＴ係数データ（ｃｏｅｆ）を量子化回路１０４に出力する。
【０００６】
量子化回路１０４は、入力されたＤＣＴ係数データ（ｃｏｅｆ）に対して、所定の量子化スケール値（ｑ＿ｓｃａｌｅ）を用いて、量子化処理を行い、量子化データ（Ｃ）を出力する。
【０００７】
逆量子化回路１０５には、量子化回路１０４から出力された量子化データ（Ｃ）のうち、参照画像データとなり得るフレームのデータ（Ｉピクチャ及びＰピクチャのＤＣＴ係数データ）が入力される。逆量子化回路１０５は、入力された量子化データ（Ｃ）に対して、その量子化データ（Ｃ）を量子化した際に用いられた量子化スケール値（ｑ＿ｓｃａｌｅ）により、逆量子化処理を行う。
【０００８】
ＩＤＣＴ回路１０６は、逆量子化回路１０５から出力されたＤＣＴ係数データ（ｃｏｅｆ）に対して、逆離散コサイン変換を行い、空間領域の画像データ（ｐｉｘｅｌ）を生成する。
【０００９】
第２の加算回路１０７には、ＩＤＣＴ回路１０６から出力された画像データ（ｐｉｘｅｌ）が入力される。また、第２の加算回路１０７には、入力された画像データがＰピクチャである場合には、その画像データの予測画像データが動き補償回路１０９から入力される。第２の加算回路１０７は、インターマクロブロックに対して処理を行う場合には、入力された画像データに対して予測画像データを加算する。第２の加算回路１０７は、イントラマクロブロックに対して処理を行う場合には、入力された画像データをそのまま出力する。第２の加算回路１０７は、出力した画像データを参照画像データとしてフレーム単位でフレームメモリ１０８に格納する。
【００１０】
動き補償回路１０９は、フレームメモリ１０８に格納された参照画像データに対して動きベクトルを参照して動き補償を行い、予測画像データを生成する。予測画像データは、第１の加算回路１０２に供給される。予測画像データのうち、参照画像となる画像データ（Ｐピクチャの予測画像データ）は、第２の加算回路１０７にも供給される。
【００１１】
以上のような画像符号化装置１０１は、量子化装置１０４から出力される量子化された量子化データ（Ｃ）を、さらに、可変長符号化して、外部へ出力する。
【００１２】
画像復号装置１１１は、図２４に示すように、逆量子化回路１１２と、ＩＤＣＴ回路１１３と、加算回路１１４と、フレームメモリ１１５と、動き補償回路１１６とを備えて構成される。
【００１３】
逆量子化回路１１２には、可変長符号の復号がされたのち、量子化データ（Ｃ）が入力される。量子化回路１１２は、入力された量子化データ（Ｃ）に対して、その量子化データ（Ｃ）を量子化した際に用いられた量子化スケール値（ｑ＿ｓｃａｌｅ）により、逆量子化処理を行い、ＤＣＴ係数データ（ｃｏｅｆ）を生成する。
【００１４】
ＩＤＣＴ回路１１３は、逆量子化回路１１２から出力されたＤＣＴ係数データ（ｃｏｅｆ）に対して、逆離散コサイン変換を行い、空間領域の画像データ（ｐｉｘｅｌ）を生成する。
【００１５】
加算回路１１４には、ＩＤＣＴ回路１１３から出力された画像データ（ｐｉｘｅｌ）が入力される。また、加算回路１１４には、入力された画像データがＰピクチャ又はＢピクチャである場合には、その画像データの予測画像データが動き補償回路１１６から入力される。加算回路１１４は、インターマクロブロックに対して処理を行う場合には、入力された画像データに対して予測画像データを加算する。加算回路１１４は、イントラマクロブロックに対して処理を行う場合には、入力された画像データをそのまま出力する。
【００１６】
フレームメモリ１１５は、加算回路１１４から出力された画像データのうち、参照画像データ（Ｉピクチャ又はＰピクチャ）をフレーム単位で格納する。
【００１７】
動き補償回路１１６は、フレームメモリ１１５に格納された参照画像データに対して動きベクトルを参照して動き補償を行い、予測画像データを生成する。予測画像データは、加算回路１１４に供給される。
【００１８】
以上のような本画像復号装置１１１は、加算回路１１４から出力された画像データを出力する。
【００１９】
ところで、ＭＰＥＧ４（ＩＳＯ／ＩＥＣ　１４４９６−２）規格では、逆量子化の演算式として下記式（１１）及び式（１２）が規定されている。
ｑ＿ｓｃａｌｅが奇数の場合
ｃｏｅｆ　＝　ｓｉｇｎ（Ｃ）　＊　｛（｜Ｃ｜　＊　２　＊　ｑ＿ｓｃａｌｅ　）＋　ｑ＿ｓｃａｌｅ　｝　　…（１１）
ｑ＿ｓｃａｌｅが偶数の場合
ｃｏｅｆ　＝　ｓｉｇｎ（Ｃ）　＊　｛（｜Ｃ｜　＊　２　＊　ｑ＿ｓｃａｌｅ　）＋　ｑ＿ｓｃａｌｅ　−　１｝　　…（１２）
【００２０】
上記式（１１），（１２）において、　“ｑ＿ｓｃａｌｅ”は量子化スケール値、“Ｃ”は量子化後のＤＣＴ係数データ、“ｃｏｅｆ”は量子化前のＤＣＴ係数データ或いは逆量子化後のＤＣＴ係数データである。また、上記式（１１），（１２）において、“Ａ　＋　Ｂ”はＡとＢとを加算した値、“Ａ　−　Ｂ”はＡからＢを減算した値、”Ａ　＊Ｂ”はＡとＢとを乗算した値、“Ａ　／　Ｂ”はＡをＢで除算して小数点以下を切り捨てた値、｜Ａ｜はＡの絶対値、ｓｉｇｎ（Ａ）は絶対値が１で正負符号がＡと同一の値、を意味している。なお、以上の変数の定義、並びに、演算式の定義は、本明細書の発明の詳細な説明の欄の他の式についても同じである。
【００２１】
ＭＰＥＧ４（ＩＳＯ／ＩＥＣ　１４４９６−２）規格に対応した符号化装置では、式（１１），（１２）に示す逆量子化の規定に対応させて、以下の丸め演算式Ａ（式（１３））又は丸め演算式Ｂ（式（１４））の２つの丸め演算式が量子化の演算式として用いられる。
【００２２】
丸め演算式Ａ
Ｃ　＝　ｓｉｇｎ（ｃｏｅｆ）　＊　｛｜　ｃｏｅｆ　｜　　／　（２　＊　ｑ＿ｓｃａｌｅ）　｝　…（１４）
丸め演算式Ｂ
Ｃ　＝　ｓｉｇｎ（ｃｏｅｆ）　＊　｛｜　ｃｏｅｆ　｜　−　（ｑ＿ｓｃａｌｅ　／　２）｝　／　（２　＊　ｑ＿ｓｃａｌｅ）　…（１５）
【００２３】
式（１４）の丸め演算式Ａは、いわゆる１／２丸め式である。つまり、ｃｏｅｆをｑ＿ｓｃａｌｅで除算したときの、余りの値の切り上げ切り下げの閾値を、量子化ステップ（２＊ｑ＿ｓｃａｌｅ）の１／２の値としている方式である。
【００２４】
式（１５）の丸め演算式Ｂは、いわゆる３／４丸め式である。つまり、ｃｏｅｆをｑ＿ｓｃａｌｅで除算したときの、余りの値の切り上げ切り下げの閾値を、量子化ステップ（２＊ｑ＿ｓｃａｌｅ）の３／４の値としている方式である。
【００２５】
丸め演算式Ａと丸め演算式Ｂとを比較した場合、丸め演算式Ｂの方が切り捨て量が多い。つまり、可変長符号化をした後の符号量が少なくなるが、量子化誤差が大きくなる。そのため、符号化装置は、丸め演算式Ａと丸め演算式Ｂとを適応的に切り換えて、ＤＣＴ係数データ（ｃｏｅｆ）の量子化を行う。例えば、発生する符号量は多くてもより誤差を少なく量子化を行いたい場合（例えば、イントラマクロブロックに対して量子化を行う場合）には丸め演算式Ａを適用し、誤差は大きくても発生する符号量を少なくしたい場合（例えば、インターマクロブロックに対して量子化を行う場合）には丸め演算式Ｂを適用する、といったように適応的に演算式を切り換えて量子化を行う。
【００２６】
なお、以上の丸め演算式Ａ及び丸め演算式Ｂは、ＭＰＥＧ４のテストモデル（ＩＳＯ／ＩＥＣ　１４４９６−５　Ｒｅｆｅｒｅｎｃｅ　Ｓｏｆｔｗａｒｅ　）で規定されている
【００２７】
【発明が解決しようとする課題】
本発明は、より効率的に量子化を行う量子化装置及び量子化方法を提供することを目的とする。
【００２８】
【課題を解決するための手段】
本発明にかかる量子化装置及び方法では、整数値である元値（Ｘ）を０又は自然数である所定の量子化スケール値（ｑ＿ｓｃａｌｅ）に基づき量子化し、出力値（Ｙ）を出力する量子化演算を、下記式（１）に基づき出力値（Ｙ）を算出する。
【００２９】
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　｛｜Ｘ｜　−　（ｑ＿ｓｃａｌｅ　＋　１）　／　２｝　／　（２　＊　ｑ＿ｓｃａｌｅ）　…（１）
なお、上記式（１）では、Ａ，Ｂが任意の整数値であるとした場合、“Ａ　＋　Ｂ”はＡとＢとを加算した値、“Ａ　−　Ｂ”はＡからＢを減算した値、”Ａ　＊　Ｂ”はＡとＢとを乗算した値、“Ａ　／　Ｂ”はＡをＢで除算して小数点以下を切り捨てた値、｜Ａ｜はＡの絶対値、ｓｉｇｎ（Ａ）は絶対値が１で正負符号がＡと同一の値、を意味している。
【００３０】
【発明の実施の形態】
以下、本発明の実施の形態として、本発明を適用した量子化装置について説明をする。
【００３１】
なお、本実施の形態の量子化装置は、いわゆるＭＰＥＧ４（ＩＳＯ／ＩＥＣ　１４４９６−２）に規定された符号化データを生成する符号化装置に適用されるものである。すなわち、本実施の形態の量子化装置は、上記式（１１），（１２）に示した逆量子化式に対応した量子化を行う装置である。
【００３２】
第１の実施の形態
まず、第１の実施の形態の量子化装置について説明をする。
【００３３】
図１に、第１の実施の形態の量子化装置１０のブロック構成図を示す。
【００３４】
量子化装置１０は、図１に示すように、量子化スケール発生部１１と、第１の演算部１２と、第２の演算部１３と、セレクタ１４と、コントローラ１５とを備えて構成されている。
【００３５】
量子化装置１０には、離散コサイン変換装置からＤＣＴ係数データ（ｃｏｅｆ）が入力される。入力されたＤＣＴ係数データ（ｃｏｅｆ）は、量子化スケール発生部１１、第１の演算部１２及び第２の演算部１３に供給される。
【００３６】
量子化スケール発生部１１は、入力されたＤＣＴ係数データ（ｃｏｅｆ）のＤＣＴブロック内における画素位置を判断し、その画素位置に対応した量子化スケール値（ｑ＿ｓｃａｌｅ）を発生する。発生される量子化スケール値は、０又は自然数の値である。量子化スケール値（ｑ＿ｓｃａｌｅ）は、第１の演算部１２及び第２の演算部１３に供給される。
【００３７】
第１の演算部１２は、上述した丸め演算式Ａに従い、入力されたＤＣＴ係数データ（ｃｏｅｆ）を量子化する。すなわち、第１の演算部１２は、入力されたＤＣＴ係数データ（ｃｏｅｆ）に対して、下記式（２１）の演算を行い、量子化データ（Ｃ）を求める。
【００３８】
丸め演算式Ａ
Ｃ　＝　ｓｉｇｎ（ｃｏｅｆ）　＊　｛｜　ｃｏｅｆ　｜　　／　（２　＊　ｑ＿ｓｃａｌｅ）　｝　…（２１）
【００３９】
丸め演算式Ａは、いわゆる１／２丸め式である。つまり、ｃｏｅｆをｑ＿ｓｃａｌｅで除算したときの、余りの値の切り上げ切り下げの閾値を、量子化ステップ（２＊ｑ＿ｓｃａｌｅ）の１／２の値としている方式である。
【００４０】
丸め演算式Ａにおけるｑ＿ｓｃａｌｅとＣとの対応関係を、具体的に図２〜図５に示す。図２はｑ＿ｓｃａｌｅ＝３の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図３はｑ＿ｓｃａｌｅ＝４の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図４はｑ＿ｓｃａｌｅ＝５の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図５はｑ＿ｓｃａｌｅ＝６の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示している。
【００４１】
第２の演算部１３は、下記式（２２）に示す丸め演算式Ｃに従い、入力されたＤＣＴ係数データ（ｃｏｅｆ）を量子化する。
【００４２】
丸め演算式Ｃ
Ｃ　＝　ｓｉｇｎ（ｃｏｅｆ）　＊　｛｜　ｃｏｅｆ　｜　−　（ｑ＿ｓｃａｌｅ　＋　１）　／　２｝　／　（２　＊　ｑ＿ｓｃａｌｅ）　…（２２）
【００４３】
式（２２）の丸め演算式Ｃは、いわゆる３／４丸め式である。つまり、ｃｏｅｆをｑ＿ｓｃａｌｅで除算したときの、余りの値の切り上げ切り下げの閾値を、量子化ステップ（２＊ｑ＿ｓｃａｌｅ）の３／４の値としている方式である。
【００４４】
ここで、式（２２）に示す丸め演算式Ｃは、同じ３／４丸め式であっても、式（１５）に示した丸め演算式Ｂとは異なった式となっている。
【００４５】
以下、丸め演算式Ｂと丸め演算式Ｃとを比較する。
【００４６】
丸め演算式Ｂにおけるｑ＿ｓｃａｌｅとＣとの対応関係を、具体的に図６〜図９に示す。図６はｑ＿ｓｃａｌｅ＝３の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図７はｑ＿ｓｃａｌｅ＝４の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図８はｑ＿ｓｃａｌｅ＝５の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図９はｑ＿ｓｃａｌｅ＝６の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示している。
【００４７】
丸め演算式Ｃにおけるｑ＿ｓｃａｌｅとＣとの対応関係を、具体的に図１０〜図１３に示す。図１０はｑ＿ｓｃａｌｅ＝３の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図１１はｑ＿ｓｃａｌｅ＝４の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図１２はｑ＿ｓｃａｌｅ＝５の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図１３はｑ＿ｓｃａｌｅ＝６の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示している。
【００４８】
丸め演算式Ｂでは、量子化スケール値（ｑ＿ｓｃａｌｅ）が偶数の場合には、切り上げの閾値が３／４ちょうどとなる。しかしながら、量子化スケール値（ｑ＿ｓｃａｌｅ）が奇数の場合には、数値として整数値のみ扱うことから、実際のステップの切り上げの閾値が３／４より小さくなってしまう。例えば、ｑ＿ｓｃａｌｅ＝３の場合には、ｃｏｅｆ＝４，１０，１６，・・・で切り上げられるので、余りの値の実際の切り上げ閾値は４／６となる。また、例えば、ｑ＿ｓｃａｌｅ＝５の場合には、ｃｏｅｆ＝６，１６，２６，・・・で切り上げられるので、余りの値の実際の切り上げ閾値は６／１０となる。
【００４９】
これに対して、丸め演算式Ｃの場合、量子化スケール値（ｑ＿ｓｃａｌｅ）が奇数の場合であっても、実際のステップの切り上げの閾値が３／４より大きくなる。例えば、ｑ＿ｓｃａｌｅ＝３の場合には、ｃｏｅｆ＝４，１０，１６，・・・で切り上げられるので、余りの値の実際の切り上げ閾値は４／６となる。また、例えば、ｑ＿ｓｃａｌｅ＝５の場合には、ｃｏｅｆ＝６，１６，２６，・・・で切り上げられるので、余りの値の実際の切り上げ閾値は６／１０となる。
【００５０】
つまり、丸め演算式Ｂと丸め演算式Ｃとを比較した場合、丸め演算式Ｃの方が切り捨てる係数量が多いため、同一の符号列に対してデータの量子化を行うと、可変長符号化後の符号量が少なくなる。すなわち、丸め演算式Ｃは、丸め演算式Ｂと同じ量子化誤差量で丸め演算を行っているにも関わらず、発生する符号量が少なくなる。
【００５１】
セレクタ１４は、コントローラ１５による切り換え制御に従い、第１の演算部１２又は第２の演算部１３のいずれか一方を選択し、選択した演算部から出力された量子化データ（Ｃ）を出力する。
【００５２】
コントローラ１５は、丸め演算式Ａと丸め演算式Ｃとを適応的に切り換えて、量子化データ（Ｃ）を出力する。例えば、発生する符号量は多くてもより誤差を少なく量子化を行いたい場合には第１の演算部１２から出力された量子化データ（Ｃ）（丸め演算式Ａで量子化されたデータ）を選択し、誤差は大きくても発生する符号量を少なくしたい場合には第２の演算部１３から出力された量子化データ（Ｃ）（丸め演算式Ｃで量子化されたデータ）を選択する。
【００５３】
以下、コントローラ１５による丸め演算式Ａと丸め演算式Ｃとの選択条件の第１〜第５の具体例を示す。
【００５４】
（第１例）
第１例では、イントラマクロブロックに含まれるＤＣＴ係数データが入力された場合には丸め演算式Ａを選択し、インターマクロブロックに含まれるＤＣＴ係数データが入力された場合には丸め演算式Ｃを選択する。イントラマクロブロックは、フレーム間の相関を利用せずにフレーム内だけで符号化された画素ブロックである。インターマクロブロックは、時間的に前又は後のフレームとの相関性を利用して圧縮符号化された画素ブロックである。動画像データを圧縮符号化した場合、イントラマクロブロックの画質が劣化すると、他のフレームも伴って画質が劣化することが多い。そのため、第１例では、イントラマクロブロックに対しては、丸め演算式Ａを用いて誤差を少なくした量子化を行い、インターマクロブロックに対しては、丸め演算式Ｃを用いて発生符号量を少なくした量子化を行う。
【００５５】
（第２例）
第２例では、符号化装置全体での圧縮率が所定の閾値（ｔｈ１）より小さい場合には丸め演算式Ａを選択し、圧縮率が所定の閾値（ｔｈ１）以上である場合には丸め演算式Ｃを選択する。圧縮率を小さくするということは、符号量を削減する要求よりも、画質を劣化させない要求の方が大きいということである。そのため、圧縮率がある一定値よりも小さい場合には、常に丸め演算式Ａを用いてより正確な丸め演算を採用する。閾値の基準となる圧縮率は、例えば、原画像信号がいわゆる４：２：０信号であれば、１マクロブロックの発生情報量（１６×１６画素）が平均で２５００ビットを超えると丸め演算式Ａを選択するといったように設定する。
【００５６】
（第３例）
第３例では、上記第１例と第２例とを組み合わせた条件で選択をする。つまり、圧縮率が所定の閾値（ｔｈ１）よりも小さい場合、又は、イントラマクロブロックに含まれるＤＣＴ係数データが入力された場合には、丸め演算式Ａを選択する。それ以外の場合には、丸め演算式Ｃを選択する。
【００５７】
（第４例）
第４例では、丸め演算式Ｃで量子化された後の量子化データ（Ｃ）の絶対値が所定の閾値（ｔｈ２）より大きい場合には丸め演算式Ａを選択し、丸め演算式Ｃで量子化された後の量子化データ（Ｃ）の絶対値が所定の閾値（ｔｈ２）以下の場合には丸め演算式Ｃを選択する。丸め演算式Ｃで量子化する場合には、丸め演算式Ａよりも発生符号量を少なくすることを目的としている。しかしながら、量子化後の量子化データ（Ｃ）がある程度以上大きい場合には、丸め演算式Ａで量子化したときと丸め演算式Ｃで量子化したときとで、可変長符号化した後の符号量が変わらないということがあり得る。この場合、丸め方式Ｃを用いても丸め方式Ａを用いても、発生する符号量が同じであるにも関わらず、量子化誤差が大きくなる。そのため、丸め演算式Ｃで量子化した後の量子化データ（Ｃ）の絶対値が所定の閾値（ｔｈ２）より大きい場合には、丸め演算式Ａを選択し、より正確な丸め演算を採用する。
【００５８】
（第５例）
第５例では、上記第１例、第２例及び第４例を組み合わせた条件で選択をする。つまり、圧縮率が所定の閾値（ｔｈ）よりも小さい場合、イントラマクロブロックに含まれるＤＣＴ係数データが入力された場合、又は、丸め演算式Ｃで量子化された後の量子化データ（Ｃ）の値が所定の閾値（ｔｈ２）より大きい場合には、丸め演算式Ａを選択する。それ以外の場合には、丸め演算式Ｃを選択する。
【００５９】
以上のように第１の実施の形態の量子化装置１０では、３／４丸めの演算式として丸め演算式Ｃを用いている。このため、第１の実施の形態の量子化装置１０では、従来の丸め演算式と同じ精度の丸め演算を行っているのみも関わらず、発生する符号量を少なくすることができる。
【００６０】
第２の実施の形態
つぎに、第２の実施の形態の量子化装置について説明をする。なお、第２の実施の形態の量子化装置の構成要素のうち、上記第１の実施の形態の量子化装置１０の構成要素と同一のものについては、図面中に同一の符号を付ける。
【００６１】
図１４に、第２の実施の形態の量子化装置２０のブロック構成図を示す。
【００６２】
量子化装置２０は、図１４に示すように、量子化スケール発生部１１と、第１の演算部１２と、第２の演算部１３と、第３の演算部２１と、セレクタ２２と、コントローラ２３とを備えて構成されている。
【００６３】
量子化装置２０には、離散コサイン変換装置からＤＣＴ係数データ（ｃｏｅｆ）が入力される。入力されたＤＣＴ係数データは、量子化スケール発生部１１、第１の演算部１２、第２の演算部１３及び第３の演算部２１に供給される。
【００６４】
量子化スケール発生部１１は、入力されたＤＣＴ係数データのＤＣＴブロック内における画素位置を判断し、その画素位置に対応した量子化スケール値（ｑ＿ｓｃａｌｅ）を発生する。発生される量子化スケール値は、０又は自然数の値である。量子化スケール値（ｑ＿ｓｃａｌｅ）は、第１の演算部１２、第２の演算部１３及び第３の演算部２１に供給される。
【００６５】
第１の演算部１２は、丸め演算式Ａに従い、入力されたＤＣＴ係数データ（ｃｏｅｆ）を量子化する。第１の演算部１２により算出された量子化データ（Ｃ）は、セレクタ２２に供給される。
【００６６】
第２の演算部１３は、丸め演算式Ｃに従い、入力されたＤＣＴ係数データ（ｃｏｅｆ）を量子化する。第２の演算部１３により算出された量子化データ（Ｃ）は、セレクタ２２に供給される。
【００６７】
第３の演算部２１は、下記式（２３）に示す丸め演算式Ｄに従い、入力されたＤＣＴ係数データ（ｃｏｅｆ）を量子化する。
【００６８】
丸め演算式Ｄ
Ｃ　＝　ｓｉｇｎ（ｃｏｅｆ）　＊　｛｜　ｃｏｅｆ　｜　−　ｑ＿ｓｃａｌｅ　｝　／　（２　＊　ｑ＿ｓｃａｌｅ）　…（２３）
【００６９】
式（２３）の丸め演算式Ｄは、いわゆる小数点以下の切捨て式である。つまり、ｃｏｅｆをｑ＿ｓｃａｌｅで除算したときの、余りの値を切り捨てている方式である。丸め演算式Ｄは、３／４丸め式である丸め演算式Ｃと比べて、発生符号量が少なくなるが、量子化誤差も大きくなる。
【００７０】
丸め演算式Ｄにおけるｑ＿ｓｃａｌｅとＣとの対応関係を、具体的に図１５〜図１８に示す。図１５はｑ＿ｓｃａｌｅ＝３の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図１６はｑ＿ｓｃａｌｅ＝４の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図１７はｑ＿ｓｃａｌｅ＝５の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示し、図１８はｑ＿ｓｃａｌｅ＝６の場合のｑ＿ｓｃａｌｅとＣとの対応関係を示している。
【００７１】
第３の演算部１４により算出された量子化データ（Ｃ）は、セレクタ２２に供給される。
【００７２】
セレクタ２２は、コントローラ２３による切り換え制御に従い、第１の演算部１２、第２の演算部１３又は第３の演算部２１のいずれか一つを選択し、選択した演算部から出力された量子化データ（Ｃ）を出力する。
【００７３】
コントローラ２３は、丸め演算式Ａ、丸め演算式Ｃ又は丸め演算式Ｄとを適応的に切り換えて、量子化データ（Ｃ）を出力する。例えば、発生する符号量は多くてもより誤差を少なく量子化を行いたい場合には第１の演算部１２から出力された量子化データ（Ｃ）（丸め演算式Ａで量子化されたデータ）を選択し、誤差は大きくても発生する符号量を非常に少なくしたい場合には第３の演算部２１から出力された量子化データ（Ｃ）（丸め演算式Ｄで量子化されたデータ）を選択し、その中間の場合には第２の演算部１３から出力された量子化データ（Ｃ）（丸め演算式Ｃで量子化されたデータ）を選択する
以下、コントローラ２３による丸め演算式Ａ、丸め演算式Ｃ及び丸め演算式Ｄの選択条件の具体例を示す。
【００７４】
（第１例）
第１例では、イントラマクロブロックに含まれるＤＣＴ係数データが入力された場合には丸め演算式Ａを選択し、インターマクロブロックに含まれるＤＣＴ係数データが入力された場合には下記のマトリクス方式で演算式を選択する。
【００７５】
マトリクス方式とは、ＤＣＴ係数データがどの周波数領域に含まれているかに応じて、選択する演算式を切り換える方式である。
【００７６】
一般的に、画像信号の場合、ＤＣＴ係数データの低域成分に対してデータの削除をすると大きな画質劣化となる可能性が大きいが、ＤＣＴ係数データの高域成分に対して情報の削除をしても大きな画質劣化とならない可能性が大きい。マトリクス方式は、以上の画像の特性を利用して、量子化方式を選択する方式である。
【００７７】
マトリクス方式での選択処理例を具体的に説明する。
【００７８】
コントローラ２３は、図１９に示すように、ＤＣＴブロック内を低域領域（図１９中“１”で示した領域）、中域領域（図１９中“２”で示した領域）、高域領域（図１９中“３”で示した領域）に分割したテーブルを備えておく。コントローラ２３は、量子化対象となるＤＣＴ係数データ（ｃｏｅｆ）が入力されると、そのＤＣＴ係数データ（ｃｏｅｆ）がどの領域に含まれているかを示す情報を図１９に示すテーブルから取得する。
【００７９】
コントローラ２３は、入力されたＤＣＴ係数データの周波数領域を示す情報を取得すると、図２０に示すフローチャートに従い処理を行う。
【００８０】
まず、コントローラ２３は、入力されたＤＣＴ係数データが低域領域に含まれているか否かを判断する（ステップＳ１）。コントローラ２３は、入力されたＤＣＴ係数データが低域係数であれば、丸め演算式Ａを選択する。
【００８１】
コントローラ２３は、ステップＳ１で判断した結果、入力されたＤＣＴ係数データが低域係数でなければ、続いて、入力されたＤＣＴ係数データが中域領域に含まれているか否かを判断する（ステップＳ２）。ステップＳ２で判断した結果、入力されたＤＣＴ係数データが中域係数であれば、続いて、丸め演算式Ｃで量子化した後の量子化データ（Ｃ）の値が所定の閾値（ｔｈ３）より大きいか否かを判断する（ステップＳ３）。つまり、丸め演算式Ｃで量子化した場合に情報量の削減の効果が現れるか否かを判断する。コントローラ２３は、丸め演算式Ｃで量子化した後の量子化データ（Ｃ）が所定の閾値（ｔｈ３）より大きい場合には、丸め演算式Ａを選択する。コントローラ２３は、丸め演算式Ｃで量子化した後の量子化データ（Ｃ）が所定の閾値（ｔｈ３）以下の場合には、丸め演算式Ｃを選択する。
【００８２】
コントローラ２３は、ステップＳ２で判断した結果、入力されたＤＣＴ係数データが中域係数でなければ、続いて、丸め演算式Ｄで量子化した後の量子化データ（Ｃ）の値が所定の閾値（ｔｈ４）より大きいか否かを判断する（ステップＳ４）。つまり、丸め演算式Ｄで量子化した場合に情報量の削減の効果が現れるか否かを判断する。コントローラ２３は、丸め演算式Ｄで量子化した後の量子化データ（Ｃ）が所定の閾値（ｔｈ４）より大きい場合には、丸め演算式Ａを選択する。コントローラ２３は、丸め演算式Ｄで量子化した後の量子化データ（Ｃ）が所定の閾値（ｔｈ４）以下の場合には、丸め演算式Ｄを選択する。
【００８３】
（第２例）
第２例は、圧縮率が所定の閾値（ｔｈ１）よりも小さい場合、又は、イントラマクロブロックに含まれるＤＣＴ係数データが入力された場合には、丸め演算式Ａを選択する。それ以外の場合には、上記のマトリクス方式を用いて選択する。
【００８４】
以上のように第２の実施の形態の量子化装置２０では、３／４丸めの演算式として丸め演算式Ｃを用いている。このため、第２の実施の形態の量子化装置２０では、従来の丸め演算式と同じ精度の丸め演算を行っているのみも関わらず、発生する符号量を少なくすることができる。また、さらに、第２の実施の形態の量子化装置２０では、切捨て演算を行う丸め演算式Ｄを用いている。このため、第２の実施の形態の量子化装置２０では、入力されるＤＣＴ係数データに応じて、１／２丸め演算、３／４丸め演算又は切捨て演算を適応的に選択することができる。従って、第２の実施の形態の量子化装置２０では、量子化の精度と符号量の削減とのバランスが取れた量子化を行うことができる。
【００８５】
（変形例）
つぎに、第２の実施の形態の量子化装置２０の変形例について説明をする。
【００８６】
コントローラ２３は、上記マトリクス方式の判断処理フローのステップＳ３及びステップＳ４で、丸め演算式Ｃ又は丸め演算式Ｄで量子化した場合に情報量の削減の効果があるか否かの判断を、量子化データ（Ｃ）と閾値（ｔｈ３，ｔｈ４）とを比較して行っている。この変形例では、この判断で用いられる閾値（ｔｈ３，ｔｈ４）を、可変長符号化後の０ラン数をモニタして、常時可変させる。
【００８７】
図２１に、第２の実施の形態の量子化装置２０を変形した量子化装置３０のブロック構成図を示す。
【００８８】
量子化装置３０には、セレクタ２２から出力される量子化データ（Ｃ）に対して、可変長符号の符号長の評価を行うＶＬＣコード長評価部３１が設けられている。
【００８９】
ＶＬＣコード長評価部３１は、セレクタ２２から出力された量子化データ（Ｃ）の０ランの数を検出し、発生する符号量に影響を与える閾値（ｔｈ３，ｔｈ３）を求める。求めた閾値（ｔｈ３，ｔｈ４）は、コントローラ２３に供給される。
【００９０】
コントローラ２３は、ＶＬＣコード長評価部３１から与えられた閾値（ｔｈ３，ｔｈ４）に基づき、上記マトリクス方式の判断処理を行う。
【００９１】
なお、第１の実施の形態の量子化装置１０でも、コントローラ２３が第３例における条件判断を行う場合には、ＶＬＣコード長評価部３１を設けて閾値（ｔｈ１）の制御を行ってもよい。
【００９２】
また、さらに、上記マトリクス方式の判断処理フローのステップＳ３及びステップＳ４で、丸め演算式Ｃ又は丸め演算式Ｄで量子化した場合に情報量の削減の効果があるか否かの判断を、可変長符号化後のデータを直接モニタして行ってもよい。この場合には、図２２に示すように、第１の演算部１２、第２の演算部１３及び第３の演算部２１のそれぞれの後段に、第１の可変長符号化（ＶＬＣ）部４１、第１の可変長符号化（ＶＬＣ）部４２及び第３の可変長符号化部４３をそれぞれ設け、セレクタ２２は、可変長符号化した後の量子化データ（Ｃ）を選択するように構成する。このように構成すれば、コントローラ２２は、可変長符号化後のデータを直接モニタして比較することができるので、より正確に、丸め演算式ごとの符号量の削減の効果の影響を判断することができる。
【００９３】
【発明の効果】
本発明にかかる量子化装置及び方法では、整数値である元値（Ｘ）を０又は自然数である所定の量子化スケール値（ｑ＿ｓｃａｌｅ）に基づき量子化し、出力値（Ｙ）を出力する量子化演算を、下記式（１）に基づき出力値（Ｙ）を算出する。
【００９４】
Ｙ　＝　ｓｉｇｎ（Ｘ）　＊　｛｜Ｘ｜　−　（ｑ＿ｓｃａｌｅ　＋　１）　／　２｝　／　（２　＊　ｑ＿ｓｃａｌｅ）　…（１）
本発明にかかる量子化装置及び方法では、効率的に量子化を行うことができる。
【図面の簡単な説明】
【図１】本発明の第１の実施の形態の量子化装置のブロック構成図である。
【図２】丸め方式Ａ及びｑ＿ｓｃａｌｅ＝３で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図３】丸め方式Ａ及びｑ＿ｓｃａｌｅ＝４で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図４】丸め方式Ａ及びｑ＿ｓｃａｌｅ＝５で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図５】丸め方式Ａ及びｑ＿ｓｃａｌｅ＝６で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図６】丸め方式Ｂ及びｑ＿ｓｃａｌｅ＝３で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図７】丸め方式Ｂ及びｑ＿ｓｃａｌｅ＝４で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図８】丸め方式Ｂ及びｑ＿ｓｃａｌｅ＝５で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図９】丸め方式Ｂ及びｑ＿ｓｃａｌｅ＝６で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１０】丸め方式Ｃ及びｑ＿ｓｃａｌｅ＝３で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１１】丸め方式Ｃ及びｑ＿ｓｃａｌｅ＝４で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１２】丸め方式Ｃ及びｑ＿ｓｃａｌｅ＝５で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１３】丸め方式Ｃ及びｑ＿ｓｃａｌｅ＝６で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１４】本発明の第２の実施の形態の量子化装置のブロック構成図である。
【図１５】丸め方式Ｄ及びｑ＿ｓｃａｌｅ＝３で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１６】丸め方式Ｄ及びｑ＿ｓｃａｌｅ＝４で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１７】丸め方式Ｄ及びｑ＿ｓｃａｌｅ＝５で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１８】丸め方式Ｄ及びｑ＿ｓｃａｌｅ＝６で量子化した場合のｑ＿ｓｃａｌｅとＣとの対応関係を示す図である。
【図１９】上記第２の実施の形態の量子化装置のコントローラに備えられているテーブルを示す図である。
【図２０】上記第２の実施の形態の量子化装置で用いられる量子化方式の選択処理を示すフローチャートである。
【図２１】上記第２の実施の形態の第１の変形例を示すブロック構成図である。
【図２２】上記第２の実施の形態の第２の変形例を示すブロック構成図である。
【図２３】ＭＰＥＧ２又はＭＰＥＧ４の符号化装置のブロック構成図である。
【図２４】ＭＰＥＧ２又はＭＰＥＧ４の復号装置のブロック構成図である。
【符号の説明】
１０，２０　量子化装置、１１　量子化スケール発生部、１２　第１の演算部、１３　第２の演算部、１４，２２　セレクタ、１５，２３　コントローラ、２１　第３の演算部
【図１】
2004048191.tif 000003
【図２】
2004048191.tif 000004
【図３】
2004048191.tif 000005
【図４】
2004048191.tif 000006
【図５】
2004048191.tif 000007
【図６】
2004048191.tif 000008
【図７】
2004048191.tif 000009
【図８】
2004048191.tif 000010
【図９】
2004048191.tif 000011
【図１０】
2004048191.tif 000012
【図１１】
2004048191.tif 000013
【図１２】
2004048191.tif 000014
【図１３】
2004048191.tif 000015
【図１４】
2004048191.tif 000016
【図１５】
2004048191.tif 000017
【図１６】
2004048191.tif 000018
【図１７】
2004048191.tif 000019
【図１８】
2004048191.tif 000020
【図１９】
2004048191.tif 000021
【図２０】
2004048191.tif 000022
【図２１】
2004048191.tif 000023
【図２２】
2004048191.tif 000024
【図２３】
2004048191.tif 000025
【図２４】
2004048191.tif 000026