【PHP実践｜実務向け】PHPで実現する滑らかなグラフ描画：3次スプライン曲線の実装

導入：なぜ「3次スプライン曲線」が必要なのか

Webアプリケーションでデータを可視化する際、単なる折れ線グラフでは角張った印象になりがちです。特に、時間の経過に伴う緩やかな変化や、統計データに基づいた予測値を表現する場合、点と点を滑らかに結ぶ「曲線」が求められます。PHPのGDライブラリを使用してグラフを描画する際、単純な線形補間（直線で結ぶ）ではなく「3次スプライン曲線」を用いることで、全ての点を通りながらも非常に滑らかで美しい曲線を描くことが可能になります。

基礎知識：3次スプライン曲線とは

3次スプライン曲線とは、与えられた複数の点の間を「3次多項式」で接続する手法です。各区間ごとに計算された多項式を繋ぎ合わせる際、隣接する区間との境界で「傾き」や「曲がり具合（2階微分）」が連続するように設計されているため、結果として非常に自然な曲線が得られます。数学的な計算量は増えますが、ビジネスレポートや分析ツールなど、視認性が重視される現場では非常に有用な技術です。

実装：スプライン曲線クラスの設計

以下に、与えられた座標配列を元にスプライン補間を行い、画像として出力する実用的なクラスの例を示します。

サンプルプログラム

以下のコードは、指定された点群を元に滑らかな曲線を描画するクラスです。

/
class CubicSpline {
private $points;
private $n;
private $q = [], $r = [], $s = [];

public function __construct(array $points) {
$this->points = $points;
$this->n = count($points) – 1;
$h = []; $b = []; $d = [];
$g = []; $u = [];

// ステップ1: 各区間の差分を計算
for ($i = 0; $i < $this->n; $i++) {
$h[$i] = $points[$i + 1][0] – $points[$i][0];
}
for ($i = 1; $i < $this->n; $i++) {
$b[$i] = 2.0 ($h[$i] + $h[$i – 1]);
$d[$i] = 3.0 (($points[$i + 1][1] – $points[$i][1]) / $h[$i] – ($points[$i][1] – $points[$i – 1][1]) / $h[$i – 1]);
}

// ステップ2: 係数の計算
$g[1] = $h[1] / $b[1];
for ($i = 2; $i < $this->n – 1; $i++) {
$g[$i] = $h[$i] / ($b[$i] – $h[$i – 1] $g[$i – 1]);
}
$u[1] = $d[1] / $b[1];
for ($i = 2; $i < $this->n; $i++) {
$u[$i] = ($d[$i] – $h[$i – 1] $u[$i – 1]) / ($b[$i] – $h[$i – 1] $g[$i – 1]);
}

// ステップ3: 3次式の係数rを算出
$this->r[0] = $this->r[$this->n] = 0.0;
for ($i = $this->n – 1; $i >= 1; $i–) {
$this->r[$i] = $u[$i] – $g[$i] $this->r[$i + 1];
}

// ステップ4: 補間用の係数q, sを算出
for ($i = 0; $i < $this->n; $i++) {
$this->q[$i] = ($points[$i + 1][1] – $points[$i][1]) / $h[$i] – $h[$i] ($this->r[$i + 1] + 2.0 $this->r[$i]) / 3.0;
$this->s[$i] = ($this->r[$i + 1] – $this->r[$i]) / (3.0 $h[$i]);
}
}

// 指定したX座標に対するY座標を算出
public function value($x) {
$i = 0;
for ($i1 = 1; $i1 < $this->n; $i1++) {
if ($x < $this->points[$i1][0]) { $i = $i1 – 1; break; }
}
$x2 = $x – $this->points[$i][0];
return $this->points[$i][1] + $x2 ($this->q[$i] + $x2 ($this->r[$i] + $this->s[$i] $x2));
}
}

// 利用例
$points = [[0, 20], [50, 80], [100, 30], [150, 90]];
$spline = new CubicSpline($points);
// この後、ループで $spline->value($x) を呼び出し、imageline() で描画します。
?>

応用・注意点：現場での運用

1. パフォーマンスの考慮: 大量のデータを描画する場合、ループ内で毎回計算を行うと処理負荷が高まります。必要に応じて計算結果をキャッシュする仕組みを導入してください。
2. 例外処理: $points のデータがX軸方向に昇順でない場合、正しく計算できません。入力データがソートされていることを保証するか、クラス内でソート処理を追加してください。
3. GDライブラリの制限: PHPのGDライブラリはアンチエイリアス処理がデフォルトで効かない場合があります。より高品質な描画を行いたい場合は、`imageantialias($image, true);` を必ずコールするようにしてください。