GeoTessCPP  2.0.0
Software to facilitate storage and retrieval of 3D information about the Earth.
 All Classes Namespaces Files Functions Variables Typedefs Friends Defines
include/HorizonLayer.h
Go to the documentation of this file.
00001 //- ****************************************************************************
00002 //- 
00003 //- Copyright 2009 Sandia Corporation. Under the terms of Contract
00004 //- DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government
00005 //- retains certain rights in this software.
00006 //- 
00007 //- BSD Open Source License.
00008 //- All rights reserved.
00009 //- 
00010 //- Redistribution and use in source and binary forms, with or without
00011 //- modification, are permitted provided that the following conditions are met:
00012 //- 
00013 //-    * Redistributions of source code must retain the above copyright notice,
00014 //-      this list of conditions and the following disclaimer.
00015 //-    * Redistributions in binary form must reproduce the above copyright
00016 //-      notice, this list of conditions and the following disclaimer in the
00017 //-      documentation and/or other materials provided with the distribution.
00018 //-    * Neither the name of Sandia National Laboratories nor the names of its
00019 //-      contributors may be used to endorse or promote products derived from
00020 //-      this software without specific prior written permission.
00021 //- 
00022 //- THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
00023 //- AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
00024 //- IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
00025 //- ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE
00026 //- LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
00027 //- CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
00028 //- SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
00029 //- INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
00030 //- CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
00031 //- ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
00032 //- POSSIBILITY OF SUCH DAMAGE.
00033 //-
00034 //- ****************************************************************************
00035 
00036 #ifndef HORIZONLAYER_H_
00037 #define HORIZONLAYER_H_
00038 
00039 // **** _SYSTEM INCLUDES_ ******************************************************
00040 
00041 #include <iostream>
00042 #include <fstream>
00043 #include <string>
00044 #include <climits>
00045 
00046 // use standard library objects
00047 using namespace std;
00048 
00049 // **** _LOCAL INCLUDES_ *******************************************************
00050 
00051 #include "Horizon.h"
00052 #include "Profile.h"
00053 #include "GeoTessPosition.h"
00054 
00055 // **** _BEGIN GEOTESS NAMESPACE_ *********************************************
00056 
00057 namespace geotess {
00058 
00059 // **** _FORWARD REFERENCES_ ***************************************************
00060 
00061 // **** _CLASS DEFINITION_ *****************************************************
00062 
00073 class GEOTESS_EXP_IMP HorizonLayer : public Horizon
00074 {
00075 
00076 private:
00077 
00082         double fraction;
00083 
00084 public:
00085 
00096         HorizonLayer(const double& fractionalPosition, const int& lyrIndex)
00097         : Horizon(lyrIndex), fraction(fractionalPosition)
00098         {};
00099 
00103         HorizonLayer(HorizonLayer& other) : Horizon(other.getLayerIndex()), fraction(other.fraction)
00104         {
00105         }
00106 
00110         HorizonLayer& operator=(const HorizonLayer& other)
00111         {
00112                 layerIndex = other.layerIndex;
00113                 fraction = other.fraction;
00114                 return *this;
00115         }
00116 
00120         virtual ~HorizonLayer() {};
00121 
00122         virtual double getValue() { return fraction; };
00123 
00124         virtual double getRadius(const double* position, Profile** profiles)
00125         {
00126                 double bottom = profiles[layerIndex]->getRadiusBottom();
00127                 if (fraction <= 0. )
00128                         return bottom;
00129                 double top = profiles[layerIndex]->getRadiusTop();
00130                 if (fraction >= 1.)
00131                         return top;
00132                 return bottom + fraction*(top-bottom);
00133         }
00134 
00135         virtual string class_name() { return "HorizonLayer"; };
00136 
00137         virtual double getRadius(GeoTessPosition& position)
00138         {
00139                 double bottom = position.getRadiusBottom(layerIndex);
00140                 if (fraction <= 0. )
00141                         return bottom;
00142                 double top = position.getRadiusTop(layerIndex);
00143                 if (fraction >= 1.)
00144                         return top;
00145                 return bottom + fraction*(top-bottom);
00146         }
00147 
00148         virtual string str()
00149         {
00150                 string s = "layer " + CPPUtils::dtos(fraction) + " " + CPPUtils::itos(layerIndex);
00151                 return s;
00152         }
00153 
00154 }; // end class DataValue
00155 
00156 } // end namespace geotess
00157 
00158 #endif /* HORIZONLAYER_H_ */