Как я могу добавить точку из базы данных в Tmap?

Question

Я хотел бы добавить точки наблюдения на карту. В частности, мои данные по очкам имеют следующую форму:

cita$lat
  [1] 41.88333 41.88333 44.95971 45.16110 45.16110 45.16110 45.16110
  [8] 45.16110 45.16110 45.16110 45.16110 45.16110 45.16110 45.16110
 [15] 45.16110 45.16110 45.16110 45.29775 45.29775 45.13333 45.16110
 [22] 45.29775 45.29775 45.29775 45.29775 45.29775 45.18885 45.19092

cita$lng
  [1] 12.500000 12.500000 10.689220  9.701432  9.701432  9.701432
  [7]  9.701432  9.701432  9.701432  9.701432  9.701432  9.701432
 [13]  9.701432  9.701432  9.701432  9.701432  9.701432 11.658382
 [19] 11.658382 10.033333  9.701432 11.658382 11.658382 11.658382
 [25] 11.658382 11.658382  9.690454  9.726434 

Я создал карту следующим образом:

library(rnaturalearth)
ita = ne_countries(country = "Italy") # United States borders
class(ita)
ita_sf = st_as_sf(ita)

tm_shape(ita_sf) +
  tm_fill() +
  tm_borders() 

Как я могу это сделать?

спасибо вы!

Answer 1

library(rnaturalearth)
library(sf)
library(tmap)

ita = ne_countries(country = "Italy") # United States borders
class(ita)
ita_sf = st_as_sf(ita)

cita_sf = st_as_sf(cita, coords = c('long', 'lat'), crs = st_crs(ita_sf)$proj4string)

tm_shape(ita_sf) +
  tm_fill() +
  tm_borders() +
tm_shape(cita_sf) +
  tm_dots()

Карта

Данные

cita <- data.frame(
  lat = c(41.88333, 41.88333, 44.95971, 45.16110, 45.16110, 45.16110, 45.16110,
          45.16110, 45.16110, 45.16110, 45.16110, 45.16110, 45.16110, 45.16110,
          45.16110, 45.16110, 45.16110, 45.29775, 45.29775, 45.13333, 45.16110,
          45.29775, 45.29775, 45.29775, 45.29775, 45.29775, 45.18885, 45.19092),
  long = c(12.500000, 12.500000, 10.689220,  9.701432,  9.701432,  9.701432,
            9.701432,  9.701432,  9.701432,  9.701432,  9.701432,  9.701432,
            9.701432,  9.701432,  9.701432,  9.701432,  9.701432, 11.658382,
            11.658382, 10.033333,  9.701432, 11.658382, 11.658382, 11.658382,
            11.658382, 11.658382,  9.690454,  9.726434)
)