我有一个简单的两列矩阵,我想将其转换为一个sf
对象,其中每一行指定一个点:
> set.seed(123);m=matrix(runif(10),ncol=2)
> m
[,1] [,2]
[1,] 0.2875775 0.0455565
[2,] 0.7883051 0.5281055
[3,] 0.4089769 0.8924190
[4,] 0.8830174 0.5514350
[5,] 0.9404673 0.4566147
最简单的方法行不通,因为apply
将点重新组合成矩阵,操作只是一个非常慢的转置函数:
> apply(m,1,st_point)
[,1] [,2] [,3] [,4] [,5]
[1,] 0.2875775 0.7883051 0.4089769 0.8830174 0.9404673
[2,] 0.0455565 0.5281055 0.8924190 0.5514350 0.4566147
如果不使用显式循环,最好的方法是使用以下代码:
> st_sfc(lapply(data.frame(t(m)),st_point))
Geometry set for 5 features
geometry type: POINT
dimension: XY
bbox: xmin: 0.2875775 ymin: 0.0455565 xmax: 0.9404673 ymax: 0.892419
epsg (SRID): NA
proj4string: NA
POINT(0.287577520124614 0.0455564993899316)
POINT(0.788305135443807 0.528105488047004)
POINT(0.4089769218117 0.892419044394046)
POINT(0.883017404004931 0.551435014465824)
POINT(0.940467284293845 0.456614735303447)
另一种选择是通过
sp
对象进行操作,但我不想这样做。我也希望只使用基本 R 解决方案,不转换为 data.table 或 tbl 等数据格式。我是否错过了简单的
as(m,"sf")
函数或类似函数?
st_as_sf()
怎么样? - RoyalTS