r - Creating Heat Map using Krigging -
i'm trying create heat map using krigging missing values.
i have following data, contains values have been measured rlevel.
i followed following link tells how use krigging. https://rpubs.com/nabilabd/118172
this following code wrote. before these steps, had removed values diedata needed values tested. values need tested refered die.data.navalues in code.
#**************************************************code***************** #step3: convert spatialpointsdataframe object coordinates(die.data) = ~x+y #step 4: 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 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 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()
when plot graph, following image values have been tested kriging method.
my question how can show heat map predicted points krig , points have. want graph show link above had posted.
description dataset: original dataset including na values have not been tested contains around 55057 points. when take out na values , use prediction grid, 390 points. majority of values rlevel within 30's range except around 100-200 points above 100.
can me out or give me guidance of how produce heatmap?
wiki
Comments
Post a Comment