билинейная интерполяция с экстрактом(): показать промежуточные этапы

Question

Я использую extract(), чтобы получить билинейные интерполированные точки из моего растрового объекта. Есть ли способ показать промежуточные этапы интерполяции? Меня интересовали бы координаты и значения ближайших четырех точек сетки и расстояния до моей точки интерполяции.

С г = растр (...) и SPDF = SpatialPointsDataFrame() мой вызов функции выглядит следующим образом:

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

Answer 1

Вы можете запустить код в raster:::.bilinearValue построчно, с raster = r; xyCoords = coordinates(spdf); layer=1, n=1