PHPCoord\CoordinateOperation\GeographicValue

Complexity

Total Complexity

7

Size/Duplication

Total Lines	50
Duplicated Lines	0 %

Test Coverage

Coverage

100%

Importance

Changes	1
Bugs	0	Features	0

5 Methods

Rating	Name	Duplication	Size	Complexity
A	__construct()	0	8	2
A	getLongitude()	0	3	1
A	getLatitude()	0	3	1
A	getHeight()	0	3	1
A	asGeocentricValue()	0	14	2

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

namespace PHPCoord\CoordinateOperation;

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

/**
 * A geographic point w/out a CRS.
 * @internal
 */
class GeographicValue
{
    private Radian $latitude;

    private Radian $longitude;

    private ?Metre $height;

    private Datum $datum;

    public function __construct(Angle $latitude, Angle $longitude, ?Length $height, Datum $datum)
    {
        $this->latitude = $latitude->asRadians();
        $this->longitude = $longitude->asRadians();
        $this->datum = $datum;

        if ($height) {
            $this->height = $height->asMetres();
        }
    }

    public function getLatitude(): Radian
    {
        return $this->latitude;
    }

    public function getLongitude(): Radian
    {
        return $this->longitude;
    }

    public function getHeight(): ?Metre
    {
        return $this->height;
    }

    public function asGeocentricValue(): GeocentricValue
    {
        $a = $this->datum->getEllipsoid()->getSemiMajorAxis()->asMetres()->getValue();
        $e2 = $this->datum->getEllipsoid()->getEccentricitySquared();
        $latitude = $this->latitude->getValue();
        $longitude = $this->longitude->getValue() - $this->datum->getPrimeMeridian()->getGreenwichLongitude()->asRadians()->getValue();
        $h = isset($this->height) ? $this->height->getValue() : 0;

The method getValue() does not exist on null.

        $nu = $a / sqrt(1 - $e2 * (sin($latitude) ** 2));
        $x = ($nu + $h) * cos($latitude) * cos($longitude);
        $y = ($nu + $h) * cos($latitude) * sin($longitude);
        $z = ((1 - $e2) * $nu + $h) * sin($latitude);

        return new GeocentricValue(new Metre($x), new Metre($y), new Metre($z), $this->datum);
    }
}