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
|
/***********************************************/
/**
* @file griddedData.h
*
* @brief Gridded values.
*
* @author Torsten Mayer-Guerr
* @date 2005-01-14
*
*/
/***********************************************/
#ifndef __GROOPS_GRIDDEDDATA__
#define __GROOPS_GRIDDEDDATA__
#include "base/importStd.h"
#include "base/vector3d.h"
#include "base/ellipsoid.h"
/** @addtogroup base */
/// @{
class GriddedDataRectangular;
/***** CLASS ***********************************/
/** @brief Point list with (multiple) data. */
class GriddedData
{
public:
Ellipsoid ellipsoid; //!< Ellipsoid for conversion to ellipsoidal coordinates
std::vector<Vector3d> points; //!< List of points
std::vector<Double> areas; //!< Area element (projected to unit sphere).
std::vector<std::vector<Double>> values; //!< data.at(dataIdx).at(pointIdx)
/// Default constructor.
GriddedData() = default;
/// Constructor with points, area elements, and multiple values for each point.
GriddedData(const Ellipsoid &ellip, const std::vector<Vector3d> &_points, const std::vector<Double> &_areas, const std::vector<std::vector<Double>> &_values) : ellipsoid(ellip), points(_points), areas(_areas), values(_values) {}
/// Constructor from GriddedDataRectangular.
GriddedData(const GriddedDataRectangular &grid) {init(grid);}
/// Create from GriddedDataRectangular.
void init(const GriddedDataRectangular &grid);
/** @brief Sort points geographically (North/West->South/East). */
void sort();
/** @brief Define points a rectangular grid?.
* if function returns FALSE, @a lambda, @a phi, @a radius contain garbage.
* if function returns TRUE, points are in same order as:
@code
UInt idx = 0;
for(UInt i=0; i<phi.size(); i++)
for(UInt k=0; k<lambda.size(); k++)
points.at(idx++) == polar(lambda.at(k), phi.at(i), r.at(i));
@endcode */
Bool isRectangle(std::vector<Angle> &lambda, std::vector<Angle> &phi, std::vector<Double> &radius) const;
/** @brief Automatically area computation of rectangular grids (overwrite areas). */
Bool computeArea();
/** @brief Is GriddedData valid?.
* Test dimensions of vectors. */
Bool isValid() const;
};
/***** CLASS ***********************************/
/** @brief Rectangular grid with (multiple) data. */
class GriddedDataRectangular
{
public:
Ellipsoid ellipsoid; //!< Ellipsoid for conversion to ellipsoidal coordinates
std::vector<Angle> longitudes; //!< Longitude (columns)
std::vector<Angle> latitudes; //!< Latitude (rows)
std::vector<Double> heights; //!< Elliposoidal height (rows)
std::vector<Matrix> values; //!< Multiple values at each point.
/// Create from GriddedData.
Bool init(const GriddedData &grid);
/** @brief Conversion from ellipsoidal to geocentric spherical polar coordinates. */
void geocentric(std::vector<Angle> &lambda, std::vector<Angle> &phi, std::vector<Double> &radius) const;
/** @brief borders of grid cell (i,k): (lat(i) - lat(i+1)) x (lon(k) - lon(k+1)). */
void cellBorders(std::vector<Double> &longitudes, std::vector<Double> &latitudes) const;
/** @brief Area elements projected on the unit sphere.
* area(i,k) = dPhi(i)*dLambda(k).
* @return total area (4pi for global grids). */
Double areaElements(std::vector<Double> &dLambda, std::vector<Double> &dPhi) const;
/** @brief Is GriddedDataRectangular valid?.
* Test dimensions of vectors. */
Bool isValid() const;
};
/// @}
/***********************************************/
#endif /* __GROOPS__ */
|
