PHPCoord\CoordinateOperation\GeographicValue

Complexity

Total Complexity

14

Size/Duplication

Total Lines	77
Duplicated Lines	0 %

Test Coverage

Coverage

91.89%

Importance

Changes	1
Bugs	0	Features	0

8 Methods

Rating	Name	Duplication	Size	Complexity
A	getLongitude()	0	3	1
A	getLatitude()	0	3	1
A	__construct()	0	7	2
A	normaliseLatitude()	0	10	3
A	normaliseLongitude()	0	10	3
A	asGeocentricValue()	0	14	2
A	getDatum()	0	3	1
A	getHeight()	0	3	1

<?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\Degree;
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 $latitude;

    private $longitude;

    private $height;

    private $datum;

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

        $this->height = $height ? $height->asMetres() : null;
    }

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

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

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

    public function getDatum(): Datum
    {
        return $this->datum;
    }

    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);
    }

    protected function normaliseLatitude(Angle $latitude): Angle
    {
        if ($latitude->asDegrees()->getValue() > 90) {
            return new Degree(90);
        }
        if ($latitude->asDegrees()->getValue() < -90) {
            return new Degree(-90);
        }

        return $latitude;
    }

    protected function normaliseLongitude(Angle $longitude): Angle
    {
        while ($longitude->asDegrees()->getValue() > 180) {
            $longitude = $longitude->subtract(new Degree(360));
        }
        while ($longitude->asDegrees()->getValue() <= -180) {
            $longitude = $longitude->add(new Degree(360));
        }

        return $longitude;
    }
}