0

次の MODIS DSR 1km 製品を使用して分析を行っています: MCD18A1.A2001001.h15v05.061.2020097222704.hdf

ただし、ピクセル サイズが変更されるため、SpatialGridDataFrame から Raster への変換に問題があります...

これが私のスクリプトです：

Filename <- "ModisProductsOriginal/MCD18A1.A2001001.h15v05.061.2020097222704.hdf"
SDSs <- getSds(Filename)
GMT_1200_DSR_sgdf <- readGDAL(SDSs$SDS4gdal[8], as.is = TRUE)
GMT_1200_DSR_rast <- raster(GMT_1200_DSR_sgdf)

警告は次のとおりです。

1: In getProjectionRef(x, OVERRIDE_PROJ_DATUM_WITH_TOWGS84 = OVERRIDE_PROJ_DATUM_WITH_TOWGS84, : 破棄されたデータム Proj4 定義では指定されていません (Clarke 1866 回転楕円体に基づく): +proj=sinu +lon_0=0 +x_0=0 +y_0=0 +ellps=clrk66 +units= m +no_defs

2: in showSRID(uprojargs, format = "PROJ", multiline = "NO", prefer_proj = prefer_proj) : Proj4 定義の Clarke 1866 楕円体に基づく破棄されたデータ Unknown

結果は次のとおりです。

GMT_1200_DSR_sgdf

GMT_1200_DSR_rast

この問題の解決策が見つかりません。これが起こるべきかどうかを知りたいのですが、そうでない場合はどうすれば解決できますか? CRS を指定しますか? 正弦投影 ("+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs") を使用すると、警告のみが表示されます。おそらく、CRS を正しく使用していない可能性があります。 .

4

1 に答える 1

0

ピクセルサイズが変わると言う理由は何ですか? CRS に関する警告のみが表示されます。この警告は、PROJ ライブラリの新しいバージョンが WGS84 以外のデータを持つ PROJ 文字列を好まないためです。

サブデータセットを含む HDF は、 での処理がはるかに簡単terraです。

library(terra)
f <- "ModisProductsOriginal/MCD18A1.A2001001.h15v05.061.2020097222704.hdf"
s <- sds(f)
s8 <- s[8] 

# or
r <- rast(f)
于 2021-04-30T19:53:00.630 に答える