Hi @nrol,
I think it’s a bug in AngularRaDec. Line 171 should probably use transformVector
instead of transformPosition
(the first one apply only the rotation while the second applies rotation + translation).
Could you please open a bug for this one too ? (two bugs in a post !! Thanks)
Two solutions in my opinion:
-
You already have the reference measurements as Orekit
ObservedMeasurements
→ Use aResidualHandler
with your propagator (see this post for an example)
Afterwards you get a Map of (Observed, Estimated) and you can compare them. -
You don’t have a list or it’s too difficult to obtain.
→ Use the list of dates in the handler instead of a list ofObservedMeasurement
In theResidualHandler
before line 90 you need to initialize theobserved
variable.
Initialize it like you would do a measurement (with sat, frame, weights, sigmas etc.) and put [0,0] for the observed values.
When observed.estimate(…) is called, Orekit will compute the estimated (i.e. theoretical) measurement. However the residuals in this case will be meaningless.