星期三, 四月 27, 2011

在R语言中进行地理信息绘图

数据可视化是将科学计算中产生的大量非直观的、抽象的或者不可见的数据,借助计算机图形学和图像处理等技术,以图形图像信息的形式,直观、形象地表达出来,并进行交互处理。地图是空间信息可视化的最主要和最常用的形式。

R语言中也有相关的扩展包可以使用,例如maps,sp等。下面的例子是仿照《R Graph Cookbook》一书中的例子,利用R中的sp空间分析包得到的。代码如下:



#首先加载空间分析包
library(sp)

#从gadm.org网站上得到中国的省区地理数据,并加载到R软件内存中
load(url("http://gadm.org/data/rda/CHN_adm1.RData"))

#从全国水资源公报上得到中国各省的人均水资源占有量,并存到变量中
gadm$water=c(1085,325,1473,3524,1079,2935,3989,2790,4147,358,2046,434
,1652,2490,451,3362,1467,871,2145,182,1000,12278,448,377,
182,1221,3135,152,4976,10000,5298,2005)

#利用空间绘图命令进行绘图
spplot(gadm,"water",col.regions = rev(terrain.colors(gadm$water)),main="中国人均水资源拥有量")

没有评论:

发表评论