R与GIS实践系列-Shapefile导入与地图显示

林场

<div id="cnblogs_post_body">1. R项目简介

R是一个开源的统计计算和图形编程语言和软件环境，广泛应用于统计分析和数据挖掘[1]。R项目遵循GNU协议，它的软件环境源码由C,Fortran和R写就。R可以运行在多个平台，包括各种Unix，Linux发行版本，Window和MacOS上, 国内有中科院、厦门大学、北京交通大学等多个镜像[2]。R是对S语言的一种实现，由Ross Ihaka 和 Robert Gentleman 所创建，R取之他们名字的首字母。R具有由用户贡献的大量的类库，能够处理各种科学计算的问题。
2.R与GIS

与地学问题相关的R类库也十分繁多，包括Spatial data,Maptools,RGdal等。 不过似乎StackOverflow的主站关于R在GIS的问题比GIS分站更为活跃。原因未知……之所以要使用R是因为R能够以Hadoop集成[3],希望利用R更高效地处理空间聚类问题。虽然Mashout也是一个很好的解决方案，但是感觉用户更加需要脚本式的问题处理工具。R的GIS的书籍比较少，基本都是以文档的形式存在，其中Applied spatial data analysis with R是一本挺好的教材。多使用Example命令和Google也是学习R的一个好途径。
3.Shapefile导入与地图显示

R中导入Shapefile是非常容易的，可以使用的包，包括RGdal,maptools,PBSmapping等。在这里我们采用maptools和sp两个类库来实现地图数据的导入和显示。在R中心选择镜像并下载。
<ol> 加载这两个类库，输入如下代码：<div class="cnblogs_code">library(maptools)library(sp)