R:raster()在从SpatialGridDF转换时更改像素大小
我正在使用以下MODIS DSR 1km产品进行一些分析:MCD18A1.A2001001.h15v05.061.2020097222704.hdf 但是,我在从SpatialGridDataFrame转换为光栅时遇到了问题,因为像素大小会发生变化 这是我的剧本:R:raster()在从SpatialGridDF转换时更改像素大小,r,pixel,raster,projection,r-modis,R,Pixel,Raster,Projection,R Modis,我正在使用以下MODIS DSR 1km产品进行一些分析:MCD18A1.A2001001.h15v05.061.2020097222704.hdf 但是,我在从SpatialGridDataFrame转换为光栅时遇到了问题,因为像素大小会发生变化 这是我的剧本: Filename <- "ModisProductsOriginal/MCD18A1.A2001001.h15v05.061.2020097222704.hdf" SDSs <- getSds(File
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)
Filename是什么让你说像素大小改变了?我只看到有关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)
图书馆(terra)
你有没有检查我上传的图片->GMT\U 1200\U DSR\U sgdf和GMT\U 1200\U DSR\U rast?但是,我发现我的问题在于使用plot()时的可视化。谢谢你的帮助。不,我没有。你可以在文本中显示图像,但我不喜欢点击东西。