r - 使用 R 在 GADM map 上绘制城市

Question

我正在尝试使用 GADM 数据并使用 R 绘制印度州/地区的详细 map 。我使用了以下代码

# Load required libraries
library(sp)
library(RColorBrewer)
# ---------------------------------------------------------------------------
# load level 2 india data downloaded from http://gadm.org/country
load("IND_adm2.RData")
ind2 = gadm
# plotting districts of a State, in this case West Bengal
wb2 = (ind2[ind2$NAME_1=="West Bengal",])
spplot(wb2,"NAME_1", main = "West Bengal Districts",
       colorkey=F, scales=list(draw=T))

生成此 map

<小时/>

我现在尝试放入几个标记(或点加文本)来显示地区总部“Purulia”，位于 lon = 86.36521 lat = 23.33208 但不知何故我无法获取 sp.layout 语法正确。我需要这个来处理一组长、纬度都已知的城镇。如果有人能帮助我解决这个问题，我将非常感激。

Answer 1

这是一个简短的示例:

library("sp")

# load level 2 india data from gadm.org
library("raster")
ind2 <- getData('GADM', country='IND', level=2)
wb2 <- ind2[ind2$NAME_1=="West Bengal",]

cities <- data.frame(name="Purulia", lon=86.36521, lat=23.33208)

spplot(wb2, "NAME_1",
       sp.layout=list("panel.points", cities$lon, cities$lat, col="red"),
       main="West Bengal Districts",
       colorkey=FALSE, scales=list(draw=TRUE))

或者如果您的城市位于 SpatialPointsDataFrame 中:

cities <- data.frame(name="Purulia", lon=86.36521, lat=23.33208)
coordinates(cities) <- ~ lon + lat
class(cities)
# [1] "SpatialPointsDataFrame"

spplot(wb2, "NAME_1",
       sp.layout=list("sp.points", cities, col="red"),
       main="West Bengal Districts",
       colorkey=FALSE, scales=list(draw=TRUE))