Search
Question: Focal stacking of images in R/EBImage
0
gravatar for phil.chapman
13 months ago by
phil.chapman120
United Kingdom
phil.chapman120 wrote:

Has anyone got experience of carrying out a focal stacking operation on an image stack in R?  This is where one has a series of images with the same field of view but different focal planes, and these are combined into a single image where everything is magically in focus. 

Thanks, Phil

ADD COMMENTlink modified 13 months ago by Bernd Fischer540 • written 13 months ago by phil.chapman120
2
gravatar for Bernd Fischer
13 months ago by
Bernd Fischer540
Germany / Heidelberg / DKFZ
Bernd Fischer540 wrote:

The easiest way to combine images from a confocal microscope is a max-projection, which can be obtained by

Img = readImage(c("img1.tif","img2.tif","img3.tif","img4.tif"))

MaxProj = apply(Img, 1:2, max)

display(MaxProj)

In many cases this is sufficient. However, other algorithms are not available in R.

 

ADD COMMENTlink written 13 months ago by Bernd Fischer540

Thanks Bernd!  So you're just taking the maximum intensity for each pixel?

ADD REPLYlink written 13 months ago by phil.chapman120

Hi Phil, this method is now provided by the new Bioconductor package MaxContrastProjection. The package also implements the novel maximum contrast projection algorithm which avoids some of the problems of the former method and yields cleaner results.

ADD REPLYlink modified 6 months ago • written 6 months ago by Andrzej Oleś670
1
gravatar for Andrzej Oleś
13 months ago by
Andrzej Oleś670
EMBL Heidelberg, Germany
Andrzej Oleś670 wrote:

Hi Phil,

focus stacking is now provided by the Bioconductor package MaxContrastProjection. It implements various projections methods of 3D image stack along the z-dimension into 2D, including the maximum intensity projection, and a novel maximum contrast projection algorithm.

An alternative would be to try some third-party tools, see e.g. this list.

ADD COMMENTlink modified 6 months ago • written 13 months ago by Andrzej Oleś670
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 158 users visited in the last hour