source: ntrip/trunk/BNC/src/rinex/reqcedit.cpp@ 6809

Last change on this file since 6809 was 6809, checked in by stuerze, 9 years ago

I/NAV - F/NAV issue (hopefully:) solved in another way
File size: 18.2 KB
Line 
1// Part of BNC, a utility for retrieving decoding and
2// converting GNSS data streams from NTRIP broadcasters.
3//
4// Copyright (C) 2007
5// German Federal Agency for Cartography and Geodesy (BKG)
6// http://www.bkg.bund.de
7// Czech Technical University Prague, Department of Geodesy
8// http://www.fsv.cvut.cz
9//
10// Email: euref-ip@bkg.bund.de
11//
12// This program is free software; you can redistribute it and/or
13// modify it under the terms of the GNU General Public License
14// as published by the Free Software Foundation, version 2.
15//
16// This program is distributed in the hope that it will be useful,
17// but WITHOUT ANY WARRANTY; without even the implied warranty of
18// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19// GNU General Public License for more details.
20//
21// You should have received a copy of the GNU General Public License
22// along with this program; if not, write to the Free Software
23// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
24
25/* -------------------------------------------------------------------------
26 * BKG NTRIP Client
27 * -------------------------------------------------------------------------
28 *
29 * Class: t_reqcEdit
30 *
31 * Purpose: Edit/Concatenate RINEX Files
32 *
33 * Author: L. Mervart
34 *
35 * Created: 11-Apr-2012
36 *
37 * Changes:
38 *
39 * -----------------------------------------------------------------------*/
40
41#include <iostream>
42#include "reqcedit.h"
43#include "bnccore.h"
44#include "bncsettings.h"
45#include "bncutils.h"
46#include "rnxobsfile.h"
47#include "rnxnavfile.h"
48
49using namespace std;
50
51// Constructor
52////////////////////////////////////////////////////////////////////////////
53t_reqcEdit::t_reqcEdit(QObject* parent) : QThread(parent) {
54
55 bncSettings settings;
56
57 _logFileName = settings.value("reqcOutLogFile").toString(); expandEnvVar(_logFileName);
58 _logFile = 0;
59 _log = 0;
60 _obsFileNames = settings.value("reqcObsFile").toString().split(",", QString::SkipEmptyParts);
61 _outObsFileName = settings.value("reqcOutObsFile").toString();
62 _navFileNames = settings.value("reqcNavFile").toString().split(",", QString::SkipEmptyParts);
63 _outNavFileName = settings.value("reqcOutNavFile").toString();
64 int version = settings.value("reqcRnxVersion").toInt();
65 if (version < 3) {
66 _rnxVersion = t_rnxObsHeader::defaultRnxObsVersion2;
67 }
68 else {
69 _rnxVersion = t_rnxObsHeader::defaultRnxObsVersion3;
70 }
71 _samplingRate = settings.value("reqcSampling").toInt();
72 _begTime = bncTime(settings.value("reqcStartDateTime").toString().toAscii().data());
73 _endTime = bncTime(settings.value("reqcEndDateTime").toString().toAscii().data());
74}
75
76// Destructor
77////////////////////////////////////////////////////////////////////////////
78t_reqcEdit::~t_reqcEdit() {
79 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
80 delete _rnxObsFiles[ii];
81 }
82 for (int ii = 0; ii < _ephs.size(); ii++) {
83 delete _ephs[ii];
84 }
85 delete _log; _log = 0;
86 delete _logFile; _logFile = 0;
87}
88
89//
90////////////////////////////////////////////////////////////////////////////
91void t_reqcEdit::run() {
92
93 // Open Log File
94 // -------------
95 _logFile = new QFile(_logFileName);
96 if (_logFile->open(QIODevice::WriteOnly | QIODevice::Text)) {
97 _log = new QTextStream();
98 _log->setDevice(_logFile);
99 }
100
101 // Log File Header
102 // ---------------
103 if (_log) {
104 *_log << QByteArray(78, '-') << endl;
105 *_log << "Concatenation of RINEX Observation and/or Navigation Files\n";
106 *_log << QByteArray(78, '-') << endl;
107
108 *_log << QByteArray("Program").leftJustified(15) << ": "
109 << BNC_CORE->pgmName() << endl;
110 *_log << QByteArray("Run by").leftJustified(15) << ": "
111 << BNC_CORE->userName() << endl;
112 *_log << QByteArray("Date").leftJustified(15) << ": "
113 << QDateTime::currentDateTime().toUTC().toString("yyyy-MM-dd hh:mm:ss") << endl;
114 *_log << QByteArray("RINEX Version").leftJustified(15) << ": "
115 << _rnxVersion << endl;
116 *_log << QByteArray("Sampling").leftJustified(15) << ": "
117 << _samplingRate << endl;
118 *_log << QByteArray("Start time").leftJustified(15) << ": "
119 << _begTime.datestr().c_str() << ' '
120 << _begTime.timestr(0).c_str() << endl;
121 *_log << QByteArray("End time").leftJustified(15) << ": "
122 << _endTime.datestr().c_str() << ' '
123 << _endTime.timestr(0).c_str() << endl;
124 *_log << QByteArray("Input Obs Files").leftJustified(15) << ": "
125 << _obsFileNames.join(",") << endl;
126 *_log << QByteArray("Input Nav Files").leftJustified(15) << ": "
127 << _navFileNames.join(",") << endl;
128 *_log << QByteArray("Output Obs File").leftJustified(15) << ": "
129 << _outObsFileName << endl;
130 *_log << QByteArray("Output Nav File").leftJustified(15) << ": "
131 << _outNavFileName << endl;
132
133 *_log << QByteArray(78, '-') << endl;
134 _log->flush();
135 }
136
137 // Handle Observation Files
138 // ------------------------
139 editObservations();
140
141 // Handle Navigations Files
142 // ------------------------
143 editEphemerides();
144
145 // Exit (thread)
146 // -------------
147 if (BNC_CORE->mode() != t_bncCore::interactive) {
148 qApp->exit(0);
149 }
150 else {
151 emit finished();
152 deleteLater();
153 }
154}
155
156// Initialize input observation files, sort them according to start time
157////////////////////////////////////////////////////////////////////////////
158void t_reqcEdit::initRnxObsFiles(const QStringList& obsFileNames,
159 QVector<t_rnxObsFile*>& rnxObsFiles,
160 QTextStream* log) {
161
162 QStringListIterator it(obsFileNames);
163 while (it.hasNext()) {
164 QString fileName = it.next();
165 if (fileName.indexOf('*') != -1 || fileName.indexOf('?') != -1) {
166 QFileInfo fileInfo(fileName);
167 QDir dir = fileInfo.dir();
168 QStringList filters; filters << fileInfo.fileName();
169 QListIterator<QFileInfo> it(dir.entryInfoList(filters));
170 while (it.hasNext()) {
171 QString filePath = it.next().filePath();
172 t_rnxObsFile* rnxObsFile = 0;
173 try {
174 rnxObsFile = new t_rnxObsFile(filePath, t_rnxObsFile::input);
175 rnxObsFiles.append(rnxObsFile);
176 }
177 catch (...) {
178 delete rnxObsFile;
179 if (log) {
180 *log << "Error in rnxObsFile " << filePath.toAscii().data() << endl;
181 }
182 }
183 }
184 }
185 else {
186 t_rnxObsFile* rnxObsFile = 0;
187 try {
188 rnxObsFile = new t_rnxObsFile(fileName, t_rnxObsFile::input);
189 rnxObsFiles.append(rnxObsFile);
190 }
191 catch (...) {
192 if (log) {
193 *log << "Error in rnxObsFile " << fileName.toAscii().data() << endl;
194 }
195 }
196 }
197 }
198 qStableSort(rnxObsFiles.begin(), rnxObsFiles.end(),
199 t_rnxObsFile::earlierStartTime);
200}
201
202//
203////////////////////////////////////////////////////////////////////////////
204void t_reqcEdit::editObservations() {
205
206 // Easy Exit
207 // ---------
208 if (_obsFileNames.isEmpty() || _outObsFileName.isEmpty()) {
209 return;
210 }
211
212 t_reqcEdit::initRnxObsFiles(_obsFileNames, _rnxObsFiles, _log);
213
214 // Initialize output observation file
215 // ----------------------------------
216 t_rnxObsFile outObsFile(_outObsFileName, t_rnxObsFile::output);
217
218 // Select observation types
219 // ------------------------
220 bncSettings settings;
221 QStringList useObsTypes = settings.value("reqcUseObsTypes").toString().split(" ", QString::SkipEmptyParts);
222
223 // Put together all observation types
224 // ----------------------------------
225 if (_rnxObsFiles.size() > 1 && useObsTypes.size() == 0) {
226 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
227 t_rnxObsFile* obsFile = _rnxObsFiles[ii];
228 for (int iSys = 0; iSys < obsFile->numSys(); iSys++) {
229 char sys = obsFile->system(iSys);
230 if (sys != ' ') {
231 for (int iType = 0; iType < obsFile->nTypes(sys); iType++) {
232 QString type = obsFile->obsType(sys, iType);
233 if (_rnxVersion < 3.0) {
234 useObsTypes << type;
235 }
236 else {
237 useObsTypes << QString(sys) + ":" + type;
238 }
239 }
240 }
241 }
242 }
243 useObsTypes.removeDuplicates();
244 }
245
246 // Loop over all input observation files
247 // -------------------------------------
248 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
249 t_rnxObsFile* obsFile = _rnxObsFiles[ii];
250 if (_log) {
251 *_log << "Processing File: " << obsFile->fileName() << " start: "
252 << obsFile->startTime().datestr().c_str() << ' '
253 << obsFile->startTime().timestr(0).c_str() << endl;
254 }
255 if (ii == 0) {
256 outObsFile.setHeader(obsFile->header(), int(_rnxVersion), &useObsTypes);
257 if (_begTime.valid() && _begTime > outObsFile.startTime()) {
258 outObsFile.setStartTime(_begTime);
259 }
260 if (_samplingRate > outObsFile.interval()) {
261 outObsFile.setInterval(_samplingRate);
262 }
263 editRnxObsHeader(outObsFile);
264 bncSettings settings;
265 QMap<QString, QString> txtMap;
266 QString runBy = settings.value("reqcRunBy").toString();
267 if (!runBy.isEmpty()) {
268 txtMap["RUN BY"] = runBy;
269 }
270 QString comment = settings.value("reqcComment").toString();
271 if (!comment.isEmpty()) {
272 txtMap["COMMENT"] = comment;
273 }
274 outObsFile.header().write(outObsFile.stream(), &txtMap);
275 }
276 t_rnxObsFile::t_rnxEpo* epo = 0;
277 try {
278 while ( (epo = obsFile->nextEpoch()) != 0) {
279 if (_begTime.valid() && epo->tt < _begTime) {
280 continue;
281 }
282 if (_endTime.valid() && epo->tt > _endTime) {
283 break;
284 }
285
286 if (_samplingRate == 0 ||
287 fmod(round(epo->tt.gpssec()), _samplingRate) == 0) {
288 applyLLI(obsFile, epo);
289 outObsFile.writeEpoch(epo);
290 }
291 else {
292 rememberLLI(obsFile, epo);
293 }
294 }
295 }
296 catch (QString str) {
297 if (_log) {
298 *_log << "Exception " << str << endl;
299 }
300 else {
301 qDebug() << str;
302 }
303 return;
304 }
305 catch (...) {
306 if (_log) {
307 *_log << "Exception unknown" << endl;
308 }
309 else {
310 qDebug() << "Exception unknown";
311 }
312 return;
313 }
314 }
315}
316
317// Change RINEX Header Content
318////////////////////////////////////////////////////////////////////////////
319void t_reqcEdit::editRnxObsHeader(t_rnxObsFile& obsFile) {
320
321 bncSettings settings;
322
323 QString oldMarkerName = settings.value("reqcOldMarkerName").toString();
324 QString newMarkerName = settings.value("reqcNewMarkerName").toString();
325 if (!newMarkerName.isEmpty()) {
326 if (oldMarkerName.isEmpty() ||
327 QRegExp(oldMarkerName).exactMatch(obsFile.markerName())) {
328 obsFile.setMarkerName(newMarkerName);
329 }
330 }
331
332 QString oldAntennaName = settings.value("reqcOldAntennaName").toString();
333 QString newAntennaName = settings.value("reqcNewAntennaName").toString();
334 if (!newAntennaName.isEmpty()) {
335 if (oldAntennaName.isEmpty() ||
336 QRegExp(oldAntennaName).exactMatch(obsFile.antennaName())) {
337 obsFile.setAntennaName(newAntennaName);
338 }
339 }
340
341 QString oldAntennaNumber = settings.value("reqcOldAntennaNumber").toString();
342 QString newAntennaNumber = settings.value("reqcNewAntennaNumber").toString();
343 if (!newAntennaNumber.isEmpty()) {
344 if (oldAntennaNumber.isEmpty() ||
345 QRegExp(oldAntennaNumber).exactMatch(obsFile.antennaNumber())) {
346 obsFile.setAntennaNumber(newAntennaNumber);
347 }
348 }
349
350
351 const ColumnVector& obsFileAntNEU = obsFile.antNEU();
352 QString oldAntennadN = settings.value("reqcOldAntennadN").toString();
353 QString newAntennadN = settings.value("reqcNewAntennadN").toString();
354 if(!newAntennadN.isEmpty()) {
355 if (oldAntennadN.isEmpty() ||
356 oldAntennadN.toDouble() == obsFileAntNEU(1)) {
357 obsFile.setAntennaN(newAntennadN.toDouble());
358 }
359 }
360 QString oldAntennadE = settings.value("reqcOldAntennadE").toString();
361 QString newAntennadE = settings.value("reqcNewAntennadE").toString();
362 if(!newAntennadE.isEmpty()) {
363 if (oldAntennadE.isEmpty() ||
364 oldAntennadE.toDouble() == obsFileAntNEU(2)) {
365 obsFile.setAntennaE(newAntennadE.toDouble());
366 }
367 }
368 QString oldAntennadU = settings.value("reqcOldAntennadU").toString();
369 QString newAntennadU = settings.value("reqcNewAntennadU").toString();
370 if(!newAntennadU.isEmpty()) {
371 if (oldAntennadU.isEmpty() ||
372 oldAntennadU.toDouble() == obsFileAntNEU(3)) {
373 obsFile.setAntennaU(newAntennadU.toDouble());
374 }
375 }
376
377 QString oldReceiverType = settings.value("reqcOldReceiverName").toString();
378 QString newReceiverType = settings.value("reqcNewReceiverName").toString();
379 if (!newReceiverType.isEmpty()) {
380 if (oldReceiverType.isEmpty() ||
381 QRegExp(oldReceiverType).exactMatch(obsFile.receiverType())) {
382 obsFile.setReceiverType(newReceiverType);
383 }
384 }
385
386 QString oldReceiverNumber = settings.value("reqcOldReceiverNumber").toString();
387 QString newReceiverNumber = settings.value("reqcNewReceiverNumber").toString();
388 if (!newReceiverNumber.isEmpty()) {
389 if (oldReceiverNumber.isEmpty() ||
390 QRegExp(oldReceiverNumber).exactMatch(obsFile.receiverNumber())) {
391 obsFile.setReceiverNumber(newReceiverNumber);
392 }
393 }
394}
395
396//
397////////////////////////////////////////////////////////////////////////////
398void t_reqcEdit::rememberLLI(const t_rnxObsFile* obsFile,
399 const t_rnxObsFile::t_rnxEpo* epo) {
400
401 if (_samplingRate == 0) {
402 return;
403 }
404
405 for (unsigned iSat = 0; iSat < epo->rnxSat.size(); iSat++) {
406 const t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iSat];
407 char sys = rnxSat.prn.system();
408 QString prn(rnxSat.prn.toString().c_str());
409
410 for (int iType = 0; iType < obsFile->nTypes(sys); iType++) {
411 QString type = obsFile->obsType(sys, iType);
412 if (!_lli[prn].contains(iType)) {
413 _lli[prn][iType] = 0;
414 }
415 if (rnxSat.obs.contains(type) && rnxSat.obs[type].lli & 1) {
416 _lli[prn][iType] |= 1;
417 }
418 }
419 }
420}
421
422//
423////////////////////////////////////////////////////////////////////////////
424void t_reqcEdit::applyLLI(const t_rnxObsFile* obsFile,
425 t_rnxObsFile::t_rnxEpo* epo) {
426
427 if (_samplingRate == 0) {
428 return;
429 }
430
431 for (unsigned iSat = 0; iSat < epo->rnxSat.size(); iSat++) {
432 t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iSat];
433 char sys = rnxSat.prn.system();
434 QString prn(rnxSat.prn.toString().c_str());
435
436 for (int iType = 0; iType < obsFile->nTypes(sys); iType++) {
437 QString type = obsFile->obsType(sys, iType);
438 if (_lli[prn].contains(iType) && _lli[prn][iType] & 1) {
439 if (rnxSat.obs.contains(type)) {
440 rnxSat.obs[type].lli |= 1;
441 }
442 }
443 }
444 }
445
446 _lli.clear();
447}
448
449/// Read All Ephemerides
450////////////////////////////////////////////////////////////////////////////
451void t_reqcEdit::readEphemerides(const QStringList& navFileNames,
452 QVector<t_eph*>& ephs) {
453
454 QStringListIterator it(navFileNames);
455 while (it.hasNext()) {
456 QString fileName = it.next();
457 if (fileName.indexOf('*') != -1 || fileName.indexOf('?') != -1) {
458 QFileInfo fileInfo(fileName);
459 QDir dir = fileInfo.dir();
460 QStringList filters; filters << fileInfo.fileName();
461 QListIterator<QFileInfo> it(dir.entryInfoList(filters));
462 while (it.hasNext()) {
463 QString filePath = it.next().filePath();
464 appendEphemerides(filePath, ephs);
465 }
466 }
467 else {
468 appendEphemerides(fileName, ephs);
469 }
470 }
471 qStableSort(ephs.begin(), ephs.end(), t_eph::earlierTime);
472}
473
474//
475////////////////////////////////////////////////////////////////////////////
476void t_reqcEdit::editEphemerides() {
477
478 // Easy Exit
479 // ---------
480 if (_navFileNames.isEmpty() || _outNavFileName.isEmpty()) {
481 return;
482 }
483
484 // Read Ephemerides
485 // ----------------
486 t_reqcEdit::readEphemerides(_navFileNames, _ephs);
487
488 // Check Satellite Systems
489 // -----------------------
490 bool haveGPS = false;
491 bool haveGlonass = false;
492 for (int ii = 0; ii < _ephs.size(); ii++) {
493 const t_eph* eph = _ephs[ii];
494 if (eph->type() == t_eph::GPS) {
495 haveGPS = true;
496 }
497 else if (eph->type() == t_eph::GLONASS) {
498 haveGlonass = true;
499 }
500 }
501
502 // Initialize output navigation file
503 // ---------------------------------
504 t_rnxNavFile outNavFile(_outNavFileName, t_rnxNavFile::output);
505
506 outNavFile.setGlonass(haveGlonass);
507
508 if ( (haveGPS && haveGlonass) || _rnxVersion >= 3.0) {
509 outNavFile.setVersion(t_rnxNavFile::defaultRnxNavVersion3);
510 }
511 else {
512 outNavFile.setVersion(t_rnxNavFile::defaultRnxNavVersion2);
513 }
514
515 bncSettings settings;
516 QMap<QString, QString> txtMap;
517 QString runBy = settings.value("reqcRunBy").toString();
518 if (!runBy.isEmpty()) {
519 txtMap["RUN BY"] = runBy;
520 }
521 QString comment = settings.value("reqcComment").toString();
522 if (!comment.isEmpty()) {
523 txtMap["COMMENT"] = comment;
524 }
525
526 outNavFile.writeHeader(&txtMap);
527
528 // Loop over all ephemerides
529 // -------------------------
530 for (int ii = 0; ii < _ephs.size(); ii++) {
531 const t_eph* eph = _ephs[ii];
532 outNavFile.writeEph(eph);
533 }
534}
535
536//
537////////////////////////////////////////////////////////////////////////////
538void t_reqcEdit::appendEphemerides(const QString& fileName,
539 QVector<t_eph*>& ephs) {
540
541 t_rnxNavFile rnxNavFile(fileName, t_rnxNavFile::input);
542 for (unsigned ii = 0; ii < rnxNavFile.ephs().size(); ii++) {
543 t_eph* eph = rnxNavFile.ephs()[ii];
544 bool isNew = true;
545 for (int iOld = 0; iOld < ephs.size(); iOld++) {
546 const t_eph* ephOld = ephs[iOld];
547 if (ephOld->prn() == eph->prn() &&
548 ephOld->TOC() == eph->TOC()) {
549 isNew = false;
550 break;
551 }
552 }
553 if (isNew) {
554 if (eph->type() == t_eph::GPS) {
555 ephs.append(new t_ephGPS(*dynamic_cast<t_ephGPS*>(eph)));
556 }
557 else if (eph->type() == t_eph::GLONASS) {
558 ephs.append(new t_ephGlo(*dynamic_cast<t_ephGlo*>(eph)));
559 }
560 else if (eph->type() == t_eph::Galileo) {
561 ephs.append(new t_ephGal(*dynamic_cast<t_ephGal*>(eph)));
562 }
563 else if (eph->type() == t_eph::QZSS) {
564 ephs.append(new t_ephGPS(*dynamic_cast<t_ephGPS*>(eph)));
565 }
566 else if (eph->type() == t_eph::SBAS) {
567 ephs.append(new t_ephSBAS(*dynamic_cast<t_ephSBAS*>(eph)));
568 }
569 else if (eph->type() == t_eph::BDS) {
570 ephs.append(new t_ephBDS(*dynamic_cast<t_ephBDS*>(eph)));
571 }
572 }
573 }
574}
Note: See TracBrowser for help on using the repository browser.