(57)【要約】【課題】  符号化効率を向上させるための技術を提供すること。【解決手段】  着目フレームの画像における着目画素ブロックに対してフレーム間符号化を行う場合、該着目フレームより１フレーム過去のフレームにおける画像において着目画素ブロックともっとも類似する画像領域に対してフィルタ処理を行う。そして、該フィルタ処理を行った画像領域と着目画素ブロックとの差分に対してフレーム間符号化を行う。【選択図】  図４
2013004992.tif 000002
【特許請求の範囲】
【請求項１】
動画像を構成する各フレームの画像について、該画像を複数の画素ブロックに分割し、該画素ブロックを整数画素精度のハイブリッド符号化を用いて符号化する符号化手段を有する画像処理装置であって、
前記符号化手段は、着目フレームの画像における着目画素ブロックに対してフレーム間符号化を行う場合、該着目フレームより１フレーム過去のフレームにおける画像において前記着目画素ブロックともっとも類似する画像領域に対してフィルタ処理を行い、該フィルタ処理を行った画像領域と前記着目画素ブロックとの差分に対してフレーム間符号化を行う
ことを特徴とする画像処理装置。
【請求項２】
前記符号化手段は、前記着目フレームの画像を撮像したときのシャッタースピードを示すシャッタースピード情報を取得し、該取得したシャッタースピード情報が示すシャッタースピードが遅いほどローパスフィルタ処理の影響が小さくなるローパスフィルタを用いて前記画像領域に対するローパスフィルタ処理を行うことを特徴とする請求項１に記載の画像処理装置。
【請求項３】
前記符号化手段は、
前記着目画素ブロック及び前記差分のそれぞれについてＭＡＤを計算し、
前記着目画素ブロックのＭＡＤが前記差分のＭＡＤよりも小さい場合、前記着目画素ブロックに対するフレーム内符号化の結果を前記着目画素ブロックの符号化結果として出力し、
前記差分のＭＡＤが前記着目画素ブロックのＭＡＤよりも小さい場合、前記差分に対するフレーム間符号化の結果を前記着目画素ブロックの符号化結果として出力する
ことを特徴とする請求項１又は２に記載の画像処理装置。
【請求項４】
前記符号化手段は、
前記１フレーム過去のフレームにおける画像中に、前記着目画素ブロックと同サイズの画像領域を複数個設定し、該設定した複数個の画像領域のそれぞれに対して、該画像領域から前記着目画素ブロックへの動きベクトルの成分に応じた指向性を有する指向性フィルタを用いてフィルタ処理を行い、
前記フィルタ処理を施した前記複数個の画像領域のうち、前記着目画素ブロックともっとも類似する画像領域を特定し、該特定した画像領域と前記着目画素ブロックとの差分に対してフレーム間符号化を行う
ことを特徴とする請求項１に記載の画像処理装置。
【請求項５】
前記符号化手段は、
前記１フレーム過去のフレームにおける画像中に、前記着目画素ブロックと同サイズの画像領域を複数個設定し、該設定した複数個の画像領域のそれぞれに対して、該画像領域から前記着目画素ブロックへの動きベクトルの大きさが整数画素単位、非整数画素単位の何れであるのかのケースに対して予め設定されたフィルタを用いてフィルタ処理を行い、
前記フィルタ処理を施した前記複数個の画像領域のうち、前記着目画素ブロックともっとも類似する画像領域を特定し、該特定した画像領域と前記着目画素ブロックとの差分に対してフレーム間符号化を行う
ことを特徴とする請求項１に記載の画像処理装置。
【請求項６】
動画像を構成する各フレームの画像について、該画像を複数の画素ブロックに分割し、該画素ブロックを整数画素精度のハイブリッド符号化を用いて符号化する符号化手段を有する画像処理装置が行う画像処理方法であって、
前記符号化手段が、着目フレームの画像における着目画素ブロックに対してフレーム間符号化を行う場合、該着目フレームより１フレーム過去のフレームにおける画像において前記着目画素ブロックともっとも類似する画像領域に対してフィルタ処理を行い、該フィルタ処理を行った画像領域と前記着目画素ブロックとの差分に対してフレーム間符号化を行う
ことを特徴とする画像処理方法。
【請求項７】
コンピュータを、請求項１乃至６の何れか１項に記載の画像処理装置の符号化手段として機能させるためのコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、符号化技術に関するものである。
【背景技術】
【０００２】
映像の符号化において、ＤＣＴと動き補償とを組み合わせた、いわゆるハイブリッド符号化方式の初期の国際標準には、Ｈ．２６１（非特許文献１）がある。Ｈ．２６１では、整数画素単位で動き補償を行っていた。
【０００３】
ＭＰＥＧ−１（非特許文献２）では、半画素精度の動きベクトルを算出して動き補償を行うように改良された。これにより、原画と参照画像の差分である残差の発生が抑制され、残差の符号化にともなう発生符号量が減少する。
【０００４】
Ｈ．２６４（非特許文献３）ではさらに、１／４画素精度の動き補償を行ようになった。しかし、動き補償における画素精度を細かくすることは、同一速度の動きであってもそれを表現する動きベクトルの情報量は増加することになり、その動きベクトルを符号化した際の発生符号量も増加することになる。つまり、高画素精度の動きベクトルを利用することにより、残差の発生符号量を削減しても、動きベクトルの発生符号量は増加してしまうために、トータルの符号量削減にはつながらないケースがある。とくに、映像の高解像度化が進むと、解像度に比例して動きベクトルの情報量も増加するため、動きベクトルの発生符号量と残差の発生符号量の双方を同時に削減する技術が重要となってくる。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】ＩＴＵ−Ｔ  Ｈ．２６１
【非特許文献２】ＩＳＯ／ＩＥＣ  １１１７２−２
【非特許文献３】ＩＴＵ−Ｔ  Ｈ．２６４（ＩＳＯ／ＩＥＣ  １４４９６  Ｐａｒｔ１０
【発明の概要】
【発明が解決しようとする課題】
【０００６】
本発明は以上の問題に鑑みてなされたものであり、符号化効率を向上させるための技術を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明の目的を達成するために、例えば、本発明の画像処理装置は、動画像を構成する各フレームの画像について、該画像を複数の画素ブロックに分割し、該画素ブロックを整数画素精度のハイブリッド符号化を用いて符号化する符号化手段を有する画像処理装置であって、前記符号化手段は、着目フレームの画像における着目画素ブロックに対してフレーム間符号化を行う場合、該着目フレームより１フレーム過去のフレームにおける画像において前記着目画素ブロックともっとも類似する画像領域に対してフィルタ処理を行い、該フィルタ処理を行った画像領域と前記着目画素ブロックとの差分に対してフレーム間符号化を行うことを特徴とする。
【発明の効果】
【０００８】
本発明の構成によれば、符号化効率を向上させることができる。
【図面の簡単な説明】
【０００９】
【図１】符号化部２００１の機能構成例を示すブロック図。
【図２】復号部２００４の機能構成例を示すブロック図。
【図３】シンタックス表現の一部を示す図。
【図４】画像処理装置の機能構成例を示すブロック図。
【図５】符号化部２００１の機能構成例を示すブロック図。
【図６】復号部２００４の機能構成例を示すブロック図。
【図７】符号化部２００１の機能構成例を示すブロック図。
【図８】復号部２００４の機能構成例を示すブロック図。
【発明を実施するための形態】
【００１０】
以下、添付図面を参照し、本発明の好適な実施形態について説明する。なお、以下説明する実施形態は、本発明を具体的に実施した場合の一例を示すもので、特許請求の範囲に記載した構成の具体的な実施例の一つである。
【００１１】
［第１の実施形態］
本実施形態では、動画像を構成する各フレームの画像を複数の画素ブロックに分割し、画素ブロックごとに、整数精度の動きベクトルを用いたハイブリッド符号化を行う画像処理装置について説明する。先ず、画像処理装置の機能構成例を図４のブロック図を用いて説明する。
【００１２】
動画像を構成する各フレームの画像は順次、符号化部２００１に入力される。符号化部２００１は、入力されたフレームの画像を複数の画素ブロックに分割し、外部から取得した情報（外部情報）を用いてそれぞれの画素ブロックに対してハイブリッド符号化を行い、画素ブロックごとの符号化データ（ストリーム）を生成する。ここで、外部情報とは、各フレームの画像サイズ、動画像フォーマット、フレームレート、などの動画像の属性、シャッタースピードなどの撮影条件、符号化条件、ユーザ設定パラメータなどを示す情報である。
【００１３】
そして符号化部２００１は、画素ブロックごとのストリームと、符号化部２００１内で生成したサイド情報と、後述するフィルタ情報と、を統合してストリーム統合部２００２に送出する。ここで、符号化部２００１内で生成されるサイド情報は、外部情報として入力された情報を用いて画素ブロックを符号化する際に必要な情報であり、既定のフォーマットにしたがって成形した情報である。
【００１４】
ストリーム統合部２００２は、外部情報を入力してサイド情報を生成すると共に、符号化部２００１が生成したサイド情報、ブロック単位のストリーム、フィルタ情報、と統合してストリームパース部２００３に出力する。ここで、ストリーム統合部２００２で生成されるサイド情報は、画像全体を表すのに必要な情報を既定のフォーマットにしたがって成形した情報である。なお、ストリーム統合部２００２からストリームパース部２００３は直接的に接続されている必要はなく、ネットーワークや光メディアなどを介して伝送するようにしてもよい。
【００１５】
ストリームパース部２００３は、ストリームを画素ブロック単位のデータに分解し、サイド情報及びフィルタ情報と共に、復号部２００４に出力する。復号部２００４は、入力したストリームを復号して各フレームの画像を構成して出力する。出力先については特に限定するものではなく、ＣＲＴや液晶画面などの表示装置に対して出力してもよいし、メモリなどの記憶装置に対して出力してもよいし、ネットワークを介して外部の装置に対して送信してもよい。
【００１６】
次に、符号化部２００１の機能構成例について、図１のブロック図を用いて説明する。本実施形態では、動画像を構成する最初のフレームと、それから一定間隔毎のフレーム又は任意の位置におけるフレームをイントラフレームとし、イントラフレームを構成する各画素ブロック（８画素×８画素）をイントラ画素ブロックと呼称する。そしてイントラ画素ブロックを符号化する場合は、周知の如く、インターフレーム（イントラフレーム以外のフレーム）における画素ブロックは参照せず、インターフレーム内での処理（フレーム内符号化）となる。また、周知の如く、フレーム毎に、画素ブロック間の動きベクトルの計算も行われる。
【００１７】
然るに、符号化部２００１は、イントラフレームの画像が入力された場合は、この画像に対してフレーム内符号化を行う。一方、符号化部２００１は、イントラフレーム以外のフレーム（インターフレーム）の画像が入力された場合は、以下に説明する各処理を実行する。即ち、以下の説明は、インターフレームの画像が入力された場合における符号化部２００１の動作についてのものである。
【００１８】
動画像を構成する各フレームの画像は順次、入力フレームバッファ１０１に格納される。以降では、入力フレームバッファ１０１に格納されているインターフレーム（着目フレーム）の画像（Ｄ１０１）に対する符号化動作について説明する。然るに以下に説明する符号化動作は、入力フレームバッファ１０１に格納されている各インターフレームの画像について行われることになる。
【００１９】
イントラデータ生成部１０２は、入力フレームバッファ１０１から着目フレームの画像を読み出し、読み出した着目フレームの画像を複数個の画素ブロックに分割する。以降では、イントラデータ生成部１０２が分割した複数個の画素ブロックのうち着目画素ブロックに対する符号化動作について説明する。然るに以下に説明する符号化動作は、着目フレームの画像を構成する各画素ブロックについて行われることになる。
【００２０】
動き探索部１０９は、入力フレームバッファ１０１に格納されている着目フレームの画像から、着目画素ブロック（Ｄ１０１Ａ）を読み出す。そして動き探索部１０９は、フレームバッファ１１０に格納されている、着目フレームよりも１フレーム前の画像において、着目画素ブロックともっとも類似する画像領域を、参照画素ブロック（Ｄ１１０）として特定する。例えば、着目フレームの画像中の着目画素ブロックの領域に対応する、着目フレームよりも１フレーム前の画像中の領域を規定倍したサイズを有する領域を探索領域として設定する。そしてこの探索領域内で、着目画素ブロックとのＳＡＤ（絶対値差分和）が最も小さい画像領域を参照画素ブロックとして特定する。ＳＡＤは２つの画像の歪み度を示す指標であり、ＳＡＤが低いほど２つの画像は似ていることになる。そして動き探索部１０９は、参照画素ブロックから着目画素ブロックへのベクトルを動きベクトル（Ｄ１０９Ａ）として求め、求めた動きベクトルを可変長符号化部１０８及びフィルタ１１３に対して送出する。
【００２１】
フィルタ１１３はフレームバッファ１１０から、着目フレームよりも１フレーム前の画像において、動きベクトルにより特定される参照画素ブロック（Ｄ１１０）を読み出す。そしてフィルタ１１３は、この読み出した参照画素ブロックを構成する各画素について、該画素の周辺の画素を用いた処理を行う。本実施形態では、参照画素ブロックに対して、下記の表１に示す成分値を有するローパスフィルタを用いたローパスフィルタ処理を行う。
【００２２】
【表１】
2013004992.tif 000003
【００２３】
表１においてｎは、フィルタ強度を示す変数であり、ｎが大きいほど、ローパスフィルタ処理の影響が小さくなるローパスフィルタとなる。このｎはフィルタ情報算出部１５０によって算出される。
【００２４】
フィルタ情報算出部１５０は、フレーム毎に、外部情報を用いてこのｎを算出する。例えば、外部情報に含まれている、着目フレームの画像を撮像したときのシャッタースピードを示すシャッタースピード情報を取得し、該取得したシャッタースピード情報が示すシャッタースピードが早いほどｎを小さく設定する。例えば以下の式１に従ってｎを決定する。
【００２５】
【数１】
2013004992.tif 000004
【００２６】
この式１においてαは定数であり、レンズの明るさなどの固有な特長によって予め定められている。また、ｓはシャッタースピード情報が示すシャッタースピードである。例えば、αを７２０とした場合、ｓ＝１／６０（ｓｅｃ）の場合にはｎ＝１２、ｓ＝１／１２０（ｓｅｃ）の場合にはｎ＝６となり、シャッタースピードが速いほど（撮像ボヤケが少ないほど）、強いローパスフィルタがかかることになる。
【００２７】
フィルタ情報算出部１５０は、上記のｎを求めると、この求めたｎをフィルタ１１３に対して送出し、更にこのｎをフィルタ情報（Ｄ１５０）として符号化してストリーム統合部２００２に送出する（サイド情報の一部）。
【００２８】
然るにフィルタ１１３は、フィルタ情報算出部１５０から送出されたｎによりローパスフィルタを確定させ、確定させたローパスフィルタを用いて、参照画素ブロックに対するローパスフィルタ処理を行う。より具体的には、表１に示すフィルタカーネルを用いて参照画素ブロックに対する畳み込み積分を行うことで、フィルタ処理済み参照画素ブロック（Ｄ１１３Ｂ）を生成する。そしてフィルタ１１３は、この生成したフィルタ処理済み参照画素ブロック（Ｄ１１３，Ｄ１１３Ｂ）を、動き補償部１１１及び残差データ生成部１１４に対して送出する。
【００２９】
残差データ生成部１１４は、入力フレームバッファ１０１に格納されている着目フレームの画像から、着目画素ブロック（Ｄ１０１Ａ）を読み出す。そして残差データ生成部１１４は、この着目画素ブロックと、フィルタ処理済み参照画素ブロックと、の差分ブロックを予測残差（Ｄ１１４）として計算する。即ち、差分ブロックにおける画素位置（ｘ、ｙ）における画素値は、着目画素ブロックの画素位置（ｘ、ｙ）における画素値とフィルタ処理済み参照画素ブロックの画素位置（ｘ、ｙ）における画素値との差分値となる。
【００３０】
モード判定部１０３は、イントラデータ生成部１０２からの着目画素ブロック、残差データ生成部１１４からの予測残差、のそれぞれのＭＡＤ（絶対差分平均）を算出する。ＭＡＤは画素ブロック内の画素値の平均と各画素値との差の絶対値和であり、この値が小さいほど、画素ブロック内の画素値の分散が少ない、つまりエントロピーが低く、ブロックを符号化した場合に符号量が小さくなりやすいことを示す。
【００３１】
モード判定部１０３は、着目画素ブロックのＭＡＤ＞予測残差のＭＡＤの場合には、動作モードをインターモードに設定すると共に、残差データ生成部１１４からの予測残差を予測データ（Ｄ１０３）としてＤＣＴ部１０４に送出する。一方、モード判定部１０３は、着目画素ブロックのＭＡＤ＜予測残差のＭＡＤの場合には、動作モードをイントラモードに設定し、イントラデータ生成部１０２からの着目画素ブロックを予測データ（Ｄ１０３）としてＤＣＴ部１０４に送出する。
【００３２】
ＤＣＴ部１０４は、予測データに対して直交変換を施して、周波数成分ごとにパワーを有するＤＣＴ係数（Ｄ１０４）に変換し、量子化部１０５に出力する。
【００３３】
量子化部１０５は、ＤＣＴ係数の周波成分ごとに異なる量子化ステップで量子化を行い、量子化データ（Ｄ１０５）を生成する。そして量子化部１０５は、この量子化データを可変長符号化部１０８及び逆量子化部１０６に対して送出する。
【００３４】
逆量子化部１０６は、量子化部１０５からの量子化データに対して、逆量子化（量子化部１０５による量子化の逆変換）を行うことで、ＤＣＴ係数（Ｄ１０６）を復元し、復元したＤＣＴ係数を逆ＤＣＴ部１０７に対して送出する。
【００３５】
逆ＤＣＴ部１０７は、逆量子化部１０６からのＤＣＴ係数に対して逆ＤＣＴ変換（ＤＣＴ部１０４によるＤＣＴ変換の逆変換）を行うことで予測データ（復号予測データ）を復元し、この復号予測データを動き補償部１１１に対して送出する。
【００３６】
動き補償部１１１は、モード判定部１０３によりイントラモードが設定されている場合、復号予測データをそのまま再構成ブロック（Ｄ１１１）としてフレームバッファ１１０に格納する。一方、動き補償部１１１は、モード判定部１０３によりインターモードが設定されている場合、復号予測データと、フィルタ１１３からのフィルタ処理済み参照画素ブロックと、を加算して再構成ブロックを生成する。この場合、再構成ブロックの画素位置（ｘ、ｙ）における画素値は、｛復号予測データにおける画素位置（ｘ、ｙ）の画素値｝＋｛フィルタ処理済み参照画素ブロックの画素位置（ｘ、ｙ）における画素値｝となる。そして動き補償部１１１は、この加算処理で生成した再構成ブロックをフレームバッファ１１０に格納する。
【００３７】
可変長符号化部１０８は、量子化部１０５による量子化結果及びモード判定部１０３が設定したモードを示す情報をハフマン符号化し、動き探索部１０９が求めた動きベクトルをＧｏｌｏｍｂ符号化することで、画素ブロックごとのストリームを生成する。そして可変長符号化部１０８は、画素ブロックごとにこのようにして生成したストリームと、サイド情報と、フィルタ情報、とを１つのストリーム（Ｄ１０８）としてストリーム統合部２００２に入力する。
【００３８】
次に、以上の符号化処理による符号化結果を復号して、各フレームの復号画像を得るための復号処理について説明する。
【００３９】
ストリームパース部２００３は、ストリーム統合部２００２からのストリームを、サイド情報、画素ブロック単位のストリーム、フィルタ情報に分解し、復号部２００４に供給する。
【００４０】
復号部２００４の機能構成例を図２のブロック図に示す。可変長復号部２０１は、着目ブロックのストリーム（Ｄ２００）を受けると、これを可変長復号し、量子化データ（Ｄ２０１）については逆量子化部２０２へと出力し、動きベクトル（Ｄ２０１Ａ）についてはフィルタ部２０６に出力する。また、モードを示す情報は、後述する各機能部に供給される。ここで、この着目画素ブロックは上記の通り、上記の着目フレーム（インターフレーム）における画素ブロックである。
【００４１】
逆量子化部２０２は、量子化データを逆量子化し、ＤＣＴ係数（Ｄ２０２）を生成し、逆ＤＣＴ部２０３へ出力する。逆ＤＣＴ部２０３は、このＤＣＴ係数に対して逆ＤＣＴを行うことで復号予測データ（Ｄ２０３）を生成し、動き補償部２０５へ出力する。
【００４２】
フィルタ部２０６は先ず、モードを示す情報を参照する。そしてこのモードを示す情報がイントラモードである場合には動作せず、インターモードである場合には、図１のフィルタ１１３と同様の処理を行う。即ち、ストリームパース部２００３からフィルタ情報としてのフィルタ強度ｎ（Ｄ２５０）を入力し、動きベクトルが指す参照画素ブロック（Ｄ２０４Ａ）をフレームバッファ２０４から取得し、表１のフィルタカーネルに基づいてフィルタ処理を行う。そしてフィルタ部２０６は、このフィルタ処理によって得られるフィルタ処理済み参照画素ブロック（Ｄ２０６）を、動き補償部２０５に送出する。
【００４３】
動き補償部２０５は先ず、モードを示す情報を参照する。そしてこのモードを示す情報がイントラモードである場合には、逆ＤＣＴ部２０３からの復号予測データを着目画素ブロックの復号結果（Ｄ２０５）としてフレームバッファ２０４に格納する。一方、モードを示す情報がインターモードである場合、動き補償部２０５は、フィルタ部２０６からのフィルタ処理済み参照画素ブロックを復号予測データに加算して再構成ブロック（Ｄ２０５）を生成し、フレームバッファ２０４に格納する。フレームバッファ２０４は、着目フレームの画像を構成する全ての画素ブロックの復号結果が格納された時点で、これらの画素ブロックをフレームデータ（Ｄ２０４）として出力する。
【００４４】
Ｈ．２６１では、整数画素精度の動き探索および動き補償において、参照画素ブロックにフィルタをかけていなかった。そのため、着目画素ブロックと参照画素ブロックの位相が小数画素単位でずれた映像を整数画素単位で動き補償する場合、整数画素単位未満でずれた画像ブロックの差分をとることになり、残差成分に不要な高周波が発生していた。これが整数画素単位の動き補償において、符号化効率が向上しない一因である。
【００４５】
本実施形態では、動きベクトルを整数精度に保ち、動きベクトルの符号量の増加を抑制しながら、フィルタ処理を行った参照画素ブロックとの間で動き補償残差成分を生成することで、残差成分の高周波の発生を抑制し、トータルの符号化効率を狙うものである。
【００４６】
なお、本実施形態では、低精度の画素単位を整数画素としているが、それに限るものではなく、例えば、２画素単位で行うことも可能である。また、本実施形態では、イントラデータ生成部１０２によるイントラデータとして画素ブロックそのものを用いているが、それに限るものではなく、例えば、Ｈ．２６４で採用されるようなフレーム内予測を用いてもよい。
【００４７】
また、本実施形態では、動き探索において、ＳＡＤを用いて動きベクトルを決定しているがこれに限るものではない。例えば、ＳＳＤ（Ｓｕｍ  ｏｆ  Ｓｑｕａｒｅ  Ｄｉｆｆｅｒｅｎｃｅ）やＳＡＴＤ（Ｓｕｍ  ｏｆ  Ａｂｓｏｌｕｔｅ  Ｔｒａｎｓｆｏｒｍ  Ｄｉｆｆｅｒｅｎｃｅ）といった画素ブロックの相違度を求める指標を利用することも可能である。
【００４８】
また、動きベクトルの決定方法は本実施形態で説明した方法に限定するものではなく、全探索とよばれる動き探索方法を採用してもよいし、全探索以外の動き探索アルゴリズムを採用してもよい。
【００４９】
また、本実施形態では、モードの決定にはＭＡＤを用いているが、画像のエントロピー示す指標であればこれに限定されない。例えば、実際にＤＣＴ、量子化、可変長符号化を行い、符号化データのサイズと復号画像の画質の比からモードを決定してもよい。
【００５０】
また、本実施形態では、可変長符号化にはハフマン符号やＧｏｌｏｍｂ符号を用いているがこれに限るものではなく、コンテキスト適応算術符号化など他の可変長符号化方法を用いてもよい。
【００５１】
また、本実施形態では、画素ブロックのサイズを８画素×８画素としているが、これに限るものではなく、４画素×４画素や１６画素×１６画素、あるいは、画素ブロックのサイズを可変としてそれらを組み合わせる構成をとってもよい。その場合のＤＣＴや量子化もブロックのサイズに応じて異なる。
【００５２】
また、本実施形態では、フィルタ情報算出部１５０は、フィルタ強度ｎをサイド情報として出力しているが、これに限るものではなく、例えば、フィルタカーネルをそのまま出力してもよい。また、フィルタカーネルを決定する情報を出力する単位も、フレーム単位に限るものではなく、スライス、ブロック等の単位毎にサイド情報として出力してもよい。
【００５３】
また、フィルタ強度ｎは、シャッタースピード情報が示すシャッタースピードｓの値が大きくなる（シャッタースピードは遅くなる）のに従って大きくなるのであれば、式１以外の式に従って求めてもよい。例えば、ｎ＝β×ｓ２（βは定数）という式に従ってフィルタ強度ｎを求めてもよい。
【００５４】
また、本実施形態ではフィルタ強度ｎは、シャッタースピード情報をもとに算出しているが、これに限るものではなく、外部情報としてユーザが指定したパラメータそのものを利用してフィルタ強度ｎを算出してもよい。
【００５５】
また、本実施形態では、フィルタ情報としてフィルタ強度ｎを符号化したが、これに限るものではなく、シャッタースピード情報を直接符号化し、復号側でも同様にフィルタ強度ｎを算出してもかまわない。また、本実施形態では、直交変換の一例としてＤＣＴ変換を挙げたが、アダマール変換やＷａｖｅｌｅｔ変換などを適用してもよい。
【００５６】
［第２の実施形態］
本実施形態では、動画像を構成する各フレームの画像を複数の画素ブロックに分割し、画素ブロックごとに、整数精度の動きベクトルを用いたハイブリッド符号化を行う画像処理装置について説明する。本実施形態と第１の実施形態との差分は、符号化部２００１及び復号部２００４である。然るに以下では、符号化部２００１及び復号部２００４について説明する。
【００５７】
本実施形態に係る符号化部２００１の機能構成例を図５に示す。図５において、図１に示した構成要件と同じ構成要件については同じ参照番号を付しており、それらに係る説明は省略する。
【００５８】
動き探索部５０９は、着目フレームよりも１フレーム前の画像から着目画素ブロックともっとも類似する画像領域を特定するために、先ず、着目画素ブロックからある方向への規定長の動きベクトルを設定する。これは、例えば、着目画素ブロックから他の画素ブロックへの動きベクトルであってもよい。そして動き探索部５０９は、この設定した動きベクトル（Ｄ５０９Ｂ）をフィルタ決定部５０１に送出する。
【００５９】
フィルタ決定部５０１は、動き探索部５０９から受けた動きベクトルの成分に応じた指向性フィルタを生成すべく、この成分に応じた指向性フィルタカーネルを生成する。この生成方法については後述する。そしてフィルタ決定部５０１は、この生成した指向性フィルタカーネル（Ｄ５０１Ｂ）と、動き探索部５０９から受けた動きベクトル（Ｄ５０９Ｂ）と、をフィルタ５１３に送出する。
【００６０】
フィルタ５１３は、フレームバッファ１１０に格納されている、着目フレームよりも１フレーム前の画像から、フィルタ決定部５０１から受けた動きベクトルによって特定される画像領域を特定する。これは、画像中における着目画素ブロックの位置から動きベクトル分だけずれた位置の画像領域を特定することで達成される。そしてフィルタ５１３は、この特定した画像領域に対して、フィルタ決定部５０１から受けた指向性フィルタカーネルを用いて畳み込み積分を行い、フィルタ処理済み参照画素ブロックを生成する。そしてフィルタ５１３は、この生成したフィルタ処理済み参照画素ブロック（Ｄ５１３Ｂ）を動き探索部５０９に送出する。
【００６１】
動き探索部５０９は、フィルタ５１３から受けたフィルタ処理済み参照画素ブロックと着目画素ブロックとのＳＡＤ（絶対値差分和）を計算する。このように、次の（１）〜（４）の処理を、方向、長さが異なる複数個の動きベクトルのそれぞれについて行う。
【００６２】
（１）  着目画素ブロックからある方向への規定長の動きベクトルを設定
（２）  動きベクトルの成分に応じた指向性フィルタカーネルを生成
（３）  動きベクトルによって特定される画像領域に対して、指向性フィルタカーネルを用いて畳み込み積分を行い、フィルタ処理済み参照画素ブロックを生成
（４）  フィルタ処理済み参照画素ブロックと着目画素ブロックとのＳＡＤ（絶対値差分和）を計算
そして、最も小さいＳＡＤが計算された画像領域を参照画素ブロックとする。また、動き探索部５０９は、最も小さいＳＡＤが計算された動きベクトル（Ｄ１０９Ａ）をフィルタ決定部５０１及び可変長符号化部１０８に対して送出する。
【００６３】
フィルタ決定部５０１は、動き探索部５０９から受けた動きベクトルの成分に応じた指向性フィルタカーネルを生成し、生成した指向性フィルタカーネル（Ｄ５０１Ｂ）と動き探索部５０９から受けた動きベクトル（Ｄ５０９Ｂ）とをフィルタ５１３に送出する。
【００６４】
フィルタ５１３は、フレームバッファ１１０に格納されている、着目フレームよりも１フレーム前の画像から、フィルタ決定部５０１から受けた動きベクトルによって特定される画像領域を特定する。そしてフィルタ５１３は、この特定した画像領域に対して、フィルタ決定部５０１から受けた指向性フィルタカーネルを用いて畳み込み積分を行い、フィルタ処理済み参照画素ブロックを生成する。そしてフィルタ５１３は、この生成したフィルタ処理済み参照画素ブロック（Ｄ１１３）を動き補償部１１１に送出する。
【００６５】
次に、フィルタ決定部５０１による、指向性フィルタカーネルの生成方法について説明する。動き探索部５０９から受けた動きベクトルの成分を（ｍｖｘ，ｍｖｙ）として表現した場合、フィルタ決定部５０１は、表２に記した成分値を有する、コンボリューションフィルタカーネルである仮指向性フィルタカーネルｋ（ｘ，ｙ）を生成する。
【００６６】
【表２】
2013004992.tif 000005
【００６７】
表２に記したマトリックスは、フィルタカーネルに指向性を持たせるマトリックスである。例えば、動きベクトルの成分が（−２，−２）の場合は、仮指向性フィルタカーネルは表３（ａ）のように斜め方向に成分をもち、（−３，０）の場合は、表３（ｂ）のように横方向に成分をもつ。
【００６８】
【表３】
2013004992.tif 000006
【００６９】
次にフィルタ決定部５０１は、この仮指向性フィルタカーネルｋ（ｘ，ｙ）を用いて、以下の式２に従って指向性フィルタカーネルｋ’（ｘ，ｙ）を生成する。
【００７０】
【数２】
2013004992.tif 000007
【００７１】
この式においてｎは第１の実施形態と同様にして求められるフィルタ強度であり、フィルタの強度を設定しつつカーネルを正規化するものである。例えば、ｎ＝１２、動きベクトルの成分が（−２，−２）の場合、指向性フィルタカーネルは表４（ａ）のように斜め方向にぼかしを加えるフィルタカーネルとなり、（−３，０）の場合、表４（ｂ）のように水平方向にぼかしをくわえるフィルタカーネルとなる。
【００７２】
【表４】
2013004992.tif 000008
【００７３】
また、動きベクトルの長さが小さい場合、例えば、（０，０）の場合、指向性フィルタカーネルは表５のようになり、このフィルタカーネルを使って畳み込み積分を行うことはフィルタリング処理を行わないことと同義となる。なお、このフィルタ強度ｎは第１の実施形態と同様、フィルタ情報（Ｄ１５０）として符号化される。
【００７４】
【表５】
2013004992.tif 000009
【００７５】
以降は第１の実施形態と同様である。即ち、本実施形態では、着目フレームより１フレーム過去のフレームの画像中に、着目画素ブロックと同サイズの画像領域を複数個設定する。そして、該設定した複数個の画像領域のそれぞれについて、該画像領域から着目画素ブロックへの動きベクトルの成分に応じた指向性を有する指向性フィルタを用いて該画像領域に対してフィルタ処理を行う。そして、フィルタ処理を施した複数個の画像領域のうち、着目画素ブロックともっとも類似する画像領域を特定し、該特定した画像領域と着目画素ブロックとの差分に対してフレーム間符号化を行う。
【００７６】
次に、以上の符号化処理による符号化結果を復号して、各フレームの復号画像を得るための復号処理について説明する。復号部２００４の機能構成例を図６のブロック図に示す。図６に示した構成要件のうち、図２に示した構成要件と同じ構成要件については同じ参照番号を付しており、その説明は省略する。
【００７７】
可変長復号部２０１により復号された動きベクトル（Ｄ２０１Ａ）はフィルタ決定部６０１に入力される。また、フィルタ情報としてのフィルタ強度ｎ（Ｄ２５０）もこのフィルタ決定部６０１に入力される。フィルタ決定部６０１は、図５のフィルタ決定部５０１の動作と同様に動作し、入力された動きベクトルから指向性フィルタカーネル（Ｄ６０１）を生成し、動きベクトル（Ｄ２０１Ａ）と共に、フィルタ部６０６へ出力する。フィルタ部６０６は、図５のフィルタ５１３と同一の動作を行う。
【００７８】
このように、動きベクトルの方向と長さを加味したフィルタリング処理を行うことにより、着目フレームと参照画素ブロックのずれにより発生する指向性をもった高周波成分をより一層抑制することができる。また、動き探索におけるＳＡＤの算出においても、フィルタをかけた参照画素ブロックを用いて算出することにより、より残差成分の発生の少ない動きベクトルを算出することができる。
【００７９】
なお、本実施形態では、フィルタ強度ｎを固定値と扱っているが、第１の実施形態と同様に可変とし、サイド情報として復号側に伝送してもよい。また、第１の実施形態で説明した様々な変形例も適宜本実施形態に適用することができる。
【００８０】
［第３の実施形態］
本実施形態では、動画像を構成する各フレームの画像を複数の画素ブロックに分割し、画素ブロックごとに、整数精度と小数精度の動きベクトルを用いたハイブリッド符号化を行う画像処理装置について説明する。以下では、第１の実施形態と異なる点のみについて説明し、それ以外については第１の実施形態に準ずるものとする。
【００８１】
本実施形態に係る符号化部２００１の機能構成例を図７に示す。図７において、図１に示した構成要件と同じ構成要件については同じ参照番号を付しており、それらに係る説明は省略する。
【００８２】
本実施形態では、動き探索部７０９とコントロール部７２０とフィルタ７１３が連動して、動きベクトル、フィルタ強度、フィルタタイプの組み合わせを試行する。そしてこれにより、最終的にＳＡＤが最も小さくなる組合わせを決定し、最終動きベクトル、最終フィルタ強度、最終ｃｔｒｌ＿ｃｏｄｅを決定する。このとき、フィルタ強度は、値として「４」、「１２」、「フィルタなしを示すＮ／Ａ識別信号」の何れかをとるものとする。また、動きベクトル、フィルタ強度、フィルタタイプの組わせは、表７で表現される四種類が存在し、それぞれにｃｔｒｌ＿ｃｏｄｅがコードが割り当てられる。
【００８３】
動き探索部７０９は、第１の実施形態及び第２の実施形態と同様にして、最小のＳＡＤをとる動きベクトルを算出するものであるが、本実施形態では、整数画素単位だけでなく、半画素単位の動きベクトルも探索対象とする。また、動きベクトルの大きさは半画素（１画素の半分）の長さを１単位とする。例えば、動きベクトル（２，０）は、ｘ軸のプラス方向に１画素分の長さを持つ動きベクトルを表すものとなる。
【００８４】
動き探索部７０９は、着目フレームよりも１フレーム前の画像から着目画素ブロックともっとも類似する画像領域を特定するために、先ず、着目画素ブロックからある方向への規定長の動きベクトルを設定する。これは、例えば、着目画素ブロックから他の画素ブロックへの動きベクトルであってもよい。そして動き探索部７０９は、この設定した動きベクトル（Ｄ７０９Ｂ）を、コントロール部７２０に入力する。
【００８５】
コントロール部７２０は、動き探索部７０９から受けた動きベクトル（ｔｍｖｘ，ｔｍｖｙ）、フィルタ強度、ｃｔｒｌ＿ｃｏｄｅの組み合わせをフィルタ７１３に出力する。
【００８６】
フィルタ７１３は、動きベクトル、フィルタ強度、ｃｔｒｌ＿ｃｏｄｅから予測画素ブロックを生成し、動き探索部７０９へ出力する。
【００８７】
動き探索部７０９は、フィルタ７１３から受けた予測画素ブロックと着目画素ブロックとのＳＡＤ（絶対値差分和）を計算する。
【００８８】
このように、動き探索部７０９は、異なる規定長の動きベクトルを再設定しながら、ＳＡＤが最も小さくなる動きベクトル、フィルタ強度、ｃｔｒｌ＿ｃｏｄｅの組み合わせを探索していく。また、ＳＡＤが最も小さくなった際の動きベクトル、フィルタ強度、ｃｔｒｌ＿ｃｏｄｅをそれぞれ、最終動きベクトル、最終フィルタ強度、最終ｃｔｒｌ＿ｃｏｄｅとする。
【００８９】
動き探索部７０９は、最終動きベクトル、最終フィルタ強度、最終ｃｔｒｌ＿ｃｏｄｅを用いて最終フィルタ処理済み参照画素ブロックを生成する。そして動き探索部７０９は、最終フィルタ処理済み参照画素ブロックと入力フレームバッファ１０１に格納されている着目画素ブロックとの差分である予測残差（Ｄ７１４）をモード判定部１０３に送出する。また、最終動きベクトル、最終フィルタ強度、最終ｃｔｒｌ＿ｃｏｄｅをコントロール部７２０へ出力する。
【００９０】
以下、動きベクトル、フィルタ強度、ｃｔｒｌ＿ｃｏｄｅの最適な組み合わせの具体的な算出手順について記述する。本実施形態では、Ｈ．２６４方式でも採用されている予測動きベクトルを用いる。着目画素ブロックの左、上、左上のそれぞれの画素ブロックの動きベクトルのｘ、ｙ成分のそれぞれの中央値をとったベクトルが予測動きベクトルとなる。本実施形態では、予測動きベクトルを（ｐｍｖｘ，ｐｍｖｙ）と表現する。
【００９１】
【数３】
2013004992.tif 000010
【００９２】
この等式が成立する場合、一つの動きベクトルに対し、ｃｔｒｌ＿ｃｏｄｅに０，１，３のケースでそれぞれＳＡＤを算出を試行する。なお、Ｌ・」は、床関数を表し、Ｌｖａｌｕｅ／２」×２は、ｖａｌｕｅ値を２単位ずつの値に変換する、つまり、整数画素精度に変換することを意味する。この等式が成立しない場合、一つの動きベクトルに対し、ｃｔｒｌ＿ｃｏｄｅを２、３のケースでそれぞれＳＡＤを算出を試行する。また、コントロール部７２０が、動きベクトル（ｔｍｖｘ，ｔｍｖｙ）、フィルタ強度、ｃｔｒｌ＿ｃｏｄｅの組み合わせをフィルタ７１３へ出力する際、ｃｔｒｌ＿ｃｏｄｅが１の場合にはフィルタ強度１２を出力する。また、ｃｔｒｌ＿ｃｏｄｅが２の場合にはフィルタ強度４を、それ以外の場合は、Ｎ／Ａ識別信号を出力する。なお、フィルタ強度はこの値に限定されず、例えば、シーケンスやフレームごとに別な値を設定してもよい。
【００９３】
フィルタ７１３内におけるフィルタ切り替え部７２３は、コントロール部７２０から受けたｃｔｒｌ＿ｃｏｄｅに基づいて、コンボリューションフィルタ７２１、補間フィルタ７２２、コピー処理部７２４の何れを動作させるのかを決定する。
【００９４】
ｃｔｒｌ＿ｃｏｄｅが０の場合には、コピー処理部７２４を動作さる。ｃｔｒｌ＿ｃｏｄｅが１または２の場合には、入力されたフィルタ強度でコンボリューションフィルタ７２１を動作させる。ｃｔｒｌ＿ｃｏｄｅが３の場合、補間フィルタ７２２を動作させる。
【００９５】
コンボリューションフィルタ７２１は、動きベクトルが示す参照画素ブロックを構成する全ての画素に対して、表１に示すフィルタカーネル（フィルタ強度を用いて確定済み）を用いて畳み込み積分を行うことでフィルタ処理済み参照画素ブロックを生成する。そしてコンボリューションフィルタ７２１は、このフィルタリングを施したフィルタ処理済み参照画素ブロック（Ｄ５１３Ｂ）を動き探索部７０９へ出力する。
コピー処理部７２４は、動きベクトルが示す参照画素ブロックをそのままコピーしたものをフィルタ処理済み参照画素ブロック（Ｄ５１３Ｂ）として、動き探索部７０９へ出力する。
【００９６】
補間フィルタ７２２は、表６に示す如く、動きベクトルのｘ成分値が整数画素単位であるか非整数画素単位であるか、ｙ成分値が整数画素単位であるか非整数画素単位であるか、の組み合わせに応じた生成方法でもって補間画像を生成する。
【００９７】
【表６】
2013004992.tif 000011
【００９８】
動きベクトルのｘ成分値＝非整数画素単位、ｙ成分値＝整数画素単位である場合、補間画像上の位置（ｘ、ｙ）における画素値を、参照画素ブロック上の位置ｄ（ｘ、ｙ）における画素値とｄ（ｘ＋１、ｙ）における画素値との平均値とする。
【００９９】
動きベクトルのｘ成分値＝整数画素単位、ｙ成分値＝非整数画素単位である場合、補間画像上の位置（ｘ、ｙ）における画素値を、参照画素ブロック上の位置ｄ（ｘ、ｙ）における画素値とｄ（ｘ、ｙ＋１）における画素値との平均値とする。
【０１００】
動きベクトルのｘ成分値＝非整数画素単位、ｙ成分値＝非整数画素単位である場合、補間画像上の位置（ｘ、ｙ）の画素値を、参照画素ブロック上の位置ｄ（ｘ、ｙ）、ｄ（ｘ＋１、ｙ）、ｄ（ｘ、ｙ＋１）、ｄ（ｘ＋１、ｙ＋１）の画素値の平均値とする。
【０１０１】
つまり、座標値が半画素を示すものであれば、隣接する画素の平均をとり、画素補間を行って値をつくりだすことを示している。そして補間フィルタ７２２は、このようにして生成された補間画像をフィルタ処理済み参照画素ブロック（Ｄ５１３Ｂ）として動き探索部７０９へ出力する。
【０１０２】
また、コントロール部７２０は、最終的動きベクトルから差分動きベクトル（ｄｍｖｘ，ｄｍｖｙ）を算出する。最終動きベクトルを（ｍｖｘ、ｍｖｙ）とすると、差分動きベクトルは以下の式で算出できる。
【０１０３】
【数４】
2013004992.tif 000012
【０１０４】
コントロール部７２０は、ｃｔｒｌ＿ｃｏｄｅと差分動きベクトルを可変長符号化部７０８に出力する。
【０１０５】
以下の表７は、シンタックスから、復号時の動きベクトルとフィルタ強度ｎを得るための表である。これは、ｃｔｒｌ＿ｃｏｄｅに応じて、動きベクトルの可変長符号化方法とフィルタの強度が連動することを示す。表７によれば、差分動きベクトルが（０，０）の場合でも、コントロール符号に応じて、フィルタの有無を切替ることができるようになっている。また、同じ整数精度の動きベクトルであっても、差分動きベクトルが（０，０）の場合はフィルタ強度ｎが大きくなり、差分動きベクトルが（０，０）でない場合に比べて、弱いフィルタをかけることになる。
【０１０６】
【表７】
2013004992.tif 000013
【０１０７】
可変長符号化部７０８は、値としてｃｔｒｌ＿ｃｏｄｅと差分動きベクトルとに対して、図３に基づいた可変長符号化を行い、モード判定結果を示すデータ（Ｄ１０３Ａ）や量子化データ（Ｄ１０５）を符号化したデータとともに、外部に出力する。
【０１０８】
図３は、画素ブロック単位の動きベクトルとコントロール符号の符号化・復号のためのシンタックス表現の一部を示す図である。なお、図３のｕ（２）、ｓｅ（ｖ）、ａｅ（ｖ）という表現はそれぞれ、２ビットフラグ、符号ありＧｏｌｍｂ符号、コンテキストアダプティブ算術圧縮を示し、Ｈ．２６４の規格書で示される表記に準拠している。Ｈ．２６４の数値は、サイズ固定ビット・Ｇｏｌｍｂ符号のみで圧縮、または、サイズ固定ビット・Ｇｏｌｍｂ符号とコンテキスト適応算術圧縮を切り替えて圧縮される。ｓｅ（ｖ）｜ａｅ（ｖ）という表現は、Ｇｏｌｍｂ符号とコンテキスト適応算術圧縮を切り替えて圧縮することを示している。コンテキスト適応算術圧縮利用は、ピクチャ単位でコンテキスト適応算術圧縮利用識別フラグを埋め込み、そのフラグを参照することにより切り替えられる。つまり、ｓｅ（ｖ）｜ａｅ（ｖ）と表記されたケースにおいて、コンテキスト適応算術圧縮利用識別フラグがＯｎであればａｅ（ｖ）を利用し、そうでない場合にはｓｅ（ｖ）を利用することになる。
【０１０９】
動き補償部１１１は最終動きベクトル、最終フィルタ強度、最終ｃｔｒｌ＿ｃｏｄｅで生成したフィルタ処理済み参照画素ブロック（Ｄ１１３）と復号予測データ（Ｄ１０７）を加算して再構築ブロック（Ｄ１１１）を生成しフレームバッファ１１０に格納する。
【０１１０】
次に、以上の符号化処理による符号化結果を復号して、各フレームの復号画像を得るための復号処理について説明する。復号部２００４の機能構成例を図８のブロック図に示す。図８に示した構成要件のうち、図２に示した構成要件と同じ構成要件については同じ参照番号を付しており、その説明は省略する。
【０１１１】
可変長復号部２０１は、図３に示されるシンタックスにしたがい、ｃｔｒｌ＿ｃｏｄｅと差分動きベクトル（ｔｍｖｘ，ｔｍｖｙ）を可変長復号し、フィルタ情報算出部８２３に出力する。フィルタ情報算出部８２３は、ｃｔｒｌ＿ｃｏｄｅと差分動きベクトル（ｔｍｖｘ，ｔｍｖｙ）から、表７に従い動きベクトル（Ｄ８２３）とフィルタ強度（Ｄ８２３Ａ）を決定し、出力する。
【０１１２】
フィルタ切替部８０６、コンボリューションフィルタ８２０、補間フィルタ８２１、コピー処理部８２２の動作はそれぞれ、第１の実施形態と同様である。即ち、それぞれは、フィルタ切り替え部７２３、コンボリューションフィルタ７２１、補間フィルタ７２２、コピー処理部７２４の、動き補償部１１１へ参照画素ブロックの出力動作と同一である。  動きベクトルの可変長符号化方法、フィルタ処理方法を個別に制御し、その制御方法を示す符号を個々に復号側へ伝送すると、個々に符号が発生し、高い符号化効率を望みにくい。しかし、本実施形態のように可変長符号化方法とフィルタ処理方法が連動して制御すると、その制御方法を一つのコントロール符号として復号側に送ることにより、符号量の発生を最小限におさえつつ、最適な処理を選択でき、高い符号化効率を得ることができる。
【０１１３】
なお、本実施形態では、画素ブロック単位でコントロール符号を埋め込んでいるが、これに限るものではなく、スライス、フレーム、シーケンス単位で埋め込みを行ってもよい。
【０１１４】
また、本実施形態ではフィルタ強度（Ｄ７２０）の値として３種類の値を用いて説明したが、この数の種類に限るものではなく、これよりも多くの種類もしくはこれよりも少ない数の種類であってもよい。
【０１１５】
［第４の実施形態］
図１，２，４，５，６，７，８に示した各部は何れも専用の回路によって構成してもよいが、汎用プロセッサを用いて構成してもよい。この場合、いくつかの機能部を１つの汎用プロセッサで実装してもよい。
【０１１６】
また、フレームバッファについてはＲＡＭやＲＯＭ等のメモリで構成し、その他の各部についてはコンピュータプログラムで構成してもよい。この場合、このメモリを有するコンピュータのＣＰＵがこのコンピュータプログラムを実行することで、このコンピュータは、上記の各実施形態に係る画像処理装置として機能することになる。
【０１１７】
（その他の実施例）
また、本発明は、以下の処理を実行することによっても実現される。即ち、上述した実施形態の機能を実現するソフトウェア（プログラム）を、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータ（またはＣＰＵやＭＰＵ等）がプログラムを読み出して実行する処理である。
【図１】
2013004992.tif 000014
【図２】
2013004992.tif 000015
【図３】
2013004992.tif 000016
【図４】
2013004992.tif 000017
【図５】
2013004992.tif 000018
【図６】
2013004992.tif 000019
【図７】
2013004992.tif 000020
【図８】
2013004992.tif 000021