Context Navigation

pppMain.cpp@ 10575

Last change on this file since 10575 was 9474, checked in by stuerze, 3 years ago
satellite antenna phase center correction is added
File size: 13.6 KB

Rev	Line
[5714]	1
	2	// Part of BNC, a utility for retrieving decoding and
	3	// converting GNSS data streams from NTRIP broadcasters.
	4	//
	5	// Copyright (C) 2007
	6	// German Federal Agency for Cartography and Geodesy (BKG)
	7	// http://www.bkg.bund.de
	8	// Czech Technical University Prague, Department of Geodesy
	9	// http://www.fsv.cvut.cz
	10	//
	11	// Email: euref-ip@bkg.bund.de
	12	//
	13	// This program is free software; you can redistribute it and/or
	14	// modify it under the terms of the GNU General Public License
	15	// as published by the Free Software Foundation, version 2.
	16	//
	17	// This program is distributed in the hope that it will be useful,
	18	// but WITHOUT ANY WARRANTY; without even the implied warranty of
	19	// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
	20	// GNU General Public License for more details.
	21	//
	22	// You should have received a copy of the GNU General Public License
	23	// along with this program; if not, write to the Free Software
	24	// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
	25
	26	/* -------------------------------------------------------------------------
	27	* BKG NTRIP Client
	28	* -------------------------------------------------------------------------
	29	*
	30	* Class: t_pppMain
	31	*
	32	* Purpose: Start of the PPP client(s)
	33	*
	34	* Author: L. Mervart
	35	*
	36	* Created: 29-Jul-2014
	37	*
[7278]	38	* Changes:
[5714]	39	*
	40	* -----------------------------------------------------------------------*/
	41
[5717]	42	#include <iostream>
	43
[5714]	44	#include "pppMain.h"
[5911]	45	#include "pppCrdFile.h"
[5714]	46	#include "bncsettings.h"
	47
[5814]	48	using namespace BNC_PPP;
[5714]	49	using namespace std;
	50
	51	// Constructor
	52	//////////////////////////////////////////////////////////////////////////////
	53	t_pppMain::t_pppMain() {
[5762]	54	_running = false;
[5714]	55	}
	56
	57	// Destructor
	58	//////////////////////////////////////////////////////////////////////////////
	59	t_pppMain::~t_pppMain() {
[5900]	60	stop();
[5975]	61	QListIterator<t_pppOptions*> iOpt(_options);
	62	while (iOpt.hasNext()) {
	63	delete iOpt.next();
	64	}
[5714]	65	}
	66
[7278]	67	//
[5714]	68	//////////////////////////////////////////////////////////////////////////////
[5946]	69	void t_pppMain::start() {
[5719]	70	if (_running) {
	71	return;
	72	}
	73
[5764]	74	try {
	75	readOptions();
[5719]	76
[5810]	77	QListIterator<t_pppOptions*> iOpt(_options);
[5764]	78	while (iOpt.hasNext()) {
[5810]	79	const t_pppOptions* opt = iOpt.next();
[5946]	80	t_pppThread* pppThread = new t_pppThread(opt);
	81	pppThread->start();
[5813]	82	_pppThreads << pppThread;
[5764]	83	_running = true;
	84	}
[5719]	85	}
[5825]	86	catch (t_except exc) {
[5765]	87	_running = true;
[5900]	88	stop();
[5764]	89	}
[5716]	90	}
	91
[7278]	92	//
[5716]	93	//////////////////////////////////////////////////////////////////////////////
[5900]	94	void t_pppMain::stop() {
[5871]	95
[5719]	96	if (!_running) {
	97	return;
	98	}
[5730]	99
[5972]	100	if (_realTime) {
	101	QListIterator<t_pppThread*> it(_pppThreads);
	102	while (it.hasNext()) {
	103	t_pppThread* pppThread = it.next();
	104	pppThread->exit();
[7856]	105	if (BNC_CORE->mode() != t_bncCore::interactive) {
	106	while(!pppThread->isFinished()) {
	107	pppThread->wait();
	108	}
	109	delete pppThread;
	110	}
[5972]	111	}
[5975]	112	_pppThreads.clear();
[5971]	113	}
[5730]	114
[5719]	115	_running = false;
[5716]	116	}
	117
[7278]	118	//
[5716]	119	//////////////////////////////////////////////////////////////////////////////
[5714]	120	void t_pppMain::readOptions() {
	121
[5810]	122	QListIterator<t_pppOptions*> iOpt(_options);
[5736]	123	while (iOpt.hasNext()) {
	124	delete iOpt.next();
	125	}
[5714]	126	_options.clear();
	127
	128	bncSettings settings;
	129
[5971]	130	_realTime = false;
[5714]	131	if (settings.value("PPP/dataSource").toString() == "Real-Time Streams") {
[5971]	132	_realTime = true;
[5714]	133	}
	134	else if (settings.value("PPP/dataSource").toString() == "RINEX Files") {
[5971]	135	_realTime = false;
[5714]	136	}
	137	else {
[7278]	138	return;
[5714]	139	}
	140
[5736]	141	QListIterator<QString> iSta(settings.value("PPP/staTable").toStringList());
	142	while (iSta.hasNext()) {
	143	QStringList hlp = iSta.next().split(",");
[5714]	144
[6964]	145	if (hlp.size() < 10) {
[5825]	146	throw t_except("pppMain: wrong option staTable");
[5714]	147	}
	148
[5810]	149	t_pppOptions* opt = new t_pppOptions();
[5714]	150
[5971]	151	opt->_realTime = _realTime;
[7278]	152	opt->_roverName = hlp[0].toStdString();
[8101]	153	opt->_aprSigCrd[0] = hlp[1].toDouble()+1e-10;
	154	opt->_aprSigCrd[1] = hlp[2].toDouble()+1e-10;
	155	opt->_aprSigCrd[2] = hlp[3].toDouble()+1e-10;
	156	opt->_noiseCrd[0] = hlp[4].toDouble()+1e-10;
	157	opt->_noiseCrd[1] = hlp[5].toDouble()+1e-10;
	158	opt->_noiseCrd[2] = hlp[6].toDouble()+1e-10;
[5913]	159	opt->_aprSigTrp = hlp[7].toDouble();
	160	opt->_noiseTrp = hlp[8].toDouble();
[6964]	161	opt->_nmeaPort = hlp[9].toInt();
[5714]	162
[5971]	163	if (_realTime) {
[7850]	164	opt->_corrMount.assign(settings.value("PPP/corrMount").toString().toStdString());
[9474]	165	opt->_isAPC = (opt->_corrMount.substr(0,4)=="SSRA");
[5714]	166	}
	167	else {
[7850]	168	opt->_rinexObs.assign(settings.value("PPP/rinexObs").toString().toStdString());
	169	opt->_rinexNav.assign(settings.value("PPP/rinexNav").toString().toStdString());
	170	opt->_corrFile.assign(settings.value("PPP/corrFile").toString().toStdString());
[9474]	171	QFileInfo tmp = QFileInfo(QString::fromStdString(opt->_corrFile));
	172	opt->_isAPC = (tmp.baseName().mid(0,4)=="SSRA");
[5714]	173	}
	174
[7850]	175	opt->_crdFile.assign(settings.value("PPP/crdFile").toString().toStdString());
	176	opt->_antexFileName.assign(settings.value("PPP/antexFile").toString().toStdString());
[5714]	177
[5997]	178	opt->_sigmaC1 = settings.value("PPP/sigmaC1").toDouble(); if (opt->_sigmaC1 <= 0.0) opt->_sigmaC1 = 2.0;
	179	opt->_sigmaL1 = settings.value("PPP/sigmaL1").toDouble(); if (opt->_sigmaL1 <= 0.0) opt->_sigmaL1 = 0.01;
[5714]	180	opt->_corrWaitTime = settings.value("PPP/corrWaitTime").toDouble();
[7945]	181	if (!_realTime \|\| opt->_corrMount.empty()) {
	182	opt->_corrWaitTime = 0;
	183	}
[5714]	184
[7261]	185	if (settings.value("PPP/lcGPS").toString() == "P1") {
	186	opt->_LCsGPS.push_back(t_lc::c1);
	187	}
	188	else if (settings.value("PPP/lcGPS").toString() == "L1") {
	189	opt->_LCsGPS.push_back(t_lc::l1);
	190	}
	191	else if (settings.value("PPP/lcGPS").toString() == "P1&L1") {
	192	opt->_LCsGPS.push_back(t_lc::c1);
	193	opt->_LCsGPS.push_back(t_lc::l1);
	194	}
	195	if (settings.value("PPP/lcGPS").toString() == "P2") {
	196	opt->_LCsGPS.push_back(t_lc::c2);
	197	}
	198	else if (settings.value("PPP/lcGPS").toString() == "L2") {
	199	opt->_LCsGPS.push_back(t_lc::l2);
	200	}
	201	else if (settings.value("PPP/lcGPS").toString() == "P2&L2") {
	202	opt->_LCsGPS.push_back(t_lc::c2);
	203	opt->_LCsGPS.push_back(t_lc::l2);
	204	}
	205	if (settings.value("PPP/lcGPS").toString() == "Pi") {
	206	opt->_LCsGPS.push_back(t_lc::c1);
	207	opt->_LCsGPS.push_back(t_lc::c2);
	208	}
	209	else if (settings.value("PPP/lcGPS").toString() == "Li") {
	210	opt->_LCsGPS.push_back(t_lc::l1);
	211	opt->_LCsGPS.push_back(t_lc::l2);
	212	}
	213	else if (settings.value("PPP/lcGPS").toString() == "Pi&Li") {
	214	opt->_LCsGPS.push_back(t_lc::c1);
	215	opt->_LCsGPS.push_back(t_lc::c2);
	216	opt->_LCsGPS.push_back(t_lc::l1);
	217	opt->_LCsGPS.push_back(t_lc::l2);
	218	}
[5715]	219	if (settings.value("PPP/lcGPS").toString() == "P3") {
[5913]	220	opt->_LCsGPS.push_back(t_lc::cIF);
[5715]	221	}
	222	else if (settings.value("PPP/lcGPS").toString() == "L3") {
[5913]	223	opt->_LCsGPS.push_back(t_lc::lIF);
[5715]	224	}
	225	else if (settings.value("PPP/lcGPS").toString() == "P3&L3") {
[5913]	226	opt->_LCsGPS.push_back(t_lc::cIF);
	227	opt->_LCsGPS.push_back(t_lc::lIF);
[5715]	228	}
	229
[7261]	230	if (settings.value("PPP/lcGLONASS").toString() == "P1") {
	231	opt->_LCsGLONASS.push_back(t_lc::c1);
	232	}
	233	else if (settings.value("PPP/lcGLONASS").toString() == "L1") {
	234	opt->_LCsGLONASS.push_back(t_lc::l1);
	235	}
	236	else if (settings.value("PPP/lcGLONASS").toString() == "P1&L1") {
	237	opt->_LCsGLONASS.push_back(t_lc::c1);
	238	opt->_LCsGLONASS.push_back(t_lc::l1);
	239	}
	240	if (settings.value("PPP/lcGLONASS").toString() == "P2") {
	241	opt->_LCsGLONASS.push_back(t_lc::c2);
	242	}
	243	else if (settings.value("PPP/lcGLONASS").toString() == "L2") {
	244	opt->_LCsGLONASS.push_back(t_lc::l2);
	245	}
	246	else if (settings.value("PPP/lcGLONASS").toString() == "P2&L2") {
	247	opt->_LCsGLONASS.push_back(t_lc::c2);
	248	opt->_LCsGLONASS.push_back(t_lc::l2);
	249	}
	250	if (settings.value("PPP/lcGLONASS").toString() == "Pi") {
	251	opt->_LCsGLONASS.push_back(t_lc::c1);
	252	opt->_LCsGLONASS.push_back(t_lc::c2);
	253	}
	254	else if (settings.value("PPP/lcGLONASS").toString() == "Li") {
	255	opt->_LCsGLONASS.push_back(t_lc::l1);
	256	opt->_LCsGLONASS.push_back(t_lc::l2);
	257	}
	258	else if (settings.value("PPP/lcGLONASS").toString() == "Pi&Li") {
	259	opt->_LCsGLONASS.push_back(t_lc::c1);
	260	opt->_LCsGLONASS.push_back(t_lc::c2);
	261	opt->_LCsGLONASS.push_back(t_lc::l1);
	262	opt->_LCsGLONASS.push_back(t_lc::l2);
	263	}
[5715]	264	if (settings.value("PPP/lcGLONASS").toString() == "P3") {
[5913]	265	opt->_LCsGLONASS.push_back(t_lc::cIF);
[5715]	266	}
	267	else if (settings.value("PPP/lcGLONASS").toString() == "L3") {
[5913]	268	opt->_LCsGLONASS.push_back(t_lc::lIF);
[5715]	269	}
	270	else if (settings.value("PPP/lcGLONASS").toString() == "P3&L3") {
[5913]	271	opt->_LCsGLONASS.push_back(t_lc::cIF);
	272	opt->_LCsGLONASS.push_back(t_lc::lIF);
[5715]	273	}
	274
[7261]	275	if (settings.value("PPP/lcGalileo").toString() == "P1") {
	276	opt->_LCsGalileo.push_back(t_lc::c1);
	277	}
	278	else if (settings.value("PPP/lcGalileo").toString() == "L1") {
	279	opt->_LCsGalileo.push_back(t_lc::l1);
	280	}
	281	else if (settings.value("PPP/lcGalileo").toString() == "P1&L1") {
	282	opt->_LCsGalileo.push_back(t_lc::c1);
	283	opt->_LCsGalileo.push_back(t_lc::l1);
	284	}
	285	if (settings.value("PPP/lcGalileo").toString() == "P5") {
	286	opt->_LCsGalileo.push_back(t_lc::c2);
	287	}
	288	else if (settings.value("PPP/lcGalileo").toString() == "L5") {
	289	opt->_LCsGalileo.push_back(t_lc::l2);
	290	}
	291	else if (settings.value("PPP/lcGalileo").toString() == "P5&L5") {
	292	opt->_LCsGalileo.push_back(t_lc::c2);
	293	opt->_LCsGalileo.push_back(t_lc::l2);
	294	}
	295	if (settings.value("PPP/lcGalileo").toString() == "Pi") {
	296	opt->_LCsGalileo.push_back(t_lc::c1);
	297	opt->_LCsGalileo.push_back(t_lc::c2);
	298	}
	299	else if (settings.value("PPP/lcGalileo").toString() == "Li") {
	300	opt->_LCsGalileo.push_back(t_lc::l1);
	301	opt->_LCsGalileo.push_back(t_lc::l2);
	302	}
	303	else if (settings.value("PPP/lcGalileo").toString() == "Pi&Li") {
	304	opt->_LCsGalileo.push_back(t_lc::c1);
	305	opt->_LCsGalileo.push_back(t_lc::c2);
	306	opt->_LCsGalileo.push_back(t_lc::l1);
	307	opt->_LCsGalileo.push_back(t_lc::l2);
	308	}
[5715]	309	if (settings.value("PPP/lcGalileo").toString() == "P3") {
[5913]	310	opt->_LCsGalileo.push_back(t_lc::cIF);
[5715]	311	}
	312	else if (settings.value("PPP/lcGalileo").toString() == "L3") {
[5913]	313	opt->_LCsGalileo.push_back(t_lc::lIF);
[5715]	314	}
	315	else if (settings.value("PPP/lcGalileo").toString() == "P3&L3") {
[5913]	316	opt->_LCsGalileo.push_back(t_lc::cIF);
	317	opt->_LCsGalileo.push_back(t_lc::lIF);
[5715]	318	}
	319
[7261]	320	if (settings.value("PPP/lcBDS").toString() == "P2") {
	321	opt->_LCsBDS.push_back(t_lc::c1);
	322	}
	323	else if (settings.value("PPP/lcBDS").toString() == "L2") {
	324	opt->_LCsBDS.push_back(t_lc::l1);
	325	}
	326	else if (settings.value("PPP/lcBDS").toString() == "P2&L2") {
	327	opt->_LCsBDS.push_back(t_lc::c1);
	328	opt->_LCsBDS.push_back(t_lc::l1);
	329	}
	330	if (settings.value("PPP/lcBDS").toString() == "P7") {
	331	opt->_LCsBDS.push_back(t_lc::c2);
	332	}
	333	else if (settings.value("PPP/lcBDS").toString() == "L7") {
	334	opt->_LCsBDS.push_back(t_lc::l2);
	335	}
	336	else if (settings.value("PPP/lcBDS").toString() == "P7&L7") {
	337	opt->_LCsBDS.push_back(t_lc::c2);
	338	opt->_LCsBDS.push_back(t_lc::l2);
	339	}
	340	if (settings.value("PPP/lcBDS").toString() == "Pi") {
	341	opt->_LCsBDS.push_back(t_lc::c1);
	342	opt->_LCsBDS.push_back(t_lc::c2);
	343	}
	344	else if (settings.value("PPP/lcBDS").toString() == "Li") {
	345	opt->_LCsBDS.push_back(t_lc::l1);
	346	opt->_LCsBDS.push_back(t_lc::l2);
	347	}
	348	else if (settings.value("PPP/lcBDS").toString() == "Pi&Li") {
	349	opt->_LCsBDS.push_back(t_lc::c1);
	350	opt->_LCsBDS.push_back(t_lc::c2);
	351	opt->_LCsBDS.push_back(t_lc::l1);
	352	opt->_LCsBDS.push_back(t_lc::l2);
	353	}
[6964]	354	if (settings.value("PPP/lcBDS").toString() == "P3") {
	355	opt->_LCsBDS.push_back(t_lc::cIF);
	356	}
	357	else if (settings.value("PPP/lcBDS").toString() == "L3") {
	358	opt->_LCsBDS.push_back(t_lc::lIF);
	359	}
	360	else if (settings.value("PPP/lcBDS").toString() == "P3&L3") {
	361	opt->_LCsBDS.push_back(t_lc::cIF);
	362	opt->_LCsBDS.push_back(t_lc::lIF);
	363	}
	364
[5911]	365	// Information from the coordinate file
	366	// ------------------------------------
[7850]	367	string crdFileName(settings.value("PPP/crdFile").toString().toStdString());
[5911]	368	if (!crdFileName.empty()) {
	369	vector<t_pppCrdFile::t_staInfo> staInfoVec;
	370	t_pppCrdFile::readCrdFile(crdFileName, staInfoVec);
	371	for (unsigned ii = 0; ii < staInfoVec.size(); ii++) {
	372	const t_pppCrdFile::t_staInfo& staInfo = staInfoVec[ii];
	373	if (staInfo._name == opt->_roverName) {
	374	opt->_xyzAprRover[0] = staInfo._xyz[0];
	375	opt->_xyzAprRover[1] = staInfo._xyz[1];
	376	opt->_xyzAprRover[2] = staInfo._xyz[2];
	377	opt->_neuEccRover[0] = staInfo._neuAnt[0];
	378	opt->_neuEccRover[1] = staInfo._neuAnt[1];
	379	opt->_neuEccRover[2] = staInfo._neuAnt[2];
	380	opt->_antNameRover = staInfo._antenna;
[7523]	381	opt->_recNameRover = staInfo._receiver;
[5911]	382	break;
	383	}
	384	}
	385	}
	386
[5997]	387	opt->_minObs = settings.value("PPP/minObs").toInt(); if (opt->_minObs < 4) opt->_minObs = 4;
[5959]	388	opt->_minEle = settings.value("PPP/minEle").toDouble() * M_PI / 180.0;
[5997]	389	opt->_maxResC1 = settings.value("PPP/maxResC1").toDouble(); if (opt->_maxResC1 <= 0.0) opt->_maxResC1 = 3.0;
	390	opt->_maxResL1 = settings.value("PPP/maxResL1").toDouble(); if (opt->_maxResL1 <= 0.0) opt->_maxResL1 = 0.03;
[5959]	391	opt->_eleWgtCode = (settings.value("PPP/eleWgtCode").toInt() != 0);
	392	opt->_eleWgtPhase = (settings.value("PPP/eleWgtPhase").toInt() != 0);
	393	opt->_seedingTime = settings.value("PPP/seedingTime").toDouble();
	394
[5793]	395	// Some default values
	396	// -------------------
[7965]	397	opt->_aprSigAmb = 1000.0;
[5920]	398	opt->_noiseClk = 1000.0;
[5793]	399
[5714]	400	_options << opt;
	401	}
	402	}
	403

Note: See TracBrowser for help on using the repository browser.

Context Navigation

source: ntrip/branches/BNC_2.12/src/pppMain.cpp@ 10575

Download in other formats: