使用Krigging创建热图

时间:2017-08-23 01:23:16

标签: r ggplot2 statistics interpolation kriging

我尝试使用Krigging创建一个好的热图来显示缺失值。 Data with RLevels

我有以下数据,其中包含已为RLevel测量的所有值。

我按照以下链接说明了如何使用krigging。 https://rpubs.com/nabilabd/118172

这是我写的以下代码。在这些步骤之前,我已经删除了我的DieData中的所有值 需要测试的值。需要测试的值在我的代码中被称为die.data.NAValues。

#**************************************************CODE***************** 
#Step3: Convert to SpatialPointsDataFrame Object 
coordinates(die.data) = ~X+Y 

#Step 4: Get the prediction Grid 
coordinates(die.data.NAValues)=~X+Y 


#Using autokride method 
kr = autoKrige(RLevel, die.data, die.data.NAValues,nmax=20) 
predicted_die_values <- kr$krige_output 
predicted_die_model <- kr$var_model 


#Get Predictions. Plot the predicted on heat map. 
g <- gstat(NULL,"RLevel",RLevel~1,die.data, model=predicted_die_model,nmax=1) 
predictedSet <- predict(g,newdata=die.data,BLUE=TRUE) 

#Plot the krigging graph 
predicted_die_values %>% as.data.frame %>% ggplot(aes(x=X,y=Y)) + geom_tile(aes(fill=v1.pred))+coord_equal() +scale_fill_gradient(low="yellow",high="red")+scale_x_continuous()+scale_y_continuous()+theme_bw() 

当我绘制图形时,我从KRIGING METHOD测试的值中得到以下图像。 enter image description here

我的问题是如何显示一个好的热图,其中包含来自KRIG的预测点和已有的点。我想让我的图表从上面发布的链接中显示出类似的内容。 enter image description here

关于我的数据集的说明:我的原始数据集(包括尚未测试的NA值)包含大约55057个点。当我拿出NA值并使用它是我的预测网格时,我得到390分。 RLevel的大部分值都在30到30的范围内,除了100-200点以上都在100以上。

任何人都可以帮助我或指导我如何制作好的热图吗?

0 个答案:

没有答案