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; } } }