define_bounds: Define boundaries with slope and intercept

In gopalpenny/anem: Analytical element models for simple aquifers

Description

Prepare boundaries with slope (m) and intercept (b) of all line objects, and add boundary ID (bID)

Usage

define_bounds(bounds_df, get_rectangular = TRUE)

Arguments

bounds_df	An sf object of straight lines (single segments), a data.frame with column names x1, y1, x2, y2, representing the endpoints, or a data.frame with columns m (slope) and b (intercept)
get_rectangular	logical that, if TRUE, forces bounds to be rectangular

Value

A data.frame containing slope (m) and intercept (b) for each line, along with original columns. If bounds is an sf object, the columns x1, y1, x2, y2 are automatically extracted using prep_wells_sf and added to the sf object before getting the slope and intercept. The function requires that the wells are labeled with a column of identifiers, pID. If it is not present, the function generates them.

Examples

bounds_df <- data.frame(bound_type=c("CH","NF","NF","NF"),
  x1=c(0,10,13,1),y1=c(0,10,9,-1),
  x2=c(10,13,1,0),y2=c(10,9,-1,0))
bounds <- define_bounds(bounds_df)


library(ggplot2)
ggplot(bounds) + geom_segment(aes(x1,y1,xend=x2,yend=y2)) + coord_equal()

bounds_df <- data.frame(bound_type=c("CH","CH","NF","NF"),
  x1=c(0,0,10,10),y1=c(0,10,10,0),
  x2=c(0,10,10,0),y2=c(10,10,0,0))
bounds <- define_bounds(bounds_df)
ggplot(bounds) + geom_segment(aes(x1,y1,xend=x2,yend=y2)) + coord_equal()

gopalpenny/anem documentation built on Dec. 20, 2020, 5:27 a.m.