全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
22014 21
2012-03-29
悬赏 2000 个论坛币 已解决
> source("map2.R")
错误于getinfo.shape(filen) : Error opening SHP file

R脚本文件:
library(maptools); x=readShapePoly('bou2_4p.shp');
plot(x);
getColor=function(mapdata,provname,provcol,othercol) {
     f=function(x,y) ifelse(x %in% y,which(y==x),0);
    colIndex=sapply(mapdata@data$NAME,f,provname);
    col=c(othercol,provcol)[colIndex+1];  return(col);
}
provname=c("北京市","天津市","上海市","重庆市");
provcol=c("red","green","yellow","purple");
plot(x,col=getColor(x,provname,provcol,"white"))

最佳答案

epoh 查看完整内容

1.如果三个文件都放在Working Directory(getwd()) x=readShapePoly('bou2_4p.shp'); 2.如果三个文件是放在 C:\....\R\win-library\2.15\maptools\shapes x=readShapePoly(system.file("shapes/bou2_4p.shp",package="maptools"))
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2012-3-29 13:13:26
耕耘使者 发表于 2012-4-1 20:07
悬赏额涨到2000元了。
1.如果三个文件都放在Working Directory(getwd())
  x=readShapePoly('bou2_4p.shp');

2.如果三个文件是放在
  C:\....\R\win-library\2.15\maptools\shapes
  x=readShapePoly(system.file("shapes/bou2_4p.shp",package="maptools"))
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-3-29 16:59:26
奇怪的是,原来能够成功。现在却不行了,估计是程序包有了改动。
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-3-29 22:30:28
继续求助
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-3-29 22:40:18
可以的哦
> source("map2.R")
载入需要的程辑包:foreign
载入需要的程辑包:sp
载入需要的程辑包:lattice
Checking rgeos availability: FALSE
        Note: when rgeos is not available, polygon geometry     computations in maptools depend on gpclib,
        which has a restricted licence. It is disabled by default;
        to enable gpclib, type gpclibPermit()
附件列表
L@DJ~713H@IP[5%9$`8L@WB.jpg

原图尺寸 39.6 KB

L@DJ~713H@IP[5%9$`8L@WB.jpg

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2012-3-30 12:10:12
以前是可以的,我的电脑上也可以。但笔记本上不行。
新装的程序包不行,可能是更新了
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群