我的工作目录的“数据”文件夹中有> 1000个netCDF文件。这些是显示每日雪覆盖程度的栅格,并在文件名中捕获了日期(yyyymmdd)。
nhtsd25e2_20120501_v01r01.nc
nhtsd25e2_20120502_v01r01.nc
nhtsd25e2_20120503_v01r01.nc
nhtsd25e2_20120504_v01r01.nc
nhtsd25e2_20120506_v01r01.nc
nhtsd25e2_20120507_v01r01.nc
nhtsd25e2_20120518_v01r01.nc
nhtsd25e2_20120520_v01r01.nc
nhtsd25e2_20120521_v01r01.nc
nhtsd25e2_20120522_v01r01.nc
当我打开一个单一文件时,发现数据字段的名称为'Merged.Snow.Cover.Extent'] >>
library(raster) snow <- raster("./data/nhtsd25e2_20120522_v01r01.nc") view(snow)
我想做的(并且不知道从哪里开始)是从“ data”文件夹中批量导入netCDF文件并创建一个RasterStack。
需要注意的是,我只想在每年的2月1日至7月1日之间导入netCDF文件以进行积雪(所以月日> =“ 0201”但<=“ 0701”,而年份无关紧要。
非常感谢您的帮助。
我的工作目录的“数据”文件夹中有> 1000个netCDF文件。这些是显示每日雪覆盖程度的栅格,并在文件名中捕获了日期(yyyymmdd)。 ...
使用list.files
获取文件名
使用list.files
获取文件名