Hi Dilek,
Yes, "--save-posteriors Lesions" will save a probability map with path *yourOutputDir*/posteriors/Lesions.mgz
1 and 2) Lesions are labeled as 99 in seg.mgz. A value of 0.04 in the stats file indicate that the lesion volume is 0.04 mm^3, so it makes sense that seg.mgz has no voxels labeled as 99. However, it means that something went wrong with the segmentation. What was the command you used to segment these DIR images? I'd probably use --lesion-mask-pattern 1 for them (although I've never tried to run SAMSEG on DIR images). Do you have a T1w scan to add for better whole-brain segmentation and consequently better lesion segmentation?
3) You can use the lesion probability map obtained with the --save-posteriors Lesions flag and use mri_binarize and mri_segstats to obtain binary segmentation and corresponding volumes at different thresholds. Note that you don't need to use the same threshold of run_samseg.
Hope it helps,
Stefano
------------------------------
Message: 2
Date: Thu, 18 Aug 2022 19:45:14 +0300
From: Dilek ?zkan <dozkan521(a)gmail.com>
Subject: [Freesurfer] samseg lesions labels and values
To: freesurfer(a)nmr.mgh.harvard.edu
Message-ID:
<CAPQSyvKJBirH9oOdKXeg=k8BUfdohYvMMoaLnscy-5X0GOMZ7w(a)mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
External Email - Use Caution
Hello dear all,
I want to use Freesurfer-Samseg for MS lesion
segmentation(cross-sectional). I have read that I need to add
"--save-posteriors Lesions" to run_samseg to obtain lesion probabilities.
1)In the stats file I can see the value of "Lesions" (0.04) But when I try
to see on freeview, it says label=99 (lesion label) has no value. So, I
can't see the lesions.
2)There are a few labels related to lesions. For example wm lesions;
(498: wmsa, 499: other_wmsa) I have no volume related to these labels
neither in my stats file nor on freeview. Should I get only the "Lesions
label" for the MS lesion mask?
3)As input, I am using only DIR ?mages. I w?ill try different thresholds. I
have read that by playing mri_binarize and mri_segstats, we can try
different thresholds without re-running anything. So, the threshold value
required for mri_binarize command refers to the same threshold input in the
initial the command of *run_samseg*?
Best,
Wish