plot.raster_cube: Plot RGB data cubes

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

Plot RGB data cubes

Description

Plot RGB raster cube

Usage

## S3 method for class 'raster_cube'
plot(
  x,
  ...,
  band = NULL,
  red = NULL,
  green = NULL,
  blue = NULL,
  tile = x[["tile"]][[1]],
  dates = NULL,
  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.
palette	An RColorBrewer palette
rev	Reverse the color order in the palette?
scale	Scale to plot map (0.4 to 1.0)
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

Note

Use scale parameter for general output control. The dates parameter indicates the date allows plotting of different dates when a single band and three dates are provided, 'sits' will plot a multi-temporal RGB image for a single band (useful in the case of SAR data). For RGB bands with multi-dates, multiple plots will be produced.

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
    )
    # plot NDVI band of the second date date of the data cube
    plot(cube, band = "NDVI", dates = sits_timeline(cube)[1])
}

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