zwdbordeaux
对不起,刚才一激动,感叹号居然超过了3个,犯了谢老大的版规。自我反省一下。
“再激动感叹号也不要超过三个!!”
jah_et
[未知用户] 呵呵,这个不绝对
hunter99
建议补充上澳门的数据,否则不好在学术或正式场合引用。
tianping
很好,很和谐,很强大!!!
yzhizhi2006
不错不错,以前我做的时候用了很麻烦的笨方法,现在好了,语句简化但用途更强大了:)
finzhang
导入GIS数据后,发现省名在R里显示是乱码,请问应该如何解决? 谢谢!
zwdbordeaux
我也遇到了与bo一样的问题。应该是语言设置的问题,当我读入x的时候发现:
Shapefile type: Polygon, (5), # of Shapes: 925
Warning message:
use readShapeSpatial:
objects other than Spatial objects defined in the sp package are deprecated
之后也就没有办法对每个省的颜色进行改变!
正在想法解决
邱怡轩
[未知用户] R是中文的还是英文的?
DavidLung-thinkeeper
没有加载到GIS数据所在的工作空间
hgzhang-qgzhang
我试着运行程序,提示
> x<-read.shape('bou2_4p.shp');
错误: 没有"read.shape"这个函数
求解答,先谢了
邱怡轩
[未知用户] 新版本的maptools包把很多原来的函数都改掉了,如果你要用的话可能得重新读一下它的帮助。
yihui
[未知用户] 用readShapePoly()
函数。如x = readShapePoly('bou2_4p.shp')
邱怡轩得更新正文了。我刚看了一下,新的函数读进来的数据的子元素直接就是NAME
那一级的,原来的x$att.data$NAME
就是现在的x$NAME
,比如plot(x, col = rainbow(length(levels(x$NAME)))[as.integer(x$NAME)])
hgzhang-qgzhang
[未知用户] 谢谢,修改后就没问题了
池振合
plot(x,fg=gray(924:0/924))
错误于plot.window(xlim = xlim, ylim = ylim, asp = asp, ...) :
图形参数"fg"的长度不对
yihui
[未知用户] 参见你头上我的回复中的例子,这个包以及sp
包都更新了,现在应该用col
参数。
诸葛小小
[未知用户] 哈哈。太感谢了。终于成功了。我用的是2.10.1
诸葛小小
这个是怎么回事?
> getID=function(mapdata,provname)
+ {
+ index=mapdata$att.data$NAME %in% provname;
+ ids=rownames(mapdata$att.data[index,]);
+ return(as.numeric(ids));
+ }
> midchina=c("河南省","山西省","湖北省","安徽省","湖南省","江西省");
> plot(x,recs=getID(x,midchina),col="green",ol="white",xlab="",
+ ylab="");
警告多于50个(用warnings()来显示第一个到第五十个)
>
邱怡轩
[未知用户] 新版本的maptools
包改动很大,我已经把文章更新了一下,参见文中的红色部分。
诸葛小小
呵呵 getColor 好像找不到了
邱怡轩
[未知用户] 什么意思?getColor()
不是写在前面吗?只是需要把里面的mapdata$att.data$NAME
改成mapdata@data$NAME
。