POST
|
Hi Shaun, Thanks for the pointer - will try it out. For the lines I am now using the resultant unsplit lines to select lines that "SHARE_A_LINE_SEGMENT_WITH" and then deleting that set. I'll try the Symmetrical Difference for the dissolved polygon "left-behinds" Not sure why ESRI chose to complicate what surely must be a 'standard' need. Regards, Zoltan
... View more
11-01-2023
09:52 PM
|
0
|
0
|
419
|
POST
|
Hi, I'm using arcpy and ArcGIS Pro 3.1.3 I have a Lines feature class with 'millions' of line features and I need to use arcpy.management.UnsplitLine() on it. The trouble is, is that this outputs the unsplit lines to another (new) FeatureClass rather than replaces the unsplit lines. How can I locate and copy from the original FeatureClass to this new FeatureClass, all the lines that were NOT 'unsplit' I have a similar problem with dissolving polygons. Thanks in advance, Zoltan
... View more
11-01-2023
04:56 AM
|
0
|
2
|
506
|
POST
|
Hi, I have over 2700 FileGDBs in which to count feature occurrences based on an attribute value. The attribute field (FeaType) has over 200 possible values. (Typically: Tree, Road, Building, etc etc.) I need to count the occurences of each feature type, per FileGDB and put them in a table with 1 row per FileGDB and 1 column per FeaType value. Putting them into a table is no issue - I am using xlsxwriter. What would be the most efficient way to count these features? I have tried iterating through each FileGDB, each FeatureClass, row by row incrementing a table entry based on the FeaType value - very slow 😞 I could try iterating through each FeaType value, then using that value to 'select' and 'getcount'. But surely there is a more efficient way? Some pointers would be great. Thanks in advance, Zoltan
... View more
05-02-2023
05:43 AM
|
0
|
3
|
719
|
POST
|
Actually, XY snap is all I need for this exercise as it is also interesting to me if the vertices are same XY but Z is different I think you've hit the nail on the head for me. My Z tolerance is not as tight as my XY tolerance, so that is why the command below is/was erratic in its hits: pg = arcpy.PointGeometry(vert, crs, True) arcpy.SelectLayerByLocation_management("hgt_feats","WITHIN_A_DISTANCE_3D",pg, 0.050, 'NEW_SELECTION') I've just removed the '_3D' and I am getting many more hits now. I've also examined some test features and the Z values do drift >5cm as I had in my search radius. Let me work with this, because I think you've just solved it for me. Big thanks, and chat soon (I'll report back on this list if the SelectLayer as used above still seems dodgy). Regards, Zoltan
... View more
01-12-2019
10:30 AM
|
0
|
0
|
1070
|
POST
|
Argis 10.4 Using Pycharm IDE and running it in the IDE
... View more
01-12-2019
08:22 AM
|
0
|
0
|
1070
|
POST
|
Hi, I'm looking for an efficient way to see which vertices in one feature touch a vertex in another feature. I will then need to edit the vertex of the "first" feature. Problem: I have a large polygon feature set that I will select some polygons from. For each polygon in this subset, I need to check if any of its vertices coincide with that of another (Polygon or polyline) feature in (two) other feature classes. If a vertex in this first subset touches a vertex in the 2nd or 3rd feature class, I will have to edit that vertex and update that polygon (in the first feature set). I have tried crawling through the vertices of the subset and using SelectLayerByLocation by coverting the XYZ of the vertex for use in the Select statement, but this yields unreliable results. I have tried using various SelectLayerByLocation with various "Intersect" or "WITHIN A DISTANCE 3D" cluses but the problem here is that I cannot (find a way to) release my latest selection without clearing the full selection (because I work with sub-selections from all 3 feature classes) I'm wondering if the way to go is to create a point file from the first subset, adding the vertex's OID and sequence information as attributes to the resulting point, and then JOINing by location the points to the 2nd and 3rd feature sets. Where there is a JOIN, they touch, where not, they don't touch. Then using this information to delete and manually rebuild the polygons in this first subset. All this complication is because I cannot find a way to see if anything exists at a particular XYZ location. Comments/thoughts welcome. Thanks in advance, Zoltan
... View more
01-12-2019
06:58 AM
|
0
|
4
|
1287
|
POST
|
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"
... View more
01-12-2019
06:37 AM
|
0
|
1
|
712
|
POST
|
No - this was my first attempt - but see my post with screenshot, above - it does work.
... View more
01-11-2019
09:53 PM
|
0
|
2
|
712
|
POST
|
Actually (psychologically) I did "override" it. If one sets and env to a SRS/CRS, one should be able to assume that all dimension and coordinate related values would be interpreted by ArcGIS to be in the CRS set by the env. By placing 'None' in a specific command, that (intuitively) explicitly says "override the env requested CRS" I know this is about symantecs of speach, but I make this point because when ESRI provide certain settings, certain expectations are also "enabled". HTH Z
... View more
01-11-2019
09:50 PM
|
0
|
1
|
1533
|
POST
|
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
... View more
01-11-2019
12:23 PM
|
2
|
3
|
1533
|
POST
|
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...........
... View more
01-11-2019
12:15 PM
|
0
|
4
|
1533
|
POST
|
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
... View more
01-11-2019
11:59 AM
|
0
|
0
|
1533
|
POST
|
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
... View more
01-11-2019
11:36 AM
|
0
|
19
|
2641
|
POST
|
Hi, Is this the correct forum to ask about arcpy.SelectLayerByLocation_management using XYZ coordinates? (If not, please send me a link) Thanks, Zoltan
... View more
01-11-2019
11:15 AM
|
0
|
2
|
535
|
Title | Kudos | Posted |
---|---|---|
2 | 01-11-2019 12:23 PM |
Online Status |
Offline
|
Date Last Visited |
11-02-2023
07:36 AM
|