火星坐标系 (GCJ-02) 与百度坐标系 (BD-09) 的转换算法(android)

    final double x_pi = 3.14159265358979324 * 3000.0 / 180.0;
    public GeoPoint bd_encrypt(GeoPoint point){
        double x = point.getLongitudeE6()/1E6;
        double y = point.getLatitudeE6()/1E6;

        double z = Math.sqrt(x*x+y*y)+0.00002*Math.sin(y*x_pi);
        double theta = Math.atan2(y, x)+0.000003*Math.cos(x*x_pi);

        GeoPoint bd_point = new GeoPoint((int) ((z*Math.sin(theta)+0.006)*1E6), (int) ((z*Math.cos(theta)+0.0065)*1E6));
        return bd_point;
    }

    public GeoPoint bd_decrypt(GeoPoint bd_point){
        double x = (double)bd_point.getLongitudeE6()/1E6 - 0.0065;
        double y = (double)bd_point.getLatitudeE6()/1E6 - 0.006;
        double z = Math.sqrt(x*x+y*y) - 0.00002*Math.sin(y*x_pi);
        double theta = Math.atan2(y, x) - 0.000003*Math.cos(x*x_pi);
        GeoPoint point = new GeoPoint((int)(z*Math.sin(theta)*1E6), (int)(z*Math.cos(theta)*1E6));
        return point;
    }