プリワーピング（Pre-warping）：信号処理・制御工学からVR/CGまで

1. プリワーピングとは何か（定義・概要）

「プリワーピング（Pre-warping）」とは、後段の処理で生じる歪みをあらかじめ逆方向に補正しておく技術の総称である。分野によって具体的な意味は異なるが、共通する本質は「予測される歪みに対する事前補正（pre-compensation）」である。

最も広く知られているのは、デジタル信号処理・制御工学における双一次変換（Bilinear Transform）の周波数プリワーピングである。双一次変換はアナログフィルタ（連続時間系）をデジタルフィルタ（離散時間系）に変換する際に用いられる標準的手法だが、変換過程で周波数軸に非線形な歪み（frequency warping）が生じる。この歪みを設計段階で事前に補正する手法がプリワーピングと呼ばれる。

また、コンピュータグラフィックス・VR分野では、レンズ歪みや3D画像ワーピングに対する事前補正処理をプリワーピングと呼ぶ。

2. どの分野で使われるか

2.1 デジタル信号処理（DSP）

双一次変換を用いたIIRデジタルフィルタ設計において最も基本的かつ重要な技術。アナログプロトタイプフィルタからデジタルフィルタを設計する際、カットオフ周波数や中心周波数が正しい位置に来るよう事前補正を行う。

主な用途:

• IIRフィルタ設計（ローパス、ハイパス、バンドパス、ノッチフィルタ）
• オーディオイコライザ（パラメトリックEQ）のバイクアッドフィルタ係数計算
• デジタルオーディオエフェクト処理

2.2 制御工学

連続時間で設計した制御器をデジタル実装する際の離散化（Tustin法 = 双一次変換法）において使用される。

主な用途:

• PID制御器の離散化
• 比例共振（PR: Proportional-Resonant）制御器の離散化（マイクログリッド・パワーエレクトロニクス分野）
• ノッチフィルタ付き補償器の離散化
• モーション制御システムのデジタル化

2.3 コンピュータグラフィックス・画像処理

3Dレンダリングやテクスチャマッピングにおいて、複雑な画像変換を効率的に分解する手法として使用される。

主な用途:

• 3D画像ワーピングの因数分解（視差処理 + テクスチャマッピングへの分解）
• イメージベースドレンダリング（IBR）
• プロジェクションマッピングの幾何補正

2.4 VR/AR（バーチャルリアリティ）

ヘッドマウントディスプレイ（HMD）のレンズが引き起こす空間的・色収差を打ち消すための画像の事前歪み補正。

主な用途:

• HMDのレンズ歪み補正
• 広角ディスプレイの歪み補正
• 色収差補正

3. 主要な手法・アルゴリズム

3.1 双一次変換の周波数プリワーピング（信号処理・制御工学）

基本原理

双一次変換の基本式:

s = (2/T) * (z - 1) / (z + 1)

ここで T はサンプリング周期。この変換により、連続時間の角周波数 ω_a と離散時間の角周波数 ω_d の間に以下の非線形関係が生じる:

ωa = (2/T) * tan(ωd * T / 2)

低周波では ω_a ≈ ω_d だが、ナイキスト周波数に近づくほど歪みが大きくなる。

プリワーピングの補正式

特定の周波数 ω₀ で連続時間と離散時間の周波数応答を一致させたい場合、変換式を以下のように修正する:

s = [ω0 / tan(ω0 * T / 2)] * (z - 1) / (z + 1)

つまり、通常の係数 K = 2/T の代わりに K = ω₀ / tan(ω₀ * T / 2) を使用する。

MATLABでの実装例

MATLABでは c2d 関数のオプションで指定可能:

opt = c2dOptions('Method', 'tustin', 'PrewarpFrequency', 3.4);
sys_d = c2d(sys_c, Ts, opt);

bilinear 関数でも周波数プリワーピングを指定できる。

設計上の注意点

• プリワーピング周波数はクロスオーバー周波数やカットオフ周波数など、最も重要な周波数に設定するのが一般的
• プリワーピング周波数がナイキスト周波数（π/T）に近すぎると、低周波域に深刻な歪みが生じ、極端な場合は離散系が不安定になる
• DC（0 Hz）と指定周波数 ω₀ の2点で周波数応答が正確に一致する

3.2 オーディオEQバイクアッドフィルタにおけるプリワーピング

Robert Bristow-Johnsonの「Audio EQ Cookbook」は、オーディオ用バイクアッドフィルタ（2次IIRフィルタ）の係数計算における標準的参考資料である。

この手法では:

1. アナログプロトタイプフィルタを設計
2. 双一次変換でデジタル化
3. 周波数プリワーピングにより、ピーク周波数・カットオフ周波数・帯域幅を正しい位置に補正

ローパス、ハイパス、バンドパス、ノッチ、ピーキングEQ、ローシェルフ、ハイシェルフなど各種フィルタタイプに対応した公式が提供されている。

3.3 PR制御器のTustin法プリワーピング離散化

比例共振（PR）制御器は共振周波数付近で非常に狭い帯域を持つため、離散化手法の選択が特に重要である。

PR制御器の離散化では:

• 前進オイラー法・後退オイラー法: 共振周波数がずれる問題がある
• Tustin法 + プリワーピング: 共振周波数が正確に維持される

マイクログリッドや再生可能エネルギーのパワーコンバータ制御において、PWM変換器の電流制御にプリワーピング付きPR制御器が適用されている。

3.4 3D画像ワーピングの因数分解（コンピュータグラフィックス）

Mark, McMillan, Bishopらの手法（1997-1998年）では、3D画像ワーピング方程式を以下の2段階に分解する:

1. プリワープ（Pre-warp）: 視差（parallax）効果のみを処理。1次元の走査線・列方向の操作で実装可能で、隣接2ピクセル間の補間のみで済む
2. テクスチャマッピング: スケーリング、回転、残りの射影変換を標準的なグラフィックスハードウェアで処理

この分解により、複雑な3Dワーピングを効率的にソフトウェア・ハードウェア実装できる。

3.5 VR/HMDにおけるプリワーピング

HMDのレンズ歪み補正には主に以下の3つのアプローチがある:

1. イメージ空間法（ポストプロセス）: レンダリング結果をテクスチャに描画し、歪みメッシュに投影して補正。最も一般的
2. オブジェクト空間法（頂点プリワーピング）: 頂点シェーダでシーンのジオメトリを事前に歪ませる。モバイルVRで性能向上に有効だが、頂点間隔が広い場合に歪みが残る問題がある
3. 仮想光学系法: 仮想カメラに補正光学系をモデル化する

Pohl, Johnson, Bolkart（2013年）の研究により、バイキュービックフィルタリングが画像品質を向上させることが示されている。

4. 実際の応用例

4.1 デジタルオーディオ機器

• パラメトリックイコライザ: スタジオ用ミキサーやDAW（Digital Audio Workstation）のEQプラグインで使用。バイクアッドフィルタ係数の計算にプリワーピングが不可欠
• デジタルクロスオーバーネットワーク: スピーカーシステムの帯域分割フィルタ
• オーディオエフェクト: リバーブ、コーラス等のフィルタ部分

4.2 電力制御システム

• マイクログリッドのPWMインバータ制御: 系統連系型電力変換器のAC電流制御にPR制御器を使用し、Tustin法+プリワーピングで離散化。THD（全高調波歪率）の低減に貢献
• 再生可能エネルギーシステム: 太陽光発電・風力発電のグリッド接続インバータ
• 能動電力フィルタ: 高調波補償制御

4.3 モーション制御

• MATLABのControl System Toolboxでは、連続時間で設計した補償器をTustin法+プリワーピングで離散化する機能が標準提供されている
• サーボモータのデジタル制御器設計

4.4 VR/ARヘッドセット

• Oculus Rift, HTC Vive, PlayStation VR等: レンズ歪み補正のプリワーピングはすべての現代的VRヘッドセットで使用されている
• 色収差補正: R/G/Bチャンネルごとに異なるプリワーピングを適用

4.5 イメージベースドレンダリング

• リアルタイム視点合成: 既存画像から新しい視点の画像をリアルタイム生成する際に、プリワープ+テクスチャマッピングの分解手法が使用される
• リリーフテクスチャマッピング: 深度情報を持つテクスチャの効率的なレンダリング

4.6 プロジェクションマッピング

• 非平面スクリーンへの投影時の幾何歪み補正
• マルチプロジェクタシステムの画像位置合わせ

5. 参考文献

信号処理・制御工学分野

1. Tustin, A. “A method of analysing the behaviour of linear systems in terms of time series.” Journal of the Institution of Electrical Engineers – Part IIA, Vol. 94, No. 1, May 1947. (DOI)
2. Bristow-Johnson, R. “Cookbook formulae for audio EQ biquad filter coefficients.” (W3C Web Audio EQ Cookbook)
3. Bristow-Johnson, R. “The Equivalence of Various Methods of Computing Biquad Coefficients for Audio Parametric Equalizers.” (PDF)
4. Teodorescu, R., Blaabjerg, F., Liserre, M., Loh, P.C. “Proportional-resonant controllers and filters for grid-connected voltage-source converters.” IEE Proceedings – Electric Power Applications, 2006.
5. Mandava, S., Gudipalli, A., Amutha Prabha, N., Rajini, G. “Control of Micro-grid by Discretized PR Controller Using Tustin Frequency Pre-wrapping Method.” Springer, 2021. (DOI)

コンピュータグラフィックス分野

6. Mark, W., McMillan, L., Bishop, G. “Post-Rendering 3D Warping.” Proceedings of the Symposium on Interactive 3D Graphics (I3D), pp. 7-16, 1997.
7. Rafferty, M., Aliaga, D., Lastra, A. VRAIS ’98, 1998.
8. McMillan, L., Bishop, G. “Plenoptic modeling: An image-based rendering system.” Computer Graphics (SIGGRAPH ’95), pp. 39-46, 1995.
9. Heckbert, P.S. “Fundamentals of Texture Mapping and Image Warping.” Master’s thesis, UCB/CSD-89/516, UC Berkeley, 1989. (Link)

VR/AR分野

10. Pohl, D., Johnson, G.S., Bolkart, T. “Improved pre-warping for wide angle, head mounted displays.” Proceedings of the 19th ACM Symposium on Virtual Reality Software and Technology (VRST ’13), 2013. (DOI)

教科書・オンラインリソース

• Wikipedia – Bilinear transform
• MathWorks – Continuous-Discrete Conversion Methods
• MathWorks – bilinear