用java实现根据轨道六根数数据计算卫星运行轨迹

模拟卫星运行轨迹有很多种实现方法，可以通过纯前端的方式实现，也可通过后端服务进行计算。

     纯前端实现方式：通过引入satellite.js库来计算卫星的轨迹，再结合Cesium来模拟显示卫星的轨迹。前端实现代码过会再贴图，这节主要介绍后端如何通过JAVA来计算卫星轨迹。

    采用JAVA计算卫星轨迹需要引入Orekit库和导入相关包和文件，具体方式可以参见上篇文章CSDNhttps://mp.csdn.net/mp_blog/creation/editor/135748302

计算轨迹核心代码

 /**
     *
     * @param semiMajorAxis 半长轴
     * @param eccentricity  偏心率
     * @param inclination     倾角
     * @param raan    升交点赤径
     * @param pa     近地点幅角
     * @param trueAnomaly   真近点角
     * @param startDate   开始时间
     * @param endDate     结束时间
     */
    public void predictSatelliteOrbit(Double semiMajorAxis, Double eccentricity, Double inclination,
            Double raan, Double pa, Double trueAnomaly, AbsoluteDate startDate, AbsoluteDate endDate){
        Frame icrfFrame = FramesFactory.getICRF();
        Orbit orbit = new KeplerianOrbit(semiMajorAxis, eccentricity, inclination, raan, pa, trueAnomaly,
                PositionAngleType.TRUE, icrfFrame, new AbsoluteDate(), Constants.WGS84_EARTH_MU);
        double stepSize = 60.0; // 步长为60s,计算每隔一分钟中的位置
        KeplerianPropagator propagator = new KeplerianPropagator(orbit);
        for(AbsoluteDate date=startDate;date.compareTo(endDate)<=0;date = date.shiftedBy(stepSize)){
            PVCoordinates pvCoordinates = propagator.getPVCoordinates(date,icrfFrame);
            Vector3D position = pvCoordinates.getPosition();
            System.out.println("日期："+date+"位置:"+position);
        }
    }

测试代码生成结果：