Bavix\Geo\Metrical

Complexity

Total Complexity

8

Size/Duplication

Total Lines	116
Duplicated Lines	0 %

Importance

Changes

0

5 Methods

Rating	Name	Duplication	Size	Complexity
A	axisComputed()	0	37	4
A	rectangle()	0	23	1
A	computedEps()	0	5	1
A	squareByHypotenuse()	0	5	1
A	square()	0	3	1

<?php

namespace Bavix\Geo;

use Bavix\Geo\Geometry\Line;
use Bavix\Geo\Geometry\Quadrangle;
use Bavix\Geo\Unit\Distance;

class Metrical
{

    /**
     * @param Coordinate $center
     * @param Distance   $unit
     *
     * @return Quadrangle
     */
    public function squareByHypotenuse(Coordinate $center, Distance $unit): Quadrangle
    {
        $item = clone $unit;
        $item->miles /= \sqrt(2);
        return $this->square($center, $item);
    }

    /**
     * @param Coordinate $center
     * @param Distance   $unit
     *
     * @return Quadrangle
     */
    public function square(Coordinate $center, Distance $unit): Quadrangle
    {
        return $this->rectangle($center, $unit, $unit);
    }

    /**
     * @param Coordinate $center
     * @param Distance   $unitX
     * @param Distance   $unitY
     *
     * @return Quadrangle
     */
    public function rectangle(Coordinate $center, Distance $unitX, Distance $unitY): Quadrangle
    {
        $vx = \rad2deg($unitX->nauticalMiles / 60.);
        $vy = \rad2deg($unitY->nauticalMiles / 60.);
        $dx = \deg2rad(\hypot($vx, $vx) / 2.);
        $dy = \deg2rad(\hypot(0, $vy));

        $computedX = $center->plus($dx, 0);
        $computedY = $center->plus(0, $dy);

        $lineX = new Line();
        $lineX->push($center);
        $lineY = clone $lineX;
        $lineX->push($computedX);
        $lineY->push($computedY);

        $latitude = $this->axisComputed($lineX, $unitX);
        $longitude = $this->axisComputed($lineY, $unitY, false);

        return Quadrangle::makeWithXY(
            $center,
            $unitX->miles * $latitude,
            $unitY->miles * $longitude
        );
    }

    /**
     * @param Line     $line
     * @param Distance $unit
     * @param bool     $isAxisX
     *
     * @return float
     */
    protected function axisComputed(Line $line, Distance $unit, bool $isAxisX = true): float
    {
        $computed = $line->pop();
        $center = $line->pop();

        $iterator = 0;
        $eps = $this->computedEps($unit);
        $computed = $computed->plus(
            $eps * $isAxisX,
            $eps * !$isAxisX
        );

        while ($iterator < 256) {

            $distance = $center->distanceTo($computed);
            $eps += $this->computedEps($distance);

            if ($distance->compareTo($unit)->greaterThanOrEqual()) {
                break;
            }

            $computed = $computed->plus(
                $eps * $isAxisX,
                $eps * !$isAxisX
            );

            $iterator++;

        }

        if ($isAxisX) {
            $result = $center->latitude->degrees - $computed->latitude->degrees;
        } else {
            $result = $center->longitude->degrees - $computed->longitude->degrees;
        }

        return \abs($result) / $unit->miles;
    }

    /**
     * @param Distance $axis
     *
     * @return float
     */
    protected function computedEps(Distance $axis): float
    {
        return \max(
            \round($axis->miles, 5) / 1e5,
            1e-7
        );
    }

}