Class LMedSRobustEasyGyroscopeCalibrator
java.lang.Object
com.irurueta.navigation.inertial.calibration.gyroscope.RobustEasyGyroscopeCalibrator
com.irurueta.navigation.inertial.calibration.gyroscope.LMedSRobustEasyGyroscopeCalibrator
- All Implemented Interfaces:
AccelerometerDependentGyroscopeCalibrator
,GyroscopeCalibrator
,GyroscopeNonLinearCalibrator
,OrderedBodyKinematicsSequenceGyroscopeCalibrator
,QualityScoredGyroscopeCalibrator
,UnknownBiasGyroscopeCalibrator
,GyroscopeBiasUncertaintySource
,GyroscopeCalibrationSource
Robustly estimates gyroscope biases, cross couplings and scaling factors
along with G-dependent cross biases introduced on the gyroscope by the
specific forces sensed by the accelerometer using LMedS robust estimator.
This calibrator assumes that the IMU is at a more or less fixed location on Earth, and evaluates sequences of measured body kinematics to perform calibration for unknown orientations on those provided sequences. Each provided sequence will be preceded by a static period where mean specific force will be measured to determine gravity (and hence partial body attitude).
Measured gyroscope angular rates is assumed to follow the model shown below:
Ωmeas = bg + (I + Mg) * Ωtrue + Gg * ftrue + wWhere: - Ωmeas is the measured gyroscope angular rates. This is a 3x1 vector. - bg is the gyroscope bias. Ideally, on a perfect gyroscope, this should be a 3x1 zero vector. - I is the 3x3 identity matrix. - Mg is the 3x3 matrix containing cross-couplings and scaling factors. Ideally, on a perfect gyroscope, this should be a 3x3 zero matrix. - Ωtrue is ground-truth gyroscope angular rates. - Gg is the G-dependent cross biases introduced by the specific forces sensed by the accelerometer. Ideally, on a perfect gyroscope, this should be a 3x3 zero matrix. - ftrue is ground-truth specific force. This is a 3x1 vector. - w is measurement noise. This is a 3x1 vector.
-
Nested Class Summary
Nested classes/interfaces inherited from class com.irurueta.navigation.inertial.calibration.gyroscope.RobustEasyGyroscopeCalibrator
RobustEasyGyroscopeCalibrator.PreliminaryResult
-
Field Summary
FieldsModifier and TypeFieldDescriptionstatic final double
Default value to be used for stop threshold.static final double
Minimum allowed stop threshold value.private double
Threshold to be used to keep the algorithm iterating in case that best estimated threshold using median of residuals is not small enough.Fields inherited from class com.irurueta.navigation.inertial.calibration.gyroscope.RobustEasyGyroscopeCalibrator
confidence, DEFAULT_CONFIDENCE, DEFAULT_ESTIMATE_G_DEPENDENT_CROSS_BIASES, DEFAULT_KEEP_COVARIANCE, DEFAULT_MAX_ITERATIONS, DEFAULT_PROGRESS_DELTA, DEFAULT_REFINE_RESULT, DEFAULT_ROBUST_METHOD, DEFAULT_USE_COMMON_Z_AXIS, inliersData, listener, MAX_CONFIDENCE, MAX_PROGRESS_DELTA, maxIterations, MIN_CONFIDENCE, MIN_ITERATIONS, MIN_PROGRESS_DELTA, preliminarySubsetSize, progressDelta, refineResult, running, sequences
-
Constructor Summary
ConstructorsConstructorDescriptionConstructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, double[] accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, double[] accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa, RobustEasyGyroscopeCalibratorListener listener) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, RobustEasyGyroscopeCalibratorListener listener) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, com.irurueta.algebra.Matrix accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, com.irurueta.algebra.Matrix accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa, RobustEasyGyroscopeCalibratorListener listener) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, RobustEasyGyroscopeCalibratorListener listener) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, double[] accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, double[] accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa, RobustEasyGyroscopeCalibratorListener listener) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, RobustEasyGyroscopeCalibratorListener listener) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, com.irurueta.algebra.Matrix accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, com.irurueta.algebra.Matrix accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa, RobustEasyGyroscopeCalibratorListener listener) Constructor.LMedSRobustEasyGyroscopeCalibrator
(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, RobustEasyGyroscopeCalibratorListener listener) Constructor. -
Method Summary
Modifier and TypeMethodDescriptionvoid
Estimates gyroscope calibration parameters containing bias, scale factors, cross-coupling errors and G-dependent coupling.com.irurueta.numerical.robust.RobustEstimatorMethod
Returns method being used for robust estimation.double
Returns threshold to be used to keep the algorithm iterating in case that best estimated threshold using median of residuals is not small enough.boolean
Indicates whether this calibrator requires quality scores for each measurement/sequence or not.void
setStopThreshold
(double stopThreshold) Sets threshold to be used to keep the algorithm iterating in case that best estimated threshold using median of residuals is not small enough.Methods inherited from class com.irurueta.navigation.inertial.calibration.gyroscope.RobustEasyGyroscopeCalibrator
attemptRefine, computeError, computePreliminarySolutions, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, create, getAccelerometerBias, getAccelerometerBias, getAccelerometerBiasAsMatrix, getAccelerometerBiasAsMatrix, getAccelerometerBiasX, getAccelerometerBiasXAsAcceleration, getAccelerometerBiasXAsAcceleration, getAccelerometerBiasY, getAccelerometerBiasYAsAcceleration, getAccelerometerBiasYAsAcceleration, getAccelerometerBiasZ, getAccelerometerBiasZAsAcceleration, getAccelerometerBiasZAsAcceleration, getAccelerometerMa, getAccelerometerMa, getAccelerometerMxy, getAccelerometerMxz, getAccelerometerMyx, getAccelerometerMyz, getAccelerometerMzx, getAccelerometerMzy, getAccelerometerSx, getAccelerometerSy, getAccelerometerSz, getConfidence, getEstimatedBiasAngularSpeedX, getEstimatedBiasAngularSpeedX, getEstimatedBiasAngularSpeedY, getEstimatedBiasAngularSpeedY, getEstimatedBiasAngularSpeedZ, getEstimatedBiasAngularSpeedZ, getEstimatedBiasAsTriad, getEstimatedBiasAsTriad, getEstimatedBiases, getEstimatedBiases, getEstimatedBiasesAsMatrix, getEstimatedBiasesAsMatrix, getEstimatedBiasStandardDeviation, getEstimatedBiasStandardDeviation, getEstimatedBiasStandardDeviationAverage, getEstimatedBiasStandardDeviationAverageAsAngularSpeed, getEstimatedBiasStandardDeviationAverageAsAngularSpeed, getEstimatedBiasStandardDeviationNorm, getEstimatedBiasStandardDeviationNormAsAngularSpeed, getEstimatedBiasStandardDeviationNormAsAngularSpeed, getEstimatedBiasX, getEstimatedBiasXStandardDeviation, getEstimatedBiasXStandardDeviationAsAngularSpeed, getEstimatedBiasXStandardDeviationAsAngularSpeed, getEstimatedBiasXVariance, getEstimatedBiasY, getEstimatedBiasYStandardDeviation, getEstimatedBiasYStandardDeviationAsAngularSpeed, getEstimatedBiasYStandardDeviationAsAngularSpeed, getEstimatedBiasYVariance, getEstimatedBiasZ, getEstimatedBiasZStandardDeviation, getEstimatedBiasZStandardDeviationAsAngularSpeed, getEstimatedBiasZStandardDeviationAsAngularSpeed, getEstimatedBiasZVariance, getEstimatedChiSq, getEstimatedCovariance, getEstimatedGg, getEstimatedMg, getEstimatedMse, getEstimatedMxy, getEstimatedMxz, getEstimatedMyx, getEstimatedMyz, getEstimatedMzx, getEstimatedMzy, getEstimatedSx, getEstimatedSy, getEstimatedSz, getInitialBias, getInitialBias, getInitialBiasAngularSpeedX, getInitialBiasAngularSpeedX, getInitialBiasAngularSpeedY, getInitialBiasAngularSpeedY, getInitialBiasAngularSpeedZ, getInitialBiasAngularSpeedZ, getInitialBiasAsMatrix, getInitialBiasAsMatrix, getInitialBiasAsTriad, getInitialBiasAsTriad, getInitialBiasX, getInitialBiasY, getInitialBiasZ, getInitialGg, getInitialGg, getInitialMg, getInitialMg, getInitialMxy, getInitialMxz, getInitialMyx, getInitialMyz, getInitialMzx, getInitialMzy, getInitialSx, getInitialSy, getInitialSz, getInliersData, getListener, getMaxIterations, getMeasurementOrSequenceType, getMinimumRequiredMeasurementsOrSequences, getPreliminarySubsetSize, getProgressDelta, getQualityScores, getSequences, isCommonAxisUsed, isCovarianceKept, isGDependentCrossBiasesEstimated, isOrderedMeasurementsOrSequencesRequired, isReady, isResultRefined, isRunning, setAccelerometerBias, setAccelerometerBias, setAccelerometerBias, setAccelerometerBias, setAccelerometerBiasX, setAccelerometerBiasX, setAccelerometerBiasY, setAccelerometerBiasY, setAccelerometerBiasZ, setAccelerometerBiasZ, setAccelerometerCrossCouplingErrors, setAccelerometerMa, setAccelerometerMxy, setAccelerometerMxz, setAccelerometerMyx, setAccelerometerMyz, setAccelerometerMzx, setAccelerometerMzy, setAccelerometerScalingFactors, setAccelerometerScalingFactorsAndCrossCouplingErrors, setAccelerometerSx, setAccelerometerSy, setAccelerometerSz, setCommonAxisUsed, setConfidence, setCovarianceKept, setGDependentCrossBiasesEstimated, setInitialBias, setInitialBias, setInitialBias, setInitialBias, setInitialBias, setInitialBiasX, setInitialBiasX, setInitialBiasY, setInitialBiasY, setInitialBiasZ, setInitialBiasZ, setInitialCrossCouplingErrors, setInitialGg, setInitialMg, setInitialMxy, setInitialMxz, setInitialMyx, setInitialMyz, setInitialMzx, setInitialMzy, setInitialScalingFactors, setInitialScalingFactorsAndCrossCouplingErrors, setInitialSx, setInitialSy, setInitialSz, setListener, setMaxIterations, setPreliminarySubsetSize, setProgressDelta, setQualityScores, setResultRefined, setSequences, setupAccelerationFixer
-
Field Details
-
DEFAULT_STOP_THRESHOLD
public static final double DEFAULT_STOP_THRESHOLDDefault value to be used for stop threshold. Stop threshold can be used to avoid keeping the algorithm unnecessarily iterating in case that best estimated threshold using median of residuals is not small enough. Once a solution is found that generates a threshold below this value, the algorithm will stop. The stop threshold can be used to prevent the LMedS algorithm iterating too many times in cases where samples have a very similar accuracy. For instance, in cases where proportion of outliers is very small (close to 0%), and samples are very accurate (i.e. 1e-6), the algorithm would iterate for a long time trying to find the best solution when indeed there is no need to do that if a reasonable threshold has already been reached. Because of this behaviour the stop threshold can be set to a value much lower than the one typically used in RANSAC, and yet the algorithm could still produce even smaller thresholds in estimated results.- See Also:
-
MIN_STOP_THRESHOLD
public static final double MIN_STOP_THRESHOLDMinimum allowed stop threshold value.- See Also:
-
stopThreshold
private double stopThresholdThreshold to be used to keep the algorithm iterating in case that best estimated threshold using median of residuals is not small enough. Once a solution is found that generates a threshold below this value, the algorithm will stop. The stop threshold can be used to prevent the LMedS algorithm iterating too many times in cases where samples have a very similar accuracy. For instance, in cases where proportion of outliers is very small (close to 0%), and samples are very accurate (i.e. 1e-6), the algorithm would iterate for a long time trying to find the best solution when indeed there is no need to do that if a reasonable threshold has already been reached. Because of this behaviour the stop threshold can be set to a value much lower than the one typically used in RANSAC, and yet the algorithm could still produce even smaller thresholds in estimated results.
-
-
Constructor Details
-
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator()Constructor. -
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg) Constructor.- Parameters:
sequences
- collection of sequences containing timestamped body kinematics measurements.initialBias
- initial gyroscope bias to be used to find a solution. This must be 3x1 and is expressed in radians per second (rad/s).initialMg
- initial gyroscope scale factors and cross coupling errors matrix. Must be 3x3.initialGg
- initial gyroscope G-dependent cross biases introduced on the gyroscope by the specific forces sensed by the accelerometer. Must be 3x3.- Throws:
IllegalArgumentException
- if any of the provided values does not have proper size.
-
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg) Constructor.- Parameters:
sequences
- collection of sequences containing timestamped body kinematics measurements.initialBias
- initial gyroscope bias to be used to find a solution. This must have length 3 and is expressed in radians per second (rad/s).initialMg
- initial gyroscope scale factors and cross coupling errors matrix. Must be 3x3.initialGg
- initial gyroscope G-dependent cross biases introduced on the gyroscope by the specific forces sensed by the accelerometer. Must be 3x3.- Throws:
IllegalArgumentException
- if any of the provided values does not have proper size.
-
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, double[] accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa) Constructor.- Parameters:
sequences
- collection of sequences containing timestamped body kinematics measurements.initialBias
- initial gyroscope bias to be used to find a solution. This must have length 3 and is expressed in radians per second (rad/s).initialMg
- initial gyroscope scale factors and cross coupling errors matrix. Must be 3x3.initialGg
- initial gyroscope G-dependent cross biases introduced on the gyroscope by the specific forces sensed by the accelerometer. Must be 3x3.accelerometerBias
- known accelerometer bias. This must have length 3 and is expressed in meters per squared second (m/s^2).accelerometerMa
- known accelerometer scale factors and cross coupling matrix. Must be 3x3.- Throws:
IllegalArgumentException
- if any of the provided values does not have proper size.
-
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, com.irurueta.algebra.Matrix accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa) Constructor.- Parameters:
sequences
- collection of sequences containing timestamped body kinematics measurements.initialBias
- initial gyroscope bias to be used to find a solution. This must be 3x1 and is expressed in radians per second (rad/s).initialMg
- initial gyroscope scale factors and cross coupling errors matrix. Must be 3x3.initialGg
- initial gyroscope G-dependent cross biases introduced on the gyroscope by the specific forces sensed by the accelerometer. Must be 3x3.accelerometerBias
- known accelerometer bias. This must be 3x1 and is expressed in meters per squared second (m/s^2).accelerometerMa
- known accelerometer scale factors and cross coupling matrix. Must be 3x3.- Throws:
IllegalArgumentException
- if any of the provided values does not have proper size.
-
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg) Constructor.- Parameters:
sequences
- collection of sequences containing timestamped body kinematics measurements.commonAxisUsed
- indicates whether z-axis is assumed to be common for accelerometer and gyroscope.estimateGDependentCrossBiases
- true if G-dependent cross biases will be estimated, false otherwise.initialBias
- initial gyroscope bias to be used to find a solution. This must be 3x1 and is expressed in radians per second (rad/s).initialMg
- initial gyroscope scale factors and cross coupling errors matrix. Must be 3x3.initialGg
- initial gyroscope G-dependent cross biases introduced on the gyroscope by the specific forces sensed by the accelerometer. Must be 3x3.- Throws:
IllegalArgumentException
- if any of the provided values does not have proper size.
-
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg) Constructor.- Parameters:
sequences
- collection of sequences containing timestamped body kinematics measurements.commonAxisUsed
- indicates whether z-axis is assumed to be common for accelerometer and gyroscope.estimateGDependentCrossBiases
- true if G-dependent cross biases will be estimated, false otherwise.initialBias
- initial gyroscope bias to be used to find a solution. This must have length 3 and is expressed in radians per second (rad/s).initialMg
- initial gyroscope scale factors and cross coupling errors matrix. Must be 3x3.initialGg
- initial gyroscope G-dependent cross biases introduced on the gyroscope by the specific forces sensed by the accelerometer. Must be 3x3.- Throws:
IllegalArgumentException
- if any of the provided values does not have proper size.
-
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, double[] initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, double[] accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa) Constructor.- Parameters:
sequences
- collection of sequences containing timestamped body kinematics measurements.commonAxisUsed
- indicates whether z-axis is assumed to be common for accelerometer and gyroscope.estimateGDependentCrossBiases
- true if G-dependent cross biases will be estimated, false otherwise.initialBias
- initial gyroscope bias to be used to find a solution. This must have length 3 and is expressed in radians per second (rad/s).initialMg
- initial gyroscope scale factors and cross coupling errors matrix. Must be 3x3.initialGg
- initial gyroscope G-dependent cross biases introduced on the gyroscope by the specific forces sensed by the accelerometer. Must be 3x3.accelerometerBias
- known accelerometer bias. This must have length 3 and is expressed in meters per squared second (m/s^2).accelerometerMa
- known accelerometer scale factors and cross coupling matrix. Must be 3x3.- Throws:
IllegalArgumentException
- if any of the provided values does not have proper size.
-
LMedSRobustEasyGyroscopeCalibrator
public LMedSRobustEasyGyroscopeCalibrator(List<BodyKinematicsSequence<StandardDeviationTimedBodyKinematics>> sequences, boolean commonAxisUsed, boolean estimateGDependentCrossBiases, com.irurueta.algebra.Matrix initialBias, com.irurueta.algebra.Matrix initialMg, com.irurueta.algebra.Matrix initialGg, com.irurueta.algebra.Matrix accelerometerBias, com.irurueta.algebra.Matrix accelerometerMa) Constructor.- Parameters:
sequences
- collection of sequences containing timestamped body kinematics measurements.commonAxisUsed
- indicates whether z-axis is assumed to be common for accelerometer and gyroscope.estimateGDependentCrossBiases
- true if G-dependent cross biases will be estimated, false otherwise.initialBias
- initial gyroscope bias to be used to find a solution. This must be 3x1 and is expressed in radians per second (rad/s).initialMg
- initial gyroscope scale factors and cross coupling errors matrix. Must be 3x3.initialGg
- initial gyroscope G-dependent cross biases introduced on the gyroscope by the specific forces sensed by the accelerometer. Must be 3x3.accelerometerBias
- known accelerometer bias. This must have length 3 and is expressed in meters per squared second (m/s^2).accelerometerMa
- known accelerometer scale factors and cross coupling matrix. Must be 3x3.- Throws:
IllegalArgumentException
- if any of the provided values does not have proper size.
-
-
Method Details
-
getStopThreshold
public double getStopThreshold()Returns threshold to be used to keep the algorithm iterating in case that best estimated threshold using median of residuals is not small enough. Once a solution is found that generates a threshold below this value, the algorithm will stop. The stop threshold can be used to prevent the LMedS algorithm to iterate too many times in cases where samples have a very similar accuracy. For instance, in cases where proportion of outliers is very small (close to 0%), and samples are very accurate (i.e. 1e-6), the algorithm would iterate for a long time trying to find the best solution when indeed there is no need to do that if a reasonable threshold has already been reached. Because of this behaviour the stop threshold can be set to a value much lower than the one typically used in RANSAC, and yet the algorithm could still produce even smaller thresholds in estimated results.- Returns:
- stop threshold to stop the algorithm prematurely when a certain accuracy has been reached.
-
setStopThreshold
public void setStopThreshold(double stopThreshold) throws com.irurueta.navigation.LockedException Sets threshold to be used to keep the algorithm iterating in case that best estimated threshold using median of residuals is not small enough. Once a solution is found that generates a threshold below this value, the algorithm will stop. The stop threshold can be used to prevent the LMedS algorithm to iterate too many times in cases where samples have a very similar accuracy. For instance, in cases where proportion of outliers is very small (close to 0%), and samples are very accurate (i.e. 1e-6), the algorithm would iterate for a long time trying to find the best solution when indeed there is no need to do that if a reasonable threshold has already been reached. Because of this behaviour the stop threshold can be set to a value much lower than the one typically used in RANSAC, and yet the algorithm could still produce even smaller thresholds in estimated results.- Parameters:
stopThreshold
- stop threshold to stop the algorithm prematurely when a certain accuracy has been reached.- Throws:
IllegalArgumentException
- if provided value is zero or negative.com.irurueta.navigation.LockedException
- if calibrator is currently running.
-
calibrate
public void calibrate() throws com.irurueta.navigation.LockedException, com.irurueta.navigation.NotReadyException, CalibrationExceptionEstimates gyroscope calibration parameters containing bias, scale factors, cross-coupling errors and G-dependent coupling.- Throws:
com.irurueta.navigation.LockedException
- if calibrator is currently running.com.irurueta.navigation.NotReadyException
- if calibrator is not ready.CalibrationException
- if estimation fails for numerical reasons.
-
getMethod
public com.irurueta.numerical.robust.RobustEstimatorMethod getMethod()Returns method being used for robust estimation.- Specified by:
getMethod
in classRobustEasyGyroscopeCalibrator
- Returns:
- method being used for robust estimation.
-
isQualityScoresRequired
public boolean isQualityScoresRequired()Indicates whether this calibrator requires quality scores for each measurement/sequence or not.- Returns:
- true if quality scores are required, false otherwise.
-