First of all,
Let me thank you for this package. It's been solving lot's of problems for me.

I have one question though. Why is ge_quiver not being distributed anymore? I had to recover it from a much older version. It works, so I'm not sure why it's not being offered with the rest of the package.

Maybe it's speed of the ge_quiver algorithm. To that, I have made a few tweaks to make it faster.. I would gladly send you the code if you want.

Darin,
I have been importing some ECI data from the OrbitTools C++ library and have stumbled upon some weird behavior.

In MATLAB I am using your functions to transform ECI to ECEF and then using the MATLAB function, ecef2geodetic to transform into normal azimuth/elevation/height coordinate space.

The weird behavior is that despite having smooth ECI data (smooth as in no jumps in the data) as inputs, I am getting discontinuous ECEF data (graph for r_ECEF(1,:) and r_ECEF(2,:) looks like saw). This is happening for GOES-8 and GOES-15 satellite data.

What is weirder is that on a macro scale (ignoring variations, including the sudden discontinuities), the position of the satellite in the mapa mundi seems to be correct.

Chad,
The two code lines you are referencing are trying to find the midnight that occurs before the provided Julian Date.
As written, the first check will be true if the Julian date occurs after midnight and 12+ hours before the time.
As written, the second check will be true if the Julian date occurs after midnight and less than 12 hours before the time.
If both checks are true, the code will use the last midnight that occurred. This is the desired functionality.
Because I am rounding down on the Julian date and subtracting a half day, there is no possible way that jd could be less than jdMin.
-Darin

Comment only

29 Aug 2013

Convert ECI to ECEF Coordinates
Transform ECI position, velocity and acceleration to ECEF position, velocity & acceleration.

Darin,
In your JD2GMST function, I have a question related to the following two lines:
jd0(jd>jdMin) = jdMin(jd>jdMin);
jd0(jd>jdMax) = jdMax(jd>jdMax);
This appears to be limiting the Julian day between a lower bound and an upper bound. Don't you want the first line to be:
jd0(jd<jdMin) = jdMin(jd<jdMin);
Otherwise, it would appear that no actual lower bounding occurs...
Chad

Comment only

27 Mar 2013

Convert ECI to ECEF Coordinates
Transform ECI position, velocity and acceleration to ECEF position, velocity & acceleration.

Darin,
I have been importing some ECI data from the OrbitTools C++ library and have stumbled upon some weird behavior.
In MATLAB I am using your functions to transform ECI to ECEF and then using the MATLAB function, ecef2geodetic to transform into normal azimuth/elevation/height coordinate space.
The weird behavior is that despite having smooth ECI data (smooth as in no jumps in the data) as inputs, I am getting discontinuous ECEF data (graph for r_ECEF(1,:) and r_ECEF(2,:) looks like saw). This is happening for GOES-8 and GOES-15 satellite data.
What is weirder is that on a macro scale (ignoring variations, including the sudden discontinuities), the position of the satellite in the mapa mundi seems to be correct.

5

28 Jun 2012

Convert ECI to ECEF Coordinates
Transform ECI position, velocity and acceleration to ECEF position, velocity & acceleration.

Drake,
J2000 is factored into the ECI to ECEF conversion routine through the call to JD2GAST.m (converts Julian Date to Greenwhich Apparent Sidereal Time). In order to do this, the effects of nutation are considered, and depend on the correct number of Julian centuries since J2000.

Comment only