ArcGIS Python Mapping Module

Arc Python to interact with Arc-map

What is ArcPy Mapping and Why is it actually used for :

Arcpy.mapping was built for the professional GIS analyst (as well as for developers). Traditionally, the scenarios listed above had to be done using ArcObjects and it often proved to be a very difficult programming environment to learn for the average GIS professional. Arcpy.mapping is a courser-grained object model, meaning that the functions are designed in a way that a single arcpy.mapping function can replace many lines of ArcObjects code. The following is a very simple example of how arcpy.mapping can be used to reference an existing map document and export it to a PDF document with only two lines of code.

To explore the Overall Arcpy mapping Module use the following PDF:

Arcpy Python Mapping Module

Code for Getting the reference to the MXD :

YourMXD= arcpy.mapping.MapDocument("path to the MXD")

To export the above MXD to PDF using arcpy:

arcpy.mapping.ExportToPdf(YourMXD,"your pdf file")

Pointer to a MXD present on the Disk

Layout View of the Map Document

Arcpy.mapping is not a replacement for ArcObjects but rather an alternative for the different scenarios it supports. ArcObjects is still necessary for finer-grain development and application customization, whereas arcpy.mapping is intended for automating the contents of existing map documents and layer files.

Creating Geodatabases using Arcpy or Model builder

Creating Geo-databases in Arc Map

What is a Geodatabase:

• A geodatabase is the common data storage and management framework for all the ArcGIS Products. It combines Geo (spatial data) with Database (data repository) to create a central data repository for all the spatial data storage and management in your ArcGIS platform.
• It can be leveraged in desktop, server, or mobile environments and allows you to store GIS data in a central location for easy access and management.

The geodatabase offers you the ability to

• Store a rich collection of spatial data in a centralized location.
• Apply sophisticated rules and relationships to the data.
• Define advanced geospatial relational models (e.g., topologies, networks).
• Maintain integrity of spatial data with a consistent, accurate database.
• Work within a multiuser access and editing environment.
• Integrate spatial data with other IT databases.
• Easily scale your storage solution.
• Support custom features and behavior.
• Leverage your spatial data to its full potential.

Arc-map geodatabase can be created using various way:

1. Manually creating through Arc-Catalog
2. Using the Model workflow to create it Programmatic-ally
3. Using the Python integrated to Create Geodatabase Programmatic

1. Model for Creating Feature Class

Model to create Geodatabase in ArcGIS

Set the workspace and OutputFC as the Model parameter

AddField will add a column to the newly created Feature class and can be populated with the user defined values

The Search, Insert or Update cursor can be used to manipulate the values present in the new Feature class

For achieving the above functionality it is advised to write a python code using the arcpy library.

To know more detailed step by step information about the creating a feature class in the Geodatabase
Click here - Create Geodatabase Feature Class - PDF

Exporting the feature class to CSV Format

Exporting Geodatabase toComma Separated Value Format

Let start creating an python addin for converting the feature class in your geodatabase to CSV format.So that everyone can understand your feature class and its values.

• If you are very new to python addin then go through- How to start with python addin 
• The input feature class and output csv file path are hard coded in the sample program
• For dynamic conversion get the path of the files through a variable and pass it as a parameter to the defined functions

CSV Fromat

def ex_fc_to_csv(fc, out_csv):

    """

    Export all vertex in a line or poly feature class to a csv with OID

   

    example

        import geom_snippets

        geom_snippets.ex_fc_to_csv(r"c:\proj\fc1.shp", r"c:\proj\fc1.csv")

    output csv looks like this

        1,56019.99998067904,69118.00001450378

        1,56159.99998080942,69026.0000144181

        1,56359.999980995686,68913.00001431286

        2,34985.00002508866,68936.00001433428

        2,35178.000025268404,68805.00001421227

    """

    import csv

    import json

    with open(out_csv, 'w') as csvfile:

        csvwriter = csv.writer(csvfile, delimiter=',', lineterminator='\n')

        with arcpy.da.SearchCursor(fc,

                                   field_names=("OID@","SHAPE@JSON")) as cursor:

            for row in cursor:

                geom = json.loads(row[1])

                for path in geom['paths']:

                    for pt in path:

                        csvwriter.writerow([row[0]] + pt)

How to check for condition in Geoprocessing Service or Model Builder

How to check for if/else conditions in Geoprocessing Service

TOOL USED:  

Arc python and Script tool

Python script will be deployed in the ARC Toolbox to create the Conversion Tool 

WORKING OF THE TOOL:

Script tool is created in the ArcMap and the required input & Output Parameter are set. The script tool source is set to the custom python code present in an Python file(.py)

Out motive is to use check for the count of a features in an feature class and performing an conditional Execution

PROCEDURE:

• Open ArcMap 10.0+ 
• Create your own Geometric Network or Open your Existing Geometric Network
• Add your Geometric network to the Table of contents in the ArcMap as Shown Below

·         Create a new Model in your myToolBox .Right Click on MyToolBox > New > Model.Following window will appear on the Screen 

•    Add your desired feature class to the Model. (Right Click on any empty space and select ADD DATA  and select your feature class in the object browser)
• ·         Add Get Count and Calculate value Model from the System Toolbox Connect the model as shown below 

• ·         Now open calculate value Model and in the Expression Textbox call the function Calculate(%RowCount%)
• ·         Here %RowCount% will act as an Inline Variable and Compiler automatically replace the %RowCount% with the value of the rowcount field
• ·         At the Code block defined the function Calculate as

def calculate(row):

if row>0:

return true

else:

return false

• ·         Make the return Data type as Boolean
• ·         The output Boolean can be set as a precondition for any model. So that the Model  get executes only when the output of calculate value is true

How to create Graphical User Interface in Arcpy

To Show GUI in your Script tool of your ArcMap

Python consist of  huge number of GUI frameworks (or toolkits) available for it, from TkInter to a number of other cross-platform solutions.

In your case we going to use the default TkInter Module to access the windows form from the python script tool.

TkInter is being proposed as the most interactive and easy to program module in Python GUI framework

1. Create a new Script tool in the Arcmap's ArcCatalog window in your desired location

2. Copy the below codes in to a notepad file and save it as python .py extension

3.Go to the script tool and Right Click > Properties > Parameter
4.Browse to your newly created python code as the source file to execute

5.All is well.Start GUI python from TkInter. Happy Programming

from Tkinter import *

import datetime

from array import *

root = Tk()

logo = PhotoImage(file="Image.Gif")

w2 = Label(root,

           image=logo).grid(row=0,column=1)

explanation = """Team Geo Zoner"""

w2 = Label(root, 

          compound = CENTER,

           fg="Red",

           bg="black",

           font="Times 55 bold",

          text=explanation)

w2.grid(row=0,column=0)

p=str(datetime.datetime.now().isoformat())

texter="Your Heading Goes Here"+str(p);

w3 = Label(root, 

          compound = CENTER,

           fg="Green",

           font="Times 20 bold",

          text=texter)

w3.grid(row=1,column=0)

Parameters = ['Number of Inputs','Number of Ouputs','Number of Processing Parameters','Number of Intermediate data']

r = 2

i=1

for c in Parameters:

    Label(text=c, relief=RIDGE,width=25).grid(row=r,column=0)

    Entry(text="Value", relief=SUNKEN,width=10).grid(row=r,column=1)

    r = r + 1

    i = i + 1

root.title("GISSTUDY")

root.mainloop()

ArcObjects to add a flag to the Geometric Network .Net

Creating a Flag in ArcMap using ArcObjects .NET

1. Add a new add-in component to the CreatingOutputs solution and name the component

AddFlagTool. Select Tool as the type of add-in, set the confi guration of the tool as shown in

2. Enter the following using directives at the top of the AddFlagTool.cs fi le’s code window:

using ESRI.ArcGIS.ArcMapUI;

using ESRI.ArcGIS.Geometry;

using ESRI.ArcGIS.Display;

using ESRI.ArcGIS.Carto;

3. You are going to write the code for handling the OnMouseDown event. In other words, you want

your code to be executed when a user clicks somewhere in the main window of the software application. Because all the events and their handlers are defined in the base class of all tools(ESRI.ArcGIS.Desktop.AddIns.Tool), you have to override the necessary handlers. So write

the following code in the AddflagTool.cs file inside the definition for the AddflagTool class and outside of any method: protected override

As soon as you type these two keywords and press the spacebar, you will see the list of all available handlers.

Find OnMouseDown in the list as shown in  below Figure and press Enter.

4.On the Mouse down event call the add flag function in the another library

5.For the Add flag function add the following code in the other library addflag function

        public static void AddFlag(IPoint pPnt, IApplication app, double snapTol)

        {

            //IProgressDialogFactory pProDFact = null;

            //IStepProgressor pStepPro = null;

            //IProgressDialog2 pProDlg = null;

            //ITrackCancel pTrkCan = null;

            IGeometricNetwork gn = null;

            IPoint snappedPoint = null;

            IFlagDisplay pFlagDisplay = null;

            INetFlag startNetFlag = null;

            INetworkAnalysisExt pNetAnalysisExt = null;

            IMap pMap = null;

            UID pID = null;

            int EID = -1;

            try

            {

                pID = new UID();

                pID.Value = "esriEditorExt.UtilityNetworkAnalysisExt";

                pNetAnalysisExt = (INetworkAnalysisExt)app.FindExtensionByCLSID(pID);

                gn = pNetAnalysisExt.CurrentNetwork;

                pMap = (app.Document as IMxDocument).FocusMap;

                startNetFlag = Globals.GetJunctionFlagWithGN(ref pPnt, ref pMap, ref gn, snapTol,out snappedPoint, out EID, out  pFlagDisplay, true) as INetFlag;

                if (startNetFlag == null)

                {

                    //startNetFlag = Globals.GetEdgeFlag(ref pPnt, ref pMap, ref gnList, snapTol, ref gnIdx, out snappedPoint, out EID, out distanceAlong, out  pFlagDisplay, true) as INetFlag;

                }

                if (app != null)

                {

                    

                    Globals.AddFlagToGN(ref pNetAnalysisExt, ref  gn, ref pFlagDisplay);

                    //  pFlagDisplay

                    pNetAnalysisExt = null;

                    pID = null;

                }

            }

            catch (Exception)

            {

                

                throw;

            }

        }

5. The above function will make a call to the getJunctionflag code.create a function block called "getjunctionflag" with the class library Global  and place the below code:

using System;

using System.Collections.Generic;

using System.Linq;

using System.Text;

using ESRI.ArcGIS.SystemUI;

using ESRI.ArcGIS.Geometry;

using ESRI.ArcGIS.GeoDatabaseUI;

using ESRI.ArcGIS.Geodatabase;

using ESRI.ArcGIS.Framework;

using ESRI.ArcGIS.esriSystem;

using ESRI.ArcGIS.Display;

using ESRI.ArcGIS.DataSourcesGDB;

using ESRI.ArcGIS.DataSourcesFile;

using ESRI.ArcGIS.CartoUI;

using ESRI.ArcGIS.Carto;

using ESRI.ArcGIS.ArcMapUI;

using ESRI.ArcGIS.ArcMap;

using ESRI.ArcGIS.NetworkAnalysis;

using ESRI.ArcGIS.NetworkAnalyst;

using ESRI.ArcGIS.EditorExt;

namespace Toolkit

{

    class Globals

    {

        public enum GNTypes

        {

            Flags = 1, Barries = 2, Results = 3

        };

        public enum flagType { EdgeFlag, JunctionFlag, EdgeBarrier, JunctionBarrier };

    

        public static IJunctionFlag GetJunctionFlagWithGN(ref IPoint point, ref IMap map, ref IGeometricNetwork gn, double snapTol, out IPoint snappedPoint, out int EID, out  IFlagDisplay pFlagDisplay, bool Flag)

        {

            //Initialize output variables

            snappedPoint = null;

            EID = -1;

            pFlagDisplay = null;

            int FCID = -1, FID = -1, subID = -1;

            IGeoDataset pDS = null;

            IPointToEID pointToEID = null;

            INetElements netElements = null;

            INetFlag junctionFlag = null;

            try

            {

                pFlagDisplay = null;

                pDS = gn.FeatureDataset as IGeoDataset;

                point.Project(pDS.SpatialReference);

                pointToEID = new PointToEIDClass() as IPointToEID;

                // find the nearest junction element to this Point

                pointToEID.GeometricNetwork = gn as IGeometricNetwork;

                pointToEID.SourceMap = map;

                pointToEID.SnapTolerance = snapTol;

                try

                {

                    pointToEID.GetNearestJunction(point, out EID, out snappedPoint);

                }

                catch (Exception ex)

                {

                }

                if (snappedPoint == null)

                    return null;

                // convert the EID to a feature class ID, feature ID, and sub ID

                netElements = gn.Network as INetElements;

                try

                {

                    netElements.QueryIDs(EID, esriElementType.esriETJunction, out FCID, out FID, out subID);

                }

                catch (Exception ex)

                {

                    return null;

                }

                //Create flag for start of trace

                junctionFlag = new JunctionFlagClass() as INetFlag;

                junctionFlag.UserClassID = FCID;

                junctionFlag.UserID = FID;

                junctionFlag.UserSubID = subID;

                if (junctionFlag is IEdgeFlag)

                {

                    pFlagDisplay = new EdgeFlagDisplayClass();

                    if (Flag)

                        pFlagDisplay.Symbol = CreateNetworkFlagBarrierSymbol(flagType.EdgeFlag) as ISymbol;

                    else

                        pFlagDisplay.Symbol = CreateNetworkFlagBarrierSymbol(flagType.EdgeBarrier) as ISymbol;

                }

                else

                {

                    pFlagDisplay = new JunctionFlagDisplayClass();

                    if (Flag)

                        pFlagDisplay.Symbol = CreateNetworkFlagBarrierSymbol(flagType.JunctionFlag) as ISymbol;

                    else

                        pFlagDisplay.Symbol = CreateNetworkFlagBarrierSymbol(flagType.JunctionBarrier) as ISymbol;

                }

                pFlagDisplay.ClientClassID = FCID;

                pFlagDisplay.FeatureClassID = FID;

                pFlagDisplay.SubID = subID;

                pFlagDisplay.Geometry = snappedPoint;

                return junctionFlag as IJunctionFlag;

            }

            catch

            {

                return null;

            }

            finally

            {

                pDS = null;

                pointToEID = null;

                netElements = null;

            }

        }

        public static ISimpleMarkerSymbol CreateNetworkFlagBarrierSymbol(flagType flgType)

        {

            ISimpleMarkerSymbol pSymbolFlag = null;

            switch (flgType)

            {

                case flagType.EdgeFlag:

                    pSymbolFlag = new SimpleMarkerSymbolClass();

                    pSymbolFlag.Style = esriSimpleMarkerStyle.esriSMSSquare;

                    pSymbolFlag.Angle = 0;

                    //pSymbolFlag.Color = 

                    pSymbolFlag.Outline = true;

                    pSymbolFlag.OutlineSize = 1;

                    //pSymbolFlag.OutlineColor = GetColor(0, 0, 0);

                    pSymbolFlag.Size = 10; //TODO: UserConfig

                    break;

                case flagType.JunctionFlag:

                    pSymbolFlag = new SimpleMarkerSymbolClass();

                    pSymbolFlag.Style = esriSimpleMarkerStyle.esriSMSCircle;

                    pSymbolFlag.Angle = 0;

                    //pSymbolFlag.Color = GetColor(0, 255, 0);

                    pSymbolFlag.Outline = true;

                    pSymbolFlag.OutlineSize = 1;

                    //pSymbolFlag.OutlineColor = GetColor(0, 0, 0);

                    pSymbolFlag.Size = 10; //TODO: UserConfig

                    break;

                case flagType.EdgeBarrier:

                    pSymbolFlag = new SimpleMarkerSymbolClass();

                    pSymbolFlag.Style = esriSimpleMarkerStyle.esriSMSDiamond;

                    pSymbolFlag.Angle = 0;

                    //pSymbolFlag.Color = GetColor(255, 0, 0);

                    pSymbolFlag.Outline = true;

                    pSymbolFlag.OutlineSize = 1;

                    //pSymbolFlag.OutlineColor = GetColor(0, 0, 0);

                    pSymbolFlag.Size = 10; //TODO: UserConfig

                    break;

                case flagType.JunctionBarrier:

                    pSymbolFlag = new SimpleMarkerSymbolClass();

                    pSymbolFlag.Style = esriSimpleMarkerStyle.esriSMSX;

                    pSymbolFlag.Angle = 0;

                    //pSymbolFlag.Color = GetColor(255, 0, 0);

                    pSymbolFlag.Outline = true;

                    pSymbolFlag.OutlineSize = 1;

                    //pSymbolFlag.OutlineColor = GetColor(0, 0, 0);

                    pSymbolFlag.Size = 10; //TODO: UserConfig

                    break;

                default:

                    pSymbolFlag = new SimpleMarkerSymbolClass();

                    pSymbolFlag.Style = esriSimpleMarkerStyle.esriSMSCircle;

                    pSymbolFlag.Angle = 0;

                    //pSymbolFlag.Color = GetColor(0, 255, 0);

                    pSymbolFlag.Outline = true;

                    pSymbolFlag.OutlineSize = 1;

                    //pSymbolFlag.OutlineColor = GetColor(0, 0, 0);

                    pSymbolFlag.Size = 10; //TODO: UserConfig

                    break;

            }

            return pSymbolFlag;

        }

        public static void AddFlagToGN(ref INetworkAnalysisExt pNetworkAnalysisExt, ref ESRI.ArcGIS.Geodatabase.IGeometricNetwork pGeomNet, ref  IFlagDisplay pFlagDsiplay)

        {

            INetworkAnalysisExtFlags pNetworkAnalysisExtFlags = null;

            try

            {

                if (pNetworkAnalysisExt.CurrentNetwork != pGeomNet)

                {

                    pNetworkAnalysisExt.CurrentNetwork = pGeomNet;

                }

                pNetworkAnalysisExtFlags = (INetworkAnalysisExtFlags)pNetworkAnalysisExt;

                if (pFlagDsiplay is IEdgeFlagDisplay)

                {

                    pNetworkAnalysisExtFlags.AddEdgeFlag(pFlagDsiplay as IEdgeFlagDisplay);

                }

                else

                {

                    pNetworkAnalysisExtFlags.AddJunctionFlag(pFlagDsiplay as IJunctionFlagDisplay);

                }

            }

            catch

            {

            }

            finally

            {

                pNetworkAnalysisExtFlags = null;

            }

        }

    }

}

6. Run code and set the default debug application as ArcMap

7. Once the ArcMap is open go to customise code and in the command tab add the tool to your window

8.click on the tool will add an flag symbol to the map(use f5 to refresh the map)

Addin Setting flow direction in Geometric Network

Addin source code for setup flow direction   

1.Create an c# object with establish flow as a function

2.Call this function from the addin on click method by passing the application reference and the type of flow direction you wish to set to your network

To know more about geometric network setting flow direction setting flow direction

Setting Flow Direction through ArcObjects Code

public static void EstablishFlow(Globals.GNFlowDirection flowDirection, IApplication app)
        {
            IProgressDialog2 progressDialogs= default(IProgressDialog2);
            IProgressDialogFactory progressDialogFactorys = null;
            IEditor editing = null;
         
            IMouseCursor appCursor = null;
            INetworkAnalysisExt netExt = null;
            UID pUID = null;
         
            List<IGeometricNetwork> geonetworklist= null;
            ITrackCancel trackCancel = null;
            Int32 int32_hWnd;
            IStepProgressor stepProgressorobj= null;
            IMxDocument mxdoc = null;
   IEditLayers eLayers = null;
   IMap pMap = null;
            try
            {




                int calcCount = 0;

                //Get editing

                editing = Globals.geteditor(ref app);

                if (editing.EditState != esriEditState.esriStateEditing)
                {
                    MessageBox.Show(Teamgeofunctions.Localizer.GetString("MustBEditg"), Teamgeofunctions.Localizer.GetString("GeoNetToolsLbl_2"));
                    return;
                }

                eLayers = editing as IEditLayers;

                //Change mouse cursor to wait - automatically changes back (ArcGIS Desktop only)
                appCursor = new MouseCursorClass();
                appCursor.SetCursor(2);

                ESRI.ArcGIS.esriSystem.IStatusBar statusBar = app.StatusBar;
                statusBar.set_Message(0, Teamgeofunctions.Localizer.GetString("GeoNetToolsWait_1"));

                //Get NA Extension in order to update the current network with the first visible network
                pUID = new UIDClass();
                pUID.Value = "esriEditorExt.UtilityNetworkAnalysisExt";
                netExt = app.FindExtensionByCLSID(pUID) as INetworkAnalysisExt;

                //Get Visible geometric networks
                pMap = editing.Map;
                geonetworklist= Globals.GetGeometricNetworksCurrentlyVisible(ref pMap);


                if (gnList.Count > 0)
                {


                    //ProgressBar
                    progressDialogFactorys= new ProgressDialogFactoryClass();

                    // Create a CancelTracker
                    trackCancel = new CancelTrackerClass();

                    // Set the properties of the Step Progressor
                    int32_hWnd = app.hWnd;
                    stepProgressorobj= progressDialogFactory.Create(trackCancel, int32_hWnd);
                    stepProgressor.MinRange = 0;
                    stepProgressor.MaxRange = gnList.Count;
                    stepProgressor.StepValue = 1;
                    stepProgressor.Message = "";
                    stepProgressor.Hide();

                    // Create the ProgressDialog. This automatically displays the dialog
                    progressDialogs= (ESRI.ArcGIS.Framework.IProgressDialog2)stepProgressor; // Explict Cast

                    // Set the properties of the ProgressDialog
                    progressDialog.CancelEnabled = false;
                    progressDialog.Description = Teamgeofunctions.Localizer.GetString("GeoNetToolsDesc_1");
                    progressDialog.Title = Teamgeofunctions.Localizer.GetString("GeoNetToolsTitle_1");
                    progressDialog.Animation = ESRI.ArcGIS.Framework.esriProgressAnimationTypes.esriProgressSpiral;

                    bool editStarted = false;
                    try
                    {// Create an edit operation enabling undo/redo
                        editing.StartOperation();
                        editStarted = true;
                    }
                    catch
                    {
                        editStarted = false;
                    }
                    IGeometricNetwork gn = null;
                    IEnumFeatureClass enumFC = null;
                    INetwork net = null;
                    IUtilityNetworkGEN unet = null;
                    IEnumNetEID edgeEIDs = null;
                    //IFeatureLayer fLayer = null;
                    try
                    {

                        for (int i = 0; i < gnList.Count; i++)
                        {

                            gn = gnList[i] as IGeometricNetwork;
                         
                            stepProgressor.Message = Teamgeofunctions.Localizer.GetString("GeoNetToolsProc_1");// +fLayer.Name;
                            //Establish flow using AncillaryRole values


                            if (flowDirection == Globals.GNFlowDirection.AncillaryRole)
                            {
                                enumFC = gn.get_ClassesByNetworkAncillaryRole(esriNetworkClassAncillaryRole.esriNCARSourceSink);
                                if (enumFC.Next() == null)
                                    MessageBox.Show(Teamgeofunctions.Localizer.GetString("GeoNetToolsError_1a") + gn.FeatureDataset.Name + Teamgeofunctions.Localizer.GetString("GeoNetToolsError_1b") + Environment.NewLine +
                                                    Teamgeofunctions.Localizer.GetString("GeoNetToolsError_1c"), Teamgeofunctions.Localizer.GetString("GeoNetToolsProc_2"));
                                else
                                {
                                    gn.EstablishFlowDirection();
                                    calcCount += 1;
                                }
                            }

                            //Establish flow direction based on digitized direction.
                            else
                            {
                                net = gn.Network;
                                unet = net as IUtilityNetworkGEN;
                                edgeEIDs = net.CreateNetBrowser(esriElementType.esriETEdge);
                                edgeEIDs.Reset(); int edgeEID;
                                for (long j = 0; j < edgeEIDs.Count; j++)
                                {
                                    edgeEID = edgeEIDs.Next();
                                    unet.SetFlowDirection(edgeEID, esriFlowDirection.esriFDWithFlow);
                                }
                                calcCount += 1;
                            }
                            stepProgressor.Step();

                        }
                    }

                    catch (Exception ex)
                    {
                        editing.AbortOperation();
                        MessageBox.Show(Teamgeofunctions.Localizer.GetString("GeoNetToolsProc_2") + "\n" + ex.Message, ex.Source);
                    }
                    finally
                    {
                        if (enumFC != null)
                            Marshal.ReleaseComObject(enumFC);

                        gn = null;
                        enumFC = null;
                        net = null;
                        unet = null;
                        edgeEIDs = null;
                        //fLayer = null;
                    }
                    if (editStarted)
                    {   // Stop the edit operation
                        if (flowDirection == Globals.GNFlowDirection.AncillaryRole)
                            editing.StopOperation(Teamgeofunctions.Localizer.GetString("GeoNetToolsProc_2"));
                        else
                            editing.StopOperation(Teamgeofunctions.Localizer.GetString("GeoNetToolsProc_3"));
                    }
                    object Missing = Type.Missing;
                    mxdoc = app.Document as IMxDocument;
                    mxdoc.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGraphics, Missing, mxdoc.ActiveView.Extent);

           
            }
            catch (Exception ex)
            {
       
                return;


            }
            finally
            {
                // delete the references to cleanup the memory
                if (progressDialogs!= null)
                    progressDialog.HideDialog();
                progressDialogs= null;
                progressDialogFactorys= null;
                editing = null;
                eLayers = null;
                appCursor = null;
                netExt = null;
                pUID = null;
                pMap = null;
                geonetworklist= null;
                trackCancel = null;

                stepProgressorobj= null;
                mxdoc = null;
            }

Custom Geoprocessing Tool

Arcobjects Geoprocessing Procedure:

1. Create a geoprocessor object and set its properties.
2. Set the appropriate values for the tool’s parameters.
3. Execute a tool using the geoprocessor object’s Execute() method.

In order to create a geoprocessor object, the IGeoProcessor2 interface can be used. This interface
is defined in the geoprocessing library of ArcObjects and is implemented by the GeoProcessor
CoClass.

IGeoProcessor2 gp = new GeoProcessorClass();
//add the result of geoprocessing as a new layer to Map
gp.AddOutputsToMap = true;
//if output of geoprocessing exists before the execution of tool
//it will be overwritten
gp.OverwriteOutput = true;

Assume that you want to create Thiessen polygons for the cities FeatureClass to create proximal
zones for all cities. The proximal zones represent full areas where any location inside the zone is
closer to its associated city than any other city. Based on the Create Thiessen Polygons reference
page, you have to provide at least an input point FeatureLayer and the path to the output
FeatureClass.
In order to create and set each parameter, the IVariantArray interface of the System library must
be used.

IVariantArray parameters = new VarArrayClass();

Each parameter has to be added to the IVariantArray interface in the exact order that is specified
on the tool’s reference page.

//in_features
parameters.Add(@"D:\DataFolder\fileGDB.gdb\cities");
//out_feature_class
parameters.Add(@"D:\DataFolder\fileGDB.gdb\citiesThiessen");
//fields_to_copy(Optional)
parameters.Add("ALL");

As it is illustrated on the Create Thiessen Polygons reference page, the third parameter is optional.
You can simply not add any value to IVariantArray or, as shown in the preceding code, you can
provide an appropriate value for the optional parameter. You can skip the optional parameter
using an empty string as input to IVariantArray’s Add() method. For example, look at the
reference page of the Buffer tool; you can see that three of the seven available parameters are
mandatory. The following code demonstrates how to skip the fourth parameter and specify the
fifth parameter:

//1-in_features
parameters.Add(@"D:\test.gdb\cities");
//2-out_feature_class
parameters.Add(@"D:\test.gdb\citiesBuffer");
//3-buffer_distance_or_field
parameters.Add("50 kilometers");
//4-line_side(Optional)
parameters.Add("");
//5-line_end_type(Optional)
parameters.Add("ROUND");
//6 &7 there is no need to provide empty string
//for the rest of parameters since you don't want to set them
After setting all the required parameters, all you need to run a tool is to call the Execute() method
of the geoprocessor object. The Execute() method solicits the name of the tool and its parameters.

gp.Execute("CreateThiessenPolygons_analysis", parameters, null);
The following code shows the complete code for this example. In order to run the code, you need to
add references to the Geoprocessing and System libraries of ArcObjects:
IGeoProcessor2 gp = new GeoProcessorClass();
//add the result of geoprocessing as a new layer to Map
gp.AddOutputsToMap = true;
//if output of geoprocessing exists before the execution of tool
//it will be overwritten
gp.OverwriteOutput = true;
IVariantArray parameters = new VarArrayClass();
//in_features
parameters.Add(@"D:\DataFolder\fileGDB.gdb\cities");
//out_feature_class
parameters.Add(@"D:\DataFolder\fileGDB.gdb\citiesThiessen");
//fields_to_copy(Optional)
parameters.Add("ALL");
//or parameters.Add("");
gp.Execute("CreateThiessenPolygons_analysis", parameters, null);


The Geoprocessing library of ArcObjects is accessible through the ESRI.ArcGIS.Geoprocessing
namespace. This library contains a few hundred types which can be used to run and manage tools
and GIS workfl ows. You learned earlier in this section that IGeoProcessor2 is the main interface of
this library and the easiest way to run a geoprocessing tool is to call its Execute() method.
However, using IGeoProcessor2 is not the only approach to run a geoprocessing tool or model.
There are some managed assemblies created by Esri to performing geoprocessing in a managed
way. A managed way means there is a native .NET assembly (the Geoprocessor assembly) that is
a wrapper for some types in the Geoprocessing library of ArcObjects, and there are other .NET
assemblies for each system toolbox.
These native .NET assemblies provide an even easier way to run a system tool. In general, the
procedure for running a tool using the Geoprocessor-managed assembly is the same as running a tool