Extending ArcGIS Pro with.net and Python: Interactive Analytics. Carlos A. Osorio-Murillo Mark Janikas

Transcription

1 Extending ArcGIS Pro with.net and Python: Interactive Analytics Carlos A. Osorio-Murillo Mark Janikas

2 Introduction ArcGIS Pro is highly customizable. From an application perspective,.net can be used to construct buttons, ribbons and tool specific user interfaces. Through the use of Python and Geoprocessing, your desired functionality can run both in demand and interactively with the.net extension you created. This session introduces the key concepts necessary to get started with extending your ArcGIS Pro Application and provides a clear case study of interactive analytics to show you the power of the strategies outlined.

3 Objectives.Net & Python Integrate.Net and Python for improving workflows Powering a Density-Based Clustering method (OPTICS) through an Add-in Quick demo how to create an Add-in Visual Studio

4 Unsupervised learning methods Find representavives groups Which criteria depends of user Partition (Kmean) Density-Based Clustering Spatially constrain (SKATER) DBSCAN OPTICS HDBSCAN

5 Density-based Clustering Three clustering algorithms Fast Automatic Detection Interactive

6 OPTICS Process Produce a special order of the database with a density-base clustering structure God for both automatic and interactive cluster analysis, including finding intrinsic clustering structure Can be represented graphically (Plot) arcpy.densitybasedclustering_stats()

7 OPTICS Process arcpy.densitybasedclustering_stats() (Generate Reachability information) Cluster detection Map clustering (Applying color exclusion schema)

8 Python Tool for Detecting cluster def execute(self, parameters, messages): import SSCluster as SSC import SSDataObject as SSDO import numpy as NUM #### Get Parameters #### inputfc = parameters[0].valueastext minpoints = int(parameters[1].valueastext) tolerance = int(parameters[2].valueastext) threshold = int(parameters[3].valueastext) #### Allow overwrite Output ### ARCPY.env.overwriteOutput = True #### Load Data #### ssdo = SSDO.SSDataObject(inputFC) ssdo.obtaindata(fields = ["REACHORDER","REACHDIST"]) #### Sort Reachability Order ### ord = ssdo.fields["reachorder"].data #### Get Reachability Distances #### reachvalues = ssdo.fields["reachdist"].data #### Sorting by Index #### data = NUM.zeros((len(reachValues), 2), dtype= float) data[:, 0] = reachvalues data[:, 1] = ord odata = data[data[:,1].argsort()]

9 Distance Threshold Clusters Not enough minimum number of points to be considered OPTICS as an extension of DBSCAN #### Function to get DBSCAN Clusters #### def getdbscan(distances, threshold, minpoints, ordervalues): clusterarr = NUM.ones(len(distances), dtype = int)*-1 indices = distances <= threshold ini = end = 0 clusterindex = 1 c = 1 for i in NUM.arange(len(indices)-1): if indices[i] and indices[i] == indices[i+1]: if c == 1: ini = i c += 1 end = i + 1 else: if (end-ini) >= minpoints-1: clusterarr[ini:end+1] = clusterindex clusterindex += 1 ini = end = 0 c = 1 return clusterarr[ordervalues]

10 Cluster shape Arcgis Pro Python Module Merged clusters OPTICS Tolerance import SSCluster as SSC #### Initialize Detect Zone Class #### zones = SSC.DetectZones() #### Get Cluster By Tolerance ### idclusters, idcolor = zones.getclusters(reachvalues, ordervalues, minpoints, tolerance)

11 Interactive Clustering detection through Add-in and Python

12 Addin Demo Detecting clusters

13 Add-in Creation using Visual Studio Resources Build-Your-First-Add-in

14 Add-in Creation using Visual Studio

15 Creating an Add-in Demo

16 Design Control WPF ComboBox <UserControl x:class="interactiveanalytics.dockpane1view" xmlns=" xmlns:x=" xmlns:mc=" xmlns:d=" xmlns:ui="clr-namespace:interactiveanalytics" xmlns:extensions="clr-namespace:arcgis.desktop.extensions;assembly=arcgis.desktop.extensions" mc:ignorable="d" d:designheight="400" d:designwidth="300" d:datacontext="binding Path=ui.Dockpane1ViewModel}">.. TextBox Slider Dockpanel

17 Add-in Events Add-In UI App Python Geoprocessing Toolbox/Tool GetLayersInMap() Time line UpdateClusters() Update Simbology() Update Dataset Async method

18 GetLayersInMap() internal async void GetLayersInMap() if (MapView.Active == null) return; Obtain Point Layers with a specific field } var currentmap = MapView.Active.Map; var layers = currentmap.getlayersasflattenedlist(); List<FeatureLayer> featurelayers = new List<FeatureLayer>(); foreach (Layer lyr in layers) if (lyr is FeatureLayer) FeatureLayer flyr = lyr as FeatureLayer; bool check = true; check = (flyr.shapetype == esrigeometrytype.esrigeometrypoint); if (check) check = (flyr.shapetype == esrigeometrytype.esrigeometrypoint); /// Layer should contain Reachability distance field int hasclusteredopticlayer = await QueuedTask.Run<int>(() => return flyr.gettable().getdefinition().findfield("reachdist"); }); await QueuedTask.Run<int>(() => return flyr.gettable().getdefinition().findfield("reachdist"); }); if (hasclusteredopticlayer >= 0) featurelayers.add(flyr); } } } _layers = featurelayers; NotifyPropertyChanged(() => Layers);

19 UpdateClusters() internal async void UpdateCluster(bool istolerance) if (MapView.Active == null) return; // Get Layer Name string inputfc = SelectedLayer.Name; int minpoints = 2; bool parsed = Int32.TryParse(MinPoints, out minpoints); // Set PYT path string tool_path ="C:\\PATH\\UpdateCluster.pyt\\UpdateClusterTool"; IReadOnlyList<string> args= null; if (istolerance) /// Arguments for executing process using Tolerance args = Geoprocessing.MakeValueArray(inputFC, minpoints, ValueSlider, -1); } else /// Arguments for executing process using Threshold args = Geoprocessing.MakeValueArray(inputFC, minpoints, -1, ValueThreshold); } string tool_path ="C:\\PATH\\UpdateCluster.pyt\\UpdateClusterTool"; Task<IGPResult> task; task = Geoprocessing.ExecuteToolAsync(tool_path, args, flags: /// Execute the Tool in the python toolbox GPExecuteToolFlags.AddToHistory); task = Geoprocessing.ExecuteToolAsync(tool_path, args, flags: GPExecuteToolFlags.AddToHistory);

20 Update Simbology await QueuedTask.Run( () => /// Search for a specific Symbol /// Other styles Arcgis/Resouces/Styles/Styles.stylx SQLite DB SymbolStyleItem symbolstyleitem = (SymbolStyleItem)style.LookupItem(StyleItemType.PointSymbol, "Circle 1_Shapes_3"); pointsymbol = (CIMPointSymbol)symbolStyleItem.Symbol; /// Cluster Ids based in Color Schema int[] ids = new int[] -1, 1, 2, 3, 4, 5, 6, 7, 8 }; /// Set Colors string[] colors = new string[] "156,156,156", "166,206,227", "31,120,180", "178,223,138", "51,160,44", "251,154,153", "227,26,28", "253,191,111", "255,127,0" }; /// Color Field String[] fields = new string[] "COLOR_ID" }; SymbolStyleItem symbolstyleitem = /// Make a reference of the point symbol CIMSymbolReference symbolpointtemplate = pointsymbol.makesymbolreference(); (SymbolStyleItem)style.LookupItem(StyleItemType.PointSymbol, "Circle 1_Shapes_3"); /// Get definition of type symbology unique values UniqueValueRendererDefinition uniquevaluerendererdef = new UniqueValueRendererDefinition(fields, symbolpointtemplate, null, symbolpointtemplate, false); /// Get Current renderer of the Selected Layer CIMUniqueValueRenderer renderer = (CIMUniqueValueRenderer)SelectedLayer.CreateRenderer(uniqueValueRendererDef); CIMUniqueValueClass[] newclasses = new CIMUniqueValueClass[colors.Count()]; UniqueValueRendererDefinition uniquevaluerendererdef = /// Get Point Symbol as string for creating other point colors string point = pointsymbol.toxml(); new UniqueValueRendererDefinition(fields, symbolpointtemplate, null, symbolpointtemplate, false); /// Create Each Color for (int i = 0; i < ids.length; i++) CIMPointSymbol npoint = CIMPointSymbol.FromXml(point); if (ids[i] == -1) npoint.setsize(4); } else npoint.setsize(6); } CIMSymbolReference symbolpointtemplatei = npoint.makesymbolreference(); newclasses[i] = new CIMUniqueValueClass(); newclasses[i].values = new CIMUniqueValue[1]; newclasses[i].values[0] = new CIMUniqueValue(); newclasses[i].values[0].fieldvalues = new string[1]; newclasses[i].values[0].fieldvalues[0] = ids[i].tostring(); newclasses[i].label = ids[i].tostring(); newclasses[i].symbol = symbolpointtemplatei; var color = colors[i].split(','); double r = Convert.ToDouble(color[0]); double g = Convert.ToDouble(color[1]); double b = Convert.ToDouble(color[2]); newclasses[i].symbol.symbol.setcolor(cimcolor.creatergbcolor(r, g, b)); } /// Add Colors into the renderer renderer.groups[0].classes = newclasses; CIMUniqueValueRenderer renderer = (CIMUniqueValueRenderer)SelectedLayer.CreateRenderer(uniqueValueRendererDef); renderer.groups[0].classes = newclasses; SelectedLayer.SetRenderer(renderer); /// Apply new renderer in the layer SelectedLayer.SetRenderer(renderer); //SelectedLayer.RecalculateRenderer(); });

21 Conclusions - Add-in can improve Geoprocessing Tools (Density-Based Clustering) - Easy way to join Python Scripts and.net through an Add-in - ArcGIS Pro SDK Extension in Visual Studio great tool to initial Add-in

22 Thank you

23