cf. - https://geocompr.robinlovelace.net - https://mhweber.github.io/AWRA_2020_R_Spatial
library(leaflet)
popup = c("Calvin University", "University of Colorado", "ITRI, University of Brighton", "CSIRO, Sydney", "Pella, IA")
leaflet() %>%
addProviderTiles("NASAGIBS.ViirsEarthAtNight2012") %>%
addMarkers(lng = c(-85.670006, -105.270546, -0.152778, 151.209900, -92.916423),
lat = c(42.963795, 40.014984, 50.827778, -33.865143, 41.408059),
popup = popup)library(sf) # classes and functions for vector data## Linking to GEOS 3.9.1, GDAL 3.2.1, PROJ 7.2.1; sf_use_s2() is TRUElibrary(terra) # classes and functions for raster data## terra version 1.4.20library(spData) # load geographic data
library(spDataLarge) # load larger geographic dataworld## Simple feature collection with 177 features and 10 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: -180 ymin: -89.9 xmax: 180 ymax: 83.64513
## Geodetic CRS: WGS 84
## # A tibble: 177 x 11
## iso_a2 name_long continent region_un subregion type area_km2 pop lifeExp
## * <chr> <chr> <chr> <chr> <chr> <chr> <dbl> <dbl> <dbl>
## 1 FJ Fiji Oceania Oceania Melanesia Sove~ 1.93e4 8.86e5 70.0
## 2 TZ Tanzania Africa Africa Eastern ~ Sove~ 9.33e5 5.22e7 64.2
## 3 EH Western ~ Africa Africa Northern~ Inde~ 9.63e4 NA NA
## 4 CA Canada North Am~ Americas Northern~ Sove~ 1.00e7 3.55e7 82.0
## 5 US United S~ North Am~ Americas Northern~ Coun~ 9.51e6 3.19e8 78.8
## 6 KZ Kazakhst~ Asia Asia Central ~ Sove~ 2.73e6 1.73e7 71.6
## 7 UZ Uzbekist~ Asia Asia Central ~ Sove~ 4.61e5 3.08e7 71.0
## 8 PG Papua Ne~ Oceania Oceania Melanesia Sove~ 4.65e5 7.76e6 65.2
## 9 ID Indonesia Asia Asia South-Ea~ Sove~ 1.82e6 2.55e8 68.9
## 10 AR Argentina South Am~ Americas South Am~ Sove~ 2.78e6 4.30e7 76.3
## # ... with 167 more rows, and 2 more variables: gdpPercap <dbl>,
## # geom <MULTIPOLYGON [°]>plot(world["geom"])
st_crs(world)## Coordinate Reference System:
## User input: EPSG:4326
## wkt:
## GEOGCRS["WGS 84",
## DATUM["World Geodetic System 1984",
## ELLIPSOID["WGS 84",6378137,298.257223563,
## LENGTHUNIT["metre",1]]],
## PRIMEM["Greenwich",0,
## ANGLEUNIT["degree",0.0174532925199433]],
## CS[ellipsoidal,2],
## AXIS["geodetic latitude (Lat)",north,
## ORDER[1],
## ANGLEUNIT["degree",0.0174532925199433]],
## AXIS["geodetic longitude (Lon)",east,
## ORDER[2],
## ANGLEUNIT["degree",0.0174532925199433]],
## USAGE[
## SCOPE["Horizontal component of 3D system."],
## AREA["World."],
## BBOX[-90,-180,90,180]],
## ID["EPSG",4326]]crs(world)## Coordinate Reference System:
## Deprecated Proj.4 representation: +proj=longlat +datum=WGS84 +no_defs
## WKT2 2019 representation:
## GEOGCRS["WGS 84",
## DATUM["World Geodetic System 1984",
## ELLIPSOID["WGS 84",6378137,298.257223563,
## LENGTHUNIT["metre",1]]],
## PRIMEM["Greenwich",0,
## ANGLEUNIT["degree",0.0174532925199433]],
## CS[ellipsoidal,2],
## AXIS["geodetic latitude (Lat)",north,
## ORDER[1],
## ANGLEUNIT["degree",0.0174532925199433]],
## AXIS["geodetic longitude (Lon)",east,
## ORDER[2],
## ANGLEUNIT["degree",0.0174532925199433]],
## USAGE[
## SCOPE["Horizontal component of 3D system."],
## AREA["World."],
## BBOX[-90,-180,90,180]],
## ID["EPSG",4326]]library(tmap)
library(rcartocolor)
library(spDataLarge)
library(terra)
# data read ---------------------------------------------------------------
cla_raster = rast(system.file("raster/srtm.tif", package = "spDataLarge"))
# plots create ------------------------------------------------------------
rast_srtm = tm_shape(cla_raster) +
tm_raster(palette = carto_pal(7, "Geyser"),
title = "Elevation (m)", style = "cont")
tmap_arrange(rast_srtm)