java地理坐标系及投影间转换代码示例

public double[] sevenParamTrans(Datum7Paras datum7Paras,double X,double Y,double Z) {
        double[] point = new double[3];
        double dx = datum7Paras.getDx();
        double dy = datum7Paras.getDy();
        double dz = datum7Paras.getDz();
        double rx = datum7Paras.getRx() * 0.0000048481373323; // 1秒=0.0000048481373323 弧度
        double ry = datum7Paras.getRy() * 0.0000048481373323;
        double rz = datum7Paras.getRz() * 0.0000048481373323;
        double m = datum7Paras.getPpm();
        double k = m / 1000000;

        point[0] = (1 + k) * (X + this.rz * Y - this.ry * Z) + this.dx;
        point[1] = (1 + k) * (-this.rz * X + Y + this.rx * Z) + this.dy;
        point[2] = (1 + k) * (this.ry * X - this.rx * Y + Z) + this.dz;
}

class Datum7Paras{
    private double dx;
    private double dy;
    private double dz;
    private double rx;
    private double ry;
    private double rz;
    private double ppm;

    public Datum7Paras(double dx, double dy, double dz, double rx, double ry, double rz, double ppm) {
        this.dx = dx;
        this.dy = dy;
        this.dz = dz;
        this.rx = rx;
        this.ry = ry;
        this.rz = rz;
        this.ppm = ppm;
    }

    public double getDx() {
        return dx;
    }

    public void setDx(double dx) {
        this.dx = dx;
    }

    public double getDy() {
        return dy;
    }

    public void setDy(double dy) {
        this.dy = dy;
    }

    public double getDz() {
        return dz;
    }

    public void setDz(double dz) {
        this.dz = dz;
    }

    public double getRx() {
        return rx;
    }

    public void setRx(double rx) {
        this.rx = rx;
    }

    public double getRy() {
        return ry;
    }

    public void setRy(double ry) {
        this.ry = ry;
    }

    public double getRz() {
        return rz;
    }

    public void setRz(double rz) {
        this.rz = rz;
    }

    public double getPpm() {
        return ppm;
    }

    public void setPpm(double ppm) {
        this.ppm = ppm;
    }
}

public double[] transform4Para(Trans4Paras transPara,double X,double Y) {
        double X1 = transPara.getDx();
        double Y1 = transPara.getDy();

        double cosAngle = Math.cos(transPara.getA());
        double sinAngle = Math.sin(transPara.getA());

        double[] point = new double[2];
        point[0] = X1 + transPara.getK() * (cosAngle * X - sinAngle * Y);
        point[1] = Y1 + transPara.getK() * (sinAngle * X + cosAngle * Y);
    }

class Trans4Paras{
    private double _dx;
    private double _dy;
    private double _a;
    private double _k;
  
    public double getDx() {
        return _dx;
    }

    public void setDx(double dx) {
        _dx = dx;
    }

    public double getDy() {
        return _dy;
    }

    public void setDy(double dy) {
        _dy = dy;
    }

    public double getA() {
        return _a;
    }

    public void setA(double a) {
        _a = a;
    }

    public double getK() {
        return _k;
    }

    public void setK(double k) {
        _k = k;
    }

    public Trans4Paras(double dx, double dy, double a, double k) {
        _dx = dx;
        _dy = dy;
        _a = a;
        _k = k;
    }

}

<dependency>
            <groupId>org.osgeo</groupId>
            <artifactId>proj4j</artifactId>
            <version>0.1.0</version>
</dependency>

import org.osgeo.proj4j.*;

public class PointTest {
    public static void main(String[] args) {
        ToCoordinate(579573.572, 4269411.894,"EPSG:32651","EPSG:4326",true);
        ToCoordinate(123.91343418163903, 38.569600138042816,"32651","4326",false);
    }

    public static void ToCoordinate(double argX,double argY,String sys1,String sys2,boolean xyTolb){
        // 创建CRS工厂
        CRSFactory crsFactory = new CRSFactory();

        // 创建UTM坐标参考系统（CRS），这里以UTM 33N为例
        CoordinateReferenceSystem utmCrs = crsFactory.createFromName("EPSG:" + sys1);

        // 创建WGS84坐标参考系统
        CoordinateReferenceSystem wgs84Crs = crsFactory.createFromName("EPSG:" + sys2);

        // 创建坐标转换工厂
        CoordinateTransformFactory transformFactory = new CoordinateTransformFactory();
        if (xyTolb){
            // 创建从UTM到WGS84的坐标转换器
            CoordinateTransform transformToWgs84 = transformFactory.createTransform(utmCrs, wgs84Crs);
            // 创建UTM坐标点
            ProjCoordinate utmPoint = new ProjCoordinate(argX, argY); // X, Y in meters

            // 转换到WGS84坐标点
            ProjCoordinate wgs84Point = new ProjCoordinate();
            transformToWgs84.transform(utmPoint, wgs84Point);

            // 打印WGS84坐标
            System.out.println("WGS84: " + wgs84Point.x + ", " + wgs84Point.y);
        }else {
            // 如果需要，可以创建从WGS84到UTM的坐标转换器
            CoordinateTransform transformToUtm = transformFactory.createTransform(wgs84Crs,utmCrs);
            // 创建UTM坐标点
            ProjCoordinate wgs84Point = new ProjCoordinate(argX, argY);
            // 转换回UTM坐标点
            ProjCoordinate utmPointBack = new ProjCoordinate();
            transformToUtm.transform(wgs84Point, utmPointBack);

            // 打印转换回的UTM坐标
            System.out.println("UTM: " + utmPointBack.x + ", " + utmPointBack.y);
        }

    }
}