WGS84 到 Google 地图位置和返回的 Java 代码

Question

搜索一些示例代码，用于将 WGS84 坐标系中的点转换为 Google 地图中的地图位置（像素位置），还支持缩放级别。

如果代码有很好的注释，那么它也可以采用其他语言。

您还可以向我指出一个开源 Java 项目:)

找到的一些资源：

OpenLayer 实现。

JOSM 项目

优秀 来自 JH LABS 的 Java 地图投影库 。 这是一个纯 java PROJ.4 端口。 从 WGS84 投影到米。 从这里开始，将米转换为平铺像素就非常简单了。

Answer 1

有人从 Google 地图中获取了 javascript 代码并将其移植到 python 中： gmerc.py

我用过这个，效果很好。

Answer 2

mapki 上的 Java 平铺实用程序代码。 com（谷歌地图开发人员的重要资源）

Answer 3

以下是 JavaSCript 中的函数...从 OpenLayers 中

function toMercator (lon, lat) {
  var x = lon * 20037508.34 / 180;
  var y = Math.log(Math.tan((90 + lat) * Math.PI / 360)) / (Math.PI / 180);
  y = y * 20037508.34 / 180;

  return [x, y];
  }

function inverseMercator (x, y) {
  var lon = (x / 20037508.34) * 180;
  var lat = (y / 20037508.34) * 180;

  lat = 180/Math.PI * (2 * Math.atan(Math.exp(lat * Math.PI / 180)) - Math.PI / 2);

  return [lon, lat];
  }

提取 转换为 Java 相当简单

Answer 4

GeoTools 有代码可以转换为从你能想象到的任何坐标系，其中也包括谷歌地图的坐标系。 它也是开源的。 然而，还应该指出的是，GeoTools 是一个大型库，因此，如果您正在寻找小型、快速且简单的东西，那么它可能不是最佳选择。

如果您还打算进行其他 GIS/坐标转换等，我会强烈推荐它。

如果您使用 GeoTools 或类似工具，您可能还会有兴趣了解 Google 地图坐标系称为 EPSG 3785。

Answer 5

我将其移植到 PHP - 这是代码，如果有人需要的话：

到 Mercator：

$lon = ($lon * 20037508.34) / 180;
$lat = log(tan((90 + $lat) * M_PI / 360)) / (M_PI / 180);
$lat = $lat * 20037508.34 / 180;

从 Mercator：

$lon = ($lon / 20037508.34) * 180;
$lat = ($lat / 20037508.34) * 180;
$lat = 180/M_PI * (2 * atan(exp($lat * M_PI / 180)) - M_PI / 2);

Answer 6

/*
 * Utility functions to transform between wgs84 and google projection coordinates
 * Derived from openmap http://openmap.bbn.com/
 */

public class MercatorTransform {
    public final static double NORTH_POLE = 90.0;
    public final static double SOUTH_POLE = -NORTH_POLE;
    public final static double DATELINE = 180.0;
    public final static double LON_RANGE = 360.0;

    final public static transient double wgs84_earthEquatorialRadiusMeters_D = 6378137.0;
    private static double latfac = wgs84_earthEquatorialRadiusMeters_D;
    private static double lonfac = wgs84_earthEquatorialRadiusMeters_D;

    final public static transient double HALF_PI_D = Math.PI / 2.0d;

    /**
     * Returns google projection coordinates from wgs84 lat,long coordinates
     */
    public static double[] forward(double lat, double lon) {

        lat = normalizeLatitude(lat);
        lon = wrapLongitude(lon);

        double latrad = Math.toRadians(lat);
        double lonrad = Math.toRadians(lon);

        double lat_m = latfac * Math.log(Math.tan(((latrad + HALF_PI_D) / 2d)));
        double lon_m = lonfac * lonrad;

        double[] x = { lon_m, lat_m };
        return x;
    }

    /**
     * Returns wgs84 lat,long coordinates from google projection coordinates
     */
    public static float[] inverse(float lon_m, float lat_m) {
        double latrad = (2d * Math.atan(Math.exp(lat_m / latfac))) - HALF_PI_D;
        double lonrad = lon_m / lonfac;

        double lat = Math.toDegrees(latrad);
        double lon = Math.toDegrees(lonrad);

        lat = normalizeLatitude(lat);
        lon = wrapLongitude(lon);
        float[] x = { (float) lat, (float) lon };

        return x;
    }

    private static double wrapLongitude(double lon) {
        if ((lon < -DATELINE) || (lon > DATELINE)) {
            lon += DATELINE;
            lon = lon % LON_RANGE;
            lon = (lon < 0) ? DATELINE + lon : -DATELINE + lon;
        }
        return lon;
    }

    private static double normalizeLatitude(double lat) {
        if (lat > NORTH_POLE) {
            lat = NORTH_POLE;
        }
        if (lat < SOUTH_POLE) {
            lat = SOUTH_POLE;
        }
        return lat;
    }

}