mirror of
https://github.com/Karaka-Management/phpOMS.git
synced 2026-01-11 17:58:41 +00:00
409 lines
9.4 KiB
PHP
409 lines
9.4 KiB
PHP
<?php
|
|
/**
|
|
* Orange Management
|
|
*
|
|
* PHP Version 7.1
|
|
*
|
|
* @category TBD
|
|
* @package TBD
|
|
* @copyright Dennis Eichhorn
|
|
* @license OMS License 1.0
|
|
* @version 1.0.0
|
|
* @link http://website.orange-management.de
|
|
*/
|
|
declare(strict_types = 1);
|
|
|
|
namespace phpOMS\Math\Geometry\Shape\D2;
|
|
|
|
/**
|
|
* Polygon class.
|
|
*
|
|
* @category Framework
|
|
* @package phpOMS\Math
|
|
* @license OMS License 1.0
|
|
* @link http://website.orange-management.de
|
|
* @since 1.0.0
|
|
*/
|
|
class Polygon implements D2ShapeInterface
|
|
{
|
|
/**
|
|
* Epsilon for float comparison.
|
|
*
|
|
* @var float
|
|
* @since 1.0.0
|
|
*/
|
|
/* public */ const EPSILON = 0.00001;
|
|
|
|
/**
|
|
* Coordinates.
|
|
*
|
|
* These coordinates define the polygon
|
|
*
|
|
* @var array[]
|
|
* @since 1.0.0
|
|
*/
|
|
private $coord = [];
|
|
|
|
/**
|
|
* Polygon perimeter.
|
|
*
|
|
* @var float
|
|
* @since 1.0.0
|
|
*/
|
|
private $perimeter = 0.0;
|
|
|
|
/**
|
|
* Polygon surface.
|
|
*
|
|
* @var float
|
|
* @since 1.0.0
|
|
*/
|
|
private $surface = 0.0;
|
|
|
|
/**
|
|
* Interior angle sum of the polygon.
|
|
*
|
|
* @var int
|
|
* @since 1.0.0
|
|
*/
|
|
private $interiorAngleSum = 0;
|
|
|
|
/**
|
|
* Exterior angle sum of the polygon.
|
|
*
|
|
* @var float
|
|
* @since 1.0.0
|
|
*/
|
|
private $exteriorAngleSum = 0.0;
|
|
|
|
/**
|
|
* Polygon barycenter.
|
|
*
|
|
* @var float[]
|
|
* @since 1.0.0
|
|
*/
|
|
private $barycenter = ['x' => 0.0, 'y' => 0.0];
|
|
|
|
/**
|
|
* Polygon edge length.
|
|
*
|
|
* @var float
|
|
* @since 1.0.0
|
|
*/
|
|
private $edgeLength = 0.0;
|
|
|
|
/**
|
|
* Polygon inner length.
|
|
*
|
|
* @var float
|
|
* @since 1.0.0
|
|
*/
|
|
private $innerLength = 0.0;
|
|
|
|
/**
|
|
* Polygon inner edge angular.
|
|
*
|
|
* @var int
|
|
* @since 1.0.0
|
|
*/
|
|
private $innerEdgeAngular = 0;
|
|
|
|
/**
|
|
* Constructor.
|
|
*
|
|
* @since 1.0.0
|
|
*/
|
|
public function __construct()
|
|
{
|
|
}
|
|
|
|
/**
|
|
* Point polygon relative position
|
|
*
|
|
* @param array $point Point location
|
|
*
|
|
* @return int
|
|
*
|
|
* @since 1.0.0
|
|
*/
|
|
public function pointInPolygon(array $point) : int
|
|
{
|
|
return self::isPointInPolygon($point, $this->coord[] = $this->coord[0]);
|
|
}
|
|
|
|
/**
|
|
* Point polygon relative position
|
|
*
|
|
* @param array $point Point location
|
|
* @param array $polygon Polygon definition
|
|
*
|
|
* @return int -1 inside polygon 0 on vertice 1 outside
|
|
*
|
|
* @link http://erich.realtimerendering.com/ptinpoly/
|
|
* @since 1.0.0
|
|
*/
|
|
public static function isPointInPolygon(array $point, array $polygon) : int
|
|
{
|
|
$length = count($polygon);
|
|
|
|
// Polygon has to start and end with same point
|
|
if ($polygon[0]['x'] !== $polygon[$length - 1]['x'] || $polygon[0]['y'] !== $polygon[$length - 1]['y']) {
|
|
$polygon[] = $polygon[0];
|
|
}
|
|
|
|
// On vertex?
|
|
if (self::isOnVertex($point, $polygon)) {
|
|
return 0;
|
|
}
|
|
|
|
// Inside or ontop?
|
|
$countIntersect = 0;
|
|
$polygon_count = count($polygon);
|
|
|
|
// todo: return based on highest possibility not by first match
|
|
for ($i = 1; $i < $polygon_count; $i++) {
|
|
$vertex1 = $polygon[$i - 1];
|
|
$vertex2 = $polygon[$i];
|
|
|
|
if (abs($vertex1['y'] - $vertex2['y']) < self::EPSILON && abs($vertex1['y'] - $point['y']) < self::EPSILON && $point['x'] > min($vertex1['x'], $vertex2['x']) && $point['x'] < max($vertex1['x'], $vertex2['x'])) {
|
|
return 0; // boundary
|
|
}
|
|
|
|
if ($point['y'] > min($vertex1['y'], $vertex2['y']) && $point['y'] <= max($vertex1['y'], $vertex2['y']) && $point['x'] <= max($vertex1['x'], $vertex2['x']) && abs($vertex1['y'] - $vertex2['y']) >= self::EPSILON) {
|
|
$xinters = ($point['y'] - $vertex1['y']) * ($vertex2['x'] - $vertex1['x']) / ($vertex2['y'] - $vertex1['y']) + $vertex1['x'];
|
|
|
|
if (abs($xinters - $point['x']) < self::EPSILON) {
|
|
return 0; // boundary
|
|
}
|
|
|
|
if (abs($vertex1['x'] - $vertex2['x']) < self::EPSILON || $point['x'] < $xinters) {
|
|
$countIntersect++;
|
|
}
|
|
}
|
|
}
|
|
|
|
if ($countIntersect % 2 !== 0) {
|
|
return -1;
|
|
}
|
|
|
|
return 1;
|
|
}
|
|
|
|
/**
|
|
* Is point on vertex?
|
|
*
|
|
* @param array $point Point location
|
|
*
|
|
* @return bool
|
|
*
|
|
* @since 1.0.0
|
|
*/
|
|
public function onVertex(array $point) : bool
|
|
{
|
|
return self::isOnVertex($point, $this->coord);
|
|
}
|
|
|
|
/**
|
|
* Is point on vertex?
|
|
*
|
|
* @param array $point Point location
|
|
* @param array $polygon Polygon definition
|
|
*
|
|
* @return bool
|
|
*
|
|
* @since 1.0.0
|
|
*/
|
|
private static function isOnVertex(array $point, array $polygon) : bool
|
|
{
|
|
foreach ($polygon as $vertex) {
|
|
if (abs($point['x'] - $vertex['x']) < self::EPSILON && abs($point['y'] - $vertex['y']) < self::EPSILON) {
|
|
return true;
|
|
}
|
|
}
|
|
|
|
return false;
|
|
}
|
|
|
|
/**
|
|
* Set polygon coordinates.
|
|
*
|
|
* @param array[] $coord Coordinates
|
|
*
|
|
* @return void
|
|
*
|
|
* @since 1.0.0
|
|
*/
|
|
public function setCoordinates($coord) /* : void */
|
|
{
|
|
$this->coord = $coord;
|
|
}
|
|
|
|
/**
|
|
* Set polygon coordinate.
|
|
*
|
|
* @param int $i Index
|
|
* @param int|float $x X coordinate
|
|
* @param int|float $y Y coordinate
|
|
*
|
|
* @return void
|
|
*
|
|
* @since 1.0.0
|
|
*/
|
|
public function setCoordinate($i, $x, $y) /* : void */
|
|
{
|
|
$this->coord[$i] = ['x' => $x, 'y' => $y];
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getInteriorAngleSum() : int
|
|
{
|
|
$this->interiorAngleSum = (count($this->coord) - 2) * 180;
|
|
|
|
return $this->interiorAngleSum;
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getExteriorAngleSum()
|
|
{
|
|
return 360;
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getInteriorAngleSumFormula()
|
|
{
|
|
return '';
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getExteriorAngleSumFormula()
|
|
{
|
|
return '';
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getSurface() : float
|
|
{
|
|
return abs($this->getSignedSurface());
|
|
}
|
|
|
|
private function getSignedSurface() : float
|
|
{
|
|
$count = count($this->coord);
|
|
$surface = 0;
|
|
|
|
for ($i = 0; $i < $count - 1; $i++) {
|
|
$surface += $this->coord[$i]['x'] * $this->coord[$i + 1]['y'] - $this->coord[$i + 1]['x'] * $this->coord[$i]['y'];
|
|
}
|
|
|
|
$surface += $this->coord[$count - 1]['x'] * $this->coord[0]['y'] - $this->coord[0]['x'] * $this->coord[$count - 1]['y'];
|
|
|
|
$surface /= 2;
|
|
|
|
return $surface;
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function setSurface($surface) /* : void */
|
|
{
|
|
$this->reset();
|
|
|
|
$this->surface = $surface;
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function reset() /* : void */
|
|
{
|
|
$this->coord = [];
|
|
$this->barycenter = ['x' => 0.0, 'y' => 0.0];
|
|
$this->perimeter = 0.0;
|
|
$this->surface = 0.0;
|
|
$this->interiorAngleSum = 0;
|
|
$this->edgeLength = 0.0;
|
|
$this->innerLength = 0.0;
|
|
$this->innerEdgeAngular = 0;
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getSurfaceFormula()
|
|
{
|
|
return '';
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getPerimeter() : float
|
|
{
|
|
$count = count($this->coord);
|
|
$this->perimeter = sqrt(($this->coord[0]['x'] - $this->coord[$count - 1]['x']) ** 2 + ($this->coord[0]['y'] - $this->coord[$count - 1]['y']) ** 2);
|
|
|
|
for ($i = 0; $i < $count - 1; $i++) {
|
|
$this->perimeter += sqrt(($this->coord[$i + 1]['x'] - $this->coord[$i]['x']) ** 2 + ($this->coord[$i + 1]['y'] - $this->coord[$i]['y']) ** 2);
|
|
}
|
|
|
|
return $this->perimeter;
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function setPerimeter($perimeter) /* : void */
|
|
{
|
|
$this->reset();
|
|
|
|
$this->perimeter = $perimeter;
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getPerimeterFormula()
|
|
{
|
|
return '';
|
|
}
|
|
|
|
/**
|
|
* {@inheritdoc}
|
|
*/
|
|
public function getBarycenter()
|
|
{
|
|
$this->barycenter['x'] = 0;
|
|
$this->barycenter['y'] = 0;
|
|
|
|
$count = count($this->coord);
|
|
|
|
for($i = 0; $i < $count - 1; $i++) {
|
|
$mult = ($this->coord[$i]['x'] * $this->coord[$i + 1]['y'] - $this->coord[$i + 1]['x'] * $this->coord[$i]['y']);
|
|
$this->barycenter['x'] += ($this->coord[$i]['x'] + $this->coord[$i + 1]['x']) * $mult;
|
|
$this->barycenter['y'] += ($this->coord[$i]['y'] + $this->coord[$i + 1]['y']) * $mult;
|
|
}
|
|
|
|
$mult = ($this->coord[$count - 1]['x'] * $this->coord[0]['y'] - $this->coord[0]['x'] * $this->coord[$count - 1]['y']);
|
|
$this->barycenter['x'] += ($this->coord[$count - 1]['x'] + $this->coord[0]['x']) * $mult;
|
|
$this->barycenter['y'] += ($this->coord[$count - 1]['y'] + $this->coord[0]['y']) * $mult;
|
|
|
|
$surface = $this->getSignedSurface();
|
|
|
|
$this->barycenter['x'] = 1 / (6 * $surface) * $this->barycenter['x'];
|
|
$this->barycenter['y'] = 1 / (6 * $surface) * $this->barycenter['y'];
|
|
|
|
return $this->barycenter;
|
|
}
|
|
}
|