我有一个单独的shapefile,其中包含美国人口最多的30个城市的多边形边界。我想循环遍历R中的shapefile,并在一个5 x 6的网格中绘制每个城市的地图,每个网格的顶部都有它们的名字。每一行都有一个地名字段和一个UniqueID (0:29),以及其他字段,比如人口、SQMILES等。我可以轻松地设置shapefile的子集并绘制单个城市。然而,我正在努力处理for循环。任何提示都将受到感谢,特别是如果您还想在一个恒定的边界框中渲染每个城市的轮廓线,以便每个城市的各自边界彼此成比例。
该数据集的分辨率为0.5x0.5度/纬度,覆盖整个地球,并以.NC文件的形式下载,我使用以下命令提取了R中的数据:f <- open.ncdf("D:/CRU/cru_ts3.21.1961.1970.tmp.dat.nc")B <- get.var.ncdf(nc=f,varid="lon")
C <- g