Tag Info

New answers tagged

0

This can be done with the spatial join tool. You need to select the JOIN_ONE_TO_ONE option. Of course, you might have to run this tool several times if your stressors are not in the same feature class (or you first merge them if they are all of the same type) JOIN_ONE_TO_ONE —If multiple join features are found that have the same spatial relationship ...


0

I found some code here Is there a python option to "join attributes by location"? that uses Shapely and Fiona that outputted the points and then rearranged the code to output the polygon instead. This seems to be working on one of the 5 grid point buffer geographies and running through the ~80,000 grid points. import fiona from shapely.geometry ...


2

You don't have one polygon, you have many. It looks like the xy coords from all polygons are just dumped into one large list to make a single polygon. Hence the polygon boundary is jumping around creating bad topology. As the above commenters say, your question is poor. But, taking a guess: if you only want one polygon merge your data together prior to ...


0

The script below assumes there are 3 feature layers: Point layer - nodes or street junctions Polyline layer - street segments. Associated table has 2 fields of integer type to store sequential number of from and to node. Polyline layer called 'connections' to store results. Only essential part of the script shown below, i.e. no error handling, parameters ...


2

To remove polygons that don't have an intersection with lines: Open ArcMap and add the polygons and lines to the table of contents. Locate the Select by Locations dialog and select the ones that do intersect. Switch the selection by going to the selection in the table of contents Right click on the layer name and select 'switch selection'. Start ...


1

In your solution you mix geometries and predicates of ogr and Shapely With Shapely only (intersection(), intersects()) from shapely.geometry import Point, LineString, mapping geom1 = LineString([(0, 0), (0, 1)]) geom2 = LineString([(1, 1), (-1,-1)]) intersection = geom1.intersection(geom2) # geo_interface -> GeoJSON mapping(intersection) {'type': ...



Top 50 recent answers are included