使用 extract() 进行双线性插值:显示中间步骤

bilinear interpolation with extract(): show intermediate steps

我使用 extract() 从光栅对象中获取双线性插值点。有没有办法显示插值的中间步骤?我会对最近的四个网格点的坐标和值以及到我的插值点的距离感兴趣。

使用 r = raster(...) 和 spdf = SpatialPointsDataFrame() 我的函数调用如下:

out <- extract(r, spdf, method="bilinear")

您可以运行逐行raster:::.bilinearValue中的代码,raster = r; xy坐标 = coordinates(spdf); layer=1n=1