One polygon in the longitude and latitude format (see
proj4string=CRS("+proj=longlat +ellps=WGS84")
).
1 2 3 4 5 6 7 8 9 10 | library(maptools)
validTractShp <- tractShp[!is.na(tractShp$BUYPOW), ] # Remove 2 tracts with NA's
plot(tractShp, col="white", border="white", axes=TRUE,
main="Dallas Census Tracts with Food Deserts")
plot(validTractShp, col="ivory2", border="white", add=TRUE)
plot(lakesShp, col="skyblue", border="skyblue",add=TRUE)
plot(hwyShp, col="cornsilk3", lwd=3, add=TRUE)
plot(foodDesertShp, border="magenta",lwd=2, add=TRUE)
plot(bndShp, border="black", add=TRUE)
box()
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.