(57)【要約】  イントラＬＭ予測画像または照度予測画像の精度を向上させるための画像予測装置を提供する。  画像予測装置は、線形予測パラメータａとｂを導出する線形予測パラメータ導出部と、上記線形予測パラメータに基づいて予測画像を生成する線形予測部を備える。上記線形予測パラメータ導出部は、第１パラメータａ１と第２パラメータａ２から線形予測パラメータａを導出する手段を備える。上記線形予測パラメータ導出部は、上記第１パラメータａ１と所定の閾値ＴＨＮを比較し、上記第１パラメータａ１が上記閾値ＴＨＮ未満もしくは上記閾値ＴＨＮ以下である場合には第１コストａ１ｃｏｓｔＮを上記第１パラメータａ１から減算する手段と、上記第１パラメータａ１と所定の閾値ＴＨＰを比較し、上記第１パラメータａ１が上記閾値ＴＨＰ未満もしくは上記閾値ＴＨＰ以下の場合には第２コストａ１ｃｏｓｔＰを上記第１パラメータａ１に加算する手段の少なくとも１つの手段を備える。
2017061189.tif 000002
【特許請求の範囲】
【請求項１】
インデックスｉに対応する画素値xiと画素値yiの組を入力として、xiからyiを予測する線形予測パラメータａとｂを導出する線形予測パラメータ導出部と、上記線形予測パラメータに基づいて予測画像を生成する線形予測部を備える画像予測装置であって、上記線形予測パラメータ導出部は、画素値xiと画素値yiの積の和XYと画素値xiの和Xと画素値yiの和Yの積に基づいて導出される第１パラメータa1と、画素値xiと画素値xiの積の和XXと画素値xiの和Xと画素値xiの和Xの積に基づいて導出される第２パラメータa2から線形予測パラメータａを導出する手段を備え、上記線形予測パラメータ導出部は、上記第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN未満もしくは上記閾値THN以下である場合には第１コストa1costNを上記第１パラメータa1から減算する手段と、上記第１パラメータa1と所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THP未満もしくは上記閾値THP以下の場合には第２コストa1costPを上記第１パラメータa1に加算する手段の少なくとも１つの手段を備えることを特徴とする画像予測装置。
【請求項２】
上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THNより小さい場合、もしくは、上記第１パラメータa1が上記所定の閾値THN以下の場合に、上記第１パラメータa1に上記第１コストa1costNを減算することを特徴とする請求項１に記載の画像予測装置。
【請求項３】
上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THPより大きい場合、もしくは、上記第１パラメータa1が上記所定の閾値THP以上の場合に、上記第１パラメータa1に上記第２コストa1costPを加算することを特徴とする請求項１に記載の画像予測装置。
【請求項４】
上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THNより小さい場合、もしくは、上記第１パラメータa1が上記所定の閾値THN以下の場合に、上記第１パラメータa1に上記第１コストa1costNを減算し、上記第１パラメータa1が上記所定の閾値THPより大きい場合、もしくは、上記第１パラメータa1が上記所定の閾値THP以上の場合に、上記第１パラメータa1に上記第２コストa1costPを加算することを特徴とする請求項１に記載の画像予測装置。
【請求項５】
上記線形予測パラメータ導出部は、上記第２パラメータa2にコストa2costを加算する手段をさらに備えることを特徴とする請求項１に記載の画像予測装置。
【請求項６】
上記線形予測パラメータ導出部は、上記第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN未満もしくは上記閾値THN以下である場合には第１コストa1costNを上記第１パラメータa1から減算し、上記第２パラメータa2にコストa2costNを加算する手段と、上記第１パラメータa1と所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THP未満もしくは上記閾値THP以下の場合には第２コストa1costPを上記第１パラメータa1に加算し上記第２パラメータa2にコストa2costPを加算する手段の少なくとも１つの手段を備えることを特徴とする請求項１に記載の画像予測装置。
【請求項７】
上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THNより小さい場合、もしくは、上記第１パラメータa1が上記所定の閾値THN以下の場合に、上記第１パラメータa1に上記第１コストa1costNを減算し、上記第２パラメータa2にコストa2costNを加算する手段を備えることを特徴とする請求項６に記載の画像予測装置。
【請求項８】
上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THPより大きい場合、もしくは、上記第１パラメータa1が上記所定の閾値THP以上の場合に、上記第１パラメータa1に上記第２コストa1costPを加算し、上記第２パラメータa2にコストa2costPを加算する手段を備えることを特徴とする請求項６に記載の画像予測装置。
【請求項９】
上記コストa2costは、上記コストa1costN以上であることを特徴とする請求項５に記載の画像予測装置。
【請求項１０】
上記線形予測パラメータ導出部は、上記閾値THNとTHPのどちらか一方もしくは両方を０とすることを特徴とする請求項１に記載の画像予測装置。
【請求項１１】
上記線形予測パラメータ導出部は、上記閾値THNとTHPのどちらか一方もしくは両方を、上記第２パラメータa2に基づいて導出することを特徴とする請求項１に記載の画像予測装置。
【請求項１２】
上記線形予測パラメータ導出部は、上記閾値THNを０、上記閾値THPを上記第２パラメータa2に基づいて導出することを特徴とする請求項１に記載の画像予測装置。
【請求項１３】
上記線形予測パラメータ導出部は、上記閾値THNを、上記第２パラメータa2を所定の定数だけ右シフトすることにより導出する、ことを特徴とする請求項６に記載の画像予測装置。
【請求項１４】
上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記画素値xiと画素値xiの積の和XXから導出することを特徴とする請求項１に記載の画像予測装置。
【請求項１５】
上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記画素値xiと画素値yiの積の和XYから導出することを特徴とする請求項１に記載の画像予測装置。
【請求項１６】
上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、画素のビットデプス値から導出することを特徴とする請求項１に記載の画像予測装置。
【請求項１７】
上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記第２パラメータa2に基づいて導出することを特徴とする請求項１に記載の画像予測装置。
【請求項１８】
上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記第２パラメータa2を所定の定数だけ右シフトすることにより、導出することを特徴とする請求項１に記載の画像予測装置。
【請求項１９】
上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記第２パラメータa1を所定の定数だけ右シフトした値と、画素のビットデプス値から導出される値との最小値から、導出することを特徴とする請求項１に記載の画像予測装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、画像復号装置、画像符号化装置、および予測画像生成装置に関する。
【背景技術】
【０００２】
動画像を効率的に伝送または記録するために、動画像を符号化することによって符号化データを生成する動画像符号化装置、および、当該符号化データを復号することによって復号画像を生成する動画像復号装置が用いられている。【０００３】
具体的な動画像符号化方式としては、例えば、Ｈ．２６４／ＭＰＥＧ−４．ＡＶＣやＨＥＶＣ（High-Efficiency Video Coding）にて提案されている方式などが挙げられる。【０００４】
このような動画像符号化方式においては、動画像を構成する画像（ピクチャ）は、画像を分割することにより得られるスライス、スライスを分割することにより得られる符号化単位（コーディングユニット（Coding Unit）と呼ばれることもある）、及び、符号化単位を分割することより得られるブロックである予測ユニット（ＰＵ）、変換ユニット（ＴＵ）からなる階層構造により管理され、ブロックごとに符号化／復号される。【０００５】
また、このような動画像符号化方式においては、通常、入力画像を符号化／復号することによって得られる局所復号画像に基づいて予測画像が生成され、当該予測画像を入力画像（原画像）から減算して得られる予測残差（「差分画像」または「残差画像」と呼ぶこともある）が符号化される。予測画像の生成方法としては、画面間予測（インター予測）、および、画面内予測（イントラ予測）が挙げられる。【０００６】
非特許文献１では、画面内予測において、既に復号した輝度成分と色差成分の関係から、線形予測パラメータを推定し、推定された線形予測パラメータを用いて対象ブロックの色差成分を予測する、方法が知られている。【０００７】
非特許文献２、非特許文献３には、さらに、線形予測パラメータの導出において、二乗誤差を最小化する最小二乗法に加え、正則化コストを加えることで、線形予測パラメータを安定させる技術が知られている。
【先行技術文献】
【非特許文献】
【０００８】
【非特許文献１】J. Kim, S.-W. Park, J.-Y. Park, B.-M. Jeon (LG), “Intra chroma prediction using inter channel correlation”, JCTVC-B021, JCT-VC 2nd Meeting: Geneva, CH, 21-28 July, 2010
【非特許文献２】T. Ikai, “3D-CE5.h related: Illumination compensation regression improvement and simplification”, JCT3V-D0061, JCT3V 4th Meeting: Incheon,KR, 20-26 Apr. 2013（2013年4月12日公開）
【非特許文献３】J. Chen, Y. Chen, M. Karczewicz, X. Li, H. Liu, L. Zhang, X. Zhao, “Coding tools investigation for next generation video coding”, ITU-T SG16 Doc. COM16-C806, Feb. 2015.
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかしながら、非特許文献１では、線形予測パラメータを安定させるような正則化コストを用いないため、導出される線形予測パラメータがノイズの影響を受けやすい、という課題があった。非特許文献２では、正則化コストを用いるが、特に線形予測式の傾きのパラメータの符号が正にも負にもなる異なる色空間（例えば輝度から色差及び色差から別の色差）を予測する色コンポーネント間予測に適した正則化コストとは限らないという課題があった。非特許文献３では、色差から別の色差を予測する場合の線形予測に用いる正則化コストは開示されているが、線形予測の傾きパラメータの符号や大きさを考慮しておらず、特に、傾きが正となることがある場合に適当ではないという課題があった。このように、従来は、色成分間の線形予測パラメータを推定し、推定した線形予測パラメータを用いて予測画像を生成する予測画像生成装置において、推定される線形予測パラメータが雑音の影響を受けやすいという課題があった。さらに、特に、線形予測式の傾きのパラメータの符号が正にも負にもなる異なる色コンポーネント間予測において推定される予測パラメータが適当ではないという課題があった。
【課題を解決するための手段】
【００１０】
本発明の一つの形態は、インデックスｉに対応する画素値xiと画素値yiの組を入力として、xiからyiを予測する線形予測パラメータａとｂを導出する線形予測パラメータ導出部と、上記線形予測パラメータに基づいて予測画像を生成する線形予測部を備える画像予測装置であって、上記線形予測パラメータ導出部は、画素値xiと画素値yiの積の和XYと画素値xiの和Xと画素値yiの和Yの積に基づいて導出される第１パラメータa1と、画素値xiと画素値xiの積の和XXと画素値xiの和Xと画素値xiの和Xの積に基づいて導出される第２パラメータa2から線形予測パラメータａを導出する手段を備え、上記線形予測パラメータ導出部は、上記第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN未満もしくは上記閾値THN以下である場合には第１コストa1costNを上記第１パラメータa1から減算する手段と、上記第１パラメータa1と所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THP未満もしくは上記閾値THP以下の場合には第２コストa1costPを上記第１パラメータa1に加算する手段の少なくとも１つの手段を備えることを特徴とする。
【発明の効果】
【００１１】
本発明によれば、推定される線形予測パラメータが雑音に強くなり、線形予測パラメータを用いて導出される予測画像の精度が向上する。線形予測を利用した画像復号装置及び画像符号化装置において符号化効率が向上する。
【図面の簡単な説明】
【００１２】
【図１】本実施形態に係るＬＭ予測部３０９３の構成を示すブロック図である。
【図２】本実施形態に係る符号化ストリームのデータの階層構造を示す図である。
【図３】参照ピクチャリストの一例を示す概念図である。
【図４】参照ピクチャの例を示す概念図である。
【図５】本実施形態に係る画像復号装置の構成を示す概略図である。
【図６】本実施形態に係るインター予測パラメータ復号部の構成を示す概略図である。
【図７】本実施形態に係るマージ予測パラメータ導出部の構成を示す概略図である。
【図８】本実施形態に係るＡＭＶＰ予測パラメータ導出部の構成を示す概略図である。
【図９】ベクトル候補の一例を示す概念図である。
【図１０】本実施形態に係るインター予測パラメータ復号制御部の構成を示す概略図である。
【図１１】本実施形態に係るインター予測画像生成部の構成を示す概略図である。
【図１２】本実施形態に係るＬＭ予測の概念図である。
【図１３】本実施形態に係る照度補償の概念図である。
【図１４】本実施形態に係る画像符号化装置の構成を示すブロック図である。
【図１５】本実施形態に係るインター予測パラメータ符号化部の構成を示す概略図である。
【図１６】本実施形態に係る照度補償で用いるテーブルを示す図である。
【図１７】本実施形態に係る線形予測パラメータａの分布を示すヒストグラムである。
【図１８】本実施形態に係るＬＭパラメータａ導出部３０９３１６の構成を示すブロック図である。
【図１９】本実施形態に係るＬＭ正則化コスト加算部３０９３１８の構成を示すブロック図である。
【図２０】本実施形態に係るＬＭ正則化コスト加算部３０９３１８の構成を示すブロック図である。
【図２１】本実施形態に係るＬＭ正則化コスト加算部３０９３１８Ａ、ＬＭ正則化コスト加算部３０９３１８Ａ２の動作を示すフローチャートである。
【図２２】本実施形態に係るＬＭ正則化コスト加算部３０９３１８Ｂ、ＬＭ正則化コスト加算部３０９３１８Ｂ２の動作を示すフローチャートである。
【図２３】本実施形態に係るＬＭ正則化コスト加算部３０９３１８Ｃ、ＬＭ正則化コスト加算部３０９３１８Ｃ２の動作を示すフローチャートである。
【図２４】本実施形態に係る照度補償部３０９３の構成を示すブロック図である。
【図２５】本発明の実施形態に係る画像伝送システムの構成を示す概略図である。
【発明を実施するための形態】
【００１３】
（第１の実施形態）
以下、図面を参照しながら本発明の実施形態について説明する。
【００１４】
図２５は、本実施形態に係る画像伝送システム１の構成を示す概略図である。
【００１５】
画像伝送システム１は、複数のレイヤ画像を符号化した符号を伝送し、伝送された符号を復号した画像を表示するシステムである。画像伝送システム１は、画像符号化装置１１、ネットワーク２１、画像復号装置３１及び画像表示装置４１を含んで構成される。
【００１６】
画像符号化装置１１には、単一レイヤもしくは複数のレイヤ画像を示す信号Ｔが入力される。レイヤ画像とは、ある解像度及びある視点で視認もしくは撮影される画像である。複数のレイヤ画像を用いて３次元画像を符号化するビュースケーラブル符号化を行う場合、複数のレイヤ画像のそれぞれは、視点画像と呼ばれる。また、空間スケーラブル符号化を行う場合、その複数のレイヤ画像は、解像度の低いベースレイヤ画像と、解像度の高い拡張レイヤ画像からなる。ＳＮＲスケーラブル符号化を行う場合、その複数のレイヤ画像は、画質の低いベースレイヤ画像と、画質の高い拡張レイヤ画像からなる。なお、画像符号化装置１１及び画像復号装置３１は、単一レイヤ画像でも良いし、ビュースケーラブル符号化、空間スケーラブル符号化、ＳＮＲスケーラブル符号化を任意に組み合わせて行っても良い。ネットワーク２１は、画像符号化装置１１が生成した符号化ストリームＴｅを画像復号装置３１に伝送する。ネットワーク２１は、インターネット（internet）、広域ネットワーク（WAN:Wide Area Network）、小規模ネットワーク（LAN:Local Area Network）又はこれらの組み合わせである。ネットワーク２１は、必ずしも双方向の通信網に限らず、地上波ディジタル放送、衛星放送等の放送波を伝送する一方向又は双方向の通信網であっても良い。また、ネットワーク２１は、ＤＶＤ（Digital Versatile Disc）、ＢＤ（Blue-ray Disc）等の符号化ストリームＴｅを記録した記憶媒体で代替されても良い。
【００１７】
画像復号装置３１は、ネットワーク２１が伝送した符号化ストリームＴｅのそれぞれを復号し、それぞれ復号した複数の復号レイヤ画像Ｔｄ（復号視点画像Ｔｄ）を生成する。
【００１８】
画像表示装置４１は、画像復号装置３１が生成した複数の復号レイヤ画像Ｔｄの全部又は一部を表示する。例えば、ビュースケーラブル符号化においては、全部の場合、３次元画像（立体画像）や自由視点画像が表示され、一部の場合、２次元画像が表示される。画像表示装置４１は、例えば、液晶ディスプレイ、有機ＥＬ（Electro-luminescence）ディスプレイ等の表示デバイスを備える。また、空間スケーラブル符号化、ＳＮＲスケーラブル符号化では、画像復号装置３１、画像表示装置４１が高い処理能力を有する場合には、画質の高い拡張レイヤ画像を表示し、より低い処理能力しか有しない場合には、拡張レイヤほど高い処理能力、表示能力を必要としないベースレイヤ画像を表示する。
【００１９】
＜符号化ストリームＴｅの構造＞
本実施形態に係る画像符号化装置１１および画像復号装置３１の詳細な説明に先立って、画像符号化装置１１によって生成され、画像復号装置３１によって復号される符号化ストリームＴｅのデータ構造について説明する。
【００２０】
図２は、符号化ストリームＴｅにおけるデータの階層構造を示す図である。符号化ストリームＴｅは、例示的に、シーケンス、およびシーケンスを構成する複数のピクチャを含む。図２の（ａ）〜（ｆ）は、それぞれ、シーケンスＳＥＱを既定するシーケンスレイヤ、ピクチャＰＩＣＴを規定するピクチャレイヤ、スライスＳを規定するスライスレイヤ、スライスデータを規定するスライスデータレイヤ、スライスデータに含まれる符号化ツリーユニットを規定する符号化ツリーレイヤ、符号化ツリーに含まれる符号化単位（Ｃｏｄｉｎｇ  Ｕｎｉｔ；ＣＵ）を規定する符号化ユニットレイヤを示す図である。
【００２１】
（シーケンスレイヤ）
シーケンスレイヤでは、処理対象のシーケンスＳＥＱ（以下、対象シーケンスとも称する）を復号するために画像復号装置３１が参照するデータの集合が規定されている。シーケンスＳＥＱは、図２の（ａ）に示すように、ビデオパラメータセット（Video Parameter Set）シーケンスパラメータセットＳＰＳ（Sequence Parameter Set）、ピクチャパラメータセットＰＰＳ（Picture Parameter Set）、ピクチャＰＩＣＴ、及び、付加拡張情報ＳＥＩ（Supplemental Enhancement Information）を含んでいる。ここで＃の後に示される値はレイヤＩＤを示す。図２では、＃０と＃１すなわちレイヤ０とレイヤ１の符号化データが存在する例を示すが、レイヤの種類およびレイヤの数はこれによらない。
【００２２】
ビデオパラメータセットＶＰＳは、複数のレイヤから構成されている動画像において、複数の動画像に共通する符号化パラメータの集合および動画像に含まれる複数のレイヤおよび個々のレイヤに関連する符号化パラメータの集合が規定されている。
【００２３】
シーケンスパラメータセットＳＰＳでは、対象シーケンスを復号するために画像復号装置３１が参照する符号化パラメータの集合が規定されている。例えば、ピクチャの幅や高さが規定される。
【００２４】
ピクチャパラメータセットＰＰＳでは、対象シーケンス内の各ピクチャを復号するために画像復号装置３１が参照する符号化パラメータの集合が規定されている。例えば、ピクチャの復号に用いられる量子化幅の基準値（pic_init_qp_minus26）や重み付き予測の適用を示すフラグ（weighted_pred_flag）が含まれる。なお、ＰＰＳは複数存在してもよい。その場合、対象シーケンス内の各ピクチャから複数のＰＰＳの何れかを選択する。
【００２５】
（ピクチャレイヤ）
ピクチャレイヤでは、処理対象のピクチャＰＩＣＴ（以下、対象ピクチャとも称する）を復号するために画像復号装置３１が参照するデータの集合が規定されている。ピクチャＰＩＣＴは、図２の（ｂ）に示すように、スライスＳ０〜ＳＮＳ−１を含んでいる（ＮＳはピクチャＰＩＣＴに含まれるスライスの総数）。
【００２６】
なお、以下、スライスＳ０〜ＳＮＳ−１のそれぞれを区別する必要が無い場合、符号の添え字を省略して記述することがある。また、以下に説明する符号化ストリームＴｅに含まれるデータであって、添え字を付している他のデータについても同様である。
【００２７】
（スライスレイヤ）
スライスレイヤでは、処理対象のスライスＳ（対象スライスとも称する）を復号するために画像復号装置３１が参照するデータの集合が規定されている。スライスＳは、図２の（ｃ）に示すように、スライスヘッダＳＨ、および、スライスデータＳＤＡＴＡを含んでいる。
【００２８】
スライスヘッダＳＨには、対象スライスの復号方法を決定するために画像復号装置３１が参照する符号化パラメータ群が含まれる。スライスタイプを指定するスライスタイプ指定情報（ｓｌｉｃｅ＿ｔｙｐｅ）は、スライスヘッダＳＨに含まれる符号化パラメータの一例である。
【００２９】
スライスタイプ指定情報により指定可能なスライスタイプとしては、（１）符号化の際にイントラ予測のみを用いるＩスライス、（２）符号化の際に単方向予測、または、イントラ予測を用いるＰスライス、（３）符号化の際に単方向予測、双方向予測、または、イントラ予測を用いるＢスライスなどが挙げられる。
【００３０】
なお、スライスヘッダＳＨには、上記シーケンスレイヤに含まれる、ピクチャパラメータセットＰＰＳへの参照（pic_parameter_set_id）を含んでいても良い。
【００３１】
（スライスデータレイヤ）
スライスデータレイヤでは、処理対象のスライスデータＳＤＡＴＡを復号するために画像復号装置３１が参照するデータの集合が規定されている。スライスデータＳＤＡＴＡは、図２の（ｄ）に示すように、符号化ツリーブロック（CTB:Coded Tree Block）を含んでいる。ＣＴＢは、スライスを構成する固定サイズ（例えば６４×６４）のブロックであり、最大符号化単位（LCU:Largest Cording Unit）と呼ぶこともある。
【００３２】
（符号化ツリーレイヤ）
符号化ツリーレイヤは、図２の（ｅ）に示すように、処理対象の符号化ツリーブロックを復号するために画像復号装置３１が参照するデータの集合が規定されている。符号化ツリーユニットは、再帰的な４分木分割により分割される。再帰的な４分木分割により得られる木構造のノードのことを符号化ツリー（coding tree）と称する。４分木の中間ノードは、符号化ツリーユニット（CTU:Coded Tree Unit）であり、符号化ツリーブロック自身も最上位のＣＴＵとして規定される。ＣＴＵは、分割フラグ（splif_flag）を含み、splif_flagが１の場合には、４つの符号化ツリーユニットＣＴＵに分割される。splif_flagが０の場合には、符号化ツリーユニットＣＴＵは４つの符号化ユニット（CU:Coded Unit）に分割される。符号化ユニットＣＵは符号化ツリーレイヤの末端ノードであり、このレイヤではこれ以上分割されない。符号化ユニットＣＵは、符号化処理の基本的な単位となる。
【００３３】
また、符号化ツリーブロックＣＴＢのサイズ６４×６４画素の場合には、符号化ユニットのサイズは、６４×６４画素、３２×３２画素、１６×１６画素、および、８×８画素の何れかをとり得る。
【００３４】
（符号化ユニットレイヤ）
符号化ユニットレイヤは、図２の（ｆ）に示すように、処理対象の符号化ユニットを復号するために画像復号装置３１が参照するデータの集合が規定されている。具体的には、符号化ユニットは、ＣＵヘッダＣＵＨ、予測ツリー、変換ツリー、ＣＵヘッダＣＵＦから構成される。ＣＵヘッダＣＵＨでは、符号化ユニットが、イントラ予測を用いるユニットであるか、インター予測を用いるユニットであるかなどが規定される。符号化ユニットは、予測ツリー（prediction tree;PT）および変換ツリー（transform tree;TT）のルートとなる。ＣＵヘッダＣＵＦは、予測ツリーと変換ツリーの間、もしくは、変換ツリーの後に含まれる。
【００３５】
予測ツリーは、符号化ユニットが１または複数の予測ブロックに分割され、各予測ブロックの位置とサイズとが規定される。別の表現でいえば、予測ブロックは、符号化ユニットを構成する１または複数の重複しない領域である。また、予測ツリーは、上述の分割により得られた１または複数の予測ブロックを含む。
【００３６】
予測処理は、この予測ブロックごとに行われる。以下、予測の単位である予測ブロックのことを、予測単位（prediction unit;PU、予測ユニット）とも称する。
【００３７】
予測ツリーにおける分割の種類は、大まかにいえば、イントラ予測の場合と、インター予測の場合との２つがある。イントラ予測とは、同一ピクチャ内の予測であり、インター予測とは、互いに異なるピクチャ間（例えば、表示時刻間、レイヤ画像間）で行われる予測処理を指す。
【００３８】
イントラ予測の場合、分割方法は、２Ｎ×２Ｎ（符号化ユニットと同一サイズ）と、Ｎ×Ｎとがある。
【００３９】
また、インター予測の場合、分割方法は、符号化データのpart_modeにより符号化され、２Ｎ×２Ｎ（符号化ユニットと同一サイズ）、２Ｎ×Ｎ、２Ｎ×ｎＵ、２Ｎ×ｎＤ、Ｎ×２Ｎ、ｎＬ×２Ｎ、ｎＲ×２Ｎ、および、Ｎ×Ｎなどがある。なお、２Ｎ×ｎＵは、２Ｎ×２Ｎの符号化ユニットを上から順に２Ｎ×０．５Ｎと２Ｎ×１．５Ｎの２領域に分割することを示す。２Ｎ×ｎＤは、２Ｎ×２Ｎの符号化ユニットを上から順に２Ｎ×１．５Ｎと２Ｎ×０．５Ｎの２領域に分割することを示す。ｎＬ×２Ｎは、２Ｎ×２Ｎの符号化ユニットを左から順に０．５Ｎ×２Ｎと１．５Ｎ×２Ｎの２領域に分割することを示す。ｎＲ×２Ｎは、２Ｎ×２Ｎの符号化ユニットを左から順に１．５Ｎ×２Ｎと０．５Ｎ×１．５Ｎの２領域に分割することを示す。分割数は１、２、４のいずれかであるため、ＣＵに含まれるＰＵは１個から４個である。これらのＰＵを順にＰＵ０、ＰＵ１、ＰＵ２、ＰＵ３と表現する。
【００４０】
また、変換ツリーにおいては、符号化ユニットが１または複数の変換ブロックに分割され、各変換ブロックの位置とサイズとが規定される。別の表現でいえば、変換ブロックは、符号化ユニットを構成する１または複数の重複しない領域のことである。また、変換ツリーは、上述の分割より得られた１または複数の変換ブロックを含む。
【００４１】
変換ツリーにおける分割には、符号化ユニットと同一のサイズの領域を変換ブロックとして割り付けるものと、上述したツリーブロックの分割と同様、再帰的な４分木分割によるものがある。
【００４２】
変換処理は、この変換ブロックごとに行われる。以下、変換の単位である変換ブロックのことを、変換単位（transform unit;TU）とも称する。
【００４３】
（予測パラメータ）
予測ユニットの予測画像は、予測ユニットに付随する予測パラメータによって導出される。予測パラメータには、イントラ予測の予測パラメータもしくはインター予測の予測パラメータがある。以下、インター予測の予測パラメータ（インター予測パラメータ）について説明する。インター予測パラメータは、予測リスト利用フラグpredFlagL0、predFlagL1と、参照ピクチャインデックスrefIdxL0、refIdxL1と、ベクトルmvL0、mvL1から構成される。予測リスト利用フラグpredFlagL0、predFlagL1は、各々Ｌ０リスト、Ｌ１リストと呼ばれる参照ピクチャリストが用いられるか否かを示すフラグであり、値が１の場合に対応する参照ピクチャリストが用いられる。なお、本明細書中「ＸＸであるか否かを示すフラグ」と記す場合、１をＸＸである場合、０をＸＸではない場合とし、論理否定、論理積などでは１を真、０を偽と扱う（以下同様）。但し、実際の装置や方法では真値、偽値として他の値を用いることもできる。２つの参照ピクチャリストが用いられる場合、つまり、predFlagL0＝１，  predFlagL1＝１の場合が、双予測に対応し、１つの参照ピクチャリストを用いる場合、すなわち(predFlagL0, predFlagL1) = (1, 0)もしくは(predFlagL0, predFlagL1) = (0, 1)の場合が単予測に対応する。なお、予測リスト利用フラグの情報は、後述のインター予測フラグinter_pred_idcで表現することもできる。通常、後述の予測画像生成部、予測パラメータメモリでは、予測リスト利用フラグが用いれ、符号化データから、どの参照ピクチャリストが用いられるか否かの情報を復号する場合にはインター予測フラグinter_pred_idcが用いられる。
【００４４】
符号化データに含まれるインター予測パラメータを導出するためのシンタックス要素には、例えば、分割モードpart_mode、マージフラグmerge_flag、マージインデックスmerge_idx、インター予測フラグinter_pred_idc、参照ピクチャインデックスrefIdxLX、予測ベクトルインデックスmvp_LX_idx、差分ベクトルmvdLXがある。
【００４５】
（参照ピクチャリストの一例）
次に、参照ピクチャリストの一例について説明する。参照ピクチャリストとは、参照ピクチャメモリ３０６（図５）に記憶された参照ピクチャからなる列である。図３は、参照ピクチャリストの一例を示す概念図である。参照ピクチャリスト６０１において、左右に一列に配列された５個の長方形は、それぞれ参照ピクチャを示す。左端から右へ順に示されている符号、Ｐ１、Ｐ２、Ｑ０、Ｐ３、Ｐ４は、それぞれの参照ピクチャを示す符号である。Ｐ１等のＰとは、視点Ｐを示し、そしてＱ０のＱとは、視点Ｐとは異なる視点Ｑを示す。Ｐ及びＱの添字は、ピクチャ順序番号ＰＯＣを示す。refIdxLXの真下の下向きの矢印は、参照ピクチャインデックスrefIdxLXが、参照ピクチャメモリ３０６において参照ピクチャＱ０を参照するインデックスであることを示す。
【００４６】
（参照ピクチャの例）
次に、ベクトルを導出する際に用いる参照ピクチャの例について説明する。図４は、参照ピクチャの例を示す概念図である。図４において、横軸は表示時刻を示し、縦軸は視点を示す。図４に示されている、縦２行、横３列（計６個）の長方形は、それぞれピクチャを示す。６個の長方形のうち、下行の左から２列目の長方形は復号対象のピクチャ（対象ピクチャ）を示し、残りの５個の長方形がそれぞれ参照ピクチャを示す。対象ピクチャから上向きの矢印で示される参照ピクチャＱ０は対象ピクチャと同表示時刻であって視点が異なるピクチャである。対象ピクチャを基準とする変位予測においては、参照ピクチャＱ０が用いられる。対象ピクチャから左向きの矢印で示される参照ピクチャＰ１は、対象ピクチャと同じ視点であって、過去のピクチャである。対象ピクチャから右向きの矢印で示される参照ピクチャＰ２は、対象ピクチャと同じ視点であって、未来のピクチャである。対象ピクチャを基準とする動き予測においては、参照ピクチャＰ１又はＰ２が用いられる。
【００４７】
（インター予測フラグと予測リスト利用フラグ）
インター予測フラグと、予測リスト利用フラグpredFlagL0、predFlagL1の関係は以下のように相互に変換可能である。そのため、インター予測パラメータとしては、予測リスト利用フラグを用いても良いし、インター予測フラグを用いてもよい。また、以下、予測リスト利用フラグを用いた判定は、インター予測フラグに置き替えても可能である。逆に、インター予測フラグを用いた判定は、予測リスト利用フラグに置き替えても可能である。
【００４８】
インター予測フラグ  ＝  （predFlagL1＜＜１）  ＋  predFlagL0
predFlagL0  ＝インター予測フラグ  ＆  １
predFlagL1  ＝インター予測フラグ  ＞＞  １
ここで、＞＞は右シフト、＜＜は左シフトである。
【００４９】
（マージ予測とＡＭＶＰ予測）
予測パラメータの復号（符号化）方法には、マージ予測（ｍｅｒｇｅ）モードとＡＭＶＰ（Adaptive Motion Vector Prediction、適応動きベクトル予測）モードがある、マージフラグmerge_flagは、これらを識別するためのフラグである。マージ予測モードでも、ＡＭＶＰモードでも、既に処理済みのブロックの予測パラメータを用いて、対象ＰＵの予測パラメータが導出される。マージ予測モードは、予測リスト利用フラグpredFlagLX（インター予測フラグinter_pred_idc）、参照ピクチャインデックスrefIdxLX、ベクトルmvLXを符号化データに含めずに、既に導出した予測パラメータをそのまま用いるモードであり、ＡＭＶＰモードは、インター予測フラグinter_pred_idc、参照ピクチャインデックスrefIdxLX、ベクトルmvLXを符号化データに含めるモードである。なおベクトルmvLXは、予測ベクトルを示す予測ベクトルインデックスmvp_LX_idxと差分ベクトル（mvdLX）として符号化される。
【００５０】
インター予測フラグinter_pred_idcは、参照ピクチャの種類および数を示すデータであり、Ｐｒｅｄ＿Ｌ０、Ｐｒｅｄ＿Ｌ１、Ｐｒｅｄ＿Ｂｉの何れかの値をとる。Ｐｒｅｄ＿Ｌ０、Ｐｒｅｄ＿Ｌ１は、各々Ｌ０リスト、Ｌ１リストと呼ばれる参照ピクチャリストに記憶された参照ピクチャが用いられることを示し、共に１枚の参照ピクチャを用いること（単予測）を示す。Ｌ０リスト、Ｌ１リストを用いた予測を各々Ｌ０予測、Ｌ１予測と呼ぶ。Ｐｒｅｄ＿Ｂｉは２枚の参照ピクチャを用いること（双予測）を示し、Ｌ０リストとＬ１リストに記憶された参照ピクチャの２つを用いることを示す。予測ベクトルインデックスmvp_LX_idxは予測ベクトルを示すインデックスであり、参照ピクチャインデックスrefIdxLXは、参照ピクチャリストに記憶された参照ピクチャを示すインデックスである。なお、ＬＸは、Ｌ０予測とＬ１予測を区別しない場合に用いられる記述方法であり、ＬＸをＬ０、Ｌ１に置き換えることでＬ０リストに対するパラメータとＬ１リストに対するパラメータを区別するする。例えば、refIdxL0はＬ０予測に用いる参照ピクチャインデックス、refIdxL1はＬ１予測に用いる参照ピクチャインデックス、refIdx（refIdxLX）は、refIdxL0とrefIdxL1を区別しない場合に用いられる表記である。
【００５１】
マージインデックスmerge_idxは、処理が完了したブロックから導出される予測パラメータ候補（マージ候補）のうち、いずれかの予測パラメータを復号対象ブロックの予測パラメータとして用いるかを示すインデックスである。
【００５２】
（動きベクトルと変位ベクトル）
ベクトルmvLXには、動きベクトルと変位ベクトル（disparity vector、視差ベクトル）がある。動きベクトルとは、あるレイヤのある表示時刻でのピクチャにおけるブロックの位置と、異なる表示時刻（例えば、隣接する離散時刻）における同一のレイヤのピクチャにおける対応するブロックの位置との間の位置のずれを示すベクトルである。変位ベクトルとは、あるレイヤのある表示時刻でのピクチャにおけるブロックの位置と、同一の表示時刻における異なるレイヤのピクチャにおける対応するブロックの位置との間の位置のずれを示すベクトルである。異なるレイヤのピクチャとしては、異なる視点のピクチャである場合、もしくは、異なる解像度のピクチャである場合などがある。特に、異なる視点のピクチャに対応する変位ベクトルを視差ベクトルと呼ぶ。以下の説明では、動きベクトルと変位ベクトルを区別しない場合には、単にベクトルmvLXと呼ぶ。ベクトルmvLXに関する予測ベクトル、差分ベクトルを、それぞれ予測ベクトルmvpLX、差分ベクトルmvdLXと呼ぶ。ベクトルmvLXおよび差分ベクトルmvdLXが、動きベクトルであるか、変位ベクトルであるかは、ベクトルに付随する参照ピクチャインデックスrefIdxLXを用いて行われる。
【００５３】
（画像復号装置の構成）
次に、本実施形態に係る画像復号装置３１の構成について説明する。図５は、本実施形態に係る画像復号装置３１の構成を示す概略図である。画像復号装置３１は、エントロピー復号部３０１、予測パラメータ復号部３０２、参照ピクチャメモリ（参照画像記憶部、フレームメモリ）３０６、予測パラメータメモリ（予測パラメータ記憶部、フレームメモリ）３０７、予測画像生成部３０８、逆量子化・逆ＤＣＴ部３１１、及び加算部３１２、残差格納部３１３（残差記録部）を含んで構成される。
【００５４】
また、予測パラメータ復号部３０２は、インター予測パラメータ復号部３０３及びイントラ予測パラメータ復号部３０４を含んで構成される。予測画像生成部３０８は、インター予測画像生成部３０９及びイントラ予測画像生成部３１０を含んで構成される。
【００５５】
エントロピー復号部３０１は、外部から入力された符号化ストリームＴｅに対してエントロピー復号を行って、個々の符号（シンタックス要素）を分離し復号する。分離された符号には、予測画像を生成するための予測情報および、差分画像を生成するための残差情報などがある。
【００５６】
エントロピー復号部３０１は、分離した符号の一部を予測パラメータ復号部３０２に出力する。分離した符号の一部とは、例えば、予測モードPredMode、分割モードpart_mode、マージフラグmerge_flag、マージインデックスmerge_idx、インター予測フラグinter_pred_idc、参照ピクチャインデックスrefIdxLX、予測ベクトルインデックmvp_LX_idx、差分ベクトルmvdLXである。どの符号を復号するか否かの制御は、予測パラメータ復号部３０２の指示に基づいて行われる。エントロピー復号部３０１は、量子化係数を逆量子化・逆ＤＣＴ部３１１に出力する。この量子化係数は、符号化処理において、残差信号に対してＤＣＴ（Discrete Cosine Transform、離散コサイン変換）を行い量子化して得られる係数である。
【００５７】
インター予測パラメータ復号部３０３は、エントロピー復号部３０１から入力された符号に基づいて、予測パラメータメモリ３０７に記憶された予測パラメータを参照してインター予測パラメータを復号する。
【００５８】
インター予測パラメータ復号部３０３は、復号したインター予測パラメータを予測画像生成部３０８に出力し、また予測パラメータメモリ３０７に記憶する。インター予測パラメータ復号部３０３の詳細については後述する。
【００５９】
イントラ予測パラメータ復号部３０４は、エントロピー復号部３０１から入力された符号に基づいて、予測パラメータメモリ３０７に記憶された予測パラメータを参照してイントラ予測パラメータを復号する。イントラ予測パラメータとは、ピクチャブロックを１つのピクチャ内で予測する処理で用いるパラメータ、例えば、イントラ予測モードＩｎｔｒａPredModeである。イントラ予測パラメータ復号部３０４は、復号したイントラ予測パラメータを予測画像生成部３０８に出力し、また予測パラメータメモリ３０７に記憶する。
【００６０】
イントラ予測パラメータ復号部３０４は、輝度と色差で異なるイントラ予測モードを導出しても良い。この場合、イントラ予測パラメータ復号部３０４は、輝度の予測パラメータとして輝度予測モードIntraPredModeY、色差の予測パラメータとして、色差予測モードIntraPredModeCを復号する。輝度予測モードIntraPredModeYは、３５モードであり、プレーナ予測（０）、ＤＣ予測（１）、方向予測（２〜３４）が対応する。色差予測モードIntraPredModeCは、プレーナ予測（０）、ＤＣ予測（１）、方向予測（２〜３４）、ＬＭモード（３５）の何れかを用いるもの。イントラ予測パラメータ復号部３０４は、IntraPredModeCは輝度モードと同じモードであるか否かを示すフラグを復号し、フラグが輝度モードと同じモードであることを示せば、IntraPredModeCにIntraPredModeYを割り当て、フラグが輝度モードと異なるモードであることを示せば、I IntraPredModeCとして、プレーナ予測（０）、ＤＣ予測（１）、方向予測（２〜３４）、ＬＭモード（３５）を復号しても良い。
【００６１】
参照ピクチャメモリ３０６は、加算部３１２が生成した参照ピクチャのブロック（参照ピクチャブロック）を、復号対象のピクチャ及びブロック毎に予め定めた位置に記憶する。
【００６２】
予測パラメータメモリ３０７は、予測パラメータを、復号対象のピクチャ及びブロック毎に予め定めた位置に記憶する。具体的には、予測パラメータメモリ３０７は、インター予測パラメータ復号部３０３が復号したインター予測パラメータ、イントラ予測パラメータ復号部３０４が復号したイントラ予測パラメータ及びエントロピー復号部３０１が分離した予測モードpredModeを記憶する。記憶されるインター予測パラメータには、例えば、予測リスト利用フラグpredFlagLX（インター予測フラグinter_pred_idc）、参照ピクチャインデックスrefIdxLX、ベクトルmvLXがある。
【００６３】
予測画像生成部３０８には、エントロピー復号部３０１から入力された予測モードpredModeが入力され、また予測パラメータ復号部３０２から予測パラメータが入力される。また、予測画像生成部３０８は、参照ピクチャメモリ３０６から参照ピクチャを読み出す。予測画像生成部３０８は、予測モードpredModeが示す予測モードで、入力された予測パラメータと読み出した参照ピクチャを用いて予測ピクチャブロックＰ（予測画像）を生成する。
【００６４】
ここで、予測モードpredModeがインター予測モードを示す場合、インター予測画像生成部３０９は、インター予測パラメータ復号部３０３から入力されたインター予測パラメータと読み出した参照ピクチャを用いてインター予測により予測ピクチャブロックＰを生成する。予測ピクチャブロックＰは予測単位ＰＵに対応する。ＰＵは、上述したように予測処理を行う単位となる複数の画素からなるピクチャの一部分、つまり１度に予測処理が行われる復号対象ブロックに相当する。
【００６５】
インター予測画像生成部３０９は、予測リスト利用フラグpredFlagLXが１である参照ピクチャリスト（Ｌ０リスト、もしくはＬ１リスト）に対し、参照ピクチャインデックスrefIdxLXで示される参照ピクチャから、復号対象ブロックを基準としてベクトルmvLXが示す位置にある参照ピクチャブロックを参照ピクチャメモリ３０６から読み出す。インター予測画像生成部３０９は、読み出した参照ピクチャブロックについて予測を行って予測ピクチャブロックＰを生成する。インター予測画像生成部３０９は、生成した予測ピクチャブロックＰを加算部３１２に出力する。
【００６６】
予測モードpredModeがイントラ予測モードを示す場合、イントラ予測画像生成部３１０は、イントラ予測パラメータ復号部３０４から入力されたイントラ予測パラメータと読み出した参照ピクチャを用いてイントラ予測を行う。具体的には、イントラ予測画像生成部３１０は、復号対象のピクチャであって、既に復号されたブロックのうち復号対象ブロックから予め定めた範囲にある参照ピクチャブロックを参照ピクチャメモリ３０６から読み出す。予め定めた範囲とは、復号対象ブロックがいわゆるラスタースキャンの順序で順次移動する場合、例えば、左、左上、上、右上の隣接ブロックのうちのいずれかであり、イントラ予測モードによって異なる。ラスタースキャンの順序とは、各ピクチャにおいて、上端から下端まで各行について、順次左端から右端まで移動させる順序である。
【００６７】
イントラ予測画像生成部３１０は、読み出した参照ピクチャブロックについてイントラ予測モードＩｎｔｒａPredModeが示す予測モードで予測を行って予測ピクチャブロックを生成する。イントラ予測画像生成部３１０は、生成した予測ピクチャブロックＰを加算部３１２に出力する。
【００６８】
イントラ予測パラメータ復号部３０４において、輝度と色差で異なるイントラ予測モードを導出する場合、イントラ予測画像生成部３１０は、輝度予測モードIntraPredModeYに応じて、プレーナ予測（０）、ＤＣ予測（１）、方向予測（２〜３４）の何れかによって輝度の予測ピクチャブロックを生成し、色差予測モードIntraPredModeCに応じて、プレーナ予測（０）、ＤＣ予測（１）、方向予測（２〜３４４）、ＬＭモード（３５）の何れかによって色差の予測ピクチャブロックを生成する。ここでは、正則化コストを用いたＬＭモードを説明する。正則化コストとは、最小２乗法による予測パラメータ導出において、目的関数にパラメータコストとして加算される項である。ＬＭモードは、対象ブロックに隣接する処理済み画像を用いて、輝度の画素値から色差の画素値を予測する線形予測パラメータを導出し、その線形予測パラメータに基づいて、処理済みの輝度のブロックから、色差のピクチャブロックを生成する（ＴＹＰＥ０）と、対象ブロックに隣接する処理済み画像を用いて、色差の画素値から別の色差の画素値を予測する線形予測パラメータを導出し、その線形予測パラメータに基づいて、処理済みの色差のブロックから、別の色差のピクチャブロックを生成する（ＴＹＰＥ１）。このようなある色コンポーネントから別の色コンポーネントを推定する線形予測パラメータを導出し、導出した線形予測パラメータを用いて予測画像を生成する予測はＬＭ予測と呼ばれる。予測の方法は、後述の照度補償部３０９３と概略は同じである。イントラ予測画像生成部３１０は、図示しないＤＣ予測部３１０１、プレーナ予測部３１０２、方向予測部３１０３、ＬＭ予測部３１０４から構成される。なお、ＬＭ予測のＬＭは、Liner Modelの略である。
【００６９】
（最小二乗法と正則化項）
入力参照信号xiと目的参照信号yiからなる２次元の数列（xi, yi）、i=0..N-1（Nは数列の個数）に対して、目的参照信号yiを入力参照信号xiから予測する以下の線形予測式を考える。なお、yiは教師信号とも呼ばれる。
【００７０】
yi = a * xi + b
このとき、線形予測の２乗誤差の和Eは以下の式で定義される。
【００７１】
E= |yi - a * xi - b| ^ 2                      式（Ｅ−１）
上記２乗誤差Ｅを最小化する線形予測パラメータa, bは、Eをa、bで偏微分して得られる式を０と置くことから以下のように導出される。
【００７２】
a = (Σxy - N*ΣxΣy) / (Σxx - NΣxΣx)       式（Ａ−１）
b = (Σy/N - a*Σx/N)
しかしながら、上記最小２乗法によるパラメータ導出では、参照信号xi, yiに雑音成分等が加わっていると、十分な数Ｎの参照信号を用意しない場合には、精度の高い線形予測パラメータa、bを推定できない場合がある。また、Ｎが小さい場合や、参照信号に含まれる入力参照画素xiと目的参照画素yiの変化が小さい場合に、オーバーフィッティングになる可能性がある。そのため線形予測パラメータを推定する場合に、誤差Ｅに適当なパラメータコスト（正則化コスト）を加えたコストを考慮して、線形予測パラメータを導出するような方式ことが知られている。一般には、正則化項として、パラメータの値を０に近くなるような正則化コスト、すなわち、以下の式のようにaとbの2乗項に重みlambda_a、lambda_bをかけた形が知られている（リッジ回帰）。
【００７３】
E= |yi - a * xi - b| ^ 2 + lambda_a * a^2 + lambda_b * b^2      式（Ｅ−２）
しかしながら、リッジ回帰の正則化項を画像に適用しても、画素間の線形回帰では、線形予測パラメータaとbが０に近い値である、という仮定自体が成り立たないため、適当な結果が得られない。
【００７４】
発明者は、画像に適用する場合には、パラメータaとパラメータbの分布や特性を調査した上で、以下のように線形パラメータaを適当な値kに近づけるような正則化項（パラメータaとkとのずれの2乗に比例する項）を用いることが適当であることを見出した。
【００７５】
E´= |yi - a * xi - b| ^ 2 + lambda_a * (a - k)^2      式（Ｅ−３）
２乗誤差の和E´をa、bで偏微分して得られる式を０と置くことで、上記式を、解くことができる。これにより、以下の線形パラメータａ、ｂの推定式が得られる。
【００７６】
a = (Σxy - N*ΣxΣy - k* lambda_a) / (Σxx - NΣxΣx - lambda_a)       式（Ａ−３）
b = (Σy/N - a*Σx/N) 
以下、線形予測パラメータaを導出する際の除算の式において、分子に相当する値を第１パラメータa1（numerator）、分母に相当する数を第２パラメータa2（denominator）と呼ぶ。ここで、Σxy - N*ΣxΣy = a1、Σxx - NΣxΣx = a2、lambda_a＝acostと置くと、最小２乗法の式（Ａ−１）は、以下のように変形できる。
【００７７】
a = a1 / a2
さらに、正則化項を含む式（Ａ−３）は、以下のように変形される。
【００７８】
a = (a1 + k*acost) / (a2 + acost)
上記は、第１パラメータa1の第２パラメータa2による除算で示される線形予測パラメータａの導出において、第１パラメータa1をk*costを加算して更新し、第２パラメータa2をacostを加算して更新することで、正則化項を考慮した線形予測パラメータａが得られることを示す。以下、本特許ではacost、及び、k*acostのような第１パラメータa1、第２パラメータa2の更新に用いられる値を正則化コスト（もしくはパラメータコスト）と呼ぶ。
【００７９】
（ＬＭ予測部３１０４）
図１は、ＬＭ予測部３１０４の構成を示すブロック図である。ＬＭ予測部３１０４は、線形予測パラメータを導出するＬＭパラメータ推定部３１０４１と、導出した線形予測パラメータを用いて実際に予測を行うＬＭ予測フィルタ部３１０４２から構成される。ＬＭパラメータ推定部３１０４１は、ＬＭ積算値導出部３１０４１２、ＬＭ加算値導出部３１０４１３、ＬＭ第１パラメータ導出部３１０４１４、ＬＭ第２パラメータ導出部３１０４１５、ＬＭパラメータａ導出部３１０４１６、ＬＭパラメータｂ導出部３１０４１７から構成される。
【００８０】
ＬＭパラメータ推定部３１０４１は、入力参照画素xiから目的参照画素yiを、yi = a * xi + bの式で予測する場合の線形予測パラメータaとbを導出する。線形予測パラメータａは傾き、線形予測パラメータｂはオフセットに相当する。
【００８１】
図１７は、ＬＭパラメータ推定部３１０４１が推定するパラメータａの分布を示すヒストグラムである。図１７（ａ）と（ｂ）は別々のピクチャから得られた分布である。パラメータａの分布はピクチャ（シーケンス）によって異なるが、ａの絶対値が０から１の範囲に比較的多く分布する。輝度と色差間の色コンポーネント予測に対応するＴＹＰＥ＝０ではパラメータａは正と負の両方に存在し、色差と別の色差間の色コンポーネント予測に対応するＴＹＹＰＥ＝１では、比較的、パラメータａが負となることが多いが、パラメータａが正となる場合も発生する。
【００８２】
ＬＭパラメータ推定部３１０４１は、図１２に示す対象ブロックの周囲のある色コンポーネントの参照ブロックの周囲に存在する画素を入力参照画素x[]、対象ブロックの周囲に存在する別の色コンポーネントの画素を目的参照画素y[]とし、これら入力参照画素x[]と目的参照画素y[]に基づいて、目的参照画素y[]を入力参照画素xから線形予測する場合のパラメータであるパラメータａ、パラメータｂを導出する。
【００８３】
ＬＭ加算値導出部３１０４１３は、目的参照画素yの和Yと入力参照画素xの和Xを、以下の式（Ｂ−２）式（Ｂ−３）により導出する。
【００８４】
ＬＭ積算値導出部３１０４１２は目的参照画素yと入力参照画素xの積の和XYと、目的参照画素xの２乗の和XXを、以下の式（Ｂ−４）〜（Ｂ−５）により導出する。
【００８５】
X = Σx[i]                                  式（Ｂ−２）
Y = Σy[i]                                  式（Ｂ−３）
XX = Σ(x[i] * x[i])                      式（Ｂ−４）
XY = Σ(y[i] * y[i])                      式（Ｂ−５）
ここで、Σは、参照領域に対する和であり、参照領域の画素を特定するインデックスi
に対する和を導出する。つまり、Σ(z[i])はiについて０から参照信号数N-1までz[i]の和をとる処理を示す。y[i]は、復号画像のインデックスiにおける画素値。x[i]は、参照画像のインデックスiにおける画素値。カウントシフト値iCountShiftは、参照領域のサイズの２の対数である。
【００８６】
iCountShift = log2 (参照領域の画素数)           式（Ｂ−６）
ＬＭ第１パラメータ導出部３１０４１４は、目的参照画素yと入力参照画素xの積の和XYと、目的参照画素yの和Yと入力参照画素xの和Xの積の差に基づいて第１パラメータa1を以下の式により導出する。
【００８７】
a1 = ( XY << iCountShift ) - (Y * X);         式（Ｂ−７）
式（Ｂ−７）に示すように、XYは、カウントシフト値iCountShiftだけ左シフト、YとXの積は積算シフト値precShiftだけ右シフトしてから差を算出する。
【００８８】
ＬＭ第２パラメータ導出部３１０４１５は、入力参照画素xの２乗の和XXと、入力参照画素xの和Xの２乗の差に基づいて第２パラメータa2を以下の式により導出する。
【００８９】
a2 = ( XX << iCountShift ) - (X * X);         式（Ｂ−８）  導出された第１パラメータa１と第２パラメータa2は、ＬＭパラメータａ導出部３１０４１６に出力される。
【００９０】
また、ＬＭ第１パラメータ導出部３１０４１４、ＬＭ第２パラメータ導出部３１０４１５は、演算に必要なビットレンジを低減させるため、以下のように式（Ｂ−７）、式（Ｂ−８）をiCountShiftだけ右シフトした式（Ｂ−７´）、（Ｂ−８´）を用いて導出しても良い。
【００９１】
aveX = X >> iCountShift
aveY = Y >> iCountShift
a1 = XY - ((aveY * aveX)<<iCountShift);         式（Ｂ−７´）
a2 = XX - ((aveX * aveX)<<iCountShift);         式（Ｂ−８´）
このとき、aveX、aveYは、XとYの平均値であり、各々XとYをiCountShiftだけ右シフトして導出する。
【００９２】
また、、ＬＭ第１パラメータ導出部３１０４１４、ＬＭ第２パラメータ導出部３１０４１５は、演算精度を向上させるために、平均値aveX、aveYを導出した際の余りremX、remYも考慮して、以下のようにXYからaveYとaveXの積に基づいた数の差にaveXとremYの積とaveYとremXの積を加えることでa1、XXからaveXとaveXの積に基づいた数の差にaveXとremXの積の２倍を加えることでa2を導出しても良い。
【００９３】
aveX = X >> iCountShift
remX = X - aveX
aveY = Y >> iCountShift
remX = Y - aveY
a1 = XY - ((aveY * aveX)<<iCountShift) + (aveX * remY) + (aveY * remX);         式（Ｂ−７´´）
a2 = XX - ((aveX * aveX)<<iCountShift) + 2*(aveX * remX);     式（Ｂ−８´´）
なお、ＬＭ積算値導出部３１０４１２は、目的参照画素yと入力参照画素xの積の和XYの導出時と、入力参照画素xの２乗の和XXにおいて、下記のように所定の積算シフト値precShiftだけ右シフトしてから加算しても良い。
【００９４】
XX = Σ(x[i] * x[i])>>precShift          式（Ｂ−４´）
XY = Σ(y[i] * y[i])>>precShift          式（Ｂ−５´）
なお、積算シフト値precShiftは以下等から導出しても良い。
【００９５】
precShift = Max(0, bitDepth - 12)        式（Ｂ−１）
なお、積算シフト値precShiftを用いる場合には以下に示すように、第１パラメータa1と第２パラメータa2、XY及びXXは、カウントシフト値iCountShiftだけ左シフト、XとXの積及びXとYの積は積算シフト値precShiftだけ右シフトしてから差を算出しても良い。
【００９６】
a1 = ( XY << iCountShift ) - (Y * X) >> precShift
a2 = ( XX << iCountShift ) - (X * X) >> precShift
一般の小数点精度の最小２乗法であればパラメータaをa1/a2で推定することができるが、本特許では、パラメータの頑健性を向上させるために、場合に応じて正則化コストを加える。また、さらに本特許では、整数演算を行うために、シフト値iShiftを導入し、パラメータaとして、iShiftだけ左シフトした値(a1 << iShift) / a2を導出する。
【００９７】
（本実施形態のＬＭパラメータ導出部）
本実施形態のＬＭパラメータ推定部３１０４１の構成は、インデックスｉ（i=0からN-1）に対応する画素値xiと画素値yiの組を入力として、xiからyiを予測する線形予測パラメータａとｂを導出する線形予測パラメータ導出部と、上記線形予測パラメータに基づいて予測画像を生成する線形予測部を備える画像予測装置であって、上記線形予測パラメータ導出部は、画素値xiと画素値yiの積の和XYと画素値xiの和Xと画素値yiの和Yの積に基づいて導出される第１パラメータa1と、画素値xiと画素値xiの積の和XXと画素値xiの和Xと画素値xiの和Xの積に基づいて導出される第２パラメータa2から線形予測パラメータａを導出する手段を備え、上記線形予測パラメータ導出部は、上記第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN未満もしくは上記閾値THN以下である場合には第１コストa1costNを上記第１パラメータa1から減算する手段と、上記第１パラメータa1と所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THP未満もしくは上記閾値THP以下の場合には第２コストa1costPを上記第１パラメータa1に加算する手段の少なくとも１つの手段を備えることを特徴とする。
【００９８】
図１８は、ＬＭパラメータａ導出部３１０４１６の構成を示すブロック図である。ＬＭパラメータａ導出部３１０４１６は、ＬＭ正則化コスト加算部３１０４１８、ＬＭ第１パラメータ正規化シフト部３１０４１６１、ＬＭ第２パラメータ正規化シフト部３１０４１６２、ＬＭ除算値導出部３１０４１６３から構成される。
【００９９】
ＬＭパラメータａ導出部３１０４１６は、図示しないＬＭ第１パラメータクリップ部により、第１パラメータa1を第２パラメータa2の大きさに応じて制限する。例えば以下の式のように、a1を-2*a2以上、a2の2倍以下にクリップしても良い。
【０１００】
a1 = Clip3(-2*a2, (127*a2)>>6, a1)          式（Ｂ−１２）
上記により、a1 / a2の値は、-2から127/64の間にクリップされる。従って、パラメータａの値であるa1 / a2 << iShiftの値も、-2<<iShiftから(127/64)<<iShiftにクリップされる。すなわち、iShift=6の場合、パラメータａは-128〜127となり、8ビット整数の範囲に入る。
【０１０１】
なお、クリップは上記に限定せず、以下の式などでも良い。
【０１０２】
a1 = Clip3(-m1*a2, m2*a2, a1)          式（Ｂ−１２´）
ここでm1、m2は正の整数。例えばm1=2から3、m2=2から3の値などが適当である。さらに整数よりも細かい単位で制限するために、以下の式でもクリップしても良い。
【０１０３】
a1 = Clip3(-(m1*a2>>m3), (m2*a2>>m43), a1)          式（Ｂ−１２´´）
ここで、m1、m2、m3、m4は正の整数。例えばm1=3、m2=3、m3=1、m4=1とすれば、a1をa2の-1.5倍〜1.5倍の間の値に制限できる。つまりa1/a2を-1.5〜1.5の間に制限できる。
【０１０４】
本実施形態のＬＭパラメータ推定部３１０４１は、ＬＭ正則化コスト加算部３１０４１８により、第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN以下である場合には０以上の第１コストa1costNを第１パラメータa1から減算するＬＭ第１パラメータ正則化コスト加算部と、第１パラメータa1と１つ以上の所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THPより以上の場合には０以上の第２コストa1costPを第１パラメータa1に加算するＬＭ第１パラメータ正則化コスト加算部の少なくとも１つの手段を備える。ここで、第１パラメータa1、第２パラメータa2に正則化コストを減算もしくは加算した更新後の値もa1、a2と表記するが、更新前と変更後の値を区別する場合には、更新前のパラメータをa1、a2、更新後のパラメータをa1´、a2´と表記する。
【０１０５】
ＬＭ正則化コスト加算部３１０４１８は、例えば、第１パラメータa1が所定の閾値THNより小さい場合、もしくは、第１パラメータa1が所定の閾値THN以下の場合に、第１パラメータa1に０以上の第１コストa1costNを減算する。その他、ＬＭ正則化コスト加算部３１０４１８の具体的な構成及びその変形例は後述する。
【０１０６】
a1 = a1 + a1costN
（ＬＭ正規化シフト部）
以下、a1 << iShift / a2に対応するパラメータａを整数処理により導出する前に、導出時に処理に必要な値域が大きくなりすぎないように、ＬＭ第１パラメータ正規化シフト部３１０４１６１で、a1の大きさをシフトにより調節しa1sを導出する。また、同様に、a2で除算する場合に参照するテーブルのサイズを制限するために、ＬＭ第２パラメータ正規化シフト部３１０４１６２によりa2の大きさをシフトにより調節しa2sを導出する。
【０１０７】
ＬＭパラメータ第２正規化シフト部３１０４１６２は、第２パラメータa2の大きさに応じて、図１６のテーブルの導出に用いた所定のビット幅ShiftA2に対して、以下の式により第２正規化シフト値iScaleShiftA2を導出する。導出された第２正規化シフト値iScaleShiftA2は、ＬＭ除算値導出部３１０４１６３に出力される。
【０１０８】
iScaleShiftA2 = Max( 0, Floor( Log2( Abs( a2 ) ) ) - (ShiftA2-1))
式（Ｂ−１４）
なお、Floor( Log2( Abs( x ) ) )は、a2を32ビットレジスタに格納した際に、ビット列の左側 Leftmost bit から見て 0 が連続している数である Number of Leading Zero (NLZ)を用いて、
Floor( Log2( Abs( x ) ) ) = 32 - NLZ(x)
により求めることができる。なお、64ビットレジスタを用いる場合には、64 - NLZ(x)により導出することができる。
【０１０９】
なお、NLZの導出は、比較的複雑な計算が必要になるため、数が少ない方が好ましい。
【０１１０】
ＬＭ第１パラメータ正規化シフト部３１０４１６１は、第２正規化シフト値iScaleShiftA2に応じて、以下の式により第１正規化シフト値iScaleShiftA1を導出する。導出された第１正規化シフト値iScaleShiftA1は、ＬＭ除算値導出部３１０４１６３に出力される。
【０１１１】
iScaleShiftA1 = Max(0, iScaleShiftA2 - offsetA1)        式（Ｂ−１３）
なお、offsetA1はここでは10以下を満たす定数とする。
【０１１２】
上記では、第１正規化シフト値を第２正規化シフト値に基づいて導出するため、第１正規化パラメータを導出する処理が容易になるという効果を奏する。
【０１１３】
なお、ＬＭ第１パラメータ正規化シフト部３１０４１６１は、第２正規化シフト値iScaleShiftA2に応じずに、以下の式で導出しても良い。
【０１１４】
iScaleShiftA1 = Max( 0, Floor( Log2( Abs( a1 ) ) ) - (ShiftA1-1))
ここでShiftA1は、例えば６〜１４などの固定値を用いても良いし、ビットデプスに依存する値（例えばbitDepth - 2）などの値を用いても良い。
【０１１５】
ＬＭ第１パラメータ正規化シフト部３１０４１６１、ＬＭ第２パラメータ正規化シフト部３１０４１６２は、第１パラメータa1を第１正規化シフト値iScaleShiftA1、第２パラメータa2を第２正規化シフト値iScaleShiftA2により右シフトし、正規化第１パラメータa1sと正規化第２パラメータa2sを導出する。
【０１１６】
a1s = a1 >> iScaleShiftA1          式（Ｂ−１５）
a2s = a2 >> iScaleShiftA2          式（Ｂ−１６）
なお、正則化コストを用いる場合には、上記式において、a1、a2の代わりに正則化コストを減算もしくは加算後のa1´、a2´を用いる。
【０１１７】
これにより、正規化第１パラメータa1sと正規化第２パラメータa2sは各々、０から2のiShiftA1乗-1、―2のShiftA2乗から2のShiftA2乗-1の間の値に正規化される。
【０１１８】
ＬＭ除算値導出部３１０４１６３は、a1s / a2sに相当に相当する数を導出する。本実施形態では、整数値として扱うために（1<< ScaleShiftA）倍した数を、テーブルを用いて導出する。
【０１１９】
ＬＭ除算値導出部３１０４１６３は、第１正規化シフト値iScaleShiftA1と第２正規化シフト値iScaleShiftA2の差に基づいて、以下の式によりパラメータａシフト値iScaleShiftAを導出する。
【０１２０】
ScaleShiftA  = ShiftA1 + iScaleShiftA2 - iScaleShiftA1 - iShift          式（Ｂ−１８）
ここで、iScaleShiftA1 = Max(0, iScaleShiftA2 - offsetA1)であるから、以下の式が得られる。
【０１２１】
ScaleShiftA  <= ShiftA1 + iScaleShiftA2 - (iScaleShiftA2 - offsetA1) - iShift
ScaleShiftA  <= ShiftA1 + offsetA1 - iShift
offsetA1は0以上、iShiftは5から8ビット、ShiftA1は14ビット〜15ビットであるから、ScaleShiftAは常に0以上になる。
【０１２２】
ＬＭ除算値導出部３１０４１６３は、正規化第２パラメータa2sに応じて定まる逆数テーブル値invTableを参照し、正規化第１パラメータa1sとの積をとりテーブルシフト値（ScaleShiftA）で右シフトすることにより、以下の式によりパラメータａを導出する。
【０１２３】
a = (a1s * invTable[a2s]) >> (ScaleShiftA)           式（Ｂ−１９）
図１６は、本実施形態で用いる逆数テーブル値invTable[]である。図１６の逆数invTable[x]は、インデックスxが０の場合に０となり、インデックスxが０以外の場合には、所定の定数（２のShiftA2）をxで割って整数化した値から導出される。すなわち、
invTable[x]  ＝  ０  （xが０の場合）                          式（Ｔ−１）
invTable[x]  ＝  Floor((2^ ShiftA2 / x/2) / x)  （xが０以外の場合）    式（Ｔ−２）
なお、Floor(x)は、小数点以下を切り捨てにより整数化する関数である。式（Ｔ−１）の代わりに、以下の式（Ｔ−２´）を用いても良い。すなわち、除数xの1/2倍を加える丸目調整を行わなくても良い。
【０１２４】
invTable[x]  ＝  Floor(M / x)  （xが０以外の場合）    式（Ｔ−２´）
逆数テーブル値invTable[]を用いることにより、a2sでの除算に相当する演算を、a2sの逆数に相当する逆数テーブル値invTable[a2s]との積と、log2(M)に相当する右シフトにより実現することができる。
【０１２５】
パラメータａの値は、第１パラメータa1と第２パラメータa2の比(a1/a2をiShiftだけ左シフトした値に相当する)。
【０１２６】
導出されたパラメータａは、ＬＭパラメータｂ導出部３１０４１７とＬＭ予測フィルタ部３１０４２に出力される。
【０１２７】
ＬＭパラメータｂ導出部３１０４１７は、目的参照画素yの和Yから、入力参照画素xの和Xにパラメータａをかけて固定シフト値iShiftだけ右シフトした値を引いた値を、参照領域の画素数で割ることにより、以下の式によりパラメータｂを導出する。
【０１２８】
b = ( Y - ( ( a * X ) >> iShift ) + ( 1 << ( iCountShift - 1 ) ) ) >> iCountShift          式（Ｂ−２０）
なお、iCountShiftの右シフトは、参照領域の画素数Ｎで割ることに相当する。
【０１２９】
（ＬＭ予測フィルタ部３１０４２）
ＬＭ予測フィルタ部３１０４２は、前述の輝度成分から色差成分を予測するＴＹＰＥ０や、色差成分から別の色差成分を予測するＴＹＰＥ１の予測の場合には、ＬＭパラメータ推定部３１０４１が導出した推定パラメータを用いて、ＬＭ予測前の入力参照画素における予測ブロックに対応する予測画像predSamples []から、ＬＭ予測後の予測画像predSamples´[]を導出する。例えば、パラメータｂを式（Ｂ−２０）により導出する場合には以下の式を用いる。
【０１３０】
predSamples´[x][y] = (a * predSamples[x][y] >> iShift)+b           式（Ｐ−１）
なお、色コンポーネントが４２０の場合など、入力参照画素predSamples[x][y]と予測画像predSamples´[]のサンプリングが異なる場合には、サンプリング後の画素をpredSamples[x][y]として用いても良い。例えば、入力参照画素predSamples[x][y]が輝度成分、予測画像predSamples´[]が色差成分の場合には、入力参照画素predSamples[x][y]として、入力参照画素の予測画像を２：１にサブサンプルしたブロック、もしくは、入力参照画素の予測画像にフィルタ処理をしてから２：１にサブサンプルしたブロックを用いても良い。
【０１３１】
（ＬＭ予測フィルタ部３１０４２´）
ＬＭ予測フィルタ部３１０４２の代わりに、ＬＭ予測フィルタ部３１０４２の別の構成であるＬＭ予測フィルタ部３１０４２´を用いる。ＬＭ予測フィルタ部３１０４２´は、ＬＭ予測前の予測画像predSamples []から、以下の式により、ＬＭ予測後の予測画像predSamples´[]を導出する。
【０１３２】
predSamples´[x][y] = (a * predSamples[x][y] + b >> iShift)           式（Ｐ−１´）
なお、ＬＭ予測フィルタ部３１０４２´のようにパラメータｂを加算してからシフトする場合には、ＬＭパラメータｂ導出部３１０４１７の代わりに、ＬＭパラメータｂ導出部３１０４１７の別の構成ＬＭパラメータｂ導出部３１０４１７´を用いても良い。この場合、目的参照画素yの和Yを固定シフト値iShiftだけ左シフトした値から、目的参照画素xの和Xにパラメータａをかけた値を引いた値を、参照画素の数で割ることにより、以下の式によりパラメータｂを導出しても良い。
【０１３３】
b = ( (Y << iShift) - ( ( a * X ) ) + ( 1 << ( iCountShift - 1 ) ) ) >> iCountShift          式（Ｂ−２０´）
（ＬＭ予測フィルタ部３１０４２Ｂ）
ＬＭ予測フィルタ部３１０４２の別の構成として、残差を利用した以下のＬＭ予測フィルタ部３１０４２Ｂを用いても良い。
【０１３４】
ＬＭ予測フィルタ部３１０４２Ｂは、ＬＭパラメータ推定部３１０４１が導出した推定パラメータを用いて、プレーナ予測、ＤＣ予測、方向予測などで導出した目的参照信号に対応する対象色成分の予測画像predSamples[x][y]と、入力参照信号に対応する参照色成
分（例えば第１色差成分のCb成分）の残差画像resSamples[x][y]から、目的参照信号に対応するＬＭ予測後の（例えば第２色差成分のCr成分の）予測画像predSamples´[]を導出する。
【０１３５】
predSamples´[x][y] = predSamples[x][y] + (a * resSamples[x][y] >> iShift)          式（Ｐ−２）
上記のように、ＬＭ予測フィルタ部３１０４２Ｂの動作は、線形予測を利用して参照色コンポーネントでの残差画像から、対象色コンポーネントでの残差画像を予測し、予測した残差画像を用いて、対象色コンポーネントでの予測画像を更新する。
【０１３６】
残差画像を利用した構成は、特に色差成分から別の色差成分を予測するＴＹＰＥ１において有効であるが、輝度成分から色差成分を予測するＴＹＰＥ０で利用することもできる。
【０１３７】
（ＬＭパラメータ推定部３１０４１Ａ）
以下、本実施形態に係るＬＭパラメータ推定部３１０４１一つの構成としてＬＭパラメータ推定部３１０４１Ａを説明する。ＬＭパラメータ推定部３１０４１Ａ（ＬＭ正則化コスト加算部３１０４１８Ａ）は、第１パラメータa1が所定の閾値THNより小さい場合、もしくは、第１パラメータa1が所定の閾値THN以下の場合に、第１パラメータa1に０以上の第１コストa1costNを減算することを特徴とする。
【０１３８】
なお、０以上の正則化コストの減算は、０以下の正則化コストの加算と等価であるから、正則化コストの減算の構成は、全て正則化コストの加算で置き換えることができる（以下、同様）。
【０１３９】
図１９は、ＬＭパラメータａ導出部３１０４１６が備えるＬＭ正則化コスト加算部３１０４１８（ＬＭ正則化コスト加算部３１０４１８Ａ、ＬＭ正則化コスト加算部３１０４１８Ｂ、ＬＭ正則化コスト加算部３１０４１８Ｃ）の構成を示すブロック図である。ＬＭ正則化コスト加算部３１０４１８は、ＬＭ正則化コスト導出部３１０４１８０、ＬＭ正則化コスト閾値判定部３１０４１８３、ＬＭ第１パラメータ正則化コスト加算部３１０４１８１から構成される。
【０１４０】
ＬＭ正則化コスト導出部３１０４１８０は、正則化コストacostを導出する。導出した正則化コストacostを、第１パラメータコストa1costNとして用いる。ＬＭ正則化コスト導出部３１０４１８０の詳細については後述する。
【０１４１】
本実施形態に係るＬＭパラメータ推定部３１０４１Ａ（ＬＭ正則化コスト加算部３１０４１８Ａ）は、第１パラメータa1の大きさと閾値との判定を行い、判定に基づいて、第１パラメータa1に正則化コストを加算する。より具体的には、ＬＭ正則化コスト閾値判定部３１０４１８３は、第１パラメータa1と閾値THNを比較し、第１パラメータa1が閾値THN未満（もしくは第１パラメータa1が閾値THN以下）の場合に、ＬＭ第１パラメータ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costNを減算する。逆に、第１パラメータa1が閾値THN以上の場合には、第１パラメータ正則化コストa1costNを減算しない。
【０１４２】
a1´ = a1 - a1costN
ＬＭ正則化コスト加算部３１０４１８Ａの動作を、図２１のＳ１１０１、Ｓ１１０２を用いてフローチャートを用いて、再度、説明する。
【０１４３】
Ｓ１１０１  ＬＭ正則化コスト閾値判定部３１０４１８３は、第１パラメータa1と閾値THNを比較する。第１パラメータa1が閾値THN未満（もしくは第１パラメータa1が閾値THN以下）の場合に、Ｓ１１０２に遷移する。
【０１４４】
Ｓ１１０２  ＬＭパラメータａ導出部３１０４１６は、ＬＭ第１パラメータ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costNを減算する。減算後のパラメータa1をa1´とも呼ぶ。
【０１４５】
ＬＭ正則化コスト加算部３１０４１８Ｂの動作はここで終了する。
【０１４６】
なお、ＬＭ正則化コスト加算部３１０４１８Ｂは、第２パラメータa2に正則化コストa2costNを加算しないが、後述のＬＭ正則化コスト加算部３１０４１８Ｂ２のように、Ｓ１２０２に続いてＳ１２０３を実行し、第２パラメータa2に第２パラメータ正則化コストa2costNを加算しても良い。
【０１４７】
Ｓ１１０３  ＬＭ第２パラメータ正則化コスト加算部３１０４１８２は、第２パラメータa2に第２パラメータ正則化コストa2costNを加算する。加算後の第２パラメータa2をa2´とも呼ぶ。
【０１４８】
なお、正則化コストの加算は、上記のように、ＬＭ第１パラメータ導出部３１０４１４および、ＬＭ第２パラメータ導出部３１０４１５で、第１パラメータa1、第２パラメータa2を導出した後に別の装置であるＬＭ正則化コスト加算部３１０４１８で行うのではなく、ＬＭ第１パラメータ導出部３１０４１４および、ＬＭ第２パラメータ導出部３１０４１５で第１パラメータa1、第２パラメータa2を導出する際に、同時に正則化コストを加算しても良い。この構成のＬＭパラメータ推定部３１０４１は、ＬＭ第１パラメータ導出部３１０４１４および、ＬＭ第２パラメータ導出部３１０４１５において、前述の式（Ｂ−７）、式（Ｂ−７´）、式（Ｂ−７´´）の代わりに、各々、以下の式（ＢＣ−７）、式（ＢＣ−７´）、式（ＢＣ−７´´）の式を用いる。
【０１４９】
a1 = (XY<<iCountShift) - X*Y - a1costN;                 式（ＢＣ−７）
a1 = XY - ((aveY * aveX)<<iCountShift) - a1costN;         式（ＢＣ−７´）
a1 = XY - ((aveY * aveX)<<iCountShift) + (aveX * remY) + (aveY * remX) - a1costN;         式（ＢＣ−７´´）
（具体例）
以下、正則化コストa1costNに具体的な値を用いた場合の例を説明する。ここでは、a1costNを第２パラメータa2に基づいて導出する（ＬＭ正則化コスト導出部３１０４１８０Ｄ、式（Ｃ−１Ｄ）参照）。例えば、a1costNがa2のk倍をLで割った以下の値である場合、
a1costN = k * a2 / L
正則化コストを加算後（減算後）の第１パラメータa1´、第２パラメータa2´を用いると、a1/a2に相当するパラメータaは以下の式で導出される。
【０１５０】
a1´/a2´ = (a1 + a1costN) / a2 = a1/a2 - k / L。
【０１５１】
ここでk=1、L=32の場合には、k / L = 0.0315であり、正則化コストの効果は、a1が閾値THN未満である場合に、パラメータa（a1/a2）を小さな値のオフセット0.0315だけ減算することと等価である。
【０１５２】
なお、ＬＭ正則化コスト閾値判定部３１０４１８３Ａは、閾値THNとして、式（Ｄ−１）のように0もしくは、式（Ｄ−１´）のようにa2の倍数に設定すると良いことが、発明者らの実験により明らかになっている。
【０１５３】
THN = 0                    式（Ｄ−１）
THN = -(a2 >> tshift2)        式（Ｄ−１´）
（説明）
発明者の知見では、最小２乗法により導出されるパラメータaは、雑音などにより、予測画像の生成に最適な値よりも、絶対値が小さくなることが多い（最適な値よりも０に近い側になる）ことが分かった。正則化コストによって、パラメータaを調整する場合には、以下の２通りがある。パラメータaを正の方向に近づけて絶対値を大きく場合には、０以上の正則化コストをパラメータa1に加算する必要があるのに対して、パラメータaを負の方向に近づけて絶対値を大きくするには、０以上の正則化コストをパラメータa1から減算する必要がある。このように、パラメータaの符号により、正則化コストを加算するか減算するかを変更する必要がある。そのためには、正則化コストで第１パラメータa1の大きさを調整する前にパラメータaの符号を判定する必要がある。この判定のために、本実施形態に係るＬＭパラメータ推定部３１０４１（変形例を含む）は、第１パラメータa1の符号、もしくは、第１パラメータa1と所定の閾値との比較結果を利用する。
【０１５４】
パラメータaは、a1/a2の整数倍(1<<iShift倍)に相当する値として導出される。第２パラメータa2は常に０以上であるから、パラメータaの符号（＝a1/a2の符号）は、第１パラメータa1の符号と等しい。従って、パラメータaの符号及び大きさは、第１パラメータa1と所定の閾値との比較で導出できる。
【０１５５】
ＬＭパラメータ推定部３１０４１は、例えばa1がTHN未満（ここではTHN=0）である場合に、a1から正則化コストa1costNだけ減算することにより、パラメータaの値を、最小２乗法で導出される値よりも０から遠い方向（ここでは負の方向）に近づけることが可能であることが理解できる。上記、構成により、ＬＭパラメータ推定部３１０４１は、雑音に強い線形予測パラメータを導出することができ、予測画像を生成する際の予測精度を向上させる効果を奏する。
【０１５６】
なお、a1costN、a1costP、a2costN、a2costPは、各々、パラメータaを概略負(Negative)の方向に近づけるa1への正則化コスト、パラメータaを概略負の方向に近づけるa2への正則化コスト、パラメータaを概略正(Positive)の方向に近づけるa1への正則化コスト、パラメータaを概略正の方向に近づけるa2への正則化コストを意味する。これらを区別しない場合には単に、acost、a1cost、a2cost、acostN、acostPなどを用いる。
【０１５７】
なお、０から遠い方向に近づけるという意味では、所定の閾値としてTHN=0が好適である。すなわち、ＬＭ正則化コスト閾値判定部３１０４１８３は、第１パラメータa1と閾値THNを比較し、第１パラメータa1が０以上の閾値THN未満（もしくは第１パラメータa1が閾値THN以下）の場合に、第１パラメータＬＭ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costNを加算する構成が好適である。
【０１５８】
しかしながら、画像によってはTHN=0に近い別の値も好適である。例えば、THN = -a2>>tshift(ここではtshiftは2〜8の定数)など。例えばtshift = 4の場合には、a1が、-a2/16未満、つまりa1/a2が-1/16未満の場合に、０から遠い方向（ここでは負の方向）に近づけることが可能である。パラメータaの値が０に近い場合に無理して拘束項をきかせないことが可能である。また、逆に、THN = a2>>tshift(ここではtshiftは2〜8の定数)とすることも可能である。例えばtshift = 4の場合には、a1が、a2/16未満、つまりa1/a2が1/16未満の場合に、０から遠い方向（ここでは負の方向）に近づけることが可能である。ここでは、パラメータaが０に非常に近い正の場合にも正則化コストによってパラメータaの値を負の方向に近づけることができる。
【０１５９】
なお、本実施形態では、ＬＭパラメータ推定部３１０４１は、パラメータaが負となる場合において、正則化コストを用いてパラメータaの値をさらに負の方向に近づけることにより、パラメータaを０から通い方向に近づけたが、パラメータaが正となる場合を考慮して、正則化コストを用いてパラメータaの値を正の方向に近づけることにより、パラメータaを０から遠い方向に近づける場合も可能である。この例を、ＬＭパラメータ推定部３１０４１Ｂとして説明する。さらに、パラメータaが負と正の場合の両方を考慮した構成をＬＭパラメータ推定部３１０４１Ｃとして説明する。
【０１６０】
（ＬＭパラメータ推定部３１０４１Ｂ）
以下、ＬＭパラメータ推定部３１０４１一つの構成としてＬＭパラメータ推定部３１０４１Ｂを説明する。ＬＭパラメータ推定部３１０４１Ｂ（ＬＭ正則化コスト加算部３１０４１８Ｂ）は、第１パラメータa1が所定の閾値THPより大きい場合、もしくは、第１パラメータa1が所定の閾値THP以上の場合に、第１パラメータa1に０以上の第１パラメータ正則化コストa1costPを加算することを特徴とする。
【０１６１】
ＬＭ正則化コスト加算部３１０４１８Ｂの動作を、図２２のフローチャートのＳ１２０１、Ｓ１２０２を用いて、再度、説明する。
【０１６２】
Ｓ１２０１  ＬＭ正則化コスト閾値判定部３１０４１８３は、第１パラメータa1と閾値THNを比較する。第１パラメータa1が閾値THPより大きい（もしくは第１パラメータa1が閾値THP以上）の場合に、Ｓ１２０２に遷移する。
【０１６３】
Ｓ１２０２  第１パラメータＬＭ正則化コスト加算部３１０４１８１は、第１パラメータＬＭ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costPを加算する。加算後のパラメータa1をa1´とも呼ぶ。
【０１６４】
ＬＭ正則化コスト加算部３１０４１８Ｂの動作はここで終了する。
【０１６５】
なお、ＬＭ正則化コスト加算部３１０４１８Ｂは、第２パラメータa2に正則化コストa2costPを加算しないが、後述のＬＭ正則化コスト加算部３１０４１８Ｂ２のように、Ｓ１２０２に続いてＳ１２０３を実行し、第２パラメータa2に第２パラメータ正則化コストa2costPを加算しても良い。
【０１６６】
Ｓ１２０３  ＬＭ第２パラメータ正則化コスト加算部３１０４１８２は、第２パラメータa2に第２パラメータ正則化コストa1costPを加算する。加算後のパラメータa2をa2´とも呼ぶ。
【０１６７】
ＬＭ正則化コスト導出部３１０４１８０は、正則化コストacostを導出する。上述の式（Ｃ−１Ａ）から式（Ｃ−１Ｄ２）の何れか一項、もしくは、その組み合わせにより導出したacostを、第１パラメータコストの第１コストa1costPとして導出する。
【０１６８】
本実施形態に係るＬＭ正則化コスト加算部３１０４１８Ｂは、第１パラメータa1の大きさと閾値との判定を行い、判定に基づいて、第１パラメータa1に正則化コストを加算する。より具体的には、ＬＭ正則化コスト閾値判定部３１０４１８３Ｂは、第１パラメータa1と０以上の閾値THPを比較し、第１パラメータa1が閾値THPより大きい（もしくは第１パラメータa1が閾値THP以上の）場合に、ＬＭ第１パラメータ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costPを加算する。逆に、第１パラメータa1が閾値THP以下（もしくは閾値THP未満）の場合には、第１パラメータ正則化コストa1costPを加算しない。
【０１６９】
a1 = a1 + a1costP
なお、正則化コストの加算は、ＬＭ第１パラメータ導出部３１０４１４および、ＬＭ第２パラメータ導出部３１０４１５で行っても良い。導出式は、既に説明した式（ＢＣ−７）、式（ＢＣ−７´）、式（ＢＣ−７´´）の式において、a1costNをa1costPに置き替えた式となる。
【０１７０】
（具体例）
以下、正則化コストa1costNに具体的な値を用いた場合の例を説明する。ここでは、a1costNを第２パラメータa2に基づいて導出する（ＬＭ正則化コスト導出部３１０４１８０Ｄ、式（Ｃ−１Ｄ）参照）。
【０１７１】
a1costP = k * a2 / L
(a1 + a1costP) / a2 = a1/a2 - k / L。
【０１７２】
k=1、L=32の場合には、k / L = 0.0315であるから、正則化コストは、a1が閾値THPより大きい場合に、パラメータa（a1/a2）を小さな値のオフセット0.0315だけ加算することと等価である。従って、例えばa1がTHPより大きい（ここではTHP=0）である場合に、正則化コストa1costPをa1に加算することにより、パラメータaの値を、最小２乗法で導出される値よりも０から遠い方向（ここでは正の方向）に近づけることが可能であり、予測画像を生成する際の予測精度を向上させる効果を奏する。
【０１７３】
なお、ＬＭ正則化コスト閾値判定部３１０４１８３Ｂは、閾値THPとして、式（Ｄ−２）のように0もしくは、式（Ｄ−２´）のようにa2の倍数に設定すると良いことが、発明者らの実験により明らかになっている。
【０１７４】
THP = 0                    式（Ｄ−２）
THP = a2 >> tshift2        式（Ｄ−２´）
（ＬＭパラメータ推定部３１０４１Ｃ）
以下、ＬＭパラメータ推定部３１０４１の一つの構成としてＬＭパラメータ推定部３１０４１Ｃを説明する。ＬＭパラメータ推定部３１０４１Ｃは、正則化コストの加算（ＬＭ正則化コスト加算部３１０４１８）において、ＬＭパラメータ推定部３１０４１ＡとＬＭパラメータ推定部３１０４１Ｂを組み合わせた構成である。ＬＭパラメータ推定部３１０４１Ｃは、第１パラメータa1が所定の閾値THNより小さい場合、もしくは、第１パラメータa1が所定の閾値THN以下の場合に、第１パラメータa1に０以上の第１コストa1costNを減算し、さらに、第１パラメータa1が所定の閾値THPより大きい場合、もしくは、第１パラメータa1が所定の閾値THP以上の場合に、第１パラメータa1に０以上の第２コストa1costPを加算することを特徴とする。
【０１７５】
ＬＭパラメータ推定部３１０４１Ｃ（ＬＭ正則化コスト加算部３１０４１８Ｃ）の動作を、図２３のフローチャートを用いて、再度、説明する。
【０１７６】
Ｓ１３０１  ＬＭ正則化コスト閾値判定部３１０４１８３は、第１パラメータa1と閾値THNを比較する。第１パラメータa1が閾値THN未満（もしくは第１パラメータa1が閾値THN以下）の場合に、Ｓ１３０２に遷移する。
【０１７７】
Ｓ１３０２  第１パラメータＬＭ正則化コスト加算部３１０４１８１は、第１パラメータＬＭ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costNを加算する。加算後のパラメータa1をa1´とも呼ぶ。
【０１７８】
ＬＭ正則化コスト加算部３１０４１８Ｃは、Ｓ１３０３をスキップしてＳ１３０４に遷移する。なお、ＬＭ正則化コスト加算部３１０４１８Ｃは、Ｓ１３０３の動作を行わないが、後述のＬＭ正則化コスト加算部３１０４１８Ｃ２のように、Ｓ１３０２に続いてＳ１３０３を実行し、第２パラメータa2にa２パラメータ正則化コストa2costNを加算する構成でも良い。
【０１７９】
Ｓ１３０３  ＬＭ第２パラメータ正則化コスト加算部３１０４１８２は、ＬＭ第２パラメータ正則化コスト加算部３１０４１８２を用いて、第２パラメータa2に第２パラメータ正則化コストa2costNを加算する。加算後のパラメータa2をa2´とも呼ぶ。
【０１８０】
Ｓ１３０４  ＬＭ正則化コスト閾値判定部３１０４１８３は、第１パラメータa1と閾値THPを比較する。第１パラメータa1が閾値THPより大きい（もしくは第１パラメータa1が閾値THP以上）の場合に、Ｓ１３０５に遷移する。
【０１８１】
Ｓ１３０５  第１パラメータＬＭ正則化コスト加算部３１０４１８１は、ＬＭ第１パラメータ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costPを加算する。加算後のパラメータa1をa1´とも呼ぶ。
【０１８２】
ＬＭ正則化コスト加算部３１０４１８Ｃの動作はここで終了する。なお、ＬＭ正則化コスト加算部３１０４１８Ｃは、Ｓ１３０６の動作を行わないが、後述のＬＭ正則化コスト加算部３１０４１８Ｃ２のように、Ｓ１３０５に続いてＳ１３０６を実行し、第２パラメータa2に第２パラメータ正則化コストa2costPを加算する構成でも良い。
【０１８３】
Ｓ１３０６  ＬＭ第２パラメータ正則化コスト加算部３１０４１８２は、第２パラメータＬＭ正則化コスト加算部３１０４１８２を用いて、第２パラメータa2に第２パラメータ正則化コストa2costPを加算する。加算後のパラメータa2をa2´とも呼ぶ。
【０１８４】
ＬＭ正則化コスト導出部３１０４１８０は、正則化コストacostを導出する。例えば、既に説明した式（Ｃ−１Ａ）から式（Ｃ−１Ｄ２）の何れか一項、もしくは、その組み合わせにより導出する。
【０１８５】
本実施形態に係るＬＭ正則化コスト加算部３１０４１８Ｃは、第１パラメータa1の大きさと閾値との判定を行い、判定に基づいて、第１パラメータa1及び第２パラメータa2に正則化コストを加算する。より具体的には、ＬＭ正則化コスト閾値判定部３１０４１８３Ｃは、第１パラメータa1と０以下の閾値THNを比較し、第１パラメータa1が閾値THNより小さい（もしくは第１パラメータa1が閾値THN以下）の場合に、第１パラメータＬＭ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1から第１パラメータ正則化コストa1costNを減算する。それ以外の場合、ＬＭ正則化コスト閾値判定部３１０４１８３Ｃは、第１パラメータa1と０以上の閾値THPを比較し、第１パラメータa1が閾値THPより大きい（もしくは第１パラメータa1が閾値THP以上）の場合に、第１パラメータＬＭ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costPを加算する。逆に、第１パラメータa1が閾値THN以上の場合には、第１パラメータ正則化コストa1costPを加算しない。
【０１８６】
ここで、減算する場合に用いる正則化コストa1costNと加算する場合の正則化コストa1costPとして、全て同じ値を用いても良い（a1costN = a1costP = acost）し、異なる値を用いても良い。例えば、加算する場合の正則化コストa1costPを減算する場合の第１パラメータ正則化コストa1costN以下として以下のように導出しても良い。
【０１８７】
a1costN = (kM*a2)>>kshiftM
a1costP = (kP*a2)>>kshiftP
ここでkM>kPまたはkshiftP > kshiftM。例えばkM=kP=1、kshiftM=5、kshiftP=6とすれば良い。
【０１８８】
以上までに説明したＬＭパラメータ推定部３１０４１Ａ、ＬＭパラメータ推定部３１０４１Ｂ、ＬＭパラメータ推定部３１０４１Ｃは、パラメータaを０から離れる側に近づけると良い、という発明者の知見に基づいていたが、さらに詳細に検討すると、パラメータaの絶対値が十分に大きい場合には、パラメータaをさらに０から遠い側に近づけると、パラメータaの絶対値が本来の最適な値よりも大きくなりすぎる場合がある、ことが分かった。これは、雑音の影響として、基本的にはパラメータaが小さくなる可能性が比較的大きいが、場合によって、雑音の影響でパラメータaが大きくなっていることがあるからであり、特に最小２乗法で求められる正則化コストなしで導出されるパラメータaが十分に大きい場合には、その確率が高くなる。以下に説明するＬＭパラメータ推定部３１０４１Ａ２、ＬＭパラメータ推定部３１０４１Ｂ２、ＬＭパラメータ推定部３１０４１Ｃ２では、構成が若干複雑になる代わりに、これらの問題を解決する。
【０１８９】
ＬＭパラメータ推定部３１０４１Ａ、ＬＭパラメータ推定部３１０４１Ｂ、ＬＭパラメータ推定部３１０４１Ｃは、ＬＭ正則化コストとして第１パラメータa1に正則化超す音を加算もしくは減算したが、ＬＭパラメータ推定部３１０４１Ａ２、ＬＭパラメータ推定部３１０４１Ｂ２、ＬＭパラメータ推定部３１０４１Ｃ２は、第１パラメータa1に正則化コストを加算もしくは減算しながら同時に第２パラメータa2に正則化コストを加算する構成になる。
【０１９０】
図２０は、ＬＭパラメータ推定部３１０４１が備えるＬＭ正則化コスト加算部３１０４１８（ＬＭ正則化コスト加算部３１０４１８Ａ２、ＬＭ正則化コスト加算部３１０４１８Ｂ２、ＬＭ正則化コスト加算部３１０４１８Ｃ２）の構成を示すブロック図である。ＬＭ正則化コスト加算部３１０４１８は、ＬＭ正則化コスト導出部３１０４１８０、ＬＭ正則化コスト閾値判定部３１０４１８３、ＬＭ第１パラメータ正則化コスト加算部３１０４１８１、ＬＭ第２パラメータ正則化コスト加算部３１０４１８２から構成される。
【０１９１】
図２０のＬＭパラメータ推定部３１０４１は、上記第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN未満もしくは上記閾値THN以下である場合には第１コストa1costNを上記第１パラメータa1から減算し、上記第２パラメータa2にコストa2costNを加算する手段と、上記第１パラメータa1と所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THP未満もしくは上記閾値THP以下の場合には第２コストa1costPを上記第１パラメータa1に加算し上記第２パラメータa2にコストa2costPを加算する手段の少なくとも１つの手段を備えることを特徴とする。
【０１９２】
（ＬＭパラメータ推定部３１０４１Ａ２）
以下、ＬＭパラメータ推定部３１０４１の一つの構成としてＬＭパラメータ推定部３１０４１Ａ２を説明する。既に説明したＬＭパラメータ推定部３１０４１Ａでは、第１パラメータa1に第１パラメータ正則化コストa1costNを減算していたが、ＬＭパラメータ推定部３１０４１Ａ２では、第１パラメータa1に正則化コストa1costNを減算するのに加えて、第２パラメータa2に第２パラメータ正則化コストa2costNを加算する。
【０１９３】
ＬＭパラメータ推定部３１０４１Ａ２（ＬＭ正則化コスト加算部３１０４１８Ａ２）の動作は、図２１のフローチャートの通りである。図２１の各ステップの詳細は、既にＬＭ正則化コスト加算部３１０４１８Ａで説明したので省略する。
【０１９４】
ＬＭ正則化コスト加算部３１０４１８Ａ２は、第１パラメータa1が所定の閾値THNより小さい場合、もしくは、第１パラメータa1が所定の閾値THN以下の場合（Ｓ１１０１）に、第１パラメータa1に０以上の第１コストa1costNを減算（Ｓ１１０２）し、第２パラメータa2に０以上の第２コストa2costNを加算する（Ｓ１１０３）。
【０１９５】
第１パラメータa1にa1costNを減算し、第２パラメータa2にa2costNだけ加算する場合、パラメータａは、-a1costN/a2costNに近づくように補正がかかる。
【０１９６】
ＬＭ正則化コスト導出部３１０４１８０は、正則化コストacostを導出する。例えば、既に説明した式（Ｃ−１Ａ）から式（Ｃ−１Ｄ２）の何れか一項、もしくは、その組み合わせにより導出する。さらに正則化コストacostから、a1costN及びa2costNを以下の式により導出する。
【０１９７】
a1costN = (acost>>kshift)      
a2costN = acost                
上記の場合、パラメータａは、-a1costN/a2costN=-(1>>kshift)に近づくように補正がかかる。具体的には、1、0.5、0.25、0.125に近づくようなkshift=0、1、2、3が適当である。
【０１９８】
なお、シフトにより計算可能な上記正則化項は計算が容易であるが、正則化コストは上記に限定されない。例えば、以下のようにすれば、k>>kshiftに近づくように補正がかかる。
【０１９９】
a1costN = (k*acost)>>kshift      
a2costN = acost                  
例えば、k=17、kshift=5とすれば、-17/32=-0.53125に近づけるような正則化項が導入できる。
【０２００】
上記のように、a2costNは、上記コストa1costN以上であることが好適である、すなわち、a1costN / a2costNが１以下となることが好適である。
【０２０１】
本実施形態に係るＬＭパラメータ推定部３１０４１Ａ２（ＬＭ正則化コスト加算部３１０４１８Ａ２）は、第１パラメータa1の大きさと閾値との判定を行い、判定に基づいて、第１パラメータa1及び第２パラメータa2に正則化コストを加算する。より具体的には、ＬＭ正則化コスト閾値判定部３１０４１８３Ａ２は、第１パラメータa1と０以上の閾値THNを比較し、第１パラメータa1が閾値THN未満（もしくは第１パラメータa1が閾値THN以下）の場合に、第１パラメータＬＭ正則化コスト加算部３１０４１８１を用いて、第１パラメータa1に第１パラメータ正則化コストa1costNを加算し、第２パラメータa2に第２パラメータ正則化コストa2costを加算する。逆に、第１パラメータa1が閾値THN以上の場合には、第１パラメータa1、第２パラメータa2に正則化コストを加算しない。
【０２０２】
上記構成のＬＭパラメータ推定部３１０４１Ａ２（ＬＭ正則化コスト加算部３１０４１８Ａ２）によれば、第１パラメータa1に第１パラメータ正則化コストa1costNを減算し、第２パラメータa2に第２パラメータ正則化コストa2costNを加算することにより、パラメータaの値を、最小２乗法で導出される値よりも-a1costN / a2costNに近い方向に近づけることが可能であり、予測画像を生成する際の予測精度を向上させる効果を奏する。ＬＭ正則化コスト加算部３１０４１８Ａは単純にパラメータaの値を、最小２乗法で導出される値よりも０から遠い方向に近づけるだけであったのに対して、ＬＭパラメータ推定部３１０４１Ａ２（ＬＭ正則化コスト加算部３１０４１８Ａ２）は、第２パラメータa2にも第２パラメータ正則化コストa2costNを加算することにより、-a1costN / a2costNに近づける。輝度色差予測に適した傾きaは、aが負の場合には、０と−１の間であるので、ＬＭ正則化コスト加算部３１０４１８Ｂのように、最小２乗法で導出される値が負の値として十分に小さい（-a1costN / a2costNよりも小さい）場合に、さらに小さくすることによる弊害をなくすことができ、ＬＭパラメータ推定部３１０４１Ａ（ＬＭ正則化コスト加算部３１０４１８Ａ）よりもさらに精度の良い線形予測パラメータを導出することができる効果を奏する。
【０２０３】
（具体例）
以下、正則化コストa1costNに具体的な値を用いた場合の例を説明する。ここでは、a1costNを第２パラメータa2に基づいて導出する（ＬＭ正則化コスト導出部３１０４１８０Ｄ、式（Ｃ−１Ｄ）参照）。
【０２０４】
発明者の実験では、正則化コストとして、例えば、例１から例４の以下の組み合わせが好ましいことが実験で得られている。
【０２０５】
例１：acost＝(a2>>4)、kshift = 1
a1costN = (k*acost)>>kshift= a2>>5
a2costN = acost = a2>>4
例２：acost＝(a2>>3)、kshift = 1
a1costN = (k*acost)>>kshift= a2>>4
a2costN = acost = a2>>3
例３：acost＝(a2>>3)、kshift = 2
a1costN = (k*acost)>>kshift= a2>>5
a2costN = acost = a2>>3
例４：acost＝(a2>>2)、kshift = 2
a1costN = (k*acost)>>kshift= a2>>4
a2costN = acost= a2>>2
なお、上記、正則化コストの組み合わせは、好適な効果を上げることが確認されているが、本発明の形態は上記の値に限定されない。
【０２０６】
（ＬＭパラメータ推定部３１０４１Ｂ２）
以下、ＬＭパラメータ推定部３１０４１の一つの構成としてＬＭパラメータ推定部３１０４１Ｂ２を説明する。既に説明したＬＭパラメータ推定部３１０４１Ｂでは、第１パラメータa1に第１パラメータ正則化コストa1costNを加算していたが、ＬＭパラメータ推定部３１０４１Ｂ２では、第１パラメータa1に第１パラメータ正則化コストa1costNを加算するのに加えて、第２パラメータa2に第２パラメータ正則化コストa2costNを加算する。
【０２０７】
ＬＭ正則化コスト加算部３１０４１８Ｂ２の動作は、図２１のフローチャートの通りである。図２１の各ステップの詳細は、既にＬＭ正則化コスト加算部３１０４１８Ｂで説明したので省略する。
【０２０８】
本実施形態に係るＬＭ正則化コスト加算部３１０４１８Ｂ２は、第１パラメータa1が所定の閾値THPより大きい場合、もしくは、第１パラメータa1が所定の閾値THP以上の場合に（Ｓ１２０１）、第１パラメータa1に０以上の第１コストa1costPを減算し（Ｓ１２０２）、第２パラメータa2に０以上の第２コストa2costPを加算する（Ｓ１２０３）。
【０２０９】
上記構成のＬＭパラメータ推定部３１０４１Ｂ２（ＬＭ正則化コスト加算部３１０４１８Ｂ２）によれば、第１パラメータa1に第１パラメータ正則化コストa1costPを加算し、第２パラメータa2に第２パラメータ正則化コストa2costPを加算することにより、パラメータaの値を、最小２乗法で導出される値よりもa1costP / a2costPに近い方向に近づけることが可能であり、予測画像を生成する際の予測精度を向上させる効果を奏する。ＬＭパラメータ推定部３１０４１Ｂ（ＬＭ正則化コスト加算部３１０４１８Ｂ）は単純にパラメータaの値を、最小２乗法で導出される値よりも０から遠い方向に近づけるだけであったのに対して、ＬＭパラメータ推定部３１０４１Ｂ２（ＬＭ正則化コスト加算部３１０４１８Ｂ２）は、第２パラメータa2にも第２パラメータ正則化コストa2costPを加算することにより、a1costP / a2costPに近づける。色コンポーネント間予測において、例えば輝度成分から色差成分を予測する場合に適した傾きaは、aが正の場合には、多くの場合０以上１未満の値であるので、ＬＭ正則化コスト加算部３１０４１８Ｂのように、最小２乗法で導出される値が十分に大きい（a1costP / a2costPよりも大きい）場合に、さらに大きくすることによる弊害をなくすことができ、ＬＭパラメータ推定部３１０４１Ｂよりもさらに精度の良い線形予測パラメータを導出することができる効果を奏する。
【０２１０】
（具体例）
以下、正則化コストa1costＰに具体的な値を用いた場合の例を説明する。ここでは、a1costNを第２パラメータa2に基づいて導出する（ＬＭ正則化コスト導出部３１０４１８０Ｄ、式（Ｃ−１Ｄ）参照）。
【０２１１】
a1costN = a2>>5
a2costN = a2>>4
上記のように、a2costPは、上記コストa1costP以上であることが好適である、すなわち、a1costP / a2costPが１以下となることが好適である。
【０２１２】
（ＬＭ正則化コスト加算部３１０４１８Ｃ２）
以下、ＬＭ正則化コスト加算部３１０４１８の一つの構成としてＬＭ正則化コスト加算部３１０４１８Ｃ２を説明する。既に説明したＬＭ正則化コスト加算部３１０４１８Ｃでは、第１パラメータa1に第１パラメータ正則化コスト（a1costN及びa1costP）を減算もしくは加算していたが、ＬＭ正則化コスト加算部３１０４１８Ｃ２では、第１パラメータa1に第１パラメータ正則化コスト（a1costN及びa1costP）を減算及び加算するのに加えて、第２パラメータa2に第２パラメータ正則化コスト（a2costN及びa2costP）を加算する。
【０２１３】
ＬＭ正則化コスト加算部３１０４１８Ｃ２の動作は、図２１のフローチャートの通りである。図２１の各ステップの詳細は、既にＬＭ正則化コスト加算部３１０４１８Ｃで説明したので省略する。
【０２１４】
本実施形態に係るＬＭ正則化コスト加算部３１０４１８Ｃ２は、第１パラメータa1が所定の閾値THNより小さい場合、もしくは、第１パラメータa1が所定の閾値THN以下の場合に（Ｓ１３０１）、第１パラメータa1に０以上の第１コストa1costNを減算し（Ｓ１３０２）第２パラメータa2に０以上の第２コストa2costNを加算する（Ｓ１３０３）、さらに、第１パラメータa1が所定の閾値THPより大きい場合、もしくは、第１パラメータa1が所定の閾値THP以上の場合に（Ｓ１３０４）、第１パラメータa1に０以上の第３コストa1costPを加算し（Ｓ１３０５）第２パラメータa2に０以上の第４コストa2costPを加算する（Ｓ１３０６）ことを特徴とする。
【０２１５】
上記構成のＬＭパラメータ推定部３１０４１Ｃ２（ＬＭ正則化コスト加算部３１０４１８Ｃ２）によれば、第１パラメータa1に第１パラメータ正則化コストa1costNを減算もしくはa1costPを加算するだけでなく、第２パラメータa2に第２パラメータ正則化コストa2costNもしくはa2costPを加算することにより、パラメータaの値を、最小２乗法で導出される値よりも絶対値としてa1costN / a2costN（もしくはa1costP/a2costP）に近い方向に近づけることが可能であり、予測画像を生成する際の予測精度を向上させる効果を奏する。輝度色差予測に適した傾きaは、aが正の場合には、多くの場合絶対値として０以上１未満の値であるので、ＬＭ正則化コスト加算部３１０４１８Ｃのように、最小２乗法で導出される値が絶対値として十分に大きい（a1costP / a2costPよりも大きい）場合に、さらに大きくすることによる弊害をなくすことができ、ＬＭパラメータ推定部３１０４１Ｃ（ＬＭ正則化コスト加算部３１０４１８Ｃ）よりもさらに精度の良い線形予測パラメータを導出することができる効果を奏する。
【０２１６】
（ＬＭ正則化コスト導出部３１０４１８０）
以下、ＬＭパラメータ推定部３１０４１の備えるＬＭ正則化コスト導出部３１０４１８０の詳細を説明する。ＬＭ正則化コスト導出部３１０４１８０は、例えば、以下の式（Ｃ−１Ａ）から式（Ｃ−１Ｄ２）の何れか一項、もしくは、その組み合わせにより導出しても良い。
【０２１７】
acost  ＝  XX >> ashift                        式（Ｃ−１Ａ）
acost  ＝  XY >> ashift                        式（Ｃ−１Ｂ）
acost  ＝  1<<(2*bitDepth)*N >> ashift          式（Ｃ−１Ｃ）
acost  ＝  1<<(2*bitDepth)<<CountShift >> ashift          式（Ｃ−１Ｃ２）
acost  ＝  k*a2 >> kshift          式（Ｃ−１Ｄ）
acost  ＝  min(k*a2 >> kshift, acost0)          式（Ｃ−１Ｄ２）
なお、式（Ｃ−１Ａ）から式（Ｃ−１Ｄ）のashiftは例えば発明者らの実験では6〜12が好適であることが確認されている。式（Ｃ−１Ｃ）のNは、入力参照信号の数（1<<iCountShift）を用いても良いし、入力参照信号の数が多くの場合8〜16であれば、８から１６の数を用いても良い。式（Ｃ−１Ｄ２）のkとkshiftについては後述する。
【０２１８】
（ＬＭ正則化コスト導出部３１０４１８０Ａ）
以下、ＬＭ正則化コスト導出部３１０４１８０の一つの構成としてＬＭ正則化コスト導出部３１０４１８０Ａを説明する。ＬＭ正則化コスト導出部３１０４１８０Ａは、入力参照信号xと入力参照信号xの積の和XXに基づいて正則化コストを導出する。この値は常に０以上になる。
【０２１９】
acost  ＝  XX >> ashift                        式（Ｃ−１Ａ）
第２パラメータa2は、既に説明した式a2 = XX - ((aveX * aveX)<<iCountShift) のように入力参照信号xと入力参照信号xの積の和XXから所定の値を減算することにより計算される。従って、入力参照信号xと入力参照信号xの積の和XXにより導出される正則化コストは、第２パラメータa2と、同じレベルの値域であることが期待でき、正則化コストとして適当である。本構成のＬＭ正則化コスト導出部３１０４１８０Ａは、入力参照画素xと目的参照画素y[]の値域が比較的等しい場合に適当な正則化コストが導出できるという効果を奏する。
【０２２０】
（ＬＭ正則化コスト導出部３１０４１８０Ｂ）
以下、ＬＭ正則化コスト導出部３１０４１８０の一つの構成としてＬＭ正則化コスト導出部３１０４１８０Ｂを説明する。ＬＭ正則化コスト導出部３１０４１８０Ｂは、入力参照信号xと目的参照信号yの積の和ＸＹに基づいて正則化コストを導出する。この値は常に０以上になる。
【０２２１】
acost  ＝  XY >> ashift                        式（Ｃ−１Ｂ）
第１パラメータa1は、既に説明した式XY - ((aveY * aveX)<<iCountShift)のように入力参照信号xと目的参照信号yの積の和ＸＹから所定の値を減算することにより（計算される。従って、XYにより導出される正則化コストは、第１パラメータa1と、同じレベルの値域であることが期待でき、正則化コストとして適当である。入力参照画素x[]と目的参照画素y[]の積の和から導出されるＸＹは、入力参照画素xの値域や、目的参照画素y[]の値域によらずに適当な正則化コストが導出できるという効果を奏する。
【０２２２】
（ＬＭ正則化コスト導出部３１０４１８０Ｃ）
以下、ＬＭ正則化コスト導出部３１０４１８０の一つの構成としてＬＭ正則化コスト導出部３１０４１８０Ｃを説明する。ＬＭ正則化コスト導出部３１０４１８０Ｃは、bitDepthに基づいて正則化コストを導出する。この値は常に０以上になる。
【０２２３】
acost  ＝  1<<(2*bitDepth)*N >> ashift          式（Ｃ−１Ｃ）
例えば、入力参照画素x[]及び目的参照画素y[]のビットデプス（ビット深度）がbitDepthである場合、x及びyはおおよそ1<<(bitDepTHN)周辺の数である場合が多い。具体的にはbitDepth=8の場合、つまり、輝度及び色差が0から255の値である場合には、輝度、色差の値として128の周辺にあることが比較的多い。従って、実際の入力参照画素x[]や目的参照画素y[]を用いずとも、入力参照画素x及び目的参照画素yとしてビット深度に依存する値（x[i]=y[i]=1<<(bitDepTHN)）を代入して、既に説明した式（Ｃ−１Ａ）、式（Ｃ−１Ｂ）を用いれば、適当な正則化コストが導出されることが期待される。ＬＭ正則化コスト導出部３１０４１８０Ｃは、この考えに基づいて、以下の式で正則化コストを導出する。
【０２２４】
acost  ＝  1<<(2*bitDepth)<<CountShift >> ashift          式（Ｃ−１Ｃ２）
ＬＭ正則化コスト導出部３１０４１８０ＣのbitDepthに基づく方法では、入力参照画素xや目的参照画素y[]を用いないため、比較的小さい演算数で正則化コストが導出できるという効果を奏する。
【０２２５】
（ＬＭ正則化コスト導出部３１０４１８０Ｄ）
以下、ＬＭ正則化コスト導出部３１０４１８０の一つの構成としてＬＭ正則化コスト導出部３１０４１８０Ｄを説明する。ＬＭ正則化コスト導出部３１０４１８０Ｄは、第２パラメータa2に基づいて正則化コストを導出する。例えば以下の式で導出する。この値は常に０以上になる。
【０２２６】
acost  ＝  k*a2 >> kshift          式（Ｃ−１Ｄ）
例えば、k=1、kshift = 2。
【０２２７】
第２パラメータa2は、入力参照画素x[]の分散に相当する。正則化コストは第１パラメータa1及び第２パラメータa2に加算もしくは減算されるが、第２パラメータa2が比較的小さい場合に、第２パラメータa2に比べて大きすぎる正則化コストacostをa2に加算する場合、結果として値が大きく変化しすぎる、課題が生じる。本実施形態では、第２パラメータa2に基づいて正則化コストを導出することにより、正則化コストが大きく効きすぎることを抑制する効果を奏する。
【０２２８】
（ＬＭ正則化コスト導出部３１０４１８０Ｄ２）
さらに、第２パラメータa2に基づいて正則化コストを導出するＬＭ正則化コスト導出部３１０４１８０Ｄは、ＸＸ、ＸＹ、bitDepthなどで導出された正則化コストacostを、a2から導出される値で上限をクリップすることで、正則化コストacostを導出しても良い（ＬＭ正則化コスト導出部３１０４１８０Ｄ２）。この値は常に０以上になる。
【０２２９】
正則化コストは第１パラメータa1及び第２パラメータa2に加算もしくは減算されるが、第２パラメータa2が比較的小さい場合に、第２パラメータa2に比べて大きすぎる正則化コストacostをa2に加算する場合、結果として値が大きく変化しすぎる、という問題が生じる。従って、正則化コストacost自体を第２パラメータa2の大きさに合わせて導出することが適当である。これにより適した正則化コストが大きく効きすぎることを抑制する効果を奏する。
【０２３０】
acost  ＝  min(k*a2 >> kshift, acost0)          式（Ｃ−１Ｄ２）
ここでacost0は、上記までに説明した、ＬＭ正則化コスト導出部３１０４１８０Ａ、ＬＭ正則化コスト導出部３１０４１８０Ｂ、ＬＭ正則化コスト導出部３１０４１８０Ｃの何れかで導出される正則化コストacostである。ＬＭ正則化コスト導出部３１０４１８０Ａは、ベースとなるacost0と、上限を定めるa2から導出される値（ここではk*a2 >> kshift）の最大値を導出することで、導出される正則化コストのついてa2の大きさに応じたクリップを行うことができる。一般に、正則化コストは、雑音などによりa1、a2の信頼度が低い場合に、推定パラメータを安定させるために加算もしくは減算するパラメータであるが、xの分散（ΣXX-ΣXX/N）に相当するa2の値が比較的小さい場合に、第２パラメータa2に比べて大きすぎる正則化コストacostをa2に加算する場合、結果として推定される線形予測パラメータに対して正則化コストが効きすぎるという問題が生じる。従って、正則化コストacostが足される項である第２パラメータa2の大きさに合わせることが適当である。これにより正則化項が大きく効きすぎることを抑制する効果を奏する。
【０２３１】
発明者らの実験では、式（Ｃ−１Ｄ２）が好適な予測精度を有するが、式（Ｃ−１Ｄ）も近い精度を実現することができる。式（Ｃ−１Ｄ）の方が正則化項の導出が容易であることから、式（Ｃ−１Ｄ）を用いることも有効である。
【０２３２】
逆量子化・逆ＤＣＴ部３１１は、エントロピー復号部３０１から入力された量子化係数を逆量子化してＤＣＴ係数を求める。逆量子化・逆ＤＣＴ部３１１は、求めたＤＣＴ係数について逆ＤＣＴ（Inverse Discrete Cosine Transform、逆離散コサイン変換）を行い、復号残差信号を算出する。逆量子化・逆ＤＣＴ部３１１は、算出した復号残差信号を加算部３１２および残差格納部３１３に出力する。
【０２３３】
加算部３１２は、インター予測画像生成部３０９及びイントラ予測画像生成部３１０から入力された予測ピクチャブロックＰと逆量子化・逆ＤＣＴ部３１１から入力された復号残差信号の信号値を画素毎に加算して、参照ピクチャブロックを生成する。加算部３１２は、生成した参照ピクチャブロックを参照ピクチャメモリ３０６に記憶し、生成した参照ピクチャブロックをピクチャ毎に統合した復号レイヤ画像Ｔｄを外部に出力する。
（インター予測パラメータ復号部の構成）
次に、インター予測パラメータ復号部３０３の構成について説明する。
【０２３４】
図６は、本実施形態に係るインター予測パラメータ復号部３０３の構成を示す概略図である。インター予測パラメータ復号部３０３は、インター予測パラメータ復号制御部３０３１、ＡＭＶＰ予測パラメータ導出部３０３２、加算部３０３５及びマージ予測パラメータ導出部３０３６を含んで構成される。
【０２３５】
インター予測パラメータ復号制御部３０３１は、インター予測に関連する符号（シンタックス要素の復号をエントロピー復号部３０１に指示し、符号化データに含まれる符号（シンタックス要素）を例えば、分割モードpart_mode、マージフラグmerge_flag、マージインデックスmerge_idx、インター予測フラグinter_pred_idc、参照ピクチャインデックスrefIdxLX、予測ベクトルインデックスmvp_LX_idx、差分ベクトルmvdLXを抽出する。
【０２３６】
インター予測パラメータ復号制御部３０３１が、あるシンタックス要素を抽出すると表現する場合は、あるシンタックス要素の復号をエントロピー復号部３０１に指示し、該当のシンタックス要素を符号化データから読み出すことを意味する。ここで、マージフラグが示す値が１、すなわち、マージ予測モードを示す場合、インター予測パラメータ復号制御部３０３１は、マージ予測に係る予測パラメータとして、マージインデックスmerge_idxを抽出する。インター予測パラメータ復号制御部３０３１は、抽出したマージインデックスmerge_idxをマージ予測パラメータ導出部３０３６に出力する。
【０２３７】
マージフラグmerge_flagが０、すなわち、ＡＭＶＰ予測モードを示す場合、インター予測パラメータ復号制御部３０３１は、エントロピー復号部３０１を用いて符号化データからＡＭＶＰ予測パラメータを抽出する。ＡＭＶＰ予測パラメータとして、例えば、インター予測フラグinter_pred_idc、参照ピクチャインデックスrefIdxLX、ベクトルインデックスmvp_LX_idx、差分ベクトルmvdLXがある。インター予測パラメータ復号制御部３０３１は、抽出したインター予測フラグinter_pred_idcから導出した予測リスト利用フラグpredFlagLXと、参照ピクチャインデックスrefIdxLXをＡＭＶＰ予測パラメータ導出部３０３２及び予測画像生成部３０８（図５）に出力し、また予測パラメータメモリ３０７（図５）に記憶する。インター予測パラメータ復号制御部３０３１は、抽出したベクトルインデックスmvp_LX_idxをＡＭＶＰ予測パラメータ導出部３０３２に出力する。インター予測パラメータ復号制御部３０３１は、抽出した差分ベクトルmvdLXを加算部３０３５に出力する。
【０２３８】
図７は、本実施形態に係るマージ予測パラメータ導出部３０３６の構成を示す概略図である。マージ予測パラメータ導出部３０３６は、マージ候補導出部３０３６１とマージ候補選択部３０３６２を備える。マージ候補格納部３０３６１１は、マージ候補導出部３０３６１から入力されたマージ候補を格納する。なお、マージ候補は、予測リスト利用フラグpredFlagLX、ベクトルmvLX、参照ピクチャインデックスrefIdxLXを含んで構成されている。マージ候補格納部３０３６１１において、格納されたマージ候補には、所定の規則に従ってインデックスが割り当てられる。マージ候補導出部３０３６１は、空間マージ候補導出部３０３６１３１と時間マージ候補導出部３０３６１３２と結合マージ候補導出部３０３６１３３とゼロマージ候補導出部３０３６１３４を含んで構成される。
【０２３９】
空間マージ候補導出部３０３６１３１は、所定の規則に従って、予測パラメータメモリ３０７が記憶している予測パラメータ（予測リスト利用フラグpredFlagLX、ベクトルmvLX、参照ピクチャインデックスrefIdxLX）を読み出し、読み出した予測パラメータをマージ候補として導出する。読み出される予測パラメータは、復号対象ブロックから予め定めた範囲内にあるブロック（例えば、復号対象ブロックの左下端、左上端、右上端にそれぞれ接するブロックの全部又は一部）のそれぞれに係る予測パラメータである。導出されたマージ候補はマージ候補格納部３０３６１１に格納される。
【０２４０】
時間マージ候補導出部３０３６１３２は、復号対象ブロックの右下の座標を含む参照画像中のブロックの予測パラメータを予測パラメータメモリ３０７から読みだしマージ候補とする。参照画像の指定方法は、例えば、スライスヘッダに置いて指定された参照ピクチャインデックスrefIdxLXでも良いし、復号対象ブロックに隣接するブロックの参照ピクチャインデックスrefIdxLXのうち最小のものを用いて指定しても良い。導出されたマージ候補はマージ候補格納部３０３６１１に格納される。
【０２４１】
結合マージ候補導出部３０３６１３３は、既に導出されマージ候補格納部３０３６１１に格納された２つの異なる導出済マージ候補のベクトルと参照ピクチャインデックスを、それぞれＬ０、Ｌ１のベクトルとして組み合わせることで結合マージ候補を導出する。導出されたマージ候補はマージ候補格納部３０３６１１に格納される。
【０２４２】
ゼロマージ候補導出部３０３６１３４は、参照ピクチャインデックスrefIdxLXが０であり、ベクトルmvLXのＸ成分、Ｙ成分が共に０であるマージ候補を導出する。導出されたマージ候補はマージ候補格納部３０３６１１に格納される。
【０２４３】
マージ候補選択部３０３６２は、マージ候補格納部３０３６１１に格納されているマージ候補のうち、インター予測パラメータ復号制御部３０３１から入力されたマージインデックスmerge_idxに対応するインデックスが割り当てられたマージ候補を、対象ＰＵのインター予測パラメータとして選択する。マージ候補選択部３０３６２は選択したマージ候補を予測パラメータメモリ３０７（図５）に記憶するとともに、予測画像生成部３０８（図５）に出力する。
【０２４４】
図８は、本実施形態に係るＡＭＶＰ予測パラメータ導出部３０３２の構成を示す概略図である。ＡＭＶＰ予測パラメータ導出部３０３２は、ベクトル候補導出部３０３３と予測ベクトル選択部３０３４を備える。ベクトル候補導出部３０３３は、参照ピクチャインデックスrefIdxに基づいて予測パラメータメモリ３０７（図５）が記憶するベクトル（動きベクトル又は変位ベクトル）をベクトル候補mvpLXとして読み出す。読み出されるベクトルは、復号対象ブロックから予め定めた範囲内にあるブロック（例えば、復号対象ブロックの左下端、左上端、右上端にそれぞれ接するブロックの全部又は一部）のそれぞれに係るベクトルである。
【０２４５】
予測ベクトル選択部３０３４は、ベクトル候補導出部３０３３が読み出したベクトル候補のうち、インター予測パラメータ復号制御部３０３１から入力されたベクトルインデックスmvp_LX_idxが示すベクトル候補を予測ベクトルmvpLXとして選択する。予測ベクトル選択部３０３４は、選択した予測ベクトルmvpLXを加算部３０３５に出力する。
【０２４６】
図９は、ベクトル候補の一例を示す概念図である。図９に示す予測ベクトルリスト６０２は、ベクトル候補導出部３０３３において導出される複数のベクトル候補からなるリストである。予測ベクトルリスト６０２において、左右に一列に配列された５個の長方形は、それぞれ予測ベクトルを示す領域を示す。左端から２番目のmvp_LX_idxの真下の下向きの矢印とその下のmvpLXは、ベクトルインデックスmvp_LX_idxが、予測パラメータメモリ３０７においてベクトルmvpLXを参照するインデックスであることを示す。
【０２４７】
候補ベクトルは、復号処理が完了したブロックであって、復号対象ブロックから予め定めた範囲のブロック（例えば、隣接ブロック）を参照し、参照したブロックに係るベクトルに基づいて生成される。なお、隣接ブロックには、対象ブロックに空間的に隣接するブロック、例えば、左ブロック、上ブロックの他、対象ブロックに時間的に隣接するブロック、例えば、対象ブロックと同じ位置で、表示時刻が異なるブロックから得られたブロックを含む。
【０２４８】
加算部３０３５は、予測ベクトル選択部３０３４から入力された予測ベクトルmvpLXとインター予測パラメータ復号制御部から入力された差分ベクトルmvdLXを加算してベクトルmvLXを算出する。加算部３０３５は、算出したベクトルmvLXを予測画像生成部３０８（図５）に出力する。
【０２４９】
図１０は、インター予測パラメータ復号制御部３０３１の構成を示す図である。インター予測パラメータ復号制御部３０３１は、追加予測フラグ復号部３０３１１、マージインデックス復号部３０３１２、ベクトル候補インデックス復号部３０３１３、及び図示しない、分割モード復号部、マージフラグ復号部、インター予測フラグ復号部、参照ピクチャインデックス復号部、ベクトル差分復号部を含んで構成される。分割モード復号部、マージフラグ復号部、マージインデックス復号部、インター予測フラグ復号部、参照ピクチャインデックス復号部、ベクトル候補インデックス復号部３０３１３、ベクトル差分復号部は各々、分割モードpart_mode、マージフラグmerge_flag、マージインデックスmerge_idx、インター予測フラグinter_pred_idc、参照ピクチャインデックスrefIdxLX、予測ベクトルインデックスmvp_LX_idx、差分ベクトルmvdLXを復号する。
【０２５０】
追加予測フラグ復号部３０３１１は、追加予測である照度予測を行うか否かを示すフラグic_enable_flagを復号し、インター予測画像生成部３０９に出力する。
【０２５１】
変位ベクトル取得部は、対象ＰＵに隣接するブロックが変位ベクトルを持つ場合には、その変位ベクトルを予測パラメータメモリ３０７から抽出し、予測パラメータメモリ３０７を参照し、対象ＰＵに隣接するブロックの予測フラグpredFlagLX、参照ピクチャインデックスrefIdxLXとベクトルmvLXを読み出す。
【０２５２】
（インター予測画像生成部３０９）
図１１は、本実施形態に係るインター予測画像生成部３０９の構成を示す概略図である。インター予測画像生成部３０９は、動き補償部３０９１、照度補償部３０９３、重み予測部３０９４を含んで構成される。
【０２５３】
（動き補償）
動き補償部３０９１は、インター予測パラメータ復号部３０３から入力された、予測リスト利用フラグpredFlagLX、参照ピクチャインデックスrefIdxLX、動きベクトルmvLXに基づいて、参照ピクチャメモリ３０６から、参照ピクチャインデックスrefIdxLXで指定された参照ピクチャの対象ブロックの位置を起点として、ベクトルmvLXだけずれた位置にあるブロックを読み出すことによって動き補償画像を生成する。ここで、ベクトルmvLXが整数ベクトルでない場合には、動き補償フィルタ（もしくは変位補償フィルタ）と呼ばれる小数位置の画素を生成するためのフィルタを施して、動き補償画像を生成する。一般に、ベクトルmvLXが動きベクトルの場合のみ、上記処理を動き補償と呼び、変位ベクトルの場合は変位補償と呼ぶ場合もあるが、本明細書では、ベクトルの種類を区別せず単に動き補償と表現する。以下、Ｌ０予測の動き補償画像をpredSamplesL0、Ｌ１予測の動き補償画像をpredSamplesL1と呼ぶ。両者を区別しない場合predSamplesLXと呼ぶ。以下、動き補償部３０９１で得られた動き変位補償画像predSamplesLXに、さらに照度補償が行われる例を説明するが、これらの出力画像もまた、動き補償画像predSamplesLXと呼ぶ。なお、以下の照度補償において、入力画像と出力画像を区別する場合には、入力画像をpredSamplesLX、出力画像をpredSamplesLX´と表現する。
【０２５４】
（照度補償）
照度補償部３０９３は、照度補償フラグic_enable_flagが１の場合に、入力された動き変位画像predSamplesLXに対して、照度補償を行う。照度補償フラグic_enable_flagが０の場合には、入力された動き補償画像predSamplesLXをそのまま出力する。照度補償部３０９３に入力される動き補償画像predSamplesLXは、動き補償部３０９１の出力画像である。照度補償は、予測画像生成の対象となる対象ブロックに隣接する隣接領域の動き補償画像の画素値と、隣接領域の復号画像の変化が、対象ブロック内の画素値と、対象ブロックの原画との変化に類似しているとの仮定に基づいて行われる。
【０２５５】
照度パラメータ推定部３０９３１は、対象ブロック（対象予測ユニット）の画素を、参照ブロックの画素から推定するための推定パラメータを求める。図１３は、照度補償を説明するための図である。図１３では、対象ブロックの周囲の画素（目的参照画素）と、対象ブロックからあるベクトルだけずれた位置にある、参照レイヤ画像上の参照ブロックの周囲の画素（入力参照画素）の位置を示している。
【０２５６】
照度パラメータ推定部３０９３１は、上記対象ブロックの周囲の目的参照画素y（y０〜yＮ−１）と上記参照ブロックの周囲の入力参照画素x（x０〜xＮ−１）から、推定パラメータ（照度変化パラメータ）を求める。
【０２５７】
（照度補償部３０９３）
図２４は、照度補償部３０９３の構成を示すブロック図である。照度補償部３０９３は、照度パラメータ推定部３０９３１と、照度補償フィルタ部３０９３２から構成される。照度パラメータ推定部３０９３１は、積算シフト値導出部３０９３１１、積算値導出部３０９３１２、加算値導出部３０９３１３、第１パラメータ導出部３０９３１４、第２パラメータ導出部３０９３１５、パラメータａ導出部３０９３１６、パラメータｂ導出部３０９３１７、正則化コスト加算部３０９３１８から構成される。
【０２５８】
照度パラメータ推定部３０９３１は、図１３に示す参照レイヤ画像上の参照ブロックの周囲の画素Ｃを入力参照画素x[]、対象ブロックの周囲の画素Ｌを目的参照画素y[]とし、これら入力参照画素x[]と目的参照画素y[]に基づいて、目的参照画素y[]を入力参照画素xから線形予測する場合のパラメータであるパラメータａ、パラメータｂを導出する。
【０２５９】
加算値導出部３０９３１３は、目的参照画素yの和Yと入力参照画素xの和Xを、以下の式（Ｂ−２）式（Ｂ−３）により導出する。
【０２６０】
積算値導出部３０９３１２は目的参照画素yと入力参照画素xの積の和XYと、入力参照画素の２乗の和XXを、以下の式（Ｂ−４）〜（Ｂ−５）により導出する。この時、積算値導出部３０９３１２は、目的参照画素yと入力参照画素xの積の和XYの導出時と、入力参照画素xの２乗の和XXの導出する。以下の和の前にX、Y、XY、XXは0に初期化しておく。
【０２６１】
X = Σx[i*2]                                  式（Ｂ−２）
Y = Σy[i*2]                                  式（Ｂ−３）
XX = Σ(x[i] * x[i])          式（Ｂ−４）
XY = Σ(x[i] * y[i])          式（Ｂ−５）
ここで、Σは、参照領域に対する和であり、参照領域の画素を特定するインデックスiに対する和を導出する。y[i]は、復号画像のインデックスiにおける画素値。x[i]は、参照画像のインデックスiにおける画素値。カウントシフト値iCountShiftは、参照領域のサイズの２の対数である。
【０２６２】
iCountShift = log2 (参照領域の画素数)           式（Ｂ−６）
第１パラメータ導出部３０９３１４は、目的参照画素yと入力参照画素xの積の和XYと、目的参照画素の和Yと入力参照画素の和Xの積の差から第１パラメータa1を以下の式により導出する。
【０２６３】
a1 = ( XY << iCountShift ) - (Y * X)        式（Ｂ−７）
式（Ｂ−７）に示すように、XYは、カウントシフト値iCountShiftだけ左シフトしてから差を算出する。
【０２６４】
第２パラメータ導出部３０９３１５は、入力参照画素の２乗の和XXと、入力参照画素の和Xの２乗の差から第２パラメータa2を以下の式により導出する。
【０２６５】
a2 = ( XX << iCountShift ) - (X * X);         式（Ｂ−８）
式（Ｂ−８）に示すように、XXは、カウントシフト値iCountShiftだけ左シフトしてから差を算出する。
【０２６６】
導出された第１パラメータa１と第２パラメータa2は、パラメータａ導出部３０９３１６に出力される。
【０２６７】
第１パラメータ導出部３０９３１４は、第１パラメータa1を第２パラメータa2の大きさに応じて制限する。例えば以下の式のように、a1を0以上、a2の2倍以下にクリップする。
【０２６８】
a1 = Clip3(0, 2*a2, a1)          式（Ｂ−１２）
第１パラメータ導出部３０９３１４により、a1 / a2の値は、0から2の間にクリップされる。従って、パラメータａの値であるa1 / a2 << iShiftの値も、0から2<<iShiftにクリップされる。すなわち、iShift=6の場合、パラメータａは0〜128となり、8ビット非負整数の範囲に入る。
【０２６９】
第１パラメータ導出部３０９３１４は、第２パラメータa2の大きさに応じて、図１６のテーブルの導出に用いた所定のビット幅ShiftA2に対して、以下の式により第２正規化シフト値iScaleShiftA2を導出する。導出された第２正規化シフト値iScaleShiftA2は、テーブルベースパラメータａ導出部３０９３１６３に出力される。
【０２７０】
iScaleShiftA2 = Max( 0, Floor( Log2( Abs( a2 ) ) ) - (ShiftA2-1))        式（Ｂ−１４）
第１パラメータ正規化シフト部３０９３１６１は、第２正規化シフト値iScaleShiftA2に応じて、以下の式により第１正規化シフト値iScaleShiftA1を導出する。導出された第１正規化シフト値iScaleShiftA1は、テーブルベースパラメータａ導出部３０９３１６３に出力される。
【０２７１】
iScaleShiftA1 = Max(0, iScaleShiftA2 - offsetA1)        式（Ｂ−１３）
なお、offsetA1はここでは10以下を満たす定数とする。
【０２７２】
上記では、第１正規化シフト値を第２正規化シフト値を用いて導出するため、第１正規化パラメータを導出する処理が容易になるという効果を奏する。
【０２７３】
第１パラメータ導出部３０９３１４、第２パラメータ導出部３０９３１５は、第１パラメータa1を第１正規化シフト値iScaleShiftA1、第２パラメータa2を第２正規化シフト値iScaleShiftA2により右シフトし、正規化第１パラメータa1sと正規化第２パラメータa2sを導出する。
【０２７４】
a1s = a1 >> iScaleShiftA1          式（Ｂ−１５）
a2s = a2 >> iScaleShiftA2          式（Ｂ−１６）
これにより、正規化第１パラメータa1sと正規化第２パラメータa2sは各々、０から2のiShiftA1乗-1、０から2のShiftA2乗-1の間の値に正規化される。
【０２７５】
パラメータａ導出部３０９３１６は、第１正規化シフト値iScaleShiftA1と第２正規化シフト値iScaleShiftA2の差に基づいて、以下の式によりパラメータａシフト値iScaleShiftAを導出する。
【０２７６】
ScaleShiftA  = ShiftA1 + iScaleShiftA2 - iScaleShiftA1 - iShift          式（Ｂ−１８）
ここで、iScaleShiftA1 = Max(0, iScaleShiftA2 - offsetA1)であるから、以下の式が得られる。
【０２７７】
ScaleShiftA  <= ShiftA1 + iScaleShiftA2 - (iScaleShiftA2 - offsetA1) - iShift
ScaleShiftA  <= ShiftA1 + offsetA1 - iShift
offsetA1は0以上、iShiftは5から8ビット、ShiftA1は14ビット〜15ビットであるから、ScaleShiftAは常に0以上になる。
【０２７８】
パラメータａ導出部３０９３１６は、正規化第２パラメータa2sに応じて定まる逆数テーブル値invTableを参照し、正規化第１パラメータa1sとの積をとりテーブルシフト値（ScaleShiftA）で右シフトすることにより、以下の式によりパラメータａを導出する。
【０２７９】
a = (a1s * invTable[a2s]) >> (ScaleShiftA)           式（Ｂ−１９）
図１６は、本実施形態で用いる逆数テーブル値invTable[]である。既に説明したように図１６の逆数invTable[x]は、インデックスxが０の場合に０となり、インデックスxが０以外の場合には、所定の定数（ここでは2の15乗）Mをxで割って整数化した値から導出される。
【０２８０】
導出されたパラメータａは、パラメータｂ導出部３０９３１７と照度補償フィルタ部３０９３２に出力される。
【０２８１】
パラメータｂ導出部３０９３１７は、目的参照画素の和Yから、入力参照画素の和Xにパラメータａをかけて固定シフト値iShiftだけ右シフトした値を引いた値を、参照領域の画素数で割ることにより、以下の式によりパラメータｂを導出する。
【０２８２】
b = ( Y - ( ( a * X ) >> iShift ) + ( 1 << ( iCountShift - 1 ) ) ) >> iCountShift          式（Ｂ−２０）
なお、iCountShiftの右シフトは、参照領域の画素数で割ることに相当する。
【０２８３】
照度補償フィルタ部３０９３２は、照度パラメータ推定部３０９３１が導出した推定パラメータを用いて、照度補償前の予測画像predSamples []から、照度補償後の予測画像predSamples´[]を導出する。例えば、パラメータｂを式（Ｂ−２０）により導出する場合には以下の式を用いる。
【０２８４】
predSamples´[x][y] = (a * predSamples[x][y] >> iShift)+b           式（Ｂ−２１）
（重み予測）
重み予測部３０９４は、入力される動き変位画像predSamplesLXに重み係数を乗算することにより予測ピクチャブロックＰ（予測画像）を生成する。入力される動き変位画像predSamplesLXは、残差予測、照度補償が行われる場合には、それらが施された画像である。参照リスト利用フラグの一方（predFlagL0もしくはpredFlagL1）が１の場合（単予測の場合）で、重み予測を用いない場合には入力された動き変位画像predSamplesLX（ＬＸはＬ０もしくはＬ１）を画素ビット数に合わせる以下の式の処理を行う。
predSamples[ x ][ y ] = Clip3( 0, ( 1 << bitDepth ) - 1, ( predSamplesLX[ x ][ y ] + offset1 ) >> shift1 )
ここで、shift1=14-bitDepth、offset1=1<<(shift1-1)である。
【０２８５】
また、参照リスト利用フラグの両者（predFlagL0もしくはpredFlagL1）が１の場合（双予測の場合）で、重み予測を用いない場合には、入力された動き変位画像predSamplesL0、predSamplesL1を平均し画素ビット数に合わせる以下の式の処理を行う。
predSamples[ x ][ y ] = Clip3( 0, ( 1 << bitDepth ) - 1, ( predSamplesL0[ x ][ y ] + predSamplesL1[ x ][ y ] + offset2 ) >> shift2 )
ここで、shift2=15-bitDepth、offset2=1<<(shift2-1)である。
【０２８６】
さらに、単予測の場合で、重み予測を行う場合には、重み予測部３０９４は、重み予測係数ｗ０とオフセットｏ０を導出し、以下の式の処理を行う。
predSamples[ x ][ y ] = Clip3( 0, ( 1 << bitDepth ) - 1, ( (predSamplesLX[ x ][ y ] * w0 + 2log2WD - 1) >> log2WD ) + o0 )
ここで、ｌｏｇ２ＷＤは所定のシフト量を示す変数である。
【０２８７】
さらに、双予測の場合で、重み予測を行う場合には、重み予測部３０９４は、重み予測係数ｗ０、ｗ１、ｏ０、ｏ１を導出し、以下の式の処理を行う。
predSamples[ x ][ y ] = Clip3( 0, ( 1 << bitDepth ) - 1, ( predSamplesL0 [ x ][ y ] * w0 + predSamplesL1[ x ][ y ] * w1 + ((o0 + o1 + 1) << log2WD) ) >> (log2WD + 1) )
（画像符号化装置の構成）
次に、本実施形態に係る画像符号化装置１１の構成について説明する。図１４は、本実施形態に係る画像符号化装置１１の構成を示すブロック図である。画像符号化装置１１は、予測画像生成部１０１、減算部１０２、ＤＣＴ・量子化部１０３、エントロピー符号化部１０４、逆量子化・逆ＤＣＴ部１０５、加算部１０６、予測パラメータメモリ（予測パラメータ記憶部、フレームメモリ）１０８、参照ピクチャメモリ（参照画像記憶部、フレームメモリ）１０９、符号化パラメータ決定部１１０、予測パラメータ符号化部１１１、残差格納部３１３（残差記録部）を含んで構成される。予測パラメータ符号化部１１１は、インター予測パラメータ符号化部１１２及びイントラ予測パラメータ符号化部１１３を含んで構成される。
【０２８８】
予測画像生成部１０１は、外部から入力されたレイヤ画像Ｔの視点毎の各ピクチャについて、そのピクチャを分割した領域であるブロック毎に予測ピクチャブロックＰを生成する。ここで、予測画像生成部１０１は、予測パラメータ符号化部１１１から入力された予測パラメータに基づいて参照ピクチャメモリ１０９から参照ピクチャブロックを読み出す。予測パラメータ符号化部１１１から入力された予測パラメータとは、例えば、動きベクトル又は変位ベクトルである。予測画像生成部１０１は、符号化対象ブロックを起点として予測された動きベクトル又は変位ベクトルが示す位置にあるブロックの参照ピクチャブロックを読み出す。予測画像生成部１０１は、読み出した参照ピクチャブロックについて複数の予測方式のうちの１つの予測方式を用いて予測ピクチャブロックＰを生成する。予測画像生成部１０１は、生成した予測ピクチャブロックＰを減算部１０２に出力する。なお、予測画像生成部１０１は、既に説明した予測画像生成部３０８と同じ動作であるため予測ピクチャブロックＰの生成の詳細は省略する。
【０２８９】
予測画像生成部１０１は、予測方式を選択するために、例えば、レイヤ画像に含まれるブロックの画素毎の信号値と予測ピクチャブロックＰの対応する画素毎の信号値との差分に基づく誤差値を最小にする予測方式を選択する。予測方式を選択する方法は、これには限られない。
【０２９０】
符号化対象のピクチャが他のレイヤに依存しないベースレイヤピクチャ（ベースビューピクチャ）である場合には、複数の予測方式とは、イントラ予測、動き予測及びマージ予測である。動き予測とは、上述のインター予測のうち、表示時刻間の予測である。マージ予測とは、既に符号化されたブロックであって、符号化対象ブロックから予め定めた範囲内にあるブロックと同一の参照ピクチャブロック及び予測パラメータを用いる予測である。符号化対象のピクチャがノンベースビューピクチャである場合には、複数の予測方式とは、イントラ予測、動き予測、マージ予測、及び変位予測である。変位予測（視差予測）とは、上述のインター予測のうち、別レイヤ画像（別視点画像）間の予測である。さらに、動き予測、マージ予測、及び変位予測である。変位予測（視差予測）に対して、追加予測（照度補償）を行う場合と行わない場合の予測がある。
【０２９１】
予測画像生成部１０１は、イントラ予測を選択した場合、予測ピクチャブロックＰを生成する際に用いたイントラ予測モードを示す予測モードpredModeを予測パラメータ符号化部１１１に出力する。
【０２９２】
予測画像生成部１０１は、動き予測を選択した場合、予測ピクチャブロックＰを生成する際に用いた動きベクトルmvLXを予測パラメータメモリ１０８に記憶し、インター予測パラメータ符号化部１１２に出力する。動きベクトルmvLXは、符号化対象ブロックの位置から予測ピクチャブロックＰを生成する際の参照ピクチャブロックの位置までのベクトルを示す。動きベクトルmvLXを示す情報には、参照ピクチャを示す情報（例えば、参照ピクチャインデックスrefIdxLX、ピクチャ順序番号ＰＯＣ）を含み、予測パラメータを表すものであっても良い。また、予測画像生成部１０１は、インター予測モードを示す予測モードpredModeを予測パラメータ符号化部１１１に出力する。
【０２９３】
予測画像生成部１０１は、変位予測を選択した場合、予測ピクチャブロックＰを生成する際に用いた変位ベクトルを予測パラメータメモリ１０８に記憶し、インター予測パラメータ符号化部１１２に出力する。変位ベクトルdvLXは、符号化対象ブロックの位置から予測ピクチャブロックＰを生成する際の参照ピクチャブロックの位置までのベクトルを示す。変位ベクトルdvLXを示す情報には、参照ピクチャを示す情報（例えば、参照ピクチャインデックスrefIdxLX、ビューＩＤview_id）を含み、予測パラメータを表すものであっても良い。また、予測画像生成部１０１は、インター予測モードを示す予測モードpredModeを予測パラメータ符号化部１１１に出力する。
【０２９４】
予測画像生成部１０１は、マージ予測を選択した場合、選択した参照ピクチャブロックを示すマージインデックスmerge_idxをインター予測パラメータ符号化部１１２に出力する。また、予測画像生成部１０１は、マージ予測モードを示す予測モードpredModeを予測パラメータ符号化部１１１に出力する。
【０２９５】
上記の、動き予測、変位予測、マージ予測において、予測画像生成部１０１は、残差予測フラグres_pred_flagが残差予測を行うことを示す場合には、既に説明したように予測画像生成部１０１に含まれる残差予測部３０９２において残差予測を行い、照度補償フラグic_enable_flagが照度補償を行うことを示す場合には、既に説明したように予測画像生成部１０１に含まれる照度補償部３０９３において照度補償予測を行う。
【０２９６】
減算部１０２は、予測画像生成部１０１から入力された予測ピクチャブロックＰの信号値を、外部から入力されたレイヤ画像Ｔの対応するブロックの信号値から画素毎に減算して、残差信号を生成する。減算部１０２は、生成した残差信号をＤＣＴ・量子化部１０３と符号化パラメータ決定部１１０に出力する。
【０２９７】
ＤＣＴ・量子化部１０３は、減算部１０２から入力された残差信号についてＤＣＴを行い、ＤＣＴ係数を算出する。ＤＣＴ・量子化部１０３は、算出したＤＣＴ係数を量子化して量子化係数を求める。ＤＣＴ・量子化部１０３は、求めた量子化係数をエントロピー符号化部１０４及び逆量子化・逆ＤＣＴ部１０５に出力する。
【０２９８】
エントロピー符号化部１０４には、ＤＣＴ・量子化部１０３から量子化係数が入力され、符号化パラメータ決定部１１０から符号化パラメータが入力される。入力される符号化パラメータには、例えば、参照ピクチャインデックスrefIdxLX、ベクトルインデックスmvp_LX_idx、差分ベクトルmvdLX、予測モードpredMode、及びマージインデックスmerge_idx等の符号がある。
【０２９９】
エントロピー符号化部１０４は、入力された量子化係数と符号化パラメータをエントロピー符号化して符号化ストリームＴｅを生成し、生成した符号化ストリームＴｅを外部に出力する。
【０３００】
逆量子化・逆ＤＣＴ部１０５は、ＤＣＴ・量子化部１０３から入力された量子化係数を逆量子化してＤＣＴ係数を求める。逆量子化・逆ＤＣＴ部１０５は、求めたＤＣＴ係数について逆ＤＣＴを行い、復号残差信号を算出する。逆量子化・逆ＤＣＴ部１０５は、算出した復号残差信号を加算部１０６に出力する。
【０３０１】
加算部１０６は、予測画像生成部１０１から入力された予測ピクチャブロックＰの信号値と逆量子化・逆ＤＣＴ部１０５から入力された復号残差信号の信号値を画素毎に加算して、参照ピクチャブロックを生成する。加算部１０６は、生成した参照ピクチャブロックを参照ピクチャメモリ１０９に記憶する。
【０３０２】
予測パラメータメモリ１０８は、予測パラメータ符号化部１１１が生成した予測パラメータを、符号化対象のピクチャ及びブロック毎に予め定めた位置に記憶する。
【０３０３】
参照ピクチャメモリ１０９は、加算部１０６が生成した参照ピクチャブロックを、符号化対象のピクチャ及びブロック毎に予め定めた位置に記憶する。
【０３０４】
符号化パラメータ決定部１１０は、符号化パラメータの複数のセットのうち、１つのセットを選択する。符号化パラメータとは、上述した予測パラメータやこの予測パラメータに関連して生成される符号化の対象となるパラメータである。予測画像生成部１０１は、これらの符号化パラメータのセットの各々を用いて予測ピクチャブロックＰを生成する。
【０３０５】
符号化パラメータ決定部１１０は、複数のセットの各々について情報量の大きさと符号化誤差を示すコスト値を算出する。コスト値は、例えば、符号量と二乗誤差に係数λを乗じた値との和である。符号量は、量子化誤差と符号化パラメータをエントロピー符号化して得られる符号化ストリームＴｅの情報量である。二乗誤差は、減算部１０２において算出された残差信号の残差値の二乗値についての画素間の総和である。係数λは、予め設定されたゼロよりも大きい実数である。符号化パラメータ決定部１１０は、算出したコスト値が最小となる符号化パラメータのセットを選択する。これにより、エントロピー符号化部１０４は、選択した符号化パラメータのセットを符号化ストリームＴｅとして外部に出力し、選択されなかった符号化パラメータのセットを出力しない。
【０３０６】
予測パラメータ符号化部１１１は、予測画像生成部１０１から入力されたパラメータに基づいて予測ピクチャを生成する際に用いる予測パラメータを導出し、導出した予測パラメータを符号化して符号化パラメータのセットを生成する。予測パラメータ符号化部１１１は、生成した符号化パラメータのセットをエントロピー符号化部１０４に出力する。
【０３０７】
予測パラメータ符号化部１１１は、生成した符号化パラメータのセットのうち符号化パラメータ決定部１１０が選択したものに対応する予測パラメータを予測パラメータメモリ１０８に記憶する。
【０３０８】
予測パラメータ符号化部１１１は、予測画像生成部１０１から入力された予測モードpredModeがインター予測モードを示す場合、インター予測パラメータ符号化部１１２を動作させる。予測パラメータ符号化部１１１は、予測モードpredModeがイントラ予測モードを示す場合、イントラ予測パラメータ符号化部１１３を動作させる。
【０３０９】
インター予測パラメータ符号化部１１２は、符号化パラメータ決定部１１０から入力された予測パラメータに基づいてインター予測パラメータを導出する。インター予測パラメータ符号化部１１２は、インター予測パラメータを導出する構成として、インター予測パラメータ復号部３０３（図５等、参照）がインター予測パラメータを導出する構成と同一の構成を含む。インター予測パラメータ符号化部１１２の構成については、後述する。
【０３１０】
イントラ予測パラメータ符号化部１１３は、符号化パラメータ決定部１１０から入力された予測モードpredModeが示すイントラ予測モードIntraPredModeをインター予測パラメータのセットとして定める。
【０３１１】
（ＬＭ予測部）
予測画像生成部１０１は、予測画像生成部３０８と同様にＬＭ予測部３１０４を備える。ＬＭ予測部３１０４は、既に説明したように、ＬＭパラメータ推定部３１０４１を備え、インデックスｉに対応する画素値xiと画素値yiの組を入力として、xiからyiを予測する線形予測パラメータａとｂを導出する線形予測パラメータ導出部と、上記線形予測パラメータに基づいて予測画像を生成する線形予測部を備える画像予測装置であって、上記線形予測パラメータ導出部は、画素値xiと画素値yiの積の和XYと画素値xiの和Xと画素値yiの和Yの積から導出される第１パラメータa1と、画素値xiと画素値xiの積の和XXと画素値xiの和Xと画素値xiの和Xの積から導出される第２パラメータa2から線形予測パラメータａを導出する手段を備え、上記線形予測パラメータ導出部は、上記第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN未満もしくは上記閾値THN以下である場合には第１コストa1costNを上記第１パラメータa1から減算する手段と、上記第１パラメータa1と所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THP未満もしくは上記閾値THP以下の場合には第２コストa1costPを上記第１パラメータa1に加算する手段の少なくとも１つの手段を備えることを特徴とする。
【０３１２】
本構成の、予測画像生成部１０１の含むＬＭパラメータ推定部３１０４１は、例えばa1がTHN未満（ここではTHN=0）である場合に、a1から正則化コストa1costNだけ減算するように閾値との比較に基づいて正則化コストを減算もしくは加算することで、パラメータaの値を、最小２乗法で導出される値よりも０から遠い方向に近づけることが可能である。上記、構成により、ＬＭパラメータ推定部３１０４１は、雑音に強い線形予測パラメータを導出することができ、予測画像を生成する際の予測精度を向上させる効果を奏する。
【０３１３】
（インター予測パラメータ符号化部の構成）
次に、インター予測パラメータ符号化部１１２の構成について説明する。インター予測パラメータ符号化部１１２は、インター予測パラメータ復号部３０３に対応する手段である。
【０３１４】
図１５は、本実施形態に係るインター予測パラメータ符号化部１１２の構成を示す概略図である。
【０３１５】
インター予測パラメータ符号化部１１２は、図示しないインター予測パラメータ符号化制御部１０３１、マージ予測パラメータ導出部３０３６、ＡＭＶＰ予測パラメータ導出部３０３２、減算部１１２３、及び予測パラメータ統合部１１２６を含んで構成される。
【０３１６】
マージ予測パラメータ導出部３０３６は、上述のマージ予測パラメータ導出部３０３６（図７参照）と同様な構成を有する。
【０３１７】
インター予測パラメータ符号化制御部１０３１は、インター予測に関連する符号（シンタックス要素の復号をエントロピー符号化部１０４に指示し、符号化データに含まれる符号（シンタックス要素）を例えば、分割モードpart_mode、マージフラグmerge_flag、マージインデックスmerge_idx、インター予測フラグinter_pred_idc、参照ピクチャインデックスrefIdxLX、予測ベクトルインデックスmvp_LX_idx、差分ベクトルmvdLXを符号化する。
【０３１８】
インター予測パラメータ符号化制御部１０３１は、追加予測フラグ符号化部１０３１１、マージインデックス符号化部１０３１２、ベクトル候補インデックス符号化部１０３１３、及び図示しない、分割モード符号化部、マージフラグ符号化部、インター予測フラグ符号化部、参照ピクチャインデックス符号化部、ベクトル差分符号化部を含んで構成される。分割モード符号化部、マージフラグ符号化部、マージインデックス符号化部、インター予測フラグ符号化部、参照ピクチャインデックス符号化部、ベクトル候補インデックス符号化部１０３１３、ベクトル差分符号化部は各々、分割モードpart_mode、マージフラグmerge_flag、マージインデックスmerge_idx、インター予測フラグinter_pred_idc、参照ピクチャインデックスrefIdxLX、予測ベクトルインデックスmvp_LX_idx、差分ベクトルmvdLXを符号化する。
【０３１９】
追加予測フラグ符号化部１０３１１は、追加予測が行われるか否かを示すために、照度補償フラグic_enable_flagを符号化する。
【０３２０】
マージ予測パラメータ導出部３０３６には、予測画像生成部１０１から入力された予測モードpredModeがマージ予測モードを示す場合、符号化パラメータ決定部１１０からマージインデックスmerge_idxが入力される。マージインデックスmerge_idxは、予測パラメータ統合部１１２６に出力される。マージ予測パラメータ導出部３０３６は、マージ候補のうちマージインデックスmerge_idxが示す参照ブロックの参照ピクチャインデックスrefIdxLX、ベクトルmvLXを予測パラメータメモリ１０８から読み出す。マージ候補とは、符号化対象となる符号化対象ブロックから予め定めた範囲にある参照ブロック（例えば、符号化対象ブロックの左下端、左上端、右上端に接する参照ブロックのうち）であって、符号化処理が完了した参照ブロックである。
【０３２１】
ＡＭＶＰ予測パラメータ導出部３０３２は、上述のＡＭＶＰ予測パラメータ導出部３０３２（図８参照）と同様な構成を有する。
【０３２２】
ＡＭＶＰ予測パラメータ導出部３０３２には、予測画像生成部１０１から入力された予測モードpredModeがインター予測モードを示す場合、符号化パラメータ決定部１１０からベクトルmvLXが入力される。ＡＭＶＰ予測パラメータ導出部３０３２は、入力されたベクトルmvLXに基づいて予測ベクトルmvpLXを導出する。ＡＭＶＰ予測パラメータ導出部３０３２は、導出した予測ベクトルmvpLXを減算部１１２３に出力する。なお、参照ピクチャインデックスrefIdx及びベクトルインデックスmvp_LX_idxは、予測パラメータ統合部１１２６に出力される。
【０３２３】
減算部１１２３は、符号化パラメータ決定部１１０から入力されたベクトルmvLXからＡＭＶＰ予測パラメータ導出部３０３２から入力された予測ベクトルmvpLXを減算して差分ベクトルmvdLXを生成する。差分ベクトルmvdLXは予測パラメータ統合部１１２６に出力する。
【０３２４】
予測画像生成部１０１から入力された予測モードpredModeがマージ予測モードを示す場合には、予測パラメータ統合部１１２６は、符号化パラメータ決定部１１０から入力されたマージインデックスmerge_idxをエントロピー符号化部１０４に出力する。
【０３２５】
予測画像生成部１０１から入力された予測モードpredModeがインター予測モードを示す場合には、予測パラメータ統合部１１２６は、次の処理を行う。
【０３２６】
予測パラメータ統合部１１２６は、符号化パラメータ決定部１１０から入力された参照ピクチャインデックスrefIdxLX及びベクトルインデックスmvp_LX_idx、減算部１１２３から入力された差分ベクトルmvdLXを統合する。予測パラメータ統合部１１２６は、統合した符号をエントロピー符号化部１０４に出力する。
【０３２７】
なお、上述した実施形態における画像符号化装置１１、画像復号装置３１の一部、例えば、エントロピー復号部３０１、予測パラメータ復号部３０２、予測画像生成部１０１、ＤＣＴ・量子化部１０３、エントロピー符号化部１０４、逆量子化・逆ＤＣＴ部１０５、符号化パラメータ決定部１１０、予測パラメータ符号化部１１１、エントロピー復号部３０１、予測パラメータ復号部３０２、予測画像生成部３０８、逆量子化・逆ＤＣＴ部３１１をコンピュータで実現するようにしても良い。その場合、この制御機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することによって実現しても良い。なお、ここでいう「コンピュータシステム」とは、画像符号化装置１１−１１ｈ、画像復号装置３１−３１ｈのいずれかに内蔵されたコンピュータシステムであって、ＯＳや周辺機器等のハードウェアを含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含んでも良い。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
【０３２８】
また、上述した実施形態における画像符号化装置１１、画像復号装置３１の一部、または全部を、ＬＳＩ（Ｌａｒｇｅ  Ｓｃａｌｅ  Ｉｎｔｅｇｒａｔｉｏｎ）等の集積回路として実現しても良い。画像符号化装置１１、画像復号装置３１の各機能ブロックは個別にプロセッサ化しても良いし、一部、または全部を集積してプロセッサ化しても良い。また、集積回路化の手法はＬＳＩに限らず専用回路、または汎用プロセッサで実現しても良い。また、半導体技術の進歩によりＬＳＩに代替する集積回路化の技術が出現した場合、当該技術による集積回路を用いても良い。
【０３２９】
以上、図面を参照してこの発明の一実施形態について詳しく説明してきたが、具体的な構成は上述のものに限られることはなく、この発明の要旨を逸脱しない範囲内において様々な設計変更等をすることが可能である。
【０３３０】
本発明の一つの形態は、インデックスｉに対応する画素値xiと画素値yiの組を入力として、xiからyiを予測する線形予測パラメータａとｂを導出する線形予測パラメータ導出部と、上記線形予測パラメータに基づいて予測画像を生成する線形予測部を備える画像予測装置であって、上記線形予測パラメータ導出部は、画素値xiと画素値yiの積の和XYと画素値xiの和Xと画素値yiの和Yの積に基づいて導出される第１パラメータa1と、画素値xiと画素値xiの積の和XXと画素値xiの和Xと画素値xiの和Xの積に基づいて導出される第２パラメータa2から線形予測パラメータａを導出する手段を備え、上記線形予測パラメータ導出部は、上記第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN未満もしくは上記閾値THN以下である場合には第１コストa1costNを上記第１パラメータa1から減算する手段と、上記第１パラメータa1と所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THP未満もしくは上記閾値THP以下の場合には第２コストa1costPを上記第１パラメータa1に加算する手段の少なくとも１つの手段を備えることを特徴とする。
【０３３１】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THNより小さい場合、もしくは、上記第１パラメータa1が上記所定の閾値THN以下の場合に、上記第１パラメータa1に上記第１コストa1costNを減算することを特徴とする。
【０３３２】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THPより大きい場合、もしくは、上記第１パラメータa1が上記所定の閾値THP以上の場合に、上記第１パラメータa1に上記第２コストa1costPを加算することを特徴とする。
【０３３３】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THNより小さい場合、もしくは、上記第１パラメータa1が上記所定の閾値THN以下の場合に、上記第１パラメータa1に上記第１コストa1costNを減算し、上記第１パラメータa1が上記所定の閾値THPより大きい場合、もしくは、上記第１パラメータa1が上記所定の閾値THP以上の場合に、上記第１パラメータa1に上記第２コストa1costPを加算することを特徴とする。
【０３３４】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第２パラメータa2にコストa2costを加算する手段をさらに備えることを特徴とする。
【０３３５】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１パラメータa1と所定の閾値THNを比較し、上記第１パラメータa1が上記閾値THN未満もしくは上記閾値THN以下である場合には第１コストa1costNを上記第１パラメータa1から減算し、上記第２パラメータa2にコストa2costNを加算する手段と、上記第１パラメータa1と所定の閾値THPを比較し、上記第１パラメータa1が上記閾値THP未満もしくは上記閾値THP以下の場合には第２コストa1costPを上記第１パラメータa1に加算し上記第２パラメータa2にコストa2costPを加算する手段の少なくとも１つの手段を備えることを特徴とする。
【０３３６】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THNより小さい場合、もしくは、上記第１パラメータa1が上記所定の閾値THN以下の場合に、上記第１パラメータa1に上記第１コストa1costNを減算し、上記第２パラメータa2にコストa2costNを加算する手段を備えることを特徴とする。
【０３３７】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１パラメータa1が上記所定の閾値THPより大きい場合、もしくは、上記第１パラメータa1が上記所定の閾値THP以上の場合に、上記第１パラメータa1に上記第２コストa1costPを加算し、上記第２パラメータa2にコストa2costNを加算する手段を備えることを特徴とする。
【０３３８】
本発明の一つの形態は、上記コストa2costは、上記コストa1costN以上であることを特徴とする。
【０３３９】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記閾値THNとTHPのどちらか一方もしくは両方を０とすることを特徴とする。
【０３４０】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記閾値THNとTHPのどちらか一方もしくは両方を、上記第２パラメータa2に基づいて導出することを特徴とする。
【０３４１】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記閾値THNを０、上記閾値THPを上記第２パラメータa2に基づいて導出することを特徴とする。
【０３４２】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記閾値THNを、上記第２パラメータa2を所定の定数だけ右シフトすることにより導出する、ことを特徴とする。
【０３４３】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記画素値xiと画素値xiの積の和XXから導出することを特徴とする。
【０３４４】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記画素値xiと画素値yiの積の和XYから導出することを特徴とする。
【０３４５】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、画素のビットデプス値から導出することを特徴とする。
【０３４６】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記第２パラメータa2に基づいて導出することを特徴とする。
【０３４７】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記第２パラメータa2を所定の定数だけ右シフトすることにより、導出することを特徴とする。
【０３４８】
本発明の一つの形態は、上記線形予測パラメータ導出部は、上記第１コストa1costNもしくは上記第２コストa1costPを、上記第２パラメータa1を所定の定数だけ右シフトした値と、画素のビットデプス値から導出される値との最小値から、導出することを特徴とする。
（関連出願の相互参照）
本出願は、2015年10月5日に出願された日本国特許出願：特願2015-197489に対して優先権の利益を主張するものであり、それを参照することにより、その内容の全てが本書に含まれる。
【産業上の利用可能性】
【０３４９】
本発明は、画像データが符号化された符号化データを復号する画像復号装置、および、画像データが符号化された符号化データを生成する画像符号化装置、および、色コンポーネント間予測を行う予測画像生成装置に好適に適用することができる。
【符号の説明】
【０３５０】
１…画像伝送システム
１１…画像符号化装置
１０１…予測画像生成部
１０２…減算部
１０３…ＤＣＴ・量子化部
１０４…エントロピー符号化部
１０５…逆量子化・逆ＤＣＴ部
１０６…加算部
１０８…予測パラメータメモリ（フレームメモリ）
１０９…参照ピクチャメモリ（フレームメモリ）
１１０…符号化パラメータ決定部
１１１…予測パラメータ符号化部
１１２…インター予測パラメータ符号化部
３０３６…マージ予測パラメータ導出部
３０３２…ＡＭＶＰ予測パラメータ導出部
１１２３…減算部
１１２６…予測パラメータ統合部
１１３…イントラ予測パラメータ符号化部
２１…ネットワーク
３１…画像復号装置
３０１…エントロピー復号部
３０２…予測パラメータ復号部
３０３…インター予測パラメータ復号部
３０３１２…マージインデックス復号部
３０３１３…ベクトル候補インデックス復号部
３０３２…ＡＭＶＰ予測パラメータ導出部
３０３５…加算部
３０３６…マージ予測パラメータ導出部
３０３６１…マージ候補導出部
３０３６１１…マージ候補格納部
３０３６１３１…空間マージ候補導出部
３０３６１３２…時間マージ候補導出部
３０３６１３３…結合マージ候補導出部
３０３６１３４…ゼロマージ候補導出部
３０３６２…マージ候補選択部
３０４…イントラ予測パラメータ復号部
３０６…参照ピクチャメモリ（フレームメモリ）
３０７…予測パラメータメモリ（フレームメモリ）
３０８…予測画像生成部
３０９…インター予測画像生成部
３０９１…動き補償部
３０９２…残差予測部
３０９３…照度補償部
３０９３１…照度パラメータ推定部
３０９３１２…積算値導出部
３０９３１３…加算値導出部
３０９３１４…第１パラメータ導出部
３０９３１５…第２パラメータ導出部
３０９３１６…パラメータａ導出部
３０９３１６１…第１パラメータ正規化シフト部
３０９３１６２…第２パラメータ正規化シフト部
３０９３１６３…テーブルベースパラメータａ導出部
３０９３１７…パラメータｂ導出部
３０９３１８…正則化コスト加算部
３０９３２…照度補償フィルタ部
３０９４…重み予測部
３１０…イントラ予測画像生成部
３１０４…ＬＭ予測部
３１０４１、３１０４１Ａ、３１０４１Ｂ、３１０４１Ｃ、３１０４１Ａ２、３１０４１Ｂ２、３１０４１Ｃ２…ＬＭパラメータ推定部
３１０４２、３１０４２Ｂ…ＬＭフィルタ部
３１０４１２…ＬＭ積算値導出部
３１０４１３…ＬＭ加算値導出部
３１０４１４…ＬＭ第１パラメータ導出部
３１０４１５…ＬＭ第２パラメータ導出部
３１０４１６…ＬＭパラメータａ導出部
３１０４１６１…ＬＭ第１パラメータ正規化シフト部
３１０４１６２…ＬＭ第２パラメータ正規化シフト部
３１０４１６３…ＬＭ除算値導出部
３１０４１７…ＬＭパラメータｂ導出部
３１０４１８、３１０４１８Ａ、３１０４１８Ｂ、３１０４１８Ｃ、３１０４１８Ａ２、３１０４１８Ｂ２、３１０４１８Ｃ２、…ＬＭ正則化コスト加算部
３１０４１８０…ＬＭ正則化コスト導出部
３１０４１８１…ＬＭ第１パラメータ正則化コスト加算部
３１０４１８２…ＬＭ第２パラメータ正則化コスト加算部
３１０４１８３…ＬＭ正則化コスト閾値判定部
３１１…逆量子化・逆ＤＣＴ部
３１２…加算部
４１…画像表示装置
【図１】
2017061189.tif 000003
【図２】
2017061189.tif 000004
【図３】
2017061189.tif 000005
【図４】
2017061189.tif 000006
【図５】
2017061189.tif 000007
【図６】
2017061189.tif 000008
【図７】
2017061189.tif 000009
【図８】
2017061189.tif 000010
【図９】
2017061189.tif 000011
【図１０】
2017061189.tif 000012
【図１１】
2017061189.tif 000013
【図１２】
2017061189.tif 000014
【図１３】
2017061189.tif 000015
【図１４】
2017061189.tif 000016
【図１５】
2017061189.tif 000017
【図１６】
2017061189.tif 000018
【図１７】
2017061189.tif 000019
【図１８】
2017061189.tif 000020
【図１９】
2017061189.tif 000021
【図２０】
2017061189.tif 000022
【図２１】
2017061189.tif 000023
【図２２】
2017061189.tif 000024
【図２３】
2017061189.tif 000025
【図２４】
2017061189.tif 000026
【図２５】
2017061189.tif 000027
【国際調査報告】
2017061189.tif 000028
2017061189.tif 000029
2017061189.tif 000030
2017061189.tif 000031