配列ディザリング

配列ディザリング（はいれつでぃざりんぐ）とは、画像をディザリングするアルゴリズムの一つである。

このアルゴリズムは通常、色深度が小さくて表示可能な色数が少ないディスプレイにおいて、なるべく滑らかな画像を表示するために使用される。たとえば、Microsoft Windowsは「VGA標準16色グラフィックモード」の時に配列ディザリング法を用いてグラフィックを表示する。画像にこのアルゴリズムを適用した場合、特徴的なクロスハッチング状のタイルパターンがハッキリと現れる。

しきい値マップ

このアルゴリズムは「しきい値マップ（ $M$ ）」を使って色数を減らす。ディスプレイに表示される画像の各ピクセルの色の値と、対応する「しきい値マップ」のセルの値を比較し、画面に表示する色を決める。その過程で、元画像の色と、そのディスプレイで利用可能なカラーパレットの色との差分が切り捨てられる。

「しきい値マップ」の大きさは任意に決めて構わないが、普通は「2の累乗かける 2の累乗」のマトリクスを使う。ちなみに、この「しきい値マップ」のことを「テーブル」という日本人がいるが、「マトリクス」が正しい（「テーブル」は「値」で構成された表、「マトリックス」は「行」と「列」で構成された表、つまり「行列」のこと）。

\mathbf {M} ={\frac {1}{4}}\times {\begin{bmatrix}0&2\\3&1\\\end{bmatrix}}

\mathbf {M} ={\frac {1}{16}}\times {\begin{bmatrix}0&8&2&10\\12&4&14&6\\3&11&1&9\\15&7&13&5\\\end{bmatrix}}

\mathbf {M} ={\frac {1}{64}}\times {\begin{bmatrix}0&48&12&60&3&51&15&63\\32&16&44&28&35&19&47&31\\8&56&4&52&11&59&7&55\\40&24&36&20&43&27&39&23\\2&50&14&62&1&49&13&61\\34&18&46&30&33&17&45&29\\10&58&6&54&9&57&5&53\\42&26&38&22&41&25&37&21\\\end{bmatrix}}

このマップは、回転や反転を行ったとしてもアルゴリズムの有効性に影響を与えない。両辺の長さがどちらも2の累乗になった、このようなしきい値マップは「インデックスマトリクス」または「ベイヤーマトリクス」とも呼ばれる^[1]。白黒2値で写真を表現するためにコダック社のブライス・ベイヤーが考案し、1973年に発表した。（ベイヤー以後にもディザリング用に独自のマトリクスを考案した者は存在するが、結局ベイヤーマトリクスが一番早くて奇麗なのでずっと使われている。ランダムな値のマトリクスを使う「ランダムディザリング法」も存在するが、綺麗にうまくいくシーンはそれほど多くない）

配列ディザリング法を使う場合は、上記のマトリクスをそのまま使うのが一般的だが、もし任意のサイズのしきい値マップを自力で作りたい場合は、簡単な法則を使って制作できる。まず、各スロットを連続した整数で埋める。次に、マップ内の2つの連続した数同士の平均距離ができるだけ大きくなるように並べ替え、テーブルが端で「折り返される」ようにする^[要出典]。

次元が2の累乗であるしきい値マップの場合、マップは以下の定式で再帰的に生成できる。

\mathbf {M} _{2n}={\frac {1}{(2n)^{2}}}\times {\begin{bmatrix}4\times \mathbf {M} _{n}&4\times \mathbf {M} _{n}+2\\4\times \mathbf {M} _{n}+3&4\times \mathbf {M} _{n}+1\end{bmatrix}}

再帰式は、ビット演算のみを使用して明示的に計算できる^[2]。

M(i, j) = bit_reverse(bit_interleave(bitwise_xor(i, j), i)) / n ^ 2

前処理したしきい値マップ

ディザリングが実行されるハードウェアの精度にもよるが、しきい値マップを、 0 から $n^{2}$ までの整数を $n$ x $n$ のマトリクスとして保持しておくという伝統的な方式よりも、前処理で浮動小数点の形式にしたものを保持していた方が便利な場合がある。

その場合、以下の式が使われる。

Mpre(i,j) = (Mint(i,j)+1) / n^2

これにより、「標準化（正規化）しきい値マトリクス」が生成される。

例えば2×2のしきい値マップは、

{\frac {1}{4}}\times {\begin{bmatrix}0&2\\3&1\\\end{bmatrix}}

であるが、これをベースとして「しきい値マップ」に前処理をかけると、

{\begin{bmatrix}0.25&0.75\\1.0&0.5\\\end{bmatrix}}

となる。加えて言うと、前処理によって各値を正規化して合計を0に平均化するようなマップを作ることもできる。この場合、しきい値マップの全ての値から0.5を減算すればよい。

Mpre(i,j) = (Mint(i,j)+1) / n^2 - 0.5

これにより、次のような前処理済ディザマップが生成される。

{\begin{bmatrix}-0.25&0.25\\0.5&0\\\end{bmatrix}}

アルゴリズム

任天堂ポケットカメラ（1998年発売）。配列ディザリング法は軽くて十分きれいなので、Z80カスタム（4MHz）を搭載した携帯型ゲーム機でも動画（縦112ドット×横128ドット、モノクロ4階調）のリアルタイム処理が可能

配列ディザリングのアルゴリズムは画像を自然に表現できるが、しかし画像のピクセルの各々においては、そのピクセルの色の値がしきい値マップの対応するセルの値で置き換えられることになるので、ピクセルの色の値のしきい値を超えた部分がカットされ、元の色とは異なる色に量子化されることになる。

配列ディザリング法のアルゴリズムを組むに際して、ほとんどの場合、すべてのピクセルにしきい値を追加して、ピクセルの値をしきい値と比較するだけで充分ディザリングの実用を満たすので、わざわざ0.5を減算して「しきい値マップ」を正規化する必要まではない。元画像のピクセルの輝度の値がベイヤーマトリクスの対応するセルの値よりも小さかった場合は、そのピクセルを「黒」で描写し、そうじゃなかった場合は「白」で描写する、みたいなアルゴリズムで充分である。ただ、正規化しなかった場合、画像の平均輝度がわずかに増加し、「ほとんど白に近いグレー」のピクセルがディザリングされずにそのまま白になる、という欠点がある。もっともこれは、グレースケールパレット（もしくは色同士の相対的な距離が（ほぼ）一定なカラーパレット）を使用する場合は問題とはならない。人間の目は明色の階調差よりも暗色の階調差の方をより正確に認識するため、「ほとんど白に近いグレー」を白飛びさせて暗色の階調差に割り振るのは、むしろ望ましいディザリング結果になるとすら言える。一方で、色数の特に少ないカラーパレットや任意の色で構成されたカラーパレットを使用してディザリングを行った場合は、元画像と全然違った見た目の画像が生成されることになるため、アルゴリズムを組む人は面倒がらずに適切な正規化を行った方が良い。

配列ディザリング法で、64色のカラーパレット（左）を使って、140 x 140 = 19600色のカラースケールを表現したもの（右）

逆の目線で言うと、アルゴリズムは、すべてのピクセルの各カラーチャンネル「 $c$ 」に対して次の変換を実行する。

c'=\mathrm {nearest\_palette\_color} {\mathopen {}}\left(c+r\times \left(M(x{\bmod {n}},y{\bmod {n}})-{\frac {1}{2}}\right){\mathclose {}}\right)

ここで、「 $M (i, j)$ 」は $i$ 番目の行と $j$ 番目の列のしきい値マップ、「 $c'$ 」は変換された色、「 $r$ 」は色空間における広がりの量である。RGBの各色が0から255までのオクテットで表された、全「 $2 3 N$ 」色の等間隔な色で構成されるカラーパレットを想定すると、一つの色は通常はこのような形で選択される。

r\approx {\frac {255}{N}}

このアルゴリズムは1つ１つのピクセルのみを操作するもので、条件分岐もないため、非常に高速で、リアルタイム変換に適している。さらに、ディザパターンは表示フレームに対して常に同じ位置にあるため、誤差拡散法よりもジッター（いわゆる「カクツキ」「チラツキ」、現代のゲームエンジンのレトロゲーム機風シェーダーでも問題になる）が発生しにくく、アニメーションに適している。誤差拡散法のディザパターンよりも反復的なタイリングであるため、配列ディザリングを使用した画像の圧縮率は高くなる。配列ディザリングは、直線が真っ直ぐに変換され、揺らぎが少ないため、イラストなど線画系の画像に適している（写真などは誤差拡散法の方が奇麗に出る）。

カラーパレットに合わせてしきい値マップを自作する場合、しきい値マップの値同士の差は、カラーパレットの色同士の最小の差と同じ範囲になるようにスケーリングすることが望ましい。

注釈

^ Bayer, Bryce (June 11–13, 1973). “An optimum method for two-level rendition of continuous-tone pictures”. IEEE International Conference on Communications 1: 11–15. オリジナルの2013-05-12時点におけるアーカイブ。.
^ Joel Yliluoma. “Arbitrary-palette positional dithering algorithm”

参照

Ordered Dithering (Graphics course project, Visgraf lab, Brazil)
Dithering algorithms (Lee Daniel Crocker, Paul Boulay and Mike Morra)

外部リンク

Matlab implementation of various dithering methods

[bayermatrix-1] Bayer, Bryce (June 11–13, 1973). “An optimum method for two-level rendition of continuous-tone pictures”. IEEE International Conference on Communications 1: 11–15. オリジナルの2013-05-12時点におけるアーカイブ。.

[2] Joel Yliluoma. “Arbitrary-palette positional dithering algorithm”

[1]

[2]

しきい値マップ

前処理したしきい値マップ

アルゴリズム

注釈

参照

関連項目

外部リンク