Có thể tạo một ô đường viền 3D không có dữ liệu liên tục trong R không?

Question

Tôi muốn tạo một đường bao của biến z với dữ liệu x, y, z. Tuy nhiên, có vẻ như chúng tôi cần cung cấp dữ liệu theo thứ tự tăng dần.

Tôi đã cố gắng sử dụng một số mã nhưng nó đã cho tôi lỗi.

Tôi đã thử các mã sau: Trial 1:

age2100 <- read.table("temp.csv",header=TRUE,sep=",") 

x <- age2100$x 
y <- age2100$y 
z <- age2100$z 

contour(x,y,z,add=TRUE,col="black") 

Tôi đã nhận lỗi sau

Error in contour.default(x, y, z, add = TRUE, col = "black") : increasing 'x' and 'y' values expected 

sau đó tôi cố gắng sử dụng ggplot2 để tạo ra các đường viền. Tôi sử dụng đoạn mã sau:

library("ggplot2") 
library("MASS") 
library("rgdal") 
library("gpclib") 
library("maptools") 
age2100 <- read.table("temp.csv",header=TRUE,sep=",") 
v <- ggplot(age2100, aes(age2100$x, age2100$y,z=age2100$z))+geom_contour() 
v 

Tôi đã nhận lỗi sau đây: Thông điệp

Cảnh báo:

Not possible to generate contour data 

hãy tìm các dữ liệu về vị trí sau https://www.dropbox.com/s/mg2bo4rcr6n3dks/temp.csv

Ai có thể cho tôi biết làm thế nào để tạo ra các dữ liệu đường viền từ biến thứ ba (z) từ temp.csv? Tôi cần phải làm như vậy nhiều lần vì vậy tôi đang cố gắng làm trên R thay vì Arcgis.

Answer 1

Dưới đây là một ví dụ về cách một nội suy sử dụng interp từ akima gói:

age2100 <- read.table("temp.csv",header=TRUE,sep=",") 

x <- age2100$x 
y <- age2100$y 
z <- age2100$z 

require(akima) 

fld <- interp(x,y,z) 

par(mar=c(5,5,1,1)) 
filled.contour(fld) 

Dưới đây là một âm mưu thay thế bằng cách sử dụng image chức năng (điều này cho phép một số linh hoạt để thêm âm mưu cấp dưới chức năng (yêu cầu hàm image.scale, được tìm thấy here):

source("image.scale.R") # http://menugget.blogspot.de/2011/08/adding-scale-to-image-plot.html 

x11(width=5, height=6) 
layout(matrix(c(1,2), nrow=1, ncol=2), widths=c(4,1), height=6, respect=TRUE) 
layout.show(2) 

par(mar=c(4,4,1,1)) 
image(fld) 
contour(fld, add=TRUE) 
points(age2100$x,age2100$y, pch=".", cex=2) 

par(mar=c(4,0,1,4)) 
image.scale(fld$z, xlab="", ylab="", xaxt="n", yaxt="n", horiz=FALSE) 
box() 
axis(4) 
mtext("text", side=4, line=2.5)