Skip to content

Problem with polygon holes and islands

2 messages · Ian Reeve, Roger Bivand

Ian Reeve
# 
Dear all,

I'm importing census tract boundaries into R from a government website where the boundaries are provided as Mapinfo Interchange Format files, using:

mypolygons <- readOGR("websitefile.mif", layer="websitefile")

Some census tracts have holes.  In these cases, the hole slot of the outer polygon is TRUE and the hole slot of the inner polygon is FALSE.  I have to reverse these manually to get these tracts to plot correctly.

Some census tracts comprise a collection of islands.  In these cases, the hole slot of the first polygon in the collection is FALSE, and in the remainder of the polygons in the collection it is TRUE.  

Again I can correct these manually (and tediously), but I'm wondering if there is either a way of preventing this happening when importing the .mif file, or a single method that can be applied to a SpatialPolygons object to fix the problems in one go.

Any leads much appreciated.

Ian Reeve
University of New England
Armidale, NSW
Roger Bivand
#
On Sat, 5 Mar 2011, Ian Reeve wrote:

            

      
      
      
    

        
If you could make a sample file available, it might be possible to fix. 
The hole slots are being assigned based on the ring direction of the 
polygon boundaries, if the files do not specify the interior or exterior 
status of the ring. This can also be fixed by running checkPolygonsHoles() 
in the maptools package on the list in the polygons slot of the 
SpatialPolygonsDataFrame object, but there may be other possibilities.

Roger