首先我先绘制了全国地图,代码如下:
download.file(url = "http://dapengde.com/r4rookies/cm.zip",
destfile = "c:/r4r/cm.zip")
unzip(zipfile = "c:/r4r/cm.zip", exdir = "c:/r4r")
require(rgdal)
cm <- readOGR(dsn = "c:/r4r/cm", layer = "bou2_4p")
plot(cm)
现在我想将单独的某一个省份地图(或者更加小的地级市县的地图)另存为一个shp文件,以及其它的小范围的道路、铁路等等数据也按地区另存为一个shp文件,方便各数据单独叠加绘制等!我现能够绘出全国地图及全国道路线路等,但想单独绘制小范围的却暂时没有头绪。请问各位见多识广前辈们,像这样的情况应该怎样进行处理呢?比如说,我想从全国地图中提取出江西省存为shp文件,从全国的铁路分布图中提取出江西省的道路存为shp文件,然后能直接调取江西省地图来进行绘制!