PHPCoord\CoordinateOperation\GeocentricValue

Complexity

Total Complexity

6

Size/Duplication

Total Lines	58
Duplicated Lines	0 %

Test Coverage

Coverage

100%

Importance

Changes	1
Bugs	0	Features	0

5 Methods

Rating	Name	Duplication	Size	Complexity
A	__construct()	0	6	1
A	getZ()	0	3	1
A	getX()	0	3	1
A	getY()	0	3	1
A	asGeographicValue()	0	23	2

<?php
/**
 * PHPCoord.
 *
 * @author Doug Wright
 */
declare(strict_types=1);

namespace PHPCoord\CoordinateOperation;

use function abs;
use function atan2;
use function cos;
use PHPCoord\Datum\Datum;
use PHPCoord\UnitOfMeasure\Angle\Radian;
use PHPCoord\UnitOfMeasure\Length\Length;
use PHPCoord\UnitOfMeasure\Length\Metre;
use function sin;
use function sqrt;

/**
 * A geocentric point w/out a CRS.
 * @internal
 */
class GeocentricValue
{
    protected const NEWTON_RAPHSON_CONVERGENCE = 1e-16;

    private Metre$x;

    private Metre $y;

    private Metre $z;

    private Datum $datum;

    public function __construct(Length $x, Length $y, Length $z, Datum $datum)
    {
        $this->x = $x->asMetres();
        $this->y = $y->asMetres();
        $this->z = $z->asMetres();
        $this->datum = $datum;
    }

    public function getX(): Metre
    {
        return $this->x;
    }

    public function getY(): Metre
    {
        return $this->y;
    }

    public function getZ(): Metre
    {
        return $this->z;
    }

    public function asGeographicValue(): GeographicValue
    {
        $a = $this->datum->getEllipsoid()->getSemiMajorAxis()->asMetres()->getValue();
        $e2 = $this->datum->getEllipsoid()->getEccentricitySquared();
        $x = $this->x->getValue();
        $y = $this->y->getValue();
        $z = $this->z->getValue();

        $longitude = atan2($y, $x);
        $p = sqrt(($x ** 2) + ($y ** 2));

        $latitude = atan2($z, ($p * (1 - $e2)));

        do {
            $phi1 = $latitude;
            $v = $a / sqrt(1 - $e2 * (sin($latitude) ** 2));
            $latitude = atan2($z + ($e2 * $v * sin($latitude)), $p);
        } while (abs($latitude - $phi1) >= self::NEWTON_RAPHSON_CONVERGENCE);

        $latitude += $this->datum->getPrimeMeridian()->getGreenwichLongitude()->asRadians()->getValue();
        $h = $p / cos($latitude) - $v;

        return new GeographicValue(new Radian($latitude), new Radian($longitude), new Metre($h), $this->datum);
    }
}