OpencachingDeutschland / oc-server3

        ];

    }

    public function geo2gk($bp, $lp)

    {

        /* 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 GK Koordinaten

           um. Geographische Länge lp und Breite bp müssen im Potsdam Datum

           gegeben sein. Berechnet werden Rechtswert rw und Hochwert hw.*/

        //Geographische Länge lp und Breite bp im Potsdam Datum

        if ($lp == '' || $bp == '') {

            return [

                0,

                0,

            ];

        }

        $lp = (float) $lp;

        $bp = (float) $bp;

        // Potsdam Datum

        // Große Halbachse a und Abplattung f

        $a = 6377397.155; // + $falseeasting;

        $f = 3.34277321e-3;

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

        // Breite in Radianten

        $br = $bp * 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 bp

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

        // Längendifferenz dl zum Bezugsmeridian lh

        $kz = round($lp / 3);

        $lh = $kz * 3;

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

        $dl2 = $dl * $dl;

        $dl4 = $dl2 * $dl2;

        $dl3 = $dl2 * $dl;

        $dl5 = $dl4 * $dl;

        // Hochwert hw und Rechtswert rw als Funktion von geographischer Breite und Länge

        $hw = ($g + $nd * $cos2 * $tan1 * $dl2 / 2 + $nd * $cos4 * $tan1 * (5 - $tan2 + 9 * $etasq)

            * $dl4 / 24);

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

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

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

        if ($nk < 0.5) {

            $hw = (int) $hw;

        } else {

            $hw = ((int) $hw) + 1;

        }

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

        if ($nk < 0.5) {

            $rw = (int) $rw;

        } else {

            $rw = (int) ($rw + 1);

        }

        return [

            $rw,

            $hw,

        ];

    }

    // return string

    public function getRD()

     *

     * @return int[]

     */

    public function geo2gk($bp, $lp)

    {

        /* 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 GK Koordinaten

           um. Geographische Länge lp und Breite bp müssen im Potsdam Datum

           gegeben sein. Berechnet werden Rechtswert rw und Hochwert hw.*/

        //Geographische Länge lp und Breite bp im Potsdam Datum

        if ($lp == '' || $bp == '') {

            return [

                0,

                0,

            ];

        }

        $lp = (float) $lp;

        $bp = (float) $bp;

        // Potsdam Datum

        // Große Halbachse a und Abplattung f

        $a = 6377397.155; // + $falseeasting;

        $f = 3.34277321e-3;

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

        // Breite in Radianten

        $br = $bp * 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 bp

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

        // Längendifferenz dl zum Bezugsmeridian lh

        $kz = round($lp / 3);

        $lh = $kz * 3;

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

        $dl2 = $dl * $dl;

        $dl4 = $dl2 * $dl2;

        $dl3 = $dl2 * $dl;

        $dl5 = $dl4 * $dl;

        // Hochwert hw und Rechtswert rw als Funktion von geographischer Breite und Länge

        $hw = ($g + $nd * $cos2 * $tan1 * $dl2 / 2 + $nd * $cos4 * $tan1 * (5 - $tan2 + 9 * $etasq)

                                                     * $dl4 / 24);

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

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

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

        if ($nk < 0.5) {

            $hw = (int) $hw;

        } else {

            $hw = ((int) $hw) + 1;

        }

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

        if ($nk < 0.5) {

            $rw = (int) $rw;

        } else {

            $rw = (int) ($rw + 1);

        }

        return [

            $rw,

            $hw,

        ];

    }

    /**

     * RD Dutch Grid