[NetTopologySuite]（1）线面相交

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

 Polygon inputGeometry = null;
         LineString analysisGeoemtry = null;
         private void button1_Click(object sender, EventArgs e)
         {
             var coords = new Coordinate[];
             // Shell Coordinates
             var coordscheck = new Coordinate[];
             double[] Xs = new double[] { , , , , , , , , , , , , ,  };
             double[] Ys = new double[] { , , , , , , ,  };
             for (var i = ; i < Xs.Length / ; i++)
             {
                 var x = Xs[i * ];
                 var y = Xs[i *  + ];
                 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[];
             var coordshole = new Coordinate[];
             for (var i = ; i < Ys.Length / ; i++)
             {
                 var x = Ys[i * ];
                 var y = Ys[i *  + ];
                 coordshole[i] = new Coordinate(x, y);
                 coordsholecheck[i] = new Coordinate(x, y);
             }
             // Hole LinearRing Arrays
             var hole = new LinearRing(coordshole);
             var holes = new ILinearRing[];
             var holeCheck = gf.CreateLinearRing(coordsholecheck);
             var holescheck = new ILinearRing[];
             holes[] = hole;
             holescheck[] = 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[] { , , ,  };
             var coords1 = new Coordinate[];
             var rnd1 = new Random();
             var coordscheck1 = new Coordinate[];
             for (var i = ; i < Lines.Length / ; i++)
             {
                 var x = Lines[i * ];
                 var y = Lines[i *  + ];
                 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 = ;
                     GraphicsHelper.DrawLineString(gc, lsr, mp);
                 }
                 if (intersectionOutput is MultiLineString)
                 {
                     MultiLineString mlsr = intersectionOutput as MultiLineString;
                     for (int i = ; 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 = ;
                             GraphicsHelper.DrawLineString(gc, lsr, mp);
                         }
                     }
                 }
             }
             gc.Dispose();
         }