MODIS (version 1.2.11)

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, ...)
# S4 method for sfc,ANY,ANY
getTile(x, tileH, tileV, ...)

Value

A MODISextent object.

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.

Author

Matteo Mattiuzzi, Florian Detsch

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 maps::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().
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, sfc, Spatial:
Except for resolution, same as for Raster*.
Other:

See Also

raster::extent(), sf::st_bbox(), maps::map(), search4map().

Examples

if (FALSE) {
# ex 1 ############
# interactive tile selection
getTile()
getTile(mode = "draw")
}

# ex 2: Spatial ############
dsn <- system.file("ex/lux.shp", package = "terra")
Up <- raster::shapefile(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)