I am trying to measure the distance between two points using the arcgis.geometry.distance function. However, I get this error: RuntimeError: Unable to complete operation. The operation cannot be performed on a non-simple geometry. (Error Code: 400) My code is: # Calculate distance for index , values in cc_df . iterrows ( ) : cap_point = values [ 'SHAPE' ] # A point object print ( type ( cap_point ) , " geometry 1 type" ) cap_stid = values [ 'Study_ID' ] # (string) for c in dist_clm : # A list of columns moi = c [ 5 : ] # The specific column moi = moi [ 5 : ] # A part of the column name used to acess the 'SHAPE' value in a dictionary r_shp = rel_dict [ cap_stid ] [ moi ] [ 'shp' ] # the 'SHAPE' column of a data frame print ( type ( r_shp ) , " geometry 2 type" ) # Printing the data type as part of testing r_wkid = rel_dict [ cap_stid ] [ moi ] [ 'shp' ] [ 'spatialReference' ] [ 'wkid' ] # The WKID (spatial refrence) d = arcgis . geometry . distance ( r_wkid , cap_point , r_shp , "METERS" , False ) ) It is basically extracting a point geometry from a dictionary and measuring the distance to another point geometry in a SPED. (Or that is what it is supposed to do) The output is: geometry 1 type geometry 2 type Unable to complete operation. The operation cannot be performed on a non-simple geometry. Are point objects not simple geometry? The documentation on this is very sparse and I can't find any examples. Here is a link to the documentation: arcgis.geometry module — arcgis 1.6.0 documentation Any insight would be helpful. Thanks arcgis 1.5.3 pandas 0.24.2

Hi Ben, For info on what Distance Units should be refer to this page: Distance—ArcGIS REST API: Services Directory | ArcGIS for Developers which branches off into: * esriSRUnitType Constants * esriSRUnit2Type Constants It seems the function fails when you supply the distance unit directly (this might constitute a defect). However, it works when you store it in a variable like so: from arcgis import geometry sr = 102100 distance_unit = 9001 geom1 = { 'x' : - 7703308.762905646 , 'y' : - 2937284.4374484867 , 'spatialReference' : { 'wkid' : 102100 } } geom2 = { 'x' : 13614373.72405075 , 'y' : 2571663.047185693 , 'spatialReference' : { 'wkid' : 102100 } } pt1 = geometry . Point ( geom1 ) pt2 = geometry . Point ( geom2 ) d = arcgis . geometry . distance ( sr , pt1 , pt2 , distance_unit , False ) print ( d ) The result is: {'distance': 22017994.663527936} Kind regards, Earl

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

How to use geometry.distance in Python API

BenjaminSperry1

I am trying to measure the distance between two points using the arcgis.geometry.distance function.

However, I get this error:

RuntimeError: Unable to complete operation. The operation cannot be performed on a non-simple geometry. (Error Code: 400)

My code is:

<SPAN class="comment token"># Calculate distance</SPAN>
<SPAN class="keyword token">for</SPAN> index<SPAN class="punctuation token">,</SPAN> values <SPAN class="keyword token">in</SPAN> cc_df<SPAN class="punctuation token">.</SPAN>iterrows<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
    cap_point <SPAN class="operator token">=</SPAN> values<SPAN class="punctuation token">[</SPAN><SPAN class="string token">'SHAPE'</SPAN><SPAN class="punctuation token">]</SPAN>  <SPAN class="comment token"># A point object</SPAN>
    <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>type<SPAN class="punctuation token">(</SPAN>cap_point<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">" geometry 1 type"</SPAN><SPAN class="punctuation token">)</SPAN>
    cap_stid <SPAN class="operator token">=</SPAN> values<SPAN class="punctuation token">[</SPAN><SPAN class="string token">'Study_ID'</SPAN><SPAN class="punctuation token">]</SPAN>  <SPAN class="comment token"># (string)</SPAN>

    <SPAN class="keyword token">for</SPAN> c <SPAN class="keyword token">in</SPAN> dist_clm<SPAN class="punctuation token">:</SPAN>  <SPAN class="comment token"># A list of columns</SPAN>
        moi <SPAN class="operator token">=</SPAN> c<SPAN class="punctuation token">[</SPAN><SPAN class="number token">5</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="punctuation token">]</SPAN>  <SPAN class="comment token"># The specific column</SPAN>
        moi <SPAN class="operator token">=</SPAN> moi<SPAN class="punctuation token">[</SPAN><SPAN class="number token">5</SPAN><SPAN class="punctuation token">:</SPAN><SPAN class="punctuation token">]</SPAN>  <SPAN class="comment token"># A part of the column name used to acess the 'SHAPE' value in a dictionary</SPAN>
        r_shp <SPAN class="operator token">=</SPAN> rel_dict<SPAN class="punctuation token">[</SPAN>cap_stid<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN>moi<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'shp'</SPAN><SPAN class="punctuation token">]</SPAN>  <SPAN class="comment token"># the 'SHAPE' column of a data frame</SPAN>
        <SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>type<SPAN class="punctuation token">(</SPAN>r_shp<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">" geometry 2 type"</SPAN><SPAN class="punctuation token">)</SPAN>  <SPAN class="comment token"># Printing the data type as part of testing</SPAN>
        r_wkid <SPAN class="operator token">=</SPAN> rel_dict<SPAN class="punctuation token">[</SPAN>cap_stid<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN>moi<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'shp'</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'spatialReference'</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="string token">'wkid'</SPAN><SPAN class="punctuation token">]</SPAN>  <SPAN class="comment token"># The WKID (spatial refrence)</SPAN>

        d <SPAN class="operator token">=</SPAN> arcgis<SPAN class="punctuation token">.</SPAN>geometry<SPAN class="punctuation token">.</SPAN>distance<SPAN class="punctuation token">(</SPAN>r_wkid<SPAN class="punctuation token">,</SPAN> cap_point<SPAN class="punctuation token">,</SPAN> r_shp<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"METERS"</SPAN><SPAN class="punctuation token">,</SPAN>
                                     <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN><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>

It is basically extracting a point geometry from a dictionary and measuring the distance to another point geometry in a SPED. (Or that is what it is supposed to do)

The output is:

<class 'arcgis.geometry._types.Point'>  geometry 1 type <class 'arcgis.geometry._types.Point'>  geometry 2 type

Unable to complete operation. The operation cannot be performed on a non-simple geometry.

Are point objects not simple geometry?

The documentation on this is very sparse and I can't find any examples.

Here is a link to the documentation: arcgis.geometry module — arcgis 1.6.0 documentation

Any insight would be helpful. Thanks

arcgis 1.5.3

pandas 0.24.2

Encontre mais publicações marcadas com

Accepted answers

EarlMedina

Hi Ben,

For info on what Distance Units should be refer to this page: Distance—ArcGIS REST API: Services Directory | ArcGIS for Developers which branches off into:

It seems the function fails when you supply the distance unit directly (this might constitute a defect). However, it works when you store it in a variable like so:

<SPAN class="keyword token">from</SPAN> arcgis <SPAN class="keyword token">import</SPAN> geometry

sr <SPAN class="operator token">=</SPAN> <SPAN class="number token">102100</SPAN>
distance_unit <SPAN class="operator token">=</SPAN> <SPAN class="number token">9001</SPAN>

geom1 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="string token">'x'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="number token">7703308.762905646</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'y'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="number token">2937284.4374484867</SPAN><SPAN class="punctuation token">,</SPAN> 
         <SPAN class="string token">'spatialReference'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="string token">'wkid'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="number token">102100</SPAN><SPAN class="punctuation token">}</SPAN><SPAN class="punctuation token">}</SPAN>
geom2 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="string token">'x'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="number token">13614373.72405075</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'y'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="number token">2571663.047185693</SPAN><SPAN class="punctuation token">,</SPAN> 
         <SPAN class="string token">'spatialReference'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="string token">'wkid'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="number token">102100</SPAN><SPAN class="punctuation token">}</SPAN><SPAN class="punctuation token">}</SPAN>

pt1 <SPAN class="operator token">=</SPAN> geometry<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>geom1<SPAN class="punctuation token">)</SPAN>
pt2 <SPAN class="operator token">=</SPAN> geometry<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>geom2<SPAN class="punctuation token">)</SPAN>

d <SPAN class="operator token">=</SPAN> arcgis<SPAN class="punctuation token">.</SPAN>geometry<SPAN class="punctuation token">.</SPAN>distance<SPAN class="punctuation token">(</SPAN>sr<SPAN class="punctuation token">,</SPAN> pt1<SPAN class="punctuation token">,</SPAN> pt2<SPAN class="punctuation token">,</SPAN> distance_unit<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>d<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>

The result is:

{'distance': 22017994.663527936}

Kind regards,

Earl

All comments

EarlMedina

Hi Ben,

For info on what Distance Units should be refer to this page: Distance—ArcGIS REST API: Services Directory | ArcGIS for Developers which branches off into:

It seems the function fails when you supply the distance unit directly (this might constitute a defect). However, it works when you store it in a variable like so:

<SPAN class="keyword token">from</SPAN> arcgis <SPAN class="keyword token">import</SPAN> geometry

sr <SPAN class="operator token">=</SPAN> <SPAN class="number token">102100</SPAN>
distance_unit <SPAN class="operator token">=</SPAN> <SPAN class="number token">9001</SPAN>

geom1 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="string token">'x'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="number token">7703308.762905646</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'y'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="number token">2937284.4374484867</SPAN><SPAN class="punctuation token">,</SPAN> 
         <SPAN class="string token">'spatialReference'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="string token">'wkid'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="number token">102100</SPAN><SPAN class="punctuation token">}</SPAN><SPAN class="punctuation token">}</SPAN>
geom2 <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="string token">'x'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="number token">13614373.72405075</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'y'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="number token">2571663.047185693</SPAN><SPAN class="punctuation token">,</SPAN> 
         <SPAN class="string token">'spatialReference'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="punctuation token">{</SPAN><SPAN class="string token">'wkid'</SPAN><SPAN class="punctuation token">:</SPAN> <SPAN class="number token">102100</SPAN><SPAN class="punctuation token">}</SPAN><SPAN class="punctuation token">}</SPAN>

pt1 <SPAN class="operator token">=</SPAN> geometry<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>geom1<SPAN class="punctuation token">)</SPAN>
pt2 <SPAN class="operator token">=</SPAN> geometry<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>geom2<SPAN class="punctuation token">)</SPAN>

d <SPAN class="operator token">=</SPAN> arcgis<SPAN class="punctuation token">.</SPAN>geometry<SPAN class="punctuation token">.</SPAN>distance<SPAN class="punctuation token">(</SPAN>sr<SPAN class="punctuation token">,</SPAN> pt1<SPAN class="punctuation token">,</SPAN> pt2<SPAN class="punctuation token">,</SPAN> distance_unit<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>d<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>

The result is:

{'distance': 22017994.663527936}

Kind regards,

Earl

BenjaminSperry1

Thanks for the speedy response Dylan.

It is a point feature so only the "Null Geometry" would apply. However, in exploring that possibility I found that deleting the distance units got rid of the error. I set geodesic to True to get meters.

arcgis.geometry.distance(r_wkid, cap_point, r_shp, geodesic=True)<SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

All the units I tried caused the error so it must be something to do with that element.

geodesic will work find for my current application but it would be nice to know why planar with units defined didn't work.

Thanks,

Ben

DylanHarwell

Not sure if this is helpful, but would your data fall into any of these categories?

https://support.esri.com/en/technical-article/000007177

Participe da Conversa

Membros conectados podem postar, seguir atualizações e mais. Novo aqui? Registre uma conta gratuita.

Links Rápidos

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