Jak obliczyć optymalny poziom powiększenia, aby wyświetlić dwa lub więcej punktów na mapie

Chcemy wyświetlić kilka znaczników na mapie statycznej i chcemy obliczyć optymalny poziom powiększenia, tak jak robi to Google Maps. Obliczyliśmy już prostokąt obwiedni i punkt środkowy mapy, ale teraz mamy trudności z obliczeniem prawidłowego poziomu powiększenia, aby wyświetlić cały prostokąt obwiedni. Czy ktoś może skierować nas w dobrym kierunku?

Aby uzyskać poziom powiększenia, musisz znać wymiary w pikselach mapy. Musisz także wykonać matematykę w sferycznych współrzędnych merkatora.

1. Przelicz szerokość, długość na sferyczny merkator x, y.
2. Uzyskaj odległość między dwoma punktami w kulistym mercatorze.
3. Równik ma rzut około 40 metrów, a płytki mają szerokość 256 pikseli, więc długość piksela tej mapy przy danym poziomie powiększenia wynosi około 256 * odległość / 40000000 * 2 ^ zbliżenie . Spróbuj powiększyć = 0, zoom = 1, zoom = 2, aż odległość będzie zbyt długa, aby pomieścić wymiary mapy w pikselach.

Oto kod C #, którego używam w Maperitive :

    public void ZoomToArea (Bounds2 mapArea, float paddingFactor)
    {
        double ry1 = Math.Log((Math.Sin(GeometryUtils.Deg2Rad(mapArea.MinY)) + 1) 
            / Math.Cos(GeometryUtils.Deg2Rad(mapArea.MinY)));
        double ry2 = Math.Log((Math.Sin(GeometryUtils.Deg2Rad(mapArea.MaxY)) + 1) 
            / Math.Cos(GeometryUtils.Deg2Rad(mapArea.MaxY)));
        double ryc = (ry1 + ry2) / 2;
        double centerY = GeometryUtils.Rad2Deg(Math.Atan(Math.Sinh(ryc)));

        double resolutionHorizontal = mapArea.DeltaX / Viewport.Width;

        double vy0 = Math.Log(Math.Tan(Math.PI*(0.25 + centerY/360)));
        double vy1 = Math.Log(Math.Tan(Math.PI*(0.25 + mapArea.MaxY/360)));
        double viewHeightHalf = Viewport.Height/2.0f;
        double zoomFactorPowered = viewHeightHalf
            / (40.7436654315252*(vy1 - vy0));
        double resolutionVertical = 360.0 / (zoomFactorPowered * 256);

        double resolution = Math.Max(resolutionHorizontal, resolutionVertical) 
            * paddingFactor;
        double zoom = Math.Log(360 / (resolution * 256), 2);
        double lon = mapArea.Center.X;
        double lat = centerY;

        CenterMapOnPoint(new PointD2(lon, lat), zoom);
    }

• mapArea: obwiednia w długich / długich współrzędnych (x = długi, y = lat)
• paddingFactor: można tego użyć, aby uzyskać efekt „120%”, o którym mówi ThomM. Wartość 1,2 da ci 120%.

Zauważ, że w moim przypadku zoommoże być liczbą rzeczywistą. W przypadku map internetowych potrzebna jest całkowita wartość powiększenia, więc powinieneś użyć czegoś takiego jak (int)Math.Floor(zoom)to.

Oczywiście ten kod dotyczy tylko projekcji Web Mercator.

Jeśli używasz OpenLayers Map.getZoomForExtent, obliczy najwyższy poziom powiększenia, który może zmieścić się w całym zakresie na mapie. Że extentmusi być w projekcji mapy. Możesz także użyć a, fill_factoraby uniknąć wyświetlania punktów na krawędzi mapy i max_zoomograniczyć możliwy zoom:

extent = extent.scale(1/fill_ratio);
var zoom = Math.min(map.getZoomForExtent(map_extent), max_zoom);
map.setCenter(extent.getCenterLonLat(), zoom);