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