问题描述
是否有直接的方法将多边形SpatVector
类对象(来自terra
库)转换为简单要素或SpatialpolygonsDataFrames
?
谢谢
解决方法
有一种直接方法将SpatVector
强制转换为sp
对象。您可以从sp
对象创建一个sf
对象。
示例数据
library(terra)
f <- system.file("ex/lux.shp",package="terra")
v <- vect(f)
v
#class : SpatVector
#geometry : polygons
#elements : 12
#extent : 5.74414,6.528252,49.44781,50.18162 (xmin,xmax,ymin,ymax)
#coord. ref. : +proj=longlat +datum=WGS84 +no_defs
#names : ID_1,NAME_1,ID_2,NAME_2,AREA
创建一个sp
对象
x <- as(v,"Spatial")
x
#class : SpatialPolygonsDataFrame
#features : 12
#extent : 5.74414,ymax)
#crs : +proj=longlat +datum=WGS84 +no_defs
#variables : 5
#names : ID_1,AREA
#min values : 1,Diekirch,1,Capellen,76
#max values : 3,Luxembourg,12,Wiltz,312
并将sp
对象强制为sf
对象
library(sf)
y <- as(x,"sf")
y
#Simple feature collection with 12 features and 5 fields
#geometry type: POLYGON
#dimension: XY
#bbox: xmin: 5.74414 ymin: 49.44781 xmax: 6.528252 ymax: 50.18162
#CRS: +proj=longlat +datum=WGS84 +no_defs
#First 10 features:
# ID_1 NAME_1 ID_2 NAME_2 AREA geometry
#1 1 Diekirch 1 Clervaux 312 POLYGON ((6.026519 50.17767...
#2 1 Diekirch 2 Diekirch 218 POLYGON ((6.178368 49.87682...
(...)
,
这是一种解决方法(https://www.gitmemory.com/issue/rspatial/terra/71/660682352)。您可以使用相同的名称,但可以使用“ sp”代替“ sf”。希望它对您有用。