Latest Contributions by XanderBakker

From the image you include I think there are 4 different triangles that can be created, and you could change the start/end points for each triangle. There is a thing to keep in mind when doing a cost distance and a subsequent least cost path analysis including a weight based on the slope. In a slope map the cells hold the maximum slope of that cell in relation to its neighbors. You can however cross a cell in different ways (uphill, downhill or going straight). To compensate you could combine an aspect map with the slopes and assume that that a certain aspect is uphill and the opposite is downhill... Below a link to a thread on hiking speeds, that might be interesting: Hiking equation in python returns wrong results‌ Kind regards, Xander ... View more

Sounds nice! To connect A to B passing by the lake, you should probably calculate path A - Lake and Lake - B separately. To enhance the analysis you could include a visibility (viewshed) analysis where the areas with visibility to the lake get a lower weight than the other pixels. This way you will force the route to have more views to the lake. If this will actually create a better path, depends strongly on the topography and landcover. Kind regards, Xander ... View more

It might help if you post a sample of your data and explain in more detail what type of chart you want to create and what need to be done with the "repeats" (redundancy) in the data. Also include the version of the software you´re using. Kind regards, Xander ... View more

Hi Paul, I'm glad it worked. Be aware that it creates the output in the same folder as the input. A second run would try to perform the same action on both the input and output files from the previous run. It might be better to create a separate workspace (folder) for the outputs. Like this: def main():     # specify the workspace     ws = r"C:\Users\PaulBritten\Desktop\tenure\All reserves"     ws_out = r"C:\Users\PaulBritten\Desktop\tenure\YourOutputFolder" # should exist     arcpy.env.workspace = ws     arcpy.env.overwriteOutput = True     lst_fc = arcpy.ListFeatureClasses()     for fc in lst_fc:         name, ext = os.path.splitext(fc)         fc_out = "{0}_merged{1}".format(name, ext)         MergePolygonsOnDate(os.path.join(ws, fc), os.path.join(ws_out, fc_out)) ... View more

Owen is right, TIN's are not supported as parameter in a geoprocessing service. In case the TIN's are pretty much static you could include them in your script, without being parameters. If you have multiple TIN's then you could provide some intelligence (for instance the extent) in your script to use those needed for the provided input featureclass. Kind regards, Xander ... View more

the 5th edition was released in june 2013: Learning Python, 5th Edition‌ ... View more

When I try the code below which copies the source to dest1 and dest2 (featureclass in a feature dataset), it works without any problem in 10.2.2. import arcpy arcpy.env.overwriteOutput = True source = r"C:\Forum\overwrite\source.gdb\myFC" dest1 = r"C:\Forum\overwrite\dest1.gdb\myFC" dest2 = r"C:\Forum\overwrite\dest2.gdb\featDS\myFC" arcpy.Copy_management(in_data=source,out_data=dest1,data_type="FeatureClass") arcpy.Copy_management(in_data=source,out_data=dest2,data_type="FeatureClass") The only way I could reproduce this error was trying to write to the dest2 featureclass without using the feature dataset in the path. Like this: dest2 = r"C:\Forum\overwrite\dest2.gdb\myFC" Are you sure you are using the output feature dataset in referencing the output featureclass? Kind regards, Xander ... View more

Since I don't know what your data looks like and I don't know where it is located, I just created an arbitrary polygon (triangle) over a set of test data pulled from a different thread. When you run your original code this should create the shapefile "X:/Users/Graphics.shp". Is this shapefile created correctly? Does it overlay with your data in the shapefile "X:/Users/ROUTE_SHAPES.shp"? Is your data frame, where you are drawing your line geometry, in the same projection? Those would be the first places to check. If the code is now working satisfactory then there may have been something with perhaps a previous shapefile that may have been conflicting... Do you remove the "Graphics.shp" or have you configured the option to overwrite the output in case it exists? I included this in my code to avoid problems. Since you are creating a tool (which I suppose you want to use more than once) and the intermediate and output results are hard coded (those will be the same each time you run the tool) you will need to include some handling of overwriting existing files. Kind regards, Xander ... View more

Maybe it's because the clip features "Graphics.shp" is not assigned a coordinate system? You can do this by changing the line where you create the polygon to this: polygon = arcpy.Polygon(array, "X:/Users/ROUTE_SHAPES.shp") ... View more

Hi Paul, Don't worry, that what the forum is all about... to help one another. Have a look at the code below. The majority is placed inside a function "MergePolygonsOnDate". The main function (lines 13 - 23) contains the logic for looping through your featureclasses. It will create a list of all the featureclasses in you workspace. Based on the input name it will create the new output names. On line 23 for each featureclass the "MergePolygonsOnDate" is called. I haven't tested the code, so if you run into problems, just post it and I'll have a look. #------------------------------------------------------------------------------- # Name:        MergePolygonsOnDate_v3.py # Purpose:     Merge polygons on date range # # Author:      xbakker # # Created:     19/08/2014 #------------------------------------------------------------------------------- import arcpy, os from datetime import datetime, timedelta def main():     # specify the workspace     ws = r"C:\Users\PaulBritten\Desktop\tenure\All reserves"     arcpy.env.workspace = ws     arcpy.env.overwriteOutput = True     lst_fc = arcpy.ListFeatureClasses()     for fc in lst_fc:         name, ext = os.path.splitext(fc)         fc_out = "{0}_merged{1}".format(name, ext)         MergePolygonsOnDate(os.path.join(ws, fc), os.path.join(ws, fc_out)) def MergePolygonsOnDate(fc, fc_out):     # change settings     fld_start = "DateFrom" # create date field and fill it with content of StartDate     fld_end = "DateTo" # create date field and fill it with content of EndDate     fld_name = "Name"     fld_outname = "OutputName" # output name     len_outname = 75     # create empty output shapefile     ws_path, fc_out_name = os.path.split(fc_out)     arcpy.CreateFeatureclass_management(ws_path, fc_out_name,         "POLYGON", fc, "SAME_AS_TEMPLATE", "SAME_AS_TEMPLATE", fc)     # add string field     arcpy.AddField_management(fc_out, fld_outname, "TEXT", "", "", len_outname)     flds = ("SHAPE@", fld_start, fld_end, fld_name)     flds_out = ("SHAPE@", fld_start, fld_name, fld_outname)     # create a unique list of dates from start date     lst_dates = list(set([row[0] for row in arcpy.da.SearchCursor(fc, (fld_start))]))     lst_dates.sort()     # insert cursor to store features to new featureclass     with arcpy.da.InsertCursor(fc_out, (flds_out)) as curs_out:         # loop through list of unique dates         for date in lst_dates:             # create an expression             expression = "{0} <= date '{2}' AND {1} >= date '{2}'".format(                 arcpy.AddFieldDelimiters(fc, fld_start),                 arcpy.AddFieldDelimiters(fc, fld_end), date)             i = 0             # create the searchcursor with the expression             with arcpy.da.SearchCursor(fc, flds, where_clause=expression) as curs:                 for row in curs:                     pol = row[0]                     name = row[3] # take name of first polygon                     if i == 0:                         # first polygon is just the polygon"                         polygon = pol                     else:                         # next polygons are added with union                         polygon = polygon.union(pol)                     i += 1                 # insert the feature into the output featureclass                 datefrom = date                 i_date = lst_dates.index(date)                 if i_date + 1 > len(lst_dates) - 1:                     dateto = ""                     outname = "{0} {1} - {2}".format(name, datefrom.strftime('%Y/%m/%d'), dateto)                 else:                     dateto = lst_dates[i_date + 1] - timedelta(days=1)                     outname = "{0} {1} - {2}".format(name, datefrom.strftime('%Y/%m/%d'), dateto.strftime('%Y/%m/%d'))                 curs_out.insertRow((polygon, row[1], name, outname)) if __name__ == '__main__':     main() ... View more

Hi Lance, As long as you know what you need, there is no problem whatsoever on summarizing on unique length. Good luck, Xander ... View more

Hi Paul, That is possible, but it might be necessary to know a little more on the folder structure where you have your input featureclasses and how you want to name the outputs. Normally a arcpy.ListFeatureclasses() can be used to get a list of featureclasses in the current workspace. In case they are in multiple folders and sub folders the arcpy.da.Walk() is better. Based on the input name, you can create an outputname (like for instance inputname_merged). The current script can be rewritten as a function to process a single featureclass. This function will be called from the loop through the input featureclasses. Kind regards, Xander ... View more