spinal-cord-7t / coil-qc-code

7T Spinal Cord Coil QC Analysis Code
0 stars 0 forks source link

Plot/figure brainstorming for cross-site comparisons #7

Open evaalonsoortiz opened 7 months ago

evaalonsoortiz commented 7 months ago

Our dataset incudes 7 sites with 3 subjects per-site (the same 3 subjects were, or are to be, scanned across all 7 sites). The goal is to compare coil designs. The purpose of this issue is to develop a consensus on what metrics we want to derive from the multi-subject data and how we want to visualize those metrics across sites.

Please comment on whether or not you agree with the "goals" and the proposed figures. Alternatively, please suggest and "goal" and/or figure for a given scan type.

Images:


TFL_B1_C3C4 ---> B1+ map acquired with the Optimal Reference Voltage

Goals:

Suggested figures: Figure 1:

Screenshot 2023-12-04 at 9 30 09 PM

Figure 1 shows us how uniform B1+ is along the SC but it does not tell us how efficient the coil is. For that we may want to include the following figure:

Figure 2:

Screen Shot 2023-12-18 at 3 54 53 PM

COILQA_SAG_LARGE (H>F phase encoding)

Images: SNR, g-factor, correlation matrix

Goals:

Suggested figures:

For SNR: same as those for TFL (B1+)

For g-factor:

Using an enlarged spinal cord (SC) ROI, compute the average g-factor along the SC report it along with the standard deviation.

Figure 3:

Screenshot 2023-12-18 at 4 58 36 PM

The drawback of Figure 3 is that the little bars for each subject might overlap quite a bit making it overall hard to read. Alternatively we can represent the g-factor along the SC for each subject as a histogram/distribution in the form of a ridge-plot:

Figure 4:

Screenshot 2023-12-18 at 5 00 40 PM

Figures 3 or 4 can be a sub-plot for a larger figure that includes other acceleration factors (3x3, 3x2).

For correlation matrices:

Compute the average correlation matrix across subjects for each site and show the average matrix for each of the 7 sites in a tiled arrangement.


COILQA_SAG_SMALL (H>F phase encoding)

Goals:

Suggested figures:

Same as Figure 3 or 4. This could be included in supporting documents.


COILQA_TRA (L>R phase encoding)

Goals:

Suggested figures:

For SNR: Compute in-plane standard-deviation and average that for all slices. Show this value for each subject, and across sites, like in Figure 1.

For g-factor: Same as Figure 3 or 4.


DREAM_LARGE

Goals:

Suggested figures:

Same as Figure 2. Note: A bias in measured B1+ where the refV is not well calibrated can also lead to differences in measured B1+ uniformity. Therefore we may want to have another version of Figure 1.


DREAM_MEDIUM + DREAM_MEDIUM_066 + DREAM_MEDIUM_HWLIMIT

Goals:

Suggested figures:

Compute the average B1+ along the spinal cord and average across subjects, for each refV. The former will be included in a tiled figure where each sub-plot represents one site.

Figure 5:

Screenshot 2023-12-18 at 6 33 04 PM

Each site has 3 violin plots (refV, 0.66refV, 1.5refV) containing the distribution of B1+ values for all subjects.

Figure 6:

Screenshot 2023-12-18 at 6 34 31 PM

GRE ---> uncombined GRE images

Goals:

Suggested figures:

A tiled figure (containing each channel's GRE magnitude) could be generated for each subject at each site and presented in supporting documentation.

jcohenadad commented 7 months ago

Thank you for the detailed figure layouts-- it makes it easier to comment.

evaalonsoortiz commented 6 months ago
jcohenadad commented 6 months ago

This would be a lot of visual data though (7 sites times R=2x2, 2x3, 3x2, 3x3). How about picking a one acceleration only and showing it across sites for one subject?

Yes, this is also what I had in mind. Showing 2x2 seems reasonable.

GRE / Shall we not show this data at all then?

I would absolutely show at least one anat image. This is what most (ie not coil experts) people will want to see. If we don't have mGRE, GRE will suffice.

evaalonsoortiz commented 6 months ago

GRE / Shall we not show this data at all then?

I would absolutely show at least one anat image. This is what most (ie not coil experts) people will want to see. If we don't have mGRE, GRE will suffice.

It would still be a very big figure: all-coil images for one subject times 7 sites. I think supplementary figs is a good place for that, do you agree/disagree?

jcohenadad commented 6 months ago

if we pick one subj, one slice, it’s only 7 images, unless i am missing something?

evaalonsoortiz commented 6 months ago

We have 1 image per coil element (not coil combined). We can coil combine but we had agreed at one point to keep them uncombined in order to assess the coil elements themselves.

evaalonsoortiz commented 5 months ago

Summary of the agreed-upon figures:

sub-xxx_acq-famp_TB1TFL

Manuscript figure: 3 panel figure with (A) single subject representative B1+ map, (B) median B1+ within spinal cord ROI vs. spinal cord level for that subject, (C) median within slice (and ROI) and then inter-quartile range (IQR) across slices for each subject (with the median across subjects) and each site

Notebook figures:

Note1: Note that for Figure 1, @jcohenadad wanted to explore how the data looks here, so this figure may evolve. We can play around with it once we have the code. _Note2: TB1TFL maps are output as a flip angle map in units of flip angle * 10 (in degrees), to convert to nT/V use the same code as in https://github.com/shimming-toolbox/rf-shimming-7t/blob/main/data_processing.ipynb_

sub-xxx_acq-refv_TB1DREAM

Manuscript figure: same as for famp_TB1TFL

Notebook figures:

Note1: DREAM outputs flip angle maps in degrees and B1+ maps that are normalized and in percent units Note2: the 3 reference voltage acquisitions have to be combined into one in the following way:

Note3: The FOV of the DREAM scans are smaller than the FOV of the TFL scans, therefore a one-to-one comparison between Figures 1 and 3, or Figures 2 and 4 can not be easily made

sub-xxx_acq-coilQaSagLarge_SNR

Manuscript figure description: 3 panel figure with (A) single subject representative SNR map, (B) median SNR within spinal cord ROI vs. spinal cord level for that subject, (C) median SNR within spinal cord ROI vs. spinal cord level averaged across subjects per site @jcohenadad I would argue here that B is not necessary and just keep A and C, let me know what you think.

Notebook figures:

jcohenadad commented 5 months ago

some comments:

jcohenadad commented 5 months ago

@jcohenadad I would argue here that B is not necessary and just keep A and C, let me know what you think.

Agreed. Single-subject statistics could go over the map, as done in the RF shimming paper, eg:

image