interpolację danych w R

Question

Załóżmy, że ma trzy przez 5 macierz R:

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

chciałbym interpolacji między tymi wartościami w celu utworzenia macierzy o rozmiarze 15 przez 25. to również tutaj określić, czy interpolacja jest liniowa, gaussowska itd. Jak mogę to zrobić?

Na przykład, jeśli mam małą matrycę jak ten

2 3 
1 3 

i chcę, aby stało się 3 na 3, to może wyglądać

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

Dawno temu Napisałem podobną zabawkę, z wyjątkiem tego, że nigdy nie byłem w pobliżu, by zdefiniować funkcję interpolacji. Jest też 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)) 
}