# NOT RUN {
library(EBImage)
img <- ReadImageData(system.file('extdata', '50.tif', package = 'nandb'))
display(normalize(img[, , 1]), method = 'raster')
brightness <- Brightness(img, tau = NA, mst = "Huang")
KmerPlot(brightness, 1.02)
KmerPlot(brightness, 1.12)
KmerPlot(brightness, 100)
KmerPlot(brightness, 1.02, log.trans = TRUE)
KmerPlot(MedianFilterB(brightness), 1.02)
# }
Run the code above in your browser using DataLab