微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

将SpatVector转换为SF或SP

如何解决将SpatVector转换为SF或SP

是否有直接的方法将多边形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”。希望它对您有用。

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。