Nothing
cshp <- function(date=NA, useGW=TRUE) {
# check input
if (!is.na(date) && !inherits(date, "Date")) {
stop("date is not of type Date")
}
if (!is.na(date) && (date < as.Date("1946-1-1") | date > as.Date("2015-6-30"))) {
stop("Specified date is out of range")
}
# where to look for the dataset
path <- paste(system.file(package = "cshapes"), "shp/cshapes.shp", sep="/")
# load the dataset
cshp.full <- readShapePoly(path, proj4string=CRS("+proj=longlat +ellps=WGS84"), IDvar="FEATUREID")
if (is.na(date)) {
cshp.full
} else if (useGW) {
cshp.full <- cshp.full[cshp.full$GWCODE>=0,]
startdate <- as.Date(paste(cshp.full$GWSYEAR, cshp.full$GWSMONTH, cshp.full$GWSDAY, sep="-"))
enddate <- as.Date(paste(cshp.full$GWEYEAR, cshp.full$GWEMONTH, cshp.full$GWEDAY, sep="-"))
cshp.part <- cshp.full[startdate <= date & enddate >= date,]
cshp.part
} else {
cshp.full <- cshp.full[cshp.full$COWCODE>=0,]
startdate <- as.Date(paste(cshp.full$COWSYEAR, cshp.full$COWSMONTH, cshp.full$COWSDAY, sep="-"))
enddate <- as.Date(paste(cshp.full$COWEYEAR, cshp.full$COWEMONTH, cshp.full$COWEDAY, sep="-"))
cshp.part <- cshp.full[startdate <= date & enddate >= date,]
cshp.part
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.