ArcGIS Engine生成等值线(C#)

转自ArcGIS Engine生成等值线(C#) - yycGIS - 博客园

本文介绍c#写的利用ArcGIS Engine生成等值线的方法。

c#写的根据雨量站的降雨量值内插出降雨量等值线的功能。做几点说明:
根据离散点生成等值线:
(1)判断等值线shapefile图层临时存放路径是否存在,如果不存在则创建,如果存在并不为空则删除重建;
(2)加【降雨量】字段并赋值,生成等值线后删除所加字段;
(3)根据所选的字段利用IDW内插生成raster,再设定间距生成contour并标注。

---------------------------------以下是源代码----------------------------------

private void button1_Click(object sender, EventArgs e)
        {            
            //定义等值线临时shapefile图层存放路径
            string ls_TempSavePath = @"C:\HNYTTEMP" ;
            if (!Directory.Exists(ls_TempSavePath))//如果目录不存在,则创建
                Directory.CreateDirectory(ls_TempSavePath);
            DirectoryInfo di = new DirectoryInfo(ls_TempSavePath);
            FileInfo[] fi = null;
            fi = di.GetFiles();
            if (fi.Length != 0)//如果目录不为空,则删除目录,然后重建
            {undefined
                Directory.Delete(ls_TempSavePath, true);
                Directory.CreateDirectory(ls_TempSavePath);
            }          
             IMap pMap  = axMapControl1.Map;  
             IInterpolationOp pInterpolationOp  = new RasterInterpolationOpClass();
  
             // Create the input point object
             IGeoDataset pInputDataset ;
             //雨量站图层【点】
             IFeatureLayer pFeatLayer = axMapControl1.get_Layer(0) as IFeatureLayer;
             // Calls function to open the point dataset from disk 
  
             IFeatureClass pFeatCla  = pFeatLayer.FeatureClass;  
            
             IFeature pFeature ;  
             IField pField  = new FieldClass();
             IFieldEdit pFieldEdit = pField as IFieldEdit;
  
             pFieldEdit .Name_2 = "降雨量";
             pFieldEdit.Type_2 = esriFieldType.esriFieldTypeDouble;
             pFieldEdit.Length_2 = 12;
             pFieldEdit.IsNullable_2 = false;
             pFieldEdit.DefaultValue_2 = 1111;
            pFeatCla.AddField(pFieldEdit);   //添加那个要用的字段进去
        for(int ii=0;ii<pFeatCla.FeatureCount(null) - 1;ii++)
            {undefined
               pFeature = pFeatCla.GetFeature(ii);
               pFeature.set_Value(pFeature.Fields.FindField("降雨量"), 22+1999/(ii+2)); 
              //要素字段赋值,这里是自己随便赋的值,仅做测试用
          pFeature.Store();
            } 
           pInputDataset = pFeatLayer as IGeoDataset;
           // Define the search radius
           IRasterRadius pRadius  = new RasterRadiusClass();
           
           object Missing = Type.Missing;
           pRadius.SetVariable(12, ref Missing);  
           //Create FeatureClassDescriptor using a value field
          IFeatureClassDescriptor pFCDescriptor = new FeatureClassDescriptorClass();
          pFCDescriptor.Create( pFeatLayer.FeatureClass, null, "降雨量"); 
          //Set cellsize for output raster in the environment
          object cellSizeProvider = 20;
    
         IRasterAnalysisEnvironment pEnv  = pInterpolationOp as IRasterAnalysisEnvironment;
         pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref  cellSizeProvider);
     
         //Perform the interpolation
         IRaster pOutRaster = pInterpolationOp.IDW(pFCDescriptor as IGeoDataset, 2, pRadius, ref Missing) as IRaster;
         
         //Add output into ArcMap as a raster layer    
         RasterLayer pOutRasLayer = new RasterLayerClass();
         pOutRasLayer.CreateFromRaster(pOutRaster);
         pOutRasLayer.Name = "栅格";
         //pMap.AddLayer(pOutRasLayer);
    
         IGeoDataset pGeoDataSet = pOutRaster as IGeoDataset;
         IWorkspaceFactory pWorkspaceFactory  = new ShapefileWorkspaceFactory();
         IWorkspace pShpWorkspace = pWorkspaceFactory.OpenFromFile(ls_TempSavePath, 0);
         ISurfaceOp2 pSurfaceOp2 = new RasterSurfaceOpClass();
         IRasterAnalysisEnvironment pRasterAnalysisEnvironment = pSurfaceOp2  as IRasterAnalysisEnvironment;
         pRasterAnalysisEnvironment.Reset();
         pRasterAnalysisEnvironment.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, ref cellSizeProvider);
         pRasterAnalysisEnvironment.OutWorkspace = pShpWorkspace;
         double  dInterval   =5;  //间距 
         IGeoDataset pOutputDataSet   = pSurfaceOp2.Contour(pGeoDataSet, dInterval,ref Missing,ref Missing);
         IFeatureClass pFeatureClass=  pOutputDataSet as IFeatureClass;
         IFeatureLayer pFeatureLayer = new FeatureLayerClass();
         pFeatureLayer.FeatureClass = pFeatureClass;
         IGeoFeatureLayer pGeoFeatureLayer  = pFeatureLayer as IGeoFeatureLayer;
         pGeoFeatureLayer.DisplayAnnotation = true;
         pGeoFeatureLayer.DisplayField = "Contour";
         pGeoFeatureLayer.Name = "降雨量等值线";
         pMap.AddLayer( pGeoFeatureLayer);
         axMapControl1.Refresh();
        
         //删除【降雨量】字段
      IFields pFields = pFeatCla.Fields;
         int lFieldnumber = pFields.FindField("降雨量");
         IField pField1 = pFields.get_Field(lFieldnumber);
         pFeatCla.DeleteField(pField);
        }
acelee

  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
ArcGIS Engine 10.1是一个地理信息系统的开发平台,可以用来构建和定制地理空间应用程序。绘制多面体(Polyhedron)是其中的一个功能。 在ArcGIS Engine 10.1中,绘制多面体可以通过几何对象类(Geometry Class)来实现。可以使用Polygon或Multipatch类来创建多边形或多面体。 要绘制多面体,首先需要确定多面体的坐标和形状。可以通过定义多边形的顶点坐标来创建一个Polygon实例,或者通过使用多边形的外观和高度属性来创建一个Multipatch实例。 对于Polygon类,需要提供一个点集合,并按照顺时针或逆时针方向定义多边形的边界。可以使用Polygon的AddPoint方法将点逐一添加到多边形中,然后调用Finish方法来完成多边形的创建。 对于Multipatch类,除了点集合外,还需要定义多面体的外观和高度属性。可以使用Multipatch的AddRing方法将多个环添加到多面体中,然后使用AddPatch方法将这些环连接起来,并使用Height属性来设置多面体的高度。 完成多面体的创建后,可以使用绘制方法将其显示在地图中。可以调用GraphicsLayer的AddGraphic方法将多边形或多面体添加到图形层中,然后在地图视图中刷新图层,即可实现多面体的绘制。 绘制多面体是ArcGIS Engine 10.1中的一个常见需求,可以通过上述的步骤来完成。同时,还可以根据具体的需求进行更复杂的定制,如添加纹理、边框等效果来美化多面体的显示效果。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值