R语言入门——读取栅格数据参数解读

Posted 2023-01-25 统计学小王子

读取栅格数据

• 引言
• 1、数据读入
• 2、参数解读
• • 2.1 class （类）
  • 2.2 dimensions（维度）
  • 2.3 resolution （像素）
  • 2.4 extent（范围）
  • 2.5 names
  • 2.6 band(波段)
• 3、总结

引言

很多非地理专业的同学会在使用天气、地理、自然灾害、温度湿度等数据时遇到一种数据格式。.nc格式。也称为栅格数据，但是读入之后确读不明白他的数据结构。下面我以2019华为杯的数据为例，为大家分析一下这种数据的数据结构。非专业解读，仅供参考。

1、数据读入

# 载入包
> library(raster)
> library(ncdf4)
> library(rasterVis)
> library(lattice)

# 数据读取
r = stack("数据/data.nc")

2、参数解读

> r
class      : RasterStack 
dimensions : 89, 180, 16020, 1985  (nrow, ncol, ncell, nlayers)
resolution : 2, 2  (x, y)
extent     : -1, 359, -89, 89  (xmin, xmax, ymin, ymax)
crs        : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0 
names      : X1854.01.01, X1854.02.01, X1854.03.01, X1854.04.01, X1854.05.01, X1854.06.01, X1854.07.01, X1854.08.01, X1854.09.01, X1854.10.01, X1854.11.01, X1854.12.01, X1855.01.01, X1855.02.01, X1855.03.01, ...

2.1 class （类）

数据的类型，栅格数据

2.2 dimensions（维度）

dimensions : 89, 180, 16020, 1985  (nrow, ncol, ncell, nlayers)

数据的维度，（行数，列数、行函数乘以列数、层数）。上述这个数据共有89行，180列，16020个数据，1985层。
每一层可以通过列表索引取出类似于矩阵格式的数据。当然取出的类还是栅格数据。下面的命令进行查看。

> Data <- r[[1]]  # 取出第一层
> class(Data)
[1] "RasterLayer"
attr(,"package")
[1] "raster"
> dim(Data)  # 第一层的维度
[1]  89 180   1

下面画一下第一层的效果图：

plot(Data)

2.3 resolution （像素）

数据的像素resolution :2, 2 (x, y)，这个是数据的像素，决定x轴y轴的经纬度的间隔。这个数x和y间隔都是2。下个参描说具体怎么用。

2.4 extent（范围）

这个的含义应该很明确，-1, 359, -89, 89 (xmin, xmax, ymin, ymax),x是从-1到359间隔为2，即：-1，1，3，5，…。

2.5 names

这个参数是每一层的名字，我猜想可以通过每一层的名字索引试试：

> Data2 <- r$X1854.01.01
> class(Data2)
[1] "RasterLayer"
attr(,"package")
[1] "raster"
> dim(Data2)

确认过眼神，确实可以。

2.6 band(波段)

这个数据没有这个参数，没有感性的认识，咱也不知道。有知道的大神欢迎评论区留言学习。

3、总结

希望对一些刚刚上手的人有用。总结不专业，有问题评论区见。谢谢支持了。