jplephem

https://pypi.python.org/pypi/jplephem
Admin
This package can load and use a Jet Propulsion Laboratory (JPL) ephemeris for predicting the...
Tokens:4,313
Snippets:74
Trust Score:10
Update:4 months ago
Show doc for...
Tokens:
Raw
### Install jplephem

Source: https://pypi.python.org/pypi/jplephem

Install the jplephem package using pip. Pip will automatically install NumPy as a dependency.

```bash
$ pip install jplephem
```

--------------------------------

### Create Ephemeris Excerpt by Date Range

Source: https://pypi.python.org/pypi/jplephem

Generate a smaller ephemeris file by specifying a start and end date range. The output file will contain only the data within this period.

```bash
python -m jplephem excerpt 2018/1/1 2018/4/1 de421.bsp excerpt421.bsp
```

--------------------------------

### Handle Negative Start Years in Excerpts

Source: https://pypi.python.org/pypi/jplephem

When specifying a negative start year for an excerpt, use the '--' argument to indicate the end of options. This prevents negative years from being misinterpreted as command-line flags.

```bash
python -m jplephem excerpt -- -800/1/1 800/1/1 de422.bsp excerpt422.bsp
```

--------------------------------

### Configure NumPy for Pretty Output

Source: https://pypi.python.org/pypi/jplephem

Set NumPy print options for better readability of vector outputs. This is useful for visualizing results.

```python
import numpy as np
np.set_printoptions(precision=3)
```

--------------------------------

### Load and Inspect SPK Kernel

Source: https://pypi.python.org/pypi/jplephem

Load an SPK ephemeris file and print its segment information. Ensure the kernel file is downloaded before opening.

```python
from jplephem.spk import SPK
kernel = SPK.open('de421.bsp')
print(kernel)
```

--------------------------------

### Load and Inspect Binary PCK File

Source: https://pypi.python.org/pypi/jplephem

Opens a binary PCK file and accesses its segments to retrieve body, frame, and data type information. Ensure the PCK file exists at the specified path.

```python
from jplephem.pck import PCK
p = PCK.open('moon_pa_de421_1900-2050.bpc')
print(p.segments[0].body)
print(p.segments[0].frame)
print(p.segments[0].data_type)
```

--------------------------------

### Display Ephemeris File Comments and Data

Source: https://pypi.python.org/pypi/jplephem

Use these commands to display the comment text or the DAF/SPK data contained within a .bsp ephemeris file.

```bash
python -m jplephem comment de421.bsp
```

```bash
python -m jplephem daf de421.bsp
```

```bash
python -m jplephem spk de421.bsp
```

```bash
python -m jplephem spk -v de421.bsp
```

--------------------------------

### Load Raw Coefficients from Segment

Source: https://pypi.python.org/pypi/jplephem

Loads the raw polynomial coefficients from an SPK segment using the load_array() method. Returns the initial epoch, interval length, and the coefficients array.

```python
>>> initial_epoch, interval_length, coefficients = segment.load_array()
>>> print(coefficients.shape)
(3, 14080, 13)
```

--------------------------------

### Create Ephemeris Excerpt by Target IDs

Source: https://pypi.python.org/pypi/jplephem

Filter the ephemeris data to include only specific target IDs within a given date range. Unrecognized targets are ignored, allowing for flexible filtering.

```bash
python -m jplephem excerpt --targets 1,2,3 2018/1/1 2018/4/1 de421.bsp excerpt421.bsp
```

--------------------------------

### Download and Excerpt Remote Ephemeris Data

Source: https://pypi.python.org/pypi/jplephem

Fetch data from a remote ephemeris file (URL) and create an excerpt. This conserves bandwidth by downloading only the necessary data blocks for the specified date range.

```bash
$ python -m jplephem excerpt 2018/1/1 2018/4/1 \
    https://naif.jpl.nasa.gov/pub/naif/generic_kernels/spk/satellites/jup365.bsp \
    excerpt.bsp
```

```bash
$ ls -lh excerpt.bsp
```

--------------------------------

### Compute Position and Velocity

Source: https://pypi.python.org/pypi/jplephem

Calculate both the position and velocity of a celestial body. Velocity is returned as a second value from the compute_and_differentiate method.

```python
position, velocity = kernel[0,4].compute_and_differentiate(2457061.5)
print(position)
print(velocity)
```

--------------------------------

### Describe SPK Segment

Source: https://pypi.python.org/pypi/jplephem

Prints a verbose description of an SPK segment, including its date range, type, frame, and source.

```python
>>> segment = kernel[3,399]
>>> print(segment.describe())
1899-07-29..2053-10-09  Type 2  Earth Barycenter (3) -> Earth (399)
  frame=1 source=DE-0421LE-0421
```

--------------------------------

### Convert Velocity to Per Second

Source: https://pypi.python.org/pypi/jplephem

Convert the computed velocity (distance per day) to velocity per second by dividing by the number of seconds in a day.

```python
velocity_per_second = velocity / 86400.0
print(velocity_per_second)
```

--------------------------------

### BaseSegment Attributes

Source: https://pypi.python.org/pypi/jplephem

Lists the available attributes for a BaseSegment object, providing details about the ephemeris data.

```python
>>> from jplephem.spk import BaseSegment
>>> help(BaseSegment)
Help on class BaseSegment in module jplephem.spk:
...
 |  segment.source - official ephemeris name, like 'DE-0430LE-0430'
 |  segment.start_second - initial epoch, as seconds from J2000
 |  segment.end_second - final epoch, as seconds from J2000
 |  segment.start_jd - start_second, converted to a Julian Date
 |  segment.end_jd - end_second, converted to a Julian Date
 |  segment.center - integer center identifier
 |  segment.target - integer target identifier
 |  segment.frame - integer frame identifier
 |  segment.data_type - integer data type identifier
 |  segment.start_i - index where segment starts
 |  segment.end_i - index where segment ends
...
```

--------------------------------

### Compute Position with Multiple Subtractions

Source: https://pypi.python.org/pypi/jplephem

Calculate the position of one body relative to another by chaining subtractions through barycenters. This demonstrates complex relative positioning.

```python
position = kernel[0,4].compute(2457061.5)
position -= kernel[0,3].compute(2457061.5)
position -= kernel[3,399].compute(2457061.5)
print(position)
```

--------------------------------

### Compute Celestial Body Position

Source: https://pypi.python.org/pypi/jplephem

Calculate the position of a celestial body relative to another for a specific Julian date. The output is in kilometers.

```python
position = kernel[0,4].compute(2457061.5)
print(position)
```

--------------------------------

### Compute Rotation and Velocity Data from PCK Segment

Source: https://pypi.python.org/pypi/jplephem

Computes both the rotation angles and their velocities for a given Julian date from a PCK segment. The first returned array contains the rotation data, and the second contains the velocity data.

```python
tdb = 2454540.34103
r, v = p.segments[0].compute(tdb, 0.0, True)
print(r)
print(v)
```

--------------------------------

### Compute Positions for Multiple Julian Dates

Source: https://pypi.python.org/pypi/jplephem

Compute celestial body positions for an array of Julian dates. This returns a vector for each component, matching the input date array length.

```python
jd = np.array([2457061.5, 2457062.5, 2457063.5, 2457064.5])
position = kernel[0,4].compute(jd)
print(position)
```

--------------------------------

### Convert Gregorian to Julian Date

Source: https://pypi.python.org/pypi/jplephem

Convert a Gregorian calendar date to its Julian date equivalent. This is necessary for ephemeris computations.

```python
from jplephem.calendar import compute_julian_date
compute_julian_date(2015, 2, 8)
```

--------------------------------

### Compute Rotation Data from PCK Segment

Source: https://pypi.python.org/pypi/jplephem

Calculates the three angles (right ascension of the pole, declination of the pole, and cumulative rotation) for a given Julian date from a PCK segment. The output is typically in radians. Set the third argument to True to also compute velocity.

```python
tdb = 2454540.34103
print(p.segments[0].compute(tdb, 0.0, False))
```

--------------------------------

### Close Ephemeris Files

Source: https://pypi.python.org/pypi/jplephem

Releases all open files and memory maps associated with an ephemeris object. This should be called when the ephemeris is no longer needed to free up resources.

```python
kernel.close()
p.close()
```
Tokens:
Raw