HDF ファイルから MODIS 土地被覆データを抽出し、ラスターに入れることができました。
library(gdalUtils) #? gdal_translate?
library(raster)
sds <- get_subdatasets(".......myfileplath/file.hdf")
# Isolate the name of the first sds
name <- sds[5]
filename <- '2009test.tif'
gdal_translate(sds[5], dst_dataset = filename)
# Load the Geotiff created into R
r <- raster(filename)
それをデータフレームに入れたいのですが、元の正弦波から再投影する+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defsと思います。
私が分析している他のデータセットと互換性のある通常の楕円体/ WGS84 に。
これは私が試したもので、うまくいくようです：
sr <- '+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0'
projected_raster <- projectRaster(r, crs = sr)
ただし、この新しい投影法で自分の土地被覆データをデータ フレームに入れると、すべての土地被覆値が NA になりました。
これは、データフレームが洞投影でどのように見えるかです (4 は土地被覆分類です)
DF <- as.data.frame(r, xy=TRUE)
head(DF)
# x y X2009test
#1 -1111718.9 1111719 4
#2 -1111255.6 1111719 4
#3 -1110792.2 1111719 4
#4 -1110328.9 1111719 4
#5 -1109865.6 1111719 4
私の再投影では、次のようになります。
reprojected_DF <- as.data.frame(projected_raster, xy=TRUE)
head(reprojected_DF)
# x y X2009test
#1 -10.173076 10.01876 NA
#2 -10.168896 10.01876 NA
#3 -10.164716 10.01876 NA
#4 -10.160536 10.01876 NA
#5 -10.156356 10.01876 NA
#6 -10.152176 10.01876 NA
#7 -10.147996 10.01876 NA
私が間違っていること、または土地被覆座標を適切に再投影する方法について何かアドバイスはありますか?
乾杯！！！！
また、NASA MODIS 再投影ツールがあったことを読みましたが、それはもはや存在しない/利用可能ではありません。それについて何か知っている人はいますか？