Latest Contributions by DuncanHornby

I can't replicate your problem so it suggests there is another issue. May be data type or default settings? When I run extract by mask from a python command line it works as expected. ... View more

I wonder rather than saying null, it's coloured or in italics. I hear what you are saying but seeing null might think I could put a point in when it's a polygon feature class or something daft like that? ... View more

OK, well done for showing your workings and its because you made the effort to show me I can make an educated guess! Your layername is "New_parval_SpatialJoin". The clue is in the name... I'm guessing you have done a spatial join with the new Add Spatial Join tool?  Well if you open your attribute table and move your cursor over the Tree_can field header, I'm going to guess its not Tree_can but something like "somelayer_spatial_Join.Tree_can"? Thats the fully qualified field name of the field because its a joined field. So what you see is what you need to use in your cursors. ... View more

You can complete this task with some python scripting which you can run directly within ArcPro from the python console. Here is an example table of my input layer, 12 rows with the values I want to put into quartile classes in the field Value and QClass is the field to record which class they are in, note the rows are all NULL as the code has yet to be run.:   The code you drop into the python console is this, you just need to change the layername and fields in the search and update cursors to match your data.:   import arcpy, numpy layername = "testdata" # Name of layer as seen in TOC lstValues = list() # list to store values in # Read values into list with arcpy.da.SearchCursor(layername, "Value") as cursor: for row in cursor: lstValues.append(row[0]) # Use numpy to compute quartiles q25 = numpy.percentile(lstValues, 25) q50 = numpy.percentile(lstValues, 50) q75 = numpy.percentile(lstValues, 75) print("Q25 = " + str(q25)) print("Q50 = " + str(q50)) print("Q75 = " + str(q75)) # Write quartile class back to layer with arcpy.da.UpdateCursor(layername, ["Value","QClass"]) as cursor: for row in cursor: v = row[0] # The value if v < q25: row[1] = 25 elif v >= q25 and v < q50: row[1] = 50 elif v >= q50 and v < q75: row[1] = 75 else: row[1] = 100 cursor.updateRow(row) Print("Finished!")     The results are shown, sorted by QClass:   ... View more

What you described work for me, my input rasters were tested as TIF and file geodatabase rasters, fed into a collects tool and that sub model connected correct with cell statistics. It would suggest the problem is in your sub-model and as you don't show it, no one can help you. ... View more

I think the short answer is no. It's a cosmetic bug in model builder. You can click on that single line and drag out the links to new positions. ... View more

I think to give people a chance to help you, you need to attach to your question the layer file with the symbology as what you show could no doubt be created multiple ways and which way you have done it will dictate the solution (if any). ... View more

If you look at the tool in toolbox you will observe that it has a script icon which means its a python script and these you can look at the source to fully understand what the tool is doing. Simply right click on tools and select > properties > Execution. If you search the source code for the lines that have arcpy.Polyline() this is where they construct the polyline geometry from an array object and if you search for the array you simply observe them adding points to it. So no fancy conversion is happening. From this you can conclude that the lines built are in the coordinate system you have provided the points in, in your case looks like WGS84.  You need to go back to the person who supplied them to you to see if they were constructed in a specific manner. ... View more

You could explore the Collect Events tool to count your stacked points then use that new layer to select only 2-count points then you can do the select by location using those selected subset? I have to say questions like yours are ALWAYS better if accompanied with a screen shot/sketch  as you have obfuscated the nature of your dataset by not saying what it is, rivers, roads, shipping lanes etc... ... View more

Thanks for the advice, this was the approach I was hoping to avoid! 😁 Interesting stuff about I/O. ... View more

If you want assistance, it will help if you state which software you require assistance in as the syntax is applicable to ArcMap or ArcPro... Edit your question, don't just simply reply to this as others then have to trawl through a pointless conversation. ALSO format your code (via the 3 dots), its difficult to read and indentation is critical in python. A question where you have put time and effort in formatting and providing all the required information gets answered a lot quicker! ... View more

That's a really useful white paper, some quality bed time reading for me! Thanks. ... View more

To loop within models you need to be using an iterator. Explore help file and look at the examples. ... View more

This page discusses how ArcPro deals with precision and tolerance: https://pro.arcgis.com/en/pro-app/latest/help/data/geodatabases/overview/the-properties-of-a-spatial-reference.htm Without talking to the developers that's as much as you are going to find out which is published. ... View more