特許 6348746 画像圧縮回路および画像圧縮方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6348746

(24)【登録日】2018年6月8日

(45)【発行日】2018年6月27日

(54)【発明の名称】画像圧縮回路および画像圧縮方法

(51)【国際特許分類】

   H04N 1/41 20060101AFI20180618BHJP

   H04N 19/124 20140101ALI20180618BHJP

   H04N 19/176 20140101ALI20180618BHJP

   H04N 19/154 20140101ALI20180618BHJP

【ＦＩ】

   H04N1/41

   H04N19/124

   H04N19/176

   H04N19/154

【請求項の数】8

【全頁数】14

(21)【出願番号】特願2014-69332(P2014-69332)

(22)【出願日】2014年3月28日

(65)【公開番号】特開2015-192372(P2015-192372A)

(43)【公開日】2015年11月2日

【審査請求日】2017年2月17日

(73)【特許権者】

【識別番号】591128453

【氏名又は名称】株式会社メガチップス

(74)【代理人】

【識別番号】100080159

【弁理士】

【氏名又は名称】渡辺 望稔

(74)【代理人】

【識別番号】100090217

【弁理士】

【氏名又は名称】三和 晴子

(72)【発明者】

【氏名】渡辺 祥則

(72)【発明者】

【氏名】橋本 順之

【審査官】 松永 隆志

(56)【参考文献】

【文献】 特表２０１０−５０３２４５（ＪＰ，Ａ）

【文献】 特表２００７−５３４２３９（ＪＰ，Ａ）

【文献】 特開平０３−１６７９８６（ＪＰ，Ａ）

【文献】 下城 拓也 ほか３名，ヒストグラムの疎性を考慮したＬｌｏｙｄ−Ｍａｘ法の効果的適用法，電子情報通信学会技術研究報告，２０１１年 ６月２３日，第111巻、第105号， p.1-6

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｈ０４Ｎ １／００− １／６４

Ｈ０４Ｎ １９／００−１９／９８

(57)【特許請求の範囲】

【請求項1】

入力されたＲＧＢ画像データをＹＵＶフォーマットに変換する画像データ変換回路と、
前記画像データ変換回路により前記ＹＵＶフォーマットに変換された画像データをそれぞれ所定の画素数からなる複数の単位領域に分割し、単位領域毎にＹ成分およびＵＶ成分の２系統に分けて、それぞれ、複数の仮想代表値を設定すると共に前記複数の仮想代表値から複数の第１の量子化候補値を算出し、各単位領域内の画素毎に前記複数の第１の量子化候補値の中から前記画像データの値にもっとも近い第１の量子化候補値を選択することで仮想量子化データを取得する第１の量子化回路と、
Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、前記第１の量子化回路により設定された前記複数の仮想代表値および前記第１の量子化回路により取得された前記仮想量子化データに基づいて前記単位領域毎にデータ圧縮の誤差を最小にする複数の最適代表値を最小二乗法により演算する最適化回路と、
Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、前記最適化回路で得られた前記複数の最適代表値に基づいて前記単位領域毎に量子化に用いられる複数の第２の量子化候補値を算出し、各単位領域内の画素毎に前記複数の第２の量子化候補値の中から前記画像データの値にもっとも近い第２の量子化候補値を選択することで最適量子化データを取得する第２の量子化回路と、
Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、前記最適化回路で演算された前記複数の最適代表値および前記第２の量子化回路で得られた前記最適量子化データを所定のデータセット構造にパッキングして出力するパッキング回路と
を備えたことを特徴とする画像圧縮回路。

【請求項2】

前記第１の量子化回路は、各単位領域内の画素毎に選択された前記第１の量子化候補値を前記複数の仮想代表値に対する係数の行列で表したものを前記仮想量子化データとし、
前記第２の量子化回路は、各単位領域内の画素毎に選択された前記第２の量子化候補値を前記複数の最適代表値に対する係数の行列で表したものを前記最適量子化データとする請求項１に記載の画像圧縮回路。

【請求項3】

前記第１の量子化回路は、前記単位領域毎に画素値の最大値および最小値を前記複数の仮想代表値として設定する請求項１または２に記載の画像圧縮回路。

【請求項4】

前記第１の量子化回路は、前記複数の仮想代表値と、前記複数の仮想代表値の間の補間値を、各単位領域に対する前記複数の第１の量子化候補値として算出する請求項３に記載の画像圧縮回路。

【請求項5】

前記最適化回路は、ヤコビ法を用いて前記複数の最適代表値を演算する請求項１〜４のいずれか一項に記載の画像圧縮回路。

【請求項6】

前記最適化回路は、整数倍された前記画像データおよび前記仮想量子化データを用いた漸化計算を複数回繰り返すことにより前記複数の最適代表値を演算する請求項５に記載の画像圧縮回路。

【請求項7】

前記第２の量子化回路は、前記複数の最適代表値と、前記複数の最適代表値の間の補間値を、各単位領域に対する前記複数の第２の量子化候補値として算出する請求項１〜６のいずれか一項に記載の画像圧縮回路。

【請求項8】

入力されたＲＧＢ画像データをＹＵＶフォーマットに変換し、
前記ＹＵＶフォーマットに変換された画像データをそれぞれ所定の画素数からなる複数の単位領域に分割し、単位領域毎にＹ成分およびＵＶ成分の２系統に分けて、それぞれ、複数の仮想代表値を設定すると共に前記複数の仮想代表値から複数の第１の量子化候補値を算出し、各単位領域内の画素毎に前記複数の第１の量子化候補値の中から前記画像データの値にもっとも近い第１の量子化候補値を選択することで仮想量子化データを取得し、
Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、設定された前記複数の仮想代表値および取得された前記仮想量子化データに基づいて前記単位領域毎にデータ圧縮の誤差を最小にする複数の最適代表値を最小二乗法により演算し、
Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、前記複数の最適代表値に基づいて前記単位領域毎に量子化に用いられる複数の第２の量子化候補値を算出し、各単位領域内の画素毎に前記複数の第２の量子化候補値の中から前記画像データの値にもっとも近い第２の量子化候補値を選択することで最適量子化データを取得し、
Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、前記複数の最適代表値および前記最適量子化データを所定のデータセット構造にパッキングして出力する
ことを特徴とする画像圧縮方法。

【発明の詳細な説明】

【技術分野】

【0001】

この発明は、画像圧縮回路および画像圧縮方法に係り、特に、テクスチャ圧縮により画像データを圧縮する画像圧縮回路および方法に関する。

【背景技術】

【0002】

例えば、特許文献１には、圧縮処理の中間で用いられる輝度変更子と圧縮処理前の輝度との間の誤差を低減するように輝度変更子および量子化色を選択する画像処理方法が開示されている。
また、特許文献２には、液晶ディスプレイのソースドライバで用いられる画像圧縮において、所定範囲の画素ブロックを輝度および色差の変化に合致した複数のサブ領域に分割し、サブ領域毎に処理前の色との間の誤差を小さくするように代表色を決定する画像処理方法が開示されている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特許第４４４４９６７号公報

【特許文献2】特開２００９−１３０６９０号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

しかしながら、特許文献１および２に開示された画像処理方法では、画像データの圧縮前後における誤差の低減はなされるものの、圧縮処理に先立って、輝度変更子および量子化色、あるいは、代表色を決定するのみであるので、誤差の低減は十分なものとは言い難く、圧縮による画質の劣化を伴うものであった。

【0005】

この発明は、このような従来の問題点を解消するためになされたもので、圧縮による画質の劣化を最小限に抑えることができる画像圧縮回路および画像圧縮方法を提供することを目的とする。

【課題を解決するための手段】

【0006】

この発明に係る画像圧縮回路は、入力されたＲＧＢ画像データをＹＵＶフォーマットに変換する画像データ変換回路と、画像データ変換回路によりＹＵＶフォーマットに変換された画像データをそれぞれ所定の画素数からなる複数の単位領域に分割し、単位領域毎にＹ成分およびＵＶ成分の２系統に分けて、それぞれ、複数の仮想代表値を設定すると共に複数の仮想代表値から複数の第１の量子化候補値を算出し、各単位領域内の画素毎に複数の第１の量子化候補値の中から画像データの値にもっとも近い第１の量子化候補値を選択することで仮想量子化データを取得する第１の量子化回路と、Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、第１の量子化回路により設定された複数の仮想代表値および第１の量子化回路により取得された仮想量子化データに基づいて単位領域毎にデータ圧縮の誤差を最小にする複数の最適代表値を最小二乗法により演算する最適化回路と、Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、最適化回路で得られた複数の最適代表値に基づいて単位領域毎に量子化に用いられる複数の第２の量子化候補値を算出し、各単位領域内の画素毎に複数の第２の量子化候補値の中から画像データの値にもっとも近い第２の量子化候補値を選択することで最適量子化データを取得する第２の量子化回路と、Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、最適化回路で演算された複数の最適代表値および第２の量子化回路で得られた最適量子化データを所定のデータセット構造にパッキングして出力するパッキング回路とを備えたものである。

【0007】

第１の量子化回路は、各単位領域内の画素毎に選択された第１の量子化候補値を複数の仮想代表値に対する係数の行列で表したものを仮想量子化データとし、第２の量子化回路は、各単位領域内の画素毎に選択された第２の量子化候補値を複数の最適代表値に対する係数の行列で表したものを最適量子化データとすることができる。
第１の量子化回路は、単位領域毎に画素値の最大値および最小値を複数の仮想代表値として設定することが好ましい。
この場合、第１の量子化回路は、複数の仮想代表値と、複数の仮想代表値の間の補間値を、各単位領域に対する複数の第１の量子化候補値として算出することができる。
最適化回路は、ヤコビ法を用いて複数の最適代表値を演算することが好ましい。
この場合、最適化回路は、整数倍された画像データおよび仮想量子化データを用いた漸化計算を複数回繰り返すことにより複数の最適代表値を演算することができる。

【0008】

また、第２の量子化回路は、複数の最適代表値と、複数の最適代表値の間の補間値を、各単位領域に対する複数の第２の量子化候補値として算出することができる。
なお、入力されたＲＧＢデータをＹＵＶフォーマットに変換し、画像データとして第１の量子化回路に出力する画像データ変換回路をさらに備えてもよい。
また、最適化回路で演算されたＹ成分の複数の最適代表値および第２の量子化回路で得られたＹ成分の最適量子化データと最適化回路で演算されたＵＶ成分の複数の最適代表値および第２の量子化回路で得られたＵＶ成分の最適量子化データを所定のデータセット構造にパッキングして出力するパッキング回路をさらに備えることもできる。

【0009】

この発明に係る画像圧縮方法は、入力されたＲＧＢ画像データをＹＵＶフォーマットに変換し、ＹＵＶフォーマットに変換された画像データをそれぞれ所定の画素数からなる複数の単位領域に分割し、単位領域毎にＹ成分およびＵＶ成分の２系統に分けて、それぞれ、複数の仮想代表値を設定すると共に複数の仮想代表値から複数の第１の量子化候補値を算出し、各単位領域内の画素毎に複数の第１の量子化候補値の中から画像データの値にもっとも近い第１の量子化候補値を選択することで仮想量子化データを取得し、Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、設定された複数の仮想代表値および取得された仮想量子化データに基づいて単位領域毎にデータ圧縮の誤差を最小にする複数の最適代表値を最小二乗法により演算し、Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、複数の最適代表値に基づいて単位領域毎に量子化に用いられる複数の第２の量子化候補値を算出し、各単位領域内の画素毎に複数の第２の量子化候補値の中から画像データの値にもっとも近い第２の量子化候補値を選択することで最適量子化データを取得し、Ｙ成分およびＵＶ成分の２系統に分けて、それぞれ、複数の最適代表値および最適量子化データを所定のデータセット構造にパッキングして出力する方法である。

【発明の効果】

【0010】

この発明によれば、複数の仮想代表値に基づいて仮想量子化データを取得した後、データ圧縮の誤差を最小にする複数の最適代表値を最小二乗法により演算し、演算された複数の最適代表値に基づいて最適量子化データを取得するので、圧縮による画質の劣化を最小限に抑えることが可能となる。

【図面の簡単な説明】

【0011】

【図1】この発明の実施の形態に係る画像圧縮回路の構成を示すブロック図である。

【図2】１つの単位領域における圧縮前のＹデータの一例を示す図である。

【図3】図２に示したＹデータを３倍したＹ成分の整数倍データを示す図である。

【図4】図３に示したデータに対して画素毎に選択された第１の量子化候補値を示す図である。

【図5】Ｙ成分の仮想量子化データを示す図である。

【図6】単位領域における画素の並びを示す図である。

【図7】ＵＶ成分の仮想量子化データを示す図である。

【図8】実施の形態における最適化回路の内部構成を示す機能ブロック図である。

【図9】圧縮されたＹ成分データを示す図である。

【図10】パッキング回路でパッキングされた圧縮データを示す図である。

【発明を実施するための形態】

【0012】

以下、この発明の実施の形態を添付図面に基づいて説明する。
図１に、実施の形態に係る画像圧縮回路の構成を示す。画像圧縮回路は、画像データ変換回路１を有し、この画像データ変換回路１に第１の量子化回路２が接続され、第１の量子化回路２に最適化回路３が接続されている。さらに、最適化回路３に第２の量子化回路４が接続され、第２の量子化回路４にパッキング回路５が接続されている。

【0013】

画像データ変換回路１は、圧縮処理の対象となるＲＧＢ画像データをＹＵＶフォーマットに変換し、変換されたＹＵＶ画像データを第１の量子化回路２に出力するものであり、汎用の変換回路を使用することができる。
一般に、画像圧縮による誤差のうち、もっとも人間の目に感じやすい誤差が輝度変化であるため、画像データ変換回路１によりＲＧＢ画像データをＹＵＶフォーマットに変換すれば、輝度を表すＹ成分のデータとそれ以外の色差に関するＵＶ成分のデータとに分けて圧縮処理を行うことで、圧縮率を維持しながら、特に輝度に関して圧縮誤差を低減することができ、人間の目から見た画質の劣化を効果的に抑制することが可能となる。

【0014】

第１の量子化回路２は、画像データ変換回路１から入力された画像データを、それぞれ所定の画素数からなる複数の単位領域に分割し、単位領域毎に複数の仮想代表値を設定して量子化を行うことで仮想量子化データを取得する回路である。具体的には、複数の仮想代表値から複数の第１の量子化候補値が算出され、各単位領域内の画素毎に複数の第１の量子化候補値の中から画像データの値にもっとも近い第１の量子化候補値を選択することで仮想量子化データが取得される。
最適化回路３は、第１の量子化回路２により設定された複数の仮想代表値および第１の量子化回路２により取得された仮想量子化データに基づき、最小二乗法により、単位領域毎にデータ圧縮の誤差を最小にする複数の最適代表値を演算する。

【0015】

第２の量子化回路４は、最適化回路３で演算された複数の最適代表値に基づいて、量子化をし直すことで最適量子化データを取得する回路である。具体的には、複数の最適代表値から複数の第２の量子化候補値が算出され、各単位領域内の画素毎に複数の第２の量子化候補値の中から画像データの値にもっとも近い第２の量子化候補値を選択することで最適量子化データが取得される。
パッキング回路５は、最適化回路３で演算された複数の最適代表値と第２の量子化回路４で得られた最適量子化データに基づき、圧縮データを所定のデータセット構造にパッキングして出力する。

【0016】

次に、具体例を挙げて、実施の形態に係る画像圧縮回路の動作について説明する。
まず、ＹＵＶデータのうち、Ｙデータについて、処理を行う場合を説明する。
Ｙデータは、画素毎に８ビットを有するデータであり、圧縮を行う単位領域を４画素×４画素の領域とする。
１つの単位領域における圧縮前のＹデータの一例を図２に示す。
第１の量子化回路２は、この単位領域内のＹデータから仮想代表値ｘ_０として画素値の最小値ｍｉｎＹと最大値ｍａｘＹを抽出する。すなわち、図２の例では、ｍｉｎＹ＝０、ｍａｘＹ＝１００となる。

【0017】

次に、最小値ｍｉｎＹおよび最大値ｍａｘＹと、これら最小値ｍｉｎＹおよび最大値ｍａｘＹの間を線形に補間した２つの補間値（２＊ｍｉｎＹ＋１＊ｍａｘＹ）／３および（１＊ｍｉｎＹ＋２＊ｍａｘＹ）／３を用いて、単位領域毎にＹデータに対する４つの第１の量子化候補値ｃ０〜ｃ３を次式（１）のように設定する。
ｃ０＝３＊minＹ
ｃ１＝２＊minＹ＋１＊maxＹ
ｃ２＝１＊minＹ＋２＊maxＹ
ｃ３＝３＊maxＹ ・・・（１）

【0018】

最小値ｍｉｎＹおよび最大値ｍａｘＹの間を線形に補間することで、圧縮誤差の低減を図っている。
なお、上記の式（１）の右辺は、圧縮前のＹデータの値を３倍したものに対応している。すなわち、式（１）の右辺を「３」で除算することで、圧縮前のＹデータの値に対応するものとなるが、計算における除算の導入による回路の煩雑化を回避するため、「３」で除算せずに、圧縮前のＹデータの値を３倍した整数倍データｂを用いて計算を進めるものとする。

【0019】

Ｙデータに対する第１の量子化候補値ｃ０〜ｃ３は、式（１）における最小値ｍｉｎＹおよび最大値ｍａｘＹの係数ｃｏｅｆＹを抜き出して以下の式（２）のような形で、第１の量子化回路２内に予め保持された係数ＬＵＴ（ルックアップテーブル）に設定することができる。
(coefＹ) for ｃ０＝(３，０)
(coefＹ) for ｃ１＝(２，１)
(coefＹ) for ｃ２＝(１，３)
(coefＹ) for ｃ３＝(０，３) ・・・（２）

【0020】

図２に示したＹデータでは、ｍｉｎＹ＝０、ｍａｘＹ＝１００であるので、
ｃ０＝０
ｃ１＝１００
ｃ２＝２００
ｃ３＝３００ ・・・（３）
となる。

【0021】

これらの第１の量子化候補値ｃ０〜ｃ３に対し、図２に示したＹデータを量子化する手順を以下に示す。
上述したように、仮想代表値ｘ_０および補間値を３倍した値で第１の量子化候補値ｃ０〜ｃ３を設定したため、図２に示したＹデータをそれぞれ３倍して図３に示すＹ成分の整数倍データｂを作成し、図４に示されるように、単位領域内の画素毎に第１の量子化候補値ｃ０〜ｃ３のうち、整数倍データｂの値にもっとも近い値を有する第１の量子化候補値を選択する。さらに、選択された第１の量子化候補値に、上記の式（２）に示した係数ｃｏｅｆＹを割り当てると、図５に示されるような画素毎のＹ成分の仮想量子化データが得られる。

【0022】

このようにして、仮想代表値ｘ_０に基づく仮想の量子化処理が行われる。処理結果として得られた画素毎の係数で表されるＹ成分の仮想量子化データは、図６に示されるような画素Ｐ１〜Ｐ１６の並びに従って、画素Ｐ１から画素Ｐ１６まで羅列した数値データとして、以下のように１６行×２列の行列Ａ_０として定義することができる。

【0023】

【数1】

【0024】

ここで、Ｙデータの最小値ｍｉｎＹおよび最大値ｍａｘＹからなる２値の仮想代表値を、２×１型のベクトルｘ_０として再定義すると共に、図３に示したように、圧縮前のＹデータの値を３倍した整数倍データを、１６×１型のベクトルｂとして再定義すると、第１の量子化回路２における仮想の量子化処理は、
Ａ_０・ｘ_０＝ｂ ・・・（４）
にもっとも近くなるように、行列Ａ_０を選定する作業に他ならない。

【0025】

このようなＹデータに対する量子化処理と同様にして、ＹＵＶデータのうち、ＵＶデータについても量子化の処理が行われる。
圧縮前のＵデータおよびＶデータも、Ｙデータと同様に、それぞれ画素毎に８ビットを有するデータであり、圧縮を行う単位領域を４画素×４画素の領域とする。

【0026】

第１の量子化回路２は、圧縮対象の単位領域内のＵデータおよびＶデータから仮想代表値ｘ_０としてそれぞれの最小値ｍｉｎＵおよびｍｉｎＶと最大値ｍａｘＵおよびｍａｘＶを抽出し、ＵおよびＶを軸とした仮想２次元平面上で、これら最小値および最大値を４点で補間することにより、単位領域毎にＵＶデータに対する４つの第１の量子化候補値ｃ０〜ｃ３を次式（５）のように設定する。
ｃ０＝(２＊minＵ，２＊minＶ)
ｃ１＝(minＵ＋maxＵ−minＶ＋maxＶ，minＵ−maxＵ＋minＶ＋maxＶ)
ｃ２＝(−minＵ＋maxＵ＋minＶ＋maxＶ，minＵ＋maxＵ＋minＶ−maxＶ)
ｃ３＝(２＊maxＵ，２＊maxＶ) ・・・（５）

【0027】

なお、上記の式（５）の右辺は、圧縮前のＵデータおよびＶデータの値を２倍したものに対応している。すなわち、式（５）の右辺を「２」で除算することで、圧縮前のＵデータおよびＶデータの値に対応するものとなるが、計算における除算の導入による回路の煩雑化を回避するため、「２」で除算せずに、圧縮前のＵデータおよびＶデータの値を２倍した整数倍データｂを用いて計算を進めるものとする。

【0028】

ＵＶデータに対する第１の量子化候補値ｃ０〜ｃ３は、式（５）におけるＵデータの最小値ｍｉｎＵおよび最大値ｍａｘＵの係数ｃｏｅｆＵとＶデータの最小値ｍｉｎＶおよび最大値ｍａｘＶの係数ｃｏｅｆＶを抜き出して以下の式（６）のような形で、第１の量子化回路２内の係数ＬＵＴに設定することができる。
(coefＵ，coefＶ) for ｃ０＝(２，０，０，０)，(０，０，２，０)
(coefＵ，coefＶ) for ｃ１＝(１，１，-１，１)，(１，-１，１，１)
(coefＵ，coefＶ) for ｃ２＝(-１，１，１，１)，(１，１，１，-１)
(coefＵ，coefＶ) for ｃ３＝(０，２，０，０)，(０，０，０，２) ・・・（６）

【0029】

これらの第１の量子化候補値ｃ０〜ｃ３を用いて、Ｙデータの際と同様に、ＵデータおよびＶデータの量子化が行われる。仮想代表値ｘ_０および補間値を２倍した値で第１の量子化候補値ｃ０〜ｃ３を設定したため、圧縮前の画素毎のＵデータおよびＶデータの値をそれぞれ２倍してＵＶ成分の整数倍データｂを生成し、整数倍データｂの値に対して、第１の量子化候補値ｃ０〜ｃ３のうちもっとも近い値を有する第１の量子化候補値を選択する。選択された第１の量子化候補値に、上記の式（６）に示した係数ｃｏｅｆＵおよびｃｏｅｆＶを割り当てると、例えば図７に示されるような画素毎のＵＶ成分の仮想量子化データが得られる。

【0030】

この場合、画素Ｐ１のＵ、画素Ｐ１のＶ、画素Ｐ２のＵ、画素Ｐ２のＶ・・・というように、画素毎にＵ成分の仮想量子化データとＶ成分の仮想量子化データを並べながら画素Ｐ１から画素Ｐ１６まで各画素の仮想量子化データを羅列すると、以下のように３２行×４列の行列Ａ_０として定義することができる。

【0031】

【数2】

【0032】

ここで、Ｕデータの最小値ｍｉｎＵおよび最大値ｍａｘＵとＶデータの最小値ｍｉｎＶおよび最大値ｍａｘＶからなる４値の仮想代表値を、４×１型のベクトルｘ_０として再定義すると共に、圧縮前のＵデータおよびＶデータの値を２倍した整数倍データを、画素Ｐ１のＵ、画素Ｐ１のＶ、画素Ｐ２のＵ、画素Ｐ２のＶ・・・の順で並べた３２×１型のベクトルｂとして再定義すると、Ｙデータの際と同様に、ＵデータおよびＶデータに対しても、第１の量子化回路２における仮想の量子化処理は、上記の式（４）にもっとも近くなるように、行列Ａ_０を選定する作業となる。

【0033】

このようにして演算された各画素の仮想量子化データＡ_０と仮想代表値ｘ_０と整数倍データｂは、圧縮による誤差のない理想的な状態においては、上記の式（４）の関係を満たすこととなるが、実際には、量子化誤差を含むため、式（４）の関係を満たなくなる。
そこで、
Ａ_０・ｘ−ｂ・・・（７）
を最小にするような代表値ｘを最適代表値ｘ_１として求め、この最適代表値ｘ_１を用いて改めて量子化処理を実施することにより、理想的な状態に近づけることができる。

【0034】

単位領域毎にデータ圧縮の誤差を最小にする最適代表値ｘ_１が、最適化回路３により演算される。
上記の式（７）を最小にすることは、Ｃ^ＴがＣの転置行列を表すものとして、
(Ａ_０・ｘ−ｂ)＊(Ａ_０・ｘ−ｂ)^Ｔ ・・・（８）
を最小にすることと同等であり、この式（８）をベクトルの各成分に関して微分したものが０になるという条件と等価となる。すなわち、
２Ａ_０^ＴＡ_０＊ｘ−２Ａ_０^Ｔｂ＝０ ・・・（９）
となり、さらに、式（９）を変形して、
Ａ_０^ＴＡ_０＊ｘ＝Ａ_０^Ｔｂ ・・・（１０）
となる。

【0035】

最適化回路３は、ヤコビ（Jacobi）法を用いて最小二乗法により式（１０）を数値的に解くことで、代表値ｘの最適化を行い、最適代表値ｘ_１を演算する。
ここで、行列Ｂ、ベクトルｚおよびｙが、
Ｂｚ＝ｙ ・・・（１１）
の関係を満たすとき、これをｚについて解いた場合に、ｋ番目のｚ_ｋと（ｋ＋１）番目のｚ_ｋ＋１に対し、以下の漸化式として表されることが知られている。
ｚ_ｋ＋１＝Ｂ_ｄｉａｇ^−１(ｙ−Ｂ_{ｎｏｎｄｉａｇ}ｚ_ｋ) ・・・（１２）
ただし、Ｂ_ｄｉａｇ^−１はＢの対角成分の逆行列、Ｂ_{ｎｏｎｄｉａｇ}はＢの非対角成分の行列を表す。

【0036】

式（１０）と式（１１）を比較すると、Ａ_０^ＴＡ_０がＢに、ｘがｚに、Ａ_０^Ｔｂがｙに相当する。従って、Ａ_０^ＴＡ_０の対角成分の逆行列および非対角成分の行列並びにＡ_０^Ｔｂの値を計算し、式（１２）の形に当てはめ、ｘを繰り返し算出することで、ヤコビ法により式（１０）をｘについて漸化的に解くことが可能となる。
最適化回路３は、図８に示されるように、直列に接続された第１の計算部３１、第２の計算部３２および第３の計算部３３を有して、漸化計算をパイプライン処理する構造を備えている。第１の計算部３１、第２の計算部３２および第３の計算部３３は、それぞれ、同一のＡ_０を用いて、式（１２）による漸化計算を１回実行する回路ブロックである。

【0037】

第１の計算部３１は、第１の量子化回路２で得られた画素毎の仮想量子化データの行列Ａ_０によるＡ_０^ＴＡ_０およびＡ_０^Ｔｂと、第１の量子化回路２で設定された仮想代表値ｘ_０とに基づいて、式（１２）による漸化計算を１回行って新たな代表値ｘａを算出する。
第２の計算部３２は、第１の量子化回路２で得られた画素毎の仮想量子化データの行列Ａ_０によるＡ_０^ＴＡ_０およびＡ_０^Ｔｂと、第１の計算部３１で算出された代表値ｘａとに基づいて、式（１２）による漸化計算をさらに１回行って新たな代表値ｘｂを算出する。
第３の計算部３３は、第１の量子化回路２で得られた画素毎の仮想量子化データの行列Ａ_０によるＡ_０^ＴＡ_０およびＡ_０^Ｔｂと、第２の計算部３２で算出された代表値ｘｂとに基づいて、式（１２）による漸化計算をさらに１回行い、算出された代表値を最適化回路３で得られた最適代表値ｘ_１として第２の量子化回路４へ出力する。
このようなパイプライン処理により算出された代表値は、仮想代表値ｘ_０から、ｘａおよびｘｂを経て最適代表値ｘ_１に至るほど、次第に精度を増して最適化されることとなる。

【0038】

なお、この実施の形態においては、Ｙデータに対する仮想量子化データの行列Ａ_０は、１６行×２列の行列であり、転置行列Ａ_０^Ｔは、２行×１６列の行列となるため、Ａ_０^ＴＡ_０は、２行×２列の行列となる。さらに、整数倍データｂは１６×１型のベクトルであるので、Ａ_０^Ｔｂは、２行×１列の行列となる。
また、ＵデータおよびＶデータに対する仮想量子化データの行列Ａ_０は、３２行×４列の行列であり、転置行列Ａ_０^Ｔは、４行×３２列の行列となるため、Ａ_０^ＴＡ_０は、４行×４列の行列となる。さらに、整数倍データｂは３２×１型のベクトルであるので、Ａ_０^Ｔｂは、４行×１列の行列となる。

【0039】

このようにして最適化回路３で演算された最適代表値ｘ_１が第２の量子化回路４に入力されると、第２の量子化回路４は、最適代表値ｘ_１を用いて改めてＹデータの量子化を実施すると共に、最適代表値ｘ_１を用いて改めてＵデータおよびＶデータの量子化を実施することでＹ成分の最適量子化データおよびＵＶ成分の最適量子化データを演算する。
第２の量子化回路４における量子化処理は、仮想代表値ｘ_０の代わりに最適代表値ｘ_１を使用する点を除いて、第１の量子化回路２で行われた仮想の量子化処理と同様である。すなわち、最適代表値ｘ_１に基づいて複数の第２の量子化候補値ｃ０〜ｃ３が算出され、各単位領域内の画素毎に複数の第２の量子化候補値ｃ０〜ｃ３の中から画像データの値にもっとも近い第２の量子化候補値を選択することで最適量子化データが取得される。

【0040】

取得された最適量子化データは、４つの第２の量子化候補値ｃ０〜ｃ３に対応する２ビットのデータで表すことができる。例えば、図９に示されるように、圧縮されたＹ成分データは、画素毎に、第２の量子化候補値ｃ０〜ｃ３にそれぞれ対応するデータ「００」、「０１」、「１０」および「１１」のいずれかとなる。
この実施の形態においては、ＹＵＶのデータを、Ｙ成分とＵＶ成分の２系統に分けてそれぞれ圧縮を行っており、双方とも同数の第２の量子化候補値ｃ０〜ｃ３が用いられるので、ＵＶ成分についても、図９に示したものと同様に、画素毎に２ビットの圧縮されたＵＶ成分データが生成される。

【0041】

このようにして第２の量子化回路４で得られた最適量子化データは、最適化回路３で演算された最適代表値ｘ_１と共にパッキング回路５に入力され、所定のデータセット構造を有する圧縮データにパッキングされる。例えば、図１０に示されるように、単位領域毎に、Ｙ成分は、最適代表値ｘ_１である最小値ｍｉｎＹおよび最大値ｍａｘＹと画素毎の圧縮されたＹ成分データＹ１〜Ｙ１６から構成され、ＵＶ成分は、最適代表値ｘ_１である最小値ｍｉｎＵおよびｍｉｎＶと最大値ｍａｘＵおよびｍａｘＶと画素毎の圧縮されたＵＶ成分データＵＶ１〜ＵＶ１６から構成されるデータセット構造となる。

【0042】

この実施の形態においては、最適代表値ｘ_１を構成するＹデータの最小値および最大値、Ｕデータの最小値および最大値、Ｖデータの最小値および最大値は、それぞれ８ビットのデータであるので、Ｙ成分とＵＶ成分を併せて４８ビットとなる。また、圧縮された成分データＹ１〜Ｙ１６およびＵＶ１〜ＵＶ１６は、それぞれ２ビットのデータであるので、併せて６４ビットとなる。すなわち、１つの単位領域の画像データが、合計１１２ビットの圧縮データで表されることとなる。

【0043】

これに対して、圧縮前のデータは、それぞれ８ビットのＹデータ、ＵデータおよびＶデータを有する画素が、単位領域内に４×４＝１６個存在するため、合計８×３×１６＝３８４ビットであり、この実施の形態により、大幅なデータ圧縮がなされたことがわかる。さらに、最適化回路３で代表値の最適化が行われているため、圧縮の誤差は最小となり、画質の劣化が最小限に抑制されたものとなっている。

【0044】

以上のようにして圧縮され、パッキングされた圧縮データは、この圧縮データを受信した回路で伸張することにより、画像データに戻すことができる。
このとき、まず、Ｙデータに関して、最適代表値ｘ_１である最小値ｍｉｎＹおよび最大値ｍａｘＹから、次の式（１３）により量子化候補値ｃ０〜ｃ３を算出する。
ｃ０＝minＹ
ｃ１＝(２＊minＹ＋１＊maxＹ)/３
ｃ２＝(１＊minＹ＋２＊maxＹ)/３
ｃ３＝maxＹ ・・・（１３）
ここでは、回路処理の煩雑化を回避するために省略した、右辺を「３」で除算する処理を行っている。
次に、各画素の圧縮されたＹ成分データに基づき、それぞれ該当する量子化候補値ｃ０〜ｃ３を参照して、式（１３）によりデータの値を変換することで、画像データへの伸張を行うことができる。

【0045】

ＵデータおよびＶデータについても、同様に、最適代表値ｘ_１である最小値ｍｉｎＵおよびｍｉｎＶと最大値ｍａｘＵおよびｍａｘＶから、次の式（１４）により量子化候補値ｃ０〜ｃ３を算出する。
ｃ０＝(minＵ，minＶ)
ｃ１＝((maxＵ＋minＵ＋maxＶ−minＶ)/２，(−maxＵ＋minＵ＋maxＶ＋minＶ)/２)
ｃ２＝((maxＵ−minＵ＋maxＶ＋minＶ)/２，(maxＵ＋minＵ−maxＶ＋minＶ)/２)
ｃ３＝(maxＵ，maxＶ)
・・・（１４）
ここでは、回路処理の煩雑化を回避するために省略した、右辺を「２」で除算する処理を行っている。
次に、各画素の圧縮されたＵＶ成分データに基づき、それぞれ該当する量子化候補値ｃ０〜ｃ３を参照して、式（１４）によりデータの値を変換することで、画像データへの伸張を行うことができる。

【0046】

図１０に示したデータ構造は、通常のテクスチャ圧縮で使用される汎用の構造であるため、圧縮データを受信した回路で伸張する場合も、特別な処理を必要とすることなく、汎用の処理回路で取り扱うことが可能となる。
以上説明したように、第１の量子化回路２で仮想代表値ｘ_０に基づいて仮想量子化データを取得した後、最適化回路３でデータ圧縮の誤差を最小にする最適代表値ｘ_１を最小二乗法により演算し、演算された複数の最適代表値ｘ_１に基づいて第２の量子化回路４で最適量子化データを取得し、さらに、パッキング回路５が、最適量子化データに基づいて圧縮データを所定のデータセット構造にパッキングするので、圧縮による画質の劣化を最小限に抑えながら、画像データの圧縮を行うことが可能となる。

【0047】

また、上記の実施の形態では、最適化回路３がヤコビ法を用いて漸化計算を効率的にパイプライン処理しているので、計算処理を繰り返しているにもかかわらず、処理速度の低下が回避されている。このため、動画処理、液晶画面のパネルセルフリフレッシュ（ＰＳＲ）機能等、画像データの処理に所定の速度が要求される用途において、特に有用な画像圧縮回路を提供することができる。
なお、最適化回路３では、第１の計算部３１、第２の計算部３２および第３の計算部３３により、漸化計算を３回行ったが、これに限るものではなく、処理結果の精度と処理に必要な時間に対応して、２回以下の漸化計算あるいは４回以上の漸化計算を行ってもよい。

【0048】

また、上記の実施の形態では、式（１）および式（５）により算出された第１の量子化候補値および第２の量子化候補値を用いたが、これに限るものではなく、他の演算式により、これら第１の量子化候補値および第２の量子化候補値を算出することもできる。
上記の実施の形態においては、圧縮を行う単位領域を４画素×４画素の領域としたが、これに限るものではなく、任意の画素数を有する単位領域を設定して圧縮処理を行うことができる。

【0049】

さらに、上記の実施の形態では、ＲＧＢ画像データをＹＵＶフォーマットに変換し、輝度を表すＹ成分のデータとそれ以外の色差に関するＵＶ成分のデータの２系統に分けて圧縮処理を行ったが、これに限るものではない。例えば、画像データをＹ成分とＵ成分とＶ成分の３系統に分けて処理してもよく、また、ＹＵＶのデータを複数系統に分けずに一度に三次元近似することもできる。さらに、ＹＵＶフォーマットに変換することなく、ＲＧＢ画像データのまま、圧縮処理を行うこともできる。

【符号の説明】

【0050】

１ 画像データ変換回路、２ 第１の量子化回路、３ 最適化回路、４ 第２の量子化回路、５ パッキング回路、３１ 第１の計算部、３２ 第２の計算部、３３ 第３の計算部、ｘ_０ 仮想代表値、ｘａ，ｘｂ 代表値、ｘ_１ 最適代表値、ｃ０〜ｃ３ 第１の量子化候補値，第２の量子化候補値、ｂ 整数倍データ、Ａ_０ 仮想量子化データ、Ａ_１ 最適量子化データ、Ｐ１〜Ｐ１６ 画素、ｍｉｎＹ，ｍｉｎＵ，ｍｉｎＶ 最小値、ｍａｘＹ，ｍａｘＵ，ｍａｘＶ 最大値、Ｙ１〜Ｙ１６ 圧縮されたＹ成分データ、ＵＶ１〜ＵＶ１６ 圧縮されたＵＶ成分データ。

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】