• 设为首页
  • 点击收藏
  • 手机版
    手机扫一扫访问
    迪恩网络手机版
  • 关注官方公众号
    微信扫一扫关注
    公众号

c#调用ArcEngine的GP工具

原作者: [db:作者] 来自: [db:来源] 收藏 邀请

转自原文c# 调用ArcEngine的GP工具AE调用GP工具

 

IAoInitialize m_AoInitialize = new AoInitializeClass();   
esriLicenseStatuslicenseStatus = esriLicenseStatus.esriLicenseUnavailable;   
licenseStatus = m_AoInitialize.Initialize(esriLicenseProductCode.esriLicenseProductCodeArcInfo);  
if (licenseStatus == esriLicenseStatus.esriLicenseCheckedOut)  
{  
    GeoprocessortGp = newGeoprocessor();  
    tGp.OverwriteOutput = true;   
    licenseStatus = m_AoInitialize.IsExtensionCodeAvailable(esriLicenseProductCode.esriLicenseProductCodeEngine,esriLicenseExtensionCode.esriLicenseExtensionCodeSpatialAnalyst);  
    licenseStatus = m_AoInitialize.CheckOutExtension(esriLicenseExtensionCode.esriLicenseExtensionCodeSpatialAnalyst);  
      
    //LICenced arcengine  
    ESRI.ArcGIS.ConversionTools.ASCIIToRastertASC = newESRI.ArcGIS.ConversionTools.ASCIIToRaster();  
    tASC.data_type = "FLOAT";  
    tASC.in_ascii_file = pWorkSpacePath + @"a1.txt";  
    tASC.out_raster = pWorkSpacePath + "a1.img";  
    //ScrollToBottom("ASCIIToRaster");  
    IGeoProcessorResulttGeoResult = (IGeoProcessorResult)tGp.Execute(tASC, null);  
    if (tGeoResult.Status == ESRI.ArcGIS.esriSystem.esriJobStatus.esriJobSucceeded)  
    {  
         ESRI.ArcGIS.SpatialAnalystTools.ReclassifytReclass = newESRI.ArcGIS.SpatialAnalystTools.Reclassify();  
        tReclass.in_raster = pWorkSpacePath + "a1.img";  
        tReclass.missing_values = "NODATA";  
        tReclass.out_raster = pWorkSpacePath + @"Out3.img";  
        tReclass.reclass_field = "VALUE";  
        tReclass.remap = "0 958 10;958 988 20;988 990 30;990 1000 40;1000 1100 50";  
       // ScrollToBottom("Reclassify");  
        tGeoResult = (IGeoProcessorResult)tGp.Execute(tReclass, null);  
        if (tGeoResult.Status == ESRI.ArcGIS.esriSystem.esriJobStatus.esriJobSucceeded)  
        {  
            ESRI.ArcGIS.ConversionTools.RasterToPolygontRTP = newESRI.ArcGIS.ConversionTools.RasterToPolygon();  
            tRTP.in_raster = pWorkSpacePath + "Out3.img";  
            tRTP.out_polygon_features = pWorkSpacePath + "a2_shp.shp";  
            tRTP.raster_field = "VALUE";  
            tRTP.simplify = "SIMPLIFY";  
            //ScrollToBottom("RasterToPolygon");  
            tGeoResult = (IGeoProcessorResult)tGp.Execute(tRTP, null);  
            if (tGeoResult.Status == ESRI.ArcGIS.esriSystem.esriJobStatus.esriJobSucceeded)  
            {  
                licenseStatus = m_AoInitialize.CheckInExtension(esriLicenseExtensionCode.esriLicenseExtensionCodeSpatialAnalyst);  
                {  
                    ESRI.ArcGIS.DataManagementTools.SmoothPolygontSP = newESRI.ArcGIS.DataManagementTools.SmoothPolygon();  
                    tSP.algorithm = "PAEK";  
                    tSP.in_features = pWorkSpacePath + "a2_shp.shp";  
                    tSP.out_feature_class = pWorkSpacePath + "a2_shp_smooth.shp";  
                    tSP.tolerance = "2";  
                    tSP.endpoint_option = "FIXED_ENDPOINT";  
                    tSP.error_option = "FLAG_ERRORS";  
                    ScrollToBottom("SmoothPolygon");  
                    tGeoResult = (IGeoProcessorResult)tGp.Execute(tSP, null);  
                }  
            }  
        }  
    }  
    // arcengine liencesd 当使用高版本的liencesd在使用中,如果需要低版本的liencesd,需要动态签出liencesd,不然将不会执行成功!  
   // writeMessage(tGp);  
}  
  
MessageBox.Show("ok");  

 

 

         Geoprocessing是ArcGIS提供的一个非常实用的工具,借由Geoprocessing工具可以方便的调用ArcToolBox中提供的各类工具,本文在ArcEngine9.2平台环境下总结了调用ArcToolBox工具的使用方法:

        1、调用ArcToolBox工具方法

         以ArcToolBox->Analysis Tools->Proximity->Buffer工具的调用为例,C#代码如下:

         using ESRI.ArcGIS.AnalysisTools;         //添加引用
         using ESRI.ArcGIS.Geoprocessor;

         Geoprocessor gp = new Geoprocessor();    //初始化Geoprocessor
         gp.OverwriteOutput = true;                     //允许运算结果覆盖现有文件

         ESRI.ArcGIS.AnalysisTools.Buffer pBuffer = new ESRI.ArcGIS.AnalysisTools.Buffer(); //定义Buffer工具
         pBuffer.in_features = pVorLineLayer;    //输入对象,既可是IFeatureLayer对象,也可是完整文件路径如“D:\\data.shp”
         pBuffer.out_feature_class = pBuffer;     //输出对象,一般是包含输出文件名的完整文件路径,如“D:\\buffer.shp”

         //设置缓冲区的大小,即可是带单位的具体数值,如0.1 Decimal Degrees;也可是输入图层中的某个字段,如“BufferLeng”
         pBuffer.buffer_distance_or_field = "BufferLeng";    
         pBuffer.dissolve_option = "ALL";     //支持融合缓冲区重叠交叉部分
         gp.Execute(pBuffer, null);                //执行缓冲区分析

         ArcToolBox中各个工具调用时需添加的引用分别如下图所示:

                     

         参考网页:http://edndoc.esri.com/arcobjects/9.2/NET/c4ff8b68-0410-435f-b8e5-682d5cea47cf.htm

        2、参数设置

        在调用ArcToolBox执行具体的分析操作时,需要设置各类输入输出参数,简单概括起来说主要分为两类:对应于Environment Settings对话框的Geoprocessor对象设置、对应于具体操作窗口的方法设置。以ArcToolBox->Analysis Tools->Overlay->Intersect为例,C#代码如下:

            Geoprocessor gp = new Geoprocessor();
            gp.OverwriteOutput = true;    //覆盖原有文件并重写

            //Environment Settings对话框参数设置,具体名称参考操作界面Help中对应参数文档 

            object obj = gp.GetEnvironmentValue("Extent");  //设置Exten,大小写无关;          

            gp.SetEnvironmentValue("Extent", "MAXOF");     //或者"113.697050 115.074770 29.969986 31.362495"

            obj = gp.GetEnvironmentValue("OutputZFlag");                    //设置Output has Z Values
            gp.SetEnvironmentValue("OutputZFlag", "DEFAULT");

            obj = gp.GetEnvironmentValue("OutputMFlag");                    //设置Output has M Values
            gp.SetEnvironmentValue("OutputMFlag", "DEFAULT");

            obj = gp.GetEnvironmentValue("OutputCoordinateSystem");  //设置Output Coordinate System
            gp.SetEnvironmentValue("OutputCoordinateSystem", Application.StartupPath + "\\zhouyang.prj");

            obj = gp.GetEnvironmentValue("QualifiedFieldNames");         //设置Maintain fully qualifid field names
            gp.SetEnvironmentValue("QualifiedFieldNames", "QUALIFIED");

            //关于Environment Settings的设置可以参考ArcMap操作界面提供的文档,如图所示:

                      

                                                            Environment Settings设置

                      

                                        参数名称参考帮助,如上图所示为Extent,其取值有三种形式

            //具体操作窗口的方法设置

            Intersect pIntersect = new Intersect();
            //多个对象的输入:用分号隔开包含完整路径的文件名
            pIntersect.in_features = pInputFeature1 + ";" + pInputFeature2;

     //多个对象的输入:使用IGpValueTableObject接口,该接口可以设置Rank(http://resources.esri.com/help/9.3/arcgisengine/dotnet/84349562-e062-44ee-8db0-9fcdcd64708b.htm

     //object inputfeature1 = @"D:\周杨\贝贝\WuhanCity\ThiessenPolygons_Line_Buffer.shp";
            //object inputfeature2 = @"D:\周杨\贝贝\wuhanCity_shp\poi Point.shp";
            //IGpValueTableObject pObject = new GpValueTableObjectClass();
            //pObject.SetColumns(2);
            //pObject.AddRow(ref inputfeature1);
            //pObject.AddRow(ref inputfeature2);
            //pIntersect.in_features = pObject;
            pIntersect.out_feature_class = pOutputFeature;
            pIntersect.join_attributes = "All";
            pIntersect.output_type = "POINT";

     gp.Execute(pIntersect, null);      //执行

     //Intersect参数设置跟弹出的Intersect对话框对应,如图所示:

    

    参考网页:http://edndoc.esri.com/arcobjects/9.2/NET/552ca115-f23b-4a74-a2c5-069c50d6cdcf.htm

  3、运行结果对象提取

  Geoprocessor对象通过Execute方法执行后将结果保存到指定输出路径下,通过也可以通过IGeoProcessorResult接口读取存储在内容中的结果对象,C#代码如下:

    //执行图层求交运算
            IGeoProcessorResult pResult = (IGeoProcessorResult)gp.Execute(pIntersect, null);

     IGPUtilities pGPUtil = new GPUtilitiesClass();
            IFeatureClass pFC;
            IQueryFilter pQF;
            pGPUtil.DecodeFeatureLayer(pResult.GetOutput(0),out pFC,out pQF);
            int count = pFC.FeatureCount(null);      //统计Feature对象个数
            IFeatureCursor pCursor = pFC.Insert(true);   //提取FeatureCursor对象
            IFeatureLayer pFeatureLayer = new FeatureLayerClass();
            pFeatureLayer.FeatureClass = pFC;
            m_mapControl.Map.AddLayer(pFeatureLayer);   //加载图层对象

参考网页:http://edndoc.esri.com/arcobjects/9.2/NET/1b14f488-84de-4e7f-8009-cfe612f8dcbe.htm

     其实总的说来,ESRI的官方帮助和各类在线帮助文档中都提供了相应的说明,可以很容易搞清楚一些内容,但是在具体的操作过程中,有时候经常得不到结果,这时候就需要关注下Environment Settings中的部分参数是否设置了,有可能没有像软件操作界面中那样进行默认设置。

 


鲜花

握手

雷人

路过

鸡蛋
该文章已有0人参与评论

请发表评论

全部评论

专题导读
热门推荐
热门话题
阅读排行榜

扫描微信二维码

查看手机版网站

随时了解更新最新资讯

139-2527-9053

在线客服(服务时间 9:00~18:00)

在线QQ客服
地址:深圳市南山区西丽大学城创智工业园
电邮:jeky_zhao#qq.com
移动电话:139-2527-9053

Powered by 互联科技 X3.4© 2001-2213 极客世界.|Sitemap