R: Plot a landsat Object

plot,landsat-method {oce} | R Documentation

## Plot a landsat Object

### Description

Plot the data within a landsat image, or information computed from the data. The second category includes possibilities such as an estimate of surface temperature and the `"terralook"` estimate of a natural-color view.

### Usage

    ## S4 method for signature 'landsat'
    plot(
      x,
      band,
      which = 1,
      decimate = TRUE,
      zlim,
      utm = FALSE,
      col = oce.colorsPalette,
      drawPalette = TRUE,
      showBandName = TRUE,
      alpha.f = 1,
      red.f = 1.7,
      green.f = 1.5,
      blue.f = 6,
      offset = c(0, -0.05, -0.2, 0),
      transform = diag(c(red.f, green.f, blue.f, alpha.f)),
      debug = getOption("oceDebug"),
      ...
    )

### Arguments

`x` | a landsat object.
`band` | If given, the name of the band. For Landsat-8 data, this may be one of: `"aerosol"`, `"blue"`, `"green"`, `"red"`, `"nir"`, `"swir1"`, `"swir2"`, `"panchromatic"`, `"cirrus"`, `"tirs1"`, or `"tirs2"`. For Landsat-7 data, this may be one of `"blue"`, `"green"`, `"red"`, `"nir"`, `"swir1"`, `"tirs1"`, `"tirs2"`, `"swir2"`, or `"panchromatic"`. For Landsat data prior to Landsat-7, this may be one of `"blue"`, `"green"`, `"red"`, `"nir"`, `"swir1"`, `"tirs1"`, `"tirs2"`, or `"swir2"`. If `band` is not given, the (`"tirs1"`) will be used if it exists in the object data, or otherwise the first band will be used. In addition to the above, using `band="temperature"` will plot an estimate of at-satellite brightness temperature, computed from the `tirs1` band, and `band="terralook"` will plot a sort of natural color by combining the `red`, `green`, `blue` and `nir` bands according to the formula provided at `https://lta.cr.usgs.gov/terralook/what_is_terralook` (a website that worked once, but failed as of Feb 2, 2017).
`which` | Desired plot type; 1=image, 2=histogram.
`decimate` | An indication of the desired decimation, passed to `imagep()` for image plots. The default yields faster plotting. Some decimation is sensible for full-size images, since no graphical displays can show 16 thousand pixels on a side.
`zlim` | Either a pair of numbers giving the limits for the colorscale, or `"histogram"` to have a flattened histogram (i.e. to maximally increase contrast throughout the domain.) If not given, the 1 and 99 percent quantiles are calculated and used as limits.
`utm` | A logical value indicating whether to use UTS (easting and northing) instead of longitude and latitude on plot.
`col` | Either a function yielding colors, taking a single integer argument with the desired number of colors, or the string `"natural"`, which combines the information in the `red`, `green` and `blue` bands and produces a natural-hue image. In the latter case, the band designation is ignored, and the object must contain the three color bands.
`drawPalette` | Indication of the type of palette to draw, if any. See `imagep()` for details.
`showBandName` | A logical indicating whether the band name is to plotted in the top margin, near the right-hand side.
`alpha.f` | Argument used if `col="natural"`, to adjust colors with `adjustcolor()`.
`red.f` | Argument used if `col="natural"`, to adjust colors with `adjustcolor()`. Higher values of `red.f` cause red hues to be emphasized (e.g. dry land).
`green.f` | Argument used if `col="natural"`, to adjust colors with `adjustcolor()`. Higher values of `green.f` emphasize green hues (e.g. forests).
`blue.f` | Argument used if `band="terralook"`, to adjust colors with `adjustcolor()`. Higher values of `blue.f` emphasize blue hues (e.g. ocean).
`offset` | Argument used if `band="terralook"`, to adjust colors with `adjustcolor()`.
`transform` | Argument used if `band="terralook"`, to adjust colors with `adjustcolor()`.
`debug` | Set to a positive value to get debugging information during processing.
`...` | optional arguments passed to plotting functions.

### Details

For Landsat-8 data, the `band` may be one of: `"aerosol"`, `"blue"`, `"green"`, `"red"`, `"nir"`, `"swir1"`, `"swir2"`, `"panchromatic"`, `"cirrus"`, `"tirs1"`, or `"tirs2"`.

For Landsat-7 data, `band` may be one of `"blue"`, `"green"`, `"red"`, `"nir"`, `"swir1"`, `"tirs1"`, `"tirs2"`, `"swir2"`, or `"panchromatic"`.

For Landsat data prior to Landsat-7, `band` may be one of `"blue"`, `"green"`, `"red"`, `"nir"`, `"swir1"`, `"tirs1"`, `"tirs2"`, or `"swir2"`.

If `band` is not given, the (`"tirs1"`) will be used if it exists in the object data, or otherwise the first band will be used.

In addition to the above there are also some pseudo-bands that can be plotted, as follows.

- Setting `band="temperature"` will plot an estimate of at-satellite brightness temperature, computed from the `tirs1` band.

- Setting `band="terralook"` will plot a sort of natural color by combining the `red`, `green`, `blue` and `nir` bands according to the formula provided at `https://lta.cr.usgs.gov/terralook/what_is_terralook` (a website that worked once, but failed as of Feb 2, 2017), namely that the `red`-band data are provided as the `red` argument of the `rgb()` function, while the `green` argument is computed as 2/3 of the `green`-band data plus 1/3 of the `nir`-band data, and the `blue` argument is computed as 2/3 of the `green`-band data minus 1/3 of the `nir`-band data. (This is not a typo: the `blue` band is not used.)

### Author(s)

Dan Kelley

### See Also

Other things related to landsat data: `[[,landsat-method`, `[[<-,landsat-method`, `landsat-class`, `landsatAdd()`, `landsatTrim()`, `landsat`, `read.landsat()`, `summary,landsat-method`

Other functions that plot oce data: `plot,adp-method`, `plot,adv-method`, `plot,amsr-method`, `plot,argo-method`, `plot,bremen-method`, `plot,cm-method`, `plot,coastline-method`, `plot,ctd-method`, `plot,gps-method`, `plot,ladp-method`, `plot,lisst-method`, `plot,lobo-method`, `plot,met-method`, `plot,odf-method`, `plot,rsk-method`, `plot,satellite-method`, `plot,sealevel-method`, `plot,section-method`, `plot,tidem-method`, `plot,topo-method`, `plot,windrose-method`, `plot,xbt-method`, `plotProfile()`, `plotScan()`, `plotTS()`, `tidem-class`

[Package oce version 1.2-0 Index]