Evaluation with the skyfield package of the azimuth and elevation of a satellite with respect to a GPS position

[dellache]

Hello,

I am trying to evaluate with the skyfield package the azimuth and elevation of a satellite with respect to a GPS position.

To do that, i use the following Python code (simple example below) :

from skyfield.api import load, wgs84, EarthSatellite

ts = load.timescale()
t = ts.now()

line1 = 'ISS (ZARYA)'
line2 = '1 25544U 98067A   24100.54831018  .00016717  00000+0  30114-3 0  9996'
line3 = '2 25544  51.6405 132.6712 0002394  85.9768  37.9614 15.50655812452923'

sat = EarthSatellite(line2, line3, line1, ts)
observer = wgs84.latlon(latitude_degrees=35, longitude_degrees=-60, elevation_m=0.0)

topocentric = observer.at(t).observe(sat).apparent()

alt, az, distance = topocentric.altaz()

print(f"Élévation (altitude) : {alt.degrees:.2f}°")
print(f"Azimut : {az.degrees:.2f}°")

But I have the following error :

    topocentric = observer.at(t).observe(sat).apparent()
                  ^^^^^^^^^^^^^^^^^^^^^^
AttributeError: 'Geocentric' object has no attribute 'observe'

Note that I use the 1.52 version of the skyfield package.

Could you help me ?

Thanks very much !

Stéphane

[dellache]

I found the solution by adapting a part of the code proposed in the post #1062 to what I am looking for.