[NetTopologySuite]（1）线面相交

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

  1 Polygon inputGeometry = null;
  2         LineString analysisGeoemtry = null;
  3         private void button1_Click(object sender, EventArgs e)
  4         {
  5             var coords = new Coordinate[7];
  6             // Shell Coordinates
  7             var coordscheck = new Coordinate[7];
  8             double[] Xs = new double[] { 30, 260, 130, 350, 250, 330, 320, 180, 250, 50, 70, 60, 30, 260 };
  9             double[] Ys = new double[] { 140, 260, 110, 150, 230, 210, 140, 260 };
 10             for (var i = 0; i < Xs.Length / 2; i++)
 11             {
 12                 var x = Xs[i * 2];
 13                 var y = Xs[i * 2 + 1];
 14                 coords[i] = new Coordinate(x, y);
 15                 coordscheck[i] = new Coordinate(x, y);
 16             }
 17             //coordscheck[19] = new Coordinate(coords[0].X, coords[0].Y);
 18             //coords[19] = new Coordinate(coords[0].X, coords[0].Y);
 19             // Shell Rings
 20             var ring = new LinearRing(coords);
 21             var gf = new GeometryFactory();
 22             var ringCheck = gf.CreateLinearRing(coordscheck);
 23             // Hole Coordinates
 24             var coordsholecheck = new Coordinate[4];
 25             var coordshole = new Coordinate[4];
 26             for (var i = 0; i < Ys.Length / 2; i++)
 27             {
 28                 var x = Ys[i * 2];
 29                 var y = Ys[i * 2 + 1];
 30                 coordshole[i] = new Coordinate(x, y);
 31                 coordsholecheck[i] = new Coordinate(x, y);
 32             }
 33             // Hole LinearRing Arrays
 34             var hole = new LinearRing(coordshole);
 35             var holes = new ILinearRing[1];
 36             var holeCheck = gf.CreateLinearRing(coordsholecheck);
 37             var holescheck = new ILinearRing[1];
 38             holes[0] = hole;
 39             holescheck[0] = holeCheck;
 40             Polygon pg = new Polygon(ring, holes);//
 41             var polygonCheck = gf.CreatePolygon(ringCheck, holescheck);
 42             var areaCheck = polygonCheck.Area;
 43             var area = pg.Area;
 44             //
 45             inputGeometry = pg;
 46 
 47             double[] Lines = new double[] { 30, 50, 320, 360 };
 48             var coords1 = new Coordinate[2];
 49             var rnd1 = new Random();
 50             var coordscheck1 = new Coordinate[2];
 51             for (var i = 0; i < Lines.Length / 2; i++)
 52             {
 53                 var x = Lines[i * 2];
 54                 var y = Lines[i * 2 + 1];
 55                 coords1[i] = new Coordinate(x, y);
 56                 coordscheck1[i] = new Coordinate(x, y);
 57             }
 58 
 59             var gf1 = new GeometryFactory();
 60             var lscheck = gf1.CreateLineString(coordscheck1);
 61             LineString ls = new LineString(coords1);
 62             analysisGeoemtry = ls;
 63             intersectionOutput = analysisGeoemtry.Intersection(inputGeometry);
 64             this.Invalidate();
 65         }
 66         IGeometry intersectionOutput = null;
 67         private void Form1_Paint(object sender, PaintEventArgs e)
 68         {
 69             Graphics gc = this.CreateGraphics();
 70             if (inputGeometry != null)
 71             {
 72                 GraphicsHelper.DrawPolygon(gc, inputGeometry, new SolidBrush(Color.Red), new Pen(new SolidBrush(Color.LightYellow)), false);
 73             }
 74             if (analysisGeoemtry != null)
 75             {
 76                 GraphicsHelper.DrawLineString(gc, analysisGeoemtry, new Pen(new SolidBrush(Color.Blue)));
 77             }
 78             if (intersectionOutput != null)
 79             {
 80                 if (intersectionOutput is LineString)
 81                 {
 82                     LineString lsr = intersectionOutput as LineString;
 83                     Pen mp = new Pen(new SolidBrush(Color.Yellow));
 84                     mp.Width = 2;
 85                     GraphicsHelper.DrawLineString(gc, lsr, mp);
 86                 }
 87                 if (intersectionOutput is MultiLineString)
 88                 {
 89                     MultiLineString mlsr = intersectionOutput as MultiLineString;
 90                     for (int i = 0; i < mlsr.NumGeometries; i++)
 91                     {
 92                         if (mlsr.Geometries[i] is LineString)
 93                         {
 94                             LineString lsr = mlsr.Geometries[i] as LineString;
 95                             Pen mp = new Pen(new SolidBrush(Color.Yellow));
 96                             mp.Width = 2;
 97                             GraphicsHelper.DrawLineString(gc, lsr, mp);
 98                         }
 99                     }
100                 }
101             }
102             gc.Dispose();
103         }