PRN-Nummer in TLE-Dateien ist die gleiche wie in RINEX-Dateien? (in GPS-Konstellation)

Ich möchte die mit Pyephem- und TLE -Dateien berechnete Position mit RINEX- Dateien und der GPSTK- Bibliothek an einem bestimmten Datum und einer bestimmten Uhrzeit vergleichen , aber ich habe nicht die Gewissheit, ob der Name des Satelliten "PRN ##" in beiden Dateien gleich ist ( TLE und RINEX)?

zB in Celestrak (TLE-Datei) für eine GPS-Konstellation lautet die ID oder der Name des Satelliten:

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.3055 20.3285

aber in den RINEX-Dateien ist nur ein Buchstabe mit einer Zahl verbunden, wie folgt:

17 01 12 00 00 0,0000000 0 21G05G16G20G21 G25 G26G27G29G31R01R02R09 R15R16R17R18R19E01E26S20S36

Dabei steht der Buchstabe für die Konstellation G für die GPS- Konstellation und R für die GLONASS- Konstellation, S für SBAS usw. Neben dem Buchstaben steht eine Zahl, z. B. G25 , die Nummer dieses Satelliten ist 25 für die RINEX-Datei.

Ich habe mit GPS-Konstellation gearbeitet. Ich teste mit etwas Code in Python, aber ich bin mir nicht sicher, ob die PRN 25 mit G25 identisch ist?, weil die Berechnung der Höhe des Satelliten mit der ID 25 an einem ganzen Tag in der Zeit, die die RINEX-Datei zeigt, unterschiedlich ist Vergleich mit TLE-Dateien ... oder vielleicht stimmt etwas in meinem Code nicht.

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)

Die Position (lon, Lat) ist die gleiche wie die der Bodenstation des RINEX-Dateispeichers. Danke.

UPDATE 1:
Ich teste mit meinem Code, ein Fahrzeug nach dem anderen, und ich kann den Vergleich der Höhen jedes Satelliten über den Tag hinweg sehen. Für RINEX-Dateien versuche ich es mit 4 Tagen, Tag des Jahres 101, 100, 80 und 10. und für diese Tage immer die TLS PRN 10 , die Ergebnisse sind ähnlich wie bei RINEX PRN 12 ; Ich sagte ähnlich, weil es nicht die gleichen Ergebnisse sind. Ich überprüfe noch.

20. April

19. April

21. März

10. Januar

Anscheinend existiert neben jeder PRN-ID für TLE eine PRN-ID für RINEX, in diesem Fall RINEX PRN 12 = TLE PRN 10. Dies ist jedoch nicht die genaue Antwort, daher lasse ich die Frage offen

UPDATE 2:

Die Positionsberechnung mit TLE ist sehr genau. aber ich habe einen Fehler bekommen, ich habe eine Position in Längengrad/Breitengrad erhalten, aber es war eine falsche Position, sicher müssen Sie die Position zum Vergleichen von der RINEX-Datei (der Beobachtungsdatei) in der Kopfzeile nehmen; Die Datei hat die Position in X, Y, Z, daher ist eine Konvertierung in die Position Längengrad-Breitengrad erforderlich.

um den Längengrad-Breitengrad zu erhalten:

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.