Координаты участка на карте

Я пытаюсь построить мои координаты с помощью R. Я уже пробовал следовать за другим сообщением (R: группировать сгруппированные координаты на карте мира; Построение координат нескольких точек на карте google в R), но я не очень хорошо разбираюсь в своих данных.

Я пытаюсь создать плоскую карту мира с координатой gps как цветные точки (каждая область имеет определенный цвет):

area         lat    long
Agullhas    -38,31  40,96
Polar       -57,59  76,51
Tasmanian   -39,47  108,93

library(RgoogleMaps)
lat <- c(-38.31, -35.50) #define our map ylim
lon <- c(40.96,37.50) #define our map xlim
center = c(mean(lat), mean(lon))  #tell what point to center on
zoom <- 2 #zoom: 1 = furthest out (entire globe), larger numbers = closer in
terrmap <- GetMap(center=center, zoom=zoom, maptype= "satallite", destfile = "satallite.png")

проблема в том, что теперь я не знаю, как добавить свои очки, и мне понравится один цвет для каждого региона.

Может ли кто-нибудь помочь мне продвигаться вперед?

другой вариант, который я пробовал:

library(maps)
library(mapdata)
library(maptools)
map(database= "world", ylim=c(-38.31, -35.5), xlim=c(40.96, 37.5), col="grey80", fill=TRUE, projection="gilbert", orientation= c(90,0,225))
lon <- c(-38.31, -35.5)  #fake longitude vector
lat <- c(40.96, 37.5)  #fake latitude vector
coord <- mapproject(lon, lat, proj="gilbert", orientation=c(90, 0, 225))  #convert points to projected lat/long
points(coord, pch=20, cex=1.2, col="red")  #plot converted points

но координаты заканчиваются в неправильном положении, и я не уверен, почему

Надеюсь, кто-то может помочь

Ответы

Ответ 1

В качестве альтернативы RgoogleMaps вы также можете использовать комбинацию ggplot2 с ggmap.

С помощью этого кода:

# loading the required packages
library(ggplot2)
library(ggmap)

# creating a sample data.frame with your lat/lon points
lon <- c(-38.31,-35.5)
lat <- c(40.96, 37.5)
df <- as.data.frame(cbind(lon,lat))

# getting the map
mapgilbert <- get_map(location = c(lon = mean(df$lon), lat = mean(df$lat)), zoom = 4,
                      maptype = "satellite", scale = 2)

# plotting the map with some points on it
ggmap(mapgilbert) +
  geom_point(data = df, aes(x = lon, y = lat, fill = "red", alpha = 0.8), size = 5, shape = 21) +
  guides(fill=FALSE, alpha=FALSE, size=FALSE)

вы получите этот результат:

Ответ 2

Другой альтернативой является пакет plotGoogleMaps, который позволяет отображать в навигаторе, позволяя увеличивать и уменьшать масштаб и т.д. Затем вы можете сделать снимок экрана своего снимка, чтобы сохранить его (хотя помните, что карты Google на законных основаниях для использования в Интернете).

 library("plotGoogleMaps")
 lat <- c(-38.31, -35.50) #define our map ylim
 lon <- c(40.96,37.50) #define our map xlim

 # make your coordinates a data frame 
 coords <- as.data.frame(cbind(lon=lon,lat=lat))

 # make it a spatial object by defining its coordinates in a reference system
 coordinates(coords) <- ~lat+lon 

 # you also need a reference system, the following should be a fine default
 proj4string(coords) <- CRS("+init=epsg:4326")
 # Note: it is a short for: 
 CRS("+init=epsg:4326")
 > CRS arguments:
 >  +init=epsg:4326 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 

 # then just plot
 a <- plotGoogleMaps(coords)
 # here `a <-` avoids that you get flooded by the html version of what you plot

И вы получите:

Ответ 3

Вот решение, использующее только Rgooglemaps, по просьбе пользователя.

# get map (from askers OP, except changed map type = "Satallite" to type = "Satellite")
library(RgoogleMaps)
lat <- c(-38.31, -35.50) #define our map ylim
lon <- c(40.96,37.50) #define our map xlim
center = c(mean(lat), mean(lon))  #tell what point to center on
zoom <- 2 #zoom: 1 = furthest out (entire globe), larger numbers = closer in
terrmap <- GetMap(center=center, zoom=zoom, type= "satellite", destfile = "satellite.png")

# plot points and save image
lat <- c(-38.31, -57.59, -39.47)
lon <- c(40.96, 76.51, 108.93)
png('map.png')
PlotOnStaticMap(terrmap, lat = lat, lon = lon, pch = 20, col = c('red', 'blue', 'green'))
dev.off()