mirrors/gnss-sdr
1
0
Fork 0
mirror of https://github.com/gnss-sdr/gnss-sdr synced 2025-02-05 21:50:09 +00:00
Code Issues Releases Wiki Activity

Update Access18 experiment with new PVT stored data structure

Browse Source
This commit is contained in:
Carles Fernandez 2018-08-31 10:03:35 +02:00
parent 34382d7ac5
commit 3fcb027138
No known key found for this signature in database
GPG Key ID: 4C583C52B0C3877D
3 changed files with 35 additions and 15 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/algorithms/PVT/libs/rtklib_solver.cc
View File

@ -599,8 +599,8 @@ bool rtklib_solver::get_PVT(const std::map<int, Gnss_Synchro>& gnss_observables_
tmp_double = pvt_sol.thres; tmp_double = pvt_sol.thres;
d_dump_file.write(reinterpret_cast<char*>(&pvt_sol.thres), sizeof(float)); d_dump_file.write(reinterpret_cast<char*>(&pvt_sol.thres), sizeof(float));
//GDOP//PDOP//HDOP//VDOP // GDOP / PDOP/ HDOP/ VDOP
d_dump_file.write(reinterpret_cast<char*>(&dop_[0]), sizeof(dop_)); d_dump_file.write(reinterpret_cast<char*>(&dop_[0]), sizeof(double) * 4);
} }
catch (const std::ifstream::failure& e) catch (const std::ifstream::failure& e)
{ {

2
src/utils/reproducibility/ieee-access18/L2-access18.conf
View File

@ -133,7 +133,7 @@ PVT.display_rate_ms=500
;# KML, GeoJSON, NMEA and RTCM output configuration ;# KML, GeoJSON, NMEA and RTCM output configuration
;#dump_filename: Log path and filename without extension. Notice that PVT will add ".dat" to the binary dump and ".kml" to GoogleEarth dump. ;#dump_filename: Log path and filename without extension. Notice that PVT will add ".dat" to the binary dump and ".kml" to GoogleEarth dump.
PVT.dump_filename=./data/PVT PVT.dump_filename=./data/access18
;#nmea_dump_filename: NMEA log path and filename ;#nmea_dump_filename: NMEA log path and filename
PVT.nmea_dump_filename=./gnss_sdr_pvt.nmea PVT.nmea_dump_filename=./gnss_sdr_pvt.nmea

44
src/utils/reproducibility/ieee-access18/plot_dump.m
View File

@ -62,26 +62,46 @@ ylabel('Navigation data bits','fontname','Times','fontsize', fontsize)
grid on grid on
fileID = fopen('data/PVT_ls_pvt.dat', 'r'); fileID = fopen('data/access18_pvt.dat', 'r');
dinfo = dir('data/PVT_ls_pvt.dat'); dinfo = dir('data/access18_pvt.dat');
filesize = dinfo.bytes; filesize = dinfo.bytes;
aux = 1; aux = 1;
while ne(ftell(fileID), filesize) while ne(ftell(fileID), filesize)
navsol.RX_time(aux) = fread(fileID, 1, 'double'); navsol.TOW_at_current_symbol_ms(aux) = fread(fileID, 1, 'uint32');
navsol.X(aux) = fread(fileID, 1, 'double'); navsol.week(aux) = fread(fileID, 1, 'uint32');
navsol.Y(aux) = fread(fileID, 1, 'double'); navsol.RX_time(aux) = fread(fileID, 1, 'double');
navsol.Z(aux) = fread(fileID, 1, 'double'); navsol.user_clock_offset(aux) = fread(fileID, 1, 'double');
navsol.user_clock(aux) = fread(fileID, 1, 'double'); navsol.X(aux) = fread(fileID, 1, 'double');
navsol.lat(aux) = fread(fileID, 1, 'double'); navsol.Y(aux) = fread(fileID, 1, 'double');
navsol.long(aux) = fread(fileID, 1, 'double'); navsol.Z(aux) = fread(fileID, 1, 'double');
navsol.height(aux) = fread(fileID, 1, 'double'); navsol.VX(aux) = fread(fileID, 1, 'double');
navsol.VY(aux) = fread(fileID, 1, 'double');
navsol.VZ(aux) = fread(fileID, 1, 'double');
navsol.varXX(aux) = fread(fileID, 1, 'double');
navsol.varYY(aux) = fread(fileID, 1, 'double');
navsol.varZZ(aux) = fread(fileID, 1, 'double');
navsol.varXY(aux) = fread(fileID, 1, 'double');
navsol.varYZ(aux) = fread(fileID, 1, 'double');
navsol.varZX(aux) = fread(fileID, 1, 'double');
navsol.latitude(aux) = fread(fileID, 1, 'double');
navsol.longitude(aux) = fread(fileID, 1, 'double');
navsol.height(aux) = fread(fileID, 1, 'double');
navsol.number_sats(aux) = fread(fileID, 1, 'uint8');
navsol.solution_status(aux) = fread(fileID, 1, 'uint8');
navsol.solution_type(aux) = fread(fileID, 1, 'uint8');
navsol.AR_ratio_factor(aux) = fread(fileID, 1, 'float');
navsol.AR_ratio_threshold(aux) = fread(fileID, 1, 'float');
navsol.GDOP(aux) = fread(fileID, 1, 'double');
navsol.PDOP(aux) = fread(fileID, 1, 'double');
navsol.HDOP(aux) = fread(fileID, 1, 'double');
navsol.VDOP(aux) = fread(fileID, 1, 'double');
aux = aux + 1; aux = aux + 1;
end end
fclose(fileID); fclose(fileID);
mean_Latitude = mean(navsol.lat); mean_Latitude = mean(navsol.latitude);
mean_Longitude = mean(navsol.long); mean_Longitude = mean(navsol.longitude);
mean_h = mean(navsol.height); mean_h = mean(navsol.height);
utmZone = findUtmZone(mean_Latitude, mean_Longitude); utmZone = findUtmZone(mean_Latitude, mean_Longitude);
[ref_X_cart, ref_Y_cart, ref_Z_cart] = geo2cart(dms2mat(deg2dms(mean_Latitude)), dms2mat(deg2dms(mean_Longitude)), mean_h, 5); [ref_X_cart, ref_Y_cart, ref_Z_cart] = geo2cart(dms2mat(deg2dms(mean_Latitude)), dms2mat(deg2dms(mean_Longitude)), mean_h, 5);