Как рассчитать оптимальный уровень масштабирования для отображения двух или более точек на карте

Мы хотим отобразить несколько маркеров на статической карте и рассчитать оптимальный уровень масштабирования, как это делают Карты Google. Мы уже вычислили ограничивающий прямоугольник и центральную точку карты, но теперь нам сложно вычислить правильный уровень масштабирования для отображения всего ограничивающего прямоугольника. Может кто-нибудь, пожалуйста, укажите нам правильное направление?

Чтобы получить уровень масштабирования, вам необходимо знать размеры вашей карты в пикселях. Вы также должны будете сделать свою математику в сферических координатах Меркатора.

1. Перевести широту, долготу в сферический меркатор x, y.
2. Получите расстояние между вашими двумя точками в сферическом меркаторе.
3. Проекция экватора составляет около 40 метров, а плитки - 256 пикселей, поэтому длина пикселя этой карты при заданном уровне масштабирования составляет примерно 256 * расстояние / 40000000 * 2 ^ увеличение . Попробуйте zoom = 0, zoom = 1, zoom = 2 до тех пор, пока расстояние не станет слишком длинным для измерений в пикселях вашей карты.

Это код C #, который я использую в 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: ограничивающий прямоугольник в координатах long / lat (x = long, y = lat)
• paddingFactor: это может быть использовано для получения эффекта "120%", на который ссылается ThomM. Значение 1,2 даст вам 120%.

Обратите внимание, что в моем случае zoomможет быть действительным числом. В случае веб-карт вам нужно целочисленное значение масштабирования, поэтому вы должны использовать что-то вроде (int)Math.Floor(zoom)этого.

Конечно, этот код относится только к проекции Web Mercator.

Если вы используете OpenLayers, то Map.getZoomForExtentрассчитаете максимальный уровень масштабирования, который может соответствовать всему экстенту на карте. В extentпотребности быть в проекции карты. Вы также можете использовать, fill_factorчтобы избежать отображения точек на краю карты и max_zoomограничить возможное увеличение:

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