[Test] // testing the Initialise method
public void Initialise()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Points");
ElementSet twoPointsElementSet = CreateElementSet("2 Points");
ElementMapper elementMapper = new ElementMapper();
IIdentifiable[] methods = SpatialAdaptedOutputFactory.GetAvailableMethods(fourPointsElementSet.ElementType, twoPointsElementSet.ElementType);
elementMapper.Initialise(methods[0], fourPointsElementSet, twoPointsElementSet);
double calculated = elementMapper.GetValueFromMappingMatrix(0, 0);
const double expected = 1;
Assert.AreEqual(expected, calculated);
}
[Test] // testing the Initialise method
public void Initialise()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Points");
ElementSet twoPointsElementSet = CreateElementSet("2 Points");
ElementMapper elementMapper = new ElementMapper();
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(fourPointsElementSet.ElementType, twoPointsElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, twoPointsElementSet);
double Calculated = elementMapper.GetValueFromMappingMatrix(0, 0);
double Expected = 1;
Assert.AreEqual(Expected, Calculated);
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_PointPolygon()
{
ElementSet gridElementSet = new ElementSet("gridElm", "G1", ElementType.XYPolygon, new SpatialReference("ref"));
ElementSet fourPointsElementSet = new ElementSet("4 points", "4P", ElementType.XYPoint, new SpatialReference("DummyID"));
Vertex v_0_20 = new Vertex(0, 20, 0);
Vertex v_0_10 = new Vertex(0, 10, 0);
Vertex v_0_0 = new Vertex(0, 0, 0);
Vertex v_0_15 = new Vertex(0, 15, 0);
Vertex v_5_15 = new Vertex(5, 15, 0);
Vertex v_10_20 = new Vertex(10, 20, 0);
Vertex v_10_15 = new Vertex(10, 15, 0);
Vertex v_10_10 = new Vertex(10, 10, 0);
Vertex v_10_0 = new Vertex(10, 0, 0);
Vertex v_15_15 = new Vertex(15, 15, 0);
Vertex v_15_5 = new Vertex(15, 5, 0);
Vertex v_20_20 = new Vertex(20, 20, 0);
Vertex v_20_10 = new Vertex(20, 10, 0);
Element square1 = new Element("square1");
Element square2 = new Element("square2");
Element square3 = new Element("square3");
square1.AddVertex(v_0_20);
square1.AddVertex(v_0_10);
square1.AddVertex(v_10_10);
square1.AddVertex(v_10_20);
square2.AddVertex(v_10_20);
square2.AddVertex(v_10_10);
square2.AddVertex(v_20_10);
square2.AddVertex(v_20_20);
square3.AddVertex(v_0_10);
square3.AddVertex(v_0_0);
square3.AddVertex(v_10_0);
square3.AddVertex(v_10_10);
gridElementSet.AddElement(square1);
gridElementSet.AddElement(square2);
gridElementSet.AddElement(square3);
Element point_5_15 = new Element("point 5, 15");
Element point_10_15 = new Element("point 10, 15");
Element point_15_15 = new Element("point 15, 15");
Element point_15_5 = new Element("point 15, 5");
point_5_15.AddVertex(v_5_15);
point_10_15.AddVertex(v_10_15);
point_15_15.AddVertex(v_15_15);
point_15_5.AddVertex(v_15_5);
fourPointsElementSet.AddElement(point_5_15);
fourPointsElementSet.AddElement(point_10_15);
fourPointsElementSet.AddElement(point_15_15);
fourPointsElementSet.AddElement(point_15_5);
ElementMapper elementMapper = new ElementMapper();
// point to polygon
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(fourPointsElementSet.ElementType, gridElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, gridElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 3));
// polygon to point
methodDescriptions = elementMapper.GetAvailableMethods(gridElementSet.ElementType, fourPointsElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), gridElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 2));
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_Pointline()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Other Points");
ElementSet lineElementSet = CreateElementSet("2 element XYline ElementSet");
ElementMapper elementMapper = new ElementMapper();
// point to polyline
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(fourPointsElementSet.ElementType, lineElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, lineElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
elementMapper.Initialise(methodDescriptions[1].ToString(), fourPointsElementSet, lineElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
// polyline to point
methodDescriptions = elementMapper.GetAvailableMethods(lineElementSet.ElementType, fourPointsElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), lineElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 1));
elementMapper.Initialise(methodDescriptions[1].ToString(), lineElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.585786437626905, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+ elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 1));
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_PointPoint()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Points");
ElementSet twoPointsElementSet = CreateElementSet("2 Points");
ElementMapper elementMapper = new ElementMapper();
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(fourPointsElementSet.ElementType, twoPointsElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, twoPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0)
+ elementMapper.GetValueFromMappingMatrix(0, 1)
+ elementMapper.GetValueFromMappingMatrix(0, 2)
+ elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+ elementMapper.GetValueFromMappingMatrix(1, 1)
+ elementMapper.GetValueFromMappingMatrix(1, 2)
+ elementMapper.GetValueFromMappingMatrix(1, 3));
elementMapper.Initialise(methodDescriptions[1].ToString(), fourPointsElementSet, twoPointsElementSet);
Assert.AreEqual(0.56310461156889, elementMapper.GetValueFromMappingMatrix(0, 0), 0.000000001);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0)
+ elementMapper.GetValueFromMappingMatrix(0, 1)
+ elementMapper.GetValueFromMappingMatrix(0, 2)
+ elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+ elementMapper.GetValueFromMappingMatrix(1, 1)
+ elementMapper.GetValueFromMappingMatrix(1, 2)
+ elementMapper.GetValueFromMappingMatrix(1, 3), 0.000000001);
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_PolygonPolygon()
{
Coordinate v1_0_10 = new Coordinate(0, 10, 0);
Coordinate v1_0_0 = new Coordinate(0, 0, 0);
Coordinate v1_10_0 = new Coordinate(10, 0, 0);
Coordinate v1_10_10 = new Coordinate(10, 10, 0);
Coordinate v1_20_0 = new Coordinate(20, 0, 0);
Coordinate v1_20_10 = new Coordinate(20, 10, 0);
Coordinate v1_5_9 = new Coordinate(5, 9, 0);
Coordinate v1_5_1 = new Coordinate(5, 1, 0);
Coordinate v1_15_5 = new Coordinate(15, 5, 0);
Element leftSquare = new Element("LeftSquare");
leftSquare.AddVertex(v1_0_10);
leftSquare.AddVertex(v1_0_0);
leftSquare.AddVertex(v1_10_0);
leftSquare.AddVertex(v1_10_10);
Element rightSquare = new Element("RightSquare");
rightSquare.AddVertex(v1_10_10);
rightSquare.AddVertex(v1_10_0);
rightSquare.AddVertex(v1_20_0);
rightSquare.AddVertex(v1_20_10);
Element triangle = new Element("Triangle");
triangle.AddVertex(v1_5_9);
triangle.AddVertex(v1_5_1);
triangle.AddVertex(v1_15_5);
ElementSet twoSquareElementSet = new ElementSet("TwoSquareElementSet", "TwoSquareElementSet", ElementType.Polygon);
ElementSet triangleElementSet = new ElementSet("TriangleElementSet", "TriangleElementSet", ElementType.Polygon);
twoSquareElementSet.AddElement(leftSquare);
twoSquareElementSet.AddElement(rightSquare);
triangleElementSet.AddElement(triangle);
ElementMapper elementMapper = new ElementMapper();
IIdentifiable[] methods = SpatialAdaptedOutputFactory.GetAvailableMethods(twoSquareElementSet.ElementType, triangleElementSet.ElementType);
elementMapper.Initialise(methods[0], triangleElementSet, twoSquareElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test1");
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0), "Test2");
elementMapper.Initialise(methods[0], twoSquareElementSet, triangleElementSet);
Assert.AreEqual(0.75, elementMapper.GetValueFromMappingMatrix(0, 0), 0.000000001, "Test3");
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 1), "Test4");
elementMapper.Initialise(methods[1], triangleElementSet, twoSquareElementSet);
Assert.AreEqual(0.3, elementMapper.GetValueFromMappingMatrix(0, 0), "Test5");
Assert.AreEqual(0.1, elementMapper.GetValueFromMappingMatrix(1, 0), "Test6");
elementMapper.Initialise(methods[1], twoSquareElementSet, triangleElementSet);
Assert.AreEqual(0.75, elementMapper.GetValueFromMappingMatrix(0, 0), 0.0000000001, "Test7");
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 1), "Test8");
Coordinate v2_0_2 = new Coordinate(0, 2, 0);
Coordinate v2_0_0 = new Coordinate(0, 0, 0);
Coordinate v2_2_0 = new Coordinate(2, 0, 0);
Coordinate v2_1_2 = new Coordinate(1, 2, 0);
Coordinate v2_1_0 = new Coordinate(1, 0, 0);
Coordinate v2_3_0 = new Coordinate(3, 0, 0);
Element leftTriangle2 = new Element("leftTriangle");
leftTriangle2.AddVertex(v2_0_2);
leftTriangle2.AddVertex(v2_0_0);
leftTriangle2.AddVertex(v2_2_0);
Element rightTriangle2 = new Element("rightTriangle");
rightTriangle2.AddVertex(v2_1_2);
rightTriangle2.AddVertex(v2_1_0);
rightTriangle2.AddVertex(v2_3_0);
ElementSet leftTriangleElementSet2 = new ElementSet("TwoSquareElementSet", "TwoSquareElementSet", ElementType.Polygon);
ElementSet rightTriangleElementSet2 = new ElementSet("TriangleElementSet", "TriangleElementSet", ElementType.Polygon);
leftTriangleElementSet2.AddElement(leftTriangle2);
rightTriangleElementSet2.AddElement(rightTriangle2);
elementMapper.Initialise(methods[0], leftTriangleElementSet2, rightTriangleElementSet2);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test9");
elementMapper.Initialise(methods[0], rightTriangleElementSet2, leftTriangleElementSet2);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test10");
elementMapper.Initialise(methods[1], leftTriangleElementSet2, rightTriangleElementSet2);
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 0), "Test11");
elementMapper.Initialise(methods[1], rightTriangleElementSet2, leftTriangleElementSet2);
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 0), "Test12");
Coordinate v3_0_2 = new Coordinate(0, 2, 0);
Coordinate v3_0_0 = new Coordinate(0, 0, 0);
Coordinate v3_2_0 = new Coordinate(2, 0, 0);
Coordinate v3_1_2 = new Coordinate(1, 2, 0);
Coordinate v3_1_0 = new Coordinate(1, 0, 0);
Coordinate v3_3_2 = new Coordinate(3, 2, 0);
Element leftTriangle3 = new Element("leftTriangle");
leftTriangle3.AddVertex(v3_0_2);
leftTriangle3.AddVertex(v3_0_0);
leftTriangle3.AddVertex(v3_2_0);
Element rightTriangle3 = new Element("rightTriangle");
rightTriangle3.AddVertex(v3_1_2);
rightTriangle3.AddVertex(v3_1_0);
rightTriangle3.AddVertex(v3_3_2);
ElementSet leftTriangleElementSet3 = new ElementSet("TwoSquareElementSet", "TwoSquareElementSet", ElementType.Polygon);
ElementSet rightTriangleElementSet3 = new ElementSet("TriangleElementSet", "TriangleElementSet", ElementType.Polygon);
leftTriangleElementSet3.AddElement(leftTriangle3);
rightTriangleElementSet3.AddElement(rightTriangle3);
elementMapper.Initialise(methods[0], leftTriangleElementSet3, rightTriangleElementSet3);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test13");
elementMapper.Initialise(methods[0], rightTriangleElementSet3, leftTriangleElementSet3);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test14");
elementMapper.Initialise(methods[1], leftTriangleElementSet3, rightTriangleElementSet3);
Assert.AreEqual(0.125, elementMapper.GetValueFromMappingMatrix(0, 0), "Test15");
elementMapper.Initialise(methods[1], rightTriangleElementSet3, leftTriangleElementSet3);
Assert.AreEqual(0.125, elementMapper.GetValueFromMappingMatrix(0, 0), "Test16");
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_PolylinePolygon()
{
ElementSet twoSquaresGrid = new ElementSet("TwoSquaresGrid", "TwoSquaresGrid", ElementType.Polygon);
Element e1 = new Element("e1");
Element e2 = new Element("e2");
e1.AddVertex(new Coordinate(1, 1, 0));
e1.AddVertex(new Coordinate(3, 1, 0));
e1.AddVertex(new Coordinate(3, 3, 0));
e1.AddVertex(new Coordinate(1, 3, 0));
e2.AddVertex(new Coordinate(3, 1, 0));
e2.AddVertex(new Coordinate(5, 1, 0));
e2.AddVertex(new Coordinate(5, 3, 0));
e2.AddVertex(new Coordinate(3, 3, 0));
twoSquaresGrid.AddElement(e1);
twoSquaresGrid.AddElement(e2);
ElementSet twoLines = new ElementSet("TwoLines", "TwoLines", ElementType.PolyLine);
Element l1 = new Element("l1");
Element l2 = new Element("l2");
l1.AddVertex(new Coordinate(0, 2.5, 0));
l1.AddVertex(new Coordinate(2, 2.5, 0));
l2.AddVertex(new Coordinate(2, 2.5, 0));
l2.AddVertex(new Coordinate(4, 1.5, 0));
twoLines.AddElement(l1);
twoLines.AddElement(l2);
// Line to Polygon
ElementMapper elementMapper = new ElementMapper();
IIdentifiable[] methods = SpatialAdaptedOutputFactory.GetAvailableMethods(twoLines.ElementType, twoSquaresGrid.ElementType);
elementMapper.Initialise(methods[0], twoLines, twoSquaresGrid);
Assert.AreEqual(1 / (1 + Math.Sqrt(1 + Math.Pow(0.5, 2))), elementMapper.GetValueFromMappingMatrix(0, 0), "Test1");
Assert.AreEqual(1 - 1 / (1 + Math.Sqrt(1 + Math.Pow(0.5, 2))), elementMapper.GetValueFromMappingMatrix(0, 1), "Test2");
Assert.AreEqual(0.0, elementMapper.GetValueFromMappingMatrix(1, 0), "Test3");
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 1), "Test4");
elementMapper.Initialise(methods[1], twoLines, twoSquaresGrid);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0), "Test5");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 1), "Test6");
Assert.AreEqual(0.0, elementMapper.GetValueFromMappingMatrix(1, 0), "Test7");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 1), "Test8");
// Polygon To PolyLine
methods = SpatialAdaptedOutputFactory.GetAvailableMethods(twoSquaresGrid.ElementType, twoLines.ElementType);
elementMapper.Initialise(methods[0], twoSquaresGrid, twoLines);
Assert.AreEqual(1.0, elementMapper.GetValueFromMappingMatrix(0, 0), "Test9");
Assert.AreEqual(0.0, elementMapper.GetValueFromMappingMatrix(0, 1), "Test10");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 0), "Test11");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 1), "Test12");
elementMapper.Initialise(methods[1], twoSquaresGrid, twoLines);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0), "Test13");
Assert.AreEqual(0.0, elementMapper.GetValueFromMappingMatrix(0, 1), "Test14");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 0), "Test15");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 1), "Test16");
}
public void UpdateMappingMatrix_PointPoint()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Points");
ElementSet twoPointsElementSet = CreateElementSet("2 Points");
ElementMapper elementMapper = new ElementMapper();
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(fourPointsElementSet.ElementType, twoPointsElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, twoPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0)
+elementMapper.GetValueFromMappingMatrix(0, 1)
+elementMapper.GetValueFromMappingMatrix(0, 2)
+elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+elementMapper.GetValueFromMappingMatrix(1, 1)
+elementMapper.GetValueFromMappingMatrix(1, 2)
+elementMapper.GetValueFromMappingMatrix(1, 3));
elementMapper.Initialise(methodDescriptions[1].ToString(), fourPointsElementSet, twoPointsElementSet);
Assert.AreEqual(0.56310461156889, elementMapper.GetValueFromMappingMatrix(0, 0),0.000000001);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0)
+elementMapper.GetValueFromMappingMatrix(0, 1)
+elementMapper.GetValueFromMappingMatrix(0, 2)
+elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+elementMapper.GetValueFromMappingMatrix(1, 1)
+elementMapper.GetValueFromMappingMatrix(1, 2)
+elementMapper.GetValueFromMappingMatrix(1, 3),0.000000001);
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_PointPolyline()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Other Points");
ElementSet lineElementSet = CreateElementSet("3 points polyline");
ElementMapper elementMapper = new ElementMapper();
// point to polyline
IIdentifiable[] methods = SpatialAdaptedOutputFactory.GetAvailableMethods(fourPointsElementSet.ElementType, lineElementSet.ElementType);
elementMapper.Initialise(methods[0], fourPointsElementSet, lineElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
elementMapper.Initialise(methods[0], fourPointsElementSet, lineElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
// polyline to point
methods = SpatialAdaptedOutputFactory.GetAvailableMethods(lineElementSet.ElementType, fourPointsElementSet.ElementType);
elementMapper.Initialise(methods[0], lineElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 1));
elementMapper.Initialise(methods[1], lineElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.585786437626905, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+ elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 1));
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_PointPoint()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Points");
ElementSet twoPointsElementSet = CreateElementSet("2 Points");
ElementMapper elementMapper = new ElementMapper();
IIdentifiable[] methods = SpatialAdaptedOutputFactory.GetAvailableMethods(fourPointsElementSet.ElementType, twoPointsElementSet.ElementType);
elementMapper.Initialise(methods[0], fourPointsElementSet, twoPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0)
+ elementMapper.GetValueFromMappingMatrix(0, 1)
+ elementMapper.GetValueFromMappingMatrix(0, 2)
+ elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+ elementMapper.GetValueFromMappingMatrix(1, 1)
+ elementMapper.GetValueFromMappingMatrix(1, 2)
+ elementMapper.GetValueFromMappingMatrix(1, 3));
elementMapper.Initialise(methods[1], fourPointsElementSet, twoPointsElementSet);
Assert.AreEqual(0.56310461156889, elementMapper.GetValueFromMappingMatrix(0, 0), 0.000000001);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0)
+ elementMapper.GetValueFromMappingMatrix(0, 1)
+ elementMapper.GetValueFromMappingMatrix(0, 2)
+ elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+ elementMapper.GetValueFromMappingMatrix(1, 1)
+ elementMapper.GetValueFromMappingMatrix(1, 2)
+ elementMapper.GetValueFromMappingMatrix(1, 3), 0.000000001);
}
public void UpdateMappingMatrix_PolylinePolygon()
{
ElementSet twoSquaresGrid = new ElementSet("TwoSquaresGrid","TwoSquaresGrid",ElementType.XYPolygon,new SpatialReference("ref"));
Element e1 = new Element("e1");
Element e2 = new Element("e2");
e1.AddVertex(new Vertex(1,1,0));
e1.AddVertex(new Vertex(3,1,0));
e1.AddVertex(new Vertex(3,3,0));
e1.AddVertex(new Vertex(1,3,0));
e2.AddVertex(new Vertex(3,1,0));
e2.AddVertex(new Vertex(5,1,0));
e2.AddVertex(new Vertex(5,3,0));
e2.AddVertex(new Vertex(3,3,0));
twoSquaresGrid.AddElement(e1);
twoSquaresGrid.AddElement(e2);
ElementSet twoLines = new ElementSet("TwoLines","TwoLines",ElementType.XYPolyLine,new SpatialReference("ref"));
Element l1 = new Element("l1");
Element l2 = new Element("l2");
l1.AddVertex(new Vertex(0,2.5,0));
l1.AddVertex(new Vertex(2,2.5,0));
l2.AddVertex(new Vertex(2,2.5,0));
l2.AddVertex(new Vertex(4,1.5,0));
twoLines.AddElement(l1);
twoLines.AddElement(l2);
// Line to Polygon
ElementMapper elementMapper = new ElementMapper();
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(twoLines.ElementType, twoSquaresGrid.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), twoLines,twoSquaresGrid);
Assert.AreEqual(1/(1+Math.Sqrt(1+Math.Pow(0.5,2))),elementMapper.GetValueFromMappingMatrix(0,0),"Test1");
Assert.AreEqual(1-1/(1+Math.Sqrt(1+Math.Pow(0.5,2))),elementMapper.GetValueFromMappingMatrix(0,1),"Test2");
Assert.AreEqual(0.0,elementMapper.GetValueFromMappingMatrix(1,0),"Test3");
Assert.AreEqual(1,elementMapper.GetValueFromMappingMatrix(1,1),"Test4");
elementMapper.Initialise(methodDescriptions[1].ToString(), twoLines,twoSquaresGrid);
Assert.AreEqual(0.5,elementMapper.GetValueFromMappingMatrix(0,0),"Test5");
Assert.AreEqual(0.5,elementMapper.GetValueFromMappingMatrix(0,1),"Test6");
Assert.AreEqual(0.0,elementMapper.GetValueFromMappingMatrix(1,0),"Test7");
Assert.AreEqual(0.5,elementMapper.GetValueFromMappingMatrix(1,1),"Test8");
// Polygon To PolyLine
methodDescriptions = elementMapper.GetAvailableMethods(twoSquaresGrid.ElementType, twoLines.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), twoSquaresGrid, twoLines);
Assert.AreEqual(1.0,elementMapper.GetValueFromMappingMatrix(0,0),"Test9");
Assert.AreEqual(0.0,elementMapper.GetValueFromMappingMatrix(0,1),"Test10");
Assert.AreEqual(0.5,elementMapper.GetValueFromMappingMatrix(1,0),"Test11");
Assert.AreEqual(0.5,elementMapper.GetValueFromMappingMatrix(1,1),"Test12");
elementMapper.Initialise(methodDescriptions[1].ToString(), twoSquaresGrid, twoLines);
Assert.AreEqual(0.5,elementMapper.GetValueFromMappingMatrix(0,0),"Test13");
Assert.AreEqual(0.0,elementMapper.GetValueFromMappingMatrix(0,1),"Test14");
Assert.AreEqual(0.5,elementMapper.GetValueFromMappingMatrix(1,0),"Test15");
Assert.AreEqual(0.5,elementMapper.GetValueFromMappingMatrix(1,1),"Test16");
}
public void UpdateMappingMatrix_PolygonPolygon()
{
Vertex v1_0_10 = new Vertex(0,10,0);
Vertex v1_0_0 = new Vertex(0,0,0);
Vertex v1_10_0 = new Vertex(10,0,0);
Vertex v1_10_10 = new Vertex(10,10,0);
Vertex v1_20_0 = new Vertex(20,0,0);
Vertex v1_20_10 = new Vertex(20,10,0);
Vertex v1_5_9 = new Vertex(5,9,0);
Vertex v1_5_1 = new Vertex(5,1,0);
Vertex v1_15_5 = new Vertex(15,5,0);
Element LeftSquare = new Element("LeftSquare");
LeftSquare.AddVertex(v1_0_10);
LeftSquare.AddVertex(v1_0_0);
LeftSquare.AddVertex(v1_10_0);
LeftSquare.AddVertex(v1_10_10);
Element RightSquare = new Element("RightSquare");
RightSquare.AddVertex(v1_10_10);
RightSquare.AddVertex(v1_10_0);
RightSquare.AddVertex(v1_20_0);
RightSquare.AddVertex(v1_20_10);
Element Triangle = new Element("Triangle");
Triangle.AddVertex(v1_5_9);
Triangle.AddVertex(v1_5_1);
Triangle.AddVertex(v1_15_5);
ElementSet TwoSquareElementSet = new ElementSet("TwoSquareElementSet","TwoSquareElementSet",ElementType.XYPolygon,new SpatialReference("ref"));
ElementSet TriangleElementSet = new ElementSet("TriangleElementSet","TriangleElementSet",ElementType.XYPolygon,new SpatialReference("ref"));
TwoSquareElementSet.AddElement(LeftSquare);
TwoSquareElementSet.AddElement(RightSquare);
TriangleElementSet.AddElement(Triangle);
ElementMapper elementMapper = new ElementMapper();
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(TwoSquareElementSet.ElementType, TriangleElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), TriangleElementSet, TwoSquareElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0),"Test1");
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0),"Test2");
elementMapper.Initialise(methodDescriptions[0].ToString(), TwoSquareElementSet, TriangleElementSet);
Assert.AreEqual(0.75, elementMapper.GetValueFromMappingMatrix(0, 0),0.000000001,"Test3");
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 1),"Test4");
elementMapper.Initialise(methodDescriptions[1].ToString(), TriangleElementSet, TwoSquareElementSet);
Assert.AreEqual(0.3, elementMapper.GetValueFromMappingMatrix(0, 0),"Test5");
Assert.AreEqual(0.1, elementMapper.GetValueFromMappingMatrix(1, 0),"Test6");
elementMapper.Initialise(methodDescriptions[1].ToString(), TwoSquareElementSet, TriangleElementSet);
Assert.AreEqual(0.75, elementMapper.GetValueFromMappingMatrix(0, 0),0.0000000001,"Test7");
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 1),"Test8");
Vertex v2_0_2 = new Vertex(0,2,0);
Vertex v2_0_0 = new Vertex(0,0,0);
Vertex v2_2_0 = new Vertex(2,0,0);
Vertex v2_1_2 = new Vertex(1,2,0);
Vertex v2_1_0 = new Vertex(1,0,0);
Vertex v2_3_0 = new Vertex(3,0,0);
Element LeftTriangle2 = new Element("LeftTriangle");
LeftTriangle2.AddVertex(v2_0_2);
LeftTriangle2.AddVertex(v2_0_0);
LeftTriangle2.AddVertex(v2_2_0);
Element RightTriangle2 = new Element("RightTriangle");
RightTriangle2.AddVertex(v2_1_2);
RightTriangle2.AddVertex(v2_1_0);
RightTriangle2.AddVertex(v2_3_0);
ElementSet LeftTriangleElementSet2 = new ElementSet("TwoSquareElementSet","TwoSquareElementSet",ElementType.XYPolygon,new SpatialReference("ref"));
ElementSet RightTriangleElementSet2 = new ElementSet("TriangleElementSet","TriangleElementSet",ElementType.XYPolygon,new SpatialReference("ref"));
LeftTriangleElementSet2.AddElement(LeftTriangle2);
RightTriangleElementSet2.AddElement(RightTriangle2);
elementMapper.Initialise(methodDescriptions[0].ToString(), LeftTriangleElementSet2, RightTriangleElementSet2);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0),"Test9");
elementMapper.Initialise(methodDescriptions[0].ToString(), RightTriangleElementSet2, LeftTriangleElementSet2);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0),"Test10");
elementMapper.Initialise(methodDescriptions[1].ToString(), LeftTriangleElementSet2, RightTriangleElementSet2);
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 0),"Test11");
elementMapper.Initialise(methodDescriptions[1].ToString(), RightTriangleElementSet2, LeftTriangleElementSet2);
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 0),"Test12");
Vertex v3_0_2 = new Vertex(0,2,0);
Vertex v3_0_0 = new Vertex(0,0,0);
Vertex v3_2_0 = new Vertex(2,0,0);
Vertex v3_1_2 = new Vertex(1,2,0);
Vertex v3_1_0 = new Vertex(1,0,0);
Vertex v3_3_2 = new Vertex(3,2,0);
Element LeftTriangle3 = new Element("LeftTriangle");
LeftTriangle3.AddVertex(v3_0_2);
LeftTriangle3.AddVertex(v3_0_0);
LeftTriangle3.AddVertex(v3_2_0);
Element RightTriangle3 = new Element("RightTriangle");
RightTriangle3.AddVertex(v3_1_2);
RightTriangle3.AddVertex(v3_1_0);
RightTriangle3.AddVertex(v3_3_2);
ElementSet LeftTriangleElementSet3 = new ElementSet("TwoSquareElementSet","TwoSquareElementSet",ElementType.XYPolygon,new SpatialReference("ref"));
ElementSet RightTriangleElementSet3 = new ElementSet("TriangleElementSet","TriangleElementSet",ElementType.XYPolygon,new SpatialReference("ref"));
LeftTriangleElementSet3.AddElement(LeftTriangle3);
RightTriangleElementSet3.AddElement(RightTriangle3);
elementMapper.Initialise(methodDescriptions[0].ToString(), LeftTriangleElementSet3, RightTriangleElementSet3);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0),"Test13");
elementMapper.Initialise(methodDescriptions[0].ToString(), RightTriangleElementSet3, LeftTriangleElementSet3);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0),"Test14");
elementMapper.Initialise(methodDescriptions[1].ToString(), LeftTriangleElementSet3, RightTriangleElementSet3);
Assert.AreEqual(0.125, elementMapper.GetValueFromMappingMatrix(0, 0),"Test15");
elementMapper.Initialise(methodDescriptions[1].ToString(), RightTriangleElementSet3, LeftTriangleElementSet3);
Assert.AreEqual(0.125, elementMapper.GetValueFromMappingMatrix(0, 0),"Test16");
}
public void UpdateMappingMatrix_PointPolyline()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Other Points");
ElementSet lineElementSet = CreateElementSet("3 points polyline");
ElementMapper elementMapper = new ElementMapper();
// point to polyline
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(fourPointsElementSet.ElementType, lineElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, lineElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, lineElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
// polyline to point
methodDescriptions = elementMapper.GetAvailableMethods(lineElementSet.ElementType, fourPointsElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), lineElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 1));
elementMapper.Initialise(methodDescriptions[1].ToString(), lineElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0.585786437626905, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0)
+ elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(3, 1));
}
public void UpdateMappingMatrix_PointPolygon()
{
ElementSet gridElementSet = new ElementSet("gridElm","G1",ElementType.XYPolygon, new SpatialReference("ref"));
ElementSet fourPointsElementSet = new ElementSet("4 points","4P",ElementType.XYPoint,new SpatialReference("DummyID"));
Vertex v_0_20 = new Vertex(0,20,0);
Vertex v_0_10 = new Vertex(0,10,0);
Vertex v_0_0 = new Vertex(0, 0,0);
Vertex v_0_15 = new Vertex(0,15,0);
Vertex v_5_15 = new Vertex(5,15,0);
Vertex v_10_20 = new Vertex(10,20,0);
Vertex v_10_15 = new Vertex(10,15,0);
Vertex v_10_10 = new Vertex(10,10,0);
Vertex v_10_0 = new Vertex(10, 0,0);
Vertex v_15_15 = new Vertex(15,15,0);
Vertex v_15_5 = new Vertex(15,5,0);
Vertex v_20_20 = new Vertex(20,20,0);
Vertex v_20_10 = new Vertex(20,10,0);
Element square1 = new Element("square1");
Element square2 = new Element("square2");
Element square3 = new Element("square3");
square1.AddVertex(v_0_20);
square1.AddVertex(v_0_10);
square1.AddVertex(v_10_10);
square1.AddVertex(v_10_20);
square2.AddVertex(v_10_20);
square2.AddVertex(v_10_10);
square2.AddVertex(v_20_10);
square2.AddVertex(v_20_20);
square3.AddVertex(v_0_10);
square3.AddVertex(v_0_0);
square3.AddVertex(v_10_0);
square3.AddVertex(v_10_10);
gridElementSet.AddElement(square1);
gridElementSet.AddElement(square2);
gridElementSet.AddElement(square3);
Element point_5_15 = new Element("point 5, 15");
Element point_10_15 = new Element("point 10, 15");
Element point_15_15 = new Element("point 15, 15");
Element point_15_5 = new Element("point 15, 5");
point_5_15.AddVertex(v_5_15);
point_10_15.AddVertex(v_10_15);
point_15_15.AddVertex(v_15_15);
point_15_5.AddVertex(v_15_5);
fourPointsElementSet.AddElement(point_5_15);
fourPointsElementSet.AddElement(point_10_15);
fourPointsElementSet.AddElement(point_15_15);
fourPointsElementSet.AddElement(point_15_5);
ElementMapper elementMapper = new ElementMapper();
// point to polygon
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(fourPointsElementSet.ElementType, gridElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, gridElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 3));
// polygon to point
methodDescriptions = elementMapper.GetAvailableMethods(gridElementSet.ElementType, fourPointsElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), gridElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 2));
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_LinePolygon()
{
ElementSet twoSquaresGrid = new ElementSet("TwoSquaresGrid", "TwoSquaresGrid", ElementType.XYPolygon, new SpatialReference("ref"));
Element e1 = new Element("e1");
Element e2 = new Element("e2");
e1.AddVertex(new Vertex(1, 1, 0));
e1.AddVertex(new Vertex(3, 1, 0));
e1.AddVertex(new Vertex(3, 3, 0));
e1.AddVertex(new Vertex(1, 3, 0));
e2.AddVertex(new Vertex(3, 1, 0));
e2.AddVertex(new Vertex(5, 1, 0));
e2.AddVertex(new Vertex(5, 3, 0));
e2.AddVertex(new Vertex(3, 3, 0));
twoSquaresGrid.AddElement(e1);
twoSquaresGrid.AddElement(e2);
ElementSet twoLines = new ElementSet("TwoLines", "TwoLines", ElementType.XYLine, new SpatialReference("ref"));
Element l1 = new Element("l1");
Element l2 = new Element("l2");
l1.AddVertex(new Vertex(0, 2.5, 0));
l1.AddVertex(new Vertex(2, 2.5, 0));
l2.AddVertex(new Vertex(2, 2.5, 0));
l2.AddVertex(new Vertex(4, 1.5, 0));
twoLines.AddElement(l1);
twoLines.AddElement(l2);
// Line to Polygon
ElementMapper elementMapper = new ElementMapper();
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(twoLines.ElementType, twoSquaresGrid.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), twoLines, twoSquaresGrid);
Assert.AreEqual(1 / (1 + Math.Sqrt(1 + Math.Pow(0.5, 2))), elementMapper.GetValueFromMappingMatrix(0, 0), "Test1");
Assert.AreEqual(1 - 1 / (1 + Math.Sqrt(1 + Math.Pow(0.5, 2))), elementMapper.GetValueFromMappingMatrix(0, 1), "Test2");
Assert.AreEqual(0.0, elementMapper.GetValueFromMappingMatrix(1, 0), "Test3");
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 1), "Test4");
elementMapper.Initialise(methodDescriptions[1].ToString(), twoLines, twoSquaresGrid);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0), "Test5");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 1), "Test6");
Assert.AreEqual(0.0, elementMapper.GetValueFromMappingMatrix(1, 0), "Test7");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 1), "Test8");
// Polygon To PolyLine
methodDescriptions = elementMapper.GetAvailableMethods(twoSquaresGrid.ElementType, twoLines.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), twoSquaresGrid, twoLines);
Assert.AreEqual(1.0, elementMapper.GetValueFromMappingMatrix(0, 0), "Test9");
Assert.AreEqual(0.0, elementMapper.GetValueFromMappingMatrix(0, 1), "Test10");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 0), "Test11");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 1), "Test12");
elementMapper.Initialise(methodDescriptions[1].ToString(), twoSquaresGrid, twoLines);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0), "Test13");
Assert.AreEqual(0.0, elementMapper.GetValueFromMappingMatrix(0, 1), "Test14");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 0), "Test15");
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(1, 1), "Test16");
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_PointPolygon()
{
ElementSet gridElementSet = new ElementSet("gridElm", "G1", ElementType.Polygon);
ElementSet fourPointsElementSet = new ElementSet("4 points", "4P", ElementType.Point);
Coordinate v_0_20 = new Coordinate(0, 20, 0);
Coordinate v_0_10 = new Coordinate(0, 10, 0);
Coordinate v_0_0 = new Coordinate(0, 0, 0);
Coordinate v_5_15 = new Coordinate(5, 15, 0);
Coordinate v_10_20 = new Coordinate(10, 20, 0);
Coordinate v_10_15 = new Coordinate(10, 15, 0);
Coordinate v_10_10 = new Coordinate(10, 10, 0);
Coordinate v_10_0 = new Coordinate(10, 0, 0);
Coordinate v_15_15 = new Coordinate(15, 15, 0);
Coordinate v_15_5 = new Coordinate(15, 5, 0);
Coordinate v_20_20 = new Coordinate(20, 20, 0);
Coordinate v_20_10 = new Coordinate(20, 10, 0);
Element square1 = new Element("square1");
Element square2 = new Element("square2");
Element square3 = new Element("square3");
square1.AddVertex(v_0_20);
square1.AddVertex(v_0_10);
square1.AddVertex(v_10_10);
square1.AddVertex(v_10_20);
square2.AddVertex(v_10_20);
square2.AddVertex(v_10_10);
square2.AddVertex(v_20_10);
square2.AddVertex(v_20_20);
square3.AddVertex(v_0_10);
square3.AddVertex(v_0_0);
square3.AddVertex(v_10_0);
square3.AddVertex(v_10_10);
gridElementSet.AddElement(square1);
gridElementSet.AddElement(square2);
gridElementSet.AddElement(square3);
Element point_5_15 = new Element("point 5, 15");
Element point_10_15 = new Element("point 10, 15");
Element point_15_15 = new Element("point 15, 15");
Element point_15_5 = new Element("point 15, 5");
point_5_15.AddVertex(v_5_15);
point_10_15.AddVertex(v_10_15);
point_15_15.AddVertex(v_15_15);
point_15_5.AddVertex(v_15_5);
fourPointsElementSet.AddElement(point_5_15);
fourPointsElementSet.AddElement(point_10_15);
fourPointsElementSet.AddElement(point_15_15);
fourPointsElementSet.AddElement(point_15_5);
ElementMapper elementMapper = new ElementMapper();
// point to polygon
IIdentifiable[] methods = SpatialAdaptedOutputFactory.GetAvailableMethods(fourPointsElementSet.ElementType, gridElementSet.ElementType);
elementMapper.Initialise(methods[0], fourPointsElementSet, gridElementSet);
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0.5, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 3));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 3));
// polygon to point
methods = SpatialAdaptedOutputFactory.GetAvailableMethods(gridElementSet.ElementType, fourPointsElementSet.ElementType);
elementMapper.Initialise(methods[0], gridElementSet, fourPointsElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(0, 2));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(1, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 0));
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(2, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(2, 2));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 0));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 1));
Assert.AreEqual(0, elementMapper.GetValueFromMappingMatrix(3, 2));
}
[Test] // testing the Initialise method
public void UpdateMappingMatrix_PolygonPolygon()
{
Vertex v1_0_10 = new Vertex(0, 10, 0);
Vertex v1_0_0 = new Vertex(0, 0, 0);
Vertex v1_10_0 = new Vertex(10, 0, 0);
Vertex v1_10_10 = new Vertex(10, 10, 0);
Vertex v1_20_0 = new Vertex(20, 0, 0);
Vertex v1_20_10 = new Vertex(20, 10, 0);
Vertex v1_5_9 = new Vertex(5, 9, 0);
Vertex v1_5_1 = new Vertex(5, 1, 0);
Vertex v1_15_5 = new Vertex(15, 5, 0);
Element LeftSquare = new Element("LeftSquare");
LeftSquare.AddVertex(v1_0_10);
LeftSquare.AddVertex(v1_0_0);
LeftSquare.AddVertex(v1_10_0);
LeftSquare.AddVertex(v1_10_10);
Element RightSquare = new Element("RightSquare");
RightSquare.AddVertex(v1_10_10);
RightSquare.AddVertex(v1_10_0);
RightSquare.AddVertex(v1_20_0);
RightSquare.AddVertex(v1_20_10);
Element Triangle = new Element("Triangle");
Triangle.AddVertex(v1_5_9);
Triangle.AddVertex(v1_5_1);
Triangle.AddVertex(v1_15_5);
ElementSet TwoSquareElementSet = new ElementSet("TwoSquareElementSet", "TwoSquareElementSet", ElementType.XYPolygon, new SpatialReference("ref"));
ElementSet TriangleElementSet = new ElementSet("TriangleElementSet", "TriangleElementSet", ElementType.XYPolygon, new SpatialReference("ref"));
TwoSquareElementSet.AddElement(LeftSquare);
TwoSquareElementSet.AddElement(RightSquare);
TriangleElementSet.AddElement(Triangle);
ElementMapper elementMapper = new ElementMapper();
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(TwoSquareElementSet.ElementType, TriangleElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), TriangleElementSet, TwoSquareElementSet);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test1");
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(1, 0), "Test2");
elementMapper.Initialise(methodDescriptions[0].ToString(), TwoSquareElementSet, TriangleElementSet);
Assert.AreEqual(0.75, elementMapper.GetValueFromMappingMatrix(0, 0), 0.000000001, "Test3");
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 1), "Test4");
elementMapper.Initialise(methodDescriptions[1].ToString(), TriangleElementSet, TwoSquareElementSet);
Assert.AreEqual(0.3, elementMapper.GetValueFromMappingMatrix(0, 0), "Test5");
Assert.AreEqual(0.1, elementMapper.GetValueFromMappingMatrix(1, 0), "Test6");
elementMapper.Initialise(methodDescriptions[1].ToString(), TwoSquareElementSet, TriangleElementSet);
Assert.AreEqual(0.75, elementMapper.GetValueFromMappingMatrix(0, 0), 0.0000000001, "Test7");
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 1), "Test8");
Vertex v2_0_2 = new Vertex(0, 2, 0);
Vertex v2_0_0 = new Vertex(0, 0, 0);
Vertex v2_2_0 = new Vertex(2, 0, 0);
Vertex v2_1_2 = new Vertex(1, 2, 0);
Vertex v2_1_0 = new Vertex(1, 0, 0);
Vertex v2_3_0 = new Vertex(3, 0, 0);
Element LeftTriangle2 = new Element("LeftTriangle");
LeftTriangle2.AddVertex(v2_0_2);
LeftTriangle2.AddVertex(v2_0_0);
LeftTriangle2.AddVertex(v2_2_0);
Element RightTriangle2 = new Element("RightTriangle");
RightTriangle2.AddVertex(v2_1_2);
RightTriangle2.AddVertex(v2_1_0);
RightTriangle2.AddVertex(v2_3_0);
ElementSet LeftTriangleElementSet2 = new ElementSet("TwoSquareElementSet", "TwoSquareElementSet", ElementType.XYPolygon, new SpatialReference("ref"));
ElementSet RightTriangleElementSet2 = new ElementSet("TriangleElementSet", "TriangleElementSet", ElementType.XYPolygon, new SpatialReference("ref"));
LeftTriangleElementSet2.AddElement(LeftTriangle2);
RightTriangleElementSet2.AddElement(RightTriangle2);
elementMapper.Initialise(methodDescriptions[0].ToString(), LeftTriangleElementSet2, RightTriangleElementSet2);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test9");
elementMapper.Initialise(methodDescriptions[0].ToString(), RightTriangleElementSet2, LeftTriangleElementSet2);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test10");
elementMapper.Initialise(methodDescriptions[1].ToString(), LeftTriangleElementSet2, RightTriangleElementSet2);
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 0), "Test11");
elementMapper.Initialise(methodDescriptions[1].ToString(), RightTriangleElementSet2, LeftTriangleElementSet2);
Assert.AreEqual(0.25, elementMapper.GetValueFromMappingMatrix(0, 0), "Test12");
Vertex v3_0_2 = new Vertex(0, 2, 0);
Vertex v3_0_0 = new Vertex(0, 0, 0);
Vertex v3_2_0 = new Vertex(2, 0, 0);
Vertex v3_1_2 = new Vertex(1, 2, 0);
Vertex v3_1_0 = new Vertex(1, 0, 0);
Vertex v3_3_2 = new Vertex(3, 2, 0);
Element LeftTriangle3 = new Element("LeftTriangle");
LeftTriangle3.AddVertex(v3_0_2);
LeftTriangle3.AddVertex(v3_0_0);
LeftTriangle3.AddVertex(v3_2_0);
Element RightTriangle3 = new Element("RightTriangle");
RightTriangle3.AddVertex(v3_1_2);
RightTriangle3.AddVertex(v3_1_0);
RightTriangle3.AddVertex(v3_3_2);
ElementSet LeftTriangleElementSet3 = new ElementSet("TwoSquareElementSet", "TwoSquareElementSet", ElementType.XYPolygon, new SpatialReference("ref"));
ElementSet RightTriangleElementSet3 = new ElementSet("TriangleElementSet", "TriangleElementSet", ElementType.XYPolygon, new SpatialReference("ref"));
LeftTriangleElementSet3.AddElement(LeftTriangle3);
RightTriangleElementSet3.AddElement(RightTriangle3);
elementMapper.Initialise(methodDescriptions[0].ToString(), LeftTriangleElementSet3, RightTriangleElementSet3);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test13");
elementMapper.Initialise(methodDescriptions[0].ToString(), RightTriangleElementSet3, LeftTriangleElementSet3);
Assert.AreEqual(1, elementMapper.GetValueFromMappingMatrix(0, 0), "Test14");
elementMapper.Initialise(methodDescriptions[1].ToString(), LeftTriangleElementSet3, RightTriangleElementSet3);
Assert.AreEqual(0.125, elementMapper.GetValueFromMappingMatrix(0, 0), "Test15");
elementMapper.Initialise(methodDescriptions[1].ToString(), RightTriangleElementSet3, LeftTriangleElementSet3);
Assert.AreEqual(0.125, elementMapper.GetValueFromMappingMatrix(0, 0), "Test16");
}
public void Initialise()
{
ElementSet fourPointsElementSet = CreateElementSet("4 Points");
ElementSet twoPointsElementSet = CreateElementSet("2 Points");
ElementMapper elementMapper = new ElementMapper();
ArrayList methodDescriptions = elementMapper.GetAvailableMethods(fourPointsElementSet.ElementType, twoPointsElementSet.ElementType);
elementMapper.Initialise(methodDescriptions[0].ToString(), fourPointsElementSet, twoPointsElementSet);
double Calculated = elementMapper.GetValueFromMappingMatrix(0, 0);
double Expected = 1;
Assert.AreEqual(Expected,Calculated);
}
