intermédiaire - no use

2018-10-06 09:56:23 +02:00
parent ebf5b0d648
commit e121c6f6dd
11 changed files with 1440 additions and 0 deletions
--- a/global/inc/cubic.php
+++ b/global/inc/cubic.php
@@ -0,0 +1,153 @@
+<?php
+
+class Plot 
+{
+  private $aCoords;
+  function __construct(&$aCoords) 
+  {
+    $this->aCoords = &$aCoords;
+  }
+
+  public function drawLine($vImage, $vColor, $iPosX = 0, $imaxX = false) 
+  {
+    $maxX = $imaxX;
+    if ($imaxX === false) $maxX = imagesx($vImage);
+
+    reset($this->aCoords);
+    list($iPrevX, $iPrevY) = each($this->aCoords);
+
+    while (list ($x, $y) = each($this->aCoords)) 
+    {
+      $laCouleur = null;
+      if (!is_array($vColor)) $laCouleur = $vColor;
+      else
+      {
+        $s = count($vColor) - 1;
+        $laCouleur = $vColor[$s]->color;
+        $pct = $x / $maxX;
+        while(($s>0)&&($pct < $vColor[$s]->pct))
+        {
+          $s -= 1;
+          $laCouleur = $vColor[$s]->color;
+        }
+      }
+      imageline($vImage, round($iPrevX), round($iPrevY), round($x), round($y), $laCouleur);
+      $iPrevX = $x;
+      $iPrevY = $y;
+    }
+  }
+
+  public function drawDots($vImage, $vColor, $iPosX = 0, $iPosY = false, $iDotSize = 1) {
+    if ($iPosY === false)
+      $iPosY = imagesy($vImage);
+    $vBorderColor = imagecolorallocate($vImage, 0, 0, 0);
+    foreach ($this->aCoords as $x => $y) {
+      imagefilledellipse($vImage, $iPosX + round($x), $iPosY - round($y), $iDotSize, $iDotSize, $vColor);
+      imageellipse($vImage, $iPosX + round($x), $iPosY - round($y), $iDotSize, $iDotSize, $vBorderColor);
+    }
+  }
+
+  public function drawAxis($vImage, $vColor, $iPosX = 0, $iPosY = false) {
+    if ($iPosY === false) $iPosY = imagesy($vImage);
+    $vImageWidth = imagesx($vImage);
+    imageline($vImage, $iPosX, $iPosY, $iPosX, 0, $vColor);
+    imageline($vImage, $iPosX, $iPosY, $vImageWidth, $iPosY, $vColor);
+    imagefilledpolygon($vImage, array($iPosX, 0, $iPosX - 3, 5, $iPosX + 3, 5), 3, $vColor);
+    imagefilledpolygon($vImage, array($vImageWidth, $iPosY, $vImageWidth - 5, $iPosY - 3, $vImageWidth - 5, $iPosY + 3), 3, $vColor);
+  }
+}
+
+class CubicSplines 
+{
+  protected $aCoords;
+  protected $aCrdX;
+  protected $aCrdY;
+  protected $aSplines = array();
+  protected $iMinX;
+  protected $iMaxX;
+  protected $iStep;
+  protected function prepareCoords(&$aCoords, $iStep, $iMinX = -1, $iMaxX = -1) {
+    $this->aCrdX = array();
+    $this->aCrdY = array();
+    $this->aCoords = array();
+    ksort($aCoords);
+    foreach ($aCoords as $x => $y) {
+      $this->aCrdX[] = $x;
+      $this->aCrdY[] = $y;
+    }
+    $this->iMinX = $iMinX;
+    $this->iMaxX = $iMaxX;
+    if ($this->iMinX == -1)
+      $this->iMinX = min($this->aCrdX);
+    if ($this->iMaxX == -1)
+      $this->iMaxX = max($this->aCrdX);
+    $this->iStep = $iStep;
+  }
+  public function setInitCoords(&$aCoords, $iStep = 1, $iMinX = -1, $iMaxX = -1) {
+    $this->aSplines = array();
+    if (count($aCoords) < 4) {
+      return false;
+    }
+    $this->prepareCoords($aCoords, $iStep, $iMinX, $iMaxX);
+    $this->buildSpline($this->aCrdX, $this->aCrdY, count($this->aCrdX));
+  }
+  public function processCoords() {
+    for ($x = $this->iMinX; $x <= $this->iMaxX; $x += $this->iStep) {
+      $this->aCoords[$x] = $this->funcInterp($x);
+    }
+    return $this->aCoords;
+  }
+  private function buildSpline($x, $y, $n) {
+    for ($i = 0; $i < $n; ++$i) {
+      $this->aSplines[$i]['x'] = $x[$i];
+      $this->aSplines[$i]['a'] = $y[$i];
+    }
+    $this->aSplines[0]['c'] = $this->aSplines[$n - 1]['c'] = 0;
+    $alpha[0] = $beta[0] = 0;
+    for ($i = 1; $i < $n - 1; ++$i) {
+      $h_i = $x[$i] - $x[$i - 1];
+      $h_i1 = $x[$i + 1] - $x[$i];
+      $A = $h_i;
+      $C = 2.0 * ($h_i + $h_i1);
+      $B = $h_i1;
+      $F = 6.0 * (($y[$i + 1] - $y[$i]) / $h_i1 - ($y[$i] - $y[$i - 1]) / $h_i);
+      $z = ($A * $alpha[$i - 1] + $C);
+      $alpha[$i] = - $B / $z;
+      $beta[$i] = ($F - $A * $beta[$i - 1]) / $z;
+    }
+    for ($i = $n - 2; $i > 0; --$i) {
+      $this->aSplines[$i]['c'] = $alpha[$i] * $this->aSplines[$i + 1]['c'] + $beta[$i];
+    }
+    for ($i = $n - 1; $i > 0; --$i) {
+      $h_i = $x[$i] - $x[$i - 1];
+      $this->aSplines[$i]['d'] = ($this->aSplines[$i]['c'] - $this->aSplines[$i - 1]['c']) / $h_i;
+      $this->aSplines[$i]['b'] = $h_i * (2.0 * $this->aSplines[$i]['c'] + $this->aSplines[$i - 1]['c']) / 6.0 + ($y[$i] - $y[$i - 1]) / $h_i;
+    }
+  }
+  private function funcInterp($x) {
+    $n = count($this->aSplines);
+    if ($x <= $this->aSplines[0]['x'])  {
+      $s = $this->aSplines[1];
+    } else {
+      if ($x >= $this->aSplines[$n - 1]['x']) {
+        $s = $this->aSplines[$n - 1];
+      } else {
+        $i = 0;
+        $j = $n - 1;
+        while ($i + 1 < $j) {
+          $k = $i + ($j - $i) / 2;
+          if ($x <= $this->aSplines[$k]['x']) {
+            $j = $k;
+          } else {
+            $i = $k;
+          }
+        }
+        $s = $this->aSplines[$j];
+      }
+    }
+    $dx = ($x - $s['x']);
+    return $s['a'] + ($s['b'] + ($s['c'] / 2.0 + $s['d'] * $dx / 6.0) * $dx) * $dx;
+  }
+}
+
+?>