[NetTopologySuite]（1）线面相交

用DotSpatial.Topology进行的测试，即使用NetTopologySuite类库进行测试：

Polygon inputGeometry = null;
        LineString analysisGeoemtry = null;
        private void button1_Click(object sender, EventArgs e)
        {
            var coords = new Coordinate[7];
            // Shell Coordinates
            var coordscheck = new Coordinate[7];
            double[] Xs = new double[] { 30, 260, 130, 350, 250, 330, 320, 180, 250, 50, 70, 60, 30, 260 };
            double[] Ys = new double[] { 140, 260, 110, 150, 230, 210, 140, 260 };
            for (var i = 0; i < Xs.Length / 2; i++)
            {
                var x = Xs[i * 2];
                var y = Xs[i * 2 + 1];
                coords[i] = new Coordinate(x, y);
                coordscheck[i] = new Coordinate(x, y);
            }
            //coordscheck[19] = new Coordinate(coords[0].X, coords[0].Y);
            //coords[19] = new Coordinate(coords[0].X, coords[0].Y);
            // Shell Rings
            var ring = new LinearRing(coords);
            var gf = new GeometryFactory();
            var ringCheck = gf.CreateLinearRing(coordscheck);
            // Hole Coordinates
            var coordsholecheck = new Coordinate[4];
            var coordshole = new Coordinate[4];
            for (var i = 0; i < Ys.Length / 2; i++)
            {
                var x = Ys[i * 2];
                var y = Ys[i * 2 + 1];
                coordshole[i] = new Coordinate(x, y);
                coordsholecheck[i] = new Coordinate(x, y);
            }
            // Hole LinearRing Arrays
            var hole = new LinearRing(coordshole);
            var holes = new ILinearRing[1];
            var holeCheck = gf.CreateLinearRing(coordsholecheck);
            var holescheck = new ILinearRing[1];
            holes[0] = hole;
            holescheck[0] = holeCheck;
            Polygon pg = new Polygon(ring, holes);//
            var polygonCheck = gf.CreatePolygon(ringCheck, holescheck);
            var areaCheck = polygonCheck.Area;
            var area = pg.Area;
            //
            inputGeometry = pg;

            double[] Lines = new double[] { 30, 50, 320, 360 };
            var coords1 = new Coordinate[2];
            var rnd1 = new Random();
            var coordscheck1 = new Coordinate[2];
            for (var i = 0; i < Lines.Length / 2; i++)
            {
                var x = Lines[i * 2];
                var y = Lines[i * 2 + 1];
                coords1[i] = new Coordinate(x, y);
                coordscheck1[i] = new Coordinate(x, y);
            }

            var gf1 = new GeometryFactory();
            var lscheck = gf1.CreateLineString(coordscheck1);
            LineString ls = new LineString(coords1);
            analysisGeoemtry = ls;
            intersectionOutput = analysisGeoemtry.Intersection(inputGeometry);
            this.Invalidate();
        }
        IGeometry intersectionOutput = null;
        private void Form1_Paint(object sender, PaintEventArgs e)
        {
            Graphics gc = this.CreateGraphics();
            if (inputGeometry != null)
            {
                GraphicsHelper.DrawPolygon(gc, inputGeometry, new SolidBrush(Color.Red), new Pen(new SolidBrush(Color.LightYellow)), false);
            }
            if (analysisGeoemtry != null)
            {
                GraphicsHelper.DrawLineString(gc, analysisGeoemtry, new Pen(new SolidBrush(Color.Blue)));
            }
            if (intersectionOutput != null)
            {
                if (intersectionOutput is LineString)
                {
                    LineString lsr = intersectionOutput as LineString;
                    Pen mp = new Pen(new SolidBrush(Color.Yellow));
                    mp.Width = 2;
                    GraphicsHelper.DrawLineString(gc, lsr, mp);
                }
                if (intersectionOutput is MultiLineString)
                {
                    MultiLineString mlsr = intersectionOutput as MultiLineString;
                    for (int i = 0; i < mlsr.NumGeometries; i++)
                    {
                        if (mlsr.Geometries[i] is LineString)
                        {
                            LineString lsr = mlsr.Geometries[i] as LineString;
                            Pen mp = new Pen(new SolidBrush(Color.Yellow));
                            mp.Width = 2;
                            GraphicsHelper.DrawLineString(gc, lsr, mp);
                        }
                    }
                }
            }
            gc.Dispose();
        }