其他
干货分享 | R-地理数据处理基础-NC格式数据读取，可视化及导出为GeoTiff文件
The following article is from TheWhoOPs Author 申泽西
目录
0.问题导入 1.示例数据 2.将NC文件中的FSDS变量导入为raster 3.数据可视化检查(图2) 4.结果写入GeoTiff文件 5.总结 6.解决本示例问题所采用的R-packages 7.致谢
0. 问题导入
目前，很多地理数据均存储为NetCDF格式（下简称NC格式），而往往内含多变量，不方便处理。其中，一种解决思路就是将NC格式文件中的多个变量数据转格式为单变量的GeoTiff文件（图1），进而进一步处理。1. 示例数据
本示例数据格式为NC格式数据，内含7个变量。
names(ncfile$var)
[1] "EDGEE" "EDGEN" "EDGES" "EDGEW" "FSDS" "LATIXY" "LONGXY"
2. 将NC文件中的FSDS变量导入为raster
input_nc = 'L:\\JianShu\\2019-12-4\\data\\clmforc.Qian.c2006.T62.Solr.1972-01.nc'varname = 'FSDS'
nc2raster = raster(input_nc, varname = varname)
3. 数据可视化检查(图2)
png("L:\\JianShu\\2019-12-4\\plot1.png",height = 15,
width = 20,
units = "cm",
res = 1000)
print(levelplot(nc2raster))
dev.off()
4. 结果写入GeoTiff文件
filename exists; use overwrite=TRUEError in .getGDALtransient(x, filename = filename, options = options, :
filename exists; use overwrite=TRUE
writeRaster(nc2raster, output, format = 'GTiff', overwrite = TRUE)
5. 总结
本篇主要解决了以下几个问题：
6. 解决本示例问题所采用的R-packages
library(raster)library(ncdf4)
library(rasterVis)
library(lattice)
7. 致谢
首先，感谢大家的持续关注，小编会继续努力，持续更新下去的！大家如果觉得有帮助啊，还麻烦大家关注点赞，也可以扩散到朋友圈，多多引导朋友加入咱们这个简书技术平台, 代码共享推动科研进程, 多谢大家啦～大家如果在使用本代码的过程有遇到问题的，可以留言评论，也可以私信我哈~~▼ ▼ ▼ ▼ ▼ ▼ ▼ ▼ ▼
转载自微信公众号TheWhoOPs
经作者授权转载
文章仅代表作者观点，与本公众号无关，版权归原作者所有
原文标题：R-地理数据处理基础-NC格式数据读取，可视化及导出为GeoTiff文件
图文排版：申茜茜
责任编辑：张英浩
审核：王波涛
终审：顾伟男 田巍 梁龙武
猜你喜欢
1. 干货分享 | R-如何基于并行计算根据extent生成shp，并将一个list的shp文件塌缩为一个shp？
2. 干货分享 | R-ggplo2 如何在空间图中绘制统计子图
4. 干货分享 | R-ggplot2中坐标轴在绘图过程中出现乱序的问题
扫描二维码，关注我们
都是成年人了，要记得告诉地小联你也“在看”哦~👇