/// <summary>
/// 显示轮廓线 张琪 20110629
/// </summary>
/// <param name="pPolygon"></param>
private void ShowCountour(IPolygon pPolygon)
{
SysCommon.CProgress vProgress = new SysCommon.CProgress("进度条");
try
{
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, ContourName);
if (!chkShowContour.Checked)
{
m_pCurrentSceneControl.SceneGraph.RefreshViewers();
return;
}
vProgress.EnableCancel = false;//设置进度条
vProgress.ShowDescription = true;
vProgress.FakeProgress = true;
vProgress.TopMost = true;
vProgress.ShowProgress();
vProgress.SetProgress("正在绘制轮廓线");
ITopologicalOperator pTopo = pPolygon as ITopologicalOperator;
IGeometry pPolyLine = new PolylineClass();
pPolyLine = pTopo.Boundary;
object StepSize = Type.Missing;;
m_SurFace.InterpolateShape(pPolyLine as IGeometry, out pPolyLine, ref StepSize);
IGroupElement pGroup = null;
//用于绘制三维效果
Cls3DMarkDraw.AddSimpleGraphic(pPolyLine as IGeometry, Cls3DMarkDraw.getRGB(30, 255, 255), 4, ContourName, m_pCurrentSceneControl.Scene, pGroup);
m_pCurrentSceneControl.SceneGraph.RefreshViewers();
vProgress.Close();
}
catch
{
vProgress.Close();
}
}
private void frm3DSection_FormClosing(object sender, FormClosingEventArgs e)
{
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, sPolyOutlineName);//清除绘制的要素
m_SurFace = null;
m_pCurrentSceneControl.CurrentTool = null;
this.Dispose();
}
private void frm3DVolumeAreaSta_FormClosing(object sender, FormClosingEventArgs e)
{
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, sRefPlaneName);
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, ContourName);
m_pCurrentSceneControl.CurrentTool = null;//释放使用工具
m_ScenePoly = null;
m_StaPolygon = null;
m_SurFace = null;
this.Dispose();
}
/// <summary>
/// 显示参考面 张琪 20110629
/// </summary>
/// <param name="pPoly">参考多边形</param>
private void ShowRefPlane(IPolygon pPoly)
{
SysCommon.CProgress vProgress = new SysCommon.CProgress("进度条");
try
{
IGeometry pBoundary = new PolylineClass();
pBoundary = null;
IGeometry pPolygon = new PolygonClass();
pPolygon = null;
object StepSize = Type.Missing;
m_SurFace.InterpolateShape(pPoly as IGeometry, out pPolygon, ref StepSize);
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, sRefPlaneName); //清除已有的参考面
if (!chkShowLWRP.Checked) //当为false时不显示参考面
{
m_pCurrentSceneControl.SceneGraph.RefreshViewers();
return;
}
vProgress.EnableCancel = false;//设置进度条
vProgress.ShowDescription = true;
vProgress.FakeProgress = true;
vProgress.TopMost = true;
vProgress.ShowProgress();
vProgress.SetProgress("正在绘制参考面");
ITopologicalOperator pTopo = pPolygon as ITopologicalOperator;
pBoundary = pTopo.Boundary as IPolyline;
//设置Z值方向显示
IZAware pZAware = pBoundary as IZAware;
pZAware.ZAware = true;
m_SurFace.InterpolateShape(pBoundary as IGeometry, out pBoundary, ref StepSize);
m_SurFace.InterpolateShape(pPoly as IGeometry, out pPolygon, ref StepSize);
IExtrude pExtrude = new GeometryEnvironmentClass();
//获取参考面的多面体
IMultiPatch pMultiPatchRefPlaneWall = pExtrude.ExtrudeAbsolute(Convert.ToDouble(txtPlaneHeight.Text), pBoundary as IGeometry) as IMultiPatch;
pMultiPatchRefPlaneWall.SpatialReference = m_pCurrentSceneControl.Scene.SpatialReference;
IZ pZ;
pZAware = pPolygon as IZAware;
pZAware.ZAware = true;
pZ = pPolygon as IZ;
pZ.SetConstantZ(Convert.ToDouble(txtPlaneHeight.Text));
//用于绘制三维效果
IGroupElement pGroup = null;
Cls3DMarkDraw.AddSimpleGraphic(pMultiPatchRefPlaneWall as IGeometry, Cls3DMarkDraw.getRGB(71, 61, 255), 1, sRefPlaneName, m_pCurrentSceneControl.Scene, pGroup);
Cls3DMarkDraw.AddSimpleGraphic(pPolygon as IGeometry, Cls3DMarkDraw.getRGB(71, 61, 255), 1, sRefPlaneName, m_pCurrentSceneControl.Scene, pGroup);
m_pCurrentSceneControl.SceneGraph.RefreshViewers();
vProgress.Close();
}
catch
{
vProgress.Close();
}
}
private void DrawGeo_Click(object sender, EventArgs e)
{
//新建三维绘制工具
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, sPolyOutlineName);
m_pCurrentSceneControl.SceneGraph.RefreshViewers();
pTool3DDrawGeo.EndDtrawd += new myEventHandler(pTool3DDrawGeo_EndDtrawd);
pTool3DDrawGeo.GeometryType = ESRI.ArcGIS.Geometry.esriGeometryType.esriGeometryLine;
pTool3DDrawGeo.pSurface = m_SurFace;
ICommand pCommad = pTool3DDrawGeo;
pCommad.OnCreate(m_pCurrentSceneControl.Object);
m_pCurrentSceneControl.CurrentTool = pCommad as ITool;
if (this.WriteLog)
{
Plugin.LogTable.Writelog("两点间剖面分析,表面集为:" + comboBoxOpen.Text);//xisheng 日志记录07.08
}
}
/// <summary>
/// 实例化绘制范围工具
/// </summary>
/// <param name="sender"></param>
/// <param name="e"></param>
private void DrawGeo_Click(object sender, EventArgs e)
{
btncount.Enabled = false;
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, sRefPlaneName);
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, ContourName);
m_pCurrentSceneControl.SceneGraph.RefreshViewers();
//pTool3DDrawGeo.BeginDrawed += new myEventHandler(pTool3DDrawGeo_BeginDrawed);
pTool3DDrawGeo.EndDtrawd += new myEventHandler(pTool3DDrawGeo_EndDtrawd);
pTool3DDrawGeo.GeometryType = ESRI.ArcGIS.Geometry.esriGeometryType.esriGeometryPolygon;
pTool3DDrawGeo.pSurface = m_SurFace;
ICommand pCommad = pTool3DDrawGeo;
pCommad.OnCreate(m_pCurrentSceneControl.Object);
m_pCurrentSceneControl.CurrentTool = pCommad as ITool;
if (this.WriteLog)
{
Plugin.LogTable.Writelog("局部填挖方分析,表面集为:" + comboBoxOpen.Text + ",基准面高度为:" + txtPlaneHeight.Text);
}
}
/// <summary>
/// 双击表示绘制对象结束
/// </summary>
public override void OnDblClick()
{
ITopologicalOperator pTopoOp;
if (m_pPointColl != null)
{
switch (m_DrawType.ToString())
{
case "esriGeometryLine":
IPointCollection pPolyLine = new PolylineClass();
pPolyLine.AddPointCollection(m_pPointColl);
pTopoOp = pPolyLine as ITopologicalOperator;
pTopoOp.Simplify();
m_Geometry = pPolyLine as IGeometry;
m_Geometry.SpatialReference = m_sceneHookHelper.Scene.SpatialReference;
Cls3DMarkDraw.DeleteAllElementsWithName(m_sceneHookHelper.Scene, sPolyOutlineName);
EndDtrawd(true); //触发结束事件
break;
case "esriGeometryPolygon":
if (m_pPointColl.PointCount < 3)
{
return;
}
IPointCollection pPolygon = new PolygonClass();
pPolygon.AddPointCollection(m_pPointColl);
pTopoOp = pPolygon as ITopologicalOperator;
pTopoOp.Simplify();
m_Geometry = pPolygon as IGeometry;
m_Geometry.SpatialReference = m_sceneHookHelper.Scene.SpatialReference;
Cls3DMarkDraw.DeleteAllElementsWithName(m_sceneHookHelper.Scene, sPolyOutlineName);
EndDtrawd(true); //触发结束事件
break;
}
m_pPointColl = null;
}
//Cls3DMarkDraw.DeleteAllElementsWithName(m_sceneHookHelper.Scene, sPolyOutlineName);
m_sceneHookHelper.SceneGraph.RefreshViewers();
}
/// <summary>
/// 获取点击地图上的点
/// </summary>
/// <param name="Button"></param>
/// <param name="Shift"></param>
/// <param name="X"></param>
/// <param name="Y"></param>
public override void OnMouseDown(int Button, int Shift, int X, int Y)
{
try
{
ESRI.ArcGIS.Geometry.IPoint pMapPoint = new ESRI.ArcGIS.Geometry.PointClass();;
ISceneGraph pSceneGraph = m_sceneHookHelper.SceneGraph;
object pOwner;
object pObject;
object before = Type.Missing;
object after = Type.Missing;
object StepSize = Type.Missing;
IDisplay3D pDisplay;
pSceneGraph.Locate(pSceneGraph.ActiveViewer, X, Y, esriScenePickMode.esriScenePickGeography, true, out pMapPoint, out pOwner, out pObject); //获取鼠标点击的位置并转化为地理坐标
if (pMapPoint == null)
{
return;
}
pMapPoint.Z = pMapPoint.Z / m_sceneHookHelper.Scene.ExaggerationFactor;
pMapPoint.Z = m_psurface.GetElevation(pMapPoint);
pMapPoint.SpatialReference = pSceneGraph.Scene.SpatialReference;
pDisplay = m_sceneHookHelper.SceneGraph as IDisplay3D;
pDisplay.FlashLocation(pMapPoint);//闪烁显示被点击的位置
IGeometry pGeom = null;
Cls3DMarkDraw.DeleteAllElementsWithName(m_sceneHookHelper.Scene, sPolyOutlineName);
//根据绘制对象类型的不同定义不同的类型
switch (m_DrawType.ToString())
{
case "esriGeometryPoint":
m_Geometry = pMapPoint;
break;
case "esriGeometryLine":
if (m_pPointColl == null)
{
m_pPointColl = new PolylineClass();
pGeom = new PolylineClass();
}
m_pPointColl.AddPoint(pMapPoint, ref before, ref after);
break;
case "esriGeometryPolygon":
if (m_pPointColl == null)
{
m_pPointColl = new PolygonClass();
pGeom = new PolygonClass();
}
m_pPointColl.AddPoint(pMapPoint, ref before, ref after);
break;
}
//BeginDrawed(true);
IGroupElement pGroup = null;
if (m_pPointColl.PointCount == 1)
{
//当为一个点时绘制点
Cls3DMarkDraw.AddSimpleGraphic(pMapPoint, Cls3DMarkDraw.getRGB(71, 61, 255), 4, sPolyOutlineName, m_sceneHookHelper.Scene, pGroup);
}
else if (m_DrawType.ToString() == "esriGeometryLine")
{
pGeom = m_pPointColl as IGeometry;
pGeom.SpatialReference = pMapPoint.SpatialReference;
m_psurface.InterpolateShape(pGeom, out pGeom, ref StepSize);
Cls3DMarkDraw.AddSimpleGraphic(pGeom, Cls3DMarkDraw.getRGB(71, 61, 255), 4, sPolyOutlineName, m_sceneHookHelper.Scene, pGroup);
m_pPointColl = pGeom as IPointCollection;
}
else
{
ITopologicalOperator pTopo = m_pPointColl as ITopologicalOperator;
pGeom = pTopo.Boundary;
pGeom.SpatialReference = pMapPoint.SpatialReference;
m_psurface.InterpolateShape(pGeom, out pGeom, ref StepSize);
Cls3DMarkDraw.AddSimpleGraphic(pGeom, Cls3DMarkDraw.getRGB(71, 61, 255), 4, sPolyOutlineName, m_sceneHookHelper.Scene, pGroup);
}
m_sceneHookHelper.SceneGraph.RefreshViewers();
}
catch
{
return;
}
}
/// <summary>
/// 获取分析图层及Surfer
/// </summary>
/// <param name="sender"></param>
/// <param name="e"></param>
private void comboBoxOpen_SelectedIndexChanged(object sender, EventArgs e)
{
try
{
if (comboBoxOpen.Items.Count > 0)
{
label.Text = "";
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, sRefPlaneName);
Cls3DMarkDraw.DeleteAllElementsWithName(m_pCurrentSceneControl.Scene, ContourName);
m_pCurrentSceneControl.SceneGraph.RefreshViewers();
m_ScenePoly = null;
m_StaPolygon = null;
m_SurFace = null;
txtArea2DSel.Text = "";
txtAreaAbove.Text = "";
txtAreaBelow.Text = "";
txtPara.Text = "";
txtPlaneHeight.Text = "";
txtVolumeAbove.Text = "";
txtVolumeBelow.Text = "";
DrawGeo.Enabled = false;
sldPlaneHeight.Enabled = false;
btncount.Enabled = false;
string LayerName = comboBoxOpen.SelectedItem.ToString();
ILayer pLayer = GetLayerByName(ref LayerName);
if (pLayer is IRasterLayer)//读取Raster数据的ISurface
{
IRasterLayer pRasterLayer = pLayer as IRasterLayer;
IGeoDataset pDataset = pRasterLayer as IGeoDataset;
ISpatialReference psr = pDataset.SpatialReference;
if (!(psr is IProjectedCoordinateSystem))
{
label.Text = "该数据为地理坐标，计算结果无效仅供参考！";
}
IRasterSurface pRasterSurface = new RasterSurfaceClass();
pRasterSurface.PutRaster(pRasterLayer.Raster, 0);
m_SurFace = pRasterSurface as ISurface;
DrawGeo.Enabled = false;
m_IsTin = false;
}
else if (pLayer is ITinLayer)//读取TIN数据的ISurface
{
ITinLayer pTinLayer = pLayer as ITinLayer;
IGeoDataset pDataset = pTinLayer as IGeoDataset;
ISpatialReference psr = pDataset.SpatialReference;
if (!(psr is IProjectedCoordinateSystem))
{
label.Text = "该数据为地理坐标，计算结果无效仅供参考！";
}
m_SurFace = pTinLayer.Dataset as ISurface;
m_IsTin = true;
DrawGeo.Enabled = true;
}
if (m_SurFace == null)
{
MessageBox.Show("请选择需要进行分析的DEM数据", "提示！");
return;
}
if (m_StaPolygon == null)
{
if (m_SurFace is ITin)
{
ITin pTin = m_SurFace as ITin;
sldPlaneHeight.Maximum = Convert.ToInt32(pTin.Extent.ZMax);
sldPlaneHeight.Minimum = Convert.ToInt32(pTin.Extent.ZMin);
}
else
{
IRasterSurface pRasterSurface = m_SurFace as IRasterSurface;
IRasterStatistics pRasterStatistic = pRasterSurface.RasterBand.Statistics;
sldPlaneHeight.Maximum = Convert.ToInt32(pRasterStatistic.Maximum);
sldPlaneHeight.Minimum = Convert.ToInt32(pRasterStatistic.Minimum);
}
double Vale = sldPlaneHeight.Maximum + sldPlaneHeight.Minimum;
Vale = Vale / 2;
txtPlaneHeight.Text = Vale.ToString();
Vale = Convert.ToDouble(txtPlaneHeight.Text);
sldPlaneHeight.Value = Convert.ToInt32(Vale);
}
else
{
double zrange;
zrange = m_StaPolygon.Envelope.ZMax - m_StaPolygon.Envelope.ZMin;
sldPlaneHeight.Maximum = Convert.ToInt32(m_StaPolygon.Envelope.ZMax + zrange);
sldPlaneHeight.Minimum = Convert.ToInt32(m_StaPolygon.Envelope.ZMin);
sldPlaneHeight.Value = Convert.ToInt32(Math.Round(m_StaPolygon.Envelope.ZMin + (zrange * 0.5)));
txtPlaneHeight.Text = sldPlaneHeight.Value.ToString();
}
sldPlaneHeight.Enabled = true;
btncount.Enabled = true;
}
}
catch
{
return;
}
}