問題タブ [ncdf4]

関連する NCDF ファイルはこちら: https://www.ncdc.noaa.gov/paleo-search/study/19419

私はNCDFファイルを持っており、次のループを使用して、最初に各ファイルをCSVファイルとして保存し、次にシェープファイルとして保存しています:

ただし、上記のコードを実行すると、次のエラーが表示されます。

私は間違って何をしていますか? 何が起こっているのか理解できず、コードのさまざまな組み合わせや順列を試しました。上記のコードでファイル「cru_drought_4.csv」を作成する必要がありますか?

参考までに、完全なコードを次に示します。

R と ncdf4 パッケージを使用して、OPeNDAP サーバーからデータを取得しようとしています。ただし、nasa eosdis サーバーにはユーザー名とパスワードが必要です。Rを使用してこの情報を渡すにはどうすればよいですか?

これが私がやろうとしていることです:

そしてエラーメッセージ：

Rsx_nc4_get_vara_double のエラー: NetCDF: 認証失敗の構文エラー、予期しない WORD_WORD、SCAN_ATTR または SCAN_DATASET または SCAN_ERROR コンテキストが必要です: HTTP^ 基本: アクセスが拒否されました。Var: nlat Ndims: 1
Start: 0 Count: 400 ncvar_get_inner(d$dimvarid$group_id, d$dimvarid$id, default_missval_ncdf4(), のエラー: C 関数 R_nc4_get_vara_double がエラーを返しました

URL https://username:password@disc2 .... を試しましたが、それも機能しませんでした。

私はGISの目的でRを使用することにかなり慣れていません。複数の次元 (x、y、z、値、時間) を持ついくつかの変数を含む netcdf ファイルがあります。これをラスターブリックに変えようとしています。データは非常に大きいため、指定された時間枠と z (深さ) からデータを取得する必要があります。これは問題ではなく、次のコードを使用して適切な次元の配列を抽出します。

私のデータの例 (深さ/z と時間の次元を削除)

ncdf から抽出された緯度と経度は不規則な間隔で、それぞれ 2 次元です (つまり、各セルの不規則な間隔の緯度と経度)。

通常、このデータからラスター ブリックを生成します。

ただし、これは不規則な間隔を考慮しておらず、ラスター セルの値が間違った位置 (経度、緯度) に配置されています。スタック オーバーフローやその他の GIS ヘルプ ソースを検索しましたが、解決策に関する同様の問題が見つからないか、適切な質問をしていません。これについてどうすればいいのか特にわかりません。netcdf からデータを抽出するときにこれを処理する必要があるのか​​、範囲を定義せずにラスター ブリックを作成した後に処理する必要があるのか​​は不明です。ラスターの経度を定義する方法を見つけようとしましたが、運がありませんでした。経度、緯度、値を 3 列のデータフレームに変換してから、raster::rasterFromXYZ 関数を使用してみました。これは、私が扱っているデータのサイズに対して十分に速く動作しません。

どんな助けでも大歓迎です