oce (version 0.9-18)

landsatTrim: Trim a landsat Image to a Geographical Region

Description

Trim a landsat image to a latitude-longitude box. This is only an approximate operation, because landsat images are provided in x-y coordinates, not longitude-latitude coordinates.

Usage

landsatTrim(x, ll, ur, box, debug = getOption("oceDebug"))

Arguments

x

A landsat object, e.g. as read by read.landsat.

ll

A list containing longitude and latitude, for the lower-left corner of the portion of the image to retain.

ur

A list containing longitude and latitude, for the upper-right corner of the portion of the image to retain.

box

A list containing x and y (each of length 2), corresponding to the values for ll and ur, such as would be produced by a call to locator(2).

debug

A flag that turns on debugging. Set to 1 to get a moderate amount of debugging information, or a higher value for more debugging.

return An object of class "landsat", with data having been

concept

satellite

Details

As of June 25, 2015, the matrices storing the image data are trimmed to indices determined by linear interpolation based on the location of the ll and ur corners within the lon-lat corners specified in the image data. (A previous version trimmed in UTM space, and in fact this may be done in future also, if a problem in lonlat/utm conversion is resolved.) An error results if there is no intersection between the trimming box and the image box.

See Also

The documentation for landsat-class explains the structure of landsat objects, and also outlines the other functions dealing with them.

Other functions dealing with satellite data: [[,landsat-method, amsr-class, g1sst-class, landsatAdd, plot,amsr-method, plot,landsat-method, plot,satellite-method, read.amsr, read.g1sst, read.landsat, satellite-class, summary,amsr-method, summary,landsat-method, summary,satellite-method