Particle Track arcpy.point failing in a loop

Question

I'm using the particle track in python and limiting it to a specific trackingtime. For each tool run, I want the rasters input (magnitude and direction) to change to the next rasters in the lists, and the starting point of the tracking (x,y) to be derived from the last polyline point of the previous run. Attached is what I have so far but I get Error999999. It works when arcpy.point is not in a loop but fails within it. Any help would be very appreciated. for idx in range(0, len(mag)): print "working on index "+str(idx) inMagnitudeRaster= mag [idx] print inMagnitudeRaster inDirectionRaster= direction [idx] print inDirectionRaster basename= inDirectionRaster.split(".tif")[0] outTrackPolylineFeatures = os.path.join (OutputFolder, basename + '.shp') outTrackFile= os.path.join(OutputFolder, basename +'.txt') ParticleTrack (Raster(inDirectionRaster), os.path.join(r'C:/Time_slices/January_newrun/Magnitude2/Rectified/Stripped',inMagnitudeRaster), sourcePoint,outTrackFile, "#", trackingTime,outTrackPolylineFeatures) arcpy.AddField_management (outTrackPolylineFeatures, "LastX", "DOUBLE") arcpy.AddField_management (outTrackPolylineFeatures, "LastY", "DOUBLE") arcpy.CalculateField_management(outTrackPolylineFeatures, "LastX", "(!SHAPE.lastpoint.X!)", "PYTHON_9.3") arcpy.CalculateField_management(outTrackPolylineFeatures, "LastY", "(!SHAPE.lastpoint.Y!)", "PYTHON_9.3") endx, endy = [i for i in arcpy.da.SearchCursor(outTrackPolylineFeatures,['LastX','LastY'])][-1] sourcePoint= arcpy.Point(float(endx),float(endy),None,None)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ particle_cursor_18_02_v5.py.zip S123 AGOL sign in error.jpg

Kris01 · Answer

Thanks Joshua. That might have been the problem. I edited the loop as I realised I only need one time step in the shapefile. So all but the first row were deleted and now the code works. for idx in range(0, len(mag)): print "working on index "+str(idx) inMagnitudeRaster= mag [idx] print inMagnitudeRaster inDirectionRaster= direction [idx] print inDirectionRaster basename= inDirectionRaster.split(".tif")[0] outTrackPolylineFeatures = os.path.join (OutputFolder, basename + '.shp') outTrackFile= os.path.join(OutputFolder, basename +'.txt') ParticleTrack (Raster(inDirectionRaster), os.path.join(r'C:/Leverhulme/Time_slices/January_newrun/Magnitude2/Rectified/Stripped',inMagnitudeRaster), sourcePoint,outTrackFile, stepLength, trackingTime,outTrackPolylineFeatures) with arcpy.da.UpdateCursor(outTrackPolylineFeatures, "FID") as cursor: for row in cursor: if row[0] > 0: cursor.deleteRow() arcpy.AddField_management (outTrackPolylineFeatures, "LastX", "DOUBLE",20 ) arcpy.AddField_management (outTrackPolylineFeatures, "LastY", "DOUBLE", 20) arcpy.CalculateField_management(outTrackPolylineFeatures, "LastX", "(!SHAPE.lastpoint.X!)", "PYTHON_9.3") arcpy.CalculateField_management(outTrackPolylineFeatures, "LastY", "(!SHAPE.lastpoint.Y!)", "PYTHON_9.3") endx, endy = [i for i in arcpy.da.SearchCursor(outTrackPolylineFeatures,['LastX','LastY'])][-1] sourcePoint= arcpy.Point(float(endx),float(endy),None,None)‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

JoshuaBixby · Answer

It helps to provide the error and traceback. I doubt the issue is with your ArcPy Point line. Looking at this line: endx, endy = [i for i in arcpy.da.SearchCursor(outTrackPolylineFeatures,['LastX','LastY'])][-1]‍ What exactly are you trying to achieve, are you trying to grab the last record of the search cursor? If so, a couple of thoughts, if your cursor only has 1 record you will run into a TypeError with your tuple unpacking. With multiple records but no SQL ORDER BY clause, you can't guarantee the last record will be the "last."

Kris01 · Answer

I did add a print statement and the correct x, y coordinates are printed. When I print the type, it's a float.