Des trous dans les polygones disparaissent après avoir effectué checkPolygonsHoles de maptools

Question

Note: Sur la suggestion de Edzer Pebesma, cette question a été crossposted à R-sig-géo, here, où il a reçu quelques bonnes réponses.

---

I rencontré le résultat inattendu suivant à l'aide checkPolygonsHoles:

# attach the worldmap as SpatialPolygonsDataFrame from the package maptools 
library(sp) 
library(maptools) 
data(wrld_simpl) 

# get a polygon with a hole 
shape_with_hole <- wrld_simpl[5,] 

# plot it (hole is left white, surrounded by blue color) 
plot(shape_with_hole, col = "blue") 

# perform checkPolygonsHoles 
shape_with_hole@polygons <- lapply(shape_with_hole@polygons, checkPolygonsHoles) 

# plot again, now holes aren't recognized as such 
plot(shape_with_hole, col = "blue") 

# and even the original SpatialPolygonsDataFrame object is changed !? 
plot(wrld_simpl[5,], col = "blue") 

Un effet secondaire irritant est que l'objet original wrld_simpl est également modifiée. Ce résultat me semble être un bug, ou ai-je raté quelque chose?

PS: l'objet shape_with_hole édité avec checkPolygonsHoles avant, continue à se comporter étrange:

# we check which polygons are marked as holes. The flags are still set 
# properly, although the `plot` function didn't recognize them: 
sapply(shape_with_hole@polygons[[1]]@Polygons, slot, "hole") 

[1] FALSE TRUE TRUE TRUE 

# load library rgdal for reprojection 
library(rgdal) 

# reproject with `spTransform`, just for testing 
shape_with_hole <- spTransform(shape_with_hole, 
CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")) 

# after reprojection all flags are set to FALSE 
sapply(shape_with_hole@polygons[[1]]@Polygons, slot, "hole") 

[1] FALSE FALSE FALSE FALSE 

Answer 1

Ce rapport clair a révélé un bogue dans le paquet R sp, qui est devenu manifeste que dans la R 3.1.x Voir here pour plus de détails. sp 1.1-0, disponible chez CRAN, a corrigé ce problème.