Eu quero converter posições x / y / z-ECEF para lla (lat / lon / alt) usando WGS84 em python com pyproj, mas parece que a conversão falha.

O código de exemplo está aqui:

import pyproj

# Example position data, should be somewhere in Germany
x = 652954.1006
y = 4774619.7919
z = -2217647.7937

ecef = pyproj.Proj(proj='geocent', ellps='WGS84', datum='WGS84')
lla = pyproj.Proj(proj='latlong', ellps='WGS84', datum='WGS84')
lon, lat, alt = pyproj.transform(ecef, lla, x, y, z, radians=True)

print lat, lon, alt

Alguém pode ver onde está o problema?

Edição: Até agora, eu acho que os cálculos estão corretos, apenas os dados que recebo do meu receptor parece estar com defeito. Alguém pode confirmar isso?