Общие формулы проекции Mercator для карт Google не работают правильно
Я создаю сервер наложения плитки для Google maps на C# и нашел несколько различных примеров кода для вычисления Y из широты. После того, как они начали работать в целом, я начал замечать некоторые случаи, когда накладки не выстраивались должным образом. Чтобы проверить это, я сделал тестовый жгут для сравнения преобразования Mercator LatToY Карты Google с формулами, которые я нашел в интернете. Как вы можете видеть ниже, они не совпадают в определенных случаях.
случае #1
Zoomed Out: проблема больше всего очевидно, когда увеличено. Вблизи, проблема едва заметна.
дело № 2
близость пункта к верхней части & дну просмотр границ: проблема хуже в середине границ просмотра, и становится лучше по краям. Такое поведение может свести на нет поведение в случае #1
Тест:
I создал страницу Google maps для отображение красных линий с помощью карты Google API встроен в преобразование Mercator, и наложите это с помощью изображения ссылочный код для выполнения Mercator преобразование. Эти преобразования представлены в виде черных линий. Сравнивать разница.
Результаты: экватор http://www.kayak411.com/Mercator/MercatorComparison%20-%20Equator.png Север Увеличен http://www.kayak411.com/Mercator/MercatorComparison%20-%20North%20Zoomed%20Out.png
Проверьте верхнюю и нижнюю линии: Северный верхний и Нижний пример http://www.kayak411.com/Mercator/MercatorComparison%20-%20North%20Zoomed%20Out%20-%20TopAndBottom.png
проблема становится визуально больше, но численно меньше, как увеличить: alt-текст http://www.kayak411.com/Mercator/MercatorComparison%20-%20North%20Zoomed%20Midway.png
и все это почти исчезает при более близких уровнях масштабирования, независимо от ориентации экрана. alt текст http://www.kayak411.com/Mercator/MercatorComparison%20-%20North%20Zoomed%20In.png
Код:
Код Стороны Клиента Google Maps:
var lat = 0;
for (lat = -80; lat <= 80; lat += 5) {
map.addOverlay(new GPolyline([new GLatLng(lat, -180), new GLatLng(lat, 0)], "#FF0033", 2));
map.addOverlay(new GPolyline([new GLatLng(lat, 0), new GLatLng(lat, 180)], "#FF0033", 2));
}
Код На Стороне Сервера:
Плитка Резак : http://mapki.com/wiki/Tile_Cutter
Вики OpenStreetMap : http://wiki.openstreetmap.org/wiki/Mercator
protected override void ImageOverlay_ComposeImage(ref Bitmap ZipCodeBitMap)
{
Graphics LinesGraphic = Graphics.FromImage(ZipCodeBitMap);
Int32 MapWidth = Convert.ToInt32(Math.Pow(2, zoom) * 255);
Point Offset =
Cartographer.Mercator2.toZoomedPixelCoords(North, West, zoom);
TrimPoint(ref Offset, MapWidth);
for (Double lat = -80; lat <= 80; lat += 5)
{
Point StartPoint = Cartographer.Mercator2.toZoomedPixelCoords(lat, -179, zoom);
Point EndPoint = Cartographer.Mercator2.toZoomedPixelCoords(lat, -1, zoom);
TrimPoint(ref StartPoint, MapWidth);
TrimPoint(ref EndPoint, MapWidth);
StartPoint.X = StartPoint.X - Offset.X;
EndPoint.X = EndPoint.X - Offset.X;
StartPoint.Y = StartPoint.Y - Offset.Y;
EndPoint.Y = EndPoint.Y - Offset.Y;
LinesGraphic.DrawLine(new Pen(Color.Black, 2),
StartPoint.X,
StartPoint.Y,
EndPoint.X,
EndPoint.Y);
LinesGraphic.DrawString(
lat.ToString(),
new Font("Verdana", 10),
new SolidBrush(Color.Black),
new Point(
Convert.ToInt32((width / 3.0) * 2.0),
StartPoint.Y));
}
}
protected void TrimPoint(ref Point point, Int32 MapWidth)
{
point.X = Math.Max(point.X, 0);
point.X = Math.Min(point.X, MapWidth - 1);
point.Y = Math.Max(point.Y, 0);
point.Y = Math.Min(point.Y, MapWidth - 1);
}
Итак, кто-нибудь когда-либо испытывал это? Осмелюсь спросить, разрешил ли это? Или просто иметь лучшую реализацию c# преобразования координат проекта Mercator?
спасибо!
2 ответов
спасибо всем за ваши предложения и помощь.
в конце концов я выяснил, что это не формула или техническая проблема, я считаю, что это проблема методологии.
вы не можете определить область просмотра в формате Lat/Lng и ожидать ее заполнения соответствующими проекциями Меркатора. Вот где происходит искажение. Вместо этого необходимо определить правильное поле просмотра в Mercator и project Mercator.
делая это, я смог правильно совпадают с Google maps.
возможно, вам придется создать несколько точек вдоль longtitude, чтобы точки были спроецированы корректно вдоль широты. В ваших примерах вы действительно проецируете только две точки в начале и конце линии и соединяете их.
проблема будет более очевидной на экваторе из-за более значительной кривизны Земли. Он будет меньше, если увеличить его по той же причине.
посмотреть http://code.google.com/apis/maps/documentation/overlays.html#Great_Circles
попробуйте создать свои полилинии Google с параметром geodsic, чтобы увидеть, имеет ли это значение. Я думаю, что это добавляет точки вдоль линии и проецирует их автоматически:
var lat = 0;
var polyOptions = {geodesic:true};
for (lat = -80; lat <= 80; lat += 5) {
map.addOverlay(new GPolyline([new GLatLng(lat, -180), new GLatLng(lat, 0)], "#FF0033", 2, polyOptions));
map.addOverlay(new GPolyline([new GLatLng(lat, 0), new GLatLng(lat, 180)], "#FF0033", 2, polyOptions));
}
Я должен был прочитать это, так как все мои измерения расстояния были неправильными в OpenLayers по аналогичным причинам:http://geographika.co.uk/watch-out-for-openlayer-distances (больше ссылки / пояснения)