public void TestPolgonWithCollapsedPointPointwise()
{
var g = reader.Read("POLYGON ((10 10, 100 100, 200 10.1, 300 100, 400 10, 10 10))");
var g2 = reader.Read("POLYGON ((10 10, 100 100, 200 10, 300 100, 400 10, 10 10))");
var gReduce = GeometryPrecisionReducer.ReducePointwise(g, pmFixed1);
AssertEqualsExactAndHasSameFactory(gReduce, g2);
}
//=======================================
private void CheckReducePointwise(string wkt, string wktExpected)
{
var g = Read(wkt);
var gExpected = Read(wktExpected);
var gReduce = GeometryPrecisionReducer.ReducePointwise(g, _pmFixed1);
AssertEqualsExactAndHasSameFactory(gExpected, gReduce);
}
/// <summary>
/// Reduces precision pointwise, then snap-rounds.
/// Note that output set may not contain non-unique linework
/// (and thus cannot be used as input to Polygonizer directly).
/// UnaryUnion is one way to make the linework unique.
/// </summary>
/// <param name="geom">A geometry containing linework to node</param>
/// <param name="scaleFactor">the precision model scale factor to use</param>
/// <returns>The noded, snap-rounded linework</returns>
public static Geometry SnapRoundWithPointwisePrecisionReduction(Geometry geom, double scaleFactor)
{
var pm = new PrecisionModel(scaleFactor);
var roundedGeom = GeometryPrecisionReducer.ReducePointwise(geom, pm);
var geomList = new List <Geometry>();
geomList.Add(roundedGeom);
var noder = new GeometryNoder(pm);
var lines = noder.Node(geomList);
return(FunctionsUtil.GetFactoryOrDefault(geom).BuildGeometry(lines.Cast <Geometry>().ToArray()));
}
/// <summary>Reduces precision pointwise, then snap-rounds.
/// <para/>
/// Note that output set may not contain non-unique linework
/// (and thus cannot be used as input to Polygonizer directly).
/// <c>UnaryUnion</c> is one way to make the linework unique.
/// </summary>
/// <param name="geom">A Geometry containing linework to node</param>
/// <param name="scaleFactor">The precision model scale factor to use</param>
/// <returns>The noded, snap-rounded linework</returns>
public static IGeometry SnapRoundLines(
IGeometry geom, double scaleFactor)
{
IPrecisionModel pm = new PrecisionModel(scaleFactor);
var roundedGeom = GeometryPrecisionReducer.ReducePointwise(geom, pm);
var geomList = new List <IGeometry>();
geomList.Add(roundedGeom);
var noder = new GeometrySnapRounder(pm);
var lines = noder.Node(geomList);
return(FunctionsUtil.GetFactoryOrDefault(geom).BuildGeometry(lines));
}
public static IGeometry SnapRound(
IGeometry geomA, IGeometry geomB,
double scaleFactor)
{
var pm = new PrecisionModel(scaleFactor);
var roundedGeomA = GeometryPrecisionReducer.ReducePointwise(geomA, pm);
var geomRound = roundedGeomA;
if (geomB != null)
{
var roundedGeomB = GeometryPrecisionReducer.ReducePointwise(geomB, pm);
geomRound = geomA.Factory.CreateGeometryCollection(new [] { roundedGeomA, roundedGeomB });
}
var noder = new GeometrySnapRounder(pm);
IGeometry snapped = noder.Node(geomRound);
return(snapped);
}
