X-Git-Url: http://git.salome-platform.org/gitweb/?a=blobdiff_plain;f=src%2FHYDRO_tests%2Ftest_HYDROData_LandCoverMap.cxx;h=1ae4bd40d0eb23d0ee5dcc620fec4370bbda463a;hb=9c947f35615e69e9e54a8c4b074dd1f2be13689c;hp=e1589b7b13e913c608eb996b7b91eacb58460908;hpb=dc6f07c898307b3b6e539d805f4e83d77e4d5215;p=modules%2Fhydro.git diff --git a/src/HYDRO_tests/test_HYDROData_LandCoverMap.cxx b/src/HYDRO_tests/test_HYDROData_LandCoverMap.cxx index e1589b7b..1ae4bd40 100644 --- a/src/HYDRO_tests/test_HYDROData_LandCoverMap.cxx +++ b/src/HYDRO_tests/test_HYDROData_LandCoverMap.cxx @@ -1,4 +1,4 @@ -// Copyright (C) 2014-2015 EDF-R&D +// Copyright (C) 2014-2015 EDF-R&D // This library is free software; you can redistribute it and/or // modify it under the terms of the GNU Lesser General Public // License as published by the Free Software Foundation; either @@ -19,60 +19,28 @@ #include #include #include +#include #include +#include #include +#include #include #include #include -#include -#include -#include -#include -#include #include +#include #include #include +#include #include #include - +#include +#include #include #include const QString REF_DATA_PATH = qgetenv( "HYDRO_REFERENCE_DATA" ); - -TopoDS_Edge Spline( const QList& theXYList, bool isClosed = false ) -{ - int n = theXYList.size()/2; - Handle(TColgp_HArray1OfPnt) aPointsArray = new TColgp_HArray1OfPnt( 1, n ); - for( int i=1; i<=n; i++ ) - { - double x = theXYList[2*i-2]; - double y = theXYList[2*i-1]; - gp_Pnt aPnt( x, y, 0 ); - aPointsArray->SetValue( i, aPnt ); - } - GeomAPI_Interpolate anInterpolator( aPointsArray, isClosed, 1E-3 ); - anInterpolator.Perform(); - bool aResult = anInterpolator.IsDone() == Standard_True; - if( aResult ) - { - Handle( Geom_BSplineCurve ) aCurve = anInterpolator.Curve(); - return BRepBuilderAPI_MakeEdge( aCurve ).Edge(); - } - else - return TopoDS_Edge(); -} - -TopoDS_Face LandCover( const QList& theXYList ) -{ - TopoDS_Edge anEdge = Spline( theXYList, true ); - if( anEdge.IsNull() ) - return TopoDS_Face(); - - TopoDS_Wire aWire = BRepBuilderAPI_MakeWire( anEdge ).Wire(); - TopoDS_Face aFace = BRepBuilderAPI_MakeFace( aWire, Standard_True ).Face(); - return aFace; -} +const QString DEF_STR_PATH = qgetenv( "HYDRO_SRC_DIR" ) + "/src/HYDROGUI/resources/def_strickler_table.txt"; void test_HYDROData_LandCoverMap::test_add_2_objects() { @@ -83,10 +51,10 @@ void test_HYDROData_LandCoverMap::test_add_2_objects() CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() ); - TopoDS_Face aLC1 = LandCover( QList() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 ); + TopoDS_Face aLC1 = Face( QList() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) ); - TopoDS_Face aLC2 = LandCover( QList() << 30 << 20 << 60 << 10 << 70 << 35 << 40 << 40 ); + TopoDS_Face aLC2 = Face( QList() << 30 << 20 << 60 << 10 << 70 << 35 << 40 << 40 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "test2" ) ); TestViewer::show( aMap->GetShape(), AIS_Shaded, true, "LandCoverMap_Add_2_Objects" ); @@ -113,12 +81,12 @@ void test_HYDROData_LandCoverMap::test_split() CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() ); - TopoDS_Face aLC = LandCover( QList() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 ); + TopoDS_Face aLC = Face( QList() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC, "test1" ) ); Handle(HYDROData_PolylineXY) aPolyline = Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) ); - TopoDS_Wire aWire = BRepBuilderAPI_MakeWire( Spline( QList() << 10 << 40 << 30 << 10 << 40 << 10 << 60 << 10 ) ).Wire(); + TopoDS_Wire aWire = Wire( QList() << 10 << 40 << 30 << 10 << 40 << 10 << 60 << 10, false ); aPolyline->SetShape( aWire ); CPPUNIT_ASSERT_EQUAL( true, aMap->Split( aPolyline ) ); @@ -148,18 +116,18 @@ void test_HYDROData_LandCoverMap::test_incomplete_split() CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() ); - TopoDS_Face aLC = LandCover( QList() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 ); + TopoDS_Face aLC = Face( QList() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC, "test1" ) ); Handle(HYDROData_PolylineXY) aPolyline = Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) ); - TopoDS_Wire aWire = BRepBuilderAPI_MakeWire( Spline( QList() << 10 << 40 << 30 << 10 << 40 << 10 ) ).Wire(); + TopoDS_Wire aWire = Wire( QList() << 10 << 40 << 30 << 10 << 40 << 10, false ); aPolyline->SetShape( aWire ); CPPUNIT_ASSERT_EQUAL( true, aMap->Split( aPolyline ) ); TestViewer::show( aMap->GetShape(), AIS_Shaded, true, "LandCoverMap_Split_2" ); - TestViewer::show( aWire, 0, true, Qt::blue ); + TestViewer::show( aWire, 0, true, Qt::green ); CPPUNIT_ASSERT_IMAGES HYDROData_LandCoverMap::Iterator anIt( aMap ); @@ -180,18 +148,18 @@ void test_HYDROData_LandCoverMap::test_merge() CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() ); - TopoDS_Face aLC1 = LandCover( QList() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 ); + TopoDS_Face aLC1 = Face( QList() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) ); - TopoDS_Face aLC2 = LandCover( QList() << 21 << 34 << 24 << 25 << 37 << 37 << 40 << 61 << - 44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 << - 128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 << - 31 << 114 ); + TopoDS_Face aLC2 = Face( QList() << 21 << 34 << 24 << 25 << 37 << 37 << 40 << 61 << + 44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 << + 128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 << + 31 << 114 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "test2" ) ); - TopoDS_Face aLC3 = LandCover( QList() << 4 << 54 << 1 << 47 << 51 << 45 << - 127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 ); + TopoDS_Face aLC3 = Face( QList() << 4 << 54 << 1 << 47 << 51 << 45 << + 127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "test3" ) ); //TestViewer::show( aMap->GetShape(), AIS_Shaded, true ); @@ -225,18 +193,18 @@ void test_HYDROData_LandCoverMap::test_remove() CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() ); - TopoDS_Face aLC1 = LandCover( QList() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 ); + TopoDS_Face aLC1 = Face( QList() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) ); - TopoDS_Face aLC2 = LandCover( QList() << 21 << 34 << 24 << 25 << 37 << 37 << 40 << 61 << - 44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 << - 128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 << - 31 << 114 ); + TopoDS_Face aLC2 = Face( QList() << 21 << 34 << 24 << 25 << 37 << 37 << 40 << 61 << + 44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 << + 128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 << + 31 << 114 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "test2" ) ); - TopoDS_Face aLC3 = LandCover( QList() << 4 << 54 << 1 << 47 << 51 << 45 << - 127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 ); + TopoDS_Face aLC3 = Face( QList() << 4 << 54 << 1 << 47 << 51 << 45 << + 127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 ); CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "test3" ) ); QString aType1, aType2; @@ -253,7 +221,7 @@ void test_HYDROData_LandCoverMap::test_remove() aDoc->Close(); } -void test_HYDROData_LandCoverMap::test_merge_faces() +void test_HYDROData_LandCoverMap::test_merge_faces_boxes() { TopoDS_Shape pp1, pp2, pp3, pp4; BRep_Builder BB; @@ -267,13 +235,13 @@ void test_HYDROData_LandCoverMap::test_merge_faces() CPPUNIT_ASSERT(!pp3.IsNull()); CPPUNIT_ASSERT(!pp4.IsNull()); - //Test mergeFaces() func; USD == true + //Test mergeFaces() func // boxes // USD == true { TopTools_ListOfShape Faces; Faces.Append(pp1); Faces.Append(pp2); TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true ); - TestViewer::show( aMergedFace, AIS_Shaded, true, "cs11" ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_11" ); CPPUNIT_ASSERT_IMAGES CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_FACE); } @@ -282,7 +250,7 @@ void test_HYDROData_LandCoverMap::test_merge_faces() Faces.Append(pp1); Faces.Append(pp3); TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true ); - TestViewer::show( aMergedFace, AIS_Shaded, true, "cs12" ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_12" ); CPPUNIT_ASSERT_IMAGES CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); } @@ -292,7 +260,7 @@ void test_HYDROData_LandCoverMap::test_merge_faces() Faces.Append(pp2); Faces.Append(pp4); TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true ); - TestViewer::show( aMergedFace, AIS_Shaded, true, "cs13" ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_13" ); CPPUNIT_ASSERT_IMAGES CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); } @@ -303,19 +271,19 @@ void test_HYDROData_LandCoverMap::test_merge_faces() Faces.Append(pp4); Faces.Append(pp2); TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true ); - TestViewer::show( aMergedFace, AIS_Shaded, true, "cs14" ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_14" ); CPPUNIT_ASSERT_IMAGES CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_FACE); } // - //Test mergeFaces() func; USD == false + //Test mergeFaces() func // boxes // USD == false { TopTools_ListOfShape Faces; Faces.Append(pp1); Faces.Append(pp2); TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false ); - TestViewer::show( aMergedFace, AIS_Shaded, true, "cs21" ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_21" ); CPPUNIT_ASSERT_IMAGES CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); } @@ -324,7 +292,7 @@ void test_HYDROData_LandCoverMap::test_merge_faces() Faces.Append(pp1); Faces.Append(pp3); TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false ); - TestViewer::show( aMergedFace, AIS_Shaded, true, "cs22" ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_22" ); CPPUNIT_ASSERT_IMAGES CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); } @@ -334,8 +302,8 @@ void test_HYDROData_LandCoverMap::test_merge_faces() Faces.Append(pp2); Faces.Append(pp4); TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false ); - TestViewer::show( aMergedFace, AIS_Shaded, true, "cs23" ); - CPPUNIT_ASSERT_IMAGES + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_23" ); + CPPUNIT_ASSERT_IMAGES CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); } { @@ -345,8 +313,278 @@ void test_HYDROData_LandCoverMap::test_merge_faces() Faces.Append(pp4); Faces.Append(pp2); TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false ); - TestViewer::show( aMergedFace, AIS_Shaded, true, "cs24" ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_boxes_24" ); + CPPUNIT_ASSERT_IMAGES + CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); + } +} + + +void test_HYDROData_LandCoverMap::test_merge_faces_circles() +{ + TopoDS_Shape ff1, ff2, ff3; + BRep_Builder BB; + BRepTools::Read(ff1, (REF_DATA_PATH + "/ff1.brep").toStdString().c_str(), BB); + BRepTools::Read(ff2, (REF_DATA_PATH + "/ff2.brep").toStdString().c_str(), BB); + BRepTools::Read(ff3, (REF_DATA_PATH + "/ff3.brep").toStdString().c_str(), BB); + + CPPUNIT_ASSERT(!ff1.IsNull()); + CPPUNIT_ASSERT(!ff2.IsNull()); + CPPUNIT_ASSERT(!ff3.IsNull()); + + //Test mergeFaces() func // circles // USD == true + { + TopTools_ListOfShape Faces; + Faces.Append(ff1); + Faces.Append(ff2); + TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_11" ); + CPPUNIT_ASSERT_IMAGES + CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_FACE); + } + { + TopTools_ListOfShape Faces; + Faces.Append(ff1); + Faces.Append(ff3); + TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_12" ); + CPPUNIT_ASSERT_IMAGES + CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); + } + { + TopTools_ListOfShape Faces; + Faces.Append(ff1); + Faces.Append(ff2); + Faces.Append(ff3); + TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_13" ); + CPPUNIT_ASSERT_IMAGES + CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); + } + { + TopTools_ListOfShape Faces; + Faces.Append(ff1); + Faces.Append(ff3); + Faces.Append(ff2); + TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, true ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_14" ); + CPPUNIT_ASSERT_IMAGES + CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); + } + // + + //Test mergeFaces() func // circles // USD == false + { + TopTools_ListOfShape Faces; + Faces.Append(ff1); + Faces.Append(ff2); + TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_21" ); CPPUNIT_ASSERT_IMAGES CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); } + { + TopTools_ListOfShape Faces; + Faces.Append(ff1); + Faces.Append(ff3); + TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_22" ); + CPPUNIT_ASSERT_IMAGES + CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); + } + { + TopTools_ListOfShape Faces; + Faces.Append(ff1); + Faces.Append(ff2); + Faces.Append(ff3); + TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_23" ); + CPPUNIT_ASSERT_IMAGES + CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); + } + { + TopTools_ListOfShape Faces; + Faces.Append(ff1); + Faces.Append(ff3); + Faces.Append(ff2); + TopoDS_Shape aMergedFace = HYDROData_LandCoverMap::MergeFaces( Faces, false ); + TestViewer::show( aMergedFace, AIS_Shaded, true, "merge_faces_circles_24" ); + CPPUNIT_ASSERT_IMAGES + CPPUNIT_ASSERT(aMergedFace.ShapeType() == TopAbs_SHELL); + } + +} + +void test_HYDROData_LandCoverMap::test_import_dbf() +{ + Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1); + + Handle(HYDROData_LandCoverMap) aMap = + Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) ); + + CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() ); + QString aFileName = REF_DATA_PATH + "/t1.dbf"; + QStringList aDBFV; + QStringList aST; + aDBFV.append("100"); + aDBFV.append("200"); + aDBFV.append("300"); + aST.append("water"); + aST.append("forest"); + aST.append("road"); + // aST.size() == aDBFV.size()!! + Handle(HYDROData_PolylineXY) aPolyline = + Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) ); + TopoDS_Wire aWire = Wire( QList() << 10 << 40 << 30 << 10 << 40 << 10 << 60 << 10 ); + aPolyline->SetShape( aWire ); + for (int i = 0; i < 3; i++) + aMap->Add(aPolyline, ""); + TopoDS_Face aLC1 = Face( QList() << 10 << 10 << 50 << 20 << 30 << 50 << 15 << 30 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) ); + + CPPUNIT_ASSERT_EQUAL( true, aMap->Split( aPolyline ) ); + QList Inds = QList() << 1 << 2 << 3; + aMap->ImportDBF(aFileName, "TESTFIELD1", aDBFV, aST, Inds); + HYDROData_LandCoverMap::Iterator anIt( aMap ); + + CPPUNIT_ASSERT_EQUAL( QString( "forest" ), anIt.StricklerType() ); + anIt.Next(); + CPPUNIT_ASSERT_EQUAL( QString( "road" ), anIt.StricklerType() ); +} + +void test_HYDROData_LandCoverMap::test_land_cover_prs_by_types() +{ + Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1); + + Handle(HYDROData_StricklerTable) aTable = + Handle(HYDROData_StricklerTable)::DownCast( aDoc->CreateObject( KIND_STRICKLER_TABLE ) ); + CPPUNIT_ASSERT_EQUAL( true, aTable->Import( DEF_STR_PATH ) ); + + Handle(HYDROData_LandCoverMap) aMap = + Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) ); + + TopoDS_Face aLC1 = Face( QList() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "Zones de champs cultivé à végétation basse" ) ); + + TopoDS_Face aLC2 = Face( QList() << 21 << 34 << 24 << 25 << 37 << 37 << 40 << 61 << + 44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 << + 128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 << + 31 << 114 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "Zones de champs cultivé à végétation haute" ) ); + + TopoDS_Face aLC3 = Face( QList() << 4 << 54 << 1 << 47 << 51 << 45 << + 127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "Zones de champs, prairies, sans cultures" ) ); + + // build presentation object + Handle(HYDROGUI_LandCoverMapPrs) aPrs = new HYDROGUI_LandCoverMapPrs( aMap ); + aPrs->SetTable( aTable ); + // show presentation in viewer + TestViewer::show( aPrs, AIS_Shaded, 4, true, "LandCoverMap_PrsByTypes" ); + // select one of faces (first) + TestViewer::select( 200, 300 ); + CPPUNIT_ASSERT_IMAGES + + aDoc->Close(); +} + +void test_HYDROData_LandCoverMap::test_land_cover_prs_by_coeff() +{ + Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1); + + Handle(HYDROData_StricklerTable) aTable = + Handle(HYDROData_StricklerTable)::DownCast( aDoc->CreateObject( KIND_STRICKLER_TABLE ) ); + CPPUNIT_ASSERT_EQUAL( true, aTable->Import( DEF_STR_PATH ) ); + + Handle(HYDROData_LandCoverMap) aMap = + Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) ); + + TopoDS_Face aLC1 = Face( QList() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "Zones de champs cultivé à végétation basse" ) ); + + TopoDS_Face aLC2 = Face( QList() << 21 << 34 << 24 << 25 << 37 << 37 << 40 << 61 << + 44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 << + 128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 << + 31 << 114 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "Zones de champs cultivé à végétation haute" ) ); + + TopoDS_Face aLC3 = Face( QList() << 4 << 54 << 1 << 47 << 51 << 45 << + 127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "Zones de champs, prairies, sans cultures" ) ); + + Handle(HYDROGUI_LandCoverMapPrs) aPrs = new HYDROGUI_LandCoverMapPrs( aMap ); + aPrs->SetTable( aTable ); + aPrs->SetColorScale( TestViewer::showColorScale( true ) ); + TestViewer::show( aPrs, AIS_Shaded, 4, true, "LandCoverMap_PrsByCoeff" ); + CPPUNIT_ASSERT_IMAGES + + TestViewer::showColorScale( false ); + aDoc->Close(); +} + +void test_HYDROData_LandCoverMap::test_dump_python() +{ + Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1); + + Handle(HYDROData_LandCoverMap) aMap = + Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) ); + aMap->SetName( "test_LCM" ); + + CPPUNIT_ASSERT_EQUAL( KIND_LAND_COVER_MAP, aMap->GetKind() ); + + TopoDS_Face aLC1 = Face( QList() << 12 << 19 << 82 << 9 << 126 << 53 << 107 << 80 << 29 << 75 ); + + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC1, "test1" ) ); + + TopoDS_Face aLC2 = Face( QList() << 21 << 34 << 24 << 25 << 37 << 37 << 40 << 61 << + 44 << 95 << 85 << 100 << 104 << 66 << 107 << 33 << + 128 << 18 << 140 << 50 << 131 << 89 << 104 << 111 << + 31 << 114 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC2, "test2" ) ); + + TopoDS_Face aLC3 = Face( QList() << 4 << 54 << 1 << 47 << 51 << 45 << + 127 << 42 << 145 << 43 << 148 << 60 << 90 << 65 ); + CPPUNIT_ASSERT_EQUAL( true, aMap->LocalPartition( aLC3, "test3" ) ); + + QString aTmpPath = QDir::tempPath() + "/lc_dump.py"; + CPPUNIT_ASSERT_EQUAL( true, aDoc->DumpToPython( aTmpPath, false ) ); + + //TODO: compare files +} + +void test_HYDROData_LandCoverMap::test_transparent_prs() +{ + Handle(HYDROData_Document) aDoc = HYDROData_Document::Document(1); + + Handle(HYDROData_StricklerTable) aTable = + Handle(HYDROData_StricklerTable)::DownCast( aDoc->CreateObject( KIND_STRICKLER_TABLE ) ); + CPPUNIT_ASSERT_EQUAL( true, aTable->Import( DEF_STR_PATH ) ); + + Handle(HYDROData_PolylineXY) aPoly = + Handle(HYDROData_PolylineXY)::DownCast( aDoc->CreateObject( KIND_POLYLINEXY ) ); + aPoly->AddSection( "", HYDROData_PolylineXY::SECTION_SPLINE, true ); + aPoly->AddPoint( 0, gp_XY( 0, 0 ) ); + aPoly->AddPoint( 0, gp_XY( 20, 0 ) ); + aPoly->AddPoint( 0, gp_XY( 10, 10 ) ); + aPoly->Update(); + + Handle(HYDROData_ImmersibleZone) aZone = + Handle(HYDROData_ImmersibleZone)::DownCast( aDoc->CreateObject( KIND_IMMERSIBLE_ZONE ) ); + aZone->SetPolyline( aPoly ); + aZone->Update(); + + Handle(HYDROData_LandCoverMap) aMap = + Handle(HYDROData_LandCoverMap)::DownCast( aDoc->CreateObject( KIND_LAND_COVER_MAP ) ); + aMap->LocalPartition( Face( QList() << 1 << 1 << 10 << 10 << 10 << 20 ), "Zones de champs cultivé à végétation haute" ); + aMap->LocalPartition( Face( QList() << 5 << 5 << 10 << 5 << 10 << 8 << 5 << 12 << 5 << 8 ), "Zones de champs cultivé à végétation haute" ); + aMap->SetName( "test_LCM" ); + + TestViewer::show( aZone->GetTopShape(), AIS_Shaded, true, "LandCoverMap_TransparentPrs" ); + + Handle(HYDROGUI_LandCoverMapPrs) aPrs = new HYDROGUI_LandCoverMapPrs( aMap ); + aPrs->SetTable( aTable ); + aPrs->SetTransparency( 0.5 ); + TestViewer::show( aPrs, AIS_Shaded, 0, true, "" ); + + CPPUNIT_ASSERT_IMAGES }