3

I have several shp files which are the components (or tiles) of a map, I have joint all these tiles together to have a large shp file. Then I realized some of the polygon is now split into 2 halves, as the upper part belongs to a tile, while the lower part belongs to another tile.

This output disturbs me a lot as I need to fill some polygon, but not part of the "split polygon", so I need to find a way to merge these split polygon back to one (preferably mechanical way, as I have 10000+ polygons, not knowing which are split, which are not).

I need to know if the two or more polygons originate from the larger polygon. And I think I can somehow use union and spRbind to create an algorithm for this task. But I doubt if I am really the first one who encounters this problem, so I wonder if there is any existing tool in R that can solve this.

Thanks.

1
  • Do you have some sample data you could post? What do you mean by 'not knowing which are split' - is there no ID column that could be used to join as @yellowcap suggests?
    – Simbamangu
    Commented Mar 20, 2012 at 12:11

1 Answer 1

4

I can only provide a partial answer to your question.

There is a function in the maptools R package called unionSpatialPolygons, see here for a manual.

The function expects an ID field, and all polygons with the same ID value will be merged. If you have an ID that identifies your polygons to be dissolved, use it. Otherwise you could try something like

unifiedPolygons <- unionSpatialPolygons(myPolygons, rep(1, length(myPolygons))

but that will probably also dissolve all other polygons that touch within you area.

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.