phpOMS/Math/Geometry/ConvexHull/MonotoneChain.php

<?php
/**
 * Orange Management
 *
 * PHP Version 7.2
 *
 * @package    phpOMS\Math\Geometry\ConvexHull
 * @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\ConvexHull;

/**
 * Andrew's monotone chain convex hull algorithm class.
 *
 * @package    phpOMS\Math\Geometry\ConvexHull
 * @license    OMS License 1.0
 * @link       http://website.orange-management.de
 * @since      1.0.0
 */
final class MonotoneChain
{
    /**
     * Constructor.
     *
     * @since  1.0.0
     * @codeCoverageIgnore
     */
    private function __construct()
    {

    }

    /**
     * Create convex hull
     *
     * @param array<int, array<string, float|int>> $points Points (Point Cloud)
     *
     * @return array<int, array<string, float|int>>
     *
     * @since  1.0.0
     */
    public static function createConvexHull(array $points) : array
    {
        if (($n = \count($points)) > 1) {
            \uasort($points, [self::class, 'sort']);

            $k      = 0;
            $result = [];

            // Lower hull
            for ($i = 0; $i < $n; ++$i) {
                while ($k >= 2 && self::cross($result[$k - 2], $result[$k - 1], $points[$i]) <= 0) {
                    --$k;
                }

                $result[$k++] = $points[$i];
            }

            // Upper hull
            for ($i = $n - 2, $t = $k + 1; $i >= 0; --$i) {
                while ($k >= $t && self::cross($result[$k - 2], $result[$k - 1], $points[$i]) <= 0) {
                    --$k;
                }

                $result[$k++] = $points[$i];
            }

            \ksort($result);

            return \array_slice($result, 0, $k - 1);
        }

        return $points;
    }

    /**
     * Counter clock wise turn?
     *
     * @param array<string, float|int> $a Point a
     * @param array<string, float|int> $b Point b
     * @param array<string, float|int> $c Point c
     *
     * @return float
     *
     * @since  1.0.0
     */
    private static function cross(array $a, array $b, array $c) : float
    {
        return ($b['x'] - $a['x']) * ($c['y'] - $a['y']) - ($b['y'] - $a['y']) * ($c['x'] - $a['x']);
    }

    /**
     * Sort by x coordinate then by z coordinate
     *
     * @param array<string, float|int> $a Point a
     * @param array<string, float|int> $b Point b
     *
     * @return float
     *
     * @since  1.0.0
     */
    private static function sort(array $a, array $b) : float
    {
        return $a['x'] === $b['x'] ? $a['y'] - $b['y'] : $a['x'] - $b['x'];
    }
}