¿El número PRN en los archivos TLE es el mismo en los archivos RINEX? (en constelación GPS)

Quiero comparar la posición calculada con archivos Pyephem y TLE , con archivos RINEX y biblioteca GPSTK en una fecha y hora específicas, pero no tengo la certeza de si el nombre del satélite "PRN ##" es el mismo en ambos archivos ( TLE y RINEX)?

por ejemplo, en Celestrak (archivo TLE) para una constelación GPS, la identificación o el nombre del satélite es:

GPS BIIF-1 (PRN 25)
1 26690C 01004A 17096.69312500 -.00000000 00000-0 00000-0 0 966 2 26690 53.0252 160.8176 0181123 256.0966 20.90285 10

pero en los archivos RINEX solo tienen una letra unida a un número, así:

17 01 12 00 00 0.0000000 0 21G05G16G20G21 G25 G26G27G29G31R01R02R09 R15R16R17R18R19E01E26S20S36

donde la letra es la constelación G para constelación GPS y R para constelación GLONASS , S es para SBAS , etc. Junto a la letra hay un número, por ejemplo, G25 , el número de ese satélite es 25 para el archivo RINEX.

Estoy trabajado con la constelación de GPS. Estoy probando con algo de código en python, pero no estoy seguro si el PRN 25 es el mismo que el G25, porque calcular la elevación del satélite con ID 25 en un día completo es diferente en el tiempo que muestra el archivo RINEX. en comparación con los archivos TLE... o puede haber algo mal en mi código.

obsHeader, obsData = gpstk.readRinex3Obs(obsfile)
navHeader, navData = gpstk.readRinex3Nav(navfile)
# setup ephemeris store to look for satellite positions
bcestore = gpstk.GPSEphemerisStore()
for navDataObj in navData:
    ephem = navDataObj.toGPSEphemeris()
    bcestore.addEphemeris(ephem)
bcestore.SearchNear()
navData.close()
#look each satellite
for obsObject in obsData: #round the data
    for satID, datumList in obsObject.obs.iteritems():
        if satID.systemGPS == satID.system and satID.id == NumSatellite:
            time=obsObject.time #satellite time from the observation
            CivilTime = gpstk.CivilTime(obsObject.time)
            #break
            eph   = bcestore.findEphemeris(satID, obsObject.time)
            svXvt = eph.svXvt(obsObject.time)
            elev = obsHeader.antennaPosition.elvAngle(svXvt.getPos())
            azim    = obsHeader.antennaPosition.azAngle(svXvt.getPos())
            
            SecondsTime.append(obsObject.time.getSecondOfDay())
            TimesSatRin.append(CivilTime)
            AziSatRin.append(azim)
            ElevSatRin.append(elev)
#for each observation day.
DateTimeRinex = CivilTime2datetime(TimesSatRin)
for t in DateTimeRinex:
    #for each Time of the rinex (for the specific Satellite)
    print 'next time: ', t

    names,elevs,azs,types = GetPositionELAZ(t,Longitude_obs,Latitude_obs,1) #just GPS with TLE files
    print names
    for i in range (0,len(names)):
        number = int(filter(str.isdigit, names[i]))
        #print number,' number sat of the TLE'
        if number == NumSatellite:
            print 'GOT',number,' number sat of the TLE'
            TimeTLE.append(t)
            ElevTLE.append(elevs)
            AzimTLE.append(azs)

La posición (lon,Lat) es la misma de la estación terrestre del repositorio de archivos RINEX. Gracias.

ACTUALIZACIÓN 1:
Estoy probando con mi código, un vehículo por uno, y puedo ver la comparación de elevaciones de cada satélite a lo largo del día. Para los archivos RINEX pruebo con 4 días, día del año 101, 100, 80 y 10. y para estos días siempre el TLS PRN 10 los resultados son similares a RINEX PRN 12 ; Dije similar, porque no son los mismos resultados. Todavía estoy comprobando.

20 de abril

19 de abril

21 de marzo

10 de enero

Aparentemente existe una ID de PRN para RINEX por cada ID de PRN para TLE, en este caso RINEX PRN 12 = TLE PRN 10. Pero esta no es la respuesta precisa, así que dejaré abierta la pregunta.

ACTUALIZACIÓN 2:

El cálculo de posición con TLE es muy preciso. pero me dio un error, me dio una posición en Longitud/Latitud, pero era una posición incorrecta, por cierto, debe tomar la posición para comparar, por el archivo RINEX (el archivo de observación), en la Cabecera; el archivo tiene la posición en X,Y,Z, por lo que es necesario convertir a posición longitud-latitud.

para obtener la Longitud-Latitud:

def Get_latitude(x, y, z):
    #radius earth
    r = 6371000.
    return np.arcsin(z/r)*(180/np.pi)

def Get_longitude(x, y, z):
    #radius earth
    r = 6371000.
    if (x > 0):
        return np.arctan(y/x)*(180./np.pi)
    elif y > 0:
        return np.arctan(y/x)*(180./np.pi) + 180.
    else:
        return np.arctan(y/x)*(180./np.pi) - 180.