36 #ifndef GEOTESSPOSITION_OBJECT_H
37 #define GEOTESSPOSITION_OBJECT_H
67 class GeoTessInterpolatorType;
151 void updatePosition2D(
int layerid,
const double*
const uVector);
170 void checkTessellation(
const int& tid)
172 if (triangle[tid] < 0)
175 triangle[tid] = grid.getTriangle(tid, 0, 0);
176 getContainingTriangle(tid);
188 void updateRadius(
int lid,
double rad)
190 if ((radius < 0.) || (rad != radius) || (lid != layerId))
194 tessid = layerTessIds[lid];
195 checkTessellation(tessid);
196 clearRadialCoefficients();
200 void clearRadialCoefficients()
202 for (
int i=0; i<(int)radialIndexes.size(); ++i)
204 radialIndexes[i].clear();
205 radialCoefficients[i].clear();
215 void updateRadialCoefficients()
219 if (radialIndexes.size() < vertices[tessid].size())
221 radialIndexes.resize(vertices[tessid].size());
222 radialCoefficients.resize(vertices[tessid].size());
227 if (radialIndexes[0].size() == 0)
229 vector<int>& v = vertices[tessid];
230 for (
int i = 0; i < (int)v.size(); ++i)
231 modlProfiles[v[i]][layerId]->setInterpolationCoefficients(radialInterpolatorType,
232 radialIndexes[i], radialCoefficients[i], radius, radiusOutOfRangeAllowed);
240 static const double TWALK_TOLERANCE;
252 void getContainingTriangle(
int tid);
262 int biggerTriangle(
int t1,
int t2)
264 const int* tv = grid.getTriangleVertexIndexes(t1);
265 double dot1 = GeoTessUtils::dot(grid.getVertex(tv[0]), grid.getVertex(tv[1]))
266 + GeoTessUtils::dot(grid.getVertex(tv[1]), grid.getVertex(tv[2]))
267 + GeoTessUtils::dot(grid.getVertex(tv[2]), grid.getVertex(tv[0]));
269 tv = grid.getTriangleVertexIndexes(t2);
270 double dot2 = GeoTessUtils::dot(grid.getVertex(tv[0]), grid.getVertex(tv[1]))
271 + GeoTessUtils::dot(grid.getVertex(tv[1]), grid.getVertex(tv[2]))
272 + GeoTessUtils::dot(grid.getVertex(tv[2]), grid.getVertex(tv[0]));
274 return dot2 > dot1 ? t1 : t2;
299 vector<double> layerRadii;
318 vector< vector< int> > vertices;
325 vector< vector< double> > linearCoefficients;
338 vector< vector< double> > hCoefficients;
354 vector<vector<int> > radialIndexes;
373 vector<vector<double> > radialCoefficients;
397 double const*
const* gridVertices;
398 int const*
const* gridTriangles;
399 int const* gridDescendants;
400 const vector<vector<Edge*> >& gridEdges;
401 const int* layerTessIds;
418 bool radiusOutOfRangeAllowed;
423 double unitVector[3];
428 int getTriangle(
int tid)
429 { checkTessellation(tid);
return triangle[tid]; };
435 virtual void update2D(
int tid) =
ABSTRACT;
509 virtual double getValue(
int attribute);
537 bool noEmptyProfiles();
559 void set(
double lat,
double lon,
double depth)
561 double uVector[3] = {0.0, 0.0, 0.0};
562 GeoTessUtils::getVectorDegrees(lat, lon, uVector);
563 double newRadius = GeoTessUtils::getEarthRadius(uVector) - depth;
564 set(uVector, newRadius);
579 void set(
const double*
const uVector,
const double& newRadius)
581 updatePosition2D(nLayers - 1, uVector);
583 int lid = getLayerId(newRadius);
585 updatePosition2D(lid, uVector);
586 updateRadius(lid, newRadius);
606 void set(
int layid,
double lat,
double lon,
double depth)
609 set(lat, lon, depth);
613 GeoTessUtils::getVectorDegrees(lat, lon, uVector);
614 updatePosition2D(layid, uVector);
615 updateRadius(layid, GeoTessUtils::getEarthRadius(uVector) - depth);
629 void set(
int layid,
const double*
const uVector,
double rad)
635 updatePosition2D(layid, uVector);
636 if (rad >= 0.0) updateRadius(layid, rad);
647 void setTop(
int layid,
const double*
const uVector)
649 updatePosition2D(layid, uVector);
650 updateRadius(layid, getRadiusTop(layid));
663 updatePosition2D(layid, uVector);
664 updateRadius(layid, getRadiusBottom(layid));
679 os << endl <<
"ERROR in GeoTessPosition::setRadius" << endl
680 <<
"Geographic position has not been specified." << endl;
683 updateRadius(layid, rad);
698 os << endl <<
"ERROR in GeoTessPosition::setRadius" << endl
699 <<
"Geographic position has not been specified." << endl;
702 updateRadius(getLayerId(rad), rad);
718 setRadius(layer, getEarthRadius()-depth);
729 setRadius(getEarthRadius()-depth);
742 os << endl <<
"ERROR in GeoTessPosition::setRadius" << endl
743 <<
"Geographic position has not been specified." << endl;
746 tessid = layerTessIds[layid];
747 checkTessellation(tessid);
748 updateRadius(layid, getRadiusTop(layid));
762 os << endl <<
"ERROR in GeoTessPosition::setRadius" << endl
763 <<
"Geographic position has not been specified." << endl;
766 tessid = layerTessIds[layid];
767 checkTessellation(tessid);
768 updateRadius(layid, getRadiusBottom(layid));
778 double getRadiusTop(
int layid);
786 double getRadiusBottom(
int layid);
795 if (earthRadius < 0.)
796 earthRadius = GeoTessUtils::getEarthRadius(unitVector);
815 { u[0] = unitVector[0]; u[1] = unitVector[1]; u[2] = unitVector[2]; };
837 const vector<int>&
getVertices()
const {
return vertices[tessid]; };
844 int getIndexOfClosestVertex()
const;
854 return grid.getVertex(getIndexOfClosestVertex());
864 {
return vertices[tessid][index]; }
873 coefficients.clear();
874 vector<int>& vtx = vertices[tessid];
875 vector<double>& hc = hCoefficients[tessid];
876 for (
int i = 0; i < (int) vtx.size(); ++i)
877 modlProfiles[vtx[i]][layerId]->getCoefficients(coefficients, radius, hc[i]);
891 maxTessLevel[layerTessIds[layid]] = maxTess;
892 triangle[layerTessIds[layid]] = -1;
893 checkTessellation(tessid);
906 {
return maxTessLevel[layerTessIds[layid]]; }
924 int getTessLevel(
const int& tId) { checkTessellation(tId);
return tessLevels[tId]; }
932 {
return getRadiusTop(layerId); };
940 {
return getRadiusBottom(layerId); };
949 return getEarthRadius() - getRadiusTop(layerId);
959 return getEarthRadius() - getRadiusBottom(layerId);
970 return getEarthRadius() - getRadiusTop(layid);
981 return getEarthRadius() - getRadiusBottom(layid);
992 return getRadiusTop(layid) - getRadiusBottom(layid);
1002 return getRadiusTop() - getRadiusBottom();
1018 {
return getEarthRadius() - radius; };
1042 for (
int i=0; i<nLayers; ++i)
1043 if (rad <= getRadiusTop(i))
1046 for (
int i=nLayers-1; i>=0; --i)
1047 if (getLayerThickness(i) > 0)
1083 errorValue = errVal;
1095 vector<int>& vtid = vertices[tessid];
1096 vector<double>& htid = hCoefficients[tessid];
1097 for (
int v = 0; v < (int) vtid.size(); ++v)
1098 if (htid[v] > 0.999999999)
1129 vector<int>& vtx = vertices[tessid];
1130 vector<double>& htid = hCoefficients[tessid];
1132 for (
int i = 0; i < (int) vtx.size(); ++i)
1133 modlProfiles[vtx[i]][layerId]->getWeights(weights, dkm, radius, htid[i]);
1145 return hCoefficients[tessid];
1158 return hCoefficients[tessid][index];
1178 return radiusOutOfRangeAllowed;
1198 if (allowed != radiusOutOfRangeAllowed)
1199 clearRadialCoefficients();
1201 radiusOutOfRangeAllowed = allowed;
1214 if (isNotReferenced())
1217 os << endl <<
"ERROR in Polygon::removeReference" << endl
1218 <<
"Reference count (" << refCount <<
") is already zero." << endl;
1245 #endif // GEOTESSPOSITION_OBJECT_H
void set(const double *const uVector, const double &newRadius)
Definition: GeoTessPosition.h:579
void copyVector(double *u)
Definition: GeoTessPosition.h:814
double getDepthBottom(int layid)
Definition: GeoTessPosition.h:979
double * getVector()
Definition: GeoTessPosition.h:806
double getEarthRadius()
Definition: GeoTessPosition.h:793
void setBottom(int layid, const double *const uVector)
Definition: GeoTessPosition.h:661
bool isRadiusOutOfRangeAllowed()
Definition: GeoTessPosition.h:1176
double getDepthTop(int layid)
Definition: GeoTessPosition.h:968
double getDepthTop()
Definition: GeoTessPosition.h:947
int getTessLevel() const
Definition: GeoTessPosition.h:914
const double * getClosestVertex() const
Definition: GeoTessPosition.h:852
Information about an interpolated point at an arbitrary position in a model.
Definition: GeoTessPosition.h:101
double getLayerThickness()
Definition: GeoTessPosition.h:1000
void setDepth(double depth)
Definition: GeoTessPosition.h:727
const vector< double > & getHorizontalCoefficients() const
Definition: GeoTessPosition.h:1143
void setTop(int layid, const double *const uVector)
Definition: GeoTessPosition.h:647
Enumeration of the interpolation algorithms supported by GeoTess including LINEAR, NATURAL_NEIGHBOR and CUBIC_SPLINE.
Definition: GeoTessInterpolatorType.h:71
Manages the geometry and topology of one or more multi-level triangular tessellations of a unit spher...
Definition: GeoTessGrid.h:163
void setTop(int layid)
Definition: GeoTessPosition.h:737
void set(int layid, double lat, double lon, double depth)
Definition: GeoTessPosition.h:606
int getLayerId(double rad)
Definition: GeoTessPosition.h:1040
int getVertexIndex()
Definition: GeoTessPosition.h:1093
double getHorizontalCoefficient(int index) const
Definition: GeoTessPosition.h:1156
void setBottom(int layid)
Definition: GeoTessPosition.h:757
int getReferenceCount()
Definition: GeoTessPosition.h:1230
void setErrorValue(double errVal)
Definition: GeoTessPosition.h:1081
void setDepth(int layer, double depth)
Definition: GeoTessPosition.h:716
int getVertex(int index)
Definition: GeoTessPosition.h:863
int getTessLevel(const int &tId)
Definition: GeoTessPosition.h:924
An exception class for all GeoTess objects.
Definition: GeoTessException.h:65
void set(int layid, const double *const uVector, double rad)
Definition: GeoTessPosition.h:629
double getErrorValue()
Definition: GeoTessPosition.h:1070
void addReference()
Definition: GeoTessPosition.h:1207
void removeReference()
Definition: GeoTessPosition.h:1212
void setRadius(int layid, double rad)
Definition: GeoTessPosition.h:674
int getNVertices()
Definition: GeoTessPosition.h:828
double getRadius()
Definition: GeoTessPosition.h:1010
#define ABSTRACT
Global constant used to make pure virtual functions readable.
Definition: CPPGlobals.h:78
double getDepth()
Definition: GeoTessPosition.h:1017
int getTriangle()
Definition: GeoTessPosition.h:822
const vector< int > & getVertices() const
Definition: GeoTessPosition.h:837
GeoTessModel * getModel()
Definition: GeoTessPosition.h:1023
GeoTessPointMap * getPointMap()
Definition: GeoTessModel.h:1051
Abstract class that manages the radii and data values that span a single layer associated with a sing...
Definition: GeoTessProfile.h:96
Top level class that manages the GeoTessMetaData, GeoTessGrid and Data that comprise a 3D Earth model...
Definition: GeoTessModel.h:119
int getTessID()
Definition: GeoTessPosition.h:1028
bool isNotReferenced()
Definition: GeoTessPosition.h:1239
void setRadiusOutOfRangeAllowed(bool allowed)
Definition: GeoTessPosition.h:1196
double getDepthBottom()
Definition: GeoTessPosition.h:957
double getRadiusBottom()
Definition: GeoTessPosition.h:939
void set(double lat, double lon, double depth)
Definition: GeoTessPosition.h:559
double getLayerThickness(int layid)
Definition: GeoTessPosition.h:990
#define GEOTESS_EXP_IMP
Definition: CPPGlobals.h:71
void getWeights(map< int, double > &weights, double dkm)
Definition: GeoTessPosition.h:1127
void setRadius(double rad)
Definition: GeoTessPosition.h:692
void getCoefficients(map< int, double > &coefficients)
Definition: GeoTessPosition.h:871
double getRadiusTop()
Definition: GeoTessPosition.h:931
int getLayerId()
Definition: GeoTessPosition.h:1056
int getMaxTessLevel(int layid)
Definition: GeoTessPosition.h:905
void setMaxTessLevel(int layid, int maxTess)
Definition: GeoTessPosition.h:889