用 R 画中国分省市地图
用R 也可以做出漂亮的依参数变化的中国地图。
主要参考(http://cos.name/author/yixuanq/)
library(maptools)
axx<-readShapePoly("bou2_4p.shp", proj4string=CRS("+proj=longlat +ellps=clrk66"))
x=axx
#plot(x,col=gray(924:0/924));#using col instead of fg!
getColor=function(mapdata,provname,provcol,othercol)
{
f=function(x,y) ifelse(x %in% y,which(y==x),0);
colIndex=sapply(mapdata$NAME,f,provname);
col=c(othercol,provcol)[colIndex+1];
return(col);
}
provname=c("北京市","天津市","河北省","山西省","内蒙古自治区",
"辽宁省","吉林省","黑龙江省","上海市","江苏省",
"浙江省","安徽省","福建省","江西省","山东省",
"河南省","湖北省","湖南省","广东省",
"广西壮族自治区","海南省","重庆市","四川省","贵州省",
"云南省","西藏自治区","陕西省","甘肃省","青海省",
"宁夏回族自治区","新疆维吾尔自治区","台湾省",
"香港特别行政区");
pop=c(1633,1115,6943,3393,2405,4298,2730,3824,1858,7625,
5060,6118,3581,4368,9367,9360,5699,6355,9449,
4768,845,2816,8127,3762,4514,284,3748,2617,
552,610,2095,2296,693);
provcol=rgb(red=1-pop/max(pop)/1,green=1-pop/max(pop)/1,blue=1/1.5);
plot(x,col=getColor(x,provname,provcol,"white"),border = "grey",xlab="",ylab="")
最新文章
- Lua 学习笔记(九)协同程序(线程thread)
- CSS3中的动画功能(二)
- MySQL Cluster搭建与测试
- Jsp字符编码过滤器
- Eclipse插件推荐:UCDetector: Unnecessary Code Detector
- 开发移动app与服务器端session的状态管理与交互
- Android之EditText
- Oracle系列之表空间
- PHP手册 – 变量和常量
- [PCB设计] 4、BAT脚本处理AD生成的GERBER文件为生产文件
- C# orderby子句
- httpie的使用
- 移动端 Retina屏border实现0.5px
- iOS将excel转plist
- GCC编绎详解
- 数值和字符串相互转换(C++ 数据类型转换技巧)
- Ionic上滑刷新
- iOS学习笔记(九)—— xml数据解析
- c++ sort用法 学习笔记
- Akka简介与Actor模型
热门文章
- DEV开发之控件XtraTabbedMdiManager
- 【leetcode刷题笔记】Convert Sorted List to Binary Search Tree
- C语言伪随机数的注意事项
- delphi的获取某坐标的颜色值
- Java -- 容器使用 Set, List, Map, Queue, Collections
- Qt中 QTableWidget用法总结
- Codeforces 461B Appleman and Tree:Tree dp
- DBGrideh 实现自动排序
- UTF-8 delphi 函数
- Java_io_02_从一个目录拷贝文件到另一个目录下