Can';t更改R中空间多边形框架的坐标
这是我的节目:Can';t更改R中空间多边形框架的坐标,r,spatial,R,Spatial,这是我的节目: library(sp) library(RColorBrewer) #get spatial data for Spain on region level con <- url("http://gadm.org/data/rda/ESP_adm1.RData") print(load(con)) close(con) # plot Spain with colors col = c("blue","blue","blue","blue","blue","blue","bl
library(sp)
library(RColorBrewer)
#get spatial data for Spain on region level
con <- url("http://gadm.org/data/rda/ESP_adm1.RData")
print(load(con))
close(con)
# plot Spain with colors
col = c("blue","blue","blue","blue","blue","blue","blue","blue","blue",
"blue","blue","blue","blue","blue","blue","blue","red","red")
spplot(gadm, "NAME_1", col.regions=col, main="Regiones de España",
colorkey=F, lwd=.4, col="black")
我尝试进行更改,例如:
gadm[14,]@polygons[[1]]@Polygons[[1]]@coords
但我得到的信息是:
没有方法将这个S4类强制为向量
也可以尝试sf
package-CRAN链接
GADM数据URL现在似乎已被弃用,但您可以从手动下载包含SpatialPolygonsDataFrame for Spanish的RDS文件
库(sf)
西班牙尝试gadm[14,]@polygons[[1]]@polygons[[1]]@coords[]
gadm[14,]@polygons[[1]]@Polygons[[1]]@coords
library(sf)
spain <- readRDS(file = "ESP_adm1.rds") # downloaded from http://gadm.org/
# convert from sp to sf object
spain_sf <- st_as_sf(spain)
# Change the coordinates of Canary Islands - the 14th geometry of spain_sf,
# itself with 20 polygons, check str(st_geometry(spain_sf)[[14]])
# Will shift them closer to mainland Spain.
# The added values are degrees because the CRS is unprojected
# (e.g. add 10 to longitudes and 6 to latitudes)
for (i in 1:length(st_geometry(spain_sf)[[14]])) {
st_geometry(spain_sf)[[14]][[i]][[1]][, 1] <- st_geometry(spain_sf)[[14]][[i]][[1]][, 1] + 10
st_geometry(spain_sf)[[14]][[i]][[1]][, 2] <- st_geometry(spain_sf)[[14]][[i]][[1]][, 2] + 6
}
# plot to check
par(omi=c(0,0,0,2)) # make some space to the right for the legend
plot(spain_sf["NAME_1"],
axes = TRUE,
graticule = st_crs(spain_sf))
# dev.off() # resets par()