From 5f2c45ad3d1b5e877f5fc04da33102e63816644a Mon Sep 17 00:00:00 2001
From: Jonathan Guinet <jonathan.guinet@c-s.fr>
Date: Wed, 12 Jul 2017 14:39:22 -0400
Subject: [PATCH] [TEST] estimateFreeParameters tests removed.
---
.../org/orekit/rugged/api/RuggedTest.java | 741 ++++++------------
1 file changed, 256 insertions(+), 485 deletions(-)
diff --git a/src/test/java/org/orekit/rugged/api/RuggedTest.java b/src/test/java/org/orekit/rugged/api/RuggedTest.java
index 600887db..47c2717e 100644
--- a/src/test/java/org/orekit/rugged/api/RuggedTest.java
+++ b/src/test/java/org/orekit/rugged/api/RuggedTest.java
@@ -26,7 +26,6 @@ import java.net.URISyntaxException;
import java.nio.MappedByteBuffer;
import java.nio.channels.FileChannel;
import java.util.ArrayList;
-import java.util.Arrays;
import java.util.Collections;
import java.util.List;
import java.util.Locale;
@@ -38,7 +37,6 @@ import org.hipparchus.analysis.differentiation.UnivariateDifferentiableFunction;
import org.hipparchus.geometry.euclidean.threed.Rotation;
import org.hipparchus.geometry.euclidean.threed.RotationConvention;
import org.hipparchus.geometry.euclidean.threed.Vector3D;
-import org.hipparchus.optim.nonlinear.vector.leastsquares.LeastSquaresOptimizer.Optimum;
import org.hipparchus.stat.descriptive.StreamingStatistics;
import org.hipparchus.stat.descriptive.rank.Percentile;
import org.hipparchus.util.FastMath;
@@ -83,7 +81,6 @@ import org.orekit.utils.IERSConventions;
import org.orekit.utils.ParameterDriver;
import org.orekit.utils.TimeStampedAngularCoordinates;
import org.orekit.utils.TimeStampedPVCoordinates;
-import org.orekit.rugged.refining.measures.SensorToGroundMapping;
public class RuggedTest {
@@ -95,7 +92,7 @@ public class RuggedTest {
@Ignore
@Test
public void testMayonVolcanoTiming()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
long t0 = System.currentTimeMillis();
int dimension = 2000;
@@ -108,11 +105,11 @@ public class RuggedTest {
// Mayon Volcano location according to Wikipedia: 13°15′24″N 123°41′6″E
GeodeticPoint summit =
- new GeodeticPoint(FastMath.toRadians(13.25667), FastMath.toRadians(123.685), 2463.0);
+ new GeodeticPoint(FastMath.toRadians(13.25667), FastMath.toRadians(123.685), 2463.0);
VolcanicConeElevationUpdater updater =
- new VolcanicConeElevationUpdater(summit,
- FastMath.toRadians(30.0), 16.0,
- FastMath.toRadians(1.0), 1201);
+ new VolcanicConeElevationUpdater(summit,
+ FastMath.toRadians(30.0), 16.0,
+ FastMath.toRadians(1.0), 1201);
final AbsoluteDate crossing = new AbsoluteDate("2012-01-06T02:27:16.139", TimeScalesFactory.getUTC());
// one line sensor
@@ -135,15 +132,15 @@ public class RuggedTest {
Propagator ephemeris = propagator.getGeneratedEphemeris();
Rugged rugged = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(lineDatation.getDate(firstLine), lineDatation.getDate(lastLine), 0.001, 5.0).
- setTrajectory(1.0 / lineDatation.getRate(0.0),
- 8, CartesianDerivativesFilter.USE_PV,AngularDerivativesFilter.USE_R,
- ephemeris).
- addLineSensor(lineSensor).
- build();
+ setDigitalElevationModel(updater, 8).
+ setAlgorithm(AlgorithmId.DUVENHAGE).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(lineDatation.getDate(firstLine), lineDatation.getDate(lastLine), 0.001, 5.0).
+ setTrajectory(1.0 / lineDatation.getRate(0.0),
+ 8, CartesianDerivativesFilter.USE_PV,AngularDerivativesFilter.USE_R,
+ ephemeris).
+ addLineSensor(lineSensor).
+ build();
try {
@@ -173,10 +170,10 @@ public class RuggedTest {
int sizeM = size / (1024 * 1024);
System.out.format(Locale.US,
"%n%n%5dx%5d:%n" +
- " Orekit initialization and DEM creation : %5.1fs%n" +
- " direct location and %3dM grid writing: %5.1fs (%.1f px/s)%n",
- lastLine - firstLine, los.getNbPixels(),
- 1.0e-3 * (t1 - t0), sizeM, 1.0e-3 * (t2 - t1), pixels / (1.0e-3 * (t2 - t1)));
+ " Orekit initialization and DEM creation : %5.1fs%n" +
+ " direct location and %3dM grid writing: %5.1fs (%.1f px/s)%n",
+ lastLine - firstLine, los.getNbPixels(),
+ 1.0e-3 * (t1 - t0), sizeM, 1.0e-3 * (t2 - t1), pixels / (1.0e-3 * (t2 - t1)));
} catch (IOException ioe) {
Assert.fail(ioe.getLocalizedMessage());
}
@@ -185,7 +182,7 @@ public class RuggedTest {
@Test
public void testLightTimeCorrection()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
int dimension = 400;
@@ -212,15 +209,15 @@ public class RuggedTest {
AbsoluteDate maxDate = lineSensor.getDate(lastLine);
RuggedBuilder builder = new RuggedBuilder().
- setAlgorithm(AlgorithmId.IGNORE_DEM_USE_ELLIPSOID).
- setEllipsoid(EllipsoidId.IERS2003, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- addLineSensor(lineSensor);
+ setAlgorithm(AlgorithmId.IGNORE_DEM_USE_ELLIPSOID).
+ setEllipsoid(EllipsoidId.IERS2003, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ addLineSensor(lineSensor);
Rugged rugged = builder.build();
Assert.assertEquals(1, rugged.getLineSensors().size());
@@ -263,7 +260,7 @@ public class RuggedTest {
@Test
public void testAberrationOfLightCorrection()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
int dimension = 400;
@@ -291,17 +288,17 @@ public class RuggedTest {
RuggedBuilder builder = new RuggedBuilder().
- setAlgorithm(AlgorithmId.IGNORE_DEM_USE_ELLIPSOID).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- addLineSensor(lineSensor).
- setLightTimeCorrection(false).
- setAberrationOfLightCorrection(true);
+ setAlgorithm(AlgorithmId.IGNORE_DEM_USE_ELLIPSOID).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ addLineSensor(lineSensor).
+ setLightTimeCorrection(false).
+ setAberrationOfLightCorrection(true);
Rugged rugged = builder.build();
GeodeticPoint[] gpWithAberrationOfLightCorrection = rugged.directLocation("line", 200);
@@ -321,7 +318,7 @@ public class RuggedTest {
@Test
public void testFlatBodyCorrection()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
int dimension = 200;
@@ -350,24 +347,24 @@ public class RuggedTest {
AbsoluteDate maxDate = lineSensor.getDate(lastLine);
TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
- FastMath.toRadians(1.0), 257);
+ new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
+ FastMath.toRadians(1.0), 257);
RuggedBuilder builder = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- addLineSensor(lineSensor);
+ setDigitalElevationModel(updater, 8).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ addLineSensor(lineSensor);
GeodeticPoint[] gpWithFlatBodyCorrection =
- builder.setAlgorithm(AlgorithmId.DUVENHAGE).build().directLocation("line", 100);
+ builder.setAlgorithm(AlgorithmId.DUVENHAGE).build().directLocation("line", 100);
GeodeticPoint[] gpWithoutFlatBodyCorrection =
- builder.setAlgorithm(AlgorithmId.DUVENHAGE_FLAT_BODY).build().directLocation("line", 100);
+ builder.setAlgorithm(AlgorithmId.DUVENHAGE_FLAT_BODY).build().directLocation("line", 100);
StreamingStatistics stats = new StreamingStatistics();
for (int i = 0; i < gpWithFlatBodyCorrection.length; ++i) {
@@ -383,7 +380,7 @@ public class RuggedTest {
@Test
public void testLocationsinglePoint()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
int dimension = 200;
@@ -412,26 +409,26 @@ public class RuggedTest {
AbsoluteDate maxDate = lineSensor.getDate(lastLine);
TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
- FastMath.toRadians(1.0), 257);
+ new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
+ FastMath.toRadians(1.0), 257);
Rugged rugged = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- addLineSensor(lineSensor).build();
+ setDigitalElevationModel(updater, 8).
+ setAlgorithm(AlgorithmId.DUVENHAGE).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ addLineSensor(lineSensor).build();
GeodeticPoint[] gpLine = rugged.directLocation("line", 100);
for (int i = 0; i < gpLine.length; ++i) {
GeodeticPoint gpPixel =
- rugged.directLocation(lineSensor.getDate(100), lineSensor.getPosition(),
- lineSensor.getLOS(lineSensor.getDate(100), i));
+ rugged.directLocation(lineSensor.getDate(100), lineSensor.getPosition(),
+ lineSensor.getLOS(lineSensor.getDate(100), i));
Assert.assertEquals(gpLine[i].getLatitude(), gpPixel.getLatitude(), 1.0e-10);
Assert.assertEquals(gpLine[i].getLongitude(), gpPixel.getLongitude(), 1.0e-10);
Assert.assertEquals(gpLine[i].getAltitude(), gpPixel.getAltitude(), 1.0e-10);
@@ -441,7 +438,7 @@ public class RuggedTest {
@Test
public void testLocationsinglePointNoCorrections()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
int dimension = 200;
@@ -470,29 +467,29 @@ public class RuggedTest {
AbsoluteDate maxDate = lineSensor.getDate(lastLine);
TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
- FastMath.toRadians(1.0), 257);
+ new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
+ FastMath.toRadians(1.0), 257);
Rugged rugged = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- setAberrationOfLightCorrection(false).
- setLightTimeCorrection(false).
- addLineSensor(lineSensor).
- build();
+ setDigitalElevationModel(updater, 8).
+ setAlgorithm(AlgorithmId.DUVENHAGE).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ setAberrationOfLightCorrection(false).
+ setLightTimeCorrection(false).
+ addLineSensor(lineSensor).
+ build();
GeodeticPoint[] gpLine = rugged.directLocation("line", 100);
for (int i = 0; i < gpLine.length; ++i) {
GeodeticPoint gpPixel =
- rugged.directLocation(lineSensor.getDate(100), lineSensor.getPosition(),
- lineSensor.getLOS(lineSensor.getDate(100), i));
+ rugged.directLocation(lineSensor.getDate(100), lineSensor.getPosition(),
+ lineSensor.getLOS(lineSensor.getDate(100), i));
Assert.assertEquals(gpLine[i].getLatitude(), gpPixel.getLatitude(), 1.0e-10);
Assert.assertEquals(gpLine[i].getLongitude(), gpPixel.getLongitude(), 1.0e-10);
Assert.assertEquals(gpLine[i].getAltitude(), gpPixel.getAltitude(), 1.0e-10);
@@ -502,7 +499,7 @@ public class RuggedTest {
@Test
public void testBasicScan()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
int dimension = 200;
@@ -531,24 +528,24 @@ public class RuggedTest {
AbsoluteDate maxDate = lineSensor.getDate(lastLine);
TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xe12ef744f224cf43l,
- FastMath.toRadians(1.0), 257);
+ new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xe12ef744f224cf43l,
+ FastMath.toRadians(1.0), 257);
RuggedBuilder builder = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- addLineSensor(lineSensor);
+ setDigitalElevationModel(updater, 8).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ addLineSensor(lineSensor);
GeodeticPoint[] gpDuvenhage =
- builder.setAlgorithm(AlgorithmId.DUVENHAGE).build().directLocation("line", 100);
+ builder.setAlgorithm(AlgorithmId.DUVENHAGE).build().directLocation("line", 100);
GeodeticPoint[] gpBasicScan =
- builder.setAlgorithm(AlgorithmId.BASIC_SLOW_EXHAUSTIVE_SCAN_FOR_TESTS_ONLY).build().directLocation("line", 100);
+ builder.setAlgorithm(AlgorithmId.BASIC_SLOW_EXHAUSTIVE_SCAN_FOR_TESTS_ONLY).build().directLocation("line", 100);
double[] data = new double[gpDuvenhage.length];
for (int i = 0; i < gpDuvenhage.length; ++i) {
@@ -565,7 +562,7 @@ public class RuggedTest {
@Ignore
@Test
public void testInverseLocationTiming()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
long t0 = System.currentTimeMillis();
int dimension = 2000;
@@ -600,19 +597,19 @@ public class RuggedTest {
AbsoluteDate maxDate = sensors.get(sensors.size() - 1).getDate(lastLine).shiftedBy(+1.0);
TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
- FastMath.toRadians(1.0), 257);
+ new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
+ FastMath.toRadians(1.0), 257);
RuggedBuilder builder = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R);
+ setDigitalElevationModel(updater, 8).
+ setAlgorithm(AlgorithmId.DUVENHAGE).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R);
builder.setLightTimeCorrection(true);
builder.setAberrationOfLightCorrection(true);
for (LineSensor lineSensor : sensors) {
@@ -662,12 +659,12 @@ public class RuggedTest {
int sizeM = (int) (size / (1024l * 1024l));
System.out.format(Locale.US,
"%n%n%5dx%5d, %d sensors:%n" +
- " Orekit initialization and DEM creation : %5.1fs%n" +
- " inverse location and %3dM grid writing: %5.1fs (%.1f px/s, %.1f%% covered)%n",
- dimension, dimension, nbSensors,
- 1.0e-3 * (t1 - t0), sizeM, 1.0e-3 * (t2 - t1),
- (badPixels + goodPixels) / (1.0e-3 * (t2 - t1)),
- (100.0 * goodPixels) / (goodPixels + badPixels));
+ " Orekit initialization and DEM creation : %5.1fs%n" +
+ " inverse location and %3dM grid writing: %5.1fs (%.1f px/s, %.1f%% covered)%n",
+ dimension, dimension, nbSensors,
+ 1.0e-3 * (t1 - t0), sizeM, 1.0e-3 * (t2 - t1),
+ (badPixels + goodPixels) / (1.0e-3 * (t2 - t1)),
+ (100.0 * goodPixels) / (goodPixels + badPixels));
} catch (IOException ioe) {
Assert.fail(ioe.getLocalizedMessage());
}
@@ -675,7 +672,7 @@ public class RuggedTest {
@Test
public void testInverseLocation()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
checkInverseLocation(2000, false, false, 4.0e-7, 5.0e-6);
checkInverseLocation(2000, false, true, 1.0e-5, 2.0e-7);
checkInverseLocation(2000, true, false, 4.0e-7, 4.0e-7);
@@ -684,16 +681,16 @@ public class RuggedTest {
@Test
public void testDateLocation()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
checkDateLocation(2000, false, false, 7.0e-7);
checkDateLocation(2000, false, true, 2.0e-5);
checkDateLocation(2000, true, false, 8.0e-7);
checkDateLocation(2000, true, true, 3.0e-6);
}
-
+
@Test
public void testLineDatation()
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
checkLineDatation(2000, 7.0e-7);
checkLineDatation(10000, 8.0e-7);
}
@@ -737,13 +734,13 @@ public class RuggedTest {
TileUpdater updater = new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l, FastMath.toRadians(1.0), 257);
RuggedBuilder builder = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(satellitePVList.get(0).getDate(), satellitePVList.get(satellitePVList.size() - 1).getDate(), 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- satellitePVList, 8, CartesianDerivativesFilter.USE_PV,
- satelliteQList, 8, AngularDerivativesFilter.USE_R);
+ setDigitalElevationModel(updater, 8).
+ setAlgorithm(AlgorithmId.DUVENHAGE).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(satellitePVList.get(0).getDate(), satellitePVList.get(satellitePVList.size() - 1).getDate(), 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ satellitePVList, 8, CartesianDerivativesFilter.USE_PV,
+ satelliteQList, 8, AngularDerivativesFilter.USE_R);
List<Vector3D> lineOfSight = new ArrayList<Vector3D>();
lineOfSight.add(new Vector3D(-0.011204, 0.181530, 1d).normalize());
@@ -875,15 +872,15 @@ public class RuggedTest {
LineSensor lineSensor = new LineSensor("QUICK_LOOK", lineDatation, offset,
new LOSBuilder(lineOfSight).build());
Rugged rugged = new RuggedBuilder().
- setAlgorithm(AlgorithmId.IGNORE_DEM_USE_ELLIPSOID).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(satellitePVList.get(0).getDate(),
- satellitePVList.get(satellitePVList.size() - 1).getDate(), 0.1, 10.0).
- setTrajectory(InertialFrameId.EME2000,
- satellitePVList, 6, CartesianDerivativesFilter.USE_P,
- satelliteQList, 8, AngularDerivativesFilter.USE_R).
- addLineSensor(lineSensor).
- build();
+ setAlgorithm(AlgorithmId.IGNORE_DEM_USE_ELLIPSOID).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(satellitePVList.get(0).getDate(),
+ satellitePVList.get(satellitePVList.size() - 1).getDate(), 0.1, 10.0).
+ setTrajectory(InertialFrameId.EME2000,
+ satellitePVList, 6, CartesianDerivativesFilter.USE_P,
+ satelliteQList, 8, AngularDerivativesFilter.USE_R).
+ addLineSensor(lineSensor).
+ build();
GeodeticPoint[] temp = rugged.directLocation("QUICK_LOOK", -250);
GeodeticPoint first = temp[0];
@@ -903,7 +900,7 @@ public class RuggedTest {
@Test
public void testInverseLocCurvedLine()
- throws RuggedException, URISyntaxException, OrekitException {
+ throws RuggedException, URISyntaxException, OrekitException {
String path = getClass().getClassLoader().getResource("orekit-data").toURI().getPath();
DataProvidersManager.getInstance().addProvider(new DirectoryCrawler(new File(path)));
@@ -931,21 +928,21 @@ public class RuggedTest {
AbsoluteDate maxDate = lineSensor.getDate(lastLine).shiftedBy(+1.0);
TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.3, 0xf0a401650191f9f6l,
- FastMath.toRadians(1.0), 257);
+ new RandomLandscapeUpdater(0.0, 9000.0, 0.3, 0xf0a401650191f9f6l,
+ FastMath.toRadians(1.0), 257);
Rugged rugged = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- addLineSensor(lineSensor).
- build();
+ setDigitalElevationModel(updater, 8).
+ setAlgorithm(AlgorithmId.DUVENHAGE).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ addLineSensor(lineSensor).
+ build();
int lineNumber = 97;
GeodeticPoint[] gp = rugged.directLocation("curved", lineNumber);
@@ -959,7 +956,7 @@ public class RuggedTest {
private void checkInverseLocation(int dimension, boolean lightTimeCorrection, boolean aberrationOfLightCorrection,
double maxLineError, double maxPixelError)
- throws RuggedException, OrekitException, URISyntaxException {
+ throws RuggedException, OrekitException, URISyntaxException {
String path = getClass().getClassLoader().getResource("orekit-data").toURI().getPath();
DataProvidersManager.getInstance().addProvider(new DirectoryCrawler(new File(path)));
@@ -987,23 +984,23 @@ public class RuggedTest {
AbsoluteDate maxDate = lineSensor.getDate(lastLine).shiftedBy(+1.0);
TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.3, 0xf0a401650191f9f6l,
- FastMath.toRadians(1.0), 257);
+ new RandomLandscapeUpdater(0.0, 9000.0, 0.3, 0xf0a401650191f9f6l,
+ FastMath.toRadians(1.0), 257);
Rugged rugged = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- setLightTimeCorrection(lightTimeCorrection).
- setAberrationOfLightCorrection(aberrationOfLightCorrection).
- addLineSensor(lineSensor).
- build();
+ setDigitalElevationModel(updater, 8).
+ setAlgorithm(AlgorithmId.DUVENHAGE).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ setLightTimeCorrection(lightTimeCorrection).
+ setAberrationOfLightCorrection(aberrationOfLightCorrection).
+ addLineSensor(lineSensor).
+ build();
double referenceLine = 0.87654 * dimension;
GeodeticPoint[] gp = rugged.directLocation("line", referenceLine);
@@ -1012,67 +1009,67 @@ public class RuggedTest {
int i = (int) FastMath.floor(p);
double d = p - i;
SensorPixel sp = rugged.inverseLocation("line",
- (1 - d) * gp[i].getLatitude() + d * gp[i + 1].getLatitude(),
- (1 - d) * gp[i].getLongitude() + d * gp[i + 1].getLongitude(),
- 0, dimension);
+ (1 - d) * gp[i].getLatitude() + d * gp[i + 1].getLatitude(),
+ (1 - d) * gp[i].getLongitude() + d * gp[i + 1].getLongitude(),
+ 0, dimension);
Assert.assertEquals(referenceLine, sp.getLineNumber(), maxLineError);
Assert.assertEquals(p, sp.getPixelNumber(), maxPixelError);
}
// point out of line (20 pixels before first pixel)
Assert.assertNull(rugged.inverseLocation("line",
- 21 * gp[0].getLatitude() - 20 * gp[1].getLatitude(),
- 21 * gp[0].getLongitude() - 20 * gp[1].getLongitude(),
- 0, dimension));
+ 21 * gp[0].getLatitude() - 20 * gp[1].getLatitude(),
+ 21 * gp[0].getLongitude() - 20 * gp[1].getLongitude(),
+ 0, dimension));
// point out of line (20 pixels after last pixel)
Assert.assertNull(rugged.inverseLocation("line",
- -20 * gp[gp.length - 2].getLatitude() + 21 * gp[gp.length - 1].getLatitude(),
- -20 * gp[gp.length - 2].getLongitude() + 21 * gp[gp.length - 1].getLongitude(),
- 0, dimension));
+ -20 * gp[gp.length - 2].getLatitude() + 21 * gp[gp.length - 1].getLatitude(),
+ -20 * gp[gp.length - 2].getLongitude() + 21 * gp[gp.length - 1].getLongitude(),
+ 0, dimension));
// point out of line (20 lines before first line)
GeodeticPoint[] gp0 = rugged.directLocation("line", 0);
GeodeticPoint[] gp1 = rugged.directLocation("line", 1);
Assert.assertNull(rugged.inverseLocation("line",
- 21 * gp0[dimension / 2].getLatitude() - 20 * gp1[dimension / 2].getLatitude(),
- 21 * gp0[dimension / 2].getLongitude() - 20 * gp1[dimension / 2].getLongitude(),
- 0, dimension));
+ 21 * gp0[dimension / 2].getLatitude() - 20 * gp1[dimension / 2].getLatitude(),
+ 21 * gp0[dimension / 2].getLongitude() - 20 * gp1[dimension / 2].getLongitude(),
+ 0, dimension));
// point out of line (20 lines after last line)
GeodeticPoint[] gp2 = rugged.directLocation("line", dimension - 2);
GeodeticPoint[] gp3 = rugged.directLocation("line", dimension - 1);
Assert.assertNull(rugged.inverseLocation("line",
- -20 * gp2[dimension / 2].getLatitude() + 21 * gp3[dimension / 2].getLatitude(),
- -20 * gp2[dimension / 2].getLongitude() + 21 * gp3[dimension / 2].getLongitude(),
- 0, dimension));
+ -20 * gp2[dimension / 2].getLatitude() + 21 * gp3[dimension / 2].getLatitude(),
+ -20 * gp2[dimension / 2].getLongitude() + 21 * gp3[dimension / 2].getLongitude(),
+ 0, dimension));
}
@Test
public void testInverseLocationDerivativesWithoutCorrections()
- throws RuggedException, OrekitException {
+ throws RuggedException, OrekitException {
doTestInverseLocationDerivatives(2000, false, false,
8.0e-9, 3.0e-10, 5.0e-12, 9.0e-8);
}
@Test
public void testInverseLocationDerivativesWithLightTimeCorrection()
- throws RuggedException, OrekitException {
+ throws RuggedException, OrekitException {
doTestInverseLocationDerivatives(2000, true, false,
3.0e-9, 9.0e-9, 2.1e-12, 9.0e-8);
}
@Test
public void testInverseLocationDerivativesWithAberrationOfLightCorrection()
- throws RuggedException, OrekitException {
+ throws RuggedException, OrekitException {
doTestInverseLocationDerivatives(2000, false, true,
4.2e-10, 3.0e-10, 3.4e-12, 7.0e-8);
}
@Test
public void testInverseLocationDerivativesWithAllCorrections()
- throws RuggedException, OrekitException {
+ throws RuggedException, OrekitException {
doTestInverseLocationDerivatives(2000, true, true,
3.0e-10, 5.0e-10, 2.0e-12, 7.0e-8);
}
@@ -1084,7 +1081,7 @@ public class RuggedTest {
double pixelTolerance,
double lineDerivativeRelativeTolerance,
double pixelDerivativeRelativeTolerance)
- throws RuggedException, OrekitException {
+ throws RuggedException, OrekitException {
try {
String path = getClass().getClassLoader().getResource("orekit-data").toURI().getPath();
@@ -1099,11 +1096,11 @@ public class RuggedTest {
// los: swath in the (YZ) plane, looking at 50° roll, 2.6" per pixel
Vector3D position = new Vector3D(1.5, 0, -0.2);
LOSBuilder losBuilder =
- TestUtils.createLOSPerfectLine(new Rotation(Vector3D.PLUS_I,
- FastMath.toRadians(50.0),
- RotationConvention.VECTOR_OPERATOR).applyTo(Vector3D.PLUS_K),
- Vector3D.PLUS_I,
- FastMath.toRadians(dimension * 2.6 / 3600.0), dimension);
+ TestUtils.createLOSPerfectLine(new Rotation(Vector3D.PLUS_I,
+ FastMath.toRadians(50.0),
+ RotationConvention.VECTOR_OPERATOR).applyTo(Vector3D.PLUS_K),
+ Vector3D.PLUS_I,
+ FastMath.toRadians(dimension * 2.6 / 3600.0), dimension);
losBuilder.addTransform(new FixedRotation("roll", Vector3D.MINUS_I, 0.0));
losBuilder.addTransform(new FixedRotation("pitch", Vector3D.MINUS_J, 0.0));
TimeDependentLOS los = losBuilder.build();
@@ -1234,247 +1231,21 @@ public class RuggedTest {
Assert.assertEquals(dXdP, result[1].getPartialDerivative(0, 1), dXdP * pixelDerivativeRelativeTolerance);
} catch (InvocationTargetException | NoSuchMethodException |
- SecurityException | IllegalAccessException |
- IllegalArgumentException | URISyntaxException |
- OrekitExceptionWrapper | RuggedExceptionWrapper e) {
+ SecurityException | IllegalAccessException |
+ IllegalArgumentException | URISyntaxException |
+ OrekitExceptionWrapper | RuggedExceptionWrapper e) {
Assert.fail(e.getLocalizedMessage());
}
}
- @Test
- public void testNoParametersSelected()
- throws OrekitException {
- try {
- Rugged perfect = createParameterEstimationContext(2000, false, false).build();
- perfect.estimateFreeParameters(Collections.singletonList(new SensorToGroundMapping("line")),
- 100, 1.0e-3);
- Assert.fail("an exception should have been thrown");
- } catch (RuggedException re) {
- Assert.assertEquals(RuggedMessages.NO_PARAMETERS_SELECTED, re.getSpecifier());
- }
- }
-
- @Test
- public void testDuplicatedParameter()
- throws OrekitException {
- try {
- TimeDependentLOS los =
- TestUtils.createLOSPerfectLine(Vector3D.PLUS_K, Vector3D.PLUS_I, 0.1, 1000).
- addTransform(new FixedRotation("pitch", Vector3D.MINUS_J, 0.0)).build();
- LineDatation lineDatation = new LinearLineDatation(AbsoluteDate.J2000_EPOCH, 50, 1.0 / 1.5e-3);
- Rugged perfect = createParameterEstimationContext(2000, false, false).
- addLineSensor(new LineSensor("other", lineDatation, Vector3D.ZERO, los)).
- build();
- // the duplicated parameter is "pitch", it will be detected even despite only "roll" is estimated
- perfect.getLineSensor("line").
- getParametersDrivers().
- filter(driver -> driver.getName().equals("roll")).
- forEach(driver -> driver.setSelected(true));
- perfect.estimateFreeParameters(Arrays.asList(new SensorToGroundMapping("line"),
- new SensorToGroundMapping("other")),
- 100, 1.0e-3);
- Assert.fail("an exception should have been thrown");
- } catch (RuggedException re) {
- Assert.assertEquals(RuggedMessages.DUPLICATED_PARAMETER_NAME, re.getSpecifier());
- Assert.assertEquals("pitch", re.getParts()[0]);
- }
- }
-
- @Test
- public void testNoReferenceMapping()
- throws OrekitException {
- try {
- Rugged perfect = createParameterEstimationContext(2000, false, false).build();
- perfect.getLineSensor("line").
- getParametersDrivers().
- filter(driver -> driver.getName().equals("roll") || driver.getName().equals("pitch")).
- forEach(driver -> driver.setSelected(true));
- perfect.estimateFreeParameters(Collections.singletonList(new SensorToGroundMapping("line")),
- 100, 1.0e-3);
- Assert.fail("an exception should have been thrown");
- } catch (RuggedException re) {
- Assert.assertEquals(RuggedMessages.NO_REFERENCE_MAPPINGS, re.getSpecifier());
- }
- }
-
- @Test
- public void testViewingModelParametersEstimationWithoutCorrections()
- throws OrekitException, RuggedException {
- doTestViewingModelParametersEstimation(2000, false, false,
- FastMath.toRadians( 0.01), 1.0e-8,
- FastMath.toRadians(-0.03), 2.0e-9,
- 6);
- }
-
- @Test
- public void testViewingModelParametersEstimationWithLightTimeCorrection()
- throws OrekitException, RuggedException {
- doTestViewingModelParametersEstimation(2000, true, false,
- FastMath.toRadians( 0.01), 1.0e-8,
- FastMath.toRadians(-0.03), 2.0e-9,
- 6);
- }
-
- @Test
- public void testViewingModelParametersEstimationWithAberrationOfLightCorrection()
- throws OrekitException, RuggedException {
- doTestViewingModelParametersEstimation(2000, false, true,
- FastMath.toRadians( 0.01), 1.0e-8,
- FastMath.toRadians(-0.03), 2.0e-9,
- 6);
- }
-
- @Test
- public void testViewingModelParametersEstimationWithAllCorrections()
- throws OrekitException, RuggedException {
- doTestViewingModelParametersEstimation(2000, true, true,
- FastMath.toRadians( 0.01), 1.0e-8,
- FastMath.toRadians(-0.03), 3.0e-9,
- 6);
- }
-
- private void doTestViewingModelParametersEstimation(int dimension,
- boolean lightTimeCorrection,
- boolean aberrationOfLightCorrection,
- double rollValue, double rollTolerance,
- double pitchValue, double pitchTolerance,
- int expectedEvaluations)
- throws OrekitException, RuggedException {
- try {
-
- // set up a perfect Rugged instance, with prescribed roll and pitch values
- Rugged perfect = createParameterEstimationContext(dimension,
- lightTimeCorrection,
- aberrationOfLightCorrection).build();
- perfect.
- getLineSensor("line").
- getParametersDrivers().
- filter(driver -> driver.getName().equals("roll")).
- findFirst().get().setValue(rollValue);
-
- perfect.
- getLineSensor("line").
- getParametersDrivers().
- filter(driver -> driver.getName().equals("pitch")).
- findFirst().get().setValue(pitchValue);
-
- // generate reference mapping
- SensorToGroundMapping mapping = new SensorToGroundMapping("line");
- LineSensor sensor = perfect.getLineSensor(mapping.getSensorName());
- for (double line = 0; line < dimension; line += 100) {
- AbsoluteDate date = sensor.getDate(line);
- for (int pixel = 0; pixel < sensor.getNbPixels(); pixel += 100) {
- GeodeticPoint gp = perfect.directLocation(date, sensor.getPosition(),
- sensor.getLOS(date, pixel));
- mapping.addMapping(new SensorPixel(line, pixel), gp);
- }
- }
-
- // set up a completely independent Rugged instance,
- // with 0.0 roll and pitch values and estimating rool and pitch
- Rugged normal = createParameterEstimationContext(dimension,
- lightTimeCorrection,
- aberrationOfLightCorrection).build();
- normal.
- getLineSensor("line").
- getParametersDrivers().
- filter(driver -> driver.getName().equals("roll") || driver.getName().equals("pitch")).
- forEach(driver -> {
- try {
- driver.setSelected(true);
- driver.setValue(0.0);
- } catch (OrekitException e) {
- throw new OrekitExceptionWrapper(e);
- }
- });
-
- // perform parameters estimation
- Optimum optimum = normal.estimateFreeParameters(Collections.singletonList(mapping), 10, 1.0e-9);
- Assert.assertEquals(expectedEvaluations, optimum.getEvaluations());
-
- // check estimated values
- double estimatedRoll = normal.getLineSensor("line").
- getParametersDrivers().
- filter(driver -> driver.getName().equals("roll")).
- findFirst().get().getValue();
- Assert.assertEquals("roll error = " + (estimatedRoll - rollValue),
- rollValue, estimatedRoll, rollTolerance);
-
- double estimatedPitch = normal.getLineSensor("line").
- getParametersDrivers().
- filter(driver -> driver.getName().equals("pitch")).
- findFirst().get().getValue();
- Assert.assertEquals("pitch error = " + (estimatedPitch - pitchValue),
- pitchValue, estimatedPitch, pitchTolerance);
-
- } catch (OrekitExceptionWrapper oew) {
- throw oew.getException();
- }
- }
-
- private RuggedBuilder createParameterEstimationContext(int dimension,
- boolean lightTimeCorrection,
- boolean aberrationOfLightCorrection)
- throws OrekitException, RuggedException {
-
- try {
-
- String path = getClass().getClassLoader().getResource("orekit-data").toURI().getPath();
- DataProvidersManager.getInstance().addProvider(new DirectoryCrawler(new File(path)));
- final BodyShape earth = TestUtils.createEarth();
- final Orbit orbit = TestUtils.createOrbit(Constants.EIGEN5C_EARTH_MU);
-
- AbsoluteDate crossing = new AbsoluteDate("2012-01-01T12:30:00.000", TimeScalesFactory.getUTC());
-
- // one line sensor
- // position: 1.5m in front (+X) and 20 cm above (-Z) of the S/C center of mass
- // los: swath in the (YZ) plane, looking at 50° roll, 2.6" per pixel
- Vector3D position = new Vector3D(1.5, 0, -0.2);
- LOSBuilder losBuilder =
- TestUtils.createLOSPerfectLine(new Rotation(Vector3D.PLUS_I,
- FastMath.toRadians(50.0),
- RotationConvention.VECTOR_OPERATOR).applyTo(Vector3D.PLUS_K),
- Vector3D.PLUS_I,
- FastMath.toRadians(dimension * 2.6 / 3600.0), dimension);
- losBuilder.addTransform(new FixedRotation("roll", Vector3D.MINUS_I, 0.0));
- losBuilder.addTransform(new FixedRotation("pitch", Vector3D.MINUS_J, 0.0));
- TimeDependentLOS los = losBuilder.build();
-
- // linear datation model: at reference time we get the middle line, and the rate is one line every 1.5ms
- LineDatation lineDatation = new LinearLineDatation(crossing, dimension / 2, 1.0 / 1.5e-3);
- int firstLine = 0;
- int lastLine = dimension;
- final LineSensor lineSensor = new LineSensor("line", lineDatation, position, los);
- AbsoluteDate minDate = lineSensor.getDate(firstLine).shiftedBy(-1.0);
- AbsoluteDate maxDate = lineSensor.getDate(lastLine).shiftedBy(+1.0);
- TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.3, 0xf0a401650191f9f6l,
- FastMath.toRadians(1.0), 257);
- return new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- setLightTimeCorrection(lightTimeCorrection).
- setAberrationOfLightCorrection(aberrationOfLightCorrection).
- addLineSensor(lineSensor);
+ /*TODO add refining tests */
- } catch (URISyntaxException e) {
- Assert.fail(e.getLocalizedMessage());
- return null; // never reached
- }
- }
private void checkDateLocation(int dimension, boolean lightTimeCorrection, boolean aberrationOfLightCorrection,
- double maxDateError)
- throws RuggedException, OrekitException, URISyntaxException {
+ double maxDateError)
+ throws RuggedException, OrekitException, URISyntaxException {
String path = getClass().getClassLoader().getResource("orekit-data").toURI().getPath();
DataProvidersManager.getInstance().addProvider(new DirectoryCrawler(new File(path)));
@@ -1502,23 +1273,23 @@ public class RuggedTest {
AbsoluteDate maxDate = lineSensor.getDate(lastLine).shiftedBy(+1.0);
TileUpdater updater =
- new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
- FastMath.toRadians(1.0), 257);
+ new RandomLandscapeUpdater(0.0, 9000.0, 0.5, 0xf0a401650191f9f6l,
+ FastMath.toRadians(1.0), 257);
Rugged rugged = new RuggedBuilder().
- setDigitalElevationModel(updater, 8).
- setAlgorithm(AlgorithmId.DUVENHAGE).
- setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
- setTimeSpan(minDate, maxDate, 0.001, 5.0).
- setTrajectory(InertialFrameId.EME2000,
- TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 8, CartesianDerivativesFilter.USE_PV,
- TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
- 2, AngularDerivativesFilter.USE_R).
- setLightTimeCorrection(lightTimeCorrection).
- setAberrationOfLightCorrection(aberrationOfLightCorrection).
- addLineSensor(lineSensor).
- build();
+ setDigitalElevationModel(updater, 8).
+ setAlgorithm(AlgorithmId.DUVENHAGE).
+ setEllipsoid(EllipsoidId.WGS84, BodyRotatingFrameId.ITRF).
+ setTimeSpan(minDate, maxDate, 0.001, 5.0).
+ setTrajectory(InertialFrameId.EME2000,
+ TestUtils.orbitToPV(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 8, CartesianDerivativesFilter.USE_PV,
+ TestUtils.orbitToQ(orbit, earth, minDate.shiftedBy(-1.0), maxDate.shiftedBy(+1.0), 0.25),
+ 2, AngularDerivativesFilter.USE_R).
+ setLightTimeCorrection(lightTimeCorrection).
+ setAberrationOfLightCorrection(aberrationOfLightCorrection).
+ addLineSensor(lineSensor).
+ build();
double referenceLine = 0.87654 * dimension;
GeodeticPoint[] gp = rugged.directLocation("line", referenceLine);
@@ -1527,9 +1298,9 @@ public class RuggedTest {
int i = (int) FastMath.floor(p);
double d = p - i;
AbsoluteDate date = rugged.dateLocation("line",
- (1 - d) * gp[i].getLatitude() + d * gp[i + 1].getLatitude(),
- (1 - d) * gp[i].getLongitude() + d * gp[i + 1].getLongitude(),
- 0, dimension);
+ (1 - d) * gp[i].getLatitude() + d * gp[i + 1].getLatitude(),
+ (1 - d) * gp[i].getLongitude() + d * gp[i + 1].getLongitude(),
+ 0, dimension);
Assert.assertEquals(0.0, date.durationFrom(lineSensor.getDate(referenceLine)), maxDateError);
}
@@ -1537,52 +1308,52 @@ public class RuggedTest {
GeodeticPoint[] gp0 = rugged.directLocation("line", 0);
GeodeticPoint[] gp1 = rugged.directLocation("line", 1);
Assert.assertNull(rugged.dateLocation("line",
- 21 * gp0[dimension / 2].getLatitude() - 20 * gp1[dimension / 2].getLatitude(),
- 21 * gp0[dimension / 2].getLongitude() - 20 * gp1[dimension / 2].getLongitude(),
- 0, dimension));
+ 21 * gp0[dimension / 2].getLatitude() - 20 * gp1[dimension / 2].getLatitude(),
+ 21 * gp0[dimension / 2].getLongitude() - 20 * gp1[dimension / 2].getLongitude(),
+ 0, dimension));
// point out of line (20 lines after last line)
GeodeticPoint[] gp2 = rugged.directLocation("line", dimension - 2);
GeodeticPoint[] gp3 = rugged.directLocation("line", dimension - 1);
Assert.assertNull(rugged.dateLocation("line",
- -20 * gp2[dimension / 2].getLatitude() + 21 * gp3[dimension / 2].getLatitude(),
- -20 * gp2[dimension / 2].getLongitude() + 21 * gp3[dimension / 2].getLongitude(),
- 0, dimension));
-
+ -20 * gp2[dimension / 2].getLatitude() + 21 * gp3[dimension / 2].getLatitude(),
+ -20 * gp2[dimension / 2].getLongitude() + 21 * gp3[dimension / 2].getLongitude(),
+ 0, dimension));
+
}
private void checkLineDatation(int dimension, double maxLineError)
- throws RuggedException, OrekitException, URISyntaxException {
-
- String path = getClass().getClassLoader().getResource("orekit-data").toURI().getPath();
- DataProvidersManager.getInstance().addProvider(new DirectoryCrawler(new File(path)));
-
- AbsoluteDate crossing = new AbsoluteDate("2012-01-01T12:30:00.000", TimeScalesFactory.getUTC());
-
- // one line sensor
- // position: 1.5m in front (+X) and 20 cm above (-Z) of the S/C center of mass
- // los: swath in the (YZ) plane, looking at 50° roll, 2.6" per pixel
- Vector3D position = new Vector3D(1.5, 0, -0.2);
- TimeDependentLOS los = TestUtils.createLOSPerfectLine(new Rotation(Vector3D.PLUS_I,
- FastMath.toRadians(50.0),
- RotationConvention.VECTOR_OPERATOR).applyTo(Vector3D.PLUS_K),
- Vector3D.PLUS_I,
- FastMath.toRadians(dimension * 2.6 / 3600.0), dimension).build();
-
- // linear datation model: at reference time we get the middle line, and the rate is one line every 1.5ms
- LineDatation lineDatation = new LinearLineDatation(crossing, dimension / 2, 1.0 / 1.5e-3);
- int firstLine = 0;
- int lastLine = dimension;
- LineSensor lineSensor = new LineSensor("line", lineDatation, position, los);
- AbsoluteDate minDate = lineSensor.getDate(firstLine).shiftedBy(-1.0);
- AbsoluteDate maxDate = lineSensor.getDate(lastLine).shiftedBy(+1.0);
-
- // Recompute the lines from the date with the appropriate shift of date
- double recomputedFirstLine = lineSensor.getLine(minDate.shiftedBy(+1.0));
- double recomputedLastLine = lineSensor.getLine(maxDate.shiftedBy(-1.0));
-
- Assert.assertEquals(firstLine, recomputedFirstLine, maxLineError);
- Assert.assertEquals(lastLine, recomputedLastLine, maxLineError);
+ throws RuggedException, OrekitException, URISyntaxException {
+
+ String path = getClass().getClassLoader().getResource("orekit-data").toURI().getPath();
+ DataProvidersManager.getInstance().addProvider(new DirectoryCrawler(new File(path)));
+
+ AbsoluteDate crossing = new AbsoluteDate("2012-01-01T12:30:00.000", TimeScalesFactory.getUTC());
+
+ // one line sensor
+ // position: 1.5m in front (+X) and 20 cm above (-Z) of the S/C center of mass
+ // los: swath in the (YZ) plane, looking at 50° roll, 2.6" per pixel
+ Vector3D position = new Vector3D(1.5, 0, -0.2);
+ TimeDependentLOS los = TestUtils.createLOSPerfectLine(new Rotation(Vector3D.PLUS_I,
+ FastMath.toRadians(50.0),
+ RotationConvention.VECTOR_OPERATOR).applyTo(Vector3D.PLUS_K),
+ Vector3D.PLUS_I,
+ FastMath.toRadians(dimension * 2.6 / 3600.0), dimension).build();
+
+ // linear datation model: at reference time we get the middle line, and the rate is one line every 1.5ms
+ LineDatation lineDatation = new LinearLineDatation(crossing, dimension / 2, 1.0 / 1.5e-3);
+ int firstLine = 0;
+ int lastLine = dimension;
+ LineSensor lineSensor = new LineSensor("line", lineDatation, position, los);
+ AbsoluteDate minDate = lineSensor.getDate(firstLine).shiftedBy(-1.0);
+ AbsoluteDate maxDate = lineSensor.getDate(lastLine).shiftedBy(+1.0);
+
+ // Recompute the lines from the date with the appropriate shift of date
+ double recomputedFirstLine = lineSensor.getLine(minDate.shiftedBy(+1.0));
+ double recomputedLastLine = lineSensor.getLine(maxDate.shiftedBy(-1.0));
+
+ Assert.assertEquals(firstLine, recomputedFirstLine, maxLineError);
+ Assert.assertEquals(lastLine, recomputedLastLine, maxLineError);
}
}
--
GitLab