time <- ncvar_get(ncdata, 'time') data <- ncvar_get(ncdata, 'lwe_thickness') # 创建输出目录 output_dir <- "F:/transfer2/grace/nc2tif" dir.create(output_dir, showWarnings = FALSE) # 循环遍历每个月的数据,并将其输出为tif格式 for (i in 1:length(time)) { # 从时间戳中获取年份...
本篇记录之前作大作业时需要将nc数据按年合成tif的代码, library(raster) library(ncdf4) # 读取nc文件 tmpfile = nc_open("./raster/tmp_2020.nc") prefile = nc_open('./raster/pre_2020.nc') …
443 0 05:23 App R语言对年数据线性拟合(每年的均值获取,计算,拟合) 1248 1 08:23 App R语言处理GIMMS_NDVI数据(nc4转tif) 2266 7 57:25 App 零基础实现基于R的多元线性回归 612 0 07:02 App R下载与安装路径要求(搬于huang老师) 500 0 15:28 App arcgis出地图-行政区划分地图简单绘制 4088 1 03...
5.需要注意的是法2用raster包写出的时候没加“.tif”,但是最后输出的文件自动带了tif后缀。terra包写出的时候就必须加“.tif”后缀。 另外,在StackExchange上有人比较了raster和terra加载nc文件,那个提问者说他用terra加载一个很大的nc用了40多分钟 while raster包只用了一两秒,然后这个问题的回答说raster读取nc文件...
NC数据的批量读写 我们将处理好的2005-2019月均雪水当量数据导出保存,以便日后使用,使用paste构造相应路径即可。 #数据导出:将月均数据导出为nc再转为tif #创建维度 y<-ncdim_def('y',units = 'm',vals = lat ) x<-ncdim_def('x',units = 'm',vals = lon) ...
nc转多波段TIFF 上面Model Builder导出的数据将所有的时间写入了一个TIF,同时时间信息丢失了,没有具体的显示,也不太方便。 QGIS查看nc数据 QGIS原生支持nc数据,直接将nc数据拖入QGIS即可打开。颜色可以从图层属性,符号化选项卡进行修改 QGIS符号化调整 波段1时间是36905,时间起始是1900-1-1,可以在EXCEL里面进行一下...
国家地球系统科学数据中心下载的全国逐月最低气温数据集,每3年数据共36个月的最低气温数据保存在一个nc文件中,现在想要将数据导出为逐月栅格数据,以15-17年数据为例。 R包加载 主要用到ncdf4包读取nc文件,raster包将读入数据转栅格数据并输出. library(raster) library(ncdf4) library(rasterVis) library(lattice)...
par()函数中的mfcol、mfrow,用于设定图像设备的布局,参数形式为c(nr, nc),子图的绘图顺序是按列还是按行就分别根据是参数指定的是mfcol还是mfrow。 par(mfcol = c(3,3)) par(mfrow = c(3,3)) mfcol和mfrow设置绘图区域的布局有一个问题,就是只能将绘图区域分成大小相等的子区域,但是我们可能希望不同的子...
学员疑问为什么重命名后的栅格波段没有改变，预览图和temp变量名称不一致？temp变量，指定的波段名称预览中波段名为月中15或16日和每个月1号的波段名不一致代码原文library(terra)library(tidyverse)library(raster)library(gimms)#读取温度cru = rast("cru_ts4.06.1901.2021.tmp.dat.nc", "tmp")#中国矢量roi =...
writeRaster(tp, filename = paste0("./tif/", tp@cpp[["names"]], ".tif")) #输出多个为单波段tif,文件名为波段名 从下图我们就可以看出,之所以修改波段名称后预览没有随之改变,并不是波段名称修改失败了,而是预览plot查看的是对应的波段时间戳,而使用time函数修改时间戳后,plot的...