Merge single linear feature with multiple polygons

Question

Hello,

The picture below is a single polygon. One of the polygons within this layer is a road (highlighted in blue) which runs through the middle of other polygons within that layer. I want the road polygon to be merged with the adjacent polygons, or alternatively, deleted, and the gap where the road was to be filled with the adjacent polygons.

It would be simple to do manually, but there are hundreds of polygons impacted and I'm wondering if there is a better way to do this.

Thanks,

Ben

Anonymous User · Answer

Hi Ben:

A colleague recommended the Eliminate GP tool, which might be an option.

Can you post a representative sample of your dataset ? (please make sure its anonymized so that no personal information is present in the attributes). There might be a way to avoid altering polygon coordinates outside of the road area, but I'd like to test that first.

-jt

Anonymous User · Answer

Hi Ben:

Here's one approach.

1.Select all roads from your input polygons.

2. reverse the selection (so only non-road polygons are selected)

3. execute polygon to raster gp tool, using an attribute that uniquely identifies each non-road polygon as the 'value_field' parameter. We'll call this output raster "PolyR". Here's an example of what "PolyR" might look like (the white is nodata):

4. use that as the input to the Euclidean Allocation tool.

5. The output should look something like this. We'll call this 'ErasedRoadsR":

6. In ErasedRoadsR, the cell values from PolyR have been 'spread' to the nearest NODATA cells.

Here are "PolyR" and "EraseRoadsR" shown together, using transparency for the latter. in EraseRoadsR each NODATA cell from POlyR is assigned the value of its closest DATA neighbor, thus 'filling in' the NODATA gaps.

7. Finally, you can use the RasterToPolygon GP tool to convert back to a polygon feature class. You can use the 'Simply' parameter to get rid of the jaggies.

Hope this helps,

-jt

BenPearse · Answer

Hi Jim,

No my road polygons don't overlap so the union wouldn't work.

The coordinates of my vertices are not sacred, although I would need to convert them back to polygons that are fairly similar to my original vector layer.

I would definitely be interested in hearing your approach, the only caveat is that I don't currently have the spatial analyst extension. I'm pushing my supervisor to get it though, so maybe this could add to my argument for it.

Cheers,

Ben