File: gnssSignalBias2Matrix.cpp

package info (click to toggle)
groops 0%2Bgit20250907%2Bds-1
  • links: PTS, VCS
  • area: non-free
  • in suites: forky, sid
  • size: 11,140 kB
  • sloc: cpp: 135,607; fortran: 1,603; makefile: 20
file content (81 lines) | stat: -rw-r--r-- 2,675 bytes parent folder | download | duplicates (2)
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
/***********************************************/
/**
* @file gnssSignalBias2Matrix.cpp
*
* @brief Computes signal biases for a given type list.
*
* @author Torsten Mayer-Guerr
* @date 2019-03-26
*/
/***********************************************/

// Latex documentation
#define DOCSTRING docstring
static const char *docstring = R"(
Computes signal biases for a given list of \configClass{types}{gnssType}.
If the type list is empty, all types contained in \configFile{inputfileSignalBias}{gnssSignalBias} are used.
The resulting \configFile{outputfileMatrix}{matrix} contains a vector with an entry for each type.
)";

/***********************************************/

#include "programs/program.h"
#include "files/fileGnssSignalBias.h"
#include "files/fileMatrix.h"
#include "files/fileStringTable.h"

/***** CLASS ***********************************/

/** @brief Computes signal biases for a given type list.
* @ingroup programsGroup */
class GnssSignalBias2Matrix
{
public:
  void run(Config &config, Parallel::CommunicatorPtr comm);
};

GROOPS_REGISTER_PROGRAM(GnssSignalBias2Matrix, SINGLEPROCESS, "Computes signal biases for a given type list", Gnss)


/***********************************************/

void GnssSignalBias2Matrix::run(Config &config, Parallel::CommunicatorPtr /*comm*/)
{
  try
  {
    FileName fileNameSignalBias, fileNameMatrix, fileNameTypes;
    std::vector<GnssType> types;

    readConfig(config, "outputfileMatrix",    fileNameMatrix,     Config::MUSTSET,  "", "");
    readConfig(config, "outputfileTypes",     fileNameTypes,      Config::OPTIONAL, "", "ASCII list of types");
    readConfig(config, "inputfileSignalBias", fileNameSignalBias, Config::MUSTSET,  "", "");
    readConfig(config, "types",               types,              Config::OPTIONAL, "", "if not set, all types in the file are used");
    if(isCreateSchema(config)) return;

    logStatus<<"read signal bias file <"<<fileNameSignalBias<<">"<<Log::endl;
    GnssSignalBias signalBias;
    readFileGnssSignalBias(fileNameSignalBias, signalBias);

    if(types.size() == 0)
      types = signalBias.types;
    Vector bias = signalBias.compute(types);

    if(!fileNameTypes.empty())
    {
      logStatus<<"write types to file <"<<fileNameTypes<<">"<<Log::endl;
      std::vector<std::string> typeNames;
      for(const auto &type : types)
        typeNames.push_back(type.str());
      writeFileStringList(fileNameTypes, typeNames);
    }

    logStatus<<"write matrix to file <"<<fileNameMatrix<<">"<<Log::endl;
    writeFileMatrix(fileNameMatrix, bias);
  }
  catch(std::exception &e)
  {
    GROOPS_RETHROW(e)
  }
}

/***********************************************/