library(MASS)
library(rbokeh)
fnmDictReal='/wk2/IdeaProjects/calf/doc/knn/dictR.csv'
nms <- c("sta","lat","lon","alt")
dictR <- read.csv(fnmDictReal,header = F)
names(dictR) <- nms
x <- dictR$lon
y <- dictR$lat
dx <- density(x)
dy <- density(y)
density<-kde2d(x,y,n=600)
filled.contour(density,color.palette=colorRampPalette(c('white','blue','yellow','red','darkred')))
评论