[論文レビュー] Fast and accurate $G^1$ fitting of clothoid curves

Enrico Bertolazzi, Marco Frego|arXiv (Cornell University)|May 28, 2013

Advanced Numerical Analysis Techniques参考文献 17被引用数 4

ひとこと要約

本稿では、1本のクロソイド曲線を用いたG1エルミート補間のための高速で安定的かつ一様に高精度なアルゴリズムを提示する。3つの非線形方程式を解析的再定式化により1変数の根を求めることに還元し、特に直線や円弧に近い極限状態（例：直線や円弧）におけるフェルマー積分の漸近展開を用いることで、あらゆる設定で高精度を確保する。適切に設計された初期推定値を用いることで、ニュートン・ラプソン法による収束が4回未満の反復で達成され、特に極限状態付近において従来手法を上回る精度と安定性を発揮する。

ABSTRACT

A new effective solution to the problem of Hermite $G^1$ interpolation with a clothoid curve is here proposed, that is a clothoid that interpolates two given points in a plane with assigned unit tangent vectors. The interpolation problem is a system of three nonlinear equations with multiple solutions which is difficult to solve also numerically. Here the solution of this system is reduced to the computation of the zeros of one single function in one variable. The location of the zero associated to the relevant solution is studied analytically: the interval containing the zero where the solution is proved to exists and to be unique is provided. A simple guess function allows to find that zero with very few iterations in all possible configurations. The computation of the clothoid curves and the solution algorithm call for the evaluation of Fresnel related integrals. Such integrals need asymptotic expansions near critical values to avoid loss of precision. This is necessary when, for example, the solution of the interpolation problem is close to a straight line or an arc of circle. A simple algorithm is presented for efficient computation of the asymptotic expansion. The reduction of the problem to a single nonlinear function in one variable and the use of asymptotic expansions make the present solution algorithm fast and robust. In particular a comparison with algorithms present in literature shows that the present algorithm requires less iterations. Moreover accuracy is maintained in all possible configurations while other algorithms have a loss of accuracy near the transition zones.

研究の動機と目的

• ケース別分解による数値的不安定性を回避する、1本のクロソイドセグメントを用いたG1エルミート補間の統一的で安定的なアルゴリズムの開発。
• 直線や円弧のようないびつな状態（退化状態）に近い場合でも、精度の低下を防ぐ一様な精度を確保すること。
• G1補間問題を1変数の非線形方程式に還元し、ニュートン・ラプソン法による効率的かつ信頼性の高い解法を可能にすること。
• 曲率や長さが0に近づくような臨界パrameter領域で、数値的安定性と精度を維持するためのフェルマー関連積分の漸近展開を提供すること。

提案手法

• クロソイドの幾何的・解析的性質を活用して、3変数のG1エルミート補間問題を1変数の非線形方程式に再定式化する。
• 解が一意に存在し保証される理論的区間を導出し、安全かつ信頼性の高い根の探査を可能にする。
• あらゆる入力設定で高速収束（通常4回未満の反復）を実現する、堅牢な初期推定関数を提案する。
• 曲率や長さが0に近づく場合（例：直線や円弧に近い状態）における精度損失を回避するため、一般化フェルマー積分の漸近展開を用いる。
• 小・大のパrameterに対して特別なルーチンを備えた、フェルマー関連積分を計算するための包括的なアルゴリズムパイプライン（疑似コードを用いて）を開発する。
• 導出された初期推定値を用いてニュートン・ラプソン法を適用し、スカラー関数のゼロを効率的かつ高精度に計算する。

実験結果

リサーチクエスチョン

• RQ1 1本のクロソイドを用いたG1エルミート補間問題を、ケース分割を避けながら1変数の非線形方程式に還元できるか？
• RQ2 直線や円弧のような退化状態に近い場合の数値的不安定性を、クロソイドフィッティングでどのように排除できるか？
• RQ3 すべての可能なエルミートデータの設定において、高速収束（少ない反復回数）を実現する初期推定関数は何か？
• RQ4 パrameterが臨界値（例：小さな曲率や長さ）に近い場合に、フェルマー積分をどのように高精度に評価できるか？
• RQ5 統一されたアルゴリズムが、極限状態でも精度と収束速度を一様に維持でき、信頼性を損なわないか？

主な発見

• 提案手法により、G1エルミート補間問題が1変数の非線形方程式に還元され、ケース分解を避けた統一的で安定した解法が実現された。
• 解は明示的に証明された明確な区間内に存在し、一意に保証される。この区間は根の探査プロセスをガイドするのに用いられる。
• 初期推定関数により、テストされたあらゆる設定（特に退化状態付近を含む）で、ニュートン・ラプソン法の収束が4回未満の反復で達成された。
• フェルマー積分の漸近展開により、クロソイドが直線や円弧に近づく場合でも一様な精度が維持され、従来手法で観察された精度損失が解消された。
• 数値的比較により、本手法は特に過渡領域において、既存のアルゴリズムを精度と反復回数の両面で上回ることが示された。
• 本アルゴリズムは、パラメータ空間全体にわたり堅牢かつ高精度であり、ゼロ曲率やゼロ長さのような極端なケースでさえ、しきい値や特別な処理を必要とせず自然に扱える。