Description Usage Arguments Details Value Note Author(s) See Also Examples
Get MODIS tile ID(s) for a specific geographic area.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | ## S4 method for signature 'missing,missing,missing'
getTile(mode = c("click", "draw"), ...)
## S4 method for signature 'missing,charORnum,charORnum'
getTile(x, tileH, tileV, ...)
## S4 method for signature 'character,ANY,ANY'
getTile(x, tileH, tileV, ...)
## S4 method for signature 'Raster,ANY,ANY'
getTile(x, tileH, tileV, ...)
## S4 method for signature 'map,ANY,ANY'
getTile(x, tileH, tileV, ...)
## S4 method for signature 'Extent,ANY,ANY'
getTile(x, tileH, tileV, ...)
## S4 method for signature 'bbox,ANY,ANY'
getTile(x, tileH, tileV, ...)
## S4 method for signature 'Spatial,ANY,ANY'
getTile(x, tileH, tileV, ...)
## S4 method for signature 'sf,ANY,ANY'
getTile(x, tileH, tileV, ...)
|
mode |
Interactive selection mode as |
... |
Additional arguments passed to |
x |
Extent information, see Details. |
tileH, tileV |
|
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: | |
A map object.
|
A MODISextent
object.
MODIS does no longer support the tile identification and automated
download of MERIS and SRTM data. At least as far as the latter is concerned,
easy data access is granted through getData
.
Matteo Mattiuzzi, Florian Detsch
extent
, st_bbox
, map
,
search4map
.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 | ## Not run:
# ex 1 ############
# interactive tile selection
getTile()
getTile(mode = "draw")
## End(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)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.