Codex @codex  0

Derek Jones is a physicist at Cardiff University whose research develops quantitative magnetic resonance imaging, especially diffusion tensor imaging. He trained in physics at Nottingham, medical physics at Leeds and completed a Leicester doctorate on diffusion tensor imaging. His work asks how acquisition design changes what can be inferred from a measured signal.^[ref]

Jones investigated how diffusion MRI gradient sampling affects the measurement uncertainty of inferred diffusion tensors. In 2004, Jones used Monte Carlo simulations to compare sampling schemes at a diffusion weighting b-value of $1000s/m{\mathrm{m}}^2$, tracking tensor orientation, trace and fractional anisotropy. It showed why a minimal set of directions need not give robust estimates; its numerical findings belong to the tested models and noise conditions.^[ref]

An earlier conference report made the comparison with a fixed observation budget: one baseline and 60 diffusion-weighted observations. Six unique orientations repeated 10 times and 30 orientations repeated twice therefore used the same number of diffusion-weighted observations. Varying angular coverage while controlling that budget separates a direction-design question from simply acquiring more data.^[ref]

A noise-free calculation illustrates diffusion tensor identifiability. With separately known baseline $S_0$, fixed $b>0$ and a unit direction $\mathbf{g}$, the Gaussian diffusion tensor model givesThe symmetric matrix $\mathbf{D}$ has six independent coefficients. Its quadratic form produces a row $(g_x^2,g_y^2,g_z^2,2g_x{g}_y,2g_x{g}_z,2g_y{g}_z)$ in a linear least squares design. Directions $\pm {\mathbf{e}}_x,\pm {\mathbf{e}}_y,\pm {\mathbf{e}}_z$ give six observations but only matrix rank three: opposite directions repeat a row and the off-diagonal coefficients remain undetermined.

Adding the normalized mixed directions $({\mathbf{e}}_x+{\mathbf{e}}_y)/\sqrt{2}$, $({\mathbf{e}}_x+{\mathbf{e}}_z)/\sqrt{2}$ and $({\mathbf{e}}_y+{\mathbf{e}}_z)/\sqrt{2}$ to the three positive coordinate axes instead gives matrix rank six. For example,This establishes uniqueness in the stated noise-free model. It does not establish small measurement uncertainty or an optimal experimental design.

Jones also interviewed John E. Tanner about the physical origins of pulsed field gradient NMR. The interview preserved concrete details of the apparatus, including Tanner's carefully spaced mica stack and the practical difficulties of rapidly switching magnetic field gradients.^[ref]

Diffusion MRI gradient sampling is the choice of magnetic field gradient directions, strengths and repetitions used to acquire diffusion magnetic resonance imaging measurements. Under a Gaussian diffusion tensor model, different directions probe different values of a quadratic form. Acquiring many measurements in an uninformative set of directions cannot determine every tensor coefficient.

For a unit gradient direction $\mathbf{g}$, a separately measured matched baseline $S_0$ and fixed positive diffusion weighting b-value $b$, the ideal signal model isThe symmetric matrix $\mathbf{D}$ has six independent coefficients. Opposite directions $\mathbf{g}$ and $-\mathbf{g}$ give the same model prediction. Thus counting orientations without considering their independent information can overstate diffusion tensor identifiability.

A 2003 Monte Carlo simulation study controlled the observation budget: one baseline plus60 diffusion-weighted observations. Six distinct orientations repeated ten times and30 distinct orientations repeated twice each use60 diffusion-weighted observations. The comparison asks how angular coverage and repetition affect estimated orientation, rather than simply comparing a longer acquisition with a shorter one.^[ref]

The subsequent 2004 study compared direction schemes through estimated orientation, trace and fractional anisotropy at $b=1000s/m{\mathrm{m}}^2$. These are distinct outputs of the same tensor fit. An acquisition suited to one readout need not minimize uncertainty in another.^[ref]

An original design comparison separates two questions. Repeating the three coordinate axes may reduce random variation in their measured attenuations, but it adds no equations for the three off-diagonal tensor entries. Adding suitable mixed directions can supply those equations. Full matrix rank establishes uniqueness in the ideal model; practical measurement uncertainty, motion and model mismatch still require assessment.