读取 shapefile
Shapefile 是共享空间数据最常见的格式之一，可以使用 rgdal 包中的 readOGR() 轻松读入 R。readOGR() 有两个重要参数：dsn 和 layer。具体传什么要看您读入的数据类型。您在视频中学到，对于 shapefile，dsn 应为包含该 shapefile 各个组成文件的目录路径，而 layer 则是该 shapefile 的文件名（不含扩展名）。
为了制作地图，您需要社区（neighborhood）边界。我们已从纽约市规划局的 Open Data Platform 下载了"Neighborhood Tabulation Areas"。下载的是 zip 压缩包，我们已将解压后的内容放在您的工作目录中。
您将使用 base R 的 dir() 函数查看工作目录内容，然后将该 shapefile 读入 R。
本练习是课程的一部分
在 R 中可视化地理空间数据
练习说明
- 不带参数使用
dir()，找出 shapefile 所在目录的名称。 - 对该目录路径调用
dir()，查看其中包含的文件。 - 现在您已经知道目录和文件名。使用
readOGR()将社区 shapefile 读入为名为neighborhoods的对象。 - 调用
summary()检查neighborhoods的内容。 - 通过绘图检查：绘制
neighborhoods。
交互式实操练习
通过完成这段示例代码来试试这个练习。
library(sp)
library(rgdal)
# Use dir() to find directory name
# Call dir() with directory name
# Read in shapefile with readOGR(): neighborhoods
# summary() of neighborhoods
# Plot neighborhoods