MODIS (version 1.2.2)

getTile: Get MODIS Tile ID(s)

Description

Get MODIS tile ID(s) for a specific geographic area.

Usage

# S4 method for missing,missing,missing
getTile(mode = c("click", "draw"), ...)
# S4 method for missing,charORnum,charORnum
getTile(x, tileH, tileV, ...)
# S4 method for character,ANY,ANY
getTile(x, tileH, tileV, ...)
# S4 method for Raster,ANY,ANY
getTile(x, tileH, tileV, ...)
# S4 method for map,ANY,ANY
getTile(x, tileH, tileV, ...)
# S4 method for Extent,ANY,ANY
getTile(x, tileH, tileV, ...)
# S4 method for bbox,ANY,ANY
getTile(x, tileH, tileV, ...)
# S4 method for Spatial,ANY,ANY
getTile(x, tileH, tileV, ...)
# S4 method for sf,ANY,ANY
getTile(x, tileH, tileV, ...)

Arguments

mode

Interactive selection mode as character. Available options are "click" (default) and "draw" that trigger interactive MODIS tile selection and free feature drawing, respectively. Triggered only if 'x' and tile IDs are omitted.

...

Additional arguments passed to MODISoptions, see Details.

x

Extent information, see Details.

tileH, tileV

numeric or character. Horizontal and vertical tile number(s) of the MODIS Sinusoidal grid (e.g., tileH = 1:5). Cropping is disabled here and full tiles (if more than one then also mosaicked) are processed instead. Ignored if x is specified.

Value

A MODISextent object.

Details

Unless stated otherwise in the following, target 'outProj' and 'pixelSize' are carried over from MODISoptions.

If 'x' is of class (see Examples for use cases)

missing:
If tile IDs (see Arguments) are also missing, a viewer window pops up that allows for interactive tile selection from the global MODIS Sinusoidal grid or, if mode = "draw", free feature drawing.
character:
The country name of a map object (see map) with pattern matching via regular expressions enabled. Alternatively, a valid file path to a single ESRI shapefile (.shp) or an image readable by raster.
Raster*:
Spatial extent, resolution, and projection of the specified Raster* are determined automatically. This information is used by runGdal to create perfectly matching files. If the Raster* comes with no valid CRS, EPSG:4326 is assumed.
Extent, bbox:
Boundary coordinates from Extent objects are generally assumed to be in EPSG:4326 as such objects have no projection information attached. The same applies for 'bbox' objects lacking CRS information.
sf, Spatial:
Except for resolution, same as for Raster*.
Other:	missing:

See Also

extent, st_bbox, map, search4map.

Examples

# NOT RUN {
# ex 1 ############
# interactive tile selection
getTile()
getTile(mode = "draw")
# }
# NOT RUN {
# ex 2: Spatial (taken from ?rgdal::readOGR) ############
dsn <- system.file("vectors/Up.tab", package = "rgdal")[1]
Up <- rgdal::readOGR(dsn, "Up")
getTile(Up)

# ex 3: sf ############
ifl <- system.file("shape/nc.shp", package = "sf")
nc <- sf::st_read(ifl, quiet = TRUE)
getTile(nc)

# ex 4: tileH,tileV ############
getTile(tileH = 18:19, tileV = 4)

# ex 5: Raster* with valid CRS ############
rst1 <- raster(xmn = 9.2, xmx = 17.47, ymn = 46.12, ymx = 49.3)
getTile(rst1)

# this also works for projected data
rst3 <- projectExtent(rst1, crs = "+init=epsg:32633")
getTile(rst3)

# ex 6: Raster* without CRS or, alternatively, Extent or bbox --> treated as EPSG:4326 ############
mat2 <- matrix(seq(180 * 360), byrow = TRUE, ncol = 360)
rst2 <- raster(mat2, xmn = -180, xmx = 180, ymn = -90, ymx = 90)
getTile(rst2)
getTile(extent(rst1))
getTile(sf::st_bbox(nc))

# ex 7: map names as returned by search4map() ############
getTile("Austria")
getTile(c("Austria", "Germany"))

# or search for specific map name patterns (use with caution):
m1 <- search4map("Per")
getTile(m1)

# or use 'map' objects directly (remember to use map(..., fill = TRUE)): 
m2 <- map("state", region = "new jersey", fill = TRUE)
getTile(m2)

# }