Petrus Hyvönen <petrus.hyvonen@gmail.com> a écrit :
Hi,
Hi Petrus,
I'm abit confused by the computeLimitVisibilityPoint results. For a
topocentric frame I would assume azimuth 0.0 is north direction? This code
snippet (python) translate it to a longitude offset to the station?
Any comments welcome.
Regards
/Petrus
ITRF2005 = FramesFactory.getITRF2005()
earth = OneAxisEllipsoid(Constants.WGS84_EARTH_EQUATORIAL_RADIUS,
...: Constants.WGS84_EARTH_FLATTENING,
...: ITRF2005)
...: staFrame = TopocentricFrame(earth, GeodeticPoint(0.0,0.0,0.0),
'test')
print
staFrame.computeLimitVisibilityPoint(Constants.WGS84_EARTH_EQUATORIAL_RADIUS+600000,
0.0,5.0)
{lat: 0 deg, lon: 148,4942265118 deg, alt: 599?999,9999999991}
staFrame.getNorth()
Out[65]: <Vector3D: {-0; -0; 1}>
I think you have found a bug. A very quick look makes me think there is an inconsistency between the pointAtDistance method and the other ones like getAzimuth/getElevation) about how the frame is oriented.
Could you file a bug report for this (and better yet analyze it and propose some fix) ?
best regards, Luc
-- _____________________________________________ Petrus Hyvönen, Uppsala, Sweden Mobile Phone/SMS:+46 73 803 19 00
---------------------------------------------------------------- This message was sent using IMP, the Internet Messaging Program.