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 にはあります。ただし、データ スロットはフラット テーブル用であり、すべてのラインにデータがあり、すべてのラインがポイントです。私が得られないのは、明らかにすべてのポイントについてこのデータを持っていないということです:それは私が最初にラスターから回復しようとしているものです! だから明らかに私が間違っていることがありますが、何がわかりません。
誰かが私が間違っていることを知っていれば、それは大きな助けになるでしょう!