深圳地方坐标系转WGS84坐标系

接到公司需求，要求将第三方供应商提供的坐标转换为wgs84大地坐标（也就是我们平时说的gps坐标），上网搜了下资料，什么大地转空间，空间转平面，网上的干货晒得不太干，难受得很，没办法，自己学习一下吧，话不多说上代码（转换四参数不能瞎给我乱改了，方法是对的，结果我也出来了，参数的话自己去找吧）。

大致流程：深圳平面坐标→wgs84平面坐标→wgs84大地坐标（调用convert方法后再调用GaussToWGSBL方法）

	/**
     * 深圳平面转为wgs84大地坐标系
     */
    private static class SZToWGS84 {

        double a = 参数自己想办法;

        double b = 参数自己想办法;

        double t = 参数自己想办法;

        double k = 参数自己想办法;

        /**
         * convert
         *
         * @return
         */
        public double[] convert(double x, double y) {
            double[] coord = {0.0, 0.0};
            double X;
            double Y;
            X = a + x * k * Math.cos(t) - y * k * Math.sin(t);
            Y = b + y * k * Math.cos(t) + x * k * Math.sin(t);
            coord[0] = X;
            coord[1] = Y;
            return coord;
        }

        /**
         * 由高斯投影坐标反算成经纬度 84坐标系
         *
         * @param X
         * @param Y
         * @return
         */
        public double[] GaussToWGSBL(double X, double Y) {

            int ProjNo;
            // 带宽
            int ZoneWide;
            double[] output = new double[2];
            // latitude0,
            double longitude1, latitude1, longitude0, X0, Y0, xval, yval;
            double e1, e2, f, a, ee, NN, T, C, M, D, R, u, fai, iPI;
            //3.1415926535898/180.0;
            iPI = 0.0174532925199433;
            a = 6378137;
            f = 1 / 298.2572236;
            // 3度带宽
            ZoneWide = 3;
            // 查找带号
            ProjNo = (int) (X / 1000000);
            // 中央经线
            longitude0 = ZoneWide * 38 * Math.PI / 180;
            X0 = ProjNo * 1000000L + 500000L;
            Y0 = 0;
            xval = X - X0;
            yval = Y - Y0; //带内大地坐标
            e2 = 2 * f - f * f;
            e1 = (1.0 - Math.sqrt(1 - e2)) / (1.0 + Math.sqrt(1 - e2));
            ee = e2 / (1 - e2);
            M = yval;
            u = M / (a * (1 - e2 / 4 - 3 * e2 * e2 / 64 - 5 * e2 * e2 * e2 / 256));
            fai = u + (3 * e1 / 2 - 27 * e1 * e1 * e1 / 32) * Math.sin(2 * u) + (21 * e1 * e1 / 16 - 55 * e1 * e1 * e1 * e1 / 32) * Math.sin(
                    4 * u) + (151 * e1 * e1 * e1 / 96) * Math.sin(6 * u) + (1097 * e1 * e1 * e1 * e1 / 512) * Math.sin(8 * u);
            C = ee * Math.cos(fai) * Math.cos(fai);
            T = Math.tan(fai) * Math.tan(fai);
            NN = a / Math.sqrt(1.0 - e2 * Math.sin(fai) * Math.sin(fai));
            R = a * (1 - e2) / Math.sqrt((1 - e2 * Math.sin(fai) * Math.sin(fai)) * (1 - e2 * Math.sin(fai) * Math.sin(fai)) * (1 - e2 * Math.sin
                    (fai) * Math.sin(fai)));
            D = xval / NN;
            //计算经度(Longitude) 纬度(Latitude)
            longitude1 = longitude0 + (D - (1 + 2 * T + C) * D * D * D / 6 + (5 - 2 * C + 28 * T - 3 * C * C + 8 * ee + 24 * T * T) * D
                    * D * D * D * D / 120) / Math.cos(fai);
            latitude1 = fai - (NN * Math.tan(fai) / R) * (D * D / 2 - (5 + 3 * T + 10 * C - 4 * C * C - 9 * ee) * D * D * D * D / 24
                    + (61 + 90 * T + 298 * C + 45 * T * T - 256 * ee - 3 * C * C) * D * D * D * D * D * D / 720);
            // 转换为度 DD
            output[0] = longitude1 / iPI;
            output[1] = latitude1 / iPI;
            return output;
        }
    }

    public static void main(String[] args) {
        SZToWGS84 szToWgs84 = new SZToWGS84();
        double[] coord = szToWgs84.convert(20880.461, 94340.704);
        System.out.println(coord[0] + "," + coord[1]);
        double[] output = szToWgs84.GaussToWGSBL(coord[1],coord[0]);
        System.out.println(output[0] + "," + output[1]);
    }

方法是对的，我自己已经在用了，我的参数第三方给的不准确，精确度大概在±10m，勉强够用了。（PS：如果是其他地方坐标系转wgs84，更改参数应该可以复用）

这里参考了一些其他博文

坐标转换流程与公式 七参数 四参数
高斯投影正算与反算的理论方法与实现代码