On 2013-03-28, at 7:29 AM, KURT PETERS wrote:
> If I'm using this SGP4 library: https://pypi.python.org/pypi/sgp4/
> which provides x,y,z of a satellite with respect to the center of
> the Earth, what do you think would be the best way to calculate the
> distance from the satellite to a lat/long point on the Earth using
> MATPLOT library basemaps?
>
> I suppose using WGS 84, I could do something like:
> # setup of basemap ('lcc' = lambert conformal conic).
> # use major and minor sphere radii from WGS84 ellipsoid.
> m =
> Basemap
> (llcrnrlon=-145.5,llcrnrlat=1.,urcrnrlon=-2.566,urcrnrlat=46.352,\
> rsphere=(6378137.00,6356752.3142),\
> resolution='l',area_thresh=1000.,projection='lcc',\
> lat_1=50.,lon_0=-107.,ax=ax)
>
> x,y = m(lon, lat)
>
> But how could I get "Z" for a particular coordinate so I could be
> this norm calculation? And... would the x,y yielded above be in the
> same coordinate frame (ie, from the center of the Earth) as the SGP4
> library in the first place?
>
> Regards,
> Kurt
The x,y you get from Basemap appear to be tailored to the particular
projection
you are using. If you change the projection, then you will probably
get different
x,y values for the same (lat,lon). This is very useful for mapping,
but it is not what you need.
You need convert your position P from geodetic coordinates (lat, lon)
(more generally, (lat, lon, elevation)) to geocentric coordinates
(x,y,z). Geocentric coordinates are also called Earth-centered-Earth-
fixed (ECEF or ECF) coordinates.
Here is a link to the formulas you need. Notice that the author uses
the terminology Latitude,Longitude,Altitude (LLA) to refer to geodetic
coordinates.
http://www.satsleuth.com/GPS_ECEF_Datum_transformation.htm
This can be difficult to interpret if you are doing it for the first
time (and even if you have done it many times before!!). Feel free to
email me directly if you have any questions. I would be glad to help.
-John
--
Own the Future-Intel® Level Up Game Demo Contest 2013
Rise to greatness in Intel's independent game demo contest.
Compete for recognition, cash, and the chance to get your game
on Steam. $5K grand prize plus 10 genre and skill prizes.
Submit your demo by 6/6/13. http://p.sf.net/sfu/intel_levelupd2d
___
Matplotlib-users mailing list
Matplotlib-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/matplotlib-users
