nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Returning the first unique result with searchcursor?

BenjaminBrewer

Looking for a way to return the next unique result with a searchcursor in arcpy.

The table is from a spatial join of areas and has been sorted through python already.

However, I would like a way to have the name of the result recorded to a text file (or other dictionary) and to only take the next result each time the script runs and to check the dictionary to determine the next row to take.

Therefore allowing me to work row by row, one at a time, where they will be exported to a GeoJSON file through quick export.

I am struggling to find a way to make this executable.

Find more posts tagged with

Accepted answers

All comments

BenjaminBrewer

Hi Joshua, yes the join is static once it has been completed. In the real world scenario this may not be the case but from what I want to achieve here then it is more than appropriate. I had thought I was going a bit over complicated and that is why I was struggling to find examples. I will be trying to implement this further, thank you for the link you provided too, gives me more to read up on and work out from there.

Thank you for your help. I will update with what is working best once I have adapted the changes.

BenjaminBrewer

I had not thought about deleting data as I wanted to keep the original dataset complete but this could be the option I need. So with the print row.FILE_NAME I can take the selected area and then export it as I wish, then as it deletes that row from the table the next time it runs it will just start from the next available row. Thank you, I will give this a try.

JoshuaBixby

The feature class is static, right? That is, the results of the spatial join aren't changing between times you extract the result for a new township. If the spatial join is static, you really want to process the whole feature class once, create a single output file with the results, and then go through that file one at a time. Trying to implement the specific steps you describe complicates the scripting and would perform poorly since you are constantly going back to the feature class.

This type of question gets asked every 6 months or so on GeoNet. GroupBy 2 columns and keep all fields has a good discussion of different approaches using ArcPy cursors and NumPy.

Adapted from one of the code snippets in the thread mentioned above, the following creates a layer and selects records. After the records are selected, you can export them however you like.

<SPAN class="keyword token">from</SPAN> itertools <SPAN class="keyword token">import</SPAN> groupby
<SPAN class="keyword token">from</SPAN> operator <SPAN class="keyword token">import</SPAN> itemgetter

fc <SPAN class="operator token">=</SPAN> <SPAN class="comment token"># path to feature class</SPAN>
case_fields <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">"DESCRIPTOR"</SPAN><SPAN class="punctuation token">]</SPAN>
max_field <SPAN class="operator token">=</SPAN> <SPAN class="string token">"Count_"</SPAN>
sql_orderby <SPAN class="operator token">=</SPAN> <SPAN class="string token">"ORDER BY {}, {} DESC"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN><SPAN class="string token">", "</SPAN><SPAN class="punctuation token">.</SPAN>join<SPAN class="punctuation token">(</SPAN>case_fields<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> max_field<SPAN class="punctuation token">)</SPAN>

fl <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>MakeFeatureLayer_management<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"fc_lyr"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>getOutput<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN>
oidField <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>fl<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>OIDFieldName
oids <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>

<SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fl<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"*"</SPAN><SPAN class="punctuation token">,</SPAN> sql_clause<SPAN class="operator token">=</SPAN><SPAN class="punctuation token">(</SPAN>None<SPAN class="punctuation token">,</SPAN> sql_orderby<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> cur<SPAN class="punctuation token">:</SPAN>
    oidField <SPAN class="operator token">=</SPAN> cur<SPAN class="punctuation token">.</SPAN>fields<SPAN class="punctuation token">.</SPAN>index<SPAN class="punctuation token">(</SPAN>oidField<SPAN class="punctuation token">)</SPAN>
    case_func <SPAN class="operator token">=</SPAN> itemgetter<SPAN class="punctuation token">(</SPAN>
        <SPAN class="operator token">*</SPAN><SPAN class="punctuation token">(</SPAN>cur<SPAN class="punctuation token">.</SPAN>fields<SPAN class="punctuation token">.</SPAN>index<SPAN class="punctuation token">(</SPAN>fld<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">for</SPAN> fld <SPAN class="keyword token">in</SPAN> case_fields<SPAN class="punctuation token">)</SPAN>
    <SPAN class="punctuation token">)</SPAN>
    <SPAN class="keyword token">for</SPAN> key<SPAN class="punctuation token">,</SPAN> group <SPAN class="keyword token">in</SPAN> groupby<SPAN class="punctuation token">(</SPAN>cur<SPAN class="punctuation token">,</SPAN> case_func<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
        oids<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>next<SPAN class="punctuation token">(</SPAN>group<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN>oidField<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
        <SPAN class="keyword token">continue</SPAN>

fl<SPAN class="punctuation token">.</SPAN>setSelectionSet<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"NEW_SELECTION"</SPAN><SPAN class="punctuation token">,</SPAN> oids<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

JB1234510

You could use the code of Rhett Zufelt and go one with writing every unique value to a list and use the max of the list to find the highest count.

<SPAN class="keyword token">for</SPAN> val <SPAN class="keyword token">in</SPAN> Uni<SPAN class="punctuation token">:</SPAN>

    MyList<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>val<SPAN class="punctuation token">)</SPAN>

MaxValue <SPAN class="operator token">=</SPAN> max<SPAN class="punctuation token">(</SPAN>MyList<SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"Max value element : "</SPAN><SPAN class="punctuation token">,</SPAN> MaxValue‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

After that you can use MaxValue to do a select by attribute on your featureclass which I would copy before. 

Let you give out the name of the selected row with another search cursor

for row in arcpy.SearchCursor(fc_Layer):     print row.FILE_NAME

Write this to a text file

Delete the selected feature and go to the next MaxValue

   

BenjaminBrewer

Hi Joshua and Rhett Zufelt‌, sorry I have had a really busy day in school. Thank you for looking over at my question.

To clarify some more of what it is I am trying to do.

I have a feature class, named TWP_Messages which is a spatial join of GPS points and the ATS grid system in Alberta.

I have sorted this feature class by a count field for the amount of GPS messages in each township. I then want to take the top result of that join and export it into a GeoJSON file which I will then use, incorporated with a script I have written to download satellite imagery from Planet Labs. I however only want to do this one township (1 cell of the ATS grid) at a time so not to be overloaded.

To achieve this I was hoping to find a way that a search cursor, or similar will take the first area from the spatial join, but as it exports this area it writes the area name in to a text file. I then ideally would like that text file to be used as a check list so when I next run the script it will then look at the text file, know that it has taken the top ranked area from the sorted join feature class and then move on to the next area in the list.

Spatial join between townships and GPS messages

So in this feature class the first time it runs it will select the row for object ID 1520, the second time the script runs I would like for it to be able to check which areas have already been selected, so the first row in the table, ignore that and then move on to the second row so it can be exported.

I hope that is easier to understand. I am sure it can be done but I cant seem to find a similar instance where something similar has been done before.

Again, thank you for all your help. I hope I have described the situation better.

JoshuaBixby

Benjamin Brewer‌, as the discussion between Rhett Zufelt‌ and I shows, there are multiple ways to interpret what you are trying to accomplish. What would be most helpful is to provide either an actual or made-up example that you and commenters can respond to directly.

rzufelt

thus the comment about not being clear about what is being asked for.

Guess I'm not sure what the 'next unique result' is. Is that the next, non-equal consecutive value in the sorted table?

the first time a value is encountered in the sorted table, etc?

OP mentions wanting to use the next result 'each time' and work row by row (assuming across the entire table).

Based on the vagueness of the question, it sounded like OP wanted to be able to perform export operations for rows with a specified value in a particular field. A set would give all the unique values loaded from that field.

Iterating through the set will give you the next value, then the next, etc. of the set itself (not necessarily the order they were in the original table), but would give a list of unique values to query on......

of course, this could all be mute depending on the actual question.

JoshuaBixby

Python sets are unordered, so the posted code snippet will provide a list of unique values for a field(s) but it won't provide the "next unique value."

rzufelt

Hard to really interpret what you are asking for, but if you really want to find the "next" unique value in a field, here is one method:

<SPAN class="keyword token">import</SPAN> arcpy

fc <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\Users\xyz\Desktop\tmp.gdb\cc4'</SPAN>
fields <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">'_type'</SPAN><SPAN class="punctuation token">]</SPAN>

<SPAN class="comment token"># load all values from the listed field into List with SearchCursor</SPAN>
<SPAN class="comment token"># then, makes a set from the list (this gives only unique values)</SPAN>

Uni <SPAN class="operator token">=</SPAN> set<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>row<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="keyword token">for</SPAN> row <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> fields<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>


<SPAN class="comment token"># since we now have a set with the unique values only, you can</SPAN>
<SPAN class="comment token"># iterate through that set now and do as you plese with each one</SPAN>



<SPAN class="keyword token">for</SPAN> val <SPAN class="keyword token">in</SPAN> Uni<SPAN class="punctuation token">:</SPAN>
    <SPAN class="comment token">#Do something with this value</SPAN>
    <SPAN class="comment token"># or, write to text file</SPAN>
    <SPAN class="comment token"># or, add to dictionary</SPAN>
    <SPAN class="comment token"># etc.</SPAN>
    ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.