Извлечение точек из растра с помощью MODIStsp_extract()
Я пытаюсь использовать пакет MODIStsp для загрузки данных MODIS и преобразования их в пригодный для использования формат (таблица R) для статистического анализа. Я загрузил данные и загрузил виртуальный растровый файл в моем сеансе R:
library(sp)
library(MODIStsp)
# loading virtual file
virtual_file = "C:/Users/hp-8570w/Desktop/internship/MODIS_data/output/Surf_Temp_Daily_1Km_v6/Time_Series/RData/Terra/LST_Day_1km/MOD11A1.RData"
data = get(load(virtual_file))
И сейчас я пытаюсь извлечь данные в некоторых конкретных точках. Для этого у меня есть набор точек (широты и долготы), хранящихся в кадре данных с 2 столбцами. Я создаю объект SpatialPoints:
points = read.csv('points.csv',sep = " ")
projection = CRS('+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181+b=6371007.181 +units=m +no_defs')
spatial_points = SpatialPoints(points, proj4string = projection)
И я хочу использовать функцию MODIStsp_extract():
test_extract = MODIStsp_extract(data,spatial_points)
Я получаю следующее сообщение об ошибке:
Error in MODIStsp_extract(raster_ts, spatial_points) :
no slot of name "data" for this object of class "SpatialPoints"
Как показывает документация ( https://www.rdocumentation.org/packages/sp/versions/1.3-1/topics/SpatialPoints), у объектов SpatialPoints нет поля "данных", но у SpatialPointsDataFrame есть. Однако слот данных предназначен для плоских таблиц, с некоторыми данными для каждой строки, каждая строка является точкой. Чего я не понимаю, так это того, что у меня, очевидно, нет этих данных для каждой точки: именно это я и пытаюсь восстановить из растра! Очевидно, что я что-то не так делаю, но не вижу, что.
Если кто-нибудь знает, что я делаю неправильно, это очень поможет!