只是尝试过滤形状文件以简化绘图
我从英国政府下载了一个形状文件：
基于此：https://www.r-bloggers.com/r-and-gis-working-with-shapefiles/ https://www.r-bloggers.com/r-and-gis-working-with-shapefiles/
我写了但不知道过滤：
setwd("~/Documents/filename")
getwd() # --double confirm real data directory
#install.packages("maptools")
library(maptools)
crswgs84=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
ukmap=readShapePoly("filename.shp",proj4string=crswgs84,verbose=TRUE)
class(ukmap)
str(ukmap@data)
str(ukmap@polygons[[1]])
ukmap@bbox
# <-- need to do some filterig
plot(ukmap) # as this will take too long and not want to plot whole UK
例如我只想要“E06000001”到“E06000020”。
（文件名是“Local_Authority_Districts_December_2016_Full_Extent_Boundaries_in_the_UK”，不确定如何将其包含在程序编码引用中）
您可以考虑使用sf包来读取 shapefile 并绘制数据。过滤sf对象与过滤数据框相同。
library(sf)
# Read the sahpefile
ukmap <- st_read("Local_Authority_Districts_December_2016_Full_Extent_Boundaries_in_the_UK.shp")
# Subset the sf object
ukmap_sub <- ukmap[ukmap$lad16cd %in% c("E06000001", "E06000020"), ]
# Plot the boundary of ukmap_sub
plot(st_geometry(ukmap_sub))
如果您更喜欢与SpatialPolygonsDataFrame，我们可以使用readOGR函数从rgdal /questions/tagged/rgdal包裹。之后，我们可以对SpatialPolygonsDataFrame就像常规数据框一样。
library(maptools)
library(rgdal)
ukmap <- readOGR(dsn = getwd(), layer = "Local_Authority_Districts_December_2016_Full_Extent_Boundaries_in_the_UK")
ukmap_sub <- ukmap[ukmap$lad16cd %in% c("E06000001", "E06000020"), ]
plot(ukmap_sub)
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
