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