使用形状文件的Choropleth贴图未在输出中显示结果

问题描述

Refering this for coding

我正在执行克里金法，并希望将结果显示在形状文件（整体图）中。预测结果已生成，但输出未显示结果。

library(raster)
boroughs <- shapefile("C:\\Downloads\\IND_adm (1)\\IND_adm3.shp")
proj4string(boroughs)
class(boroughs)
plot(boroughs)

boroughoutline <- fortify(boroughs)

krig<-krige(formula=Reach ~ 1,locations=rd,newdata=Grd,model=model.variog)
krig.output=as.data.frame(krig)
names(krig.output)[1:3]<-c("long","lat","var1.pred")

plot<-ggplot(data=krig.output,aes(x=long,y=lat))
layer1<-c(geom_tile(data=krig.output,aes(fill=var1.pred)))
layer2<-c(geom_path(data=boroughoutline,aes(long,lat,group=group),colour = "grey40",size=1))
plot+layer1+layer2+scale_fill_gradient(low="#FEEBE2",high="#7A0177")+coord_equal()

在输出中唯一显示的没有颜色的形状文件

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）