• 综合主站
  • R绘制中国地图,并展示流行病学数据

从图片制作shp过程在GIS专业里属于地图数字化的一部分,一般使用ArcGIS等专业GIS软件,可以很好解决多边形边界咬合分离、添加投影、地图几何配准的问题,而且提取多边形外边界,线转面等都很容易实现。可以考虑shp制作在ArcGIS等专业软件来实现。
在地图上标识名称的代码可能有问题
tmp = coordinates(Changsha)
tmp = as.data.frame(tmp)

Error in data.frame(list(c(113.00043, 112.9974, 112.99068, 112.98441, :
arguments imply differing number of rows: 627, 791, 470, 201

我的R环境是R 3.1.0
取坐标可能要用
Changsha$CENTROID_Y
Changsha$CENTROID_X
其实R做地图并没什么优势,写一大堆的代码,还不如geoda或Arcgis的一个菜单管用。甚至也不如SAS命令来的简单,喜欢钻研R的可以当做消遣看看,实用主义者不建议用R这么复杂。
这篇文章确实是新鲜出炉,今年所写,上半年还能下载地图。
后来地理中心开始更新数据: http://ngcc.sbsm.gov.cn/article/xwzx/xmjs/201406/20140600004739.shtml
同时 ,也有网友询问,官方也解释了由于数据更新,现在不能下载: http://ngcc.sbsm.gov.cn/article/khly/lyzx/
如果有官方最新地图,还是用官方的吧。
gadm的地图不知道在国界细节上有没有出入。
6 天 后
很好的文章,说理清楚明白,代码可重复性也很高,赞一个...
5 天 后
没研究出imageJ提取地区界的方法,可有人知道?
4 个月 后
为什么我会出现这样的错误
> mysh = fortify(Changsha, region = "NAME99")
Error: isTRUE(gpclibPermitStatus()) is not TRUE
我运行到这一步就出现这样的错误了
我知道了,是"NAME99"的错误,为什么能解释下么?我把他删除也可以做成
4 个月 后
问题同上,哪位大神可以来棒棒忙吗?
4 个月 后
错误原因是maptools包依赖的gpclib不是开源软件,故无法使用。解决方法是安装先编译安装最新版geos软件(开源),然后在R中安装rgeos包,再重新安装maptools。希望对你有帮助!
3 个月 后
大神好,通过你的描述,我现在可以使用R绘画基本的地理地图,非常非常感谢。
有个问题想请教下,我按照Kriging包的例子,画了一个Kriging预测图,这个预测图没有相应的省份,如何能够给Kriging预测图加上省份边界呀,谢谢。
4 个月 后
[未知用户] 知不知道如何画一个英文中国地图,把这几个省用灰色表示出来即可Heilongjiang ,Jiangsu, Zhejiang, Jiangxi, Shandong,
Hubei ,Sichuan, Gansu 。。。。*(所有省名均用英文标注出来)
6 个月 后
[未知用户] 是的,按照你的方法已经解决
1 个月 后
请问我根据地图(地图shp数据就是来源于文中提到的国家基础地理信息中心的数据。)获取了各省的几何中心的地理坐标,如何转成经纬度呢?
1 个月 后
你好~你的这个问题解决了么?有一个上路新手。。
3 个月 后
请问出现这种情况是什么原因,如何解决?
Error in getinfo.shape(filen) : Error opening SHP file