1 // Copyright (C) 2014-2015 EDF-R&D
2 // This library is free software; you can redistribute it and/or
3 // modify it under the terms of the GNU Lesser General Public
4 // License as published by the Free Software Foundation; either
5 // version 2.1 of the License, or (at your option) any later version.
7 // This library is distributed in the hope that it will be useful,
8 // but WITHOUT ANY WARRANTY; without even the implied warranty of
9 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
10 // Lesser General Public License for more details.
12 // You should have received a copy of the GNU Lesser General Public
13 // License along with this library; if not, write to the Free Software
14 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
16 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
19 #ifndef HYDROData_DTM_HeaderFile
20 #define HYDROData_DTM_HeaderFile
22 #include "HYDROData_Bathymetry.h"
24 #include <gp_Pnt2d.hxx>
26 class Handle_HYDROData_Profile;
27 class Handle_Geom2d_BSplineCurve;
28 class Handle_Geom2d_Curve;
33 DEFINE_STANDARD_HANDLE( HYDROData_DTM, HYDROData_Bathymetry )
35 /**\class HYDROData_DTM
36 * \brief Class that represents the Digital Terrain Model
38 class HYDROData_DTM : public HYDROData_Bathymetry
42 * Enumeration of tags corresponding to the persistent object parameters.
46 DataTag_First = HYDROData_Bathymetry::DataTag_First + 100, ///< first tag, to reserve
53 DEFINE_STANDARD_RTTI( HYDROData_DTM );
55 HYDRODATA_EXPORT HYDROData_SequenceOfObjects GetProfiles() const;
56 HYDRODATA_EXPORT void SetProfiles( const HYDROData_SequenceOfObjects& );
58 HYDRODATA_EXPORT double GetDDZ() const;
59 HYDRODATA_EXPORT void SetDDZ( double );
61 HYDRODATA_EXPORT double GetSpatialStep() const;
62 HYDRODATA_EXPORT void SetSpatialStep( double );
64 HYDRODATA_EXPORT virtual void Update();
69 PointUZ( double u=0, double z=0 ) { U=u; Z=z; }
73 class CurveUZ : public std::vector<PointUZ>
76 CurveUZ( double theXcurv );
81 CurveUZ operator + ( const CurveUZ& ) const;
82 CurveUZ operator * ( double ) const;
90 void reserve( int theNbPoints );
91 void push_back( const gp_Pnt& thePnt, const gp_Dir& theTangent );
93 TopoDS_Edge createEdge( bool is3d ) const;
96 std::vector<gp_Pnt> myPoints;
97 std::vector<gp_Dir> myDirs;
101 friend class HYDROData_Iterator;
102 friend class test_HYDROData_DTM;
104 HYDRODATA_EXPORT HYDROData_DTM();
105 virtual HYDRODATA_EXPORT ~HYDROData_DTM();
107 static Handle_Geom2d_BSplineCurve CreateHydraulicAxis(
108 const std::vector<Handle_HYDROData_Profile>& theProfiles,
109 std::vector<double>& theDistances );
111 static std::vector<Handle_Geom2d_Curve> ProfileToParametric( const Handle_HYDROData_Profile& theProfile,
112 double& theUMin, double& theUMax );
114 static void GetProperties( const Handle_HYDROData_Profile& theProfile,
115 gp_Pnt& theLowestPoint, gp_Vec2d& theDir,
117 double& theZMin, double& theZMax );
119 static void ProfileDiscretization( const Handle_HYDROData_Profile& theProfile,
120 double theXCurv, double theMinZ, double theMaxZ, double theDDZ,
121 CurveUZ& theMidPointCurve,
122 CurveUZ& theWidthCurve,
123 double theTolerance = 1E-6 );
125 static void CurveTo3D( const Handle_Geom2d_BSplineCurve& theHydraulicAxis,
126 const CurveUZ& theMidCurve, const CurveUZ& theWidthCurve,
127 AltitudePoints& thePoints,
128 Bank* theLeftBank = 0,
129 Bank* theRightBank = 0,
132 static void Interpolate( const CurveUZ& theCurveA, const CurveUZ& theCurveB,
133 int theNbSteps, std::vector<CurveUZ>& theInterpolation,
136 static AltitudePoints Interpolate( const Handle_Geom2d_BSplineCurve& theHydraulicAxis,
137 const Handle_HYDROData_Profile& theProfileA,
139 const Handle_HYDROData_Profile& theProfileB,
141 double theDDZ, int theNbSteps, bool isAddSecond,
142 Bank* theLeftBank = 0,
143 Bank* theRightBank = 0 );
145 static AltitudePoints Interpolate( const std::vector<Handle_HYDROData_Profile>& theProfiles,
146 double theDDZ, double theSpatialStep,
147 Bank* theLeftBank = 0,
148 Bank* theRightBank = 0 );