shimming-toolbox / b0-fieldmap-realistic-simulation

Optimization of shim coil layout placement using simulated B0 and B1 maps constructed from segmented MRI scans of the head and full spine region. This is a collaboration between NeuroPoly Lab (Montreal) and MGH Martinos Center (Boston).
MIT License
0 stars 0 forks source link

Verify and brush in voxels on sides (around and bottom) of torso that weren't labelled as tissue by TotalSegmentator #16

Open mathieuboudreau opened 1 week ago

mathieuboudreau commented 1 week ago

eg

Screenshot 2024-11-04 at 10 19 58 AM

mathieuboudreau commented 1 week ago

If these aren't filled in, they will be assigned as air and will create pockets of air when padding via "edge" mode (copies the values of the boundaries for the padding)

mathieuboudreau commented 1 week ago

Could also be due to this kind of "roundness" near the bottom of the torso / bumpy; maybe manually smoothing and stretching it to the boundary would improve the red-blue gradient in some maps.

Screenshot 2024-11-04 at 10 29 58 AM Screenshot 2024-11-04 at 10 29 44 AM

mathieuboudreau commented 1 week ago

Hmm, a "good" subject also has this and wasn't as much of an issue far from there,

Screenshot 2024-11-04 at 10 33 01 AM

mathieuboudreau commented 1 week ago

Using the 3D view I think I've gotten a better insight.

All my saggital slices are cut through the trachea; thus, between the two lungs.

I checked for 4 subjects, 2 "good" (i.e. little gradient in torso) and 2 "bad" (i.e. stronger gradient in the torso).

From this small cohort, it seems that the "good ones had 1) lung segmentation that were more symmetrical in terms of "fullness", whereas the "bad" had the left (from the patient's POV) one much thinner than the right. Maybe this is what caused the gradient in B0?

"Good" sub 1

Screenshot 2024-11-04 at 11 02 34 AM Screenshot 2024-11-04 at 10 51 30 AM Screenshot 2024-11-04 at 10 51 26 AM Screenshot 2024-11-04 at 10 51 22 AM

"Bad" sub 2

Screenshot 2024-11-04 at 11 03 22 AM Screenshot 2024-11-04 at 10 51 22 AM Screenshot 2024-11-04 at 10 50 41 AM Screenshot 2024-11-04 at 11 05 18 AM

"Good" sub 3

Screenshot 2024-11-04 at 11 06 11 AM

Screenshot 2024-11-04 at 10 58 15 AM Screenshot 2024-11-04 at 10 58 13 AM Screenshot 2024-11-04 at 10 58 09 AM

"Bad" sub 4

Screenshot 2024-11-04 at 11 07 24 AM Screenshot 2024-11-04 at 10 57 30 AM Screenshot 2024-11-04 at 10 57 28 AM

Screenshot 2024-11-04 at 10 57 32 AM Screenshot 2024-11-04 at 11 08 34 AM

mathieuboudreau commented 1 week ago

And then, because we're demodulating for the tissue B0 values, the subjects that have "stronger b0 gradients" will have more (higher) red values, driving the mean value up, and thus making the neck/head of the subjects darker blue, eg

Screenshot 2024-11-04 at 9 58 36 AM

instead of red-ish from the subjects that don't have this lung-induced gradient in the torso.

@jcohenadad maybe I should write two sentences on this in the abstract results section? It's not that the "head/neck" have more deviation towards negative B0 (blue) values for some subjects, but that the mean B0 changed due to higher lung asymmetry which then induced more B0 variation between/in the torso. Does this make sense? Pinging @evaalonsoortiz too

mathieuboudreau commented 1 week ago

Here's a comparison of the B0 map for an axial slice for a "good" (sub 3 above) and a "bad" subject (sub 4 above)

"Good" sub 3

Screenshot 2024-11-04 at 11 18 44 AM

"Bad" sub 4

Screenshot 2024-11-04 at 11 17 46 AM

They both seem to "touch" the walls of the FOV as much. Not sure if I can see another explanation for this apparent B0 gradient between the lungs, @evaalonsoortiz anything you see here?

mathieuboudreau commented 1 week ago

So if I go even lower, the furthest I can go before the "good" subject's left lung starts to disapear is this:

Screenshot 2024-11-04 at 11 28 27 AM

Whereas the "bad" subject, it's here:

Screenshot 2024-11-04 at 11 28 08 AM

So maybe it's more about the symmetry in length than thickness/thinness; for the "bad ones" there are a lot of slices where there's left lung but no/little right lung, whereas the "good" subject's are pretty symmetric

evaalonsoortiz commented 1 week ago

I agree with your explanations (both wrt to the left/right symmetry issue along the length of the spinal cord, and within the axial plane). Essentially, I think what this is all suggesting is that we should explore the possibility of acquiring anatomical scans with FOVs that encompass the whole torso (capturing the entirety of the lungs and their unique shapes seems important). If that is not doable, then perhaps in the future we should search for a large CT-based dataset that we can use for B0 simulations.

mathieuboudreau commented 1 week ago

The other problem is that for all images, the lower torso has really, really low SNR. So, the lower lung segmentations performed more or less ok there, decent but not phenomenal,

Screenshot 2024-11-04 at 12 13 49 PM Screenshot 2024-11-04 at 12 13 42 PM