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
|
/***********************************************/
/**
* @file goceXml2Gradiometer.cpp
*
* @brief Read ESA XML GOCE Data.
*
* @author Torsten Mayer-Guerr
* @date 2010-10-19
*
*/
/***********************************************/
// Latex documentation
#define DOCSTRING docstring
static const char *docstring = R"(
Read ESA XML GOCE Data.
The \config{outputfileGradiometer} is written as \file{instrument file (GRADIOMETER)}{instrument}.
)";
/***********************************************/
#include "programs/program.h"
#include "parser/xml.h"
#include "files/fileInstrument.h"
/***** CLASS ***********************************/
/** @brief Read ESA XML GOCE Data.
* @ingroup programsConversionGroup */
class GoceXml2Gradiometer
{
void readFileGoceGradiometer(const FileName &fileName, GradiometerArc &arc);
public:
void run(Config &config, Parallel::CommunicatorPtr comm);
};
GROOPS_REGISTER_PROGRAM(GoceXml2Gradiometer, SINGLEPROCESS, "read ESA XML GOCE Data", Conversion, Instrument)
/***********************************************/
void GoceXml2Gradiometer::run(Config &config, Parallel::CommunicatorPtr /*comm*/)
{
try
{
FileName outName;
std::vector<FileName> fileName;
readConfig(config, "outputfileGradiometer", outName, Config::MUSTSET, "", "");
readConfig(config, "inputfile", fileName, Config::MUSTSET, "", "");
if(isCreateSchema(config)) return;
GradiometerArc arc;
for(UInt i=0; i<fileName.size(); i++)
{
logStatus<<"read file <"<<fileName.at(i)<<">"<<Log::endl;
readFileGoceGradiometer(fileName.at(i), arc);
}
logStatus<<"write gradiometer file <"<<outName<<">"<<Log::endl;
InstrumentFile::write(outName, arc);
Arc::printStatistics(arc);
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
void GoceXml2Gradiometer::readFileGoceGradiometer(const FileName &fileName, GradiometerArc &arc)
{
try
{
GradiometerEpoch epoch;
InFile file(fileName);
XmlNodePtr rootNode = XmlNode::read(file);
XmlNodePtr dataNode = getChild(rootNode, "Data_Block", TRUE);
XmlNodePtr ggtNode = getChild(dataNode, "EGG_GGT_DS", TRUE);
UInt epochCount = childCount(ggtNode, "EGG_GGT_1i", TRUE);
for(UInt i=0; i<epochCount; i++)
{
XmlNodePtr epochNode = getChild(ggtNode, "EGG_GGT_1i", TRUE);
Double t = 0;
readXml(epochNode, "Tt_GPS", t, TRUE);
epoch.time = seconds2time(t) + date2time(1980,1,6);
XmlNodePtr tnsNode = getChild(epochNode, "Gravity_Grad_Tensor", TRUE);
std::string tnsStr;
readXml(tnsNode, "U_G", tnsStr, TRUE);
std::stringstream ss(tnsStr);
ss>>epoch.gravityGradient.xx();
ss>>epoch.gravityGradient.yy();
ss>>epoch.gravityGradient.zz();
ss>>epoch.gravityGradient.xy();
ss>>epoch.gravityGradient.xz();
ss>>epoch.gravityGradient.yz();
arc.push_back(epoch);
}
}
catch(std::exception &e)
{
GROOPS_RETHROW(e)
}
}
/***********************************************/
/***********************************************/
|
