using System;

using System.Collections.Generic;

using NetTopologySuite.Operation.Distance;

using System.Linq;

public class Program

{

    public static void Main()

    {

        var gf32632 = NetTopologySuite.NtsGeometryServices.Instance.CreateGeometryFactory(32632);

        var journeyPatternCoordinates = new List<NetTopologySuite.Geometries.Coordinate>()

            {

                new(100, 100),

                new(110, 105),

                new(140, 95),

                new(230, 160),

                new(180, 120),

                new(160, 50),

                new(150, -20),

                new(140, -70),

                new(10, -120),

                new(-80, -230),

                new(-170, -290)

            };

        var serviceLinkCoords = new List<NetTopologySuite.Geometries.Coordinate>();

        var pntStart = gf32632.CreatePoint(new NetTopologySuite.Geometries.Coordinate(133, 88));

        var pntEnd = gf32632.CreatePoint(new NetTopologySuite.Geometries.Coordinate(155, 41));

​

        double minDistanceStart = double.MaxValue;

        int breakIndexStart = 0;

        // scorro tutte le coordinate fino a trovare la coppia con distanza minima rispetto alla fermata di PARTENZA della tratta

        // Ris: C1--C2---C3-C4 break 4

        for (int k = 1; k < journeyPatternCoordinates.Count; k++)

        {

            var prev = journeyPatternCoordinates[k- 1];

            var curr = journeyPatternCoordinates[k];

​

            var pairLine = gf32632.CreateLineString(

                new NetTopologySuite.Geometries.Coordinate[]

                {

                    prev,

                    curr

                });

            var distOp = new DistanceOp(pairLine, pntStart);

            var distance = distOp.Distance();

            var closestPt = distOp.NearestPoints();

​

            if (distance < minDistanceStart)

            {

                minDistanceStart = distance;

                breakIndexStart = k;

            }

        }

        Console.WriteLine("breakIndexStart = " + breakIndexStart);

​

​

        double minDistanceEnd = double.MaxValue;

        int breakIndexEnd = 0;

        // scorro tutte le coordinate fino a trovare la coppia con distanza minima rispetto alla fermata di ARRIVO della tratta

        // Ris C4----C5--C6-----C7 break 7

        for (int k = breakIndexEnd + 1; k < journeyPatternCoordinates.Count; k++)

        {

            var prev = journeyPatternCoordinates[k - 1];

            var curr = journeyPatternCoordinates[k];

​

            var pairLine = gf32632.CreateLineString(

                new NetTopologySuite.Geometries.Coordinate[]

                {

                    prev,

                    curr

                });

            var distOp = new DistanceOp(pairLine, pntEnd);

            var distance = distOp.Distance();

            var closestPt = distOp.NearestPoints();

​

            if (distance < minDistanceEnd)

            {

                minDistanceEnd = distance;

                breakIndexEnd = k;

            }

        }

        Console.WriteLine("breakIndexEnd = " + breakIndexEnd);

​

        // Ris PC--C5--C6---PD

        serviceLinkCoords = new List<NetTopologySuite.Geometries.Coordinate>() { pntStart.Coordinate };

        serviceLinkCoords.AddRange(journeyPatternCoordinates.Skip(breakIndexStart).Take(breakIndexEnd - breakIndexStart));

        serviceLinkCoords.Add(pntEnd.Coordinate);

        foreach(var s in serviceLinkCoords)

            Console.WriteLine(s);

        //

    }

}