Finding Ellipses in a polygon feature class?

Steve_Salas

I recently had difficulty running the geoprocessing tool Find Identical (arcpy.FindIdentical_management) in ArcMap - it was crashing with the unhelpful "999999: Error executing function."  I eventually eliminated possibilities down to two very small ellipse polygons (as opposed to our typical multi-vertex polygons).  Once i eliminated these two polygons, the tool ran fine on the remaining 200,000+ features.

I didn't find anything particularly helpful in the ST_Geometry functions (Oracle) that would help me query for an ellipse polygon next time I need to determine if I have some in my data.  Even sde.st_astext(shape) spits out a significant list of coordinate pairs that I would guess would be the ellipse if generalized to multiple vertices.

A rather clumsy method in geoprocessing is to run the Feature Vertices to Points tool ("All" option) and summarize the output on the original feature identifier - ellipses appear to return only two output points (start and end point) opposed to the 3 or more points from a valid polygon composed of vertices.

Is there a better way to select records that use ellipse (and circle) types for their geometry?  I'm not a ArcObjects programmer but that might be an alternative there if not in SQL or ArcPy?  Thanks for any ideas!    -S.

JoshuaBixby

Although Esri's Geometry objects—Common Data Types | ArcGIS for Developers support circular arcs, elliptical arcs, and Bezier curves, there is a catch when it comes to enterprise geodatabases.  Esri's curved shapes are not always stored natively in the spatial data column.  Depending on the DBMS and spatial data type, Esri's curved shapes are stored in a separate field in a separate table, and what is stored in the base table's spatial data column is a linear approximation of the curved shape.  This is why your sde.ST_AsText is returning the long list of coordinates you are seeing.

Given the variablitiy of how Esri's curved shapes are handled in different DMBSs and spatial data types, I recommend using an ArcPy cursor to retrieve the true Esri shape and then test it for circular arcs, elliptical arcs, and Bezier curves.

Something like the following should work:

fc <SPAN class="operator token">=</SPAN> <SPAN class="comment token"># feature class to check for curved polygons</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>fc<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="string token">"OID@"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"SHAPE@JSON"</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> cur<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; oid_list <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>oid <SPAN class="keyword token">for</SPAN> oid<SPAN class="punctuation token">,</SPAN>json <SPAN class="keyword token">in</SPAN> cur <SPAN class="keyword token">if</SPAN> json <SPAN class="operator token">and</SPAN> <SPAN class="string token">"curveRing"</SPAN> <SPAN class="keyword token">in</SPAN> json<SPAN class="punctuation token">]</SPAN>
<SPAN class="keyword token">print</SPAN><SPAN class="punctuation token">(</SPAN>oid_list<SPAN class="punctuation token">)</SPAN><SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

If oid_list is empty, there are no curved polygons.