Creates or adds a stars object based on the input object or a set of parameters.
It accepts either an object of class stars
or lsp
.
In the first case, the output is created based on
the window
parameter.
In the second case, the output converts the lsp
object into
a stars
object.
lsp_add_stars(x = NULL, window = NULL, metadata = TRUE)
# Default S3 method
lsp_add_stars(x = NULL, window = NULL, metadata = TRUE)
# S3 method for class 'lsp'
lsp_add_stars(x = NULL, window = NULL, metadata = TRUE)
Object of class stars
or lsp
.
For stars
, window
or window_size
can be used.
Specifies areas for analysis. It can be either: NULL
, a numeric value, or an sf
object. If window=NULL
calculations are performed for a whole area. If the window
argument is numeric, it is a length of the side of a square-shaped block of cells. Expressed in the numbers of cells, it defines the extent of a local pattern. If an sf
object is provided, each feature (row) defines the extent of a local pattern. The sf
object should have one attribute (otherwise, the first attribute is used as an id).
Logical. Only when x`` is of class
lsp. If
TRUE, the output object will have metadata ("id" and "na_prop"). If
FALSE`, the output object will not have metadata ("id" and "na_prop").
A stars
object converted from the input object or a provided set of parameters
library(stars)
landform = read_stars(system.file("raster/landforms.tif", package = "motif"))
plot(landform)
landform_lsp = lsp_add_stars(landform, window = 100)
plot(landform_lsp)
lc_cove = lsp_signature(landform, type = "cove", window = 200, normalization = "pdf")
lc_cove_lsp = lsp_add_stars(lc_cove)
plot(lc_cove_lsp)
plot(lc_cove_lsp["na_prop"])
# \donttest{
# larger data example
library(stars)
landform = read_stars(system.file("raster/landform.tif", package = "motif"))
plot(landform)
#> downsample set to 10
#> Error in image.stars(x, ..., axes = axes, breaks = breaks, col = col, key.pos = key.pos, key.width = key.width, key.length = key.length, main = NULL): rgb plotting not supported on this device
landform_lsp = lsp_add_stars(landform, window = 100)
plot(landform_lsp)
lc_cove = lsp_signature(landform, type = "cove", window = 200, normalization = "pdf")
lc_cove_lsp = lsp_add_stars(lc_cove)
plot(lc_cove_lsp)
plot(lc_cove_lsp["na_prop"])
# }