Calculer l'aire de la section transversale pour faire varier la hauteur

Question

J'essaie de comprendre comment calculer l'aire d'une section transversale d'une rivière.

Pour la section transversale, j'ai la profondeur à tous les 25 cm sur la rivière de 5 m de large.

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) 

Si quelqu'un a des suggestions sur la façon dont cela pourrait être fait, il est très apprécié.

Answer 1

Nous pouvons utiliser le package sf pour créer un polygone présentant la section, puis calculer la surface. Notez que pour créer un polygone, il est nécessaire de fournir trois autres points comme c(0, 0), c(500, 0) et c(0, 0) lors de la création de la matrice m.

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