Latest Contributions by JohannesBierer

Hello, i've slightly adapted a script of Bera - Stackexchange - https://gis.stackexchange.com/questions/308705/split-a-line-at-a-specific-elevation It should set a point every 10 cm hight change on a line. The script is somehow working but not very efficient, is there a better way to do it? Originally the script was for one best point on a line ... How not to write single bestpoints to fclass? Any ideas? # erzeugt Punktshape alle x cm Hoehe entlang Linie import arcpy import numpy as np arcpy.env.overwriteOutput = True arcpy.env.workspace = r"Path\Berechnung_4.gdb" fc = r'Path_to_FGDB\Gew_Morre' inDGM = r"Path_to_FGDB\dgm01m" step = 1 # 0.1 in script - Change, as small as possible without making the execution time unacceptably long arcpy.CheckOutExtension("3D") Gewaesser_3D = "Gewaesser_3D" arcpy.InterpolateShape_3d(inDGM, fc, Gewaesser_3D) arcpy.AddSurfaceInformation_3d(Gewaesser_3D, inDGM, "Z_MIN;Z_MAX", "BILINEAR") arcpy.CheckInExtension("3D") with arcpy.da.SearchCursor(Gewaesser_3D,['SHAPE@','Z_MIN','Z_MAX']) as cursor: for row in cursor: print("valueMin = {} and valueMAx = {}".format(row[1], row[2])) arr = np.arange(row[1], row[2], 0.1) # hier Hoehenstufen definieren! print(arr) for r in arr: geom = row[0] points = [] position = 0 while position < geom.length: points.append(geom.positionAlongLine(position)) position+=step bestpoint = min(points, key=lambda x: abs(x.centroid.Z-r)) strR = round(r, 5) print (strR) ReplaceRstrR = str(strR).replace(".", "_") arcpy.CopyFeatures_management(in_features=bestpoint, out_feature_class='bestpoint{}'.format(ReplaceRstrR)) featureclasses = arcpy.ListFeatureClasses("bestpoint*") arcpy.Merge_management(featureclasses, "merge_ergebnis_bruchhof_10cm")   ... View more

Hello, I've problems with the "MULTIPOINT" option in arcpy.ASCII3DToFeatureClass_3d, it creates stranges results in the outgoing fc. See attached image in red there are squares create with the wrong hight information out of lidar vegetation (first pulse) data. In the image I analyse vegetation hight over 3m. If I use the "POINT" option the vegetation is correctly shown in the output raster data (black in the image). What's wrong, can't you use "MULTIPOINT" with lidar data with a resolution of 8 points/m²? Is it a bug? It takes very long time to process the data with the "POINT" option ...   ... View more

Hello, is there a size limit creating a terrain. I want to add about 3000 tiles of lidar masspoint feature classes? Or is there a other possibility to create a seamless ground point raster out of tiles of lidar xyz files? So far I use this tools: 1. arcpy.ASCII3DToFeatureClass_3d 2. arcpy.AddFeatureClassToTerrain_3d Best wishes johannes ... View more