Контур - Участок - по возрастанию

Question

У меня есть матрица, содержащая положение (X, Y, высота). Я добавляю столбец в матрицу, которую я называю «индексом». Я создаю вектор X и Y из матрицы. Оба они включают индексный столбец. Затем я сортирую по возрастанию значение вектора X и Y, который я только что сделал. Затем я создаю матрицу Z, содержащую высоту, и я отношу ее к позиции с использованием индекса. Затем я пытаюсь использовать контур команды (я хочу построить контурный график), и я получаю ошибку, говоря, что X и Y должны быть в порядке возрастания ... которые я только что сделал !!! Что я сделал не так?

noeud<-read.table("position.out") 
Matrice_Noeud<-matrix(ncol = ncol(noeud), nrow=nrow(noeud)) 
for (i in 1:nrow(noeud)) { 
    for (j in 1:ncol(noeud)) { 
    Matrice_Noeud[i,j]<-noeud[i,j] 
    } 
} 

Matrice_Noeud <- cbind(Matrice_Noeud, c(seq(1,nrow(noeud),1))) 


x<-data.frame(x=Matrice_Noeud[,1],Index=Matrice_Noeud[,4]) 
y<-data.frame(y=Matrice_Noeud[,2],Index=Matrice_Noeud[,4]) 

X<-x[order(x$x),] 
Y<-y[order(y$y),] 


Z<-matrix(NA, ncol=nrow(noeud),nrow=nrow(noeud)) 
for (x_i in 1:nrow(noeud)) { 
    for (y_i in 1:nrow(noeud)) { 
    if (Y$Index[y_i]==X$Index[x_i]) { 
     niveau<-which(Matrice_Noeud[,4]==Y$Index[y_i]) 
     Z[x_i,y_i]<-Matrice_Noeud[niveau,3] 
    } 


    } 
    } 

Xx<-array(X[,1]) 
Yy<-array(Y[,1]) 
Zz<-data.frame(Z) 

contour(Xx,Yy,Zz) 

Answer 1

ОК, так как я начал это делать, я это сделал.

#### making example data 
## assumptions: length(unique(x))=19, length(unique(y))=12, nrow(data)=121 
## (They mean the number of grid points is 19 * 12 = 228, but z.value is only 121.) 
xyz.f <- function(m, n) - m + (n - 7)^2 + 16   # make z from x and y (it means nothing special) 
xyz <- cbind(xyz <- expand.grid(x = round(seq(11,15,,19), 2), y = round(seq(6,10,,12), 2)), 
       z = apply(xyz, 1, function(k) xyz.f(k[1], k[2]))) 
set.seed(1); ind <- sample(19*12, 121)    # decide to use the 121 z of 19*12 
noeud <- as.matrix(xyz[ind,])      # example data maked out 

#### making contour()'s arguments 
Xx <- sort(unique(noeud[,1])) 
Yy <- sort(unique(noeud[,2])) # nrow(noeud); length(Xx); length(Yy) # OK (121, 19, 12) 

Zz <- matrix(NA, ncol=length(Yy), nrow=length(Xx))   # make 19 x 12 Z matrix (empty) 
# In each row, calculate x (y) value is what number in Xx (Yy) (= the position in Z matrix) 
X0 <- as.numeric(factor(noeud[,1])) # (edit) using Mr.Tufte's code in R help mailing. 
Y0 <- as.numeric(factor(noeud[,2])) 
apply(cbind(X0, Y0, noeud[,3]), 1, function (a) Zz[ a[1], a[2] ] <<- a[3]) 
## contour()'s arguments (Xx, Yy, Zz) maked out  
contour(Xx, Yy, Zz, xlab="including NAs") # length(Zz); length(Zz[!is.na(Zz)]) # OK (228,121) 

#### interpolating 
## I know few packages having interpolation functions. 
library(akima)      # use cubic spline interpolation methods of H. Akima 
NOEUD <- interp(noeud[,1], noeud[,2], noeud[,3]) 

#### results 
par.old <- par(no.readonly=T); par(mfrow=c(1,3), mar=c(4,0,1,0)) 
contour(Xx, Yy, Zz, xlab="including NAs", yaxt="n") # the including NAs data 
contour(NOEUD, xlab="Akima interpolation", yaxt="n") # the Akima interpolation data 
contour(Xx, Yy, matrix(xyz[,3], nrow=19), xlab="origin", yaxt="n") # the origin data 
# (edit) I noticed some interp()'s arguments make a difference (default: linear=T, extrap=F). 
contour(interp(noeud[,1], noeud[,2], noeud[,3], linear=T, extrap=F), xlab="Akima interp() default") 
contour(interp(noeud[,1], noeud[,2], noeud[,3], linear=F, extrap=F), xlab="interp(linear=F)") 
contour(interp(noeud[,1], noeud[,2], noeud[,3], linear=F, extrap=T), xlab="interp(linear=F, extrap=T)") 
par(par.old) 

### supplement (using the same data, output is about the same) 
noeud2 <- data.frame(x=noeud[,1], y=noeud[,2], z=noeud[,3])  # equal to the including NAs data 
NOEUD2 <- cbind(expand.grid(x=NOEUD$x, y=NOEUD$y), z=c(NOEUD$z)) # equal to the Akima interpolation data 
ggplot2::ggplot(noeud2, aes(x, y, z = z)) + geom_contour() 
lattice::contourplot(z ~ x * y, NOEUD2)