Comments moved.

[modules/hydro.git] / src / HYDROData / HYDROData_CalculationCase.cxx

diff --git a/src/HYDROData/HYDROData_CalculationCase.cxx b/src/HYDROData/HYDROData_CalculationCase.cxx
index 9d6b58f460957305786c6056ca4d85439d3c886c..5a5b3ad2dc78fe86434a6d11b93dcf041520135c 100644 (file)
--- a/src/HYDROData/HYDROData_CalculationCase.cxx
+++ b/src/HYDROData/HYDROData_CalculationCase.cxx
@@ -2,6 +2,7 @@
 #include "HYDROData_CalculationCase.h"
 
 #include "HYDROData_ArtificialObject.h"
+#include "HYDROData_Bathymetry.h"
 #include "HYDROData_Document.h"
 #include "HYDROData_Iterator.h"
 #include "HYDROData_NaturalObject.h"
@@ -19,7 +20,8 @@
 #include <TopExp.hxx>
 #include <TopTools_ListOfShape.hxx>
 #include <TopTools_ListIteratorOfListOfShape.hxx>
-
+#include <BRepTopAdaptor_FClass2d.hxx>
+#include <TopAbs.hxx>
 #define CALCULATION_REGIONS_PREF GetName() + "_Reg"
 #define CALCULATION_ZONES_PREF GetName() + "_Zone"
 
@@ -518,4 +520,126 @@ TopoDS_Shell HYDROData_CalculationCase::GetShell()
   return aShell;
 }
 
+double HYDROData_CalculationCase::GetAltitudeForPoint( const gp_XY& thePoint ) const
+{
+  double aResAltitude = HYDROData_Bathymetry::GetInvalidAltitude();
+
+  Handle(HYDROData_Zone) aZone = GetZoneFromPoint( thePoint );
+  if ( aZone.IsNull() )
+    return aResAltitude;
+
+  HYDROData_Zone::MergeBathymetriesType aZoneMergeType = aZone->GetMergeType();
+  if ( !aZone->IsMergingNeed() )
+  {
+    aZoneMergeType = HYDROData_Zone::Merge_UNKNOWN;
+  }
+  else if ( aZoneMergeType == HYDROData_Zone::Merge_UNKNOWN )
+  {
+    return aResAltitude;
+  }
+
+  if ( aZoneMergeType == HYDROData_Zone::Merge_Object )
+  {
+    Handle(HYDROData_Bathymetry) aMergeBathymetry = aZone->GetMergeBathymetry();
+    if ( !aMergeBathymetry.IsNull() )
+      aResAltitude = aMergeBathymetry->GetAltitudeForPoint( thePoint );
+  }
+  else
+  {
+    HYDROData_SequenceOfObjects aZoneObjects = aZone->GetGeometryObjects();
+    HYDROData_SequenceOfObjects::Iterator anIter( aZoneObjects );
+    for ( ; anIter.More(); anIter.Next() )
+    {
+      Handle(HYDROData_Object) aZoneObj =
+        Handle(HYDROData_Object)::DownCast( anIter.Value() );
+      if ( aZoneObj.IsNull() )
+        continue;
+
+      Handle(HYDROData_Bathymetry) anObjBathymetry = aZoneObj->GetBathymetry();
+      if ( anObjBathymetry.IsNull() )
+        continue;
 
+      double aPointAltitude = anObjBathymetry->GetAltitudeForPoint( thePoint );
+      if ( ValuesEquals( aPointAltitude, HYDROData_Bathymetry::GetInvalidAltitude() ) )
+        continue;
+
+      if ( aZoneMergeType == HYDROData_Zone::Merge_UNKNOWN )
+      {
+        aResAltitude = aPointAltitude;
+        break;
+      }
+      else if ( aZoneMergeType == HYDROData_Zone::Merge_ZMIN )
+      {
+        if ( ValuesEquals( aResAltitude, HYDROData_Bathymetry::GetInvalidAltitude() ) ||
+             aResAltitude > aPointAltitude )
+        {
+          aResAltitude = aPointAltitude;
+        }
+      }
+      else if ( aZoneMergeType == HYDROData_Zone::Merge_ZMAX )
+      {
+        if ( ValuesEquals( aResAltitude, HYDROData_Bathymetry::GetInvalidAltitude() ) ||
+             aResAltitude < aPointAltitude )
+        {
+          aResAltitude = aPointAltitude;
+        }
+      }
+    }
+  }
+
+  return aResAltitude;
+}
+
+Handle(HYDROData_Zone) HYDROData_CalculationCase::GetZoneFromPoint( const gp_XY& thePoint ) const
+{
+  Handle(HYDROData_Zone) aResZone;
+
+  HYDROData_SequenceOfObjects aRegions = GetRegions();
+
+  HYDROData_SequenceOfObjects::Iterator anIter( aRegions );
+  for ( ; anIter.More() && aResZone.IsNull(); anIter.Next() )
+  {
+    Handle(HYDROData_Region) aRegion =
+      Handle(HYDROData_Region)::DownCast( anIter.Value() );
+    if ( aRegion.IsNull() )
+      continue;
+
+    HYDROData_SequenceOfObjects aZones = aRegion->GetZones();
+    HYDROData_SequenceOfObjects::Iterator aZonesIter( aZones );
+    for ( ; aZonesIter.More() && aResZone.IsNull(); aZonesIter.Next() )
+    {
+      Handle(HYDROData_Zone) aRegZone =
+        Handle(HYDROData_Zone)::DownCast( aZonesIter.Value() );
+      if ( aRegZone.IsNull() )
+        continue;
+
+      PointClassification aPointRelation = GetPointClassification( thePoint, aRegZone );
+      if ( aPointRelation != POINT_OUT )
+        aResZone = aRegZone; // We found the desired zone
+    }
+  }
+
+  return aResZone;
+}
+HYDROData_CalculationCase::PointClassification HYDROData_CalculationCase::GetPointClassification(
+  const gp_XY&                  thePoint,
+  const Handle(HYDROData_Zone)& theZone ) const
+{
+  PointClassification aRes = POINT_OUT;
+  if ( theZone.IsNull() )
+    return aRes;
+
+  TopoDS_Face aZoneFace = TopoDS::Face( theZone->GetShape() );
+  if ( aZoneFace.IsNull() )
+    return aRes;
+
+  BRepTopAdaptor_FClass2d aClassifier( aZoneFace, Precision::Confusion() );
+  TopAbs_State State = aClassifier.Perform( gp_Pnt2d(thePoint), Standard_False );
+  if (State == TopAbs_OUT)
+       aRes =  POINT_OUT;
+  else if(State == TopAbs_IN)
+    aRes =  POINT_IN;
+  else if(State == TopAbs_ON)
+    aRes =  POINT_ON;    
+  return aRes;
+}