OpencachingDeutschland / oc-server3

    }

    // array(zone, letter, north, east)

    public function getUTM()

    {

        /* 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 wandelt geographische Koordinaten in UTM Koordinaten

           um. Geographische Länge lw und Breite bw müssen im WGS84 Datum

           gegeben sein. Berechnet werden UTM Zone, Ostwert ew und Nordwert nw.*/

        //Geographische Länge lw und Breite bw im WGS84 Datum

        if ($this->nLon == 0 || $this->nLat == 0) {

            return [

                'zone' => '',

                'letter' => '',

                'north' => 'N ' . 0,

                'east' => 'E ' . 0,

            ];

        }

        if ($this->nLon <= - 180 || $this->nLon > 180 || $this->nLat <= - 80 || $this->nLat >= 84) {

            // Werte nicht im Bereich des UTM Systems -180° <= nLon < +180°, -80° < nLat < 84° N

            return [

                '',

                '',

                0,

                0,

            ];

        }

        $lw = (float) $this->nLon;

        $bw = (float) $this->nLat;

        //WGS84 Datum

        //Große Halbachse a und Abplattung f

        $a = 6378137.000;

        $f = 3.35281068e-3;

        $b_sel = 'CDEFGHJKLMNPQRSTUVWXX';

        //Polkrümmungshalbmesser c

        $c = $a / (1 - $f);

        //Quadrat der zweiten numerischen Exzentrizität

        $ex2 = (2 * $f - $f * $f) / ((1 - $f) * (1 - $f));

        $ex4 = $ex2 * $ex2;

        $ex6 = $ex4 * $ex2;

        $ex8 = $ex4 * $ex4;

        //Koeffizienten zur Berechnung der Meridianbogenlänge

        $e0 = $c * (pi() / 180) * (1 - 3 * $ex2 / 4 + 45 * $ex4 / 64 - 175 * $ex6 / 256 + 11025 * $ex8 / 16384);

        $e2 = $c * (-3 * $ex2 / 8 + 15 * $ex4 / 32 - 525 * $ex6 / 1024 + 2205 * $ex8 / 4096);

        $e4 = $c * (15 * $ex4 / 256 - 105 * $ex6 / 1024 + 2205 * $ex8 / 16384);

        $e6 = $c * (-35 * $ex6 / 3072 + 315 * $ex8 / 12288);

        //Längenzone lz und Breitenzone (Band) bz

        $lzn = intval(($lw + 180) / 6) + 1;

        if ($bw >= 56.0 && $bw < 64.0 && $lw >= 3.0 && $lw < 12.0) {

            $lzn = 32;

        }

        // Special zones for Svalbard.

        if ($bw >= 72.0 && $bw < 84.0) {

            if ($lw >= 0.0 && $lw < 9.0) {

                $lzn = 31;

            } elseif ($lw >= 9.0 && $lw < 21.0) {

                $lzn = 33;

            } elseif ($lw >= 21.0 && $lw < 33.0) {

                $lzn = 35;

            } elseif ($lw >= 33.0 && $lw < 42.0) {

                $lzn = 37;

            }

        }

        $lz = "$lzn";

        if ($lzn < 10) {

            $lz = '0' . $lzn;

        }

        $bd = (int) (1 + ($bw + 80) / 8);

        $bz = substr($b_sel, $bd - 1, 1);

        //Geographische Breite in Radianten br

        $br = $bw * pi() / 180;

        $tan1 = tan($br);

        $tan2 = $tan1 * $tan1;

        $tan4 = $tan2 * $tan2;

        $cos1 = cos($br);

        $cos2 = $cos1 * $cos1;

        $cos4 = $cos2 * $cos2;

        $cos3 = $cos2 * $cos1;

        $cos5 = $cos4 * $cos1;

        $etasq = $ex2 * $cos2;

        //Querkrümmungshalbmesser nd

        $nd = $c / sqrt(1 + $etasq);

        //Meridianbogenlänge g aus gegebener geographischer Breite bw

        $g = ($e0 * $bw) + ($e2 * sin(2 * $br)) + ($e4 * sin(4 * $br)) + ($e6 * sin(6 * $br));

        //Längendifferenz dl zum Bezugsmeridian lh

        $lh = ($lzn - 30) * 6 - 3;

        $dl = ($lw - $lh) * pi() / 180;

        $dl2 = $dl * $dl;

        $dl4 = $dl2 * $dl2;

        $dl3 = $dl2 * $dl;

        $dl5 = $dl4 * $dl;

        //Maßstabsfaktor auf dem Bezugsmeridian bei UTM Koordinaten m = 0.9996

        //Nordwert nw und Ostwert ew als Funktion von geographischer Breite und Länge

        if ($bw < 0) {

            $nw = 10e6 + 0.9996 * ($g + $nd * $cos2 * $tan1 * $dl2 / 2 +

                    $nd * $cos4 * $tan1 * (5 - $tan2 + 9 * $etasq) * $dl4 / 24);

        } else {

            $nw = 0.9996 * ($g + $nd * $cos2 * $tan1 * $dl2 / 2 +

                    $nd * $cos4 * $tan1 * (5 - $tan2 + 9 * $etasq) * $dl4 / 24);

        }

        $ew = 0.9996 * ($nd * $cos1 * $dl + $nd * $cos3 * (1 - $tan2 + $etasq) * $dl3 / 6 +

                $nd * $cos5 * (5 - 18 * $tan2 + $tan4) * $dl5 / 120) + 500000;

        $nk = $nw - (int) $nw;

        if ($nk < 0.5) {

            $nw = '' . (int) $nw;

        } else {

            $nw = '' . ((int) $nw + 1);

        }

        while (strlen($nw) < 7) {

            $nw = '0' . $nw;

        }

        $nk = $ew - (int) $ew;

        if ($nk < 0.5) {

            $ew = '0' . (int) $ew;

        } else {

            $ew = '0' . intval($ew + 1);

        }

        return [

            'zone' => $lz,

            'letter' => $bz,

            'north' => 'N ' . floor($nw),

            'east' => 'E ' . floor($ew),

        ];

    }

    // return string

    public function getGK()

     *

     * @return array|string[]

     */

    public function getUTM()

    : 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 wandelt geographische Koordinaten in UTM Koordinaten

           um. Geographische Länge lw und Breite bw müssen im WGS84 Datum

           gegeben sein. Berechnet werden UTM Zone, Ostwert ew und Nordwert nw.*/

        //Geographische Länge lw und Breite bw im WGS84 Datum

        if ($this->nLon == 0 || $this->nLat == 0) {

            return [

                'zone' => '',

                'letter' => '',

                'north' => 'N ' . 0,

                'east' => 'E ' . 0,

            ];

        }

        if ($this->nLon <= - 180 || $this->nLon > 180 || $this->nLat <= - 80 || $this->nLat >= 84) {

            // Werte nicht im Bereich des UTM Systems -180° <= nLon < +180°, -80° < nLat < 84° N

            return [

                '',

                '',

                0,

                0,

            ];

        }

        $lw = (float) $this->nLon;

        $bw = (float) $this->nLat;

        //WGS84 Datum

        //Große Halbachse a und Abplattung f

        $a = 6378137.000;

        $f = 3.35281068e-3;

        $b_sel = 'CDEFGHJKLMNPQRSTUVWXX';

        //Polkrümmungshalbmesser c

        $c = $a / (1 - $f);

        //Quadrat der zweiten numerischen Exzentrizität

        $ex2 = (2 * $f - $f * $f) / ((1 - $f) * (1 - $f));

        $ex4 = $ex2 * $ex2;

        $ex6 = $ex4 * $ex2;

        $ex8 = $ex4 * $ex4;

        //Koeffizienten zur Berechnung der Meridianbogenlänge

        $e0 = $c * (pi() / 180) * (1 - 3 * $ex2 / 4 + 45 * $ex4 / 64 - 175 * $ex6 / 256 + 11025 * $ex8 / 16384);

        $e2 = $c * (- 3 * $ex2 / 8 + 15 * $ex4 / 32 - 525 * $ex6 / 1024 + 2205 * $ex8 / 4096);

        $e4 = $c * (15 * $ex4 / 256 - 105 * $ex6 / 1024 + 2205 * $ex8 / 16384);

        $e6 = $c * (- 35 * $ex6 / 3072 + 315 * $ex8 / 12288);

        //Längenzone lz und Breitenzone (Band) bz

        $lzn = intval(($lw + 180) / 6) + 1;

        if ($bw >= 56.0 && $bw < 64.0 && $lw >= 3.0 && $lw < 12.0) {

            $lzn = 32;

        }

        // Special zones for Svalbard.

        if ($bw >= 72.0 && $bw < 84.0) {

            if ($lw >= 0.0 && $lw < 9.0) {

                $lzn = 31;

            } elseif ($lw >= 9.0 && $lw < 21.0) {

                $lzn = 33;

            } elseif ($lw >= 21.0 && $lw < 33.0) {

                $lzn = 35;

            } elseif ($lw >= 33.0 && $lw < 42.0) {

                $lzn = 37;

            }

        }

        $lz = "$lzn";

        if ($lzn < 10) {

            $lz = '0' . $lzn;

        }

        $bd = (int) (1 + ($bw + 80) / 8);

        $bz = substr($b_sel, $bd - 1, 1);

        //Geographische Breite in Radianten br

        $br = $bw * pi() / 180;

        $tan1 = tan($br);

        $tan2 = $tan1 * $tan1;

        $tan4 = $tan2 * $tan2;

        $cos1 = cos($br);

        $cos2 = $cos1 * $cos1;

        $cos4 = $cos2 * $cos2;

        $cos3 = $cos2 * $cos1;

        $cos5 = $cos4 * $cos1;

        $etasq = $ex2 * $cos2;

        //Querkrümmungshalbmesser nd

        $nd = $c / sqrt(1 + $etasq);

        //Meridianbogenlänge g aus gegebener geographischer Breite bw

        $g = ($e0 * $bw) + ($e2 * sin(2 * $br)) + ($e4 * sin(4 * $br)) + ($e6 * sin(6 * $br));

        //Längendifferenz dl zum Bezugsmeridian lh

        $lh = ($lzn - 30) * 6 - 3;

        $dl = ($lw - $lh) * pi() / 180;

        $dl2 = $dl * $dl;

        $dl4 = $dl2 * $dl2;

        $dl3 = $dl2 * $dl;

        $dl5 = $dl4 * $dl;

        //Maßstabsfaktor auf dem Bezugsmeridian bei UTM Koordinaten m = 0.9996

        //Nordwert nw und Ostwert ew als Funktion von geographischer Breite und Länge

        if ($bw < 0) {

            $nw = 10e6 + 0.9996 * ($g + $nd * $cos2 * $tan1 * $dl2 / 2 +

                                   $nd * $cos4 * $tan1 * (5 - $tan2 + 9 * $etasq) * $dl4 / 24);

        } else {

            $nw = 0.9996 * ($g + $nd * $cos2 * $tan1 * $dl2 / 2 +

                            $nd * $cos4 * $tan1 * (5 - $tan2 + 9 * $etasq) * $dl4 / 24);

        }

        $ew = 0.9996 * ($nd * $cos1 * $dl + $nd * $cos3 * (1 - $tan2 + $etasq) * $dl3 / 6 +

                        $nd * $cos5 * (5 - 18 * $tan2 + $tan4) * $dl5 / 120) + 500000;

        $nk = $nw - (int) $nw;

        if ($nk < 0.5) {

            $nw = '' . (int) $nw;

        } else {

            $nw = '' . ((int) $nw + 1);

        }

        while (strlen($nw) < 7) {

            $nw = '0' . $nw;

        }

        $nk = $ew - (int) $ew;

        if ($nk < 0.5) {

            $ew = '0' . (int) $ew;

        } else {

            $ew = '0' . intval($ew + 1);

        }

        return [

            'zone' => $lz,

            'letter' => $bz,

            'north' => 'N' . floor($nw),

            'east' => 'E' . floor($ew),

        ];

    }

    /**

     * Gauß Krüger