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

Last change on this file since 7280 was 7280, checked in by weber, 9 years ago

Avoid error message through opening log file only if file name is specified
File size: 19.9 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
77// Destructor
78////////////////////////////////////////////////////////////////////////////
79t_reqcEdit::~t_reqcEdit() {
80 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
81 delete _rnxObsFiles[ii];
82 }
83 for (int ii = 0; ii < _ephs.size(); ii++) {
84 delete _ephs[ii];
85 }
86 delete _log; _log = 0;
87 delete _logFile; _logFile = 0;
88}
89
90//
91////////////////////////////////////////////////////////////////////////////
92void t_reqcEdit::run() {
93
94 // Open Log File
95 // -------------
96 if (!_logFileName.isEmpty()) {
97 _logFile = new QFile(_logFileName);
98 if (_logFile->open(QIODevice::WriteOnly | QIODevice::Text)) {
99 _log = new QTextStream();
100 _log->setDevice(_logFile);
101 }
102 }
103
104 // Log File Header
105 // ---------------
106 if (_log) {
107 *_log << QByteArray(78, '-') << endl;
108 *_log << "Concatenation of RINEX Observation and/or Navigation Files\n";
109 *_log << QByteArray(78, '-') << endl;
110
111 *_log << QByteArray("Program").leftJustified(15) << ": "
112 << BNC_CORE->pgmName() << endl;
113 *_log << QByteArray("Run by").leftJustified(15) << ": "
114 << BNC_CORE->userName() << endl;
115 *_log << QByteArray("Date").leftJustified(15) << ": "
116 << QDateTime::currentDateTime().toUTC().toString("yyyy-MM-dd hh:mm:ss") << endl;
117 *_log << QByteArray("RINEX Version").leftJustified(15) << ": "
118 << _rnxVersion << endl;
119 *_log << QByteArray("Sampling").leftJustified(15) << ": "
120 << _samplingRate << endl;
121 *_log << QByteArray("Start time").leftJustified(15) << ": "
122 << _begTime.datestr().c_str() << ' '
123 << _begTime.timestr(0).c_str() << endl;
124 *_log << QByteArray("End time").leftJustified(15) << ": "
125 << _endTime.datestr().c_str() << ' '
126 << _endTime.timestr(0).c_str() << endl;
127 *_log << QByteArray("Input Obs Files").leftJustified(15) << ": "
128 << _obsFileNames.join(",") << endl;
129 *_log << QByteArray("Input Nav Files").leftJustified(15) << ": "
130 << _navFileNames.join(",") << endl;
131 *_log << QByteArray("Output Obs File").leftJustified(15) << ": "
132 << _outObsFileName << endl;
133 *_log << QByteArray("Output Nav File").leftJustified(15) << ": "
134 << _outNavFileName << endl;
135
136 *_log << QByteArray(78, '-') << endl;
137 _log->flush();
138 }
139
140 // Handle Observation Files
141 // ------------------------
142 editObservations();
143
144 // Handle Navigations Files
145 // ------------------------
146 editEphemerides();
147
148 // Exit (thread)
149 // -------------
150 if (BNC_CORE->mode() != t_bncCore::interactive) {
151 qApp->exit(0);
152 }
153 else {
154 emit finished();
155 deleteLater();
156 }
157}
158
159// Initialize input observation files, sort them according to start time
160////////////////////////////////////////////////////////////////////////////
161void t_reqcEdit::initRnxObsFiles(const QStringList& obsFileNames,
162 QVector<t_rnxObsFile*>& rnxObsFiles,
163 QTextStream* log) {
164
165 QStringListIterator it(obsFileNames);
166 while (it.hasNext()) {
167 QString fileName = it.next();
168 if (fileName.indexOf('*') != -1 || fileName.indexOf('?') != -1) {
169 QFileInfo fileInfo(fileName);
170 QDir dir = fileInfo.dir();
171 QStringList filters; filters << fileInfo.fileName();
172 QListIterator<QFileInfo> it(dir.entryInfoList(filters));
173 while (it.hasNext()) {
174 QString filePath = it.next().filePath();
175 t_rnxObsFile* rnxObsFile = 0;
176 try {
177 rnxObsFile = new t_rnxObsFile(filePath, t_rnxObsFile::input);
178 rnxObsFiles.append(rnxObsFile);
179 }
180 catch (...) {
181 delete rnxObsFile;
182 if (log) {
183 *log << "Error in rnxObsFile " << filePath.toAscii().data() << endl;
184 }
185 }
186 }
187 }
188 else {
189 t_rnxObsFile* rnxObsFile = 0;
190 try {
191 rnxObsFile = new t_rnxObsFile(fileName, t_rnxObsFile::input);
192 rnxObsFiles.append(rnxObsFile);
193 }
194 catch (...) {
195 if (log) {
196 *log << "Error in rnxObsFile " << fileName.toAscii().data() << endl;
197 }
198 }
199 }
200 }
201 qStableSort(rnxObsFiles.begin(), rnxObsFiles.end(),
202 t_rnxObsFile::earlierStartTime);
203}
204
205//
206////////////////////////////////////////////////////////////////////////////
207void t_reqcEdit::editObservations() {
208
209 // Easy Exit
210 // ---------
211 if (_obsFileNames.isEmpty() || _outObsFileName.isEmpty()) {
212 return;
213 }
214
215 t_reqcEdit::initRnxObsFiles(_obsFileNames, _rnxObsFiles, _log);
216
217 // Initialize output observation file
218 // ----------------------------------
219 t_rnxObsFile outObsFile(_outObsFileName, t_rnxObsFile::output);
220
221 // Select observation types
222 // ------------------------
223 bncSettings settings;
224 QStringList useObsTypes = settings.value("reqcUseObsTypes").toString().split(" ", QString::SkipEmptyParts);
225
226 // Put together all observation types
227 // ----------------------------------
228 if (_rnxObsFiles.size() > 1 && useObsTypes.size() == 0) {
229 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
230 t_rnxObsFile* obsFile = _rnxObsFiles[ii];
231 for (int iSys = 0; iSys < obsFile->numSys(); iSys++) {
232 char sys = obsFile->system(iSys);
233 if (sys != ' ') {
234 for (int iType = 0; iType < obsFile->nTypes(sys); iType++) {
235 QString type = obsFile->obsType(sys, iType);
236 if (_rnxVersion < 3.0) {
237 useObsTypes << type;
238 }
239 else {
240 useObsTypes << QString(sys) + ":" + type;
241 }
242 }
243 }
244 }
245 }
246 useObsTypes.removeDuplicates();
247 }
248
249 // Put together all phase shifts
250 // -----------------------------
251 QStringList phaseShifts;
252 if (_rnxVersion >= 3.0 && _rnxObsFiles.size() > 1) {
253 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
254 t_rnxObsFile* obsFile = _rnxObsFiles[ii];
255 phaseShifts << obsFile->phaseShifts();
256 }
257 phaseShifts.removeDuplicates();
258 }
259
260 // Put together all GLONASS biases
261 // -------------------------------
262 QStringList gloBiases;
263 if (_rnxVersion >= 3.0 && _rnxObsFiles.size() > 1) {
264 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
265 t_rnxObsFile* obsFile = _rnxObsFiles[ii];
266 if (ii == 0 && obsFile->numGloBiases() == 4) {
267 break;
268 }
269 else {
270 gloBiases << obsFile->gloBiases();
271 }
272 }
273 gloBiases.removeDuplicates();
274 }
275
276 // Put together all GLONASS slots
277 // -----------------------------
278 QStringList gloSlots;
279 if (_rnxVersion >= 3.0 && _rnxObsFiles.size() > 1) {
280 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
281 t_rnxObsFile* obsFile = _rnxObsFiles[ii];
282 if (ii == 0 &&
283 obsFile->numGloSlots() == signed(t_prn::MAXPRN_GLONASS)) {
284 break;
285 }
286 else {
287 gloSlots << obsFile->gloSlots();
288 }
289 }
290 gloSlots.removeDuplicates();
291 }
292
293 // Loop over all input observation files
294 // -------------------------------------
295 for (int ii = 0; ii < _rnxObsFiles.size(); ii++) {
296 t_rnxObsFile* obsFile = _rnxObsFiles[ii];
297 if (_log) {
298 *_log << "Processing File: " << obsFile->fileName() << " start: "
299 << obsFile->startTime().datestr().c_str() << ' '
300 << obsFile->startTime().timestr(0).c_str() << endl;
301 }
302 if (ii == 0) {
303 outObsFile.setHeader(obsFile->header(), int(_rnxVersion), &useObsTypes,
304 &phaseShifts, &gloBiases, &gloSlots);
305 if (_begTime.valid() && _begTime > outObsFile.startTime()) {
306 outObsFile.setStartTime(_begTime);
307 }
308 if (_samplingRate > outObsFile.interval()) {
309 outObsFile.setInterval(_samplingRate);
310 }
311 editRnxObsHeader(outObsFile);
312 bncSettings settings;
313 QMap<QString, QString> txtMap;
314 QString runBy = settings.value("reqcRunBy").toString();
315 if (!runBy.isEmpty()) {
316 txtMap["RUN BY"] = runBy;
317 }
318 QString comment = settings.value("reqcComment").toString();
319 if (!comment.isEmpty()) {
320 txtMap["COMMENT"] = comment;
321 }
322 outObsFile.header().write(outObsFile.stream(), &txtMap);
323 }
324 t_rnxObsFile::t_rnxEpo* epo = 0;
325 try {
326 while ( (epo = obsFile->nextEpoch()) != 0) {
327 if (_begTime.valid() && epo->tt < _begTime) {
328 continue;
329 }
330 if (_endTime.valid() && epo->tt > _endTime) {
331 break;
332 }
333
334 if (_samplingRate == 0 ||
335 fmod(round(epo->tt.gpssec()), _samplingRate) == 0) {
336 applyLLI(obsFile, epo);
337 outObsFile.writeEpoch(epo);
338 }
339 else {
340 rememberLLI(obsFile, epo);
341 }
342 }
343 }
344 catch (QString str) {
345 if (_log) {
346 *_log << "Exception " << str << endl;
347 }
348 else {
349 qDebug() << str;
350 }
351 return;
352 }
353 catch (...) {
354 if (_log) {
355 *_log << "Exception unknown" << endl;
356 }
357 else {
358 qDebug() << "Exception unknown";
359 }
360 return;
361 }
362 }
363}
364
365// Change RINEX Header Content
366////////////////////////////////////////////////////////////////////////////
367void t_reqcEdit::editRnxObsHeader(t_rnxObsFile& obsFile) {
368
369 bncSettings settings;
370
371 QString oldMarkerName = settings.value("reqcOldMarkerName").toString();
372 QString newMarkerName = settings.value("reqcNewMarkerName").toString();
373 if (!newMarkerName.isEmpty()) {
374 if (oldMarkerName.isEmpty() ||
375 QRegExp(oldMarkerName).exactMatch(obsFile.markerName())) {
376 obsFile.setMarkerName(newMarkerName);
377 }
378 }
379
380 QString oldAntennaName = settings.value("reqcOldAntennaName").toString();
381 QString newAntennaName = settings.value("reqcNewAntennaName").toString();
382 if (!newAntennaName.isEmpty()) {
383 if (oldAntennaName.isEmpty() ||
384 QRegExp(oldAntennaName).exactMatch(obsFile.antennaName())) {
385 obsFile.setAntennaName(newAntennaName);
386 }
387 }
388
389 QString oldAntennaNumber = settings.value("reqcOldAntennaNumber").toString();
390 QString newAntennaNumber = settings.value("reqcNewAntennaNumber").toString();
391 if (!newAntennaNumber.isEmpty()) {
392 if (oldAntennaNumber.isEmpty() ||
393 QRegExp(oldAntennaNumber).exactMatch(obsFile.antennaNumber())) {
394 obsFile.setAntennaNumber(newAntennaNumber);
395 }
396 }
397
398
399 const ColumnVector& obsFileAntNEU = obsFile.antNEU();
400 QString oldAntennadN = settings.value("reqcOldAntennadN").toString();
401 QString newAntennadN = settings.value("reqcNewAntennadN").toString();
402 if(!newAntennadN.isEmpty()) {
403 if (oldAntennadN.isEmpty() ||
404 oldAntennadN.toDouble() == obsFileAntNEU(1)) {
405 obsFile.setAntennaN(newAntennadN.toDouble());
406 }
407 }
408 QString oldAntennadE = settings.value("reqcOldAntennadE").toString();
409 QString newAntennadE = settings.value("reqcNewAntennadE").toString();
410 if(!newAntennadE.isEmpty()) {
411 if (oldAntennadE.isEmpty() ||
412 oldAntennadE.toDouble() == obsFileAntNEU(2)) {
413 obsFile.setAntennaE(newAntennadE.toDouble());
414 }
415 }
416 QString oldAntennadU = settings.value("reqcOldAntennadU").toString();
417 QString newAntennadU = settings.value("reqcNewAntennadU").toString();
418 if(!newAntennadU.isEmpty()) {
419 if (oldAntennadU.isEmpty() ||
420 oldAntennadU.toDouble() == obsFileAntNEU(3)) {
421 obsFile.setAntennaU(newAntennadU.toDouble());
422 }
423 }
424
425 QString oldReceiverType = settings.value("reqcOldReceiverName").toString();
426 QString newReceiverType = settings.value("reqcNewReceiverName").toString();
427 if (!newReceiverType.isEmpty()) {
428 if (oldReceiverType.isEmpty() ||
429 QRegExp(oldReceiverType).exactMatch(obsFile.receiverType())) {
430 obsFile.setReceiverType(newReceiverType);
431 }
432 }
433
434 QString oldReceiverNumber = settings.value("reqcOldReceiverNumber").toString();
435 QString newReceiverNumber = settings.value("reqcNewReceiverNumber").toString();
436 if (!newReceiverNumber.isEmpty()) {
437 if (oldReceiverNumber.isEmpty() ||
438 QRegExp(oldReceiverNumber).exactMatch(obsFile.receiverNumber())) {
439 obsFile.setReceiverNumber(newReceiverNumber);
440 }
441 }
442}
443
444//
445////////////////////////////////////////////////////////////////////////////
446void t_reqcEdit::rememberLLI(const t_rnxObsFile* obsFile,
447 const t_rnxObsFile::t_rnxEpo* epo) {
448
449 if (_samplingRate == 0) {
450 return;
451 }
452
453 for (unsigned iSat = 0; iSat < epo->rnxSat.size(); iSat++) {
454 const t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iSat];
455 char sys = rnxSat.prn.system();
456 QString prn(rnxSat.prn.toString().c_str());
457
458 for (int iType = 0; iType < obsFile->nTypes(sys); iType++) {
459 QString type = obsFile->obsType(sys, iType);
460 if (!_lli[prn].contains(iType)) {
461 _lli[prn][iType] = 0;
462 }
463 if (rnxSat.obs.contains(type) && rnxSat.obs[type].lli & 1) {
464 _lli[prn][iType] |= 1;
465 }
466 }
467 }
468}
469
470//
471////////////////////////////////////////////////////////////////////////////
472void t_reqcEdit::applyLLI(const t_rnxObsFile* obsFile,
473 t_rnxObsFile::t_rnxEpo* epo) {
474
475 if (_samplingRate == 0) {
476 return;
477 }
478
479 for (unsigned iSat = 0; iSat < epo->rnxSat.size(); iSat++) {
480 t_rnxObsFile::t_rnxSat& rnxSat = epo->rnxSat[iSat];
481 char sys = rnxSat.prn.system();
482 QString prn(rnxSat.prn.toString().c_str());
483
484 for (int iType = 0; iType < obsFile->nTypes(sys); iType++) {
485 QString type = obsFile->obsType(sys, iType);
486 if (_lli[prn].contains(iType) && _lli[prn][iType] & 1) {
487 if (rnxSat.obs.contains(type)) {
488 rnxSat.obs[type].lli |= 1;
489 }
490 }
491 }
492 }
493
494 _lli.clear();
495}
496
497/// Read All Ephemerides
498////////////////////////////////////////////////////////////////////////////
499void t_reqcEdit::readEphemerides(const QStringList& navFileNames,
500 QVector<t_eph*>& ephs) {
501
502 QStringListIterator it(navFileNames);
503 while (it.hasNext()) {
504 QString fileName = it.next();
505 if (fileName.indexOf('*') != -1 || fileName.indexOf('?') != -1) {
506 QFileInfo fileInfo(fileName);
507 QDir dir = fileInfo.dir();
508 QStringList filters; filters << fileInfo.fileName();
509 QListIterator<QFileInfo> it(dir.entryInfoList(filters));
510 while (it.hasNext()) {
511 QString filePath = it.next().filePath();
512 appendEphemerides(filePath, ephs);
513 }
514 }
515 else {
516 appendEphemerides(fileName, ephs);
517 }
518 }
519 qStableSort(ephs.begin(), ephs.end(), t_eph::earlierTime);
520}
521
522//
523////////////////////////////////////////////////////////////////////////////
524void t_reqcEdit::editEphemerides() {
525
526 // Easy Exit
527 // ---------
528 if (_navFileNames.isEmpty() || _outNavFileName.isEmpty()) {
529 return;
530 }
531
532 // Read Ephemerides
533 // ----------------
534 t_reqcEdit::readEphemerides(_navFileNames, _ephs);
535
536 // Check Satellite Systems
537 // -----------------------
538 bool haveGPS = false;
539 bool haveGlonass = false;
540 for (int ii = 0; ii < _ephs.size(); ii++) {
541 const t_eph* eph = _ephs[ii];
542 if (eph->type() == t_eph::GPS) {
543 haveGPS = true;
544 }
545 else if (eph->type() == t_eph::GLONASS) {
546 haveGlonass = true;
547 }
548 }
549
550 // Initialize output navigation file
551 // ---------------------------------
552 t_rnxNavFile outNavFile(_outNavFileName, t_rnxNavFile::output);
553
554 outNavFile.setGlonass(haveGlonass);
555
556 if ( (haveGPS && haveGlonass) || _rnxVersion >= 3.0) {
557 outNavFile.setVersion(t_rnxNavFile::defaultRnxNavVersion3);
558 }
559 else {
560 outNavFile.setVersion(t_rnxNavFile::defaultRnxNavVersion2);
561 }
562
563 bncSettings settings;
564 QMap<QString, QString> txtMap;
565 QString runBy = settings.value("reqcRunBy").toString();
566 if (!runBy.isEmpty()) {
567 txtMap["RUN BY"] = runBy;
568 }
569 QString comment = settings.value("reqcComment").toString();
570 if (!comment.isEmpty()) {
571 txtMap["COMMENT"] = comment;
572 }
573
574 outNavFile.writeHeader(&txtMap);
575
576 // Loop over all ephemerides
577 // -------------------------
578 for (int ii = 0; ii < _ephs.size(); ii++) {
579 const t_eph* eph = _ephs[ii];
580 bncTime begTime = _begTime;
581 bncTime endTime = _endTime;
582 if (eph->type() == t_eph::BDS) {
583 begTime += 14;
584 endTime += 14;
585 }
586 if (begTime.valid() && eph->TOC() < begTime) {
587 continue;
588 }
589 if (endTime.valid() && eph->TOC() > endTime) {
590 break;
591 }
592 outNavFile.writeEph(eph);
593 }
594}
595
596//
597////////////////////////////////////////////////////////////////////////////
598void t_reqcEdit::appendEphemerides(const QString& fileName,
599 QVector<t_eph*>& ephs) {
600
601 t_rnxNavFile rnxNavFile(fileName, t_rnxNavFile::input);
602 for (unsigned ii = 0; ii < rnxNavFile.ephs().size(); ii++) {
603 t_eph* eph = rnxNavFile.ephs()[ii];
604 bool isNew = true;
605 for (int iOld = 0; iOld < ephs.size(); iOld++) {
606 const t_eph* ephOld = ephs[iOld];
607 if (ephOld->prn() == eph->prn() &&
608 ephOld->TOC() == eph->TOC()) {
609 isNew = false;
610 break;
611 }
612 }
613 if (isNew) {
614 if (eph->type() == t_eph::GPS) {
615 ephs.append(new t_ephGPS(*dynamic_cast<t_ephGPS*>(eph)));
616 }
617 else if (eph->type() == t_eph::GLONASS) {
618 ephs.append(new t_ephGlo(*dynamic_cast<t_ephGlo*>(eph)));
619 }
620 else if (eph->type() == t_eph::Galileo) {
621 ephs.append(new t_ephGal(*dynamic_cast<t_ephGal*>(eph)));
622 }
623 else if (eph->type() == t_eph::QZSS) {
624 ephs.append(new t_ephGPS(*dynamic_cast<t_ephGPS*>(eph)));
625 }
626 else if (eph->type() == t_eph::SBAS) {
627 ephs.append(new t_ephSBAS(*dynamic_cast<t_ephSBAS*>(eph)));
628 }
629 else if (eph->type() == t_eph::BDS) {
630 ephs.append(new t_ephBDS(*dynamic_cast<t_ephBDS*>(eph)));
631 }
632 }
633 }
634}
Note: See TracBrowser for help on using the repository browser.