src/HYDRO_tests/test_HYDROData_LandCoverMap.cxx

   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.
   6 //
   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.
  11 //
  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
  15 //
  16 // See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
  17 //
  18
  19 #include <test_HYDROData_LandCoverMap.h>
  20 #include <HYDROData_CalculationCase.h>
  21 #include <HYDROData_Document.h>
  22 #include <HYDROData_LandCoverMap.h>
  23 #include <HYDROData_ImmersibleZone.h>
  24 #include <HYDROData_PolylineXY.h>
  25 #include <HYDROData_StricklerTable.h>
  26 #include <HYDROData_Tool.h>
  27 #include <HYDROData_ShapeFile.h>
  28 #include <HYDROGUI_LandCoverMapPrs.h>
  29 #include <TopoDS_Edge.hxx>
  30 #include <TopoDS_Wire.hxx>
  31 #include <TopoDS_Face.hxx>
  32 #include <TopoDS.hxx>
  33 #include <TestViewer.h>
  34 #include <TestShape.h>
  35 #include <TopTools_ListOfShape.hxx>
  36 #include <AIS_DisplayMode.hxx>
  37 #include <Aspect_ColorScale.hxx>
  38 #include <BRepMesh_IncrementalMesh.hxx>
  39 #include <QString>
  40 #include <QColor>
  41 #include <QMap>
  42 #include <QDir>
  43 #include <BRepTools.hxx>
  44 #include <BRep_Builder.hxx>
  45 #include <BRepCheck_Analyzer.hxx>
  46
  47 const QString REF_DATA_PATH = qgetenv( "HYDRO_REFERENCE_DATA" );
  48 const QString DEF_STR_PATH = qgetenv( "HYDRO_SRC_DIR" ) + "/src/HYDROGUI/resources/def_strickler_table.txt";
  49
  50 void test_HYDROData_LandCoverMap::test_add_2_objects()
  51 {
  52   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
  53
  54   Handle(HYDROData_LandCoverMap) aMap =
  55     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
  56
  57   CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() );
  58
  59   TopoDS_Face aLC1 = Face( QList<double>() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 );
  60   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) );
  61
  62   TopoDS_Face aLC2 = Face( QList<double>() << 30 << 20 << 60 << 10 << 70 << 35 << 40 << 40 );
  63   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "test2" ) );
  64
  65   TestViewer::show( aMap->GetShape(), AIS_Shaded, true, "LandCoverMap_Add_2_Objects" );
  66   CPPUNIT_ASSERT_IMAGES
  67
  68   HYDROData_LandCoverMap::Explorer anIt( aMap );
  69   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
  70   CPPUNIT_ASSERT_EQUAL( QString( "test1" ), anIt.StricklerType() );
  71   anIt.Next();
  72   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
  73   CPPUNIT_ASSERT_EQUAL( QString( "test2" ), anIt.StricklerType() );
  74   anIt.Next();
  75   CPPUNIT_ASSERT_EQUAL( false, anIt.More() );
  76
  77   aDoc->Close();
  78 }
  79
  80 void test_HYDROData_LandCoverMap::test_split()
  81 {
  82   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
  83
  84   Handle(HYDROData_LandCoverMap) aMap =
  85     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
  86
  87   CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() );
  88
  89   TopoDS_Face aLC = Face( QList<double>() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 );
  90   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC, "test1" ) );
  91
  92   Handle(HYDROData_PolylineXY) aPolyline =
  93     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
  94   TopoDS_Wire aWire = Wire( QList<double>() << 10 << 40 << 30 << 10 << 40 << 10 << 60 << 10, false );
  95   aPolyline->SetShape( aWire );
  96
  97   CPPUNIT_ASSERT_EQUAL( true, aMap->Split( aPolyline ) );
  98
  99   TestViewer::show( aMap->GetShape(), AIS_Shaded, true, "LandCoverMap_Split_1" );
 100   TestViewer::show( aWire, 0, true, Qt::blue );
 101   CPPUNIT_ASSERT_IMAGES
 102
 103   HYDROData_LandCoverMap::Explorer anIt( aMap );
 104   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
 105   CPPUNIT_ASSERT_EQUAL( QString( "test1" ), anIt.StricklerType() );
 106   anIt.Next();
 107   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
 108   CPPUNIT_ASSERT_EQUAL( QString( "test1" ), anIt.StricklerType() );
 109   anIt.Next();
 110   CPPUNIT_ASSERT_EQUAL( false, anIt.More() );
 111
 112   aDoc->Close();
 113 }
 114
 115 void test_HYDROData_LandCoverMap::test_incomplete_split()
 116 {
 117   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 118
 119   Handle(HYDROData_LandCoverMap) aMap =
 120     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 121
 122   CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() );
 123
 124   TopoDS_Face aLC = Face( QList<double>() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 );
 125   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC, "test1" ) );
 126
 127   Handle(HYDROData_PolylineXY) aPolyline =
 128     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 129   TopoDS_Wire aWire = Wire( QList<double>() << 10 << 40 << 30 << 10 << 40 << 10, false );
 130   aPolyline->SetShape( aWire );
 131
 132   CPPUNIT_ASSERT_EQUAL( true, aMap->Split( aPolyline ) );
 133
 134   TestViewer::show( aMap->GetShape(), AIS_Shaded, true, "LandCoverMap_Split_2" );
 135   TestViewer::show( aWire, 0, true, Qt::green );
 136   CPPUNIT_ASSERT_IMAGES
 137
 138   HYDROData_LandCoverMap::Explorer anIt( aMap );
 139   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
 140   CPPUNIT_ASSERT_EQUAL( QString( "test1" ), anIt.StricklerType() );
 141   anIt.Next();
 142   CPPUNIT_ASSERT_EQUAL( false, anIt.More() );
 143
 144   aDoc->Close();
 145 }
 146
 147 void test_HYDROData_LandCoverMap::test_merge()
 148 {
 149   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 150
 151   Handle(HYDROData_LandCoverMap) aMap =
 152     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 153
 154   CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() );
 155
 156   TopoDS_Face aLC1 = Face( QList<double>() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 );
 157
 158   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) );
 159
 160   TopoDS_Face aLC2 = Face( QList<double>() << 21 << 34 << 24 << 25 << 37   << 37 << 40  << 61 <<
 161                                               44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 <<
 162                                              128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 <<
 163                                               31 << 114 );
 164   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "test2" ) );
 165
 166   TopoDS_Face aLC3 = Face( QList<double>() << 4 << 54 << 1   << 47 << 51  << 45 <<
 167                                             127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 );
 168   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "test3" ) );
 169
 170   //TestViewer::show( aMap->GetShape(), AIS_Shaded, true );
 171   //TestViewer::AssertEqual( "LandCoverMap_Before_Merge" );
 172
 173   //BRepTools::Write( aMap->GetShape(), "c:\\map.brep" );
 174
 175   QString aType1, aType2;
 176   gp_Pnt2d aPnt1( 25, 35 );
 177   gp_Pnt2d aPnt2( 45, 55 );
 178   //gp_Pnt2d aPnt3( 45, 10 );
 179   //gp_Pnt2d aPnt4( 45, 80 );
 180   TopTools_ListOfShape aList;
 181   aList.Append( aMap->FindByPoint( aPnt1, aType1 ) );
 182   aList.Append( aMap->FindByPoint( aPnt2, aType2 ) );
 183   CPPUNIT_ASSERT_EQUAL( true, aMap->Merge( aList, "new" ) );
 184
 185   TestViewer::show( aMap->GetShape(), AIS_Shaded, true, "LandCoverMap_Merge_1" );
 186   //TestViewer::show( BRepBuilderAPI_MakeEdge( gp_Pnt(aPnt1.X(), aPnt1.Y(), 0), gp_Pnt(aPnt2.X(), aPnt2.Y(), 0) ).Edge(), QColor( Qt::blue ), AIS_Shaded );
 187   CPPUNIT_ASSERT_IMAGES
 188
 189   aDoc->Close();
 190 }
 191
 192 void test_HYDROData_LandCoverMap::test_merge_refs_691()
 193 {
 194   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 195
 196   Handle(HYDROData_PolylineXY) Polyline_4 =
 197     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 198   Polyline_4->SetName( "test_Polyline_4" );
 199   Polyline_4->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 200   Polyline_4->AddPoint( 0, gp_XY( -242.741935, -16.129032 ) );
 201   Polyline_4->AddPoint( 0, gp_XY( -339.516129, -95.161290 ) );
 202   Polyline_4->AddPoint( 0, gp_XY( -204.032258, -167.740000 ) );
 203   Polyline_4->AddPoint( 0, gp_XY( -86.290000, -158.060000 ) );
 204   Polyline_4->AddPoint( 0, gp_XY( -83.064516, -37.096774 ) );
 205   Polyline_4->AddPoint( 0, gp_XY( -116.935484, 17.741935 ) );
 206   Polyline_4->AddSection( "Section_2", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 207   Polyline_4->AddPoint( 1, gp_XY( -239.516129, -295.161290 ) );
 208   Polyline_4->AddPoint( 1, gp_XY( -204.030000, -167.740000 ) );
 209   Polyline_4->AddPoint( 1, gp_XY( -86.290000, -158.060000 ) );
 210   Polyline_4->AddPoint( 1, gp_XY( -2.419355, -288.709677 ) );
 211   Polyline_4->AddPoint( 1, gp_XY( -165.322581, -351.612903 ) );
 212   Polyline_4->Update();
 213
 214   Handle(HYDROData_LandCoverMap) aMap =
 215     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 216
 217   CPPUNIT_ASSERT_EQUAL( true, aMap->Add( Polyline_4, "test" ) );
 218
 219   TopTools_ListOfShape aList;
 220   HYDROData_LandCoverMap::Explorer anIt( aMap );
 221   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
 222   CPPUNIT_ASSERT_EQUAL( QString( "test" ), anIt.StricklerType() );
 223   aList.Append( anIt.Face() );
 224   anIt.Next();
 225   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
 226   CPPUNIT_ASSERT_EQUAL( QString( "test" ), anIt.StricklerType() );
 227   aList.Append( anIt.Face() );
 228   anIt.Next();
 229   CPPUNIT_ASSERT_EQUAL( false, anIt.More() );
 230
 231   BRepCheck_Analyzer aCheck1( aMap->GetShape() );
 232   CPPUNIT_ASSERT_EQUAL( true, (bool)aCheck1.IsValid() );
 233
 234   CPPUNIT_ASSERT_EQUAL( true, aMap->Merge( aList, "test" ) );
 235
 236   BRepCheck_Analyzer aCheck2( aMap->GetShape() );
 237   CPPUNIT_ASSERT_EQUAL( true, (bool)aCheck2.IsValid() );
 238
 239   // we try construct the triangulation to check that it is correct
 240   BRepMesh_IncrementalMesh anIncMesh( aMap->GetShape(), 0.01 );
 241   CPPUNIT_ASSERT_EQUAL( true, (bool)anIncMesh.IsDone() );
 242   BRepCheck_Analyzer aCheck3( aMap->GetShape() );
 243   CPPUNIT_ASSERT_EQUAL( true, (bool)aCheck3.IsValid() );
 244
 245   Handle(HYDROGUI_LandCoverMapPrs) aPrs = new HYDROGUI_LandCoverMapPrs( aMap );
 246   TestViewer::show( aPrs, AIS_Shaded, 0, true, "LandCoverMap_Merge_2" );
 247   CPPUNIT_ASSERT_IMAGES
 248
 249   aDoc->Close();
 250 }
 251
 252 void test_HYDROData_LandCoverMap::test_remove()
 253 {
 254   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 255
 256   Handle(HYDROData_LandCoverMap) aMap =
 257     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 258
 259   CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() );
 260
 261   TopoDS_Face aLC1 = Face( QList<double>() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 );
 262
 263   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) );
 264
 265   TopoDS_Face aLC2 = Face( QList<double>() << 21 << 34 << 24 << 25 << 37   << 37 << 40  << 61 <<
 266                                               44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 <<
 267                                              128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 <<
 268                                               31 << 114 );
 269   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "test2" ) );
 270
 271   TopoDS_Face aLC3 = Face( QList<double>() << 4 << 54 << 1   << 47 << 51  << 45 <<
 272                                             127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 );
 273   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "test3" ) );
 274
 275   QString aType1, aType2;
 276   gp_Pnt2d aPnt1( 25, 35 );
 277   gp_Pnt2d aPnt2( 45, 55 );
 278   TopTools_ListOfShape aList;
 279   aList.Append( aMap->FindByPoint( aPnt1, aType1 ) );
 280   aList.Append( aMap->FindByPoint( aPnt2, aType2 ) );
 281   CPPUNIT_ASSERT_EQUAL( true, aMap->Remove( aList ) );
 282
 283   TestViewer::show( aMap->GetShape(), AIS_Shaded, true, "LandCoverMap_Remove_1" );
 284   CPPUNIT_ASSERT_IMAGES
 285
 286   aDoc->Close();
 287 }
 288
 289 void test_HYDROData_LandCoverMap::test_merge_faces_boxes()
 290 {
 291   TopoDS_Shape pp1, pp2, pp3, pp4;
 292   BRep_Builder BB;
 293   BRepTools::Read(pp1, (REF_DATA_PATH + "/pp1.brep").toStdString().c_str(), BB);
 294   BRepTools::Read(pp2, (REF_DATA_PATH + "/pp2.brep").toStdString().c_str(), BB);
 295   BRepTools::Read(pp3, (REF_DATA_PATH + "/pp3.brep").toStdString().c_str(), BB);
 296   BRepTools::Read(pp4, (REF_DATA_PATH + "/pp4.brep").toStdString().c_str(), BB);
 297
 298   CPPUNIT_ASSERT(!pp1.IsNull());
 299   CPPUNIT_ASSERT(!pp2.IsNull());
 300   CPPUNIT_ASSERT(!pp3.IsNull());
 301   CPPUNIT_ASSERT(!pp4.IsNull());
 302
 303   //Test mergeFaces() func // boxes // USD == true
 304   {
 305     TopTools_ListOfShape Faces;
 306     Faces.Append(pp1);
 307     Faces.Append(pp2);
 308     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true );
 309     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_11" );
 310     CPPUNIT_ASSERT_IMAGES
 311     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_FACE);
 312     BRepCheck_Analyzer aBCA(aMergedFace);
 313     CPPUNIT_ASSERT(aBCA.IsValid());
 314   }
 315   {
 316     TopTools_ListOfShape Faces;
 317     Faces.Append(pp1);
 318     Faces.Append(pp3);
 319     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true );
 320     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_12" );
 321     CPPUNIT_ASSERT_IMAGES
 322     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 323     BRepCheck_Analyzer aBCA(aMergedFace);
 324     CPPUNIT_ASSERT(aBCA.IsValid());
 325   }
 326   {
 327     TopTools_ListOfShape Faces;
 328     Faces.Append(pp1);
 329     Faces.Append(pp2);
 330     Faces.Append(pp4);
 331     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true );
 332     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_13" );
 333     CPPUNIT_ASSERT_IMAGES
 334     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 335     BRepCheck_Analyzer aBCA(aMergedFace);
 336     CPPUNIT_ASSERT(aBCA.IsValid());
 337   }
 338   {
 339     TopTools_ListOfShape Faces;
 340     Faces.Append(pp1);
 341     Faces.Append(pp3);
 342     Faces.Append(pp4);
 343     Faces.Append(pp2);
 344     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true );
 345     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_14" );
 346     CPPUNIT_ASSERT_IMAGES
 347     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_FACE);
 348     BRepCheck_Analyzer aBCA(aMergedFace);
 349     CPPUNIT_ASSERT(aBCA.IsValid());
 350   }
 351   //
 352
 353   //Test mergeFaces() func // boxes // USD == false
 354   {
 355     TopTools_ListOfShape Faces;
 356     Faces.Append(pp1);
 357     Faces.Append(pp2);
 358     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false );
 359     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_21" );
 360     CPPUNIT_ASSERT_IMAGES
 361     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL);
 362     BRepCheck_Analyzer aBCA(aMergedFace);
 363     CPPUNIT_ASSERT(aBCA.IsValid());
 364   }
 365   {
 366     TopTools_ListOfShape Faces;
 367     Faces.Append(pp1);
 368     Faces.Append(pp3);
 369     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false );
 370     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_22" );
 371     CPPUNIT_ASSERT_IMAGES
 372     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 373     BRepCheck_Analyzer aBCA(aMergedFace);
 374     CPPUNIT_ASSERT(aBCA.IsValid());
 375   }
 376   {
 377     TopTools_ListOfShape Faces;
 378     Faces.Append(pp1);
 379     Faces.Append(pp2);
 380     Faces.Append(pp4);
 381     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false );
 382     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_23" );
 383     CPPUNIT_ASSERT_IMAGES
 384     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 385     BRepCheck_Analyzer aBCA(aMergedFace);
 386     CPPUNIT_ASSERT(aBCA.IsValid());
 387   }
 388   {
 389     TopTools_ListOfShape Faces;
 390     Faces.Append(pp1);
 391     Faces.Append(pp3);
 392     Faces.Append(pp4);
 393     Faces.Append(pp2);
 394     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false );
 395     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_24" );
 396     CPPUNIT_ASSERT_IMAGES
 397     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL);
 398     BRepCheck_Analyzer aBCA(aMergedFace);
 399     CPPUNIT_ASSERT(aBCA.IsValid());
 400   }
 401 }
 402
 403
 404 void test_HYDROData_LandCoverMap::test_merge_faces_circles()
 405 {
 406   TopoDS_Shape ff1, ff2, ff3;
 407   BRep_Builder BB;
 408   BRepTools::Read(ff1, (REF_DATA_PATH + "/ff1.brep").toStdString().c_str(), BB);
 409   BRepTools::Read(ff2, (REF_DATA_PATH + "/ff2.brep").toStdString().c_str(), BB);
 410   BRepTools::Read(ff3, (REF_DATA_PATH + "/ff3.brep").toStdString().c_str(), BB);
 411
 412   CPPUNIT_ASSERT(!ff1.IsNull());
 413   CPPUNIT_ASSERT(!ff2.IsNull());
 414   CPPUNIT_ASSERT(!ff3.IsNull());
 415
 416   //Test mergeFaces() func // circles // USD == true
 417   {
 418     TopTools_ListOfShape Faces;
 419     Faces.Append(ff1);
 420     Faces.Append(ff2);
 421     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true );
 422     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_11" );
 423     CPPUNIT_ASSERT_IMAGES
 424     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_FACE);
 425     BRepCheck_Analyzer aBCA(aMergedFace);
 426     CPPUNIT_ASSERT(aBCA.IsValid());
 427   }
 428   {
 429     TopTools_ListOfShape Faces;
 430     Faces.Append(ff1);
 431     Faces.Append(ff3);
 432     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true );
 433     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_12" );
 434     CPPUNIT_ASSERT_IMAGES
 435     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 436     BRepCheck_Analyzer aBCA(aMergedFace);
 437     CPPUNIT_ASSERT(aBCA.IsValid());
 438   }
 439   {
 440     TopTools_ListOfShape Faces;
 441     Faces.Append(ff1);
 442     Faces.Append(ff2);
 443     Faces.Append(ff3);
 444     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true );
 445     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_13" );
 446     CPPUNIT_ASSERT_IMAGES
 447     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 448     BRepCheck_Analyzer aBCA(aMergedFace);
 449     CPPUNIT_ASSERT(aBCA.IsValid());
 450   }
 451   {
 452     TopTools_ListOfShape Faces;
 453     Faces.Append(ff1);
 454     Faces.Append(ff3);
 455     Faces.Append(ff2);
 456     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true );
 457     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_14" );
 458     CPPUNIT_ASSERT_IMAGES
 459     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 460     BRepCheck_Analyzer aBCA(aMergedFace);
 461     CPPUNIT_ASSERT(aBCA.IsValid());
 462   }
 463   //
 464
 465   //Test mergeFaces() func // circles // USD == false
 466   {
 467     TopTools_ListOfShape Faces;
 468     Faces.Append(ff1);
 469     Faces.Append(ff2);
 470     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false );
 471     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_21" );
 472     CPPUNIT_ASSERT_IMAGES
 473     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL);
 474     BRepCheck_Analyzer aBCA(aMergedFace);
 475     CPPUNIT_ASSERT(aBCA.IsValid());
 476   }
 477   {
 478     TopTools_ListOfShape Faces;
 479     Faces.Append(ff1);
 480     Faces.Append(ff3);
 481     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false );
 482     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_22" );
 483     CPPUNIT_ASSERT_IMAGES
 484     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 485     BRepCheck_Analyzer aBCA(aMergedFace);
 486     CPPUNIT_ASSERT(aBCA.IsValid());
 487   }
 488   {
 489     TopTools_ListOfShape Faces;
 490     Faces.Append(ff1);
 491     Faces.Append(ff2);
 492     Faces.Append(ff3);
 493     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false );
 494     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_23" );
 495     CPPUNIT_ASSERT_IMAGES
 496     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 497     BRepCheck_Analyzer aBCA(aMergedFace);
 498     CPPUNIT_ASSERT(aBCA.IsValid());
 499   }
 500   {
 501     TopTools_ListOfShape Faces;
 502     Faces.Append(ff1);
 503     Faces.Append(ff3);
 504     Faces.Append(ff2);
 505     TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false );
 506     TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_24" );
 507     CPPUNIT_ASSERT_IMAGES
 508     CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_COMPOUND);
 509     BRepCheck_Analyzer aBCA(aMergedFace);
 510     CPPUNIT_ASSERT(aBCA.IsValid());
 511   }
 512 }
 513
 514 void test_HYDROData_LandCoverMap::test_import_dbf()
 515 {
 516   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 517
 518   Handle(HYDROData_LandCoverMap) aMap =
 519     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 520
 521   CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() );
 522   QString aFileName = REF_DATA_PATH + "/t1.dbf";
 523   QStringList aDBFV;
 524   QStringList aST;
 525   aDBFV.append("100");
 526   aDBFV.append("200");
 527   aDBFV.append("300");
 528   aST.append("water");
 529   aST.append("forest");
 530   aST.append("road");
 531   // aST.size() == aDBFV.size()!!
 532   Handle(HYDROData_PolylineXY) aPolyline =
 533     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 534   TopoDS_Wire aWire = Wire( QList<double>() << 10 << 40 << 30 << 10 << 40 << 10 << 60 << 10 );
 535   aPolyline->SetShape( aWire );
 536   for (int i = 0; i < 3; i++)
 537     aMap->Add(aPolyline, "");
 538   TopoDS_Face aLC1 = Face( QList<double>() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 );
 539   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) );
 540
 541   CPPUNIT_ASSERT_EQUAL( true, aMap->Split( aPolyline ) );
 542   QList<int> Inds = QList<int>() << 1 << 2 << 3;
 543   aMap->ImportDBF(aFileName, "TESTFIELD1", aDBFV, aST, Inds);
 544   HYDROData_LandCoverMap::Explorer anIt( aMap );
 545   CPPUNIT_ASSERT_EQUAL( QString( "forest" ), anIt.StricklerType() );
 546   anIt.Next();
 547   CPPUNIT_ASSERT_EQUAL( QString( "road" ), anIt.StricklerType() );
 548
 549   aDoc->Close();
 550 }
 551
 552 void test_HYDROData_LandCoverMap::test_land_cover_prs_by_types()
 553 {
 554   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 555
 556   Handle(HYDROData_StricklerTable) aTable =
 557     Handle(HYDROData_StricklerTable)::DownCast( aDoc->CreateObject( KIND_STRICKLER_TABLE ) );
 558   CPPUNIT_ASSERT_EQUAL( true, aTable->Import( DEF_STR_PATH ) );
 559
 560   Handle(HYDROData_LandCoverMap) aMap =
 561     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 562
 563   TopoDS_Face aLC1 = Face( QList<double>() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 );
 564   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "Zones de champs cultivé à végétation basse" ) );
 565
 566   TopoDS_Face aLC2 = Face( QList<double>() << 21 << 34 << 24 << 25 << 37   << 37 << 40  << 61 <<
 567                                               44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 <<
 568                                              128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 <<
 569                                               31 << 114 );
 570   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "Zones de champs cultivé à végétation haute" ) );
 571
 572   TopoDS_Face aLC3 = Face( QList<double>() << 4 << 54 << 1   << 47 << 51  << 45 <<
 573                                             127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 );
 574   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "Zones de champs, prairies, sans cultures" ) );
 575
 576   // build presentation object
 577   Handle(HYDROGUI_LandCoverMapPrs) aPrs = new HYDROGUI_LandCoverMapPrs( aMap );
 578   aPrs->SetTable( aTable );
 579   // show presentation in viewer
 580   TestViewer::show( aPrs, AIS_Shaded, 4, true, "LandCoverMap_PrsByTypes" );
 581   // select one of faces (first)
 582   TestViewer::select( 200, 300 );
 583   CPPUNIT_ASSERT_IMAGES
 584
 585   aDoc->Close();
 586 }
 587
 588 void test_HYDROData_LandCoverMap::test_land_cover_prs_by_coeff()
 589 {
 590   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 591
 592   Handle(HYDROData_StricklerTable) aTable =
 593     Handle(HYDROData_StricklerTable)::DownCast( aDoc->CreateObject( KIND_STRICKLER_TABLE ) );
 594   CPPUNIT_ASSERT_EQUAL( true, aTable->Import( DEF_STR_PATH ) );
 595
 596   Handle(HYDROData_LandCoverMap) aMap =
 597     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 598
 599   TopoDS_Face aLC1 = Face( QList<double>() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 );
 600   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "Zones de champs cultivé à végétation basse" ) );
 601
 602   TopoDS_Face aLC2 = Face( QList<double>() << 21 << 34 << 24 << 25 << 37   << 37 << 40  << 61 <<
 603                                               44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 <<
 604                                              128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 <<
 605                                               31 << 114 );
 606   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "Zones de champs cultivé à végétation haute" ) );
 607
 608   TopoDS_Face aLC3 = Face( QList<double>() << 4 << 54 << 1   << 47 << 51  << 45 <<
 609                                             127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 );
 610   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "Zones de champs, prairies, sans cultures" ) );
 611
 612   Handle(HYDROGUI_LandCoverMapPrs) aPrs = new HYDROGUI_LandCoverMapPrs( aMap );
 613   aPrs->SetTable( aTable );
 614   TestViewer::showColorScale( true );
 615   aPrs->SetColorScale( TestViewer::colorScale() );
 616   TestViewer::show( aPrs, AIS_Shaded, 4, true, "LandCoverMap_PrsByCoeff" );
 617   CPPUNIT_ASSERT_IMAGES
 618
 619   TestViewer::showColorScale( false );
 620   aDoc->Close();
 621 }
 622
 623 void test_HYDROData_LandCoverMap::test_dump_python()
 624 {
 625   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 626
 627   Handle(HYDROData_StricklerTable) aTable =
 628     Handle(HYDROData_StricklerTable)::DownCast( aDoc->CreateObject( KIND_STRICKLER_TABLE ) );
 629   CPPUNIT_ASSERT_EQUAL( true, aTable->Import( DEF_STR_PATH ) );
 630   aTable->SetName( "DefStr" );
 631
 632   Handle(HYDROData_LandCoverMap) aMap =
 633     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 634   aMap->SetName( "test_LCM" );
 635
 636   CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() );
 637
 638   TopoDS_Face aLC1 = Face( QList<double>() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 );
 639
 640   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) );
 641
 642   TopoDS_Face aLC2 = Face( QList<double>() << 21 << 34 << 24 << 25 << 37   << 37 << 40  << 61 <<
 643                                               44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 <<
 644                                              128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 <<
 645                                               31 << 114 );
 646   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "test2" ) );
 647
 648   TopoDS_Face aLC3 = Face( QList<double>() << 4 << 54 << 1   << 47 << 51  << 45 <<
 649                                             127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 );
 650   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "test3" ) );
 651
 652   QString aBaseName = "lc_dump.py";
 653   QString aBaseShp = "lc_dump.shp";
 654   QString aBaseDbf = "lc_dump.dbf";
 655   QString aTmpPath = QDir::tempPath() + "/" + aBaseName;
 656   CPPUNIT_ASSERT_EQUAL( true, aDoc->DumpToPython( aTmpPath, false ) );
 657
 658   //TODO: CPPUNIT_ASSERT_SCRIPTS_EQUAL( aBaseName, true, false, 20 );
 659   //TODO: CPPUNIT_ASSERT_FILES_EQUAL( aBaseShp );
 660   //TODO: CPPUNIT_ASSERT_FILES_EQUAL( aBaseDbf );
 661
 662   aDoc->Close();
 663 }
 664
 665 void test_HYDROData_LandCoverMap::test_transparent_prs()
 666 {
 667   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 668
 669   Handle(HYDROData_StricklerTable) aTable =
 670     Handle(HYDROData_StricklerTable)::DownCast( aDoc->CreateObject( KIND_STRICKLER_TABLE ) );
 671   CPPUNIT_ASSERT_EQUAL( true, aTable->Import( DEF_STR_PATH ) );
 672
 673   Handle(HYDROData_PolylineXY) aPoly =
 674     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 675   aPoly->AddSection( "", HYDROData_PolylineXY::SECTION_SPLINE, true );
 676   aPoly->AddPoint( 0, gp_XY( 0, 0 ) );
 677   aPoly->AddPoint( 0, gp_XY( 20, 0 ) );
 678   aPoly->AddPoint( 0, gp_XY( 10, 10 ) );
 679   aPoly->Update();
 680
 681   Handle(HYDROData_ImmersibleZone) aZone =
 682     Handle(HYDROData_ImmersibleZone)::DownCast( aDoc->CreateObject( KIND_IMMERSIBLE_ZONE ) );
 683   aZone->SetPolyline( aPoly );
 684   aZone->Update();
 685
 686   Handle(HYDROData_LandCoverMap) aMap =
 687     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 688   aMap->LocalPartition( Face( QList<double>() << 1 << 1 << 10 << 10 << 10 << 20 ), "Zones de champs cultivé à végétation haute" );
 689   aMap->LocalPartition( Face( QList<double>() << 5 << 5 << 10 << 5 << 10 << 8 << 5 << 12 << 5 << 8 ), "Zones de champs cultivé à végétation haute" );
 690   aMap->SetName( "test_LCM" );
 691
 692   TestViewer::show( aZone->GetTopShape(), AIS_Shaded, true, "LandCoverMap_TransparentPrs" );
 693
 694   aMap->SetTransparency( 0.5 );
 695   CPPUNIT_ASSERT_DOUBLES_EQUAL( 0.5, aMap->GetTransparency(), 1E-6 );
 696
 697   Handle(HYDROGUI_LandCoverMapPrs) aPrs = new HYDROGUI_LandCoverMapPrs( aMap );
 698   aPrs->SetTable( aTable );
 699   TestViewer::show( aPrs, AIS_Shaded, 0, true, "" );
 700
 701   CPPUNIT_ASSERT_IMAGES
 702   aDoc->Close();
 703 }
 704
 705 void test_HYDROData_LandCoverMap::test_add_triangles()
 706 {
 707   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 708
 709   Handle(HYDROData_PolylineXY) Polyline_1 =
 710     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 711
 712   Polyline_1->SetName( "Polyline_1" );
 713   Polyline_1->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 714   Polyline_1->AddPoint( 0, gp_XY( -264.84, 323.44 ) );
 715   Polyline_1->AddPoint( 0, gp_XY( 254.45, 301.19 ) );
 716   Polyline_1->AddPoint( 0, gp_XY( -291.54, -47.48 ) );
 717   Polyline_1->Update();
 718
 719   Handle(HYDROData_PolylineXY) Polyline_2 =
 720     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 721   Polyline_2->SetName( "Polyline_2" );
 722   Polyline_2->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, false );
 723   Polyline_2->AddPoint( 0, gp_XY( -42.28, -351.63 ) );
 724   Polyline_2->AddPoint( 0, gp_XY( 218.84, -146.88 ) );
 725   Polyline_2->AddPoint( 0, gp_XY( 413.20, -235.91 ) );
 726   Polyline_2->AddPoint( 0, gp_XY( 466.62, -388.72 ) );
 727   Polyline_2->Update();
 728
 729   Handle(HYDROData_PolylineXY) Polyline_3 =
 730     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 731   Polyline_3->SetName( "Polyline_3" );
 732   Polyline_3->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 733   Polyline_3->AddPoint( 0, gp_XY( -108.31, 213.65 ) );
 734   Polyline_3->AddPoint( 0, gp_XY( 127.60, 382.79 ) );
 735   Polyline_3->AddPoint( 0, gp_XY( 192.88, -91.99 ) );
 736   Polyline_3->Update();
 737
 738   Handle(HYDROData_PolylineXY) Polyline_4 =
 739     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 740   Polyline_4->SetName( "Polyline_4" );
 741   Polyline_4->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 742   Polyline_4->AddPoint( 0, gp_XY( 215.52, 344.61 ) );
 743   Polyline_4->AddPoint( 0, gp_XY( 68.25, 258.52 ) );
 744   Polyline_4->AddPoint( 0, gp_XY( 200.58, 154.65 ) );
 745   Polyline_4->Update();
 746
 747   Handle(HYDROData_PolylineXY) Polyline_5 =
 748     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 749   Polyline_5->SetName( "Polyline_5" );
 750   Polyline_5->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, false );
 751   Polyline_5->AddPoint( 0, gp_XY( 16.15, -58.08 ) );
 752   Polyline_5->AddPoint( 0, gp_XY( 116.47, 84.21 ) );
 753   Polyline_5->AddPoint( 0, gp_XY( 266.59, 29.43 ) );
 754   Polyline_5->Update();
 755
 756   Handle(HYDROData_LandCoverMap) aMap =
 757     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 758
 759   CPPUNIT_ASSERT_EQUAL( true, aMap->Add( Polyline_1, "" ) );
 760   CPPUNIT_ASSERT_EQUAL( true, aMap->Add( Polyline_3, "" ) );
 761   CPPUNIT_ASSERT_EQUAL( true, aMap->Add( Polyline_4, "" ) );
 762
 763   Handle(HYDROGUI_LandCoverMapPrs) aPrs = new HYDROGUI_LandCoverMapPrs( aMap );
 764   TopoDS_Shape aSh = aMap->GetShape();
 765   TestViewer::show( aPrs, AIS_Shaded, 0, true, "LandCover_Triangles" );
 766   CPPUNIT_ASSERT_IMAGES
 767
 768   CPPUNIT_ASSERT_EQUAL( true, aMap->Split( Polyline_5 ) );
 769   Handle(HYDROGUI_LandCoverMapPrs) aPrs2 = new HYDROGUI_LandCoverMapPrs( aMap );
 770   TestViewer::show( aPrs2, AIS_Shaded, 0, true, "LandCover_Triangles_Split" );
 771   CPPUNIT_ASSERT_IMAGES
 772
 773   aDoc->Close();
 774 }
 775
 776 void test_HYDROData_LandCoverMap::test_assign_to_calc_case()
 777 {
 778   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 779
 780   Handle(HYDROData_LandCoverMap) aMap =
 781     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 782   Handle(HYDROData_CalculationCase) aCalcCase =
 783     Handle(HYDROData_CalculationCase)::DownCast( aDoc->CreateObject( KIND_CALCULATION ) );
 784
 785   aCalcCase->SetLandCoverMap( aMap );
 786   CPPUNIT_ASSERT_EQUAL( aMap->myLab, aCalcCase->GetLandCoverMap()->myLab );
 787
 788   aDoc->Close();
 789 }
 790
 791 void test_HYDROData_LandCoverMap::test_shp_import_cyp()
 792 {
 793   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 794   QString aFileName = REF_DATA_PATH + "/cyprus_natural.shp";
 795   HYDROData_ShapeFile anImporter;
 796   QStringList PolygonList;
 797   TopTools_SequenceOfShape PolygonFaces;
 798   int Type = -1;
 799   int aStat = anImporter.ImportPolygons(aFileName, PolygonList, PolygonFaces, Type);
 800   CPPUNIT_ASSERT(aStat == 1);
 801   CPPUNIT_ASSERT_EQUAL(5, Type);
 802   CPPUNIT_ASSERT_EQUAL(268, PolygonFaces.Length());
 803
 804   Handle(HYDROData_LandCoverMap) LCM = Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 805   HYDROData_MapOfFaceToStricklerType aMapFace2ST;
 806
 807   for ( int i = 1; i <= PolygonFaces.Length(); i++ )
 808   {
 809     TopoDS_Shape aShape = PolygonFaces(i);
 810     if ( aShape.IsNull() )
 811       continue;
 812     aMapFace2ST.Add( TopoDS::Face( aShape ), "" );
 813   }
 814
 815   LCM->StoreLandCovers(aMapFace2ST);
 816   TestViewer::show( LCM->GetShape(), AIS_Shaded, true, "cyprus_natural_all" );
 817   CPPUNIT_ASSERT_IMAGES
 818   aDoc->Close();
 819 }
 820
 821 void test_HYDROData_LandCoverMap::test_shp_import_clc_dec()
 822 {
 823   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 824   QString aFileName = REF_DATA_PATH + "/CLC_decoupe.shp";
 825   HYDROData_ShapeFile anImporter;
 826   QStringList PolygonList;
 827   TopTools_SequenceOfShape PolygonFaces;
 828   int Type = -1;
 829   int aStat = anImporter.ImportPolygons(aFileName, PolygonList, PolygonFaces, Type);
 830   CPPUNIT_ASSERT(aStat == 1);
 831   CPPUNIT_ASSERT_EQUAL(5, Type);
 832   CPPUNIT_ASSERT_EQUAL(625, PolygonFaces.Length());
 833
 834   Handle(HYDROData_LandCoverMap) LCM = Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 835   HYDROData_MapOfFaceToStricklerType aMapFace2ST;
 836
 837   for ( int i = 150; i <= 350; i++ )
 838   {
 839     TopoDS_Shape aShape = PolygonFaces(i);
 840     if ( aShape.IsNull() )
 841       continue;
 842     aMapFace2ST.Add( TopoDS::Face( aShape ), "" );
 843   }
 844
 845   LCM->StoreLandCovers(aMapFace2ST);
 846   TestViewer::show( LCM->GetShape(), AIS_Shaded, true, "clc_dec_150_350" );
 847   CPPUNIT_ASSERT_IMAGES
 848   aDoc->Close();
 849 }
 850
 851 void test_HYDROData_LandCoverMap::test_shp_import_nld_areas()
 852 {
 853   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 854   QString aFileName = REF_DATA_PATH + "/NLD_water_areas_dcw.shp";
 855   HYDROData_ShapeFile anImporter;
 856   QStringList PolygonList;
 857   TopTools_SequenceOfShape PolygonFaces;
 858   int Type = -1;
 859   int aStat = anImporter.ImportPolygons(aFileName, PolygonList, PolygonFaces, Type);
 860   CPPUNIT_ASSERT(aStat == 1);
 861   CPPUNIT_ASSERT_EQUAL(5, Type);
 862   CPPUNIT_ASSERT_EQUAL(127, PolygonFaces.Length());
 863
 864   Handle(HYDROData_LandCoverMap) LCM = Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 865   HYDROData_MapOfFaceToStricklerType aMapFace2ST;
 866
 867   for ( int i = 1; i <= PolygonFaces.Length(); i++ )
 868   {
 869     TopoDS_Shape aShape = PolygonFaces(i);
 870     if ( aShape.IsNull() )
 871       continue;
 872     aMapFace2ST.Add( TopoDS::Face( aShape ), "" );
 873   }
 874
 875   LCM->StoreLandCovers(aMapFace2ST);
 876   TestViewer::show( LCM->GetShape(), AIS_Shaded, true, "NLD_water_areas_dcw" );
 877   CPPUNIT_ASSERT_IMAGES
 878   aDoc->Close();
 879 }
 880
 881 void test_HYDROData_LandCoverMap::test_shp_loop_back()
 882 {
 883   QString aFFileName = REF_DATA_PATH + "/CLC06-cut_1.shp";
 884   QString aSFileName = REF_DATA_PATH + "/CLC06-cut_1_res.shp";
 885   //This test verify only some geom data (shp+shx) without any dbf reading/writing
 886   {
 887     Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 888     HYDROData_ShapeFile anImporter;
 889     QStringList PolygonList;
 890     TopTools_SequenceOfShape PolygonFaces;
 891     int Type = -1;
 892
 893     //import LCM from file (#2-4; #12-14 polygons)
 894     CPPUNIT_ASSERT( anImporter.ImportPolygons(aFFileName, PolygonList, PolygonFaces, Type));
 895     Handle(HYDROData_LandCoverMap) LCM = Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 896     HYDROData_MapOfFaceToStricklerType aMapFace2ST;
 897
 898     aMapFace2ST.Add( TopoDS::Face( PolygonFaces(2) ), "" );
 899     aMapFace2ST.Add( TopoDS::Face( PolygonFaces(3) ), "" );
 900     aMapFace2ST.Add( TopoDS::Face( PolygonFaces(4) ), "" );
 901     aMapFace2ST.Add( TopoDS::Face( PolygonFaces(12) ), "" );
 902     aMapFace2ST.Add( TopoDS::Face( PolygonFaces(13) ), "" );
 903     aMapFace2ST.Add( TopoDS::Face( PolygonFaces(14) ), "" );
 904
 905     LCM->StoreLandCovers(aMapFace2ST);
 906     LCM->SetName("lcm_1");
 907
 908     TopoDS_Shape aSh = LCM->GetShape();
 909     BRepCheck_Analyzer aBCA(aSh);
 910     CPPUNIT_ASSERT(aBCA.IsValid());
 911
 912     //export lcm
 913     HYDROData_ShapeFile anExporter;
 914     QStringList aNonExpList;
 915     anExporter.Export(aSFileName, LCM, aNonExpList);
 916     CPPUNIT_ASSERT (aNonExpList.empty());
 917
 918     aDoc->Close();
 919   }
 920   {
 921     Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
 922     HYDROData_ShapeFile anImporter;
 923     QStringList PolygonList;
 924     TopTools_SequenceOfShape PolygonFaces;
 925     int Type = -1;
 926
 927     //import all
 928     CPPUNIT_ASSERT( anImporter.ImportPolygons(aSFileName, PolygonList, PolygonFaces, Type));
 929     Handle(HYDROData_LandCoverMap) LCM = Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 930     HYDROData_MapOfFaceToStricklerType aMapFace2ST;
 931     CPPUNIT_ASSERT_EQUAL(PolygonFaces.Length(), 6);
 932     for ( int i = 1; i <= PolygonFaces.Length(); i++ )
 933       aMapFace2ST.Add( TopoDS::Face( PolygonFaces(i) ), "" );
 934
 935     LCM->StoreLandCovers(aMapFace2ST);
 936     LCM->SetName("lcm_2");
 937
 938     TopoDS_Shape OutSh = LCM->GetShape();
 939     BRepCheck_Analyzer aBCA(OutSh);
 940     CPPUNIT_ASSERT(aBCA.IsValid());
 941
 942     TestViewer::show( OutSh, AIS_Shaded, true, "CLC06-cut_1_res" );
 943     CPPUNIT_ASSERT_IMAGES
 944
 945     aDoc->Close();
 946   }
 947   CPPUNIT_ASSERT_EQUAL(0, remove(aSFileName.toStdString().c_str()));
 948 }
 949
 950 void test_HYDROData_LandCoverMap::test_add_2_section_poly()
 951 {
 952   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( 1 );
 953
 954   Handle(HYDROData_PolylineXY) Polyline_12 =
 955     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 956   Polyline_12->SetName( "Polyline_12" );
 957
 958   Polyline_12->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 959   Polyline_12->AddPoint( 0, gp_XY( -148.70, 145.68 ) );
 960   Polyline_12->AddPoint( 0, gp_XY( -81.71, 210.97 ) );
 961   Polyline_12->AddPoint( 0, gp_XY( -21.09, 122.79 ) );
 962   Polyline_12->AddPoint( 0, gp_XY( -106.30, 102.44 ) );
 963   Polyline_12->AddSection( "Section_2", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 964   Polyline_12->AddPoint( 1, gp_XY( -7.09, 227.08 ) );
 965   Polyline_12->AddPoint( 1, gp_XY( 46.32, 228.78 ) );
 966   Polyline_12->AddPoint( 1, gp_XY( 69.64, 165.61 ) );
 967   Polyline_12->AddPoint( 1, gp_XY( 13.68, 156.28 ) );
 968
 969   Polyline_12->Update();
 970
 971   Handle(HYDROData_LandCoverMap) LCM = Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
 972   CPPUNIT_ASSERT_EQUAL( true, LCM->Add( Polyline_12, "" ) );
 973
 974   TestViewer::show( LCM->GetShape(), AIS_Shaded, true, "lcm_poly_2_sections" );
 975   CPPUNIT_ASSERT_IMAGES
 976
 977   aDoc->Close();
 978 }
 979
 980 void test_HYDROData_LandCoverMap::test_split_2()
 981 {
 982   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document( 1 );
 983
 984   Handle(HYDROData_PolylineXY) Polyline_1 =
 985     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 986   Polyline_1->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 987   Polyline_1->AddPoint( 0, gp_XY( -264.84, 323.44 ) );
 988   Polyline_1->AddPoint( 0, gp_XY( 254.45, 301.19 ) );
 989   Polyline_1->AddPoint( 0, gp_XY( -291.54, -47.48 ) );
 990   Polyline_1->Update();
 991
 992   Handle(HYDROData_PolylineXY) Polyline_6 =
 993     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
 994   Polyline_6->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
 995   Polyline_6->AddPoint( 0, gp_XY( -155.15, 199.47 ) );
 996   Polyline_6->AddPoint( 0, gp_XY( -71.19, 289.12 ) );
 997   Polyline_6->AddPoint( 0, gp_XY( 22.01, 211.57 ) );
 998   Polyline_6->AddPoint( 0, gp_XY( -84.00, 134.73 ) );
 999   Polyline_6->Update();
1000
1001   Handle(HYDROData_PolylineXY) Polyline_7 =
1002     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
1003   Polyline_7->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
1004   Polyline_7->AddPoint( 0, gp_XY( -50.56, 232.20 ) );
1005   Polyline_7->AddPoint( 0, gp_XY( 40.51, 286.98 ) );
1006   Polyline_7->AddPoint( 0, gp_XY( 118.77, 256.39 ) );
1007   Polyline_7->AddPoint( 0, gp_XY( 72.52, 114.10 ) );
1008   Polyline_7->Update();
1009
1010   Handle(HYDROData_PolylineXY) Polyline_8 =
1011     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
1012   Polyline_8->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, true );
1013   Polyline_8->AddPoint( 0, gp_XY( -44.16, 190.93 ) );
1014   Polyline_8->AddPoint( 0, gp_XY( 9.91, 226.51 ) );
1015   Polyline_8->AddPoint( 0, gp_XY( 35.53, 175.99 ) );
1016   Polyline_8->AddPoint( 0, gp_XY( -9.30, 157.50 ) );
1017   Polyline_8->Update();
1018
1019   Handle(HYDROData_PolylineXY) Polyline_9 =
1020     Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) );
1021   Polyline_9->AddSection( "Section_1", HYDROData_PolylineXY::SECTION_POLYLINE, false );
1022   Polyline_9->AddPoint( 0, gp_XY( -103.92, 288.40 ) );
1023   Polyline_9->AddPoint( 0, gp_XY( -100.09, 285.22 ) );
1024   Polyline_9->AddPoint( 0, gp_XY( -96.10, 281.90 ) );
1025   Polyline_9->AddPoint( 0, gp_XY( -90.43, 277.19 ) );
1026   Polyline_9->AddPoint( 0, gp_XY( -53.41, 246.43 ) );
1027   Polyline_9->AddPoint( 0, gp_XY( 10.62, 245.01 ) );
1028
1029
1030   Polyline_9->Update();
1031
1032
1033   Handle(HYDROData_LandCoverMap) LCM = Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
1034   CPPUNIT_ASSERT_EQUAL( true, LCM->Add( Polyline_1, "" ) );
1035   CPPUNIT_ASSERT_EQUAL( true, LCM->Add( Polyline_6, "" ) );
1036   CPPUNIT_ASSERT_EQUAL( true, LCM->Add( Polyline_7, "" ) );
1037
1038   CPPUNIT_ASSERT_EQUAL( true, LCM->Split( Polyline_9 ) );
1039
1040   //Handle(HYDROGUI_LandCoverMapPrs) aPrs1 = new HYDROGUI_LandCoverMapPrs( LCM );
1041   TestViewer::show( LCM->GetShape(), AIS_Shaded, /*0,*/ true, "LCM_split_2a" );
1042   CPPUNIT_ASSERT_IMAGES
1043
1044   CPPUNIT_ASSERT_EQUAL( true, LCM->Split( Polyline_8 ) );
1045
1046   // Handle(HYDROGUI_LandCoverMapPrs) aPrs2 = new HYDROGUI_LandCoverMapPrs( LCM );
1047   TestViewer::show( LCM->GetShape(), AIS_Shaded, /*0,*/ true, "LCM_split_2b" );
1048   CPPUNIT_ASSERT_IMAGES
1049
1050   aDoc->Close();
1051 }
1052
1053 void test_HYDROData_LandCoverMap::test_export_telemac()
1054 {
1055   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
1056
1057   Handle(HYDROData_StricklerTable) aTable =
1058     Handle(HYDROData_StricklerTable)::DownCast( aDoc->CreateObject( KIND_STRICKLER_TABLE ) );
1059   CPPUNIT_ASSERT_EQUAL( true, aTable->Import( DEF_STR_PATH ) );
1060
1061   Handle(HYDROData_LandCoverMap) aMap =
1062     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
1063
1064   TopoDS_Face aLC1 = Face( QList<double>() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 );
1065   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "Zones de champs cultivé à végétation basse" ) );
1066
1067   TopoDS_Face aLC2 = Face( QList<double>() << 21 << 34 << 24 << 25 << 37   << 37 << 40  << 61 <<
1068                                               44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 <<
1069                                              128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 <<
1070                                               31 << 114 );
1071   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "Zones de champs cultivé à végétation haute" ) );
1072
1073   TopoDS_Face aLC3 = Face( QList<double>() << 4 << 54 << 1   << 47 << 51  << 45 <<
1074                                             127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 );
1075   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "Zones de champs, prairies, sans cultures" ) );
1076
1077   QString aTmpFileName = "test.telemac";
1078   QString aTmpPath = QDir::tempPath() + "/" + aTmpFileName;
1079   CPPUNIT_ASSERT_EQUAL( true, aMap->ExportTelemac( aTmpPath, 1E-4, aTable ) );
1080   CPPUNIT_ASSERT_SCRIPTS_EQUAL( aTmpFileName, true, true, 0 );
1081
1082   aDoc->Close();
1083 }
1084
1085 void test_HYDROData_LandCoverMap::test_copy()
1086 {
1087   Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1);
1088
1089   Handle(HYDROData_LandCoverMap) aMap =
1090     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
1091   aMap->SetName( "map_1" );
1092
1093   TopoDS_Face aLC1 = Face( QList<double>() << 10 << 10 << 30 << 10 << 20 << 20 );
1094   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "Zones de champs cultivé à végétation basse" ) );
1095
1096   TopoDS_Face aLC2 = Face( QList<double>() << 110 << 10 << 130 << 10 << 120 << 20 );
1097   CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "Zones de champs cultivé à végétation haute" ) );
1098
1099   Handle(HYDROData_LandCoverMap) aMap2 =
1100     Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) );
1101   aMap->CopyTo( aMap2, true );
1102
1103   CPPUNIT_ASSERT_EQUAL( QString( "map_2" ), aMap2->GetName() );
1104   HYDROData_LandCoverMap::Explorer anIt( aMap2 );
1105   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
1106   CPPUNIT_ASSERT_EQUAL( QString( "Zones de champs cultivé à végétation basse" ), anIt.StricklerType() );
1107   CPPUNIT_ASSERT( anIt.Face().TShape()!=aLC1.TShape() );
1108   anIt.Next();
1109   CPPUNIT_ASSERT_EQUAL( true, anIt.More() );
1110   CPPUNIT_ASSERT_EQUAL( QString( "Zones de champs cultivé à végétation haute" ), anIt.StricklerType() );
1111   CPPUNIT_ASSERT( anIt.Face().TShape()!=aLC2.TShape() );
1112   anIt.Next();
1113   CPPUNIT_ASSERT_EQUAL( false, anIt.More() );
1114
1115   TopoDS_Shape aShape1 = aMap->GetShape();
1116   TopoDS_Shape aShape2 = aMap2->GetShape();
1117   gp_Trsf aTr;
1118   aTr.SetTranslation( gp_Vec( 0, 50, 0 ) );
1119   aShape2.Move( TopLoc_Location( aTr ) );
1120   TestViewer::show( aShape1, AIS_Shaded, /*0,*/ true, "LCM_copied" );
1121   TestViewer::show( aShape2, AIS_Shaded, /*0,*/ true, "" );
1122   CPPUNIT_ASSERT_IMAGES
1123
1124   aDoc->Close();
1125 }