1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
|
/***********************************************/
/**
* @file gnssGlonassFrequencyNumberUpdate.cpp
*
* @brief Update/set GLONASS frequency number in transmitter info files.
*
* @author Sebastian Strasser
* @date 2019-08-29
*/
/***********************************************/
// Latex documentation
#define DOCSTRING docstring
static const char *docstring = R"(
Update/set GLONASS frequency number in \configFile{inputfileTransmitterInfo}{platform} files.
PRN/SVN to frequency number source: \url{http://semisys.gfz-potsdam.de/semisys/api/?symname=2002&format=json&satellite=GLO}.
See also \program{GnssAntex2AntennaDefinition}.
)";
/***********************************************/
#include "programs/program.h"
#include "files/fileGnssReceiverDefinition.h"
#include "files/filePlatform.h"
#include "files/fileMatrix.h"
/***** CLASS ***********************************/
/** @brief Update/set GLONASS frequency number in transmitter info files.
* @ingroup programsGroup */
class GnssGlonassFrequencyNumberUpdate
{
public:
void run(Config &config, Parallel::CommunicatorPtr comm);
};
GROOPS_REGISTER_PROGRAM(GnssGlonassFrequencyNumberUpdate, SINGLEPROCESS, "Update/set GLONASS frequency number in transmitter info files.", Gnss)
/***********************************************/
void GnssGlonassFrequencyNumberUpdate::run(Config &config, Parallel::CommunicatorPtr /*comm*/)
{
try
{
FileName fileNameOutTransmitterInfo, fileNameInTransmitterInfo, fileNameInPrnSvn2FrequencyNumber;
std::vector<std::string> prnList;
std::string variableNamePrn;
readConfig(config, "outputfileTransmitterInfo", fileNameOutTransmitterInfo, Config::OPTIONAL, "", "templated for PRN list (variableNamePrn)");
readConfig(config, "inputfileTransmitterInfo", fileNameInTransmitterInfo, Config::MUSTSET, "", "templated for PRN list (variableNamePrn)");
readConfig(config, "inputfilePrn2FrequencyNumber", fileNameInPrnSvn2FrequencyNumber, Config::MUSTSET, "", "GROOPS matrix with columns: GLONASS PRN, SVN, mjdStart, mjdEnd, frequencyNumber");
readConfig(config, "prn", prnList, Config::OPTIONAL, "", "PRN (e.g. R01) for transmitter info files");
readConfig(config, "variableNamePrn", variableNamePrn, Config::OPTIONAL, "prn", "variable name for PRN in transmitter info files");
if(isCreateSchema(config)) return;
std::vector<Platform> transmitterInfos(prnList.size());
if(!fileNameInTransmitterInfo.empty())
{
VariableList varList;
varList.setVariable(variableNamePrn, "R**");
logStatus<<"read transmitter infos from <"<<fileNameInTransmitterInfo(varList)<<">"<< Log::endl;
for(UInt idPrn=0; idPrn<prnList.size(); idPrn++)
{
varList.setVariable(variableNamePrn, prnList.at(idPrn));
readFilePlatform(fileNameInTransmitterInfo(varList), transmitterInfos.at(idPrn));
}
}
logStatus<<"read GLONASS PRN/SVN to frequency number matrix from <"<<fileNameInPrnSvn2FrequencyNumber<<">"<< Log::endl;
Matrix prnSvn2FreqNo;
readFileMatrix(fileNameInPrnSvn2FrequencyNumber, prnSvn2FreqNo);
//---------------------------------------------------------------
std::vector<std::vector<std::shared_ptr<PlatformGnssReceiver>>> receiverInfos(prnList.size());
for(UInt i=0; i<prnSvn2FreqNo.rows(); i++)
{
const std::string prn = prnSvn2FreqNo(i,0)%"R%02i"s;
const std::string svn = prnSvn2FreqNo(i,1)%"R%03i"s;
const Time timeStart = mjd2time(prnSvn2FreqNo(i,2));
const Time timeEnd = mjd2time(prnSvn2FreqNo(i,3));
const std::string freqNo = prnSvn2FreqNo(i,4)%"%i"s;
Bool found = FALSE;
for(UInt idTrans=0; idTrans<transmitterInfos.size(); idTrans++)
if(transmitterInfos.at(idTrans).markerNumber == prn)
{
// new receiver list entry
auto recv = transmitterInfos.at(idTrans).findEquipment<PlatformGnssReceiver>(0.5*(timeStart+timeEnd));
if(!recv)
continue;
found = TRUE;
auto info = std::make_shared<PlatformGnssReceiver>(*recv);
info->timeStart = timeStart;
info->timeEnd = timeEnd;
info->version = freqNo;
receiverInfos.at(idTrans).push_back(info);
}
if(!found)
logWarning<<prn<<": no transmitter info entry found for time period "<<timeStart.dateTimeStr()<<" to "<<timeEnd.dateTimeStr()<<Log::endl;
}
//---------------------------------------------------------------
if(!fileNameOutTransmitterInfo.empty())
{
VariableList varList;
varList.setVariable(variableNamePrn, "R**");
logStatus<<"write transmitter infos to <"<<fileNameOutTransmitterInfo(varList)<<">"<< Log::endl;
for(UInt idPrn=0; idPrn<prnList.size(); idPrn++)
{
auto &eq = transmitterInfos.at(idPrn).equipments;
// remove old receivers
eq.erase(std::remove_if(eq.begin(), eq.end(), [](auto &x) {return std::dynamic_pointer_cast<PlatformGnssReceiver>(x);}), eq.end());
// insert new receivers
std::sort(receiverInfos.at(idPrn).begin(), receiverInfos.at(idPrn).end(), [](auto &info1, auto &info2){return info1->timeStart < info2->timeStart;});
eq.insert(eq.end(), receiverInfos.at(idPrn).begin(), receiverInfos.at(idPrn).end());
varList.setVariable(variableNamePrn, prnList.at(idPrn));
writeFilePlatform(fileNameOutTransmitterInfo(varList), transmitterInfos.at(idPrn));
}
}
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
|
