ポイントからラインストリングを効率的に作成するにはどうすればよいですか?
2つの別々のデータフレームにgeomPOINTがあります。私がやりたいのは、ポイントを線で接続することです(後でマップ上で)。そのため、これらのデータフレームからポイントのペアごとにLinestringを作成します。私はそれをこのようにしました:
coordsCust <- table %>%
st_as_sf(coords = c("lonCust","latCust"), crs = 4326)
coordsApp <- table %>%
st_as_sf(coords = c("lonApp","latApp"), crs = 4326) %>%
st_geometry()
およびラインストリング:
lines <- st_sfc(mapply(function(a,b){
st_cast(st_union(a,b),"LINESTRING")},
coordsCust$geometry, coordsApp$geometry, SIMPLIFY=FALSE))
このコードは機能します。ポイントのペアごとに行ごとにラインストリングを作成できます。
LINESTRING (14.035 51.65182, 14.33418 53.53346)
LINESTRING (20.42767 49.98073, 16.62978 52.31037)
LINESTRING (20.18762 50.03337, 16.62978 52.31037)
LINESTRING (19.04625 49.79234, 16.62978 52.31037)
LINESTRING (21.35808 50.92382, 16.62978 52.31037)
問題は、30 000行の場合、このソリューションの動作が非常に遅く、約21秒になることです。ポイントからラインストリングを作成する他の方法はありますか?はるかに高速に動作するものはありますか?私はウェブ上でいくつかの解決策を探しましたが、無駄でした。sfを行列に変換して使用pmap
する方法について読んだことがありますが、ここで実装する方法がわかりません。
更新: sfheaders :: sf_linestring関数を使用する場合は、両方のデータセットのジオメトリを結合する必要があります。私はこのようにします:
df <- cbind(coordsCust,coordsApp)
最終的なデータフレーム(私はその最も重要な部分を示しました)を以下に示します。
残念ながら、sf_linestringはこのデータフレームでは正しく機能しません。画面に表示されているように、各行のPOINT間に個別にラインストリングを作成する必要があります。
回答
2 SymbolixAU
例のデータセットがないと、質問に完全に答えることは困難です。しかし、data.frameを「長い」形式にsfheaders
することができれば、これを瞬時に行うことができます
n <- 30000
df <- data.frame(
x = rnorm(n)
, y = rnorm(n)
)
df$id <- rep(1:(n/2), each = 2)
sfheaders::sf_linestring(
obj = df
, x = "x"
, y = "y"
, linestring_id = "id"
)
# Simple feature collection with 15000 features and 1 field
# geometry type: LINESTRING
# dimension: XY
# bbox: xmin: -4.297631 ymin: -4.118291 xmax: 3.782847 ymax: 4.053399
# CRS: NA
# First 10 features:
# id geometry
# 1 1 LINESTRING (0.2780517 0.243...
# 2 2 LINESTRING (0.4261505 2.503...
# 3 3 LINESTRING (0.8662821 -0.11...
# 4 4 LINESTRING (-0.5335952 -0.1...
# 5 5 LINESTRING (1.154309 -1.352...
# 6 6 LINESTRING (0.05512324 -0.4...
# 7 7 LINESTRING (1.945868 -0.744...
# 8 8 LINESTRING (0.0427066 -0.08...
# 9 9 LINESTRING (0.06738045 0.41...
# 10 10 LINESTRING (0.4128964 -0.04...