I’m trying to analyze a series of mosaic images which are very big in size, so I’m attaching a sample area below (this is a combination of 4 layers (blue for DAPI, green cells, red cells, and white marker area)
where I have two areas I’m looking at: one positive and the other negative for the white marker stain (basically the rest of the epithelium that’s not white). In each of these areas, I’m looking to quantify the number of cells double positive for the red and green stain out of the total number of green positive cells.
(For my first pipeline, I did roughly the following:
Identify primary objects by nuclei>threshold green cells>mask the nuclei with green cells>threshold red cells>mask the nuclei/green cells layer with the red cell layer >threshold the white area>mask the nuclei/green/red cells with the white area, count the objects.) However, one of the issues I ran into was that the red marker stain is not uniform throughout the tissue, i.e. it is fainter in the center vs the edges (picture attached)
thus thresholding the entire image for the red positive cells does not work. Since the strength of the red stain roughly correlates with the white area, I decided to threshold the red areas separately, based on the white area masks…).
So, one of my questions is if there’s any other way I can threshold the red cells? I think another option I’ll try is to
a) adjust the contrast in ImageJ to enhance the contrast or b) since this is a mosaic picture for which I have all of the original unstitched images, would it be possible to threshold the red channel for each of the pictures separately (I have roughly 150 images that make up the mosaic image)?
I’m attaching 2 batches of images of the new pipeline I made (the images used for analysis are mosaics of very big size, so I’m using a representative part as an example), where I’m listing the other questions I have. I’m copying them here as well:
- The white area is not uniformly stained (step 2 in the pipeline), i.e. not all of the cells in that area are stained positive for the marker. Is there a way to fill in the gaps for that layer in include the entire area based on its border? Alternatively, is there a way to hand select this area?
- I’m trying to count the cells in the epithelium area only, but unfortunately it turns out some nuclei co-localize with green and red in the surrounding areas (step 11 in the pipeline). Is there a way to exclude the area from the counting ? Perhaps manually?
Thanks so much for any imput, and please let me know if you have any other clarifications!