Saya ingin mengukur biomassa di lumpur pantai. Saya hanya dapat mengakses poin di dalam poligon. Adakah metode yang tersedia yang memungkinkan saya memperkirakan nilai titik di luar poligon, berdasarkan nilai titik di dalam poligon?
set.seed(5)
x <- rnorm(50, -1.841, 0.01)
set.seed(50)
y <- rnorm(50, 55.663, 0.01)
xy <- data.frame(x,y, values=rnorm(50))
coordinates(xy) <- c("x", "y")
proj4string(xy) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")
plot(xy)
makePolygons <- function(coordsx, coordsy){
coords <- matrix(c(c(coordsx, coordsy)), ncol=2)
p <- Polygon(coords)
p <- Polygons(list(p), ID = "p")
myPoly <- SpatialPolygons(list(p))
spdf = SpatialPolygonsDataFrame(myPoly, data.frame(variable1 = c(2),
variable2 = c(3), row.names = c("p")))
proj4string(spdf) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")
print("polygon is in longlat!!!")
spdf
}
myPoly <- makePolygons(coordsx=c(-1.841960, -1.843464, -1.888623, -1.841960),
coordsy=c(55.633696, 55.68178, 55.63841, 55.633696))
plot(myPoly, add=T)
Jawaban:
Tebakan terbaik saya adalah mengubah poligon menjadi titik grid dan memperkirakan nilai di setiap titik yang tidak tumpang tindih dengan poin yang berisi data. Ada tutorial yang cukup rapi di sini .
sumber