# set.seed(2)
# dd <- data.frame(lon = runif(100, min = 110, max = 117),
# lat = runif(100, min = 31, max = 37),
# value = runif(100, min = -10, max = 10))
# grd <- points2grid(dd, extent = c(110, 117, 31, 37), gridsize = 0.5)
# grd <- na.omit(grd)
# awa(grd$value, grd$lat) # area weighted average
Run the code above in your browser using DataLab