diff --git a/R/stack_rasters.R b/R/stack_rasters.R
index f97a68e..a9a3c66 100644
--- a/R/stack_rasters.R
+++ b/R/stack_rasters.R
@@ -29,7 +29,7 @@
 #' If `resolution` or `extent` are provided, they override the values from the
 #' reference raster.
 #' @param resampling_method The method to use when resampling to different
-#' resolutions.
+#' resolutions in the VRT.
 #' @param band_names Either a character vector of band names, or a function that
 #' when given a character vector of band names, returns a character vector of
 #' the same length containing new band names.
@@ -61,7 +61,6 @@ stack_rasters <- function(rasters,
                           resampling_method = "bilinear",
                           band_names,
                           gdalwarp_options = c(
-                            "-r", "bilinear",
                             "-multi",
                             "-overwrite",
                             "-co", "COMPRESS=DEFLATE",
diff --git a/man/rsi-package.Rd b/man/rsi-package.Rd
index d71e288..41b9227 100644
--- a/man/rsi-package.Rd
+++ b/man/rsi-package.Rd
@@ -4,11 +4,11 @@
 \name{rsi-package}
 \alias{rsi}
 \alias{rsi-package}
-\title{rsi: Utilities for Retrieving Satellite Imagery, Calculating Spectral Indices, and Wrangling the Outputs}
+\title{rsi: Efficiently Retrieve and Process Satellite Imagery}
 \description{
 \if{html}{\figure{logo.png}{options: style='float: right' alt='logo' width='120'}}
 
-Tools for downloading spatial data from 'spatiotemporal' asset catalogs ('STAC') and computing standard spectral indices from raster data. Methods focus on interoperability with the broader spatial ecosystem; function arguments and outputs use classes from 'sf' and 'terra', and data downloading functions support complex 'CQL2' queries using 'rstac'.
+Downloads spatial data from spatiotemporal asset catalogs ('STAC'), computes standard spectral indices from the Awesome Spectral Indices project (Montero et al. (2023) \doi{10.1038/s41597-023-02096-0}) against raster data, and glues the outputs together into predictor bricks. Methods focus on interoperability with the broader spatial ecosystem; function arguments and outputs use classes from 'sf' and 'terra', and data downloading functions support complex 'CQL2' queries using 'rstac'.
 }
 \seealso{
 Useful links:
diff --git a/man/stack_rasters.Rd b/man/stack_rasters.Rd
index b2d45d7..236560b 100644
--- a/man/stack_rasters.Rd
+++ b/man/stack_rasters.Rd
@@ -13,8 +13,8 @@ stack_rasters(
   reference_raster = 1,
   resampling_method = "bilinear",
   band_names,
-  gdalwarp_options = c("-r", "bilinear", "-multi", "-overwrite", "-co",
-    "COMPRESS=DEFLATE", "-co", "PREDICTOR=2", "-co", "NUM_THREADS=ALL_CPUS"),
+  gdalwarp_options = c("-multi", "-overwrite", "-co", "COMPRESS=DEFLATE", "-co",
+    "PREDICTOR=2", "-co", "NUM_THREADS=ALL_CPUS"),
   gdal_config_options = c(VSI_CACHE = "TRUE", GDAL_CACHEMAX = "30\%", VSI_CACHE_SIZE =
     "10000000", GDAL_HTTP_MULTIPLEX = "YES", GDAL_INGESTED_BYTES_AT_OPEN = "32000",
     GDAL_DISABLE_READDIR_ON_OPEN = "EMPTY_DIR", GDAL_HTTP_VERSION = "2",
@@ -49,7 +49,7 @@ If \code{resolution} or \code{extent} are provided, they override the values fro
 reference raster.}
 
 \item{resampling_method}{The method to use when resampling to different
-resolutions.}
+resolutions in the VRT.}
 
 \item{band_names}{Either a character vector of band names, or a function that
 when given a character vector of band names, returns a character vector of