dotspatial mysql_Dotspatial 要素重叠分析

private void toolStripButton30_Click(object sender, EventArgs e)

{

//面状重叠分析

if (mapMain.Layers.Count == 0)

{

return;

}

//重叠分析

//遍历要素，显示面积

FeatureSet fs = null;

fs = Lzq_LayerManager.getFeatureSetByName(layerNamePolygon, mapMain); //(FeatureSet) map1.Layers[0].DataSet;

//MessageBox.Show("图层类型：" + fs.FeatureType.ToString());

if (fs.FeatureType != FeatureType.Polygon)

{

return;

}

DotSpatial.Topology.Geometry gm1, gm2;

for (int i = 0; i < fs.Features.Count; i++)

{

gm1 = (DotSpatial.Topology.Geometry) (fs.Features[i].BasicGeometry);

for (int j = i + 1; j < fs.Features.Count; j++)

{

gm2 = (DotSpatial.Topology.Geometry) (fs.Features[j].BasicGeometry);

if (gm1.IsSimple && gm2.IsSimple)

{

try

{

if (gm1.Overlaps(gm2)) //重叠

{

//try

//{

Geometry p = (Geometry)gm1.Intersection(gm2);

DotSpatial.Topology.CoordinateArrays.RemoveRepeatedPoints(p.Coordinates);

if (p.Area > 0.001)

{

MessageBox.Show("\ngm1——lb:" +

fs.Features[i].DataRow["林班号"].ToString() + " xb:" +

fs.Features[i].DataRow["小班号"].ToString() + "——gm1 Area" +

gm1.Area.ToString() +

"\ngm2——lb:" +

fs.Features[j].DataRow["林班号"].ToString() + " xb:" +

fs.Features[j].DataRow["小班号"].ToString() + "——gm2 Area" +

gm2.Area.ToString() +

"\n重叠 面积:" +

p.Area);

}

//}

//catch

//{

//}

}

}

catch

{

}

}

//}

//catch (Exception exc)

//{

// // MessageBox.Show(exc.Message);

//}

}

}

fs = null;

}