如何在不使用Google地圖（圖像）的情況下在地圖上繪制數據？

Question

我有一個csv文件，其變量名為“緯度”，“經度”，“ PM10濃度”。 您可以在此處下載數據。 我想根據其緯度和經度在韓國地圖上繪制PM10數據。 我也想將它們顯示為具有不同大小和顏色的氣泡。

在此示例之后，我已經在Google Map上繪制了PM10數據。 但是現在我想不使用Google地圖而是通過創建空間對象或以其他任何方式來執行此操作。

我嘗試編寫一些代碼，但是我已經下載了韓國行政區域（GADM）的空間數據。 但是我不確定這種方法是對還是錯。

library(rgdal)
library(ggplot2)
library(maptools)

map<-readOGR('D:/BACKUP/R/GSTAT/R File/shape file korea map',layer ='KOR_adm2')
summary(kmap)

EPSG<-make_EPSG()
EPSG[grepl("WGS 84$", EPSG$note), ]

kmap84<-spTransform(kmap, CRS("+init=epsg:4326"))
kmaps<-fortify(kmap84)

我不知道下一步該怎么做。

Answer 1

這是一個例子：

library(raster)
library(ggplot2)
download.file("https://docs.google.com/uc?id=0ByY3OAw62EShakxJZkplOXZ0RGM&export=download", tf <- tempfile(fileext = ".csv"))
df <- read.csv(tf, row.names = 1)
skorea <- getData("GADM", country = "South Korea", level = 2)
skorea <- fortify(skorea)
ggplot() +  
  geom_map(data = skorea, map = skorea, aes(x = long, y = lat, map_id = id, group = group), 
           fill = NA, colour = "black") + 
  geom_point(data = df, aes(x = LON, y = LAT, size = PM10), colour = "red", alpha = .5) + 
  scale_size(range = c(1, 5))