# soil data
data("worldMap")
worldMap = terra::unwrap(worldMap)
oldpar=map.new(worldMap)
plot(worldMap, border='red', lwd=3, add=TRUE)
plot(worldMap[worldMap$NAME == 'Brazil',],
add=TRUE, col='green')
par(oldpar)
Run the code above in your browser using DataLab