dữ liệu Interpolating trong R

Question

Giả sử tôi có một 3 5 ma trận trong R:

4 5 5 6 8 
3 4 4 5 6 
2 3 3 3 4 

Tôi muốn suy ở giữa những giá trị này để tạo ra một ma trận kích thước 15 bởi 25. Tôi cũng muốn xác định nếu nội suy là tuyến tính, gaussian, vv Làm cách nào tôi có thể làm điều này?

Ví dụ, nếu tôi có một ma trận nhỏ như thế này

2 3 
1 3 

và tôi muốn nó trở thành 3 bằng 3, sau đó nó có thể trông như

2 2.5 3 
    1.5 2.2 3 
    1 2  3 

Answer 1

app <- function(x, n) approx(x, n=n)$y # Or whatever interpolation that you want 

apply(t(apply(x, 1, function(x) app(x, nc))), 2, function(x) app(x, nr)) 
    [,1] [,2] [,3] 
[1,] 2.0 2.50 3 
[2,] 1.5 2.25 3 
[3,] 1.0 2.00 3 

Answer 2

thời gian dài trước đây Tôi đã viết một món đồ chơi tương tự, ngoại trừ tôi không bao giờ có xung quanh để xác định chức năng nội suy. Ngoài ra còn có raster::disaggregate.

zexpand<-function(inarray, fact=2, interp=FALSE, ...) { 
# do same analysis of fact to allow one or two values, fact >=1 required, etc. 
fact<-as.integer(round(fact)) 
switch(as.character(length(fact)), 
      '1' = xfact<-yfact<-fact, 
      '2'= {xfact<-fact[1]; yfact<-fact[2]}, 
      {xfact<-fact[1]; yfact<-fact[2];warning(' fact is too long. First two values used.')}) 
if (xfact < 1) { stop('fact[1] must be > 0') } 
if (yfact < 1) { stop('fact[2] must be > 0') } 
bigtmp <- matrix(rep(t(inarray), each=xfact), nrow(inarray), ncol(inarray)*xfact, byr=T) #does column expansion 
bigx <- t(matrix(rep((bigtmp),each=yfact),ncol(bigtmp),nrow(bigtmp)*yfact,byr=T)) 
# the interpolation would go here. Or use interp.loess on output (won't 
# handle complex data). Also, look at fields::Tps which probably does 
# a much better job anyway. Just do separately on Re and Im data 
return(invisible(bigx)) 
}