Examples with Attitude gaiasky.util.gaia.Attitude used on opensource projects

Example 1 with Attitude

use of gaiasky.util.gaia.Attitude in project gaiasky by langurmonkey.

the class GaiaAttitudeWriter method attitudeZAxisToFile.

public static void attitudeZAxisToFile(OutputType type) {
    Calendar iniCal = GregorianCalendar.getInstance();
    // 2014-01-15T15:43:04
    iniCal.set(2014, Calendar.FEBRUARY, 15, 15, 43, 4);
    Date ini = iniCal.getTime();
    // 2019-06-20T06:20:05
    iniCal.set(2019, Calendar.JULY, 20, 5, 40, 0);
    Date end = iniCal.getTime();
    long tini = ini.getTime();
    long tend = end.getTime();
    Vector3d v1 = new Vector3d();
    Vector3d v2 = new Vector3d();
    Vector3d sph1 = new Vector3d();
    Vector3d sph2 = new Vector3d();
    String filenameEcl = "Gaia-ecliptic-Z-" + System.currentTimeMillis() + ".csv";
    String filenameEq = "Gaia-equatorial-Z-" + System.currentTimeMillis() + ".csv";
    File fecl = new File(System.getProperty("java.io.tmpdir") + File.separator + filenameEcl);
    File feq = new File(System.getProperty("java.io.tmpdir") + File.separator + filenameEq);
    try {
        Writer writerEcl = new BufferedWriter(new OutputStreamWriter(new FileOutputStream(fecl), StandardCharsets.UTF_8));
        Writer writerEq = new BufferedWriter(new OutputStreamWriter(new FileOutputStream(feq), StandardCharsets.UTF_8));
        if (type.equals(OutputType.UP_VECTOR)) {
            writerEcl.append("#EclLong[deg], EclLat[deg], t[ms-January_1_1970_00:00:00_GMT], current[date], attitudeFile\n");
            writerEq.append("#alpha[deg], delta[deg], t[ms-January_1_1970_00:00:00_GMT], current[date], attitudeFile\n");
        } else if (type.equals(OutputType.FOV_VECTORS)) {
            writerEcl.append("#EclLong-fov1[deg], EclLat-fov1[deg], EclLong-fov2[deg], EclLat-fov2[deg],t[ms-January_1_1970_00:00:00_GMT], current[date], attitudeFile\n");
            writerEq.append("#alpha-fov1[deg], delta-fov1[deg], alpha-fov2[deg], delta-fov2[deg],t[ms-January_1_1970_00:00:00_GMT], current[date], attitudeFile\n");
        }
        // 10 minutes
        long step = 10 * 60000;
        for (long t = tini; t <= tend; t += step) {
            Date current = new Date(t);
            try {
                Attitude att = GaiaAttitudeServer.instance.getAttitude(current);
                Quaterniond quat = att.getQuaternion();
                if (type.equals(OutputType.UP_VECTOR)) {
                    // Up vector in Gaia coordinates
                    v1.set(0, 1, 0);
                    // Equatorial
                    v1.rotateVectorByQuaternion(quat);
                    Coordinates.cartesianToSpherical(v1, sph1);
                    writerEq.append(String.valueOf(Math.toDegrees(sph1.x))).append(", ").append(String.valueOf(Math.toDegrees(sph1.y))).append(", ").append(String.valueOf(t)).append(", ").append(String.valueOf(current)).append(", ").append(GaiaAttitudeServer.instance.getCurrentAttitudeName()).append("\n");
                    // Ecliptic
                    v1.mul(Coordinates.eqToEcl());
                    Coordinates.cartesianToSpherical(v1, sph1);
                    writerEcl.append(String.valueOf(Math.toDegrees(sph1.x))).append(", ").append(String.valueOf(Math.toDegrees(sph1.y))).append(", ").append(String.valueOf(t)).append(", ").append(String.valueOf(current)).append(", ").append(GaiaAttitudeServer.instance.getCurrentAttitudeName()).append("\n");
                } else if (type.equals(OutputType.FOV_VECTORS)) {
                    // Fov 1 and Fov 2
                    v1.set(0, 0, 1).rotate(BAM_2, 0, 1, 0);
                    v2.set(0, 0, 1).rotate(-BAM_2, 0, 1, 0);
                    // Equatorial
                    v1.rotateVectorByQuaternion(quat);
                    v2.rotateVectorByQuaternion(quat);
                    Coordinates.cartesianToSpherical(v1, sph1);
                    Coordinates.cartesianToSpherical(v2, sph2);
                    writerEq.append(String.valueOf(Math.toDegrees(sph1.x))).append(", ").append(String.valueOf(Math.toDegrees(sph1.y))).append(", ").append(String.valueOf(Math.toDegrees(sph2.x))).append(", ").append(String.valueOf(Math.toDegrees(sph2.y))).append(", ").append(String.valueOf(t)).append(", ").append(String.valueOf(current)).append(", ").append(GaiaAttitudeServer.instance.getCurrentAttitudeName()).append("\n");
                    // Ecliptic
                    v1.mul(Coordinates.eqToEcl());
                    v2.mul(Coordinates.eqToEcl());
                    Coordinates.cartesianToSpherical(v1, sph1);
                    Coordinates.cartesianToSpherical(v2, sph2);
                    writerEcl.append(String.valueOf(Math.toDegrees(sph1.x))).append(", ").append(String.valueOf(Math.toDegrees(sph1.y))).append(", ").append(String.valueOf(Math.toDegrees(sph2.x))).append(", ").append(String.valueOf(Math.toDegrees(sph2.y))).append(", ").append(String.valueOf(t)).append(", ").append(String.valueOf(current)).append(", ").append(GaiaAttitudeServer.instance.getCurrentAttitudeName()).append("\n");
                }
            } catch (Exception e) {
                System.out.println("Error: t=" + current);
            }
        }
        writerEcl.close();
        writerEq.close();
    } catch (IOException e) {
        e.printStackTrace();
    }
}