Salome HOME
refs #430: incorrect coordinates in dump polyline
[modules/hydro.git] / src / HYDROData / HYDROData_ObstacleAltitude.cxx
1
2 #include "HYDROData_ObstacleAltitude.h"
3
4 #include "HYDROData_Document.h"
5 #include "HYDROData_Object.h"
6 #include "HYDROData_Projection.h"
7
8 #include <Precision.hxx>
9
10 #include <TopoDS_Shape.hxx>
11
12 IMPLEMENT_STANDARD_HANDLE(HYDROData_ObstacleAltitude, HYDROData_IAltitudeObject)
13 IMPLEMENT_STANDARD_RTTIEXT(HYDROData_ObstacleAltitude, HYDROData_IAltitudeObject)
14
15 HYDROData_ObstacleAltitude::HYDROData_ObstacleAltitude()
16 : HYDROData_IAltitudeObject()
17 {
18 }
19
20 HYDROData_ObstacleAltitude::~HYDROData_ObstacleAltitude()
21 {
22 }
23
24 double HYDROData_ObstacleAltitude::GetAltitudeForPoint( const gp_XY& thePoint ) const
25 {
26   double aResAltitude = GetInvalidAltitude();
27
28   Handle(HYDROData_Object) anObject =
29     Handle(HYDROData_Object)::DownCast( GetFatherObject() );
30   if ( anObject.IsNull() )
31     return aResAltitude;
32
33   TopoDS_Shape anObjectShape3D = anObject->GetShape3D();
34   if ( anObjectShape3D.IsNull() )
35     return aResAltitude;
36
37   HYDROData_Make3dMesh aMesher3D( anObjectShape3D, Precision::Intersection() );
38
39   gp_Pnt aHighestPoint;
40   if ( aMesher3D.GetHighestOriginal( thePoint.X(), thePoint.Y(), aHighestPoint ) )
41     aResAltitude = aHighestPoint.Z();
42
43   return aResAltitude;
44 }
45
46
47
48