#install.packages('data.table')
#install.packages('MASS')
#install.packages('rbokeh')
library(data.table)
library(MASS)
library(rbokeh)
# fnmDictReal='/wk2/prj/calf4dl/doc/knn/stas.csv'
fnmDictReal='/wk2/prj/calf4dl/doc/knn/dictA.csv'
dictR <- fread(fnmDictReal)
nms <- c("sta","lat","lon","alt")
names(dictR) <- nms
DT0 <- dictR[lat>18.0]
DT1 <- DT0[lat<55.0]
DT2 <- DT1[lon>60.0]
DT <- DT2[lon<140.0]
length(DT$sta)
length(dictR$sta)
dx <- density(DT$lon)
dy <- density(DT$lat)
dxy<- kde2d(DT$lon,DT$lat,n=600)
filled.contour(dxy,color.palette=colorRampPalette(c('white','blue','yellow','red','darkred')))
评论