栅格的合并（马赛克）改变分辨率答案

【问题标题】：Merge (mosaic) of rasters changes resolution栅格的合并（马赛克）改变分辨率
【发布时间】：2021-08-05 05:37:46
【问题描述】：

我正在使用我开发的 R 脚本合并两个 MODIS DSR 磁贴，这些是产品： https://drive.google.com/drive/folders/1RG3JkXlbaotBax-h5lEMT7lEn-ObwWsD?usp=sharing

所以，我从同一日期 (2019180) 打开两个产品（平铺 h15v05 和平铺 h16v05），然后打开每个 SDS 并将它们合并在一起（来自 h15v05 的 00h 和来自 h16v05 的 00h 等等......）

两种产品的 Panoply 可视化（使用合并选项）：

紫色方块是分隔两个图块的分割线的位置。

通过我的代码，我获得了一个具有不同分辨率（以及不同的最小/最大值）像素的图，但我不明白为什么：

我怀疑得到的结果是由于：

1- 从 Sinusoidal CRS 更改为 longlat WGS84 CRS；

2- 使用重采样（方法 ngb）处理马赛克。

我的代码很广泛，但这里有一些部分：

# Open scientific dataset as raster
SDSs <- sds(HDFfile)
SDS <- SDSs[SDSnumber]
crs(SDS) <- crs("+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs")
SDSreprojected <- project(SDS, DesiredCRS)
SDSasRaster <- as(SDSreprojected, "Raster")

# Resample SDS based on a reference SDS (SDS GMT_1200_DSR of a first product), I need to do this to be able to use mosaic
SDSresampled <- resample(SDSasRaster,ResampleReference_Raster,method='ngb')

# Create mosaic of same SDS, but first convert stack to list to use mosaic
ListWith_SameSDS_OfGroupFiles <- as.list(StackWith_SameSDS_OfGroupFiles)
ListWith_SameSDS_OfGroupFiles.mosaicargs <- ListWith_SameSDS_OfGroupFiles
ListWith_SameSDS_OfGroupFiles.mosaicargs$fun <- mean
SDSmosaic <- do.call(mosaic, ListWith_SameSDS_OfGroupFiles.mosaicargs)

# Save SDSs mosaic stack to netCDF
writeRaster(StackWith_AllMosaicSDSs_OfGroupFiles, NetCDFpath, overwrite=TRUE, format="CDF", varname= "DSR", varunit="w/m2", longname="Downward Shortwave Radiation", xname="Longitude", yname="Latitude", zname="TimeGMT", zunit="GMT")

有谁知道结果之间不匹配的原因可能是什么？

打印(ResampleReference_Raster)

class      : RasterLayer 
dimensions : 1441, 897, 1292577  (nrow, ncol, ncell)
resolution : 0.01791556, 0.006942043  (x, y)
extent     : -39.16222, -23.09196, 29.99652, 40  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : memory
names      : MCD18A1.A2019180.h15v05.061.2020343034815 
values     : 227.5543, 970.2346  (min, max)

打印（SDSasRaster）

class      : RasterLayer 
dimensions : 1399, 961, 1344439  (nrow, ncol, ncell)
resolution : 0.01515284, 0.007149989  (x, y)
extent     : -26.10815, -11.54627, 29.99717, 40  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : memory
names      : MCD18A1.A2019180.h16v05.061.2020343040755 
values     : 0, 0  (min, max)

打印（SDSmosaic）

class      : RasterLayer 
dimensions : 1441, 897, 1292577  (nrow, ncol, ncell)
resolution : 0.01791556, 0.006942043  (x, y)
extent     : -39.16222, -23.09196, 29.99652, 40  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : memory
names      : layer 
values     : 0, 62.7663  (min, max)

另外，脚本忽略了一些岛屿（右下）...

【问题讨论】：

我怀疑是你打给resample() 的电话，不幸的是，我不能确定，因为你的代码在你发布的时候是不可重现的。您介意发布print(ResampleReference_Raster) 和print(SDSasRaster) 的输出吗？这将显示栅格的分辨率和范围。
关于投影，如果您尝试重新投影栅格，我相信正确的函数是projectRaster()。作为最佳实践，您可能需要考虑使用 WGS84 的语法 crs(SDS)<-crs("+init=epsg:4326") 使用其 EPSG 代码定义投影字符串
我将打印添加到帖子中，ResampleReference_Raster 和 SDSasRaster 具有不同的分辨率，因为我使用了重新采样。投影似乎是正确的，唯一的区别是 Panoply 中像素的倾斜度。
我认为问题在于程度。我使用的参考范围忽略了一些岛屿。我需要栅格具有相同的范围，因为我想堆叠它们并进行马赛克，我怎样才能获得包含两个图块的范围？

标签： r merge raster resolution mosaic

【解决方案1】：

对不起，我没有早点回复。所以我认为你是对的，这个问题是你重新采样的程度。我认为您可以通过创建一个虚拟栅格来解决此问题，该栅格具有您要重新采样的栅格范围，但具有您要镶嵌到的栅格的分辨率。尝试：

dummy<-raster(ext = SDSasRaster@extent, resolution=ResampledReference_Raster@res, crs=SDSasRaster@crs)
SDS2<-resample(SDSasRaster, dummy, method="ngb")
Final<-moasic(SDS2, ResampledReference_Raster, fun=mean)

【讨论】：