Latest Contributions by XanderBakker

There are several ways to do this. The nice thing about having spatial information is that you can treat it as such. Since you are evaluating difference in height for the "same" location you could process parts of the file (selected by location). Now this will be a little less slow, however creating a spatial selection might be slow anyway when processing 13 million points (in a shapefile which is I believe slower than a file geodatabase that is properly indexed). Now what I wonder is this; the 13 million points sound as LiDAR data. If you have a LAS file you can create a LAS dataset and use the LAS Point Statistics as Raster tool (ArcGIS Help (10.2, 10.2.1, and 10.2.2) ) which has an option to extract the Z_RANGE. If you don't have the LAS file, you can create one using LAStools converting a text file to a LAS file. ... View more

You could change the expression to: def FindLabel([LOT] ,[SUBDIVISION]):   return "{0}\n{1}".format([LOT], [SUBDIVISION]) It suppresses the error. However, this will return None if there is no data: The nice part of using the python expression is that you can enhance the result, replacing the None for a different text def FindLabel([LOT] , [SUBDIVISION]):   return "{0}\n{1}".format([LOT] if [LOT] != None else "No LOT",  [SUBDIVISION] if [SUBDIVISION] != None else "No SUBDIVISION") ... will result in: ... or adding some formatting to the text: def FindLabel([LOT] , [SUBDIVISION]):   return "{0}\n{1}".format([LOT] if [LOT] != None else "<CLR red='255'>No LOT</CLR>",  [SUBDIVISION] if [SUBDIVISION] != None else "<CLR red='255'>No SUBDIVISION</CLR>") will give: ... View more

OK, just in case you're interested, but I'm sure Dan will come up with a better solution: import arcpy import os import math txt = r"D:\Xander\GeoNet\PointCloud\sample.txt" fc_pnt = r"D:\Xander\GeoNet\PointCloud\gdb\test.gdb\points05" fc_lines = r"D:\Xander\GeoNet\PointCloud\gdb\test.gdb\lines05" max_dist = 25 arcpy.env.overwriteOutput = True # sr = arcpy.SpatialReference() # create empty fc pnt_ws, pnt_name = os.path.split(fc_pnt) arcpy.CreateFeatureclass_management(pnt_ws, pnt_name, "POINT", "#", "DISABLED", "ENABLED") fld_X = "X" fld_Y = "Y" fld_Z = "Z" fld_A = "Att" fld_dist = "Dist" arcpy.AddField_management(fc_pnt, fld_X, "Double") arcpy.AddField_management(fc_pnt, fld_Y, "Double") arcpy.AddField_management(fc_pnt, fld_Z, "Double") arcpy.AddField_management(fc_pnt, fld_A, "Integer") arcpy.AddField_management(fc_pnt, fld_dist, "Double") flds = ("SHAPE@", fld_X, fld_Y, fld_Z, fld_A, fld_dist) lst_lines = [] with arcpy.da.InsertCursor(fc_pnt, flds) as curs:     # open txt file     arr_line = arcpy.Array()     i = 0     with open(txt, 'r') as f:         for r in f.readlines():             i += 1             r = r.replace('\n', '')             lst = r.split(' ')             if i % 1000 == 0:                 print "Processing line: {0}".format(i)             if i == 1:                 x = float(lst[3])                 y = float(lst[5])                 xp = x                 yp = y             else:                 xp = x                 yp = y                 x = float(lst[3])                 y = float(lst[5])             z = float(lst[10])             a = int(lst[12])             dist = math.hypot(xp-x, yp-y)             pnt = arcpy.Point(x, y, z)             pnt_g = arcpy.PointGeometry(pnt)             curs.insertRow((pnt_g, x, y, z, a, dist, )) ##            if dist < max_dist: ##                # add to line ##                arr_line.add(pnt) ##            else: ##                # add previous line to list ##                if arr_line.count > 1: ##                    polyline = arcpy.Polyline(arr_line) ##                    lst_lines.append(polyline) ## ##                # create a new line ##                arr_line = arcpy.Array() ##                arr_line.add(pnt)             if a == 1:                 # add to line                 arr_line.add(pnt)             else:                 # add previous line to list                 if arr_line.count > 1:                     polyline = arcpy.Polyline(arr_line)                     lst_lines.append(polyline)                 # create a new line                 arr_line = arcpy.Array()                 arr_line.add(pnt) arcpy.CopyFeatures_management(lst_lines, fc_lines) ... View more

It is not necessary to create a double loop, but to provide a sample that does not have that much changed see the script below. I introduced a dictionary "dct" holding the input name and the output name as key, value pair in the second loop the fc is split into the workspace and the featureclass name if the old name is in the dictionary (in the keys) then it is renamed import arcpy import os workspace = "D:\GIS\ZOne1\Zone1A" feature_classes = [] walk = arcpy.da.Walk(workspace, datatype="FeatureClass", type="Polygon") dct = {'featureclassName': 'featureclassName_renameTest',       'featureclassName1': 'featureclassName1_renameTest',       'another input name': 'corresponding output name'} for dirpath, dirnames, filenames in walk:     for filename in filenames:         feature_classes.append(os.path.join(dirpath, filename)) for fc in feature_classes:     fc_ws, fc_name = os.path.split(fc)     if fc_name in dct:         arcpy.Rename_management(fc, os.path.join(fc_ws, dct[fc_name])) ... View more

Hi, First of all you should really use syntax highlighting: Posting Code blocks in the new GeoNet To recursively loop through folder and sub folder and detect the layerfiles (.lyr) you can use the arcpy.da.Walk functionality (asuming you have a recent version of ArcGIS). import arcpy import os workspace = r"D:\Xander\GeoNet" walk = arcpy.da.Walk(workspace, datatype="Layer") for dirpath, dirnames, filenames in walk:     for filename in filenames:         print (os.path.join(dirpath, filename)) This prints each .lyr file it finds: D:\Xander\GeoNet\StackedLabels\aLayer.lyr D:\Xander\GeoNet\StackedLabels\subfolder\subsubfolder\subsubsubfolder\AnotherLayer.lyr In case you want to add all the layer found to an existing MXD you should open the MXD before the loop and add put the add layerfile logic inside the loop. Adding a lot of layer to an MXD will probably make it very slow. Your code might look like this (in case you run this inside a session of ArcMap and want to add the layers to the current session of ArcMap, hence the keyword "CURRENT"): import arcpy import os # suppose you want to add it to the current MXD (open MXD) mxd = arcpy.mapping.MapDocument("CURRENT") dataFrame = arcpy.mapping.ListDataFrames(mxd, "*")[0] # base folder workspace = r"D:\Xander\GeoNet" walk = arcpy.da.Walk(workspace, datatype="Layer") for dirpath, dirnames, filenames in walk:     for filename in filenames:         layerfile = os.path.join(dirpath, filename)         addlayer = arcpy.mapping.Layer(layerfile)         arcpy.mapping.AddLayer(dataFrame, addlayer, "BOTTOM") arcpy.RefreshTOC() arcpy.RefreshActiveView() del addlayer, mxd In case you want to run it as a standalone script and/or add the layers to a different existing MXD, you should replace the keyword "CURRENT" on line 5 by a reference (string containing path and filename) to the MXD. ... View more