R/plot_bpip.R

In dKvale/bpip: Interface to EPA's downwash model BPIP

#' Plot buildings from a BPIP dataframe
#'
#' Plot a BPIP building dataframe.
#' @param data Data frame of bpip parameters. Run 'bld_tbl()' for an example.
#' @keywords plot bpip buildings downwash
#' @export
#' @examples
#' buildings <- new_bpip()
#' 
#' plot_bpip(data = buildings)
# 
#

plot_bpip <- function(data) {
  
  source_xy <- data[1, ]$source_xy[[1]]

  coord_msg <- TRUE
  
  for(i in 1:nrow(data)) {
    
    if(is.null(data[i, "bld_xcoords"]) | 
       is.null(data[i, "bld_ycoords"]) | 
       length(unlist(data[i, "bld_xcoords"], ",")) < 3 |
       length(unlist(data[i, "bld_ycoords"], ",")) < 3 ) {
      
      if(coord_msg) {
        print("Building vertices were calculated based on a rectanglular footprint. To create a custom building shape provide 3 or more x,y coordinates.")
        
        coord_msg <- FALSE
      }
      
      coords <- bld_coords(source_xy           = data[i, ]$source_xy[[1]],
                           dist_from_source    = data[i, ]$dist_from_source,
                           angle_from_source   = data[i, ]$angle_from_source,
                           width_x             = data[i, ]$width_x,
                           length_y            = data[i, ]$length_y,
                           bld_rotation        = data[i, ]$bld_rotation,
                           angle_units         = data[i, ]$angle_units,
                           show_plot           = FALSE)
      
    
      data[i, ]$bld_xcoords  <- list(coords$x_coords)
      data[i, ]$bld_ycoords  <- list(coords$y_coords)
    }
  }
  
  xcoords <- unlist(data$bld_xcoords)
  ycoords <- unlist(data$bld_ycoords)
  
  bld_range <- max(abs(source_xy[2] - ycoords), abs(source_xy[1] - xcoords))
  
  graphics::plot(xcoords,
                 ycoords,
                 col  = "steelblue",
                 xlab = paste("South"),
                 ylab = paste("West"),
                 xlim = c(source_xy[1] - bld_range, source_xy[1] + bld_range),
                 ylim = c(source_xy[2] - bld_range, source_xy[2] + bld_range))
  
  for(i in 1:nrow(data)) {
    graphics::polygon(unlist(data[i, ]$bld_xcoords), unlist(data[i,]$bld_ycoords), col = "steelblue") 
    }
  
  graphics::points(source_xy[1], source_xy[2], pch =13, col = "orange", cex=1.8)
}