Questions tagged [skyfield]

Skyfield high precision scientific astronomy library for Python and NumPy

The Skyfield project provides a modern alternative to the elderly PyEphem astronomy library for Python. Instead of requiring a custom C extension to be built or installed, Skyfield does all of its computation in pure Python and requires only the NumPy high performance numeric library in order to do its work — making Skyfield perfect for developers who can install NumPy directly, or who use scientific distributions that install NumPy automatically.

93 questions
0
votes
1 answer

Skyfield satellite position wrong?

I am trying to find the position of the MODIS AQUA satellite at arbitrary times. I found the skyfield package for this, however comparisons with NASA's LAADS page are wildly different. Am I doing something wrong with the Skyfield package, or is…
0
votes
0 answers

Computing apparent magnitude of Algol (Python)

I'm trying to use Python to compute the apparent magnitude of the star Algol. Wikipedia: Algol's magnitude is usually near-constant at 2.1, but regularly dips to 3.4 every 2.86 days during the roughly 10-hour-long partial eclipses.) Where's the…
velkyvont
  • 23
  • 3
0
votes
0 answers

Calculating Satellite orbit prpagation, plotting and minimal distance

I want to programm a Code which calculates an Orbit for 2 Satellites to recreate an event. It should read in 2 TLE files, then create Earthsatellite Objects with the skyfield package in python. After that it should search the TLE files for a…
G3nki4
  • 11
0
votes
1 answer

Determine Planet Conjunction Times per Year using Skyfield & Scipy

Expanding on the solution provided here: Determine coordinates at conjunction times ... I coded the following to give me all the conjunctions of 5 planets (and the sun) for any given year (within the ephemeris, of course), sorted by date. CORRECTED…
Aendie
  • 321
  • 1
  • 13
0
votes
1 answer

How do you calculate a satellite's position in GCRF from an RA/DEC measurement in Skyfield?

I have a measurement of the RA and Dec for an Earth orbiting satellite, as measured from a sensor on the Earth's surface. I'm trying to calculate the satellite's position vector in the GCRF reference frame (so a vector from the centre of the…
Michael
  • 31
  • 4
0
votes
1 answer

Skyfield timescale object does not include leapseconds during creation of a UTC date?

Say I have a unix epoch time stamp 1664203669. When generating a UTC date from this time stamp using datetime I get the following date. import datetime datetime.datetime.utcfromtimestamp(1664203669).strftime("%Y-%m-%d %H:%M:%S") Out[25]: '2022-09-26…
KDecker
  • 6,928
  • 8
  • 40
  • 81
0
votes
1 answer

Skyfield returning incorrect latitude and longitude coordinates?

I'm trying to make a program which takes latitude and longitude coordinates of the ISS (from the Skyfield API) and plots it onto a scatterplot (using matplotlib) to show its estimated orbit path. The problem is that when it's plotted onto the…
bisho
  • 3
  • 1
0
votes
1 answer

Iterate trough a converted datetime pandas dataframe with a external function

https://rhodesmill.org/skyfield/positions.html#azimuth-and-altitude-from-a-geographic-position Hi I have function that generates a sun-shot azimuth on a specific date and time on a specific place, using the package skyfield for astronomical…
L Blauw
  • 1
  • 1
0
votes
1 answer

Discrepancy between PyEphem and Skyfield HA/dec results

I'm refactoring some old code that used PyEphem to use Skyfield, and I'm getting a slight difference in the results of the GHA/dec of a body. def sf(year): from skyfield.api import N, W, wgs84 from skyfield.api import load from…
MerseyViking
  • 389
  • 3
  • 19
0
votes
1 answer

Satellite velocity from km/s to arcsec

I can't change the current satellite velocity from km/s units to arcsec units. This is a snippet of my code: from astropy import units as u def check_the_satellites_current_velocity(self, available_satellites): satellites_velocity = {} for…
0
votes
0 answers

Night flight time calculation: IndexError: invalid index to scalar variable

I'm using vsimakhin's night time calculator (https://github.com/vsimakhin/night-time-calculator) to calculate total night time flown during a flight given the departure and arrival locations and times. I am new to Python so i'm not sure what the…
Bruno
  • 511
  • 2
  • 6
  • 19
0
votes
0 answers

AttributeError: 'list' object has no attribute 'center' skyfiled wgs84.geographic_position_of ()

I'm trying to determine the latitude-longitude position on Earth of a position using the skifield python package from skyfield.api import wgs84 position = [7033.45648791, 178.93192295, 0.00290722] ground =…
0
votes
1 answer

I am having problems with skyfield python

I want to find out the heliocentric position of jupiter, and the corresponding zodiac but i am having basic problems. from skyfield.api import sun, jupiter p = sun(utc=(1980, 1, 1)).observe(jupiter) lat, lon, d = p.ecliptic_latlon() …
Agent π
  • 33
  • 1
  • 8
0
votes
2 answers

How to find Planet Position using skyfield

I want to find planet position using skyfield in the form of Ecliptic longitude,latitude, speed and distance in geocentric
0
votes
1 answer

Degree of Solar Longitude vs Earth to calculate Equinox and Solstice ,etc

I try build some function to calculate Solar longitude, commonly abbreviated as Ls, is the ecliptic longitude of the sun, i.e. the position of the sun on the celestial sphere along the ecliptic. For example when input date 07/Nov/2021 it will return…