plot.vector_cube: Plot RGB vector data cubes

In sits: Satellite Image Time Series Analysis for Earth Observation Data Cubes

Plot RGB vector data cubes

Description

Plot RGB raster cube

Usage

## S3 method for class 'vector_cube'
plot(
  x,
  ...,
  band = NULL,
  red = NULL,
  green = NULL,
  blue = NULL,
  tile = x[["tile"]][[1]],
  dates = NULL,
  seg_color = "yellow",
  line_width = 0.3,
  palette = "RdYlGn",
  rev = FALSE,
  scale = 1,
  first_quantile = 0.02,
  last_quantile = 0.98,
  max_cog_size = 1024
)

Arguments

x	Object of class "raster_cube".
...	Further specifications for plot.
band	Band for plotting grey images.
red	Band for red color.
green	Band for green color.
blue	Band for blue color.
tile	Tile to be plotted.
dates	Dates to be plotted.
seg_color	Color to show the segment boundaries
line_width	Line width to plot the segments boundary (in pixels)
palette	An RColorBrewer palette
rev	Reverse the color order in the palette?
scale	Scale to plot map (0.4 to 1.5)
first_quantile	First quantile for stretching images
last_quantile	Last quantile for stretching images
max_cog_size	Maximum size of COG overviews (lines or columns)

Value

A plot object with an RGB image or a B/W image on a color scale using the palette

Note

The following optional parameters are available to allow for detailed control over the plot output:

• graticules_labels_size: size of coordinates labels (default = 0.7)

• legend_title_size: relative size of legend title (default = 0.7)

• legend_text_size: relative size of legend text (default = 0.7)

• legend_bg_color: color of legend background (default = "white")

• legend_bg_alpha: legend opacity (default = 0.3)

Author(s)

Gilberto Camara, [e-mail]

Examples

if (sits_run_examples()) {
    # create a data cube from local files
    data_dir <- system.file("extdata/raster/mod13q1", package = "sits")
    cube <- sits_cube(
        source = "BDC",
        collection = "MOD13Q1-6.1",
        data_dir = data_dir
    )
    # Segment the cube
    segments <- sits_segment(
        cube = cube,
        output_dir = tempdir(),
        multicores = 2,
        memsize = 4
    )
    # plot NDVI band of the second date date of the data cube
    plot(segments, band = "NDVI", date = sits_timeline(cube)[1])
}

sits documentation built on Sept. 11, 2024, 6:36 p.m.