GeoTessCPP  2.2.3
Software to facilitate storage and retrieval of 3D information about the Earth.
GeoTessPositionLinear.h
Go to the documentation of this file.
1 //- ****************************************************************************
2 //-
3 //- Copyright 2009 National Technology & Engineering Solutions of Sandia, LLC
4 //- (NTESS). Under the terms of Contract DE-NA0003525 with NTESS, the U.S.
5 //- Government retains certain rights in this software.
6 //-
7 //- BSD Open Source License
8 //- All rights reserved.
9 //-
10 //- Redistribution and use in source and binary forms, with or without
11 //- modification, are permitted provided that the following conditions are met:
12 //-
13 //- 1. Redistributions of source code must retain the above copyright notice,
14 //- this list of conditions and the following disclaimer.
15 //-
16 //- 2. Redistributions in binary form must reproduce the above copyright
17 //- notice, this list of conditions and the following disclaimer in the
18 //- documentation and/or other materials provided with the distribution.
19 //-
20 //- 3. Neither the name of the copyright holder nor the names of its
21 //- contributors may be used to endorse or promote products derived from
22 //- this software without specific prior written permission.
23 //-
24 //- THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
25 //- AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
26 //- IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
27 //- ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
28 //- LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
29 //- CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
30 //- SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
31 //- INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
32 //- CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
33 //- ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
34 //- POSSIBILITY OF SUCH DAMAGE.
35 //-
36 //- ****************************************************************************
37 
38 #ifndef GEOTESSPOSITIONLINEAR_OBJECT_H
39 #define GEOTESSPOSITIONLINEAR_OBJECT_H
40 
41 // **** _SYSTEM INCLUDES_ ******************************************************
42 
43 #include <iostream>
44 #include <string>
45 #include <fstream>
46 #include <vector>
47 
48 // use standard library objects
49 using namespace std;
50 
51 // **** _LOCAL INCLUDES_ *******************************************************
52 
53 #include "CPPUtils.h"
54 #include "GeoTessUtils.h"
55 #include "GeoTessPosition.h"
56 #include "GeoTessGrid.h"
58 
59 // **** _BEGIN GEOTESS NAMESPACE_ **********************************************
60 
61 namespace geotess
62 {
63 
64 // **** _FORWARD REFERENCES_ ***************************************************
65 
66 class GeoTessModel;
67 
68 // **** _CLASS DEFINITION_ *****************************************************
69 
79 {
80 protected:
81 
88  virtual void update2D(int tid)
89  {
90  vector<int>& vt = vertices[tid];
91  vector<double>& hc = hCoefficients[tid];
92 
93  vt.clear();
94  hc.clear();
95 
96  const int* trngl = grid.getTriangleVertexIndexes(getTriangle(tid));
97 
98  // if the interpolation point falls on a grid node set the list of interpolation
99  // vertices to include only the identified vertex, and the interpolation coefficient = 1.
100  if (GeoTessUtils::dot(unitVector, grid.getVertex(trngl[0])) > cos(1e-7))
101  { vt.push_back(trngl[0]); hc.push_back(1.0); return; }
102 
103  if (GeoTessUtils::dot(unitVector, grid.getVertex(trngl[1])) > cos(1e-7))
104  { vt.push_back(trngl[1]); hc.push_back(1.0); return; }
105 
106  if (GeoTessUtils::dot(unitVector, grid.getVertex(trngl[2])) > cos(1e-7))
107  { vt.push_back(trngl[2]); hc.push_back(1.0); return; }
108 
109  // interpolation point is not coincident with a grid node.
110  vector<double>& lc = linearCoefficients[tid];
111  vt.push_back(trngl[0]);
112  hc.push_back(lc[0]);
113  vt.push_back(trngl[1]);
114  hc.push_back(lc[1]);
115  vt.push_back(trngl[2]);
116  hc.push_back(lc[2]);
117 }
118 
119 public:
120 
125  const GeoTessInterpolatorType& radialType);
126 
131 
136  virtual const GeoTessInterpolatorType& getInterpolatorType() const { return GeoTessInterpolatorType::LINEAR; }
137 
138  virtual long getMemory() { return (long) sizeof(GeoTessPositionLinear) + memory(); }
139 
140 };
141 // end class GeoTessModelLinear
142 
143 }// end namespace geotess
144 
145 #endif // GEOTESSPOSITIONLINEAR_OBJECT_H
geotess
Definition: ArrayReuse.h:57
geotess::GeoTessPositionLinear::update2D
virtual void update2D(int tid)
Definition: GeoTessPositionLinear.h:88
geotess::GeoTessPositionLinear::getInterpolatorType
virtual const GeoTessInterpolatorType & getInterpolatorType() const
Definition: GeoTessPositionLinear.h:136
GeoTessUtils.h
geotess::GeoTessPositionLinear::getMemory
virtual long getMemory()
Definition: GeoTessPositionLinear.h:138
geotess::GeoTessPositionLinear::~GeoTessPositionLinear
virtual ~GeoTessPositionLinear()
GeoTessInterpolatorType.h
GeoTessGrid.h
GEOTESS_EXP_IMP
#define GEOTESS_EXP_IMP
Definition: CPPGlobals.h:73
geotess::GeoTessPositionLinear::GeoTessPositionLinear
GeoTessPositionLinear(GeoTessModel *model, const GeoTessInterpolatorType &radialType)
GeoTessPosition.h
geotess::GeoTessPositionLinear
Implements linear interpolation in geographic dimensions of a grid.
Definition: GeoTessPositionLinear.h:79
geotess::GeoTessModel
Top level class that manages the GeoTessMetaData, GeoTessGrid and GeoTessData that comprise a 3D Eart...
Definition: GeoTessModel.h:123
geotess::GeoTessPosition
Information about an interpolated point at an arbitrary position in a model.
Definition: GeoTessPosition.h:104
CPPUtils.h
geotess::GeoTessInterpolatorType
Enumeration of the interpolation algorithms supported by GeoTess including LINEAR,...
Definition: GeoTessInterpolatorType.h:74