5

我有不均匀定位的图像样本，并且想插入到规则网格，因为（除其他外）大多数图像图形函数都需要规则网格。我注意到有一些 MatLab 函数（例如，请参见随机像素的图像插值）显然可以做到这一点，但找不到可以做到的 R 包。
这是一个简单的例子。 

#make up some 2D func
y<-matrix(rep(1:10,10) -.5 + runif(100),nrow=10)
x<-matrix(rep(1:10,10) -.5 + runif(100),nrow=10)
inmat<-sin(x) + cos(y)

所以 的值inmat在随机位置。我想要某种outmat<-interpolate(inmat,x,y,gridx,gridy)函数，其中inmat,x和y都是矩阵或所有向量（展开矩阵）。

我还看到 SciPy 有http://docs.scipy.org/doc/scipy/reference/generated/scipy.interpolate.interp2d.html可以做到这一点。包中是否有这样的功能，R还是我需要移植SciPy或MatLab编码？

4

1 回答 1

1

链接页面提供了指向R执行克里金法或其他插值函数的大量软件包的指针。

我发布我的个人选择作为答案只是为了结束这个问题。

我发现akima::interp这是一个简单的函数，可以对任意样本位置集合进行 2D 插值。
这并不意味着它对每个人都是最好的，我猜那些使用地理数据的人可能更喜欢那些设计用来处理特定地理调查相关文件类型和纬度/经度坐标系的软件包。

于 2013-10-01T11:20:38.657 回答