ArcGISEngine实现栅格数据裁剪

ArcGIS空间分析扩展模块提供了对栅格数据进行裁剪抽取操作的功能，在ArcGISEngine中可以调用相关的GP工具来对栅格数据操作。

•  调用空间分析ArcToolbox中Extract     By Polygon工具

privatevoid ExtractByPolygon()

        {

            IVariantArray parameters = newVarArrayClass();

            string sDocPath = @"C:\filerastergdb.gdb\raster";      //地图文档路径

            string sDocPath1 = @"C:\filerasteroutput.gdb\outraster";

           parameters.Add(sDocPath);

           

            string ss1 = "'545160.9444639276.889';'545216.242 4639266.57';'545227.355 4639308.375'";;

           parameters.Add(ss1);

           parameters.Add(sDocPath1);

           parameters.Add("INSIDE");

            IGeoProcessor2 gp = newGeoProcessorClass();

           gp.OverwriteOutput = true;

            IGeoProcessorResult pResult = null;

            try

            {

                pResult =gp.Execute("ExtractByPolygon", parameters, null);

                IWorkspaceFactory pFWSF = new ESRI.ArcGIS.DataSourcesGDB.FileGDBWorkspaceFactoryClass();

               ESRI.ArcGIS.Geodatabase.IWorkspace pWS = pFWSF.OpenFromFile(@"C:\filerasteroutput.gdb", 0) ;

               ESRI.ArcGIS.Geodatabase.IRasterWorkspaceEx pRasterWs = pWS asESRI.ArcGIS.Geodatabase.IRasterWorkspaceEx;

                IRasterDataset pRasterDataset =pRasterWs.OpenRasterDataset("outraster");

                IRasterLayer pRasterLayer = newRasterLayerClass();

               pRasterLayer.CreateFromDataset(pRasterDataset);

                axMapControl1.AddLayer(pRasterLayer);

                string mess = null;

                for (int k = 0; k < gp.MessageCount; k++)

                {

                    mess+= gp.GetMessage(k) + "\n";

                }

                MessageBox.Show(mess);

            }

            catch (Exception err)

            {

                string msg = err.Message;

                string mess = null;

                for (int k = 0; k < gp.MessageCount; k++)

                {

                    mess+= gp.GetMessage(k) + "\n";

                }

                MessageBox.Show(mess);

            }

        }

 

• 调用IExtractionOp接口对栅格数据裁剪

        publicvoid ExtractionRasterData(IGeoDataset pGeoDataset, ESRI.ArcGIS.Geometry.IPolygon pPolygone, bool binside)

        {

           ESRI.ArcGIS.SpatialAnalyst.IExtractionOp pExtractionOp = new ESRI.ArcGIS.SpatialAnalyst.RasterExtractionOpClass();

           ESRI.ArcGIS.Geodatabase.IGeoDataset pGeoOutput = pExtractionOp.Polygon(pGeoDataset, pPolygone,binside);

        }