Latest Contributions by XanderBakker

You can do something with a Stacked bar in ArcMap: But this is only posible if for each point the soils have the same order, and there is a column for each soil type and it contains the thickness of the layer. Since that is not realistic, you could do this with extrusión, but it requires multiple points on the same location containing data on soilt ype, depth from and depth to. Something similar as in slide 11 of this presentation (in Dutch, sorry): Bodemdaling model Provincie Utrecht - Grontmij . I created a featureclass with points for every grid cell and extracted all the values and created a related table holding the soil layer for each point. ... View more

If it is a matter of speeding up the copy process of large files, you may want to look at this article: Faster Python File Copy - Blumetech's Tech Blog If you want to record the time it took to copy the file, you can record the current time just before and just after and determine the timedelta or read the post by Owen Earley‌: https://community.esri.com/thread/119781#448946 ... View more

You can use some python to obtain the points you want: import arcpy, os fc_in = r"C:\Forum\Pasture\gdb\Contours.gdb\contours2" fc_out = r"C:\Forum\Pasture\gdb\Contours.gdb\contour_pnts03" fld_contour = "Contour" interval = 1 # create empty output featureclass sr = arcpy.Describe(fc_in).spatialReference out_ws, out_name = os.path.split(fc_out) arcpy.CreateFeatureclass_management(out_ws, out_name, "POINT", fc_in, "DISABLED", "DISABLED", sr) # add field with contour value to output fc arcpy.AddField_management(fc_out, fld_contour, "DOUBLE") cnt = 0 with arcpy.da.InsertCursor(fc_out, ("SHAPE@", fld_contour)) as curs_out:     with arcpy.da.SearchCursor(fc_in, ("SHAPE@", fld_contour)) as curs_in:         for row_in in curs_in:             polyline = row_in[0]             contour = row_in[1]             max_len = polyline.length             d = 0             while d < max_len:                 pnt = polyline.positionAlongLine(d, False)                 curs_out.insertRow((pnt.firstPoint, contour, ))                 cnt += 1                 if cnt % 100 == 0:                     print "Processing point: {0}".format(cnt)                 # pnt.firstPoint                 d += interval             curs_out.insertRow((polyline.lastPoint, contour, ))             cnt += 1 Change the lines 3 and 4 to point to your input featureclass and the output you want to create. ... View more

... and if you want to create an output that contains the attributes of the input, you can try this: import arcpy, os def main():     arcpy.env.overwriteOutput = True     # inputs     fc_in = r'D:\Xander\EPM\Datos pruebas\gdb\datos erase.gdb\points03'     width = 0.005 # with of the output rectangle     height = 0.005 # height of the output rectangle     # output     fc_out = r'D:\Xander\EPM\Datos pruebas\gdb\datos erase.gdb\rectangles3'     # create empty output featureclass     sr = arcpy.Describe(fc_in).spatialReference     out_ws, out_name = os.path.split(fc_out)     arcpy.CreateFeatureclass_management(out_ws, out_name, "POLYGON", fc_in, "DISABLED", "DISABLED", sr)     # make fields list     flds_use = correctFieldList(arcpy.ListFields(fc_in))     # create insert cursor     with arcpy.da.InsertCursor(fc_out, flds_use) as curs_out:         # loop through input features         with arcpy.da.SearchCursor(fc_in, flds_use) as curs_in:             for row_in in curs_in:                 polygon = createRectangleFromPoint(row_in[0].firstPoint, width, height, sr)                 row_out = buildRow(row_in, flds_use, polygon)                 curs_out.insertRow(row_out) def createRectangleFromPoint(pnt, width, height, sr):     arrPnts = arcpy.Array()     pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)     arrPnts.add(pnt2)     pnt2 = arcpy.Point(pnt.X - width, pnt.Y + height)     arrPnts.add(pnt2)     pnt2 = arcpy.Point(pnt.X + width, pnt.Y + height)     arrPnts.add(pnt2)     pnt2 = arcpy.Point(pnt.X + width, pnt.Y - height)     arrPnts.add(pnt2)     pnt2 = arcpy.Point(pnt.X - width, pnt.Y - height)     arrPnts.add(pnt2)     return arcpy.Polygon(arrPnts, sr) def correctFieldList(flds):     flds_use = ['Shape@']     fldtypes_not = ['Geometry', 'Guid', 'OID']     for fld in flds:         if not fld.type in fldtypes_not:             flds_use.append(fld.name)     return flds_use def buildRow(row_in, flds_use, polygon):     try:         lst_in = list(row_in)         lst_in[0] = polygon         return tuple(lst_in)     except Exception, e:         print("buildRow exception: {0}".format(e))         return None if __name__ == '__main__':     main() ... View more