## ----style, echo=FALSE, results='asis'------------------------------------- BiocStyle::latex() ## ----eval=FALSE------------------------------------------------------------ # if (!requireNamespace("BiocManager", quietly=TRUE)) # install.packages("BiocManager") # BiocManager::install("MaxContrastProjection") ## -------------------------------------------------------------------------- library(MaxContrastProjection) library(EBImage) data(cells) dim(cells) ## -------------------------------------------------------------------------- max_contrast = contrastProjection(imageStack = cells, w_x = 15, w_y = 15, smoothing = 5, brushShape = "box") ## -------------------------------------------------------------------------- max_contrast_large = contrastProjection(imageStack = cells[30:95, 55:115,], w_x = 15, w_y = 15, smoothing = 5, brushShape = "box") max_contrast_small = contrastProjection(imageStack = cells[30:95, 55:115,], w_x = 3, w_y = 3, smoothing = 5, brushShape = "box") ## ----echo=FALSE------------------------------------------------------------ display(max_contrast_large / max(cells), method = "raster") ## ----echo=FALSE------------------------------------------------------------ display(max_contrast_small / max(cells), method = "raster") ## -------------------------------------------------------------------------- contrastStack = getContrastStack(imageStack = cells, w_x = 15, w_y = 15, brushShape = "box") indexMap = getIndexMap(contrastStack = contrastStack, smoothing = 5) max_contrast_fromMap = projection_fromMap(imageStack = cells, indexMap = indexMap) ## -------------------------------------------------------------------------- max_intensity_proj = intensityProjection(imageStack = cells, projType = "max") min_intensity_proj = intensityProjection(imageStack = cells, projType = "min") mean_intensity_proj = intensityProjection(imageStack = cells, projType = "mean") median_intensity_proj = intensityProjection(imageStack = cells, projType = "median") sd_intensity_proj = intensityProjection(imageStack = cells, projType = "sd") sum_intensity_proj = intensityProjection(imageStack = cells, projType = "sum") ## ----sessioninfo, results='asis'------------------------------------------- toLatex(sessionInfo())