MODIS (version 1.1.7)

getTile: Get MODIS Tile ID(s)

Description

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

Usage

getTile(x, tileH = NULL, tileV = NULL, mode = c("click", "draw"), ...)

Arguments

x

Extent information, see Details. Ignored if tileH and tileV are specified.

tileH, tileV

numeric or character. Horizontal and vertical tile number(s) of the MODIS Sinusoidal grid (e.g., tileH = 1:5). If specified, no cropping is performed and the full tile(s) (if more than one then also mosaicked) is (are) processed!

mode

Interactive selection mode as character. Available options are "click" (default) and "draw" that trigger interactive MODIS tile selection and free feature drawing, respectively. Ignored if any of 'x' or 'tileH,tileV' is NOT missing.

...

Additional arguments passed to MODISoptions. Here, only 'outProj' and 'pixelSize' are relevant, and this only if 'x' is an object of class character, map, Extent or bbox.

Value

A MODISextent object.

Details

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

missing:
If 'tileH,tileV' are specified, 'x' will be ignored. If no such tile indices are provided and 'x' is missing, a viewer window pops up that allows interactive tile selection from the global MODIS Sinusoidal grid or, if mode = "draw", free feature drawing.
character:
Either the country name of a map object (see map) or a valid file path of a raster image or ESRI shapefile (shp). The former approach also supports pattern matching via regular expressions.
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 and bbox objects are assumed to be in EPSG:4326 as such objects have no projection information attached.
Other: missing:

See Also

extent, st_bbox, map, search4map.

Examples

Run this code
# NOT RUN {
# ex 1 ############
# interactive tile selection
getTile()
# }
# 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)

# }

Run the code above in your browser using DataCamp Workspace