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; } } }

Random Selection of Preexisting Points w/ Definition Query

RobertCelestin1

Hi everyone,

I want to use definition queries with this code, taken from this ESRI page, in order to randomly select a subset of preexisting points by count:

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">def</SPAN> SelectRandomByCount <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> count<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
 <SPAN class="keyword token">import</SPAN> random
 layerCount <SPAN class="operator token">=</SPAN> int <SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>GetCount_management <SPAN class="punctuation token">(</SPAN>layer<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><SPAN class="punctuation token">)</SPAN>
 <SPAN class="keyword token">if</SPAN> layerCount <SPAN class="operator token"><</SPAN> count<SPAN class="punctuation token">:</SPAN>
 <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"input count is greater than layer count"</SPAN>
 <SPAN class="keyword token">return</SPAN>
 oids <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>oid <SPAN class="keyword token">for</SPAN> oid<SPAN class="punctuation token">,</SPAN> <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"OID@"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>
 oidFldName <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>OIDFieldName
 delimOidFld <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>AddFieldDelimiters <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> oidFldName<SPAN class="punctuation token">)</SPAN>
 randOids <SPAN class="operator token">=</SPAN> random<SPAN class="punctuation token">.</SPAN>sample <SPAN class="punctuation token">(</SPAN>oids<SPAN class="punctuation token">,</SPAN> count<SPAN class="punctuation token">)</SPAN>
 oidsStr <SPAN class="operator token">=</SPAN> <SPAN class="string token">", "</SPAN><SPAN class="punctuation token">.</SPAN>join <SPAN class="punctuation token">(</SPAN>map <SPAN class="punctuation token">(</SPAN>str<SPAN class="punctuation token">,</SPAN> randOids<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
 sql <SPAN class="operator token">=</SPAN> <SPAN class="string token">"{0} IN ({1})"</SPAN><SPAN class="punctuation token">.</SPAN>format <SPAN class="punctuation token">(</SPAN>delimOidFld<SPAN class="punctuation token">,</SPAN> oidsStr<SPAN class="punctuation token">)</SPAN>
 arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByAttribute_management <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> sql<SPAN class="punctuation token">)</SPAN>
__________________________________________________________________
SelectRandomByPercent <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"layer"</SPAN><SPAN class="punctuation token">,</SPAN> num<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>

The code works with a shapefile, but does not work for me with a personal geodatabase.
My desired scenario with using a definition query is as follows:

I have 27,000 points in one layer and 511 polygons in another. The points are draw over the polygons
and are only contained within the boundaries of the polygons.

Phase One:
1. I want to first filter out a subset of points, using a simple definition query (e.g. "ID = #") that
is set via the layer properties
2. I want to randomly select a max number of points from the layer that has the definition query set,
using the script above
3. Once that selection is made, I will assign the selected points a unique name in a field i've created.

Phase Two:
1. I want to set a definition query and repeat Phase One - Step 1. with the addition of excluding
the points assigned a name.

Can anyone help with this? I've attached a representation of the points and polygons. Thank you

Best,
Robert

Water2.JPG

Water.JPG

table example.png

Index example.png

Encuentra más publicaciones etiquetadas con

Accepted answers

PriscillaCole

Hi Robert,

I tested your code on a personal geodatabase and I found where it's getting stuck. The arcpy field delimiters function isn't returning the corrector delimiter. Personal geodatabases need square brackets wrapping their fields for SQL queries, but instead it's wrapping with double quotes. I added in a test to detect if the datasource is coming from a .mdb.

NOTE - I plan to report to ESRI about the field delimiter bug. This arcpy function shouldn't be behaving in this way.

<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</SPAN> random

<SPAN class="keyword token">def</SPAN> SelectRandomByCount <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> count<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;layerCount <SPAN class="operator token">=</SPAN> int <SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>GetCount_management <SPAN class="punctuation token">(</SPAN>layer<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><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">if</SPAN> layerCount <SPAN class="operator token"><</SPAN> count<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"input count is greater than layer count"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;oids <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>oid <SPAN class="keyword token">for</SPAN> oid<SPAN class="punctuation token">,</SPAN> <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"OID@"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;desc <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">if</SPAN> <SPAN class="string token">".mdb"</SPAN> <SPAN class="keyword token">in</SPAN> desc<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">:</SPAN> 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;name <SPAN class="operator token">=</SPAN> <SPAN class="string token">"["</SPAN> <SPAN class="operator token">+</SPAN> desc<SPAN class="punctuation token">.</SPAN>OIDFieldName <SPAN class="operator token">+</SPAN> <SPAN class="string token">"]"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN> 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;name <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>AddFieldDelimiters <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> desc<SPAN class="punctuation token">.</SPAN>OIDFieldName<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;randOids <SPAN class="operator token">=</SPAN> random<SPAN class="punctuation token">.</SPAN>sample <SPAN class="punctuation token">(</SPAN>oids<SPAN class="punctuation token">,</SPAN> count<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;oidsStr <SPAN class="operator token">=</SPAN> <SPAN class="string token">", "</SPAN><SPAN class="punctuation token">.</SPAN>join <SPAN class="punctuation token">(</SPAN>map <SPAN class="punctuation token">(</SPAN>str<SPAN class="punctuation token">,</SPAN> randOids<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;sql <SPAN class="operator token">=</SPAN> <SPAN class="string token">"{0} IN ({1})"</SPAN><SPAN class="punctuation token">.</SPAN>format <SPAN class="punctuation token">(</SPAN>name<SPAN class="punctuation token">,</SPAN> oidsStr<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByAttribute_management <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> sql<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>

All comments

JoshuaBixby

You mention definition query several times, but you are posting Python code. Since definition queries are written in SQL, are you using the Python code to return a list of OIDs to base the definition query on? Reading over Phase One, you talk about setting a definition query in Step 2 but then making a selection in Step 3.

Instead of describing the process you think you need to follow, can you summarize what you want the outcome to be? My impression is that you want to randomly group features in a feature class and assign each group a new label in a new field, is that what you are trying to do?

RandyBurton

The code works with a shapefile, but does not work for me with a personal geodatabase.

Can you explain "does not work" a bit more. Are you getting an error message; if so, what? Are you working in the Python window in ArcMap, or as a stand-alone script?

The code appears to work with a file geodatabase. I was successful in selecting a random group of points in a feature layer. If I understand what you want to do, it seems that line 14 in your code could be an update cursor that uses the "sql" where clause to update a field.

RobertCelestin1

Hi, thanks for replying.

Sorry for the confusion. I am not using the python code to return an OID to use in a definition query. The definition query is set before the python code is used, in order to first confine the set of points I am working with. The first image that I attached to my post is an example of that confined set of points. The second image is an example of what results from using the python code on the points.

You are correct in your impression of what i'd like the outcome to be. However, I do not want to select/assign all points to groups at the same time. To properly reiterate my desired outcome: I would like to randomly select some points (no greater than a max number that I define) from a specific set of points (a subset of the 27,000). Then assign that selection a label in a new field. After this, I need to repeat the process of random selection on the same initial specific set of points, but exclude those points that have already been assigned a label.

Whenever I go through the process of random selection, I only wanted to work with a subset of points from the layer of 27,000. This is why I thought to use a definition query.

PriscillaCole

Hi Robert,

NOTE - I plan to report to ESRI about the field delimiter bug. This arcpy function shouldn't be behaving in this way.

<SPAN class="keyword token">import</SPAN> arcpy<SPAN class="punctuation token">,</SPAN> random

<SPAN class="keyword token">def</SPAN> SelectRandomByCount <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> count<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;layerCount <SPAN class="operator token">=</SPAN> int <SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>GetCount_management <SPAN class="punctuation token">(</SPAN>layer<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><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">if</SPAN> layerCount <SPAN class="operator token"><</SPAN> count<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">print</SPAN> <SPAN class="string token">"input count is greater than layer count"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;oids <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>oid <SPAN class="keyword token">for</SPAN> oid<SPAN class="punctuation token">,</SPAN> <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"OID@"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;desc <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">if</SPAN> <SPAN class="string token">".mdb"</SPAN> <SPAN class="keyword token">in</SPAN> desc<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">:</SPAN> 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;name <SPAN class="operator token">=</SPAN> <SPAN class="string token">"["</SPAN> <SPAN class="operator token">+</SPAN> desc<SPAN class="punctuation token">.</SPAN>OIDFieldName <SPAN class="operator token">+</SPAN> <SPAN class="string token">"]"</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN> 
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;name <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>AddFieldDelimiters <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> desc<SPAN class="punctuation token">.</SPAN>OIDFieldName<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;randOids <SPAN class="operator token">=</SPAN> random<SPAN class="punctuation token">.</SPAN>sample <SPAN class="punctuation token">(</SPAN>oids<SPAN class="punctuation token">,</SPAN> count<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;oidsStr <SPAN class="operator token">=</SPAN> <SPAN class="string token">", "</SPAN><SPAN class="punctuation token">.</SPAN>join <SPAN class="punctuation token">(</SPAN>map <SPAN class="punctuation token">(</SPAN>str<SPAN class="punctuation token">,</SPAN> randOids<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;sql <SPAN class="operator token">=</SPAN> <SPAN class="string token">"{0} IN ({1})"</SPAN><SPAN class="punctuation token">.</SPAN>format <SPAN class="punctuation token">(</SPAN>name<SPAN class="punctuation token">,</SPAN> oidsStr<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByAttribute_management <SPAN class="punctuation token">(</SPAN>layer<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> sql<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>

SteveLynch

if you have a Geostatistical Analyst license you could use the SubsetFeatures geoprocessing tool.

RobertCelestin1

Hi Priscilla,

Thank you for the modification, it worked quite well including setting a def query before executing. And thanks for contacting ESRI about this.

Best,

Robert

Únete a la Conversación

Los miembros registrados pueden publicar, seguir actualizaciones y más. ¿Nuevo aquí? Regístrate gratis.

Enlaces Rápidos

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