/// <summary>
/// Tests execution for multiple reference systems by using points.
/// </summary>
/// <param name="metadataPreservation">Indicates whether the metadata should be preserved.</param>
private void TestExecuteForReferenceSystems(Boolean metadataPreservation)
{
// projected to projected
Coordinate sourceCoordinate = new Coordinate(10, 10);
GeoCoordinate sourceGeoCoordinate = ProjectedCoordinateReferenceSystems.HD72_EOV.Projection.Reverse(sourceCoordinate);
GeoCoordinate expectedGeoCoordinate = GeographicTransformations.HD72_WGS84_V1.Forward(sourceGeoCoordinate);
Coordinate expectedCoordinate = ProjectedCoordinateReferenceSystems.WGS84_WorldMercator.Projection.Forward(expectedGeoCoordinate);
IPoint sourcePoint = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreatePoint(10, 10);
IDictionary <OperationParameter, Object> parameters = new Dictionary <OperationParameter, Object>();
parameters.Add(ReferenceOperationParameters.TargetReferenceSystem, ProjectedCoordinateReferenceSystems.WGS84_WorldMercator);
parameters.Add(CommonOperationParameters.MetadataPreservation, metadataPreservation);
ReferenceTransformation transformation = new ReferenceTransformation(sourcePoint, parameters);
transformation.Execute();
IPoint resultPoint = transformation.Result as IPoint;
Assert.AreEqual(expectedCoordinate, resultPoint.Coordinate);
Assert.AreEqual(sourcePoint.Metadata, resultPoint.Metadata);
// geographic to projected
sourceGeoCoordinate = new GeoCoordinate(Angle.FromDegree(45), Angle.FromDegree(45));
expectedGeoCoordinate = GeographicTransformations.HD72_WGS84_V1.Forward(sourceGeoCoordinate);
expectedCoordinate = ProjectedCoordinateReferenceSystems.WGS84_WorldMercator.Projection.Forward(expectedGeoCoordinate);
sourcePoint = new GeometryFactory(Geographic2DCoordinateReferenceSystems.HD72).CreatePoint(45, 45);
parameters = new Dictionary <OperationParameter, Object>();
parameters.Add(ReferenceOperationParameters.TargetReferenceSystem, ProjectedCoordinateReferenceSystems.WGS84_WorldMercator);
transformation = new ReferenceTransformation(sourcePoint, parameters);
transformation.Execute();
resultPoint = transformation.Result as IPoint;
Assert.AreEqual(expectedCoordinate, resultPoint.Coordinate);
// projected to geographic
sourceCoordinate = new Coordinate(10, 10);
sourceGeoCoordinate = ProjectedCoordinateReferenceSystems.HD72_EOV.Projection.Reverse(sourceCoordinate);
expectedGeoCoordinate = GeographicTransformations.HD72_WGS84_V1.Forward(sourceGeoCoordinate);
sourcePoint = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreatePoint(sourceCoordinate);
parameters = new Dictionary <OperationParameter, Object>();
parameters.Add(ReferenceOperationParameters.TargetReferenceSystem, Geographic2DCoordinateReferenceSystems.WGS84);
transformation = new ReferenceTransformation(sourcePoint, parameters);
transformation.Execute();
resultPoint = transformation.Result as IPoint;
Assert.AreEqual(expectedGeoCoordinate.Latitude.GetValue(UnitsOfMeasurement.Degree), resultPoint.Coordinate.X);
Assert.AreEqual(expectedGeoCoordinate.Longitude.GetValue(UnitsOfMeasurement.Degree), resultPoint.Coordinate.Y);
// projected to projected (reverse)
sourceCoordinate = new Coordinate(10, 10);
sourceGeoCoordinate = ProjectedCoordinateReferenceSystems.WGS84_WorldMercator.Projection.Reverse(sourceCoordinate);
expectedGeoCoordinate = GeographicTransformations.HD72_WGS84_V1.Reverse(sourceGeoCoordinate);
expectedCoordinate = ProjectedCoordinateReferenceSystems.HD72_EOV.Projection.Forward(expectedGeoCoordinate);
sourcePoint = new GeometryFactory(ProjectedCoordinateReferenceSystems.WGS84_WorldMercator).CreatePoint(10, 10);
parameters = new Dictionary <OperationParameter, Object>();
parameters.Add(ReferenceOperationParameters.TargetReferenceSystem, ProjectedCoordinateReferenceSystems.HD72_EOV);
transformation = new ReferenceTransformation(sourcePoint, parameters);
transformation.Execute();
resultPoint = transformation.Result as IPoint;
Assert.AreEqual(expectedCoordinate, resultPoint.Coordinate);
// same reference system
sourcePoint = new GeometryFactory(ProjectedCoordinateReferenceSystems.WGS84_WorldMercator).CreatePoint(10, 10);
parameters = new Dictionary <OperationParameter, Object>();
parameters.Add(ReferenceOperationParameters.TargetReferenceSystem, ProjectedCoordinateReferenceSystems.WGS84_WorldMercator);
transformation = new ReferenceTransformation(sourcePoint, parameters);
transformation.Execute();
resultPoint = transformation.Result as IPoint;
Assert.AreEqual(sourcePoint.Coordinate, resultPoint.Coordinate);
// no source reference system
sourcePoint = new GeometryFactory().CreatePoint(10, 10);
parameters = new Dictionary <OperationParameter, Object>();
parameters.Add(ReferenceOperationParameters.TargetReferenceSystem, ProjectedCoordinateReferenceSystems.WGS84_WorldMercator);
transformation = new ReferenceTransformation(sourcePoint, parameters);
transformation.Execute();
resultPoint = transformation.Result as IPoint;
Assert.AreEqual(sourcePoint.Coordinate, resultPoint.Coordinate);
// no target reference system
parameters = new Dictionary <OperationParameter, Object>();
parameters.Add(ReferenceOperationParameters.TargetReferenceSystem, null);
Assert.Throws <ArgumentException>(() => transformation = new ReferenceTransformation(sourcePoint, parameters));
}
/// <summary>
/// Tests execution for multiple geometry types.
/// </summary>
/// <param name="metadataPreservation">Indicates whether the metadata should be preserved.</param>
private void TestExecuteForGeometries(Boolean metadataPreservation)
{
Random random = new Random(0);
Coordinate[] sourceCoordinates = Enumerable.Range(0, 100).Select(value => random.Next(0, 10000)).Select(value => new Coordinate(value, value)).ToArray();
Coordinate[] expectedCoordinates = new Coordinate[sourceCoordinates.Length];
for (Int32 i = 0; i < sourceCoordinates.Length; i++)
{
GeoCoordinate sourceGeoCoordinate = ProjectedCoordinateReferenceSystems.HD72_EOV.Projection.Reverse(sourceCoordinates[i]);
GeoCoordinate expectedGeoCoordinate = GeographicTransformations.HD72_WGS84_V1.Forward(sourceGeoCoordinate);
expectedCoordinates[i] = ProjectedCoordinateReferenceSystems.WGS84_WorldMercator.Projection.Forward(expectedGeoCoordinate);
}
// lines
ILine sourceLine = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreateLine(sourceCoordinates[0], sourceCoordinates[1]);
IDictionary <OperationParameter, Object> parameters = new Dictionary <OperationParameter, Object>();
parameters.Add(ReferenceOperationParameters.TargetReferenceSystem, ProjectedCoordinateReferenceSystems.WGS84_WorldMercator);
parameters.Add(CommonOperationParameters.MetadataPreservation, metadataPreservation);
ReferenceTransformation transformation = new ReferenceTransformation(sourceLine, parameters);
transformation.Execute();
Assert.IsInstanceOf <ILine>(transformation.Result);
Assert.AreEqual(sourceLine.Metadata, transformation.Result.Metadata);
ILine resultLine = transformation.Result as ILine;
for (Int32 i = 0; i < resultLine.Count; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultLine.GetCoordinate(i));
}
// test case 2: line strings
ILineString sourceLineString = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreateLineString(sourceCoordinates);
transformation = new ReferenceTransformation(sourceLineString, parameters);
transformation.Execute();
Assert.IsInstanceOf <ILineString>(transformation.Result);
Assert.AreEqual(sourceLineString.Metadata, transformation.Result.Metadata);
ILineString resultLineString = transformation.Result as ILineString;
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultLineString.GetCoordinate(i));
}
// linear rings
ILinearRing sourceLinearRing = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreateLinearRing(sourceCoordinates);
transformation = new ReferenceTransformation(sourceLinearRing, parameters);
transformation.Execute();
Assert.IsInstanceOf <ILinearRing>(transformation.Result);
Assert.AreEqual(sourceLinearRing.Metadata, transformation.Result.Metadata);
ILinearRing resultLinearRing = transformation.Result as ILinearRing;
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultLinearRing.GetCoordinate(i));
}
// polygon without hole
IPolygon sourcePolygon = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreatePolygon(sourceCoordinates);
transformation = new ReferenceTransformation(sourcePolygon, parameters);
transformation.Execute();
Assert.IsInstanceOf <IPolygon>(transformation.Result);
Assert.AreEqual(sourcePolygon.Metadata, transformation.Result.Metadata);
IPolygon resultPolygon = transformation.Result as IPolygon;
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultPolygon.Shell.GetCoordinate(i));
}
// polygon with hole
sourcePolygon = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreatePolygon(sourceCoordinates, new Coordinate[][] { sourceCoordinates });
transformation = new ReferenceTransformation(sourcePolygon, parameters);
transformation.Execute();
Assert.IsInstanceOf <IPolygon>(transformation.Result);
Assert.AreEqual(sourcePolygon.Metadata, transformation.Result.Metadata);
resultPolygon = transformation.Result as IPolygon;
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultPolygon.Shell.GetCoordinate(i));
}
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultPolygon.Holes[0].GetCoordinate(i));
}
// multi point
IMultiPoint sourceMultiPoint = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreateMultiPoint(sourceCoordinates);
transformation = new ReferenceTransformation(sourceMultiPoint, parameters);
transformation.Execute();
Assert.IsInstanceOf <IMultiPoint>(transformation.Result);
Assert.AreEqual(sourceMultiPoint.Metadata, transformation.Result.Metadata);
IMultiPoint resultMultiPoint = transformation.Result as IMultiPoint;
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultMultiPoint[i].Coordinate);
}
// multi linestring
IMultiLineString sourceMultiLineString = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreateMultiLineString(Enumerable.Repeat(sourceLineString, 5));
transformation = new ReferenceTransformation(sourceMultiLineString, parameters);
transformation.Execute();
Assert.IsInstanceOf <IMultiLineString>(transformation.Result);
Assert.AreEqual(sourceMultiLineString.Metadata, transformation.Result.Metadata);
IMultiLineString resultMultiLineString = transformation.Result as IMultiLineString;
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultMultiLineString[0].GetCoordinate(i));
}
// multi polygon
IMultiPolygon sourceMultiPolygon = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreateMultiPolygon(Enumerable.Repeat(sourcePolygon, 5));
transformation = new ReferenceTransformation(sourceMultiPolygon, parameters);
transformation.Execute();
Assert.IsInstanceOf <IMultiPolygon>(transformation.Result);
Assert.AreEqual(sourceMultiPolygon.Metadata, transformation.Result.Metadata);
IMultiPolygon resultMultiPolygon = transformation.Result as IMultiPolygon;
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], resultMultiLineString[0].GetCoordinate(i));
}
// geometry collection
IGeometryCollection <IGeometry> sourceCollection = new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV).CreateGeometryCollection(new IGeometry[] { sourceLine, sourceLineString, sourceLinearRing, sourcePolygon });
transformation = new ReferenceTransformation(sourceCollection, parameters);
transformation.Execute();
Assert.IsInstanceOf <IGeometryCollection <IGeometry> >(transformation.Result);
Assert.AreEqual(sourceCollection.Metadata, transformation.Result.Metadata);
IGeometryCollection <IGeometry> resultCollection = transformation.Result as IGeometryCollection <IGeometry>;
for (Int32 i = 0; i < (resultCollection[0] as ILine).Count; i++)
{
Assert.AreEqual(expectedCoordinates[i], (resultCollection[0] as ILine).GetCoordinate(i));
}
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], (resultCollection[1] as ILineString).GetCoordinate(i));
}
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], (resultCollection[2] as ILinearRing).GetCoordinate(i));
}
for (Int32 i = 0; i < expectedCoordinates.Length; i++)
{
Assert.AreEqual(expectedCoordinates[i], (resultCollection[3] as IPolygon).Shell.GetCoordinate(i));
}
// not supported geometry
Mock <IGeometry> geometryMock = new Mock <IGeometry>(MockBehavior.Loose);
geometryMock.Setup(geometry => geometry.Factory).Returns(() => new GeometryFactory(ProjectedCoordinateReferenceSystems.HD72_EOV));
geometryMock.Setup(geometry => geometry.ReferenceSystem).Returns(() => ProjectedCoordinateReferenceSystems.HD72_EOV);
transformation = new ReferenceTransformation(geometryMock.Object, parameters);
Assert.Throws <InvalidOperationException>(() => transformation.Execute());
}
