小编Orp*_*eus的帖子

我有一个名为"seoul032823"的81个观察的每小时PM10数据集.你可以从这里下载.我在这个数据集上进行了普通克里金法,并得到了克里金预测的空间图.我还可以在国家地图上显示观察数据点.但我不能在国家地图上重叠克里格空间预测图.

我想做什么:我想在南韩地图(不是整个韩国)上重叠我的空间预测图.我感兴趣的领域是北纬37.2N到37.7N和经度126.6E到127.2E.这意味着我需要从韩国地图中裁剪这个区域并在此处重叠预测图.我还需要根据浓度值显示原始观测数据点,这些点将遵循空间地图的颜色.例如,我想要这种类型的地图:

我的R代码用于克里金法,并在韩国地图上显示数据点:

library(sp)
library(gstat)
library(automap)
library(rgdal)
library(e1071)
library(dplyr)
library(lattice)

seoul032823 <- read.csv ("seoul032823.csv")

#plotting the pm10 data on Korea Map
library(ggplot2)
library(raster)

seoul032823 <- read.csv ("seoul032823.csv")
skorea<- getData("GADM", country= "KOR", level=1)
plot(skorea)

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=seoul032823, aes(x=LON, y=LAT), 
             colour= "red", alpha=0.7,na.rm=T) +
  #scale_size(range=c(2,4))+
  labs(title= "PM10 Concentration in Seoul Area at South Korea",
       x="Longitude", y= "Latitude", size="PM10(microgm/m3)")+
  theme(title= element_text(hjust = 0.5,vjust = 1,face= c("bold")))

# Reprojection
coordinates(seoul032823) <- ~LON+LAT …

我有一个数据帧包含名为ID,经度(LON),纬度(LAT)的变量,您可以从这里下载.我使用ggplot2包含以下代码的包在相同颜色的国家/地区地图上绘制了一些经度 - 纬度:

library(ggplot2)
library(raster)
read.csv("station.csv")
skorea<- getData("GADM", country= "KOR", level=1)
plot(skorea)
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=station, aes(x=LON, y=LAT), 
             colour= "red", alpha=1,na.rm=T) +
  scale_size(range=c(2,7))+
  labs(title= "coordinate data of seoul",
       x="Longitude", y= "Latitude")+
  theme(title= element_text(hjust = 0.5,vjust = 1,face= c("bold")))

我得到了以下情节

现在我希望其中一些点将根据其ID在不同颜色的地图上绘制.例如,我希望这些ID(111141,111142,111241,111281,111301,131141,131144,131161)将为蓝色,其余ID将保持红色.我怎样才能做到这一点？