/// <summary>
/// Creates a new MultiPoint geometry from a MultiPoint shape
/// </summary>
/// <param name="shape">The shape to convert</param>
/// <param name="factory">The geometry factory to use.</param>
/// <param name="copyAttributes">A value indicating whether or not to copy the <see cref="Data.Shape.Attributes"/> to <see cref="GeoAPIGeometry.UserData"/></param>
/// <returns>The geometry representing the converted shape.</returns>
/// <returns></returns>
private static GeoAPIGeometry FromMultiPointShape(Data.Shape shape, GeoAPIGeometryFactory factory, bool copyAttributes)
{
var coords = new List <GeoAPICoordinate>();
foreach (var part in shape.Range.Parts)
{
var i = part.StartIndex;
foreach (var vertex in part)
{
var c = new GeoAPICoordinate(vertex.X, vertex.Y);
coords.Add(c);
//if (shape.HasM()) c.M = shape.M[i];
if (shape.HasZ())
{
c.Z = shape.Z[i];
}
i++;
}
}
var ret = factory.CreateMultiPoint(coords.ToArray());
if (copyAttributes)
{
ret.UserData = shape.Attributes;
}
return(ret);
}
/// <summary>
/// Get the point for this shape if this is a point shape.
/// </summary>
/// <param name="shape">The shape to convert</param>
/// <param name="factory">The geometry factory to use.</param>
/// <param name="copyAttributes">A value indicating whether or not to copy the <see cref="Data.Shape.Attributes"/> to <see cref="GeoAPIGeometry.UserData"/></param>
/// <returns>The geometry representing the converted shape.</returns>
private static GeoAPIGeometry FromPointShape(Data.Shape shape, GeoAPIGeometryFactory factory, bool copyAttributes)
{
var part = shape.Range.Parts[0];
var i = part.StartIndex;
var c = new GeoAPICoordinate(part.Vertices[0], part.Vertices[1]);
//if (shape.HasM()) c.M = shape.M[i]
if (shape.HasZ())
{
c.Z = shape.Z[i];
}
var ret = factory.CreatePoint(c);
if (copyAttributes)
{
ret.UserData = shape.Attributes;
}
return(ret);
}
public static GeoAPIGeometry ToGeoAPI(this Data.Shape shape, GeoAPI.Geometries.IGeometryFactory factory = null, bool copyAttributes = false)
{
var useFactory = factory ?? GeoAPI.GeometryServiceProvider.Instance.CreateGeometryFactory();
if (shape.Range.FeatureType == FeatureType.Polygon)
{
return(FromPolygonShape(shape, useFactory, copyAttributes));
}
if (shape.Range.FeatureType == FeatureType.Line)
{
return(FromLineShape(shape, useFactory, copyAttributes));
}
if (shape.Range.FeatureType == FeatureType.MultiPoint)
{
return(FromMultiPointShape(shape, useFactory, copyAttributes));
}
if (shape.Range.FeatureType == FeatureType.Point)
{
return(FromPointShape(shape, useFactory, copyAttributes));
}
return(null);
}
/// <summary>
/// Gets the line for the specified index
/// </summary>
/// <param name="shape">The shape to convert</param>
/// <param name="factory">The geometry factory to use.</param>
/// <param name="copyAttributes">A value indicating whether or not to copy the <see cref="Data.Shape.Attributes"/> to <see cref="GeoAPIGeometry.UserData"/></param>
/// <returns>The geometry representing the converted shape.</returns>
private static GeoAPIGeometry FromLineShape(Data.Shape shape, GeoAPIGeometryFactory factory, bool copyAttributes)
{
var lines = new List <GeoAPILineString>();
foreach (var part in shape.Range.Parts)
{
var i = part.StartIndex;
var coords = new List <GeoAPICoordinate>();
foreach (var d in part)
{
var c = new GeoAPICoordinate(d.X, d.Y);
coords.Add(c);
//if (shape.HasM()) c.M = M[i];
if (shape.HasZ())
{
c.Z = shape.Z[i];
}
i++;
}
lines.Add(factory.CreateLineString(coords.ToArray()));
}
if (lines.Count == 1)
{
if (copyAttributes)
{
lines[0].UserData = shape.Attributes;
}
return(lines[0]);
}
var ret = factory.CreateMultiLineString(lines.ToArray());
if (copyAttributes)
{
ret.UserData = shape.Attributes;
}
return(ret);
}
private static bool HasZ(this Data.Shape self)
{
return(self.Z != null && self.Z.Length > 0);
}
/// <summary>
/// Creates a Polygon or MultiPolygon from this Polygon shape.
/// </summary>
/// <param name="shape">The shape to convert</param>
/// <param name="factory">The geometry factory to use.</param>
/// <param name="copyAttributes">A value indicating whether or not to copy the <see cref="Data.Shape.Attributes"/> to <see cref="GeoAPIGeometry.UserData"/></param>
/// <returns>The geometry representing the converted shape.</returns>
private static GeoAPIGeometry FromPolygonShape(Data.Shape shape, GeoAPIGeometryFactory factory, bool copyAttributes)
{
var shells = new List <GeoAPILinearRing>();
var holes = new List <GeoAPILinearRing>();
foreach (var part in shape.Range.Parts)
{
var coords = new List <GeoAPICoordinate>();
var i = part.StartIndex;
foreach (var d in part)
{
var c = new GeoAPICoordinate(d.X, d.Y);
//if (shape.HasM()) c.M = M[i];
if (shape.HasZ())
{
c.Z = shape.Z[i];
}
i++;
coords.Add(c);
}
var ring = factory.CreateLinearRing(coords.ToArray());
if (shape.Range.Parts.Count == 1)
{
shells.Add(ring);
}
else
{
if (NetTopologySuite.Algorithm.CGAlgorithms.IsCCW(ring.Coordinates))
{
holes.Add(ring);
}
else
{
shells.Add(ring);
}
}
}
//// Now we have a list of all shells and all holes
var holesForShells = new List <GeoAPILinearRing> [shells.Count];
for (var i = 0; i < shells.Count; i++)
{
holesForShells[i] = new List <GeoAPILinearRing>();
}
// Find holes
foreach (var testRing in holes)
{
GeoAPILinearRing minShell = null;
GeoAPIEnvelope minEnv = null;
var testEnv = testRing.EnvelopeInternal;
var testPt = testRing.Coordinates[0];
for (int j = 0; j < shells.Count; j++)
{
var tryRing = shells[j];
var tryEnv = tryRing.EnvelopeInternal;
if (minShell != null)
{
minEnv = minShell.EnvelopeInternal;
}
var isContained = tryEnv.Contains(testEnv) &&
(NetTopologySuite.Algorithm.CGAlgorithms.IsPointInRing(testPt, tryRing.Coordinates) ||
(PointInList(testPt, tryRing.Coordinates)));
// Check if this new containing ring is smaller than the current minimum ring
if (isContained)
{
if (minShell == null || minEnv.Contains(tryEnv))
{
minShell = tryRing;
}
holesForShells[j].Add(testRing);
}
}
}
var polygons = new GeoAPIPolygon[shells.Count];
for (var i = 0; i < shells.Count; i++)
{
polygons[i] = factory.CreatePolygon(shells[i], holesForShells[i].ToArray());
}
if (polygons.Length == 1)
{
if (copyAttributes)
{
polygons[0].UserData = shape.Attributes;
}
return(polygons[0]);
}
// It's a multi part
var ret = factory.CreateMultiPolygon(polygons);
if (copyAttributes)
{
ret.UserData = shape.Attributes;
}
return(ret);
}
