2つのデータフレームがあります(ddd, env0)。これらを緯度経度の近さに基づいて空間結合させたいと思っています。ただし、近い緯度経度は同じstormから選ばれるようにしたいです。
head(ddd)
storm lat long Mesh3rd
1 270 34.88750 137.2562 52372260
2 270 34.88750 137.2688 52372261
3 270 34.89583 137.2562 52372270
4 270 34.89583 137.2688 52372271
5 270 34.89583 137.2812 52372272
6 270 34.89583 137.2937 52372273
head(env0)
storm time track_time lat long track_env track_pres track_wind
1 40 1 55286.00 24.99216 139.5000 1010.416 1008.19 13.22159
2 40 2 55286.25 24.99216 140.6250 1008.091 1005.41 14.09173
3 40 3 55286.50 24.43054 142.3125 1008.779 1006.15 15.27882
4 40 4 55286.75 24.99216 144.0000 1006.264 1003.04 14.82848
5 40 5 55287.00 24.99216 144.5625 1007.882 1004.36 14.56058
6 40 6 55287.25 24.43054 145.6875 1005.603 1002.48 15.59767
試したコードは以下の通りです。
temp <- ddd
temp3 <- split(ddd, f = dd$storm)
coordinates(temp) <- ~long+lat
mylist <- split(temp, f = temp$storm)
names(env0)<-c("storm","time","track_time","lat","long","track_env",
"track_pres","track_wind")
u<-unique(temp$storm)
envv<-filter(env0, storm %in% u)
temp1 <- envv
temp4 <- split(envv, f = envv$storm)
coordinates(temp1) <- ~long+lat
mylist1 <- split(temp1, f = temp1$storm)
test<-mylist
for(i in 1:length(mylist)){
test[[i]]<-pointDistance(mylist[[i]],mylist1[[i]],lonlat = TRUE, allpairs = TRUE)
}
nearest1<-temp3[[1]]
for(i in 1:length(test)){
nearest1[[i]] <- apply(test[[i]], 1, which.min)
}
sptsd<-mylist[[1]]
for(i in 1:length(mylist)){
sptsd[[i]]<- cbind(temp3[[i]], temp4[[i]][nearest1[[i]],1:8])
}