h-a-graham

 library(raster)
#> Loading required package: sp
#> Warning: no function found corresponding to methods exports from 'raster' for:
#> 'area'
 
 circ_matrix <- function(n){
   if (!n %% 2){
     stop("n must be an odd number.")
   } 

h-a-graham

library(tibble)
library(purrr)

t1 <- tibble(`C1 Template` = sample(c(166:180), 10),
             `C2 Template` = sample(c(120:130), 10))

t1
#> # A tibble: 10 × 2
#>    `C1 Template` `C2 Template`

h-a-graham

#' download drive files from dribble with path column
#'
#' @param x a dribble with path column
#' @param .overwrite logical. Should files be overwritten
#' 
#' @return The original input dribble
#' @noRd
drive_down_files <- function(x, .overwrite = TRUE) {
  

h-a-graham

#https://twitter.com/AnalyticalEdge/status/1582206425585324034

library(ggplot2)
library(stars)
library(sf)
library(dplyr)
library(plotly)

.d <- dim(volcano)

h-a-graham

#' Buffer a point along a terrain surface.
#'
#' @param dtm A SpatRaster Digital Terrain Model
#' @param pnts A SpatVector of points
#' @param .dist Numeric - the buffer distance.
#' @param .res Numeric - approximate resampling resolution of the (internally
#' calculated) cost raster.
#' @param .smooth Logical default FALSE. Should the output buffer be smoothed using `smoothr::smooth`?
#' @param .method The smoothing method to use - see `smoothr::smooth` for details.
#' @param ... Passed to `smoothr::smooth` for eg. to control smoothness.

h-a-graham

#' @title Get administritive outlines for a country
#' @description using the geoBoundaires API, get the administritive polygon(s)
#' for a country
#' @param country character vector: a country name
#' @param admin_level character vector: the admin level to download
#' @param quiet logical, should st_read be quiet?
#' @return sf object of the outlines
#' @details check out the documentation for the geoboundaries API at:
#' geoBoundaries.org
#'

h-a-graham

# Description: This script is used to download HLS data from NASA's Earthdata
# STAC API and build a cloud-free composite image.

# ----- Functions -----

# band mapping for HLS SL data
#' @return a named character vector of band mappings for HLS SL data
hlssl_band_mapping <- function() {
  c(
    B01 = "A", B02 = "B", B03 = "G", B04 = "R",

h-a-graham

#' Project an sf/sfc object to a generic projected coordinate system
#' @param x an sf or sfc object
#' @param proj a character vector. The projection to use. One of "laea", "aeqd",
#' "utm", "pconic", or "eqdc".
#' @param ellps a character vector. The ellipsoid to use. Select from
#' `sf_proj_info(type = "ellps")`.
#' @param opts a character vector. Additional proj options to pass to the
#' proj string. see details for more information.
#' @return an sf or sfc object
#' @details For further info about the available "generic" projects see:

h-a-graham

as_the_crow_flies <- function(city1, city2,
                              country1 = NULL, country2 = NULL,
                              units = "km", quiet = FALSE) {
  gcl <- tidygeocoder::geo(
    city = c(city1, city2), country = c(country1, country2),
    method = "osm", quiet = TRUE, progress_bar = FALSE, verbose = FALSE
  )
  l1 <- gcl[1, ]
  l2 <- gcl[2, ]

h-a-graham

library(terra)
#> terra 1.8.0

# create example landcover raster
r <- rast(
  extent = c(0, 15000, 0, 15000),
  resolution = 30,
  crs = "+proj=laea"
)