hunzikp · mwip · Feb 12, 2019 · Jun 28, 2019 · Jul 1, 2019
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -2,9 +2,10 @@ Package: velox
 Type: Package
 Title: Fast Raster Manipulation and Extraction
 Version: 0.2.0.9002
-Date: 2018-02-13
+Date: 2019-02-12
 Author: Philipp Hunziker
 Maintainer: Philipp Hunziker <hunzikp@gmail.com>
+Contributor: Matthias Weigand
 BugReports: https://github.com/hunzikp/velox/issues
 Description: C++ accelerated raster manipulation and extraction.
 Depends: R (>= 3.0.2)
@@ -19,8 +20,8 @@ Imports:
     sf (>= 0.5-5)
 LinkingTo: Rcpp, BH (>= 1.65.0-1)
 NeedsCompilation: yes
-SystemRequirements:	C++11
+SystemRequirements: C++11
 Packaged: 2017-11-30 12:15 UTC; hunzikp
 LazyData: TRUE
-RoxygenNote: 5.0.1
+RoxygenNote: 6.1.1
 Suggests: testthat
diff --git a/R/RcppExports.R b/R/RcppExports.R
@@ -21,6 +21,10 @@ meanfocal_cpp <- function(rasterband, weights, wrow, wcol, band) {
     .Call('_velox_meanfocal_cpp', PACKAGE = 'velox', rasterband, weights, wrow, wcol, band)
 }
 
+stddevfocal_cpp <- function(rasterband, wrow, wcol, band) {
+    .Call('_velox_stddevfocal_cpp', PACKAGE = 'velox', rasterband, wrow, wcol, band)
+}
+
 hittest_cpp <- function(rasterbands, dim, extent, res, polyX, polyY, polyCorners) {
     .Call('_velox_hittest_cpp', PACKAGE = 'velox', rasterbands, dim, extent, res, polyX, polyY, polyCorners)
 }

diff --git a/R/velox_focal.R b/R/velox_focal.R
@@ -39,6 +39,47 @@ VeloxRaster$methods(medianFocal = function(wrow, wcol, bands=1) {
 
 
 
+#' @title Standard Deviation focal
+#'
+#' @name VeloxRaster_stddevFocal
+#'
+#' @description
+#' Applies a standard deviation filter of dimension \code{wcol x wrow} to a VeloxRaster.
+#'
+#' @details
+#' Padding is currently not implemented.
+#'
+#' @param wrow y dimension of filter. Must be uneven integer.
+#' @param wcol x dimension of filter. Must be uneven integer.
+#' @param bands Numeric vector indicating bands where filter is applied.
+#'
+#' @return Void.
+#'
+#' @examples
+#' ## Make VeloxRaster with two bands
+#' mat1 <- matrix(1:100, 10, 10)
+#' mat2 <- matrix(100:1, 10, 10)
+#' vx <- velox(list(mat1, mat2), extent=c(0,1,0,1), res=c(0.1,0.1),
+#'             crs="+proj=longlat +datum=WGS84 +no_defs")
+#' ## Median focal
+#' vx$stddevFocal(wrow=5, wcol=5, bands=c(1,2))
+#'
+NULL
+VeloxRaster$methods(stddevFocal = function(wrow, wcol, bands=1) {
+  "See \\code{\\link{VeloxRaster_stddevFocal}}."
+  if (any(!(bands %in% 1:nbands))) {
+    stop(paste("VeloxRaster only has", nbands, "bands."))
+  }
+  if (wrow < 0 | wcol < 0 | (wrow %% 2) == 0 | (wcol %% 2) == 0) {
+    stop(paste("wrow and wcol must be positive uneven integers."))
+  }
+  for (i in bands) {
+    rasterbands[[i]] <<- stddevfocal_cpp(rasterband=rasterbands[[i]], wrow=wrow, wcol=wcol, band=i)
+  }
+})
+
+
+
 #' @title Sum focal
 #'
 #' @name VeloxRaster_sumFocal

diff --git a/man/BoostBoxGrid-class.Rd b/man/BoostBoxGrid-class.Rd
diff --git a/man/BoostFactory.Rd b/man/BoostFactory.Rd
diff --git a/man/BoostGeometries-class.Rd b/man/BoostGeometries-class.Rd
diff --git a/man/BoostGrid-class.Rd b/man/BoostGrid-class.Rd
diff --git a/man/BoostMultiLines-class.Rd b/man/BoostMultiLines-class.Rd
diff --git a/man/BoostMultiPoints-class.Rd b/man/BoostMultiPoints-class.Rd
diff --git a/man/BoostMultiPolygons-class.Rd b/man/BoostMultiPolygons-class.Rd
diff --git a/man/BoostPointGrid-class.Rd b/man/BoostPointGrid-class.Rd
diff --git a/man/BoxGrid.Rd b/man/BoxGrid.Rd
diff --git a/man/MultiLineCollection.Rd b/man/MultiLineCollection.Rd
diff --git a/man/MultiPointCollection.Rd b/man/MultiPointCollection.Rd
diff --git a/man/MultiPolygonCollection.Rd b/man/MultiPolygonCollection.Rd
diff --git a/man/PointGrid.Rd b/man/PointGrid.Rd
diff --git a/man/VeloxRaster-class.Rd b/man/VeloxRaster-class.Rd
diff --git a/man/VeloxRaster_aggregate.Rd b/man/VeloxRaster_aggregate.Rd
diff --git a/man/VeloxRaster_as.RasterBrick.Rd b/man/VeloxRaster_as.RasterBrick.Rd
diff --git a/man/VeloxRaster_as.RasterLayer.Rd b/man/VeloxRaster_as.RasterLayer.Rd
diff --git a/man/VeloxRaster_as.RasterStack.Rd b/man/VeloxRaster_as.RasterStack.Rd
diff --git a/man/VeloxRaster_as.matrix.Rd b/man/VeloxRaster_as.matrix.Rd
diff --git a/man/VeloxRaster_col2im.Rd b/man/VeloxRaster_col2im.Rd
diff --git a/man/VeloxRaster_crop.Rd b/man/VeloxRaster_crop.Rd
diff --git a/man/VeloxRaster_drop.Rd b/man/VeloxRaster_drop.Rd
diff --git a/man/VeloxRaster_extract.Rd b/man/VeloxRaster_extract.Rd
diff --git a/man/VeloxRaster_extract_points.Rd b/man/VeloxRaster_extract_points.Rd
diff --git a/man/VeloxRaster_getCoordinates.Rd b/man/VeloxRaster_getCoordinates.Rd
diff --git a/man/VeloxRaster_get_data_type.Rd b/man/VeloxRaster_get_data_type.Rd
diff --git a/man/VeloxRaster_im2col.Rd b/man/VeloxRaster_im2col.Rd
diff --git a/man/VeloxRaster_meanFocal.Rd b/man/VeloxRaster_meanFocal.Rd
diff --git a/man/VeloxRaster_medianFocal.Rd b/man/VeloxRaster_medianFocal.Rd
diff --git a/man/VeloxRaster_rasterize.Rd b/man/VeloxRaster_rasterize.Rd
diff --git a/man/VeloxRaster_stddevFocal.Rd b/man/VeloxRaster_stddevFocal.Rd
diff --git a/man/VeloxRaster_sumFocal.Rd b/man/VeloxRaster_sumFocal.Rd
diff --git a/man/VeloxRaster_write.Rd b/man/VeloxRaster_write.Rd
diff --git a/man/bg_intersects.generic.Rd b/man/bg_intersects.generic.Rd
diff --git a/man/boost.Rd b/man/boost.Rd
diff --git a/man/length.BoostGeometries.Rd b/man/length.BoostGeometries.Rd
diff --git a/man/plot.BoostGeometries.Rd b/man/plot.BoostGeometries.Rd
diff --git a/man/subset.BoostGeometries.Rd b/man/subset.BoostGeometries.Rd