# Thu Sep 10 10:01:24 2020 -
# 字符编码：UTF-8
# R 版本：R x64 4.0.2
# cgh163email@163.com
#个人笔记不负责任
# —— 拎了个梨🍐
rm(list=ls());gc()
.rs.restartR()
require(raster)
# Thu Sep 10 10:01:29 2020 --bind绑定，附加到某个对象----------------------------
# 将Spatial *对象绑定（附加）到单个对象中。 所有对象必须具有相同的向量类型基类（SpatialPoints，SpatialLines或SpatialPolygons）
#
p <- readRDS(system.file("external/lux.rds", package = "raster"))
names(p)
# [1] "ID_1" "NAME_1" "ID_2" "NAME_2" "AREA"
mersch <- p[p$NAME_2 == 'Mersch',]
plot(mersch)
diekirch <- p[p$NAME_2 == 'Diekirch',]
plot(diekirch)
remich <- p[p$NAME_2 == 'Remich',]
plot(remich)
remich$NAME_1 <- NULL
x <- bind(mersch, diekirch, remich)
plot(x)
data.frame(x)
# 并列展示：
par(mfrow = c(2, 2)) # 一页多图
plot(mersch, main = '1')
plot(diekirch, main = '2')
plot(remich, main = '3')
plot(x, main = '附加、绑定在一起')
# Thu Sep 10 10:27:59 2020 --DIY--用sp对象--------------------------
rm(list = ls());gc()
load('广州各区sf对象.Rdata')
# map.gz.xian <- sf::st_read('https://geo.datav.aliyun.com/areas_v2/bound/440100_full.json')
map1 <- as(map.gz.xian[2, ], 'Spatial')
str(map1)
map2 <- as(map.gz.xian[4, ], 'Spatial')
str(map2)
map.all <- cbind(map1, map2)
par(mfrow = c(3, 1)) # 一页多图
plot(map1, main = '1')
plot(map2, main = '2')
plot(map.all, main = '合并')
官方实例
DIY实例
