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
|
/***********************************************/
/**
* @file orbit2GroopsAscii.cpp
*
* @brief write orbit positions to TVGOGO ASCII format.
*
* @author Norbert Zehentner
* @date 2013-08-06
*
*/
/***********************************************/
// Latex documentation
#define DOCSTRING docstring
static const char *docstring = R"(
Convert groops orbits and corresponding covariance information to ASCII format.
The format is used to publish TUG orbits. It contains a two line header
with a short description of the orbit defined in \config{firstLine}.
The orbit is rotated to the Earth fixed frame (TRF) with \configClass{earthRotation}{earthRotationType} and given as one line per epoch.
The epoch lines contained time [MJD GPS time], position x, y and z [m], and the epoch covariance xx, yy, zz, xy, xz and yz [$m^2$].
See also \program{GroopsAscii2Orbit}.
)";
/***********************************************/
#include "programs/program.h"
#include "inputOutput/file.h"
#include "files/fileInstrument.h"
#include "classes/earthRotation/earthRotation.h"
/***** CLASS ***********************************/
/** @brief write orbit positions into TVGOGO ASCII file.
* @ingroup programsConversionGroup */
class Orbit2GroopsAscii
{
public:
void run(Config &config, Parallel::CommunicatorPtr comm);
};
GROOPS_REGISTER_PROGRAM(Orbit2GroopsAscii, SINGLEPROCESS, "write orbit positions to TUG ASCII format", Conversion, Orbit, Covariance, Instrument)
GROOPS_RENAMED_PROGRAM(Orbit2Ascii, Orbit2GroopsAscii, date2time(2020, 6, 14))
/***********************************************/
void Orbit2GroopsAscii::run(Config &config, Parallel::CommunicatorPtr /*comm*/)
{
try
{
FileName outName, inName, inNameCov;
std::string textLine1, textDatum;
EarthRotationPtr earthRotation;
readConfig(config, "outputfile", outName, Config::MUSTSET, "", "");
readConfig(config, "inputfileOrbit", inName, Config::MUSTSET, "", "");
readConfig(config, "inputfileCovariance", inNameCov, Config::MUSTSET, "", "");
readConfig(config, "earthRotation", earthRotation, Config::MUSTSET, "", "");
readConfig(config, "firstLine", textLine1, Config::OPTIONAL, "", "Text for first line");
if(isCreateSchema(config)) return;
logStatus<<"read orbit file <"<<inName<<">"<<Log::endl;
OrbitArc orbit = InstrumentFile::read(inName);
Covariance3dArc cov = InstrumentFile::read(inNameCov);
Arc::checkSynchronized({orbit, cov});
logStatus<<"write file <"<<outName<<">"<<Log::endl;
OutFile file(outName);
file<<"# "<<textLine1<<textDatum<<std::endl;
file<<"# MJD GPS time [days] x [m] y [m] z [m] xx [m^2] yy [m^2] zz [m^2] xy [m^2] xz [m^2] yz [m^2]"<<std::endl;
Single::forEach(orbit.size(), [&](UInt i)
{
const Rotary3d rot = earthRotation->rotaryMatrix(orbit.at(i).time);
const Vector3d position = rot.rotate(orbit.at(i).position);
Tensor3d cv;
if(cov.size())
cv = rot.rotate(cov.at(i).covariance);
file<<orbit.at(i).time%"%21.15f"s;
file<<position.x()%" %15.4f"s<<position.y()%" %15.4f"s<<position.z()%" %15.4f"s;
file<<cv.xx()%" %15.8e"s<<cv.yy()%" %15.8e"s<<cv.zz()%" %15.8e"s<<cv.xy()%" %15.8e"s<<cv.xz()%" %15.8e"s<<cv.yz()%" %15.8e"s;
file<<std::endl;
});
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
|
