Sentinel-3 RPS experiment Python package

In [1]:
# import packages
import numpy as np
import sentinel3 as s3

Parameters

In [2]:
t0 = np.datetime64('2018-08-31T00:00:00')
# GPS satellite
sv_name = 'G01'
sv_sp3_import = True
sv_sp3_file = 'igs20165'
sv_sp3_dir = '../data/test/'
# Sentinel-3 orbit file
S3_orbit_file = 'S3A_OPER_AUX_POEORB_POD__20180925T071034_V20180830T215942_20180831T235942_DGNS'
S3_orbit_dir = '../data/test/'
# Sentinel-3 rinex file
rinex_freq = 'C1C'
rinex_file = 'S3A_OPER_AUX_GNSSRD_POD__20180907T024226_V20180830T235942_20180831T235941'
rinex_dir = '../data/test/'
# output
output_dir = '../output/'

Reading data

In [3]:
rinex_data, S3_orbit, gps_orbit = s3.read_data(t0, sv_name, output_dir,
                rinex_dir, rinex_file, 
                S3_orbit_dir, S3_orbit_file,
                sv_sp3_dir, sv_sp3_file,
                read_EOF=True, read_DBL=False, read_nc=True, read_sp3=True)
Reading netcdf Rinex file: ../output/S3A_OPER_AUX_GNSSRD_POD__20180907T024226_V20180830T235942_20180831T235941.nc
List of sats in view: ['G01' 'G02' 'G03' 'G05' 'G06' 'G07' 'G08' 'G09' 'G10' 'G11' 'G12' 'G13'
 'G14' 'G15' 'G16' 'G17' 'G18' 'G19' 'G20' 'G21' 'G22' 'G23' 'G24' 'G25'
 'G26' 'G27' 'G28' 'G29' 'G30' 'G31' 'G32']
Reading Sentinel-3 orbit file: ../data/test/S3A_OPER_AUX_POEORB_POD__20180925T071034_V20180830T215942_20180831T235942_DGNS.EOF
Reading sp3 orbit file: ../data/test/igs20165.sp3