Java代码如何实现经纬度转XY

实现经纬度转XY，需要用到投影算法，这里给出一个Java代码示例，使用的是墨卡托投影算法：

public class LatLon2Meters {
    private static final double EARTH_RADIUS = 6378137; // 地球半径
    private static final double MAX_LATITUDE = 85.0511287798; // 最大纬度

    public static double[] latLonToMeters(double latitude, double longitude) {
        double[] xy = new double[2];
        xy[0] = longitude * EARTH_RADIUS * Math.PI / 180; // X轴
        xy[1] = Math.log(Math.tan((90 + latitude) * Math.PI / 360)) / (Math.PI / 180);
        xy[1] = xy[1] * EARTH_RADIUS * Math.PI / 180; // Y轴
        return xy;
    }
}

其中，`latitude`和`longitude`分别代表要转换的经纬度，返回的`xy`是一个长度为2的数组，表示转换后的XY坐标。请注意，由于墨卡托投影算法是一个近似算法，因此转换后的坐标并不是精确的，但是可以满足大部分应用场景的需求。

相关问题

java代码实现cgcs2000大地坐标系xy值转化为对应经纬度

Java代码示例如下：

import org.apache.commons.math3.geometry.euclidean.threed.Vector3D;

public class CoordinateConverter {
    private static final double a = 6378137;  // 长半轴
    private static final double b = 6356752.3142;  // 短半轴
    private static final double f = (a - b) / a;  // 扁率
    private static final double e1 = Math.sqrt(2 * f - f * f);  // 第一偏心率
    private static final double e2 = Math.sqrt((a * a - b * b) / (b * b));  // 第二偏心率

    public static Vector3D convertXYToLatLon(double x, double y) {
        // 投影坐标转换为大地坐标
        double Bf0 = x / a;
        double Nf0 = (Math.pow(a, 2) / b) / Math.sqrt(1 + Math.pow(e2, 2) * Math.pow(Math.cos(Bf0), 2));
        double Mf0 = a * (1 - Math.pow(e1, 2)) / Math.pow(Math.sqrt(1 - Math.pow(e1, 2) * Math.pow(Math.sin(Bf0), 2)), 3);
        double nf = Math.sqrt(1 + Math.pow(e2, 2) * Math.pow(Math.cos(Bf0), 2));
        double etaf = e2 * Math.cos(Bf0);
        double T1f = Math.pow(Math.tan(Bf0), 2);
        double C1f = Math.pow(e1, 2) * Math.pow(Math.cos(Bf0), 2);
        double R = a * (1 - Math.pow(e1, 2))
                / Math.pow(Math.sqrt(1 - Math.pow(e1, 2) * Math.pow(Math.sin(Bf0), 2)), 2);

        double latitude = Bf0 - (Nf0 * Math.tan(Bf0) / Mf0)
                * (Math.pow(x, 2) / (2 * R))
                + (Nf0 * Math.tan(Bf0) / (24 * Mf0 * Mf0 * Mf0))
                * ((5 + 3 * T1f + 10 * C1f - 4 * C1f * C1f - 9 * etaf * etaf)
                * Math.pow(x, 4) / 24
                - (61 + 90 * T1f + 298 * C1f + 45 * T1f * T1f - 252 * etaf * etaf - 3 * C1f * C1f)
                * Math.pow(x, 6) / 720);
        double longitude = (x / (nf * Mf0 * Math.cos(Bf0))
                - (1 + 2 * T1f + C1f) * Math.pow(x, 3) / (6 * nf * nf * nf * Mf0 * Math.pow(Math.cos(Bf0), 3))
                + (5 + 28 * T1f + 24 * T1f * T1f + 6 * C1f + 8 * etaf * etaf)
                * Math.pow(x, 5) / (120 * Math.pow(nf, 5) * Mf0 * Math.pow(Math.cos(Bf0), 5)))
                / Math.cos(Bf0);

        return new Vector3D(Math.toDegrees(latitude), Math.toDegrees(longitude), y);
    }

    public static void main(String[] args) {
        double x = 1000000;  // x坐标
        double y = 2000000;  // y坐标

        Vector3D latLon = convertXYToLatLon(x, y);
        System.out.printf("经度：%f, 纬度：%f, 高程：%f", latLon.getY(), latLon.getX(), latLon.getZ());
    }
}

该代码实现了cgcs2000大地坐标系的xy值转化为对应经纬度。在代码中，首先定义了大地坐标系相关参数，如长半轴、短半轴、扁率、偏心率等。然后通过投影坐标转换为大地坐标的公式，计算出经度、纬度和高程，并将结果返回。在main方法中，我们可以指定xy坐标，然后调用转换函数，将结果输出为经度、纬度和高程。

地理信息 54经纬度坐标124.793263,46.26888194 转 54坐标

54坐标系是中国新的统一的大地坐标系，其投影方式是横轴为高斯-克吕格投影，纵轴为墨卡托投影。54坐标系是一种平面直角坐标系，以米为单位，与WGS84坐标系之间可以通过转换公式进行转换。下面是将WGS84坐标系中的经纬度坐标转换为54坐标系的Java代码示例：

public class CoordinateTransform {
    private static final double A = 6378137.0;   // 长半轴
    private static final double B = 6356752.3142;   // 短半轴
    private static final double L0 = 121.0;   // 中央经线

    private static double rad(double degree) {
        return degree * Math.PI / 180.0;
    }

    public static double[] WGS84To54(double lng, double lat) {
        double L = rad(lng);
        double B = rad(lat);
        double L0 = rad(CoordinateTransform.L0);
        double e2 = 1 - (B / A) * (B / A);
        double N = A / Math.sqrt(1 - e2 * Math.sin(B) * Math.sin(B));
        double T = Math.tan(B) * Math.tan(B);
        double C = e2 * Math.cos(B) * Math.cos(B);
        double A1 = (L - L0) * Math.cos(B);
        double A2 = (L - L0) * (L - L0) * Math.cos(B) * Math.cos(B) * Math.cos(B) * (1 - T + C) / 6.0;
        double A3 = (L - L0) * (L - L0) * (L - L0) * Math.cos(B) * Math.cos(B) * Math.cos(B) * Math.cos(B) * (5 - 18 * T + T * T + 72 * C - 58 * e2) / 120.0;
        double X = 1000000 + N * A1 + 500000 * Math.sin(B) * Math.cos(B) * (A2 + A3);
        double Y = 1000000 + N * B * Math.cos(B) * (1 + (A1 * A1) / 2.0 + (1 - T + C) * (A2 * A2) / 24.0 + (5 - 4 * T + 14 * C + 13 * C * C - 28 * T * C) * (A3 * A3) / 720.0);
        return new double[]{X, Y};
    }
}

使用方法：

double lng = 124.793263;
double lat = 46.26888194;
double[] xy = CoordinateTransform.WGS84To54(lng, lat);
System.out.println("54坐标系 X 坐标：" + xy[0]);
System.out.println("54坐标系 Y 坐标：" + xy[1]);

其中，`A` 和 `B` 分别为 WGS84 椭球体的长半轴和短半轴，`L0` 为中央经线经度，`e2` 为椭球体第一偏心率的平方，`N` 为子午线曲率半径，`T` 和 `C` 分别为一、二次纬度变形参数，`X` 和 `Y` 为计算得到的54坐标系的 X 和 Y 坐标。