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

Help using PointGeometry with SelectLayerByLocation

ZoltanSzecsei

Hi,

The following code fails:

pg = arcpy.PointGeometry(vert,'None',True)
arcpy.SelectLayerByLocation_management("hgt_feats", "WITHIN_A_DISTANCE_3D", pg, 0.005, 'NEW_SELECTION')

with:

File "Z:/Projects/Development/python/ToGround/ToGround.py", line 140, in <module>
arcpy.SelectLayerByLocation_management("hgt_feats", "WITHIN_A_DISTANCE_3D", pg, 0.005, 'NEW_SELECTION')
File "C:\Program Files (x86)\ArcGIS\Desktop10.6\ArcPy\arcpy\management.py", line 7873, in SelectLayerByLocation
raise e
RuntimeError: Object: Error in executing tool

I'm trying to see if there is an underlying polygon vertex at a specific XYZ location.

I do not need to select the underlying polygon, I just need to know if there is one at that location

Any help would be appreciated.

Regards and thanks,

Zoltan

Trouver plus de publications taguées avec

Accepted answers

DarrenWiens2

Try specifying a spatial reference on your PointGeometry, even if it is unprojected.

All comments

bixb0012

In your new thread, explain what isn't working reliably, the geoprocessing tools themselves or the spatial analysis part.

ZoltanSzecsei

OK - it seems not to work reliably, si I need to rethink my approach.

I'll start a new thread called "Efficient way to check if a vertex coincides with another feature's vertex"

ZoltanSzecsei

No - this was my first attempt - but see my post with screenshot, above - it does work.

ZoltanSzecsei

Yes it does Notice i==1 in watch window.

bixb0012

Should it work for Zoltan with the PointGeometry, I think so, and it appears to be running without errors.

I can't say I am a fan of the changes to Select Layer by Attributes and Selection Layer By Location. The tools are called "Select Layer," I don't see what is wrong with requiring an actual layer for the in_layer parameter, which is the historical behavior. The newfangled, "convenient" behavior both confuses the semantics of the tools as well as changes some of the behaviors of the tools when using ArcPy geometries as the select_features. I find it humorous and annoying how "backwards compatibility" is so often used as the excuse when my organization asks about changing ArcPy behavior, but it doesn't seem to stop Esri from changing long-standing tool behavior when they want.

DanPatterson_Retired

That makes sense (like the interactive featurelayer creation in tools for Pro

... but it does the arcpy's PointGeometry Zoltan is using work?

bixb0012

With ArcGIS Pro, not sure whether from 1.0 or a later release, you can use Select Layer By Location with only ArcPy geometries. The Select Layer By Location—Data Management toolbox | ArcGIS Pro alludes to it rather than stating it directly:

Usage
If the input is a feature class or dataset path, this tool will automatically create and return a new layer with the result of the tool applied.

<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> <SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> 
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> SR <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>SpatialReference<SPAN class="punctuation token">(</SPAN><SPAN class="number token">3857</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> polys <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>FromWKT<SPAN class="punctuation token">(</SPAN><SPAN class="string token">'POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))'</SPAN><SPAN class="punctuation token">,</SPAN> SR<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>FromWKT<SPAN class="punctuation token">(</SPAN><SPAN class="string token">'POLYGON((11 11, 21 11, 21 21, 11 21, 11 11))'</SPAN><SPAN class="punctuation token">,</SPAN> SR<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>FromWKT<SPAN class="punctuation token">(</SPAN><SPAN class="string token">'POLYGON((0 11, 5 11, 5 16, 0 16, 0 11))'</SPAN><SPAN class="punctuation token">,</SPAN> SR<SPAN class="punctuation token">)</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> <SPAN class="punctuation token">]</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> pt <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>FromWKT<SPAN class="punctuation token">(</SPAN><SPAN class="string token">'POINT(10.5 10.5)'</SPAN><SPAN class="punctuation token">,</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>SpatialReference<SPAN class="punctuation token">(</SPAN><SPAN class="number token">3857</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> 
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> fl<SPAN class="punctuation token">,</SPAN> <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>SelectLayerByLocation_management<SPAN class="punctuation token">(</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; polys<SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">'WITHIN_A_DISTANCE'</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; pt<SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="number token">1.0</SPAN><SPAN class="punctuation token">,</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN>&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="string token">'NEW_SELECTION'</SPAN>
<SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN><SPAN class="punctuation token">.</SPAN> <SPAN class="punctuation token">)</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator token">></SPAN> fl<SPAN class="punctuation token">.</SPAN>getSelectionSet<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="punctuation token">{</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">2</SPAN><SPAN class="punctuation token">}</SPAN>
<SPAN class="operator token">>></SPAN><SPAN class="operator 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>

The IDs returned from getSelectionSet() are offset by 1 from the index of the list passed to it.

In terms of ArcGIS Desktop/ArcMap, this new functionality was not supposed to be back ported, but it appears to have been back ported because I can run this code with 10.6.1. That said, I think using an ArcPy geometry as the select_features has worked for several versions.

DanPatterson_Retired

Select by Location according to the help files, requires 2 feature layers as input. (A hidden feature of arcpy? )

Zoltan, do you have the example where an arcpy geometry object and a feature layer made it work?

bixb0012

You didn't override it with None, you just created a point geometry without one. In your original code, I assume vert is an ArcPy Point, and ArcPy Points don't have spatial references so there is nothing to override.

ZoltanSzecsei

AAaargh! Yes! (Darren says: Try specifying a spatial reference on your PointGeometry)

I override it with 'None', so that negates my env setttings.

I've now added the Spatial Reference into the arcpy.PointGeometry() and it is now working.

Big thanks Darren 🙂

Zoltan

ZoltanSzecsei

I have - and have also set arcpy.env.outputCoordinateSystem = arcpy.SpatialReference(blah blah)

It is a projected system with metric coords

Dan's put me on to using geometry methods - now that I take a deeper look, it is a better way to solve my particular issue - I'm just busy wrapping my head around using disjoint/intersect/whatever...........

DarrenWiens2

Try specifying a spatial reference on your PointGeometry, even if it is unprojected.

ZoltanSzecsei

Hi Dan

I'll take a deeper look at the Point Geometry classes - but I did find plenty of examples using Geometry with SelectLayerByLocation, so I was wondering why my forehead has square indentations and my keyboard broken 🙂

Thanks and keep well,

Zoltan

DanPatterson_Retired

PointGeometry—ArcPy classes | ArcGIS Desktop intersects perhaps, if you want to work with arcpy geometry objects

The Select Layer By Location—Data Management toolbox | ArcGIS Desktop are using feature layers

Rejoindre la conversation

Les membres connectés peuvent publier, suivre les mises à jour, et plus encore. Nouveau ici ? Inscrivez-vous gratuitement.

Liens rapides

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