satkit 0.8.5__tar.gz → 0.9.0__tar.gz

{satkit-0.8.5 → satkit-0.9.0}/CHANGES.md

@@ -1,9 +1,10 @@
 # Changelog
 
+
 ## 0.5.5 - 2025-01-27
 
 ### Low-Precision Ephemeris
- - Coefficients for low-precision planetary ephemerides did not match paper or JPL website referenced in documentation.  Not sure where original numbers came from.  Some plantes (e.g., Mercury) matched.  Others (e.g., Mars) did not, although numbers for Mars did match a google search. Very strange ... regardless, update so they match.  Coefficients for years 1800 to 2050 AD are correct and remain unchanged.
+ - Coefficients for low-precision planetary ephemerides did not match paper or JPL website referenced in documentation.  Not sure where original nu mbers came from.  Some plantes (e.g., Mercury) matched.  Others (e.g., Mars) did not, although numbers for Mars did match a google search. Very strange ... regardless, update so they match.  Coefficients for years 1800 to 2050 AD are correct and remain unchanged.
  - Add Pluto as a planet for low-precision ephemerides (it is included in reference paper, but not JPL website)
 
 ### Two-Line Element Sets
@@ -143,3 +144,14 @@
   - Add "x", "y", "z", "w" property getters to quaternion in python
   - Add pickle serialize/deserialize tests in python testing
   - allow for duration division by float and by other duration in python
+
+-----------------
+
+
+## Rust 0.9.0, Python 0.9.0
+  - Support Orbital Mean-Element Messages (OMM) in JSON format
+  - Add OMM documentation, tests, and python example
+  - Structured output of SGP4 propagator in rust
+  - add "as_datetime" function in python for satkit.time (for consistent nomenclature)
+  - Rename time-interval boundary nomenclature from `start/stop` to `begin/end` across Rust + Python APIs (including propagation functions and related settings/results).
+

{satkit-0.8.5 → satkit-0.9.0}/Cargo.toml

@@ -1,6 +1,6 @@
 [package]
 name = "satkit"
-version = "0.8.3"
+version = "0.9.0"
 edition = "2021"
 description = "Satellite Toolkit"
 readme = "README.md"
@@ -12,7 +12,7 @@ keywords = ["satellite", "orbit", "ephemeris", "tle", "astrodynamics"]
 categories = ["aerospace", "algorithms", "mathematics", "science"]
 
 
-# See more keys and their definitions at 
+# See more keys and their definitions at
 # https://doc.rust-lang.org/cargo/reference/manifest.html
 [lib]
 name = "satkit"
@@ -34,6 +34,7 @@ ureq = "3.1.2"
 json = "0.12.4"
 process_path = "0.1.4"
 serde = { version = "1.0", features = ["derive"] }
+serde_json = "1.0"
 serde-pickle = "1.2.0"
 itertools = "0.14.0"
 anyhow = "1"

{satkit-0.8.5 → satkit-0.9.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: satkit
-Version: 0.8.5
+Version: 0.9.0
 Summary: Satellite Orbital Dynamics Toolkit
 Author-email: Steven Michael <ssmichael@gmail.com>
 Maintainer-email: Steven Michael <ssmichael@gmail.com>
@@ -40,6 +40,9 @@ Description-Content-Type: text/markdown
 License-File: LICENSE
 Requires-Dist: numpy>=1.0.0
 Requires-Dist: satkit-data>=0.7.0
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Requires-Dist: xmltodict; extra == "test"
 Dynamic: license-file
 
 # Satellite Toolkit (satkit)
@@ -48,7 +51,7 @@ Dynamic: license-file
 
 Satkit provides robust, high-performance satellite orbital mechanics calculations with a clean, intuitive API. Built from the ground up in Rust for maximum performance and memory safety, it offers complete Python bindings for all functionality, making advanced orbital mechanics accessible to both systems programmers and data scientists.
 
------ 
+-----
 
 ![Build Passing?](https://github.com/ssmichael1/satkit/actions/workflows/build.yml/badge.svg)
 ![Wheel Passing?](https://github.com/ssmichael1/satkit/actions/workflows/wheels.yml/badge.svg)
@@ -62,7 +65,7 @@ Satkit provides robust, high-performance satellite orbital mechanics calculation
 ![PyPI - Status](https://img.shields.io/pypi/status/satkit)
 ![PyPI - Downloads](https://img.shields.io/pypi/dm/satkit)
 ![Read the Docs](https://img.shields.io/readthedocs/satellite-toolkit)
- 
+
 ------
 
 ## Language Bindings
@@ -84,7 +87,7 @@ High-precision conversions between multiple reference frames with full support f
   - **CIRF** - Celestial Intermediate Reference Frame (IAU-2006 intermediate)
   - **TIRF** - Terrestrial Intermediate Reference Frame (Earth-rotation intermediate)
   - **Geodetic** - Latitude, longitude, altitude with WGS-84 ellipsoid
-  
+
 ### Orbit Propagation
 Multiple propagation methods optimized for different accuracy and performance requirements:
   - **Numerical Integration**: High-precision propagation using adaptive Runge-Kutta 9(8) methods with dense output
@@ -95,27 +98,29 @@ Multiple propagation methods optimized for different accuracy and performance re
     - Full AFSPC and improved mode support
     - TLE fitting from high-precision states with drag estimation
     - Batch processing for multiple satellites
+    - Support TLE and Orbital Mean-Element Messages
   - **Keplerian**: Fast analytical two-body propagation for preliminary analysis
-  
-### Force Models
+
+### Numerical Integration Force Models
 Comprehensive perturbation modeling for high-fidelity orbit propagation:
   - **Earth Gravity**: Spherical harmonic models up to degree/order 360
     - Multiple models: JGM2, JGM3, EGM96, ITU GRACE16
     - Efficient computation with configurable truncation order
     - Gravity gradient support for state transition matrix
   - **Third-Body Gravity**: Solar and lunar perturbations using JPL ephemerides
-  - **Atmospheric Drag**: NRLMSISE-00 density model with space weather integration
+  - **Atmospheric Drag**:
+    - NRLMSISE-00 density model with space weather integration
     - Automatic space weather data updates (F10.7, Ap index)
     - Configurable ballistic coefficients
   - **Solar Radiation Pressure**: Cannon-ball model with shadow function
-  
-### Ephemerides
+
+### Solar System Ephemerides
 Access to high-precision solar system body positions:
   - **JPL DE440/DE441**: State-of-the-art planetary ephemerides
     - Chebyshev polynomial interpolation for accuracy
     - Support for all major planets, sun, moon, and solar system barycenter
   - **Low-Precision Models**: Fast analytical models for sun and moon when high precision isn't required
-  
+
 ### Time Systems
 Comprehensive support for all standard astronomical time scales:
   - **UTC** - Coordinated Universal Time with leap second handling
@@ -125,7 +130,7 @@ Comprehensive support for all standard astronomical time scales:
   - **UT1** - Universal Time with Earth orientation corrections
   - **GPS** - GPS Time
   - Automatic conversion between all time scales with microsecond precision
-  
+
 ### Geodetic Utilities
 - **Geodesic Calculations**: Accurate distance and azimuth between ground locations using Vincenty's formulae
 - **Coordinate Conversions**: ITRF ↔ Geodetic ↔ East-North-Up ↔ North-East-Down
@@ -178,12 +183,12 @@ The library requires external data files for various calculations. These are aut
   - DE440/DE441 planetary ephemerides (~100 MB)
   - Provides positions of sun, moon, planets, and solar system barycenter
   - Valid for years 1550-2650 CE
-  
+
 - **Gravity Models** ([ICGEM](http://icgem.gfz-potsdam.de/home))
   - JGM2, JGM3, EGM96, ITU GRACE16 spherical harmonic coefficients
   - International Centre for Global Earth Models standardized format
   - Up to degree/order 360 for high-fidelity propagation
-  
+
 - **IERS Nutation Tables** ([IERS Conventions](https://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn36.html))
   - IAU-2006 nutation series coefficients
   - Required for GCRF ↔ ITRF transformations
@@ -196,7 +201,7 @@ The library requires external data files for various calculations. These are aut
   - Ap geomagnetic index
   - Critical for atmospheric density modeling and drag calculations
   - Updated daily by NOAA Space Weather Prediction Center
-  
+
 - **Earth Orientation Parameters** ([Celestrak](https://celestrak.org/SpaceData/))
   - Polar motion (x, y)
   - UT1-UTC time difference
@@ -218,15 +223,15 @@ The library includes comprehensive test suites ensuring correctness of calculati
 - **JPL Ephemerides**: Validated against JPL-provided test vectors for Chebyshev polynomial interpolation
   - Over 10,000 test cases covering all planets and time ranges
   - Accuracy verified to within JPL's published tolerances (sub-meter precision)
-  
+
 - **SGP4**: Verified using official test vectors from the original C++ distribution
   - All test cases from Vallado's SGP4 implementation
   - Includes edge cases and error conditions
-  
+
 - **Coordinate Transformations**: Cross-validated against multiple reference implementations
   - SOFA library comparisons for IAU-2006 transformations
   - Vallado test cases for GCRF ↔ ITRF conversions
-  
+
 - **Numerical Propagation**: Validated against high-precision commercial tools
   - Orbit fits to GPS SP3 ephemerides
   - Multi-day propagations with sub-meter accuracy
@@ -266,7 +271,7 @@ pip install satkit
 
 Pre-built binary wheels are available for:
 - **Windows**: AMD64
-- **macOS**: Intel (x86_64) and Apple Silicon (ARM64)  
+- **macOS**: Intel (x86_64) and Apple Silicon (ARM64)
 - **Linux**: x86_64 and ARM64 (aarch64)
 - **Python versions**: 3.8, 3.9, 3.10, 3.11, 3.12, 3.13, 3.14
 
@@ -336,9 +341,9 @@ settings.gravity_order = 8
 
 # Propagate for 1 day
 result = sk.propagate(
-    state0, 
-    time0, 
-    stop=time0 + sk.duration.from_days(1),
+    state0,
+    time0,
+    end=time0 + sk.duration.from_days(1),
     propsettings=settings
 )
 
@@ -375,13 +380,13 @@ use satkit::{Instant, SolarSystem, jplephem};
 fn main() -> Result<(), Box<dyn std::error::Error>> {
     // Create time instant
     let time = Instant::from_datetime(2024, 1, 1, 0, 0, 0.0)?;
-    
+
     // Get Moon position and velocity in GCRF
     let (pos, vel) = jplephem::geocentric_state(SolarSystem::Moon, &time)?;
-    
+
     println!("Moon position: {:?}", pos);
     println!("Moon velocity: {:?}", vel);
-    
+
     Ok(())
 }
 ```

{satkit-0.8.5 → satkit-0.9.0}/README.md

@@ -4,7 +4,7 @@
 
 Satkit provides robust, high-performance satellite orbital mechanics calculations with a clean, intuitive API. Built from the ground up in Rust for maximum performance and memory safety, it offers complete Python bindings for all functionality, making advanced orbital mechanics accessible to both systems programmers and data scientists.
 
------ 
+-----
 
 ![Build Passing?](https://github.com/ssmichael1/satkit/actions/workflows/build.yml/badge.svg)
 ![Wheel Passing?](https://github.com/ssmichael1/satkit/actions/workflows/wheels.yml/badge.svg)
@@ -18,7 +18,7 @@ Satkit provides robust, high-performance satellite orbital mechanics calculation
 ![PyPI - Status](https://img.shields.io/pypi/status/satkit)
 ![PyPI - Downloads](https://img.shields.io/pypi/dm/satkit)
 ![Read the Docs](https://img.shields.io/readthedocs/satellite-toolkit)
- 
+
 ------
 
 ## Language Bindings
@@ -40,7 +40,7 @@ High-precision conversions between multiple reference frames with full support f
   - **CIRF** - Celestial Intermediate Reference Frame (IAU-2006 intermediate)
   - **TIRF** - Terrestrial Intermediate Reference Frame (Earth-rotation intermediate)
   - **Geodetic** - Latitude, longitude, altitude with WGS-84 ellipsoid
-  
+
 ### Orbit Propagation
 Multiple propagation methods optimized for different accuracy and performance requirements:
   - **Numerical Integration**: High-precision propagation using adaptive Runge-Kutta 9(8) methods with dense output
@@ -51,27 +51,29 @@ Multiple propagation methods optimized for different accuracy and performance re
     - Full AFSPC and improved mode support
     - TLE fitting from high-precision states with drag estimation
     - Batch processing for multiple satellites
+    - Support TLE and Orbital Mean-Element Messages
   - **Keplerian**: Fast analytical two-body propagation for preliminary analysis
-  
-### Force Models
+
+### Numerical Integration Force Models
 Comprehensive perturbation modeling for high-fidelity orbit propagation:
   - **Earth Gravity**: Spherical harmonic models up to degree/order 360
     - Multiple models: JGM2, JGM3, EGM96, ITU GRACE16
     - Efficient computation with configurable truncation order
     - Gravity gradient support for state transition matrix
   - **Third-Body Gravity**: Solar and lunar perturbations using JPL ephemerides
-  - **Atmospheric Drag**: NRLMSISE-00 density model with space weather integration
+  - **Atmospheric Drag**:
+    - NRLMSISE-00 density model with space weather integration
     - Automatic space weather data updates (F10.7, Ap index)
     - Configurable ballistic coefficients
   - **Solar Radiation Pressure**: Cannon-ball model with shadow function
-  
-### Ephemerides
+
+### Solar System Ephemerides
 Access to high-precision solar system body positions:
   - **JPL DE440/DE441**: State-of-the-art planetary ephemerides
     - Chebyshev polynomial interpolation for accuracy
     - Support for all major planets, sun, moon, and solar system barycenter
   - **Low-Precision Models**: Fast analytical models for sun and moon when high precision isn't required
-  
+
 ### Time Systems
 Comprehensive support for all standard astronomical time scales:
   - **UTC** - Coordinated Universal Time with leap second handling
@@ -81,7 +83,7 @@ Comprehensive support for all standard astronomical time scales:
   - **UT1** - Universal Time with Earth orientation corrections
   - **GPS** - GPS Time
   - Automatic conversion between all time scales with microsecond precision
-  
+
 ### Geodetic Utilities
 - **Geodesic Calculations**: Accurate distance and azimuth between ground locations using Vincenty's formulae
 - **Coordinate Conversions**: ITRF ↔ Geodetic ↔ East-North-Up ↔ North-East-Down
@@ -134,12 +136,12 @@ The library requires external data files for various calculations. These are aut
   - DE440/DE441 planetary ephemerides (~100 MB)
   - Provides positions of sun, moon, planets, and solar system barycenter
   - Valid for years 1550-2650 CE
-  
+
 - **Gravity Models** ([ICGEM](http://icgem.gfz-potsdam.de/home))
   - JGM2, JGM3, EGM96, ITU GRACE16 spherical harmonic coefficients
   - International Centre for Global Earth Models standardized format
   - Up to degree/order 360 for high-fidelity propagation
-  
+
 - **IERS Nutation Tables** ([IERS Conventions](https://www.iers.org/IERS/EN/Publications/TechnicalNotes/tn36.html))
   - IAU-2006 nutation series coefficients
   - Required for GCRF ↔ ITRF transformations
@@ -152,7 +154,7 @@ The library requires external data files for various calculations. These are aut
   - Ap geomagnetic index
   - Critical for atmospheric density modeling and drag calculations
   - Updated daily by NOAA Space Weather Prediction Center
-  
+
 - **Earth Orientation Parameters** ([Celestrak](https://celestrak.org/SpaceData/))
   - Polar motion (x, y)
   - UT1-UTC time difference
@@ -174,15 +176,15 @@ The library includes comprehensive test suites ensuring correctness of calculati
 - **JPL Ephemerides**: Validated against JPL-provided test vectors for Chebyshev polynomial interpolation
   - Over 10,000 test cases covering all planets and time ranges
   - Accuracy verified to within JPL's published tolerances (sub-meter precision)
-  
+
 - **SGP4**: Verified using official test vectors from the original C++ distribution
   - All test cases from Vallado's SGP4 implementation
   - Includes edge cases and error conditions
-  
+
 - **Coordinate Transformations**: Cross-validated against multiple reference implementations
   - SOFA library comparisons for IAU-2006 transformations
   - Vallado test cases for GCRF ↔ ITRF conversions
-  
+
 - **Numerical Propagation**: Validated against high-precision commercial tools
   - Orbit fits to GPS SP3 ephemerides
   - Multi-day propagations with sub-meter accuracy
@@ -222,7 +224,7 @@ pip install satkit
 
 Pre-built binary wheels are available for:
 - **Windows**: AMD64
-- **macOS**: Intel (x86_64) and Apple Silicon (ARM64)  
+- **macOS**: Intel (x86_64) and Apple Silicon (ARM64)
 - **Linux**: x86_64 and ARM64 (aarch64)
 - **Python versions**: 3.8, 3.9, 3.10, 3.11, 3.12, 3.13, 3.14
 
@@ -292,9 +294,9 @@ settings.gravity_order = 8
 
 # Propagate for 1 day
 result = sk.propagate(
-    state0, 
-    time0, 
-    stop=time0 + sk.duration.from_days(1),
+    state0,
+    time0,
+    end=time0 + sk.duration.from_days(1),
     propsettings=settings
 )
 
@@ -331,13 +333,13 @@ use satkit::{Instant, SolarSystem, jplephem};
 fn main() -> Result<(), Box<dyn std::error::Error>> {
     // Create time instant
     let time = Instant::from_datetime(2024, 1, 1, 0, 0, 0.0)?;
-    
+
     // Get Moon position and velocity in GCRF
     let (pos, vel) = jplephem::geocentric_state(SolarSystem::Moon, &time)?;
-    
+
     println!("Moon position: {:?}", pos);
     println!("Moon velocity: {:?}", vel);
-    
+
     Ok(())
 }
 ```

{satkit-0.8.5 → satkit-0.9.0}/pyproject.toml

@@ -9,7 +9,7 @@ requires-python = ">= 3.8"
 authors = [{ name = "Steven Michael", email = "ssmichael@gmail.com" }]
 maintainers = [{ name = "Steven Michael", email = "ssmichael@gmail.com" }]
 readme = "README.md"
-version = "0.8.5"
+version = "0.9.0"
 license = { file = "LICENSE" }
 description = "Satellite Orbital Dynamics Toolkit"
 keywords = [
@@ -37,6 +37,8 @@ classifiers = [
     "Programming Language :: Python :: 3.13",
     "Programming Language :: Python :: 3.14",
 ]
+[project.optional-dependencies]
+test = ["pytest", "xmltodict"]
 
 
 [tool.setuptools.packages.find]
@@ -69,7 +71,7 @@ before-build = [
 ]
 environment = { PATH = "$PATH:$HOME/.cargo/bin", MACOSX_DEPLOYMENT_TARGET = "10.12" }
 # Skip trying to test arm64 builds on Intel Macs
-test-requires = "pytest"
+test-requires = "pytest xmltodict"
 test-command = [
     "SATKIT_DATA={package}/astro-data SATKIT_TESTVEC_ROOT={package}/satkit-testvecs pytest {project}/python/test/test.py",
 ]
@@ -83,7 +85,7 @@ before-build = [
     "curl https://sh.rustup.rs -sSf | sh -s -- --default-toolchain=stable --profile=minimal -y",
     "rustup show",
 ]
-test-requires = "pytest"
+test-requires = "pytest xmltodict"
 test-command = [
     "SATKIT_DATA={package}/astro-data SATKIT_TESTVEC_ROOT={package}/satkit-testvecs pytest {project}/python/test/test.py",
 ]
@@ -98,7 +100,7 @@ before-build = [
     "rustup override set stable-x86_64-pc-windows-msvc",
     "rustup show",
 ]
-test-requires = "pytest"
+test-requires = "pytest xmltodict"
 before-test = [
     "set SATKIT_DATA={package}\\astro-data",
     "set SATKIT_TESTVEC_ROOT={package}\\satkit-testvecs",

satkit-0.9.0/python/docs/source/TLE.md

@@ -0,0 +1,118 @@
+# SGP4, Two-Line Element Sets (TLEs), and Orbital Mean-Element Messages (OMMs)
+
+## Satellite Catalog
+
+The United States maintains a public catalog of Earth-orbiting objects, including active satellites, rocket bodies, and debris.
+
+There are multiple places to access the catalog, including:
+* <https://www.celestrak.org/>
+* <https://www.space-track.org/>
+
+## SGP4
+
+**SGP4 (Simplified General Perturbations No. 4)** is an analytical orbital propagation model created in the **1960s–1970s by NORAD** to efficiently predict the motion of **Earth-orbiting satellites** from **Two-Line Element (TLE)** data.
+
+It was developed to support U.S. space surveillance as a fast, closed-form alternative to numerical integration, modeling Earth’s oblateness (J2–J4), atmospheric drag via the TLE *B\** term, and key secular and periodic perturbations.
+
+Today, SGP4 is the standard propagator for TLEs published by organizations like NORAD and CelesTrak, and is widely used for satellite tracking, visualization, conjunction screening, and mission planning—though its accuracy is fundamentally limited by TLE quality and simplifying assumptions.
+
+## Ephemeris Representation
+
+### TLE
+A **Two-Line Element Set (TLE)** is a compact, legacy format (originally constrained by punch-card era line lengths) for describing satellite orbits. Despite its age, it remains widely used because it is easy to publish and often provides sufficient accuracy for many applications.
+
+In addition to the familiar (mean) Keplerian elements, TLEs include parameters related to perturbations (e.g., atmospheric drag via $B^*$, and derivatives of mean motion).
+
+TLEs are designed to be used with **SGP4**, an analytic model that produces orbital state vectors (position and velocity) from the augmented mean elements in the TLE. The perturbations modeled include Earth oblateness (which produces precession) and drag.
+
+TLEs are often preceded by an additional “line 0” containing the satellite name. For an overview, see <https://en.wikipedia.org/wiki/Two-line_element_set>.
+
+### Orbital Mean-Element Messages
+**Orbital Mean-Element Messages (OMMs)** are a more modern way of representing mean-element ephemerides. They are described by a CCSDS standard (<https://ccsds.org/Pubs/502x0b3e1.pdf>), although real-world sources may not adhere to the standard perfectly.
+
+OMMs are commonly published as:
+* JSON
+* XML
+* KVN (key–value notation)
+
+The satkit Python interface does not load OMM files directly. Instead, it expects you to provide the decoded OMM as a Python `dict` (for example, parsed from JSON or XML). The interface supports OMM dictionary layouts produced by CelesTrak and Space-Track.
+
+
+
+
+## Example Usage
+
+### SGP4 state computation from TLE
+
+```python
+import satkit as sk
+
+# The two-line element set
+# Let's pick a random Starlink satellite
+# The lines below were downloaded from https://www.celestrak.org
+tle_lines = [
+    '0 STARLINK-30477',
+    '1 57912U 23146X   24099.49439401  .00006757  00000+0  51475-3 0  9997',
+    '2 57912  43.0018 157.5807 0001420 272.5369  87.5310 15.02537576 31746'
+]
+
+# Create a TLE object
+starlink30477 = sk.TLE.from_lines(tle_lines)
+
+# We want the orbital state at April 9 2024, 12:00pm UTC
+thetime = sk.time(2024, 4, 9, 12, 0, 0)
+
+# The state is output in the "TEME" frame, which is an approximate inertial
+# frame that does not include precession or nutation
+# pTEME is geocentric position in meters
+# vTEME is geocentric velocity in meters / second
+# for now we will ignore the velocity
+pTEME, _vTEME = sk.sgp4(starlink30477, thetime)
+
+# Suppose we want currrent latitude, longitude, and altitude of satellite:
+# we need to rotate into an Earth-fixed frame, the ITRF
+# We use a "quaternion" to represent the rotation.  Quaternion rotations
+# in the satkit toolbox can be represented as multiplications of a 3-vector
+pITRF = sk.frametransform.qteme2itrf(thetime) * pTEME
+
+# Now lets make a "ITRFCoord" object to extract geodetic coordinates
+coord = sk.itrfcoord(pITRF)
+
+# Get the latitude, longitude, and
+# altitude (height above ellipsoid, or hae) of the satellite
+print(coord)
+
+# this should produce:
+# ITRFCoord(lat:  29.3890 deg, lon: 170.8051 deg, hae: 560.11 km)
+
+```
+
+### SGP4 State Computation from OMM representation of International Space Station (ISS)
+
+```python
+
+import satkit as sk
+import requests
+import json
+
+# Query the current ephemeris for the International Space Station (ISS)
+# from celestrak.org
+url = 'https://celestrak.org/NORAD/elements/gp.php?CATNR=25544&FORMAT=json'
+with requests.get(url) as response:
+    omm = response.json()
+
+# Get a representative time from the output
+epoch = sk.time(omm[0]['EPOCH'])
+# create a list of times .. once every 10 minutes
+time_array = [epoch + sk.duration(minutes=i*10) for i in range(6)]
+
+# TEME (inertial) output from SGP4
+pTEME, _vTEME = sk.sgp4(omm[0], time_array)
+
+# Rotate to Earth-fixed
+pITRF = [sk.frametransform.qteme2itrf(t) * p for t, p in zip(time_array, pTEME)]
+
+# Geodetic coordinates of space station at given times
+coord = [sk.itrfcoord(x) for x in pITRF]
+
+```

{satkit-0.8.5 → satkit-0.9.0}/python/docs/source/examples/Eclipse.ipynb

@@ -44,7 +44,7 @@
     "# Project the moon to the surface of the Earth along the Sun vector\n",
     "# This is projecting onto the surface of a sphere.  However, the Earth is actually\n",
     "# an oblate spheroid.  We can account for this by scaling the z axis by the flattening\n",
-    "# of the Earth.  \n",
+    "# of the Earth.\n",
     "scalefac = 1.0 / (1.0 - sk.consts.wgs84_f)\n",
     "moon_itrf_scaled = moon_itrf\n",
     "moon_itrf_scaled[:,2] = moon_itrf_scaled[:,2] * scalefac\n",
@@ -56,7 +56,7 @@
     "# via the law of Cosines and the quadratic equation\n",
     "lcostheta = np.sum(sun_itrf_scaled_hat * moon_itrf_scaled, axis=1)\n",
     "sqrtterm = lcostheta**2 - np.sum(moon_itrf_scaled**2, axis=1) + sk.consts.earth_radius**2\n",
-    "# Valid indices (where the projection hits the earth) are where the term under \n",
+    "# Valid indices (where the projection hits the earth) are where the term under\n",
     "# the square root is positive\n",
     "vidx = np.argwhere(sqrtterm > 0).flatten()\n",
     "\n",
@@ -171,13 +171,13 @@
     "\n",
     "    if len(eidx) > 0:\n",
     "        data[\"total\"] = {\n",
-    "            \"start\": timearr[eidx[0][0]].datetime(),  # type: ignore\n",
-    "            \"stop\": timearr[eidx[-1][0]].datetime(), # type: ignore\n",
+    "            \"begin\": timearr[eidx[0][0]].datetime(),  # type: ignore\n",
+    "            \"end\": timearr[eidx[-1][0]].datetime(), # type: ignore\n",
     "            \"duration_seconds\": (timearr[eidx[-1][0]] - timearr[eidx[0][0]]).seconds, # type: ignore\n",
     "        }\n",
     "        data[\"partial\"] = {\n",
-    "            \"start\": timearr[pidx[0][0]].datetime(),\n",
-    "            \"stop\": timearr[pidx[-1][0]].datetime(),\n",
+    "            \"begin\": timearr[pidx[0][0]].datetime(),\n",
+    "            \"end\": timearr[pidx[-1][0]].datetime(),\n",
     "            \"duration_seconds\": (timearr[pidx[-1][0]] - timearr[pidx[0][0]]).seconds, # type: ignore\n",
     "            \"peak\": None,\n",
     "            \"minangle_deg\": None,\n",
@@ -208,10 +208,10 @@
     "        max_frac_diam_occluded = 1 - (sun_extent_rad + mintheta - moon_extent_rad) / (\n",
     "            2 * sun_extent_rad\n",
     "        )\n",
-    "    \n",
+    "\n",
     "        data[\"partial\"] = {\n",
-    "            \"start\": timearr[pidx[0][0]].datetime(), # type: ignore\n",
-    "            \"stop\": timearr[pidx[-1][0]].datetime(), # type: ignore\n",
+    "            \"begin\": timearr[pidx[0][0]].datetime(), # type: ignore\n",
+    "            \"end\": timearr[pidx[-1][0]].datetime(), # type: ignore\n",
     "            \"peak\": timearr[idx].datetime(), # type: ignore\n",
     "            \"duration_seconds\": (timearr[pidx[-1][0]] - timearr[pidx[0][0]]).seconds, # type: ignore\n",
     "            \"minangle_deg\": np.min(theta) * 180.0 / m.pi,\n",
@@ -266,11 +266,11 @@
     "    [\n",
     "        {\n",
     "            \"Location\": loc[\"name\"],\n",
-    "            \"Total Eclipse Start\": loc[\"stats\"][\"total\"][\"start\"] if loc[\"stats\"][\"total\"] else None,\n",
-    "            \"Total Eclipse Stop\": loc[\"stats\"][\"total\"][\"stop\"] if loc[\"stats\"][\"total\"] else None,\n",
+    "            \"Total Eclipse Begin\": loc[\"stats\"][\"total\"][\"begin\"] if loc[\"stats\"][\"total\"] else None,\n",
+    "            \"Total Eclipse End\": loc[\"stats\"][\"total\"][\"end\"] if loc[\"stats\"][\"total\"] else None,\n",
     "            \"Total Eclipse Duration (s)\": loc[\"stats\"][\"total\"][\"duration_seconds\"] if loc[\"stats\"][\"total\"] else None,\n",
-    "            \"Partial Eclipse Start\": loc[\"stats\"][\"partial\"][\"start\"] if loc[\"stats\"][\"partial\"] else None,\n",
-    "            \"Partial Eclipse Stop\": loc[\"stats\"][\"partial\"][\"stop\"] if loc[\"stats\"][\"partial\"] else None,\n",
+    "            \"Partial Eclipse Begin\": loc[\"stats\"][\"partial\"][\"begin\"] if loc[\"stats\"][\"partial\"] else None,\n",
+    "            \"Partial Eclipse End\": loc[\"stats\"][\"partial\"][\"end\"] if loc[\"stats\"][\"partial\"] else None,\n",
     "            \"Partial Eclipse Peak\": loc[\"stats\"][\"partial\"][\"peak\"] if loc[\"stats\"][\"partial\"] else None,\n",
     "            \"Partial Eclipse Duration (s)\": loc[\"stats\"][\"partial\"][\"duration_seconds\"] if loc[\"stats\"][\"partial\"] else None,\n",
     "            \"Min Seperation (deg)\": loc[\"stats\"][\"partial\"][\"minangle_deg\"] if loc[\"stats\"][\"partial\"] else None,\n",
@@ -280,10 +280,10 @@
     "        for loc in locations\n",
     "    ]\n",
     ")\n",
-    "df.style.format({\"Total Eclipse Start\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
-    "                 \"Total Eclipse Stop\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
-    "                 \"Partial Eclipse Start\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
-    "                 \"Partial Eclipse Stop\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
+    "df.style.format({\"Total Eclipse Begin\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
+    "                 \"Total Eclipse End\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
+    "                 \"Partial Eclipse Begin\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
+    "                 \"Partial Eclipse End\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
     "                 \"Partial Eclipse Peak\": lambda x: x.strftime(\"%H:%M:%S\") if not pd.isnull(x) else 'N/A',\n",
     "                 \"Total Eclipse Duration (s)\": lambda x: f\"{x:.0f}\" if not pd.isnull(x) else 'N/A',\n",
     "                 \"Partial Eclipse Duration (s)\": lambda x: f\"{x:.0f}\" if not pd.isnull(x) else 'N/A',\n",

{satkit-0.8.5 → satkit-0.9.0}/python/docs/source/examples/High Precision Propagation.ipynb

@@ -136,8 +136,8 @@
     "\n",
     "    res = sk.propagate(\n",
     "        np.concatenate((pgcrf[0, :], v[0:3])),\n",
-    "        start=tstart,\n",
-    "        stop=tend,\n",
+    "        begin=tstart,\n",
+    "        end=tend,\n",
     "        propsettings=settings,\n",
     "        satproperties=satprops,\n",
     "    )\n",
@@ -157,8 +157,8 @@
     "satprops = sk.satproperties_static(craoverm = r.x[3])\n",
     "res = sk.propagate(\n",
     "    np.concatenate((pgcrf[0, :], r.x[0:3])),\n",
-    "    start=timearr[0],\n",
-    "    stop=timearr[-1],\n",
+    "    begin=timearr[0],\n",
+    "    end=timearr[-1],\n",
     "    propsettings=settings,\n",
     "    satproperties=satprops,\n",
     ")\n",