jsor / geokit Goto Github PK

$oBBox = BoundingBox::fromCoordinates([-118.45982654727163,33.858738223292775,-118.27313345274219,34.19377442763222]);
$fShrinkDistance = Distance::fromString('12km');
$oShrinkBBox = $oBBox->shrink($fShrinkDistance);

$oShrinkBox coords (note NE x is < SW x)
    [$oShrinkBox] => Geokit\BoundingBox Object
        (
            [southWest:Geokit\BoundingBox:private] => Geokit\Position Object
                (
                    [x:Geokit\Position:private] => -118.32986907162
                    [y:Geokit\Position:private] => 33.96665666694
                )

            [northEast:Geokit\BoundingBox:private] => Geokit\Position Object
                (
                    [x:Geokit\Position:private] => -118.40360502789
                    [y:Geokit\Position:private] => 34.085855983985
                )
        )

        /**
         * @see http://janmatuschek.de/LatitudeLongitudeBoundingCoordinates
         */
        private static function transformBoundingBox(BoundingBox $bbox, float $distanceInMeters): BoundingBox
        {
            $latSW = deg2rad($bbox->southWest()->latitude());
            $lngSW = deg2rad($bbox->southWest()->longitude());

            $latNE = deg2rad($bbox->northEast()->latitude());
            $lngNE = deg2rad($bbox->northEast()->longitude());

            $angularDistance = $distanceInMeters / Earth::RADIUS;

            $minLat = $latSW - $angularDistance;
            $maxLat = $latNE + $angularDistance;
            $fMinLat = min($minLat,$maxLat);
            $fMaxLat = max($minLat,$maxLat);
            $minLat = $fMinLat;
            $maxLat = $fMaxLat;

            $deltaLonSW = asin(sin($angularDistance) / cos($latSW));
            $deltaLonNE = asin(sin($angularDistance) / cos($latNE));

            $minLon = $lngSW - $deltaLonSW;
            $maxLon = $lngNE + $deltaLonNE;

            // asin can be neg Mark E.
            $fMinLon = min($minLon,$maxLon);
            $fMaxLon = max($minLon,$maxLon);
            $minLon = $fMinLon;
            $maxLon = $fMaxLon;

            $positionSW = Position::fromXY(rad2deg($minLon), rad2deg($minLat));
            $positionNE = Position::fromXY(rad2deg($maxLon), rad2deg($maxLat));

            // Check if we're shrinking too much
            if ($positionSW->latitude() > $positionNE->latitude()) {
                $center = $bbox->center();

                return BoundingBox::fromCornerPositions($center, $center);
            }

            return BoundingBox::fromCornerPositions($positionSW, $positionNE);
        }