さまざまな高さの断面の面積を計算します

Question

私は川の断面の面積を計算する方法を理解しようとしています。

断面については、5m幅の川の上に25cmごとに深さがあります。

x_profile <- seq(0, 500, 25) 
y_profile = c(50, 73, 64, 59, 60, 64, 82, 78, 79, 76, 72, 68, 63, 65, 62, 61, 56, 50, 44, 39, 25) 

誰かがこれをrでどのように行うことができるかについてのご意見がありましたら、高く評価されます。

Answer 1

sfパッケージを使用して、断面を示すポリゴンを作成し、面積を計算することができます。ポリゴンを作成するには、マトリクスmを作成するときに、c(0, 0),c(500, 0)、c(0, 0)という3つのポイントを追加する必要があります。

x_profile <- seq(0, 500, 25) 
y_profile <- c(50, 73, 64, 59, 60, 64, 82, 78, 79, 76, 72, 
       68, 63, 65, 62, 61, 56, 50, 44, 39, 25) 

library(sf) 

# Create matrix with coordinates 
m <- matrix(c(0, x_profile, 500, 0, 0, -y_profile, 0, 0), 
      byrow = FALSE, ncol = 2) 

# Create a polygon 
poly <- st_polygon(list(m)) 

# View the polygon 
plot(poly) 

# Calcualte the area 
st_area(poly) 
31312.5