Latest Contributions by PeterWilson

Hi Dan Unfortunately not yet, I started working through your code and got stuck in understanding your reasoning behind reshaping the array from two dimensional to three dimensional lines 31 to 36. I still have quiet a bit to learn about NumPy and haven't needed to work with anything more than a two dimensional array. ... View more

Hi Xander Thanks for the solution, this will definitely work and the approach is clean and simple, not over complicated. ... View more

Hi Xander I'm looking to mosaic the smaller areas, as you mentioned the following wouldn't be a problem if I mosaiced the tiles into a single image. The group mosaic's are a lot smaller in file such and easier to load within the AutoCAD software. ... View more

Hi Dan Thanks for the advice, I'll have look graph groupings let you know what I finally use to solve the following. ... View more

Hi Joshua Thanks for the following, sorry I should have been more clear in my post. I'm trying to re-write an existing Python module based on a previous post Extract Polyline Geometry Z and M - Equal Area Slope where I'd like to replace the use of ArcPy "SplitLine_Management" with NumPy arrays and Pandas DataFrames. I wan't to remove the need for and Advanced Licence for the end user. ... View more

Hi Dan I actually need the distance between each vertices. I'm trying to re-write the Python module to calculate the "Equal Area Slope"  using NumPy and Pandas and not ArcPy "SplitLine_Management" based on my previous post: Extract Polyline Geometry Z and M - Equal Area Slope. ... View more

Hi Dan, Thanks this will get me started. I'll post shortly what I come up with and if I get stuck. ... View more

Hi Joshua I've attached a copy of the polylines unsplit as well as a copy of the polylines split with the output fields and values that I'm trying to achieve. ... View more

Hi Dan Thanks for the following. In order to get the following solved before my deadline. I extracted the Z values from the DTM using the InterpolateShape_3d. I then split the polylines using SplitLine_management to split the polyline in segments. I added the the following fields: ["STElev", "EndElev", "STElevation", "EndElevation"]. I then calculated the Start Elevation and End Elevation using Calculate Geometry. Thereafter  I calculated the "STElevation" by subtracting each elevation by the initial STElev (i.e. 355.9268 - 355.8495), "EndElevation" (i.e. 355.6268 - 355.9268). I then used a Search Cursor and calculated the area for each trapezoid using the following formula: ''' Created on 22 Jun 2017 @author: PeterW ''' # import modules and site packages import arcpy # set environment settings arcpy.env.overwriteOutput = True def equal_area_slope(lfp_3d, lfp_3d_split): """Determine the Equal Area Slope for each Drainage Area's Longest Flowpath""" slope_area_dict = {} with arcpy.da.SearchCursor(lfp_3d_split, ["DrainID", "STElevation", "EndElevation", "SHAPE@LENGTH"]) as scur: # @UndefinedVariable for row in scur: trap_area = ((row[1]+row[2])/2)*row[3] slope_area_dict.setdefault(row[0], []).append(trap_area) with arcpy.da.UpdateCursor(lfp_3d, ["DrainID", "EASlope", "SHAPE@LENGTH"]) as upcur: # @UndefinedVariable for row in upcur: if row[0] in slope_area_dict: slope_area = sum(slope_area_dict[row[0]]) right_ordinate = (slope_area*2)/row[2] ea_slope = (right_ordinate/row[2])*100 row[1] = ea_slope upcur.updateRow(row) if __name__ == "__main__": lfp_3d = r"E:\Projects\2016\G113665\geoHydro\Model02\Results.gdb\LongestFlowPath_3D" lfp_3d_split = r"E:\Projects\2016\G113665\geoHydro\Model02\Results.gdb\LongestFlowPath_3D_split" equal_area_slope(lfp_3d, lfp_3d_split)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ I've started re-writing the following using NumPy arrays. I've converted the 3D polyline into a NumPy array. The array contains the "DrainID" X Coordinate, Y Coordinate, Z Coordinate for each point. How do I go about replicating the following using NumPy arrays? I trying to figure out the most Pythonic way of creating an array of "DrainID", "STElevation", "EndElevation" and "SegmentLength" Current Code: '' Created on 01 Jul 2017 @author: PeterW ''' # import modules and site-packages import numpy as np import pandas as pd from pathlib import Path import arcpy # set environment settings arcpy.env.overwriteOutput = True # check out extensions arcpy.CheckOutExtension("3D") def longestflowpath_3d(dtm, lfp_2d): """Generate LongestFlowPath 3D from LongestFlowPath 2D and DTM""" output_gdb = str(Path(lfp_2d).parents[0]) lfp_3d = "{0}\\{1}".format(output_gdb, "lfp_3d") lfp_3d = arcpy.InterpolateShape_3d(dtm, lfp_2d, lfp_3d, vertices_only=True) arcpy.FlipLine_edit(lfp_3d) return lfp_3d def area_under_profile(lfp_3d): """Determine Area Under Profile""" arr = arcpy.da.FeatureClassToNumPyArray(lfp_3d, ["DrainID", "SHAPE@X", "SHAPE@Y", "SHAPE@Z"], explode_to_points=True) # @UndefinedVariable print(arr) def equal_area_slope(dtm, lfp_2d): """Determine the Equal Area Slope for each Drainage Area's Longest Flowpath""" lfp_3d = longestflowpath_3d(dtm, lfp_2d) area_under_profile(lfp_3d) if __name__ == "__main__": dtm = r"E:\Projects\2016\G113665\geoHydro\DTM2\raw1" lfp_2d = r"E:\Projects\2016\G113665\geoHydro\Model02\Results.gdb\LongestFlowPath_2D_Orig" equal_area_slope(dtm, lfp_2d)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Hi Dan,  Thanks for the following, worked like a charm. MST - Minimum Spanning Tree Peano Curve - Space-Filling Curve ... View more

Hi Curtis The following has worked perfectly, with a few minor manual adjustments to the clusters. Thanks so much for the following. Same group for some reason, not sure why. ... View more