Closed lcolladotor closed 3 years ago
We will be using dotdotdot for segmenting signal as this is IF.
Ok! Thanks for the clarification =)
The background suppression step in dotdotdot
removes artifacts with low intensity in the MAP2 and GFAP channel (which actually has to be masked out from Abeta and pTau channels). So I adapted new segmentation code, made the functions segmentIF
and refineIF
. The refineIF allows the user to manually set thresholds on Abeta, GFAP, MAP2, and pTau channels if the user feels the automatic thresholds generated from segmentIF do not work well. The final outputs of these functions are for every capture area, described below.
VIFAD1_V10A27-004_A1_segmentation.mat
= final segmentation of (capture area A1, experiment/slide 1) all channels stored in mat format(reuired input for countSpots and spotspotcheck functions).VIFAD1_V10A27-004_A1_test.png
= a test field picked by @shkwon17 that represents the entire capture area, displays the segmentation results. The image has 4 rows of tiles (raw, segmented, masked, size thresholded) and the columns correspond to different channels (Abeta, DAPI, GFAP, Lipofuscin, MAP2, pTau).VIFAD1_V10A27-004_A1_Abeta.png
= segmented binary image of the entire Abeta channelVIFAD1_V10A27-004_A1_pTau.png
= segmented binary image of the entire pTau channel
Once the images from #3 are ready and we have run spaceranger in #4, Maddy can use
VistoSeg
to compute the number of cells per spot as shown at http://research.libd.org/VistoSeg/step-4-gui-to-count-nuclei-in-a-visium-spot.html.