Oc\Repository\CoordinatesRepository::wgs2pot()

Complexity

Conditions	6
Paths	4

Size

Total Lines

83

Duplication

Lines	83
Ratio	100 %

Importance

Changes

0

<?php
    public function wgs2pot($bw, $lw)

This method seems to be duplicated in your project.

    : array {
        /* Copyright (c) 2006, HELMUT H. HEIMEIER
           Permission is hereby granted, free of charge, to any person obtaining a
           copy of this software and associated documentation files (the "Software"),
           to deal in the Software without restriction, including without limitation
           the rights to use, copy, modify, merge, publish, distribute, sublicense,
           and/or sell copies of the Software, and to permit persons to whom the
           Software is furnished to do so, subject to the following conditions:
           The above copyright notice and this permission notice shall be included
           in all copies or substantial portions of the Software.*/

        /* Die Funktion verschiebt das Kartenbezugssystem (map datum) vom
           WGS84 Datum (World Geodetic System 84) zum in Deutschland
           gebräuchlichen Potsdam-Datum. Geographische Länge lw und Breite
           bw gemessen in grad auf dem WGS84 Ellipsoid müssen
           gegeben sein. Ausgegeben werden geographische Länge lp
           und Breite bp (in grad) auf dem Bessel-Ellipsoid.
           Bei der Transformation werden die Ellipsoidachsen parallel
           verschoben um dx = -606 m, dy = -23 m und dz = -413 m.*/

        // Geographische Länge lw und Breite bw im WGS84 Datum
        if ($lw == '' || $bw == '') {
            return [
                0,
                0,
            ];
        }
        $lw = (float) $lw;
        $bw = (float) $bw;

        // Quellsystem WGS 84 Datum
        // Große Halbachse a und Abplattung fq
        $a = 6378137.000;
        $fq = 3.35281066e-3;

        // Zielsystem Potsdam Datum
        // Abplattung f
        $f = $fq - 1.003748e-5;

        // Parameter für datum shift
        $dx = - 606;
        $dy = - 23;
        $dz = - 413;

        // Quadrat der ersten numerischen Exzentrizität in Quell- und Zielsystem
        $e2q = (2 * $fq - $fq * $fq);
        $e2 = (2 * $f - $f * $f);

        // Breite und Länge in Radianten
        $b1 = $bw * (pi() / 180);
        $l1 = $lw * (pi() / 180);

        // Querkrümmungshalbmesser nd
        $nd = $a / sqrt(1 - $e2q * sin($b1) * sin($b1));

        // Kartesische Koordinaten des Quellsystems WGS84
        $xw = $nd * cos($b1) * cos($l1);
        $yw = $nd * cos($b1) * sin($l1);
        $zw = (1 - $e2q) * $nd * sin($b1);

        // Kartesische Koordinaten des Zielsystems (datum shift) Potsdam
        $x = $xw + $dx;
        $y = $yw + $dy;
        $z = $zw + $dz;

        // Berechnung von Breite und Länge im Zielsystem
        $rb = sqrt($x * $x + $y * $y);
        $b2 = (180 / pi()) * atan(($z / $rb) / (1 - $e2));

        if ($x > 0) {
            $l2 = (180 / pi()) * atan($y / $x);
        } elseif ($x < 0 && $y > 0) {
            $l2 = (180 / pi()) * atan($y / $x) + 180;
        } else {
            $l2 = (180 / pi()) * atan($y / $x) - 180;
        }

        return [
            $l2,
            $b2,
        ];
    }