2d Bilineare interpolation in R

Habe ich einige mehrdimensionale Datenstrukturen, die ich brauche, um einige der interpolation auf. Ich bin wirklich finden es schwer zu bekommen ein Beispiel, wie dies in R!

Als ein Beispiel, wenn ich das array old:

old <- array(runif(10*12), dim=c(12,10))

str(old)
num [1:12, 1:10] 0.763 0.429 0.792 0.923 0.476 ...

was ich tun möchte, ist das array new:

new <- interp2d(old, newx=6, newy=5)

also ich will zu ändern, die Dimensionen des Arrays, so dass neue ist ein 6*5 array der old Daten -- in diesem Beispiel wäre es wichtig zu erhalten die Ergebnisse in der old raster. Die Zeile oben ist ein Beispiel, was ich gerne tun würde, ich weiß nicht, wie und hoffte, jemand würde??! Danke!

  • Ein ??interpolate zeigt (für mich) stats::approx (in base-R), spatstat::interp.im und fBasics::linearInterp und e1071::interpolate. Vielleicht ist einer dieser arbeiten konnte? (Ich hätte Interesse an einem base-R-Lösung ohne extra-Pakete, aber stats::approx sieht ein wenig schwerfällig, dass gibt es eine Liste).
  • Dies ist sehr ähnlich zu stackoverflow.com/questions/9171904/...
InformationsquelleAutor Alex Archibald | 2012-02-10
  1. 2

    Hier ist ein Anfang, aber nicht ich bin mir sicher, dass es das macht, was Sie wollen (erhalten Summen). In der Tat, ich bin nicht sicher, ob ich sehen, wie können Sie erhalten genau die Summen, wenn man die Verrohung in dieser Art und Weise ... es sei denn, Sie möchte lieber "bin" als "interpolieren"?

    library(sos)
findFn("{bilinear interpolation}")

set.seed(101)
old <- array(runif(10*12), dim=c(12,10))

library(fields)

interp2d <- function(old, newx, newy) {
  interp.surface.grid(list(x=seq(nrow(old)),y=seq(ncol(old)),z=old),
                      list(x=seq(1,nrow(old),length=newx),
                           y=seq(1,ncol(old),length=newy)))$z
}

newmat <- interp2d(old, newx=6, newy=5)
    • Danke Ben!!! Dies scheint zu tun, was ich brauche. Seine eine Schande, es erfordert etwas anderes als das Basis-Paket (D. H. Felder), aber es funktioniert, also bin ich glücklich!
    InformationsquelleAutor Ben Bolker
Schreibe einen Kommentar