首页 > 解决方案 > 如何创建每周光栅马赛克?

问题描述

我正在处理一些哥白尼数据(土壤水分产品,时间分辨率 - 1 天)。我已经下载了 2021 年的所有可用数据(365 幅图像的集合)。是否可以创建从连续 7 天创建的马赛克,然后在接下来的 7 天和接下来的 7 天形成等等?例如,1 月 1 日至 7 日,1 月 8 日至 14 日……老实说,我不知道该怎么做,所以只编写文件路径如下的代码:

library(raster)
library(rgdal)
library(maptools)

setwd("E:/igik/ssm")
getwd()
path<-getwd()

export<-file.path("E:/igik/ssm_mosaic")

r1 <- raster("E:/igik/ssm/c_gls_SSM1km-ssm_201901010000_CUSTOM_S1CSAR_V1.1.1.tiff")

r1
 
class      : RasterLayer 
dimensions : 771, 1280, 986880  (nrow, ncol, ncell)
resolution : 0.008928571, 0.008928571  (x, y)
extent     : 13.3125, 24.74107, 48.41071, 55.29464  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +no_defs 
source     : c_gls_SSM1km-ssm_201901010000_CUSTOM_S1CSAR_V1.1.1.tiff 
names      : c_gls_SSM1km.ssm_201901010000_CUSTOM_S1CSAR_V1.1.1 
values     : 0, 255  (min, max)

带有图像的文件夹(“E:/igik/ssm_mosaic”)包含 365 个单波段 tiff(1 个图像 - 1 天)。文件命名如下:

只有日期在变化。

标签: rimagegisraster

解决方案


推荐阅读