問題タブ [terra]

terraパッケージを使用して、ラスター スタックから抽出した値を Spatial オブジェクトの data.frame に追加したい

出力ファイルからわかるように、以前はパッケージの引数valを使用して取得していたものはありません。抽出されたdata.frameにフィールドを追加するにはどうすればよいですか?sprasterval

次の質問に対する回答を複製しようとしていますが、R Terra area コマンドでエラーが発生しています。

平射極図法を使用して衛星データからグリッド セル領域を推定する

私が実行している正確なコードは次のとおりです。

印刷 r 収量

area() コマンドの行に次のエラーが表示されます。

デフォルトの R インスタンスでこれを試し、クリーンなインスタンスを最初からセットアップしました。両方のインストールで同じエラーが発生します。どちらのシステムも macOS です。area コマンドを cellSize コマンドに変更しようとすると、R が完全に終了して再起動します。

各セットアップの sessionInfo() は次のとおりです。

システム 2: MacOS Big Sur のクリーン インストール

システム 1: macOS Catalina

私は R で作業しています。Google Earth Engine (NAIP 2018 NDVI) からダウンロードした解像度 60cm のラスター (.tif) タイルのディレクトリがあります。(ポリゴン データには人間を対象とする必要があるため、Google Earth Engine ではなく PC で解析を実行しています。) 52 個のタイルのサイズはそれぞれ 1.2 ～ 3.8 GB です。また、これらのラスターからゾーン平均を計算しようとしている 982 個のポリゴンもあります。私のコード (以下) は terra パッケージを使用しており、タイルを非常に大きな特異なラスターにモザイク化する代わりに、VRT (仮想ラスター) ファイルを作成することを選択しました。

このコードを Xeon Gold 6134 @ 3.2GHZ で実行しており、128GB の RAM を搭載しています。terraOptions() を何に設定しても、R はプロセッサまたは RAM の潜在能力のかなりの部分を使用することさえできません。

このコードでは、982 個のポリゴンすべてを実行するのに 11.8 日かかります。まだ試していないかもしれない特定のトリック/ツールを誰かが教えてくれたら、非常にありがたいです (たとえば、すべての terraOptions で作業してみました。raster パッケージと exact_extract パッケージを試しましたexact_extract()。 SpatRaster/VRT とsfポリゴン オブジェクトを入力として使用しているので、うまくいきません。非常に大きな特異なラスターをモザイク化するのを避けるためです。)

ありがとうございました。（データが大きすぎるか、人間の被験者に関連しているため、データを共有できないことをお詫びします...）ループされていないコードは次のとおりです。

編集：それぞれ1.2〜3.8GBの52タイル。150GB の総ディレクトリ サイズという最初の見積もりは、ArcGIS で圧縮されたサイズであったため、正しくありませんでした。

show(ras) は以下を生成します:

そのピクセル内にあるポイントを持つピクセルの重心の座標を抽出しようとしています。ここに再現可能な例があります

抽出された xy 座標は元の座標からわずかにずれているだけだと思いますが、代わりに値がすべて混同されているようです。コードに間違いがありますか、またはこれらの値を取得する別の方法はありますか?

app()家族の機能について簡単な質問があります

(n 個の値の) ベクトルをラスターで除算し、n 個の結果のラスタ​​ーを取得したい

同じことが SpatRasterDataset で行われることを知りたいです。それは可能ですか？

次に例を示します。