r 语言 处理栅格数据

r 语言 处理栅格数据

作为新手，记录自己学到的一些东西，也希望能对需要的人有稍许帮助

常用的包

• raster #处理栅格数据
• rasterVis #栅格数据可视化
• RColorBrewer #颜色设置
• ncdf4 #nc数据的读取

常用的数据类型

• tiff
• nc

数据导入和简单的可视化

• tiff数据读取

  • 读取一个tiff文件 raster() 函数

    ## 导入使用的包
    library(raster)
    library(rasterVis)
    library(RColorBrewer)
    library(ncdf4)
    
    n_data_path <- dir('数据存放地址',full.names = TRUE)
    n_deposition <- raster(n_data_path[19])
    plot(n_deposition) # 使用基础的plot函数
    levelplot(n_deposition, margin = list(FUN = 'median'),
        par.settings=RdBuTheme) # 使用levelplot函数
    

    
    

  • 读取多个tiff文件 stack() 函数

    n_depo_all <- stack(n_data_path[c(1,6,11,16)])
    plot(n_depo_all) # 使用基础的plot函数
    levelplot(n_depo_all)  # 使用levelplot函数
    

  
  

• nc数据读取
  同样使用 raster() 和 stack() 函数

  sif_nc <- nc_open(dir('数据存放地址',full.names = TRUE)