mirrors/gnss-sdr
1
0
Fork 0
mirror of https://github.com/gnss-sdr/gnss-sdr synced 2025-10-31 07:13:03 +00:00
Code Issues Releases Wiki Activity

Fix some merging problems

Browse Source
This commit is contained in:
Carles Fernandez
2018-03-10 12:56:30 +01:00
parent 074fb8b9f0
commit e648f38f6f
3 changed files with 105 additions and 107 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
conf/gnss-sdr_GLONASS_L1_CA_GPS_L1_CA_ibyte.conf
View File

@@ -6,7 +6,7 @@ Receiver.sources_count=2
;######### SIGNAL_SOURCE CONFIG ############ ;######### SIGNAL_SOURCE CONFIG ############
SignalSource0.implementation=File_Signal_Source SignalSource0.implementation=File_Signal_Source
SignalSource0.filename=/media/dmiralles/Seagate Backup Plus Drive/GNSS Data/NT1065_L1_20160923_fs6625e6_if60e3_schar.bin ; <- PUT YOUR FILE HERE SignalSource0.filename=/archive/NT1065_L1_20160923_fs6625e6_if60e3_schar.bin ; <- PUT YOUR FILE HERE
SignalSource0.item_type=ibyte SignalSource0.item_type=ibyte
SignalSource0.sampling_frequency=6625000 SignalSource0.sampling_frequency=6625000
SignalSource0.samples=0 SignalSource0.samples=0
@@ -14,7 +14,7 @@ SignalSource0.dump=false;
SignalSource0.dump_filename=/archive/signal_glonass.bin SignalSource0.dump_filename=/archive/signal_glonass.bin
SignalSource1.implementation=File_Signal_Source SignalSource1.implementation=File_Signal_Source
SignalSource1.filename=/media/dmiralles/Seagate Backup Plus Drive/GNSS Data/NT1065_GLONASS_L1_20160923_fs6625e6_if0e3_schar.bin ; <- PUT YOUR FILE HERE SignalSource1.filename=/archive/NT1065_GLONASS_L1_20160923_fs6625e6_if0e3_schar.bin ; <- PUT YOUR FILE HERE
SignalSource1.item_type=ibyte SignalSource1.item_type=ibyte
SignalSource1.sampling_frequency=6625000 SignalSource1.sampling_frequency=6625000
SignalSource1.samples=0 SignalSource1.samples=0

2
conf/gnss-sdr_GLONASS_L1_CA_ibyte.conf
View File

@@ -5,7 +5,7 @@ GNSS-SDR.internal_fs_sps=6625000
;######### SIGNAL_SOURCE CONFIG ############ ;######### SIGNAL_SOURCE CONFIG ############
SignalSource.implementation=File_Signal_Source SignalSource.implementation=File_Signal_Source
SignalSource.filename=/media/dmiralles/Seagate Backup Plus Drive/GNSS Data/NT1065_GLONASS_L1_20160924_fs6625e6_if0e3_schar.bin ; <- PUT YOUR FILE HERE SignalSource.filename=/archive/NT1065_GLONASS_L1_20160924_fs6625e6_if0e3_schar.bin ; <- PUT YOUR FILE HERE
SignalSource.item_type=ibyte SignalSource.item_type=ibyte
SignalSource.sampling_frequency=6625000 SignalSource.sampling_frequency=6625000
SignalSource.samples=0 SignalSource.samples=0

206
src/algorithms/PVT/libs/rinex_printer.cc
View File

@@ -47,8 +47,6 @@
using google::LogMessage; using google::LogMessage;
DEFINE_string(RINEX_version, "3.02", "Specifies the RINEX version (2.11 or 3.02)");
Rinex_Printer::Rinex_Printer(int conf_version) Rinex_Printer::Rinex_Printer(int conf_version)
{ {
@@ -3960,31 +3958,31 @@ void Rinex_Printer::rinex_obs_header(std::fstream& out, const Gps_CNAV_Ephemeris
out << line << std::endl; out << line << std::endl;
// -------- Line MARKER NAME // -------- Line MARKER NAME
line.clear(); line.clear();
line += Rinex_Printer::leftJustify("DEFAULT MARKER NAME", 60); // put a flag or a property, line += Rinex_Printer::leftJustify("DEFAULT MARKER NAME", 60); // put a flag or a property,
line += Rinex_Printer::leftJustify("MARKER NAME", 20); line += Rinex_Printer::leftJustify("MARKER NAME", 20);
Rinex_Printer::lengthCheck(line); Rinex_Printer::lengthCheck(line);
out << line << std::endl; out << line << std::endl;
// -------- Line MARKER NUMBER / TYPE // -------- Line MARKER NUMBER / TYPE
if (version == 2) if (version == 2)
{ {
line.clear(); line.clear();
line += Rinex_Printer::leftJustify("GROUND_CRAFT", 20); // put a flag or a property line += Rinex_Printer::leftJustify("GROUND_CRAFT", 20); // put a flag or a property
line += std::string(40, ' '); line += std::string(40, ' ');
line += Rinex_Printer::leftJustify("MARKER NUMBER", 20); line += Rinex_Printer::leftJustify("MARKER NUMBER", 20);
Rinex_Printer::lengthCheck(line); Rinex_Printer::lengthCheck(line);
out << line << std::endl; out << line << std::endl;
} }
if (version == 3) if (version == 3)
{ {
line.clear(); line.clear();
line += Rinex_Printer::leftJustify("GROUND_CRAFT", 20); // put a flag or a property line += Rinex_Printer::leftJustify("GROUND_CRAFT", 20); // put a flag or a property
line += std::string(40, ' '); line += std::string(40, ' ');
line += Rinex_Printer::leftJustify("MARKER TYPE", 20); line += Rinex_Printer::leftJustify("MARKER TYPE", 20);
Rinex_Printer::lengthCheck(line); Rinex_Printer::lengthCheck(line);
out << line << std::endl; out << line << std::endl;
} }
// -------- Line OBSERVER / AGENCY // -------- Line OBSERVER / AGENCY
line.clear(); line.clear();
@@ -4279,7 +4277,7 @@ void Rinex_Printer::rinex_obs_header(std::fstream& out, const Galileo_Ephemeris&
// -------- Line MARKER NAME // -------- Line MARKER NAME
line.clear(); line.clear();
line += Rinex_Printer::leftJustify("DEFAULT MARKER NAME", 60); // put a flag or a property, line += Rinex_Printer::leftJustify("DEFAULT MARKER NAME", 60); // put a flag or a property,
line += Rinex_Printer::leftJustify("MARKER NAME", 20); line += Rinex_Printer::leftJustify("MARKER NAME", 20);
Rinex_Printer::lengthCheck(line); Rinex_Printer::lengthCheck(line);
out << line << std::endl; out << line << std::endl;
@@ -6011,20 +6009,20 @@ void Rinex_Printer::log_rinex_obs(std::fstream& out, const Glonass_Gnav_Ephemeri
//double utc_t = nav_msg.utc_time(nav_msg.sv_clock_correction(obs_time)); //double utc_t = nav_msg.utc_time(nav_msg.sv_clock_correction(obs_time));
//double gps_t = eph.sv_clock_correction(obs_time); //double gps_t = eph.sv_clock_correction(obs_time);
std::string month (timestring, 4, 2); std::string month(timestring, 4, 2);
std::string day (timestring, 6, 2); std::string day(timestring, 6, 2);
std::string hour (timestring, 9, 2); std::string hour(timestring, 9, 2);
std::string minutes (timestring, 11, 2); std::string minutes(timestring, 11, 2);
double utc_sec = modf (obs_time , &int_sec) + p_glonass_time.time_of_day().seconds() ; double utc_sec = modf(obs_time, &int_sec) + p_glonass_time.time_of_day().seconds();
if (version == 2) if (version == 2)
{ {
line.clear(); line.clear();
std::string year (timestring, 2, 2); std::string year(timestring, 2, 2);
line += std::string(1, ' '); line += std::string(1, ' ');
line += year; line += year;
line += std::string(1, ' '); line += std::string(1, ' ');
if (month.compare(0, 1 , "0") == 0) if (month.compare(0, 1, "0") == 0)
{ {
line += std::string(1, ' '); line += std::string(1, ' ');
line += month.substr(1, 1); line += month.substr(1, 1);
@@ -6034,7 +6032,7 @@ void Rinex_Printer::log_rinex_obs(std::fstream& out, const Glonass_Gnav_Ephemeri
line += month; line += month;
} }
line += std::string(1, ' '); line += std::string(1, ' ');
if (day.compare(0, 1 , "0") == 0) if (day.compare(0, 1, "0") == 0)
{ {
line += std::string(1, ' '); line += std::string(1, ' ');
line += day.substr(1, 1); line += day.substr(1, 1);
@@ -6152,9 +6150,9 @@ void Rinex_Printer::log_rinex_obs(std::fstream& out, const Glonass_Gnav_Ephemeri
line += std::string(1, ' '); line += std::string(1, ' ');
// Add extra 0 if seconds are < 10 // Add extra 0 if seconds are < 10
if (utc_sec < 10) if (utc_sec < 10)
{ {
line += std::string(1, '0'); line += std::string(1, '0');
} }
line += Rinex_Printer::asString(utc_sec, 7); line += Rinex_Printer::asString(utc_sec, 7);
line += std::string(2, ' '); line += std::string(2, ' ');
// Epoch flag 0: OK 1: power failure between previous and current epoch <1: Special event // Epoch flag 0: OK 1: power failure between previous and current epoch <1: Special event
@@ -8121,10 +8119,10 @@ void Rinex_Printer::to_date_time(int gps_week, int gps_tow, int& year, int& mont
} }
// find month // find month
for (int m = 1; true ; m++) for (int m = 1; true; m++)
{ {
int secs_in_month_m = days_per_month[m-1]*secs_per_day; int secs_in_month_m = days_per_month[m - 1] * secs_per_day;
if (is_leap_year && m == 2 ) // consider February of leap year if (is_leap_year && m == 2) // consider February of leap year
{ {
secs_in_month_m += secs_per_day; secs_in_month_m += secs_per_day;
} }
@@ -8142,14 +8140,14 @@ void Rinex_Printer::to_date_time(int gps_week, int gps_tow, int& year, int& mont
//std::cout << "month: m=" << m << " secs_in_month_m="<< secs_in_month_m << " remaining_secs="<< remaining_secs << std::endl; //std::cout << "month: m=" << m << " secs_in_month_m="<< secs_in_month_m << " remaining_secs="<< remaining_secs << std::endl;
} }
day = remaining_secs/secs_per_day+1; day = remaining_secs / secs_per_day + 1;
remaining_secs = remaining_secs%secs_per_day; remaining_secs = remaining_secs % secs_per_day;
hour = remaining_secs/(60*60); hour = remaining_secs / (60 * 60);
remaining_secs = remaining_secs%(60*60); remaining_secs = remaining_secs % (60 * 60);
minute = remaining_secs/60; minute = remaining_secs / 60;
second = remaining_secs%60; second = remaining_secs % 60;
} }
@@ -8296,82 +8294,82 @@ boost::posix_time::ptime Rinex_Printer::compute_Galileo_time(const Galileo_Ephem
boost::posix_time::ptime Rinex_Printer::compute_UTC_time(const Glonass_Gnav_Ephemeris& eph, const double obs_time) boost::posix_time::ptime Rinex_Printer::compute_UTC_time(const Glonass_Gnav_Ephemeris& eph, const double obs_time)
{ {
double tod = 0.0; double tod = 0.0;
double glot2utc = 3*3600; double glot2utc = 3 * 3600;
double obs_time_utc = 0.0, obs_time_glot = 0.0; double obs_time_utc = 0.0, obs_time_glot = 0.0;
int i = 0; int i = 0;
// Get observation time in nearly GLONASS time. Correction for leap seconds done at the end // Get observation time in nearly GLONASS time. Correction for leap seconds done at the end
obs_time_glot = obs_time + glot2utc; obs_time_glot = obs_time + glot2utc;
// Get seconds of day in glonass time // Get seconds of day in glonass time
tod = fmod (obs_time_glot, 86400); tod = fmod(obs_time_glot, 86400);
// Form date and time duration types // Form date and time duration types
boost::posix_time::time_duration t1(0, 0, tod); boost::posix_time::time_duration t1(0, 0, tod);
boost::gregorian::date d1(eph.d_yr, 1, 1); boost::gregorian::date d1(eph.d_yr, 1, 1);
boost::gregorian::days d2(eph.d_N_T-1); boost::gregorian::days d2(eph.d_N_T - 1);
boost::posix_time::ptime glo_time(d1 + d2, t1); boost::posix_time::ptime glo_time(d1 + d2, t1);
// Convert to utc // Convert to utc
boost::posix_time::time_duration t2(0, 0, glot2utc); boost::posix_time::time_duration t2(0, 0, glot2utc);
boost::posix_time::ptime utc_time = glo_time - t2; boost::posix_time::ptime utc_time = glo_time - t2;
// Adjust for leap second correction // Adjust for leap second correction
for (i = 0; GLONASS_LEAP_SECONDS[i][0]>0; i++) for (i = 0; GLONASS_LEAP_SECONDS[i][0] > 0; i++)
{ {
boost::posix_time::time_duration t3(GLONASS_LEAP_SECONDS[i][3], GLONASS_LEAP_SECONDS[i][4], GLONASS_LEAP_SECONDS[i][5]); boost::posix_time::time_duration t3(GLONASS_LEAP_SECONDS[i][3], GLONASS_LEAP_SECONDS[i][4], GLONASS_LEAP_SECONDS[i][5]);
boost::gregorian::date d3(GLONASS_LEAP_SECONDS[i][0], GLONASS_LEAP_SECONDS[i][1], GLONASS_LEAP_SECONDS[i][2]); boost::gregorian::date d3(GLONASS_LEAP_SECONDS[i][0], GLONASS_LEAP_SECONDS[i][1], GLONASS_LEAP_SECONDS[i][2]);
boost::posix_time::ptime ls_time(d3, t3); boost::posix_time::ptime ls_time(d3, t3);
if (utc_time >= ls_time) if (utc_time >= ls_time)
{ {
// We subtract the leap second when going from gpst to utc, values store as negatives // We subtract the leap second when going from gpst to utc, values store as negatives
utc_time = utc_time + boost::posix_time::time_duration(0,0,GLONASS_LEAP_SECONDS[i][6]); utc_time = utc_time + boost::posix_time::time_duration(0, 0, GLONASS_LEAP_SECONDS[i][6]);
break; break;
} }
} }
return utc_time; return utc_time;
} }
double Rinex_Printer::get_leap_second(const Glonass_Gnav_Ephemeris& eph, const double gps_obs_time) double Rinex_Printer::get_leap_second(const Glonass_Gnav_Ephemeris& eph, const double gps_obs_time)
{ {
double tod = 0.0; double tod = 0.0;
double glot2utc = 3*3600; double glot2utc = 3 * 3600;
double obs_time_glot = 0.0; double obs_time_glot = 0.0;
int i = 0; int i = 0;
double leap_second = 0; double leap_second = 0;
// Get observation time in nearly GLONASS time. Correction for leap seconds done at the end // Get observation time in nearly GLONASS time. Correction for leap seconds done at the end
obs_time_glot = gps_obs_time + glot2utc; obs_time_glot = gps_obs_time + glot2utc;
// Get seconds of day in glonass time // Get seconds of day in glonass time
tod = fmod (obs_time_glot, 86400); tod = fmod(obs_time_glot, 86400);
// Form date and time duration types // Form date and time duration types
boost::posix_time::time_duration t1(0, 0, tod); boost::posix_time::time_duration t1(0, 0, tod);
boost::gregorian::date d1(eph.d_yr, 1, 1); boost::gregorian::date d1(eph.d_yr, 1, 1);
boost::gregorian::days d2(eph.d_N_T-1); boost::gregorian::days d2(eph.d_N_T - 1);
boost::posix_time::ptime glo_time(d1 + d2, t1); boost::posix_time::ptime glo_time(d1 + d2, t1);
// Convert to utc // Convert to utc
boost::posix_time::time_duration t2(0, 0, glot2utc); boost::posix_time::time_duration t2(0, 0, glot2utc);
boost::posix_time::ptime utc_time = glo_time - t2; boost::posix_time::ptime utc_time = glo_time - t2;
// Adjust for leap second correction // Adjust for leap second correction
for (i = 0; GLONASS_LEAP_SECONDS[i][0]>0; i++) for (i = 0; GLONASS_LEAP_SECONDS[i][0] > 0; i++)
{ {
boost::posix_time::time_duration t3(GLONASS_LEAP_SECONDS[i][3], GLONASS_LEAP_SECONDS[i][4], GLONASS_LEAP_SECONDS[i][5]); boost::posix_time::time_duration t3(GLONASS_LEAP_SECONDS[i][3], GLONASS_LEAP_SECONDS[i][4], GLONASS_LEAP_SECONDS[i][5]);
boost::gregorian::date d3(GLONASS_LEAP_SECONDS[i][0], GLONASS_LEAP_SECONDS[i][1], GLONASS_LEAP_SECONDS[i][2]); boost::gregorian::date d3(GLONASS_LEAP_SECONDS[i][0], GLONASS_LEAP_SECONDS[i][1], GLONASS_LEAP_SECONDS[i][2]);
boost::posix_time::ptime ls_time(d3, t3); boost::posix_time::ptime ls_time(d3, t3);
if (utc_time >= ls_time) if (utc_time >= ls_time)
{ {
// We subtract the leap second when going from gpst to utc // We subtract the leap second when going from gpst to utc
leap_second = fabs(GLONASS_LEAP_SECONDS[i][6]); leap_second = fabs(GLONASS_LEAP_SECONDS[i][6]);
break; break;
} }
} }
return leap_second; return leap_second;
} }