Latest Contributions by BL

Hi I'm using ArcGIS 10.2 Advanced/32-bit Python 2.7.3. I'm needing to create multiple graphs using Arcpy (4 feature classes per graph: 2 vertical line, 2 scatter plot). My code works great using a template for 2 vertical lines. However, I modified my python code and the graph template to add 2 additional xy/scatterplot series and now the template is not being followed. I know for sure my data is being input into Graph()/MakeGraph correctly. The only thing I can think of why the template is not following my design for 4 feature classes is that I am adding series out of order (in python) in relation to how the graph series are originally ordered/read from the .tee. I would love to just manually edit a few parameters in the .tee file to fix the issue. Anyone know where to find a list of parameters for these .tee graph templates or ESRI's graph() class? I prefer not to switch to matplotlib or Grapher over this matter. Thanks for any help you may offer. For example: Current order of input in template.tee: object Series3: TPointSeries ... end object Series4: TPointSeries ... end object Series1: TLineSeries ... end object Series2: TLineSeries ... end The way input is currently being input into arcpy.Graph() using python's 'zip' function: for yr0_points, yr1_points, yr0_datums, yr1_datums in zip(yr0_points_list, yr1_points_list, yr0_datums_list, yr1_datums_list): ... # Create the graph graph = arcpy.Graph() # Add a vertical bar series to the graph graph.addSeriesLineVertical(yr0_lyr, "Elevation", "NEAR_DIST", "") #relates to object Series1: TLineSeries graph.addSeriesLineVertical(yr1_lyr, "Elevation", "NEAR_DIST", "")  #relates to object Series2: TLineSeries graph.addSeriesScatterPlot(yr0_datum_lyr, "Elevation","", "Descriptio") #relates to object Series3: TPointSeries graph.addSeriesScatterPlot(yr1_datum_lyr, "Elevation","", "Descriptio") #relates to object Series4: TPointSeries ... arcpy.MakeGraph_management(graph_grf, graph, "temp_graph") arcpy.SaveGraph_management("temp_graph", output_graph_dir, "MAINTAIN_ASPECT_RATIO", 1209, 815) ... View more

I have a polygon feature class in which features are symbolized by unique values (region code). Now I want to create colored map labels for my features using the same colors from the symbology tab. I'm using the label method 'define different classes/label features differently' (Label Field: region code). Is there a way to import/assign the same RGB values of the features (via symbology tab) to the text symbol for each label class? Right now, the text symbols for all my text classes in the Labels tab are defaulted to black and I have to go and manually change RGB values for each class which can be a little time consuming. Btw, I'm using ArcGIS 10 SP5 (Arcinfo). Thank you. ... View more

I have 60 DEMs in a particular county created from aerial LiDAR that I'm trying to create watersheds for. One of the preliminary things I need to do is fill sinks, but first I need to know the fill limit for sinks in each DEM. To solve this problem, I read ESRI's section on creating depressionless DEMs and how to find sink depth. The problem is ESRI left me hanging as to which statistic is considerd the z limit.(http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//009z0000005q000000.htm). Below is the sinkdepth operation output for the first DEM. Which one of the following sinkdepth output statistics is the z limit value for Fill tool? <HistMin>1.1444091796875e-005</HistMin> <HistMax>0.2707290649414063</HistMax>  <Approximate>0</Approximate> <Metadata> <MDI key="STATISTICS_MINIMUM">1.1444091796875e-005</MDI> <MDI key="STATISTICS_MAXIMUM">0.27072906494141</MDI> <MDI key="STATISTICS_MEAN">0.03058447859501</MDI> <MDI key="STATISTICS_STDDEV">0.020851988436709</MDI> The code below is what i used to create the sink depth raster including .xml statistics of each DEM. It takes all unfilled DEMs in a folder and outputs DEMs for sinkdepth: # --------------------------------------------------------------------------- # find sink depth.py # Created on: 2012-04-27 00:02:16.00000 # Author: T. Lupher # Description: # --------------------------------------------------------------------------- # Import arcpy module import arcpy, sys, os from arcpy import env from arcpy.sa import * # Set workspace env.workspace = "D:/DEM/LIDAR/" # Allow overwrites arcpy.env.overwriteOutput = True # Check out any necessary licenses arcpy.CheckOutExtension("spatial") # Create list of rasters rasters = arcpy.ListRasters("*", "TIF") rasterCount = len(rasters)     counter = 0 i = 1 try:     while counter < rasterCount:         for item in rasters:                         # input/output directories             inputPath = "D:/DEM/"             flowdirFolder = "flowdirection/"             sinkFolder = "sinks/"             #pourPoints = "D:/DEM/pourpoints/alltraplocations.shp"             contribFolder = "sink_contributing_areas/"             zonestatFolder = "zonalstatistics/"             zonalfillFolder = "zonalfill/"             sinkdepthFolder = "sinkdepth/"             # input/output variables             outDropRaster = ""             flow = inputPath+flowdirFolder+"flowdir_"+str(i)+".tif"             sink = inputPath+sinkFolder+"sink_"+str(i)+".tif"             contrib = inputPath+contribFolder+"conarea_"+str(i)+".tif"             zstat = inputPath+zonestatFolder+"zonestat_"+str(i)+".tif"             zfill = inputPath+zonalfillFolder+"zonefill_"+str(i)+".tif"             sinkdepth = inputPath+sinkdepthFolder+"sinkdepth_"+str(i)+".tif"                                 # Process: Flow Direction             outFlow = FlowDirection(item, "NORMAL")             outFlow.save(flow)             # Process: Sink             outSink = Sink(flow)             outSink.save(sink)                         # Process: Watershed             outWatershed = Watershed(flow, sink, "VALUE")             outWatershed.save(contrib)             # Process: Zonal Statistics             outZstat = ZonalStatistics(contrib, "VALUE", item, "MINIMUM", "DATA")             outZstat.save(zstat)             # Process: Zonal Fill             outZfill = ZonalFill(contrib, item)             outZstat.save(zfill)             # Process: Minus             outSinkDepth = Minus(outZfill, outZstat)             outSinkDepth.save(sinkdepth)                         counter += 1             i += 1 except:     # If an error occurred while running a tool, then print the messages.     print arcpy.GetMessages()  ... View more