Nothing
classify.map <-
function(map, suit.classes, output.name = NULL, load.map = FALSE) {
if (class(suit.classes) != "CBI")
stop("The function classify.map needs an object of class 'CBI'!")
if (is.null(output.name))
output.name <- paste(map, "_class", sep = "")
breaks <- suit.classes$intervals[2:4]
print(breaks)
settings_text <- " 0 thru x1x = 1 unsuitable\n x2x thru x3x = 2 marginal\n x4x thru x5x = 3 suitable\n x6x thru 1 = 4 optimal"
settings_text <- sub("x1x", breaks[1], settings_text)
settings_text <- sub("x2x", breaks[1], settings_text)
settings_text <- sub("x3x", breaks[2], settings_text)
settings_text <- sub("x4x", breaks[2], settings_text)
settings_text <- sub("x5x", breaks[3], settings_text)
settings_text <- sub("x6x", breaks[3], settings_text)
write(settings_text, file = "classes.txt")
execGRASS("r.reclass", input = map, output = output.name, rules = "classes.txt", flags = c("overwrite", "quiet"))
cat(paste(output.name, " raster map was sucessfully classify!\n", sep = ""))
if (load.map) {
map <- list()
map[[output.name]] <- raster(readRAST(output.name))
class(map) <- c("hsm", "classified")
return(map)
}
}
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.