1

我正在尝试使用克里金法为缺失值创建一个好的热图。 具有 RLevel 的数据

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

我按照以下链接说明如何使用克里金法。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 方法测试的值中得到以下图像。 在此处输入图像描述

我的问题是如何用 KRIG 的预测点和已有的点显示一个好的热图。我希望我的图表从我发布的上面的链接中显示类似的内容。在此处输入图像描述

关于我的数据集的描述：我的原始数据集（包括未经测试的 NA 值）包含大约 55057 个点。当我取出 NA 值并使用它作为我的预测网格时，我得到 390 分。RLevel 的大多数值都在 30 的范围内，除了大约 100-200 点高于 100。

谁能帮助我或指导我如何制作好的热图？

4

0 回答 0