PHPCoord\CoordinateOperation\GTXGrid

Complexity

Total Complexity

9

Size/Duplication

Total Lines	84
Duplicated Lines	0 %

Test Coverage

Coverage

92%

Importance

Changes	2
Bugs	0	Features	0

4 Methods

Rating	Name	Duplication	Size	Complexity
A	__construct()	0	15	2
A	getRecord()	0	17	1
A	getHeader()	0	20	3
A	getValues()	0	14	3

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

namespace PHPCoord\CoordinateOperation;

use PHPCoord\UnitOfMeasure\Length\Metre;
use SplFileObject;

use function in_array;
use function substr;
use function unpack;

/**
 * @see https://vdatum.noaa.gov/docs/gtx_info.html for documentation
 */
class GTXGrid extends GeographicGeoidHeightGrid
{
    use BilinearInterpolation;

    private int $headerLength;
    private string $shiftDataType;

    public function __construct($filename)
    {
        $this->gridFile = new SplFileObject($filename);
        $this->storageOrder = self::STORAGE_ORDER_INCREASING_LONGITUDE_INCREASING_LATIITUDE;

        $header = $this->getHeader();
        $this->startX = $header['xlonsw'];
        $this->startY = $header['xlatsw'];
        $this->numberOfColumns = $header['nlon'];
        $this->numberOfRows = $header['nlat'];
        $this->columnGridInterval = $header['dlon'];
        $this->rowGridInterval = $header['dlat'];

        if ($this->startX > 180) { // normalise if necessary
            $this->startX -= 360;
        }
    }

    /**
     * @return Metre[]
     */
    public function getValues($x, $y): array
    {
        if ($x < $this->startX) { // normalise if necessary
            $x += 360;
        }

        $shift = $this->interpolate($x, $y)[0];

        // These are in millimeters for some reason... :/
        if (in_array($this->gridFile->getBasename(), ['vertconc.gtx', 'vertcone.gtx', 'vertconw.gtx'], true)) {
            $shift /= 1000;
        }

        return [new Metre($shift)];
    }

    private function getRecord(int $longitudeIndex, int $latitudeIndex): GridValues

The method getRecord() is not used, and could be removed.

    {
        $recordId = match ($this->storageOrder) {
            self::STORAGE_ORDER_INCREASING_LATITUDE_INCREASING_LONGITUDE => $longitudeIndex * $this->numberOfRows + $latitudeIndex,
            self::STORAGE_ORDER_INCREASING_LONGITUDE_DECREASING_LATIITUDE => ($this->numberOfRows - $latitudeIndex - 1) * $this->numberOfColumns + $longitudeIndex,
            self::STORAGE_ORDER_INCREASING_LONGITUDE_INCREASING_LATIITUDE => $latitudeIndex * $this->numberOfColumns + $longitudeIndex,
        };

        $offset = $this->headerLength + $recordId * 4;
        $this->gridFile->fseek($offset);
        $rawRow = $this->gridFile->fread(4);
        $data = unpack("{$this->shiftDataType}shift", $rawRow);

        return new GridValues(
            $longitudeIndex * $this->columnGridInterval + $this->startX,
            $latitudeIndex * $this->rowGridInterval + $this->startY,
            [$data['shift']]
        );
    }

    private function getHeader(): array
    {
        $this->gridFile->fseek(0);
        $rawHeader = $this->gridFile->fread(44);
        $ikind = substr($rawHeader, 40, 4);
        if (unpack('Nikind', $ikind)['ikind'] === 1) { // big endian
            $this->headerLength = 44;
            $this->shiftDataType = 'G';
            $data = unpack('Exlatsw/Exlonsw/Edlat/Edlon/Nnlat/Nnlon', $rawHeader);
        } elseif (unpack('Vikind', $ikind)['ikind'] === 1) { // little endian
            $this->headerLength = 44;
            $this->shiftDataType = 'g';
            $data = unpack('exlatsw/exlonsw/edlat/edlon/Vnlat/Vnlon', $rawHeader);
        } else { // not all files (e.g. NZ) have this endian check column, assume big endian
            $this->headerLength = 40;
            $this->shiftDataType = 'G';
            $data = unpack('Exlatsw/Exlonsw/Edlat/Edlon/Nnlat/Nnlon', $rawHeader);
        }

        return $data;
    }
}