Tôi muốn đo sinh khối trên bãi bồi ven biển. Tôi chỉ có thể truy cập các điểm bên trong đa giác. Có phương pháp nào khả dụng cho phép tôi ước tính các giá trị của các điểm bên ngoài đa giác, dựa trên các giá trị của các điểm bên trong đa giác không?
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)
Thảm thực vật nằm trên bãi bồi trong khu vực ngập triều - một số khu vực quá nguy hiểm để đến. Sinh khối thực vật là một biến dự đoán
—
luciano
Bạn có một vấn đề quan trọng, vì rất có khả năng không thể tiếp cận và sinh khối có liên quan. Điều đó làm cho nó không hợp lệ để ngoại suy dữ liệu thu được từ những nơi có thể truy cập đến tất cả các nơi. Đối với một cách tiếp cận hợp lệ, bạn phải tìm một cách nào đó - ngay cả khi chỉ là một cách thay thế - để đo một phần đại diện của các khu vực không thể tiếp cận. Kriging (và hầu hết các quy trình đường viền khác) sẽ giải quyết vấn đề một cách tuyệt vời và phần mềm sẽ vui vẻ cung cấp cho bạn kết quả rất chi tiết, cực kỳ sai.
—
whuber
Cách tiếp cận của tôi sẽ liên quan đến ước tính sinh khối trên mặt đất với các giá trị NDVI, có lẽ dựa trên dữ liệu Landsat. Sử dụng hồi quy để dự đoán sinh khối từ NDVI trong khu vực nguy hiểm.
—
Aaron
@whuber mặc dù các khu vực trong đa giác màu đỏ có thể truy cập được, chúng không bao giờ được sử dụng bởi mọi người.
—
luciano
Hiện có vấn đề gì? Làm thế nào mà thay đổi bản chất của nghiên cứu của bạn hoặc quy trình lấy mẫu?
—
whuber