library(sp)
# Load the NGS data:
data(geochm)
coordinates(geochm) <- ~LONGITUDE+LATITUDE
proj4string(geochm) <- CRS("+proj=longlat +ellps=clrk66 +datum=NAD27 +no_defs")
## Not run:
# require(plotKML)
# data(SAGA_pal)
# # replace the missing values with half the detection limit:
# geochm$PB_ICP40 <- ifelse(geochm$PB_ICP40 < 0, 2, geochm$PB_ICP40)
# shape = "http://maps.google.com/mapfiles/kml/pal2/icon18.png"
# kml(geochm, shape = shape, colour = log1p(PB_ICP40), labels = "",
# colour_scale = SAGA_pal[[1]], kmz = TRUE)
# ## End(Not run)
Run the code above in your browser using DataLab