public static IShapeCollection AdditiveSimplifyByAngle(IShapeCollection shapes, double threshold)
{
return(Process(shapes, threshold, (i, t) => ShapeUtility.AdditiveSimplifyByAngle(i, t)));
//if (shapes.First().Type == ShapeType.Polygon || shapes.First().Type == ShapeType.PolygonZ)
//{
// var result = new ShapeCollection<Polygon>();
// foreach (ISimplePoints feature in shapes)
// {
// //var temp = ShapeUtility.SimplifyByAngle(feature.Points, threshold);
// var temp = feature.Parts.Select((i, index) => ShapeUtility.AdditiveSimplifyByAngle(feature.GetPart(index), threshold)).ToArray();
// if (temp == null || temp[0].Length == 0)
// continue;
// //if (temp.Length > 2)
// //{
// result.Add(new Polygon(temp));
// //}
// }
// return result;
//}
//else if (shapes.First().Type == ShapeType.PolyLine || shapes.First().Type == ShapeType.PolyLineZ || shapes.First().Type == ShapeType.PolyLineM)
//{
// var result = new ShapeCollection<PolyLine>();
// foreach (ISimplePoints feature in shapes)
// {
// //var temp = ShapeUtility.SimplifyByAngle(feature.Points, threshold);
// var temp = feature.Parts.Select((i, index) => ShapeUtility.AdditiveSimplifyByAngle(feature.GetPart(index), threshold)).ToArray();
// if (temp == null || temp.Length == 0 || temp.Sum(i => i.Length) == 0)
// continue;
// //if (temp.Length > 1)
// //{
// result.Add(new PolyLine(temp));
// //}
// }
// return result;
//}
//else
//{
// throw new NotImplementedException();
//}
}
public static IShapeCollection SimplifyByAngle(IShapeCollection shapes, double threshold)
{
return(Process(shapes, threshold, (i, t) => ShapeUtility.SimplifyByAngle(i, t)));
//if (shapes.First().Type == ShapeType.Polygon || shapes.First().Type == ShapeType.PolygonZ)
//{
// var result = new ShapeCollection<Polygon>();
// foreach (ISimplePoints feature in shapes)
// {
// var temp = ShapeUtility.SimplifyByAngle(feature.Points, threshold);
// if (temp == null)
// continue;
// if (temp.Length > 2)
// {
// result.Add(new Polygon(temp));
// }
// }
// return result;
//}
//else if (shapes.First().Type == ShapeType.PolyLine || shapes.First().Type == ShapeType.PolyLineZ || shapes.First().Type == ShapeType.PolyLineM)
//{
// var result = new ShapeCollection<PolyLine>();
// foreach (ISimplePoints feature in shapes)
// {
// var temp = ShapeUtility.SimplifyByAngle(feature.Points, threshold);
// if (temp == null)
// continue;
// if (temp.Length > 2)
// {
// result.Add(new PolyLine(temp));
// }
// }
// return result;
//}
//else
//{
// throw new NotImplementedException();
//}
}
public static IShapeCollection AdditiveSimplifyByArea(IShapeCollection shapes, double threshold)
{
return(Process(shapes, threshold, (i, t) => ShapeUtility.AdditiveSimplifyByArea(i, t)));
//ShapeCollection<Polygon> result = new ShapeCollection<Polygon>();
//foreach (ISimplePoints feature in shapes)
//{
// var points = feature.Parts.Select((i, index) => ShapeUtility.AdditiveSimplifyByArea(feature.GetPart(index), threshold)).ToArray();
// if (points == null)
// continue;
// result.Add(new Polygon(points));
//}
//return result;
}
public static IShapeCollection SimplifyByArea(IShapeCollection shapes, double threshold)
{
return(Process(shapes, threshold, (i, t) => ShapeUtility.SimplifyByArea(i, t)));
//ShapeCollection<PolyLine> result = new ShapeCollection<PolyLine>();
//foreach (ISimplePoints feature in shapes)
//{
// var temp = ShapeUtility.SimplifyByArea(feature.Points, threshold);
// if (temp == null)
// continue;
// if (temp.Length > 2)
// {
// result.Add(new PolyLine(temp));
// }
//}
//return result;
}
public static IShapeCollection AdditiveSimplifyByAreaPlus(IShapeCollection shapes, double threshold)
{
return(Process(shapes, threshold, (i, t) => ShapeUtility.AdditiveSimplifyByAreaPlus(i, t)));
}
static void SuperSlim(string file)
{
var shapes = IRI.Ket.ShapefileFormat.Shapefile.Read(file);
ShapeCollection <Polygon> polygons = new ShapeCollection <Polygon>();
//var stat = new FeatureStatistics[120];
int index = 0;
foreach (Polygon feature in shapes)
{
index = 0;
//UpdateFeatureStat(ref stat[index++], feature.Points, -1);
var temp = ShapeUtility.SimplifyByArea(feature.Points, 100);
//UpdateFeatureStat(ref stat[index++], temp, 100);
if (temp == null)
{
continue;
}
for (int i = 1; i < 30; i++)
{
var threshold = i * 100 + 100;
temp = ShapeUtility.SimplifyByArea(temp, threshold);
//UpdateFeatureStat(ref stat[index++], temp, threshold);
}
for (int i = 1; i < 11; i++)
{
var threshold = (i * 500) + 5000;
temp = ShapeUtility.SimplifyByArea(temp, threshold);
//UpdateFeatureStat(ref stat[index++], temp, threshold);
}
for (int i = 1; i < 11; i++)
{
var threshold = (i * 1000) + 10000;
temp = ShapeUtility.SimplifyByArea(temp, threshold);
//UpdateFeatureStat(ref stat[index++], temp, threshold);
}
for (int i = 1; i < 17; i++)
{
var threshold = (i * 5000) + 20000;
temp = ShapeUtility.SimplifyByArea(temp, threshold);
//UpdateFeatureStat(ref stat[index++], temp, threshold);
}
for (int i = 1; i < 25; i++)
{
var threshold = (i * 20000) + 100000;
temp = ShapeUtility.SimplifyByArea(temp, threshold);
//UpdateFeatureStat(ref stat[index++], temp, threshold);
}
if (temp.Length > 0)
{
polygons.Add(new Polygon(temp));
}
}
IRI.Ket.ShapefileFormat.Writer.ShpWriter.Write(@"E:\Data\0. Test\Large Data\94.04.10\IRI1OstanMercatorSuperSlim.shp", polygons, true);
//Debug.Print("\n****************\n");
//foreach (var item in stat)
//{
// Debug.Print(item.ToString());
//}
//Debug.Print("\n****************\n");
//foreach (var item in stat)
//{
// Debug.Print(item.ToCSV());
//}
}