plot.uncertainty_vector_cube: Plot uncertainty vector cubes

In e-sensing/sits: Satellite Image Time Series Analysis for Earth Observation Data Cubes

Plot uncertainty vector cubes

Description

plots a probability cube using stars

Usage

## S3 method for class 'uncertainty_vector_cube'
plot(
  x,
  ...,
  tile = x$tile[[1]],
  palette = "RdYlGn",
  rev = TRUE,
  tmap_options = NULL
)

Arguments

x	Object of class "probs_vector_cube".
...	Further specifications for plot.
tile	Tile to be plotted.
palette	RColorBrewer palette
rev	Reverse order of colors in palette?
tmap_options	Named list with optional tmap parameters max_cells (default: 1e+06) scale (default: 1.0) graticules_labels_size (default: 0.7) legend_title_size (default: 1.0) legend_text_size (default: 1.0) legend_bg_color (default: "white") legend_bg_alpha (default: 0.5)

Value

A plot containing probabilities associated to each class for each pixel.

Author(s)

Gilberto Camara, [e-mail]

Examples

if (sits_run_examples()) {
    # create a random forest model
    rfor_model <- sits_train(samples_modis_ndvi, sits_rfor())
    # create a data cube from local files
    data_dir <- system.file("extdata/raster/mod13q1", package = "sits")
    cube <- sits_cube(
        source = "BDC",
        collection = "MOD13Q1-6",
        data_dir = data_dir
    )
    # segment the image
    segments <- sits_segment(
        cube = cube,
        seg_fn = sits_slic(step = 5,
                           compactness = 1,
                           dist_fun = "euclidean",
                           avg_fun = "median",
                           iter = 20,
                           minarea = 10,
                           verbose = FALSE),
        output_dir = tempdir()
    )
    # classify a data cube
    probs_vector_cube <- sits_classify(
        data = segments,
        ml_model = rfor_model,
        output_dir = tempdir()
    )
    # measure uncertainty
    uncert_vector_cube <- sits_uncertainty(
        cube = probs_vector_cube,
        type = "margin",
        output_dir = tempdir()
    )
    # plot the resulting uncertainty cube
    plot(uncert_vector_cube)
}

e-sensing/sits documentation built on Jan. 28, 2024, 6:05 a.m.