private IModelPoint Polygon2ModelPoint(IPolygon polygon, out Gvitech.CityMaker.Resource.IModel gvModel)
{
if (polygon == null || !polygon.IsValid)
{
gvModel = null;
return(null);
}
IVector3 vector = polygon.QueryNormal();
if (vector.Z < 0.0)
{
polygon.ExteriorRing.ReverseOrientation();
}
IModelPoint result = null;
IGeometryConvertor geometryConvertor = new GeometryConvertorClass();
ITriMesh triMesh = geometryConvertor.PolygonToTriMesh(polygon);
if (triMesh != null)
{
IMultiTriMesh multiTriMesh = ((IGeometryFactory) new GeometryFactoryClass()).CreateGeometry(gviGeometryType.gviGeometryMultiTrimesh, polygon.VertexAttribute) as IMultiTriMesh;
multiTriMesh.AddTriMesh(triMesh);
bool flag = geometryConvertor.TriMeshToModelPoint(multiTriMesh, out gvModel, out result);
if (flag && gvModel != null)
{
this.Model2Water(ref gvModel);
}
}
else
{
result = null;
gvModel = null;
}
return(result);
}
private void btnConstructMultiTriMesh_Click(object sender, EventArgs e)
{
try
{
if (gfactory == null)
gfactory = new GeometryFactory();
multiTriMesh = gfactory.CreateGeometry(gviGeometryType.gviGeometryMultiTrimesh, gviVertexAttribute.gviVertexAttributeZ) as IMultiTriMesh;
if (multiTriMesh == null)
return;
triMesh = gfactory.CreateGeometry(gviGeometryType.gviGeometryTriMesh, gviVertexAttribute.gviVertexAttributeZ) as ITriMesh;
if (triMesh == null)
return;
ArrayList vertices = new ArrayList(); //ITopoNode数组
point = gfactory.CreatePoint(gviVertexAttribute.gviVertexAttributeZ);
for (int i = 0; i < 17; ++i)
{
point.SetCoords(i, i, i, 0, 0);
vertices.Add(triMesh.AddPoint(point));
}
triMesh.AddTriangle(vertices[0] as ITopoNode, vertices[1] as ITopoNode, vertices[2] as ITopoNode);
triMesh.AddTriangle(vertices[2] as ITopoNode, vertices[1] as ITopoNode, vertices[3] as ITopoNode);
triMesh.AddTriangle(vertices[3] as ITopoNode, vertices[1] as ITopoNode, vertices[4] as ITopoNode);
triMesh.AddTriangle(vertices[3] as ITopoNode, vertices[4] as ITopoNode, vertices[5] as ITopoNode);
triMesh.AddTriangle(vertices[5] as ITopoNode, vertices[4] as ITopoNode, vertices[6] as ITopoNode);
triMesh.AddTriangle(vertices[5] as ITopoNode, vertices[6] as ITopoNode, vertices[7] as ITopoNode);
triMesh.AddTriangle(vertices[7] as ITopoNode, vertices[6] as ITopoNode, vertices[8] as ITopoNode);
triMesh.AddTriangle(vertices[7] as ITopoNode, vertices[8] as ITopoNode, vertices[9] as ITopoNode);
triMesh.AddTriangle(vertices[7] as ITopoNode, vertices[9] as ITopoNode, vertices[10] as ITopoNode);
triMesh.AddTriangle(vertices[10] as ITopoNode, vertices[9] as ITopoNode, vertices[11] as ITopoNode);
triMesh.AddTriangle(vertices[10] as ITopoNode, vertices[11] as ITopoNode, vertices[12] as ITopoNode);
triMesh.AddTriangle(vertices[13] as ITopoNode, vertices[10] as ITopoNode, vertices[12] as ITopoNode);
triMesh.AddTriangle(vertices[13] as ITopoNode, vertices[14] as ITopoNode, vertices[10] as ITopoNode);
triMesh.AddTriangle(vertices[13] as ITopoNode, vertices[15] as ITopoNode, vertices[14] as ITopoNode);
triMesh.AddTriangle(vertices[15] as ITopoNode, vertices[16] as ITopoNode, vertices[14] as ITopoNode);
triMesh.AddTriangle(vertices[16] as ITopoNode, vertices[5] as ITopoNode, vertices[14] as ITopoNode);
triMesh.AddTriangle(vertices[3] as ITopoNode, vertices[5] as ITopoNode, vertices[16] as ITopoNode);
triMesh.AddTriangle(vertices[2] as ITopoNode, vertices[3] as ITopoNode, vertices[16] as ITopoNode);
multiTriMesh.AddGeometry(triMesh);
MultiSurface geo = new MultiSurface();
// Geometry属性
geo.Dimension = multiTriMesh.Dimension;
if (multiTriMesh.Envelope != null)
{
geo.MaxX = multiTriMesh.Envelope.MaxX;
geo.MaxY = multiTriMesh.Envelope.MaxY;
geo.MaxZ = multiTriMesh.Envelope.MaxZ;
geo.MinX = multiTriMesh.Envelope.MinX;
geo.MinY = multiTriMesh.Envelope.MinY;
geo.MinZ = multiTriMesh.Envelope.MinZ;
}
geo.GeometryType = multiTriMesh.GeometryType;
geo.IsEmpty = multiTriMesh.IsEmpty;
geo.IsValid = multiTriMesh.IsValid;
geo.VertexAttribute = multiTriMesh.VertexAttribute;
geo.HasId = multiTriMesh.HasId();
geo.HasM = multiTriMesh.HasM();
geo.HasZ = multiTriMesh.HasZ();
// GeometryCollection属性
geo.GeometryCount = multiTriMesh.GeometryCount;
geo.IsOverlap = multiTriMesh.IsOverlap;
// MultiSurface属性
geo.Area = multiTriMesh.GetArea();
this.propertyGrid1.SelectedObject = geo;
}
catch (System.Exception ex)
{
System.Diagnostics.Trace.WriteLine(ex.Message);
}
}
void axRenderControl1_RcMouseClickSelect(IPickResult PickResult, IPoint IntersectPoint, gviModKeyMask Mask, gviMouseSelectMode EventSender)
{
try
{
this.axRenderControl1.FeatureManager.UnhighlightAll();
if (PickResult != null)
{
if (PickResult.Type == gviObjectType.gviObjectFeatureLayer)
{
IFeatureLayerPickResult flpr = PickResult as IFeatureLayerPickResult;
int fid = flpr.FeatureId;
this.axRenderControl1.FeatureManager.HighlightFeature(_featureClass, fid, System.Drawing.Color.Yellow);
//////////////////////////////////////////////////////////////////////////
//
// GeometryConvert的代码添加在这里
//
//////////////////////////////////////////////////////////////////////////
fidList.Clear();
fidList.Add(fid);
IRowBuffer rowGC = _featureClass.GetRow(fidList[0]);
int nPose = rowGC.FieldIndex("Geometry");
if (nPose == -1)
{
MessageBox.Show("不存在Geometry列");
return;
}
// 获取ModelPoint
IModelPoint modelPointGC = null;
if (rowGC != null)
{
nPose = rowGC.FieldIndex("Geometry");
IGeometry geo = rowGC.GetValue(nPose) as IGeometry;
if (geo.GeometryType == gviGeometryType.gviGeometryModelPoint)
{
modelPointGC = geo as IModelPoint;
}
}
this.Text = "拾取成功";
// 获取Model
string modelName = modelPointGC.ModelName;
IModel modelGC = (_featureClass.FeatureDataSet as IResourceManager).GetModel(modelName);
// 1、获取MultiTriMesh
IGeometryConvertor gc = new GeometryConvertor();
IMultiTriMesh multiTM = gc.ModelPointToTriMesh(modelGC, modelPointGC, false);
this.Text = "ModelToTriMesh完成";
if (this.cbCreateRenderTriMesh.Checked)
{
// 创建RenderTriMesh在三维上显示
for (int i = 0; i < multiTM.GeometryCount; i++)
{
ITriMesh tm = multiTM.GetGeometry(i) as ITriMesh;
IRenderTriMesh rtm = this.axRenderControl1.ObjectManager.CreateRenderTriMesh(tm, null, rootId);
//随机颜色填充TriMesh
Random randObj = new Random(i);
int aColor = randObj.Next(0, 255);
int gColor = randObj.Next(0, 255);
int rColor = randObj.Next(0, 255);
ISurfaceSymbol ss = new SurfaceSymbol();
ss.Color = System.Drawing.Color.FromArgb(rColor, gColor, aColor);
rtm.Symbol = ss;
rTMeshList.Add(rtm);
}
}
if (this.cbCreateRenderPolygon.Checked)
{
// 2、获取投影MultiPolygon
IMultiPolygon multiPolygon = gc.ProjectTriMeshToPolygon(multiTM, 1.0);
this.Text = "MultiTriMeshToFootprint完成。面积:" + multiPolygon.GetArea();
// 创建RenderPolygon在三维上显示
for (int i = 0; i < multiPolygon.GeometryCount; i++)
{
IPolygon polygon = multiPolygon.GetGeometry(i) as IPolygon;
IRenderPolygon rpolygon = this.axRenderControl1.ObjectManager.CreateRenderPolygon(polygon, null, rootId);
rPolygonList.Add(rpolygon);
}
}
if (this.cbCreateRenderPolygon2.Checked)
{
// 3、获取切割MultiPolygon
double heightSpec = 0.0, heightIntersect = 0.0;
Double.TryParse(IntersectPoint.Z.ToString(), out heightIntersect);
fidList.Clear();
fidList.Add(fid);
IEnvelope box = _featureClass.GetFeaturesEnvelope(fidList.ToArray(), "Geometry");
heightSpec = System.Math.Abs(heightIntersect - box.MinZ);
// 注意：CutTriMeshToPolygon方法的第三个参数为空间分辨率，应该选择合适值。
// 值过大会导致结果不精确，值过小会导致转换时间过长。使用时应设置大小合适的值来平衡精度和效率问题。
IMultiPolygon multiPolygon2 = gc.CutTriMeshToPolygon(multiTM, heightSpec, 0.5);
this.Text = "MultiTriMeshToFootprint2完成。面积:" + multiPolygon2.GetArea();
// 创建RenderPolygon在三维上显示
for (int i = 0; i < multiPolygon2.GeometryCount; i++)
{
IPolygon polygon2 = multiPolygon2.GetGeometry(i) as IPolygon;
IRenderPolygon rpolygon2 = this.axRenderControl1.ObjectManager.CreateRenderPolygon(polygon2, null, rootId);
rPolygonList2.Add(rpolygon2);
}
}
if (this.cbCreateRenderMulPoint.Checked)
{
// 4、获取MultiPoint
IMultiPoint multiPoint = gc.MultiTriMeshToMultiPoint(multiTM, 3.0);
this.Text = "MultiTriMeshToMultiPoint完成。MultiPoint个数为：" + multiPoint.GeometryCount;
//创建RenderPoint在三维上显示
IRenderMultiPoint rpoint = this.axRenderControl1.ObjectManager.CreateRenderMultiPoint(multiPoint, null, rootId);
rPointList.Add(rpoint);
}
}
}
}
catch (System.Exception ex)
{
if (ex.GetType().Name.Equals("UnauthorizedAccessException"))
{
MessageBox.Show("需要标准runtime授权");
}
else
{
MessageBox.Show(ex.Message);
}
}
}
void axRenderControl1_RcObjectEditFinish()
{
this.axRenderControl1.InteractMode = gviInteractMode.gviInteractNormal;
polygonDraw = currentGeometry as IPolygon;
//抬高一点
for (int i = 0; i < polygonDraw.ExteriorRing.PointCount; i++)
{
IPoint pointOnExr = polygonDraw.ExteriorRing.GetPoint(i).Clone() as IPoint;
pointOnExr.Z += 1;
polygonDraw.ExteriorRing.UpdatePoint(i, pointOnExr);
}
renderpolygonDraw.SetFdeGeometry(polygonDraw);
renderpolygonDraw.VisibleMask = gviViewportMask.gviViewNone;
double height = 0.0;
try
{
height = double.Parse(this.numHeight.Value.ToString());
}
catch (System.Exception ex)
{
MessageBox.Show("请检查挖洞高度");
return;
}
//构造底面polygon
IPolygon polygonBottom = geoFactory.CreateGeometry(gviGeometryType.gviGeometryPolygon, gviVertexAttribute.gviVertexAttributeZ) as IPolygon;
for (int i = 0; i < polygonDraw.ExteriorRing.PointCount; i++)
{
IPoint pointOnExr = polygonDraw.ExteriorRing.GetPoint(i).Clone() as IPoint;
maxZ = pointOnExr.Z;
pointOnExr.Z -= height;
minZ = pointOnExr.Z;
polygonBottom.ExteriorRing.AppendPoint(pointOnExr);
}
polygonBottom.Close();
PolygonList.Add(polygonBottom);
//CreateRenderPolygon(polygonBottom);
//构造侧面polygon
for (int i = 0; i < polygonDraw.ExteriorRing.PointCount; i++)
{
IPoint pointOnExr = polygonDraw.ExteriorRing.GetPoint(i).Clone() as IPoint;
pointOnExr.Z -= height;
IPoint pointOnExr2 = null;
if (i == polygonDraw.ExteriorRing.PointCount - 1)
{
pointOnExr2 = polygonDraw.ExteriorRing.GetPoint(0).Clone() as IPoint;
}
else
{
pointOnExr2 = polygonDraw.ExteriorRing.GetPoint(i + 1).Clone() as IPoint;
}
pointOnExr2.Z -= height;
IPolygon polygonTemp = geoFactory.CreateGeometry(gviGeometryType.gviGeometryPolygon, gviVertexAttribute.gviVertexAttributeZ) as IPolygon;
polygonTemp.ExteriorRing.AppendPoint(polygonDraw.ExteriorRing.GetPoint(i));
polygonTemp.ExteriorRing.AppendPoint(pointOnExr);
polygonTemp.ExteriorRing.AppendPoint(pointOnExr2);
if (i == polygonDraw.ExteriorRing.PointCount - 1)
{
polygonTemp.ExteriorRing.AppendPoint(polygonDraw.ExteriorRing.GetPoint(0));
}
else
{
polygonTemp.ExteriorRing.AppendPoint(polygonDraw.ExteriorRing.GetPoint(i + 1));
}
polygonTemp.Close();
PolygonList.Add(polygonTemp);
//CreateRenderPolygon(polygonTemp);
}
List <IRowBuffer> list = new List <IRowBuffer>();
ISpatialFilter filter = new SpatialFilter();
filter.Geometry = polygonDraw;
filter.SpatialRel = gviSpatialRel.gviSpatialRelEnvelope;
filter.GeometryField = "Geometry";
IRowBuffer row = null;
IFdeCursor cursor = __fc.Search(filter, false);
while ((row = cursor.NextRow()) != null)
{
list.Add(row);
}
foreach (IRowBuffer r in list)
{
__fid = (Int32)r.GetValue(0);
int geometryIndex = -1;
geometryIndex = r.FieldIndex("Geometry");
int nameIndex = -1;
nameIndex = r.FieldIndex("Name");
if (geometryIndex != -1)
{
// 获取ModelPoint
IModelPoint mp = r.GetValue(geometryIndex) as IModelPoint;
// 获取Model
string modelName = mp.ModelName;
IModel m = (__fc.FeatureDataSet as IResourceManager).GetModel(modelName);
Logger.WriteMsg(LogLevel.Message, "---------开始ModelPointToTriMesh----------", DateTime.Now);
// 获取MultiTriMesh
IMultiTriMesh multiTM = gc.ModelPointToTriMesh(m, mp, false);
Logger.WriteMsg(LogLevel.Message, string.Format("---------完成ModelPointToTriMesh:{0}----------", multiTM.GeometryCount), DateTime.Now);
// 获取Name
string strName = r.GetValue(nameIndex).ToString();
for (int i = 0; i < multiTM.GeometryCount; i++)
{
ITriMesh tm = multiTM.GetGeometry(i) as ITriMesh;
//if (!tm.IsClosed)
// continue;
// 生成剖面
for (int p = 0; p < PolygonList.Count; p++)
{
IPolygon curPolygon = PolygonList[p];
Logger.WriteMsg(LogLevel.Message, string.Format("TM:{0} PG:{1}", i, p), DateTime.Now);
if (gc.GetSolidProfile(tm, curPolygon, out multiPolygon))
{
Logger.WriteMsg(LogLevel.Message, string.Format("TM:{0} PG:{1} TRUE", i, p), DateTime.Now);
if (multiPolygon != null && multiPolygon.GeometryCount > 0)
{
for (int j = 0; j < multiPolygon.GeometryCount; j++)
{
IPolygon tm2 = multiPolygon.GetGeometry(j) as IPolygon;
//RenderPolygonList.Add(CreateRenderPolygon(tm2));
ICurveSymbol cs = new CurveSymbol();
cs.Color = Color.FromArgb(0, Color.White);
ISurfaceSymbol ss = new SurfaceSymbol();
ss.BoundarySymbol = cs;
ss.Color = m.GetGroup(0).GetPrimitive(0).Material.DiffuseColor;
IRenderPolygon tmPolygon = this.axRenderControl1.ObjectManager.CreateRenderPolygon(tm2, ss, rootId);
tmPolygon.SetClientData("Name", strName);
RenderPolygonList.Add(tmPolygon);
#region 用当前选中的模型材质，生成截面模型
/*
* IModel outModel = null;
* IModelPoint outMP = null;
* gc.PolygonToModelPoint(tm2, out outModel, out outMP);
*
* if (outModel != null && outMP != null)
* {
*
* IDrawGroup dgroup = outModel.GetGroup(0);
* IDrawPrimitive dpri = dgroup.GetPrimitive(0);
* IDrawMaterial dmar = dpri.Material;
* dmar.CullMode = gviCullFaceMode.gviCullNone;
*
* //处理贴图
* //IPropertySet psOut = new PropertySet();
* string imgName = m.GetGroup(0).GetPrimitive(0).Material.TextureName;
* if (imgName != "")
* {
* dmar.WrapModeS = gviTextureWrapMode.gviTextureWrapRepeat;
* dmar.WrapModeT = gviTextureWrapMode.gviTextureWrapRepeat;
*
* //string tmpImgPath = (strMediaPath + @"\1\") + imgName + ".dds";
* //IImage imageOut = (__fc.FeatureDataSet as IResourceManager).GetImage(imgName);
* //imageOut.WriteFile(tmpImgPath);
* //psOut.SetProperty(imgName, imageOut);
* //dmar.TextureName = imgName;
*
* IImage image = (__fc.FeatureDataSet as IResourceManager).GetImage(imgName);
* this.axRenderControl1.ObjectManager.AddImage(imgName, image);
* dmar.TextureName = imgName;
*
* if (dpri.VertexArray.Length == 12)
* {
* IFloatArray texcoords = new FloatArray();
* texcoords.Append(0);
* texcoords.Append(0);
* texcoords.Append(1.0f);
* texcoords.Append(0);
* texcoords.Append(1.0f);
* texcoords.Append(1.0f);
* texcoords.Append(0);
* texcoords.Append(1.0f);
* dpri.TexcoordArray = texcoords;
* }
* else
* {
* IFloatArray texcoords = new FloatArray();
* texcoords.Append(0);
* texcoords.Append(0);
* for (int v = 3; v < dpri.VertexArray.Length - 2; )
* {
* texcoords.Append(Math.Abs(dpri.VertexArray.Get(v) - dpri.VertexArray.Get(0)));
* texcoords.Append(Math.Abs(dpri.VertexArray.Get(v+2) - dpri.VertexArray.Get(2)));
* v += 3;
* }
* dpri.TexcoordArray = texcoords;
* }
* }
* else
* {
* dmar.DepthBias = m.GetGroup(0).GetPrimitive(0).Material.DepthBias;
* dmar.DiffuseColor = m.GetGroup(0).GetPrimitive(0).Material.DiffuseColor;
* dmar.EnableBlend = m.GetGroup(0).GetPrimitive(0).Material.EnableBlend;
* dmar.EnableLight = m.GetGroup(0).GetPrimitive(0).Material.EnableLight;
* dmar.SpecularColor = m.GetGroup(0).GetPrimitive(0).Material.SpecularColor;
* dmar.WrapModeS = m.GetGroup(0).GetPrimitive(0).Material.WrapModeS;
* dmar.WrapModeT = m.GetGroup(0).GetPrimitive(0).Material.WrapModeT;
* }
*
* dpri.Material = dmar;
* dgroup.SetPrimitive(0, dpri);
* outModel.SetGroup(0, dgroup);
*
* //处理模型
* //string osgNameOut = __fid + "_" + p + "_" + j;
* //string tmpOSGPath = (strMediaPath + @"\1\" + osgNameOut + ".osg");
* //outModel.WriteFile(tmpOSGPath, psOut);
* //outMP.ModelName = tmpOSGPath;
*
*
* string osgName = __fid + "_" + p + "_" + j;
* this.axRenderControl1.ObjectManager.AddModel(osgName, outModel);
* outMP.ModelName = osgName;
*
* //可视化
* outMP.ModelEnvelope = outModel.Envelope;
* IRenderModelPoint outRenderMP = this.axRenderControl1.ObjectManager.CreateRenderModelPoint(outMP, null, rootId);
* //this.axRenderControl1.RefreshModel(null, tmpOSGPath);
* outRenderMP.ShowOutline = true;
* //outRenderMP.VisibleMask = gviViewportMask.gviViewNone;
* RenderMPProfileList.Add(outRenderMP);
* }
*/
#endregion
}
}
}
} //遍历PolygonList结束
} //遍历multiTM结束
//去掉内环模型
__fl.VisibleMask = gviViewportMask.gviViewNone;
IModel modelInterior = null;
IModelPoint mpInterior = null;
IRenderModelPoint rmpInterior = null;
IModel modelExterior = null;
IModelPoint mpExterior = null;
IRenderModelPoint rmpExterior = null;
IMultiPolygon mltiPolygon = geoFactory.CreateGeometry(gviGeometryType.gviGeometryMultiPolygon, gviVertexAttribute.gviVertexAttributeZ) as IMultiPolygon;
mltiPolygon.AddGeometry(polygonDraw);
Logger.WriteMsg(LogLevel.Message, "---------开始SplitModelPointByPolygon2DWithZ----------", DateTime.Now);
if (gc.SplitModelPointByPolygon2DWithZ(mltiPolygon, m, mp, minZ, maxZ, out modelInterior, out mpInterior, out modelExterior, out mpExterior))
{
Logger.WriteMsg(LogLevel.Message, "SplitModelPointByPolygon2DWithZ返回值为TRUE", DateTime.Now);
if (modelExterior != null && mpExterior != null)
{
this.axRenderControl1.ObjectManager.AddModel(mpExterior.ModelName, modelExterior);
string[] imagenames = modelExterior.GetImageNames();
for (int q = 0; q < imagenames.Length; q++)
{
IImage image = (__fc.FeatureDataSet as IResourceManager).GetImage(imagenames[q]);
this.axRenderControl1.ObjectManager.AddImage(imagenames[q], image);
}
rmpExterior = this.axRenderControl1.ObjectManager.CreateRenderModelPoint(mpExterior, null, rootId);
RenderMPExteriorList.Add(rmpExterior);
}
if (modelInterior != null && mpInterior != null)
{
this.axRenderControl1.ObjectManager.AddModel(mpInterior.ModelName, modelInterior);
string[] imagenames = modelInterior.GetImageNames();
for (int q = 0; q < imagenames.Length; q++)
{
IImage image = (__fc.FeatureDataSet as IResourceManager).GetImage(imagenames[q]);
this.axRenderControl1.ObjectManager.AddImage(imagenames[q], image);
}
rmpInterior = this.axRenderControl1.ObjectManager.CreateRenderModelPoint(mpInterior, null, rootId);
rmpInterior.VisibleMask = gviViewportMask.gviViewNone;
RenderMPInteriorList.Add(rmpInterior);
}
}
}
}//遍历RowBufferList结束
//MessageBox.Show("执行成功");
}
private void btnAnalyse_Click(object sender, EventArgs e)
{
try
{
ClearRes();
WaitForm.Start("正在分析...", "请稍后");
if (this._rLine == null)
{
return;
}
List <DF3DFeatureClass> list = DF3DFeatureClassManager.Instance.GetAllFeatureClass();
if (list == null || list.Count == 0)
{
return;
}
IPolyline line = this._rLine.GetFdeGeometry() as IPolyline;
ILine l = this._geoFactory.CreateGeometry(gviGeometryType.gviGeometryLine, gviVertexAttribute.gviVertexAttributeZ) as ILine;
l.StartPoint = line.StartPoint;
l.EndPoint = line.EndPoint;
ISpatialFilter filter = new SpatialFilter();
filter.SpatialRel = gviSpatialRel.gviSpatialRelEnvelope;
filter.GeometryField = "Geometry";
List <InterPt> listInterPts = new List <InterPt>();
foreach (DF3DFeatureClass dffc in list)
{
IFeatureClass fc = dffc.GetFeatureClass();
if (fc == null)
{
continue;
}
IFeatureLayer fl = dffc.GetFeatureLayer();
if (fl != null)
{
if (fl.VisibleMask == gviViewportMask.gviViewNone)
{
continue;
}
}
int indexGeo = fc.GetFields().IndexOf("Geometry");
if (indexGeo == -1)
{
continue;
}
filter.Geometry = l;
IFdeCursor cursor = null;
IRowBuffer row = null;
try
{
cursor = fc.Search(filter, false);
while ((row = cursor.NextRow()) != null)
{
if (!row.IsNull(indexGeo))
{
IGeometry geo = row.GetValue(indexGeo) as IGeometry;
if (geo.GeometryType == gviGeometryType.gviGeometryModelPoint)
{
IModelPoint modelPoint = geo as IModelPoint;
IModel model = (fc.FeatureDataSet as IResourceManager).GetModel(modelPoint.ModelName);
IGeometryConvertor gc = new GeometryConvertor();
IMultiTriMesh triMesh = gc.ModelPointToTriMesh(model, modelPoint, false);
if (triMesh != null)
{
IVector3 v3 = triMesh.LineSegmentIntersect(l);
if (v3 != null)
{
IPoint pttemp = this._geoFactory.CreatePoint(gviVertexAttribute.gviVertexAttributeZ);
pttemp.X = v3.X;
pttemp.Y = v3.Y;
pttemp.Z = v3.Z;
double dis = (pttemp.X - line.StartPoint.X) * (pttemp.X - line.StartPoint.X) +
(pttemp.Y - line.StartPoint.Y) * (pttemp.Y - line.StartPoint.Y) +
(pttemp.Z - line.StartPoint.Z) * (pttemp.Z - line.StartPoint.Z);
InterPt ip = new InterPt();
ip.pt = pttemp;
ip.dis = dis;
listInterPts.Add(ip);
}
}
}
}
}
}
catch (Exception ex)
{
}
finally
{
if (row != null)
{
System.Runtime.InteropServices.Marshal.ReleaseComObject(row);
row = null;
}
if (cursor != null)
{
System.Runtime.InteropServices.Marshal.ReleaseComObject(cursor);
cursor = null;
}
}
}
if (listInterPts.Count > 0)
{
this._rLine.VisibleMask = gviViewportMask.gviViewNone;
listInterPts.Sort(new CmpInterPt());
IPoint pt = listInterPts[0].pt;
IPolyline lineStart = this._geoFactory.CreateGeometry(gviGeometryType.gviGeometryPolyline, gviVertexAttribute.gviVertexAttributeZ) as IPolyline;
lineStart.AppendPoint(line.StartPoint);
lineStart.AppendPoint(pt);
ICurveSymbol csS = new CurveSymbolClass();
csS.Color = 0xff0000ff;
csS.Width = -2;
IRenderPolyline rLineStart = this._3DControl.ObjectManager.CreateRenderPolyline(lineStart, csS, this._3DControl.ProjectTree.RootID);
this._listRes.Add(rLineStart.Guid);
IPolyline lineEnd = this._geoFactory.CreateGeometry(gviGeometryType.gviGeometryPolyline, gviVertexAttribute.gviVertexAttributeZ) as IPolyline;
lineEnd.AppendPoint(line.EndPoint);
lineEnd.AppendPoint(pt);
ICurveSymbol csE = new CurveSymbolClass();
csE.Color = 0xffff0000;
csE.Width = -2;
IRenderPolyline rLineEnd = this._3DControl.ObjectManager.CreateRenderPolyline(lineEnd, csE, this._3DControl.ProjectTree.RootID);
this._listRes.Add(rLineEnd.Guid);
}
else
{
this._rLine.VisibleMask = gviViewportMask.gviViewNone;
ICurveSymbol csS = new CurveSymbolClass();
csS.Color = 0xff0000ff;
csS.Width = -2;
IRenderPolyline rLine = this._3DControl.ObjectManager.CreateRenderPolyline(line, csS, this._3DControl.ProjectTree.RootID);
this._listRes.Add(rLine.Guid);
}
}
catch (Exception ex)
{
}
finally
{
WaitForm.Stop();
}
}
public void SelectFeaturesFromBaseLyr(IPolyline polyLine)
{
btn_analyse.Enabled = false;
btn_createLine.Enabled = false;
this.btnFlyToSourcePoint.Enabled = false;
this.btnFlyToTargetPoint.Enabled = false;
IFdeCursor cursor = null;
try
{
this.dataGridView1.Rows.Clear();
axRenderControl1.FeatureManager.UnhighlightAll();
IRowBuffer row = null;
int index = 0;
List <IRowBuffer> list = new List <IRowBuffer>();
this.Text = "开始粗查询";
foreach (IFeatureClass fc in fcMap.Keys)
{
ISpatialFilter filter = new SpatialFilter();
filter.Geometry = polyline;
filter.SpatialRel = gviSpatialRel.gviSpatialRelEnvelope;
filter.GeometryField = "Geometry";
cursor = fc.Search(filter, false);
while ((row = cursor.NextRow()) != null)
{
list.Add(row);
}
this.Text = "开始细查询";
foreach (IRowBuffer r in list)
{
index++;
int geometryIndex = -1;
geometryIndex = r.FieldIndex("Geometry");
if (geometryIndex != -1)
{
IModelPoint modelPoint = r.GetValue(geometryIndex) as IModelPoint;
IModel model = (fc.FeatureDataSet as IResourceManager).GetModel(modelPoint.ModelName);
IGeometryConvertor gc = new GeometryConvertor();
this.Text = "正在计算第" + index.ToString() + "个IMultiTriMesh是否与线相交";
IMultiTriMesh triMesh = gc.ModelPointToTriMesh(model, modelPoint, false);
if (triMesh != null)
{
ILine l = geoFactory.CreateGeometry(gviGeometryType.gviGeometryLine, gviVertexAttribute.gviVertexAttributeZ) as ILine;
l.StartPoint = polyline.StartPoint;
l.EndPoint = polyline.EndPoint;
IVector3 v3 = triMesh.LineSegmentIntersect(l);
if (v3 != null)
{
string fid = "";
string fName = "";
string groupId = "";
IEnvelope env = null;
int fidPos = r.FieldIndex(fc.FidFieldName);
axRenderControl1.FeatureManager.HighlightFeature(fc, int.Parse(r.GetValue(fidPos).ToString()), System.Drawing.Color.Red);
for (int i = 0; i < r.FieldCount; i++)
{
string fieldName = r.Fields.Get(i).Name;
if (r.Fields.Get(i).Name == "oid")
{
fid = r.GetValue(i).ToString();
}
else if (r.Fields.Get(i).Name == "groupid")
{
groupId = r.GetValue(i).ToString();
}
else if (r.Fields.Get(i).Name == "Name")
{
fName = r.GetValue(i).ToString();
}
else if (r.Fields.Get(i).Name == "Geometry")
{
IGeometry geometry = r.GetValue(i) as IModelPoint;
env = geometry.Envelope;
}
}
RowObject ro = new RowObject()
{
FID = fid, GroupId = groupId, Name = fName, FeatureClass = fc, Envelop = env
};
if (!rowMap.ContainsKey(ro.FID))
{
rowMap.Add(ro.FID, ro);
}
}
}
}
} // end of foreach (IRowBuffer r in list)
} // end of foreach (IFeatureClass fc in fcMap.Keys)
this.Text = "通视分析完成！";
LoadGridView();
}
catch (Exception ex)
{
if (ex.GetType().Name.Equals("UnauthorizedAccessException"))
{
MessageBox.Show("需要标准runtime授权");
}
else
{
MessageBox.Show(ex.Message);
}
}
finally
{
this.btn_analyse.Enabled = true;
this.btn_createLine.Enabled = true;
this.btnFlyToSourcePoint.Enabled = true;
this.btnFlyToTargetPoint.Enabled = true;
if (cursor != null)
{
//Marshal.ReleaseComObject(cursor);
cursor = null;
}
}
}
