Question

假设我在R中有一个3乘5的矩阵：

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

我想在这些值之间进行插值以创建一个15乘25的矩阵。我还想指定插值是线性的，高斯的等等。我该怎么做？

例如，如果我有一个像这样的小矩阵

2 3
1 3

我希望它变为3乘3，然后它可能看起来像

  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

很久以前我写了一个类似的玩具，除了我从来没有定义插值函数。还有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))
}

R中的插值数据

2 个答案: