An official website of the United States government
This Jupyter notebook demonstrates using the sscws Python package to access satellite location and (modeled) magnetic field information. This notebook contains the following sections:
Install the prerequisite software from the Python Package Index if it is not already installed.
) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.#%pip install sscws
Execute some preliminary code that is necessary before the code that follows.
import matplotlib as mpl
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt
from packaging import version
from sscws.sscws import SscWs
ssc = SscWs()
The following code demonstrates how to get the list of available observatories.
result = ssc.get_observatories()
for observatory in result['Observatory'][:5]:
print('{:15s} {:20.20s} {:25s}'.format(observatory['Id'], observatory['Name'], observatory['StartTime'].isoformat()))
print('...')
ace ACE 1997-08-25T17:48:00+00:00 active Active 1989-09-29T00:00:00+00:00 aec AE-C 1973-12-17T08:01:00+00:00 aed AE-D 1975-10-17T00:00:00+00:00 aee AE-E 1975-11-20T21:04:00+00:00 ...
The following code demonstrates how to get the list of available ground stations.
result = ssc.get_ground_stations()
for ground_station in result['GroundStation'][:5]:
location = ground_station['Location']
print('{:5s} {:20.20s} {:7.2f} {:7.2f}'.format(ground_station['Id'], ground_station['Name'],location['Latitude'], location['Longitude']))
print('...')
SPA South Pole -89.99 -13.32 MCM McMurdo -77.85 166.70 SBF Scott Base -77.85 166.75 SPL Siple -76.00 -84.00 HBA Halley Bay -75.52 -26.60 ...
The following code gets location information for the International Space Station (ISS) spacecraft and prints a few X coordinate values.
result = ssc.get_locations(['iss'],
['2020-01-01T00:00:00Z',
'2020-01-01T01:00:00Z'])
data = result['Data'][0]
coords = data['Coordinates'][0]
print(coords['X'][:5])
[-1473.76002814 -1035.57214121 -592.61584877 -146.91796194 299.48049963]
And the following code plots the ISS location information.
fig = plt.figure()
if version.parse(mpl.__version__) < version.parse('3.4'):
ax = fig.gca(projection='3d')
else:
ax = Axes3D(fig, auto_add_to_figure=False)
fig.add_axes(ax)
ax.set_xlabel('X (km)')
ax.set_ylabel('Y (km)')
ax.set_zlabel('Z (km)')
title = data['Id'] + ' Orbit (' + \
coords['CoordinateSystem'].value.upper() + ')'
ax.plot(coords['X'], coords['Y'], coords['Z'], label=title)
ax.legend()
plt.show()
View the sscws API for additonal features. Additional notebook examples are also available.
