raster 包

在 R 中可视化地理空间数据

Charlotte Wickham

Assistant Professor at Oregon State University

数据框不适合存储空间数据

head(preds)

        lon      lat predicted_price
1 -123.3168 44.52539        258936.2
2 -123.3168 44.52740        257258.4
3 -123.3168 44.52940        255543.1
4 -123.3168 44.53141        253791.0
5 -123.3168 44.53342        252002.4
6 -123.3168 44.53542        250178.7
  • 无 CRS 信息
  • 存储低效
  • 显示低效
在 R 中可视化地理空间数据

更好的栅格数据结构

  • 数据矩阵 + 网格信息 + CRS

栅格结构

在 R 中可视化地理空间数据

raster 包

  • sp 提供了一些栅格类：
    • SpatialGrid、SpatialPixels、SpatialGridDataFrame、SpatialPixelsDataFrame
  • 但 raster 更好：
    • 更易导入栅格
    • 大栅格不读入内存
    • 提供栅格运算函数
  • 也使用 S4，并在适当时提供相同函数
在 R 中可视化地理空间数据

raster 为 sp 对象提供打印方法

library(sp)
countries_spdf

An object of class "SpatialPolygonsDataFrame"
Slot "data":
                      name iso_a3 population         gdp                  region
1              Afghanistan    AFG   28400000    22270.00  
2                   Angola    AGO   12799293   110300.00  
3                  Albania    ALB    3639453    21810.00  
...
Slot "proj4string":
CRS arguments:
 +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0
在 R 中可视化地理空间数据

raster 为 sp 对象提供打印方法

library(raster)
countries_spdf

class       : SpatialPolygonsDataFrame 
features    : 177 
extent      : -180, 180, -90, 83.64513  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 ...
variables   : 6
names       :        name, iso_a3, population,         gdp, 
min values  : Afghanistan,    -99,        140,       16.00, 
max values  :    Zimbabwe,    ZWE, 1338612970, 15094000.00, ...
在 R 中可视化地理空间数据

¡Vamos a practicar!

在 R 中可视化地理空间数据

Preparing Video For Download...