SIENA - Analysis of Structural Brain MRI DataSIENA - Structural Image Evaluation, using Normalisation, of Atrophy - v2.6intro - tools used - SIENA - SIENAX - voxelwise SIENA statistics | ![]() |
SIENA is a package for both single-time-point ("cross-sectional") and two-time-point ("longitudinal") analysis of brain change, in particular, the estimation of atrophy (volumetric loss of brain tissue). SIENA has been used in many clinical studies.
siena estimates percentage brain volume change (PBVC) betweem two input images, taken of the same subject, at different points in time. It calls a series of FSL programs to strip the non-brain tissue from the two images, register the two brains (under the constraint that the skulls are used to hold the scaling constant during the registration) and analyse the brain change between the two time points. It is also possible to project the voxelwise atrophy measures into standard space in a way that allows for multi-subject voxelwise statistical testing.
sienax estimages total brain tissue volume, from a single image, normalised for skull size. It calls a series of FSL programs: It first strips non-brain tissue, and then uses the brain and skull images to estimate the scaling between the subject's image and standard space. It then runs tissue segmentation to estimate the volume of brain tissue, and multiplies this by the estimated scaling factor, to reduce head-size-related variability between subjects.
For more detail on SIENA and technical reports, see the SIENA web page.
If you use SIENA in your research, please make sure that you reference the following articles. You may alternatively wish to use the brief descriptive methods text and expanded list of references given below.
Two-timepoint percentage brain volume change was estimated with SIENA [Smith 2002], part of FSL [Smith 2004].
[Smith 2002] S.M. Smith, Y. Zhang, M. Jenkinson, J. Chen, P.M. Matthews, A. Federico, and N. De Stefano.
[Smith 2004] S.M. Smith, M. Jenkinson, M.W. Woolrich, C.F. Beckmann, T.E.J. Behrens, H. Johansen-Berg, P.R. Bannister, M. De Luca, I. Drobnjak, D.E. Flitney, R. Niazy, J. Saunders, J. Vickers, Y. Zhang, N. De Stefano, J.M. Brady, and P.M. Matthews.
or
Brain tissue volume, normalised for subject head size, was estimated with SIENAX [Smith 2002], part of FSL [Smith 2004].
Accurate, robust and automated longitudinal and cross-sectional brain change analysis.
NeuroImage, 17(1):479-489, 2002.
Advances in functional and structural MR image analysis and implementation as FSL.
NeuroImage, 23(S1):208-219, 2004.
SIENA. Two-timepoint percentage brain volume change was estimated with SIENA [Smith 2001, Smith 2002], part of FSL [Smith 2004]. SIENA starts by extracting brain and skull images from the two-timepoint whole-head input data [Smith 2002b]. The two brain images are then aligned to each other [Jenkinson 2001, Jenkinson 2002] (using the skull images to constrain the registration scaling); both brain images are resampled into the space halfway between the two. Next, tissue-type segmentation is carried out [Zhang 2001] in order to find brain/non-brain edge points, and then perpendicular edge displacement (between the two timepoints) is estimated at these edge points. Finally, the mean edge displacement is converted into a (global) estimate of percentage brain volume change between the two timepoints.
SIENAX. Brain tissue volume, normalised for subject head size, was estimated with SIENAX [Smith 2001, Smith 2002], part of FSL [Smith 2004]. SIENAX starts by extracting brain and skull images from the single whole-head input data [Smith 2002b]. The brain image is then affine-registered to MNI152 space [Jenkinson 2001, Jenkinson 2002] (using the skull image to determine the registration scaling); this is primarily in order to obtain the volumetric scaling factor, to be used as a normalisation for head size. Next, tissue-type segmentation with partial volume estimation is carried out [Zhang 2001] in order to calculate total volume of brain tissue (including separate estimates of volumes of grey matter, white matter, peripheral grey matter and ventricular CSF).
Voxelwise multi-subject SIENA statistics. First, SIENA was
run separately for each subject. Next, for each subject, the edge
displacement image (encoding, at brain/non-brain edge points, the
outwards or inwards edge change between the two timepoints) was
dilated, transformed into MNI152 space, and masked by a standard
MNI152-space brain edge image. In this way the edge displacement
values were warped onto the standard brain edge [Bartsch 2004]. Next,
the resulting images from all subjects were fed into voxelwise
statistical analysis to test for.....
[Smith 2001] S.M. Smith, N. De Stefano, M. Jenkinson, and P.M. Matthews.
[Smith 2002] S.M. Smith, Y. Zhang, M. Jenkinson, J. Chen, P.M. Matthews, A. Federico, and N. De Stefano.
[Smith 2004] S.M. Smith, M. Jenkinson, M.W. Woolrich, C.F. Beckmann, T.E.J. Behrens, H. Johansen-Berg, P.R. Bannister, M. De Luca, I. Drobnjak, D.E. Flitney, R. Niazy, J. Saunders, J. Vickers, Y. Zhang, N. De Stefano, J.M. Brady, and P.M. Matthews.
[Smith 2002b] S.M. Smith.
[Jenkinson 2001] M. Jenkinson and S.M. Smith.
[Jenkinson 2002] M. Jenkinson, P.R. Bannister, J.M. Brady, and S.M. Smith.
[Zhang 2001] Y. Zhang, M. Brady, and S. Smith.
[Bartsch 2004] A.J. Bartsch, N. Bendszus, N. De Stefano, G. Homola, and S. Smith.
Normalised accurate measurement of longitudinal brain change.
Journal of Computer Assisted Tomography, 25(3):466-475, May/June 2001.
Accurate, robust and automated longitudinal and cross-sectional brain change analysis.
NeuroImage, 17(1):479-489, 2002.
Advances in functional and structural MR image analysis and
implementation as FSL.
NeuroImage, 23(S1):208-219, 2004.
Fast robust automated brain extraction.
Human Brain Mapping, 17(3):143-155, November 2002.
A global optimisation method for robust affine registration of brain images.
Medical Image Analysis, 5(2):143-156, June 2001.
Improved optimisation for the robust and accurate linear registration and motion correction of brain images.
NeuroImage, 17(2):825-841, 2002.
Segmentation of brain MR images through a hidden Markov random field model and the expectation maximization algorithm.
IEEE Trans. on Medical Imaging, 20(1):45-57, 2001.
Extending SIENA for a multi-subject statistical analysis of sample-specific cerebral edge shifts: Substantiation of early brain regeneration through abstinence from alcoholism.
In Tenth Int. Conf. on Functional Mapping of the Human Brain, 2004.
bet - Brain Extraction Tool. This automatically removes all non-brain tissue from the image. It can optionally output the binary brain mask that was derived during this process, and output an estimate of the external surface of the skull, for use as a scaling constraint in later registration.
pairreg, a script supplied with FLIRT - FMRIB's Linear Image Registration Tool. This script calls FLIRT with a special optimisation schedule, to register two brain images whilst at the same time using two skull images to hold the scaling constant (in case the brain has shrunk over time, or the scanner calibration has changed). The script first calls FLIRT to register the brains as fully as possible. This registration is then applied to the skull images, but only the scaling and skew are allowed to change. This is then applied to the brain images, and a final pass optimally rotates and translates the brains to get the best final registration.
fast - FMRIB's Automated Segmentation Tool. This program automatically segments a brain-only image into different tissue types (normally background, grey matter, white matter, CSF and other). It also corrects for bias field. It is used in various ways in the SIENA scripts. Note that both siena and sienax allow you to choose between segmentation of grey matter and white matter as separate classes or a single class. It is important to choose the right option here, depending on whether there is or is not reasonable grey-white contrast in the image.
A default SIENA analysis is run by typing:
siena <input1> <input2>
The input filenames must not contain directory names - i.e. all must be done within a single directory.
Other options are:
-o <output-dir> : set output directory (the default output is <input1>_to_<input2>_siena)
-d : debug (don't delete intermediate files)
-B "bet options" : if you want to change the BET defaults,
put BET options inside double-quotes after using the -B flag. For
example, to increase the size of brain estimation, use: -B "-f
0.3"
-2 : two-class segmentation (don't segment grey and white matter separately) - use this if there is poor grey/white contrast
-t2: tell FAST that the input images are T2-weighted and not T1
-m : use standard-space masking as well as BET (e.g. if it is proving hard to get reliable brain segmentation from BET, for example if eyes are hard to segment out) - register to standard space in order to use a pre-defined standard-space brain mask
-t <t>: ignore from t (mm) upwards in MNI152/Talairach space - if you need to ignore the top part of the head (e.g. if some subjects have the top missing and you need consistency across subjects)
-b <b>: ignore from b (mm) downwards in MNI152/Talairach space; b should probably be -ve
-S "siena_diff options" : if you want to send options to the
siena_diff program (that estimates change between two aligned
images), put these options in double-quotes after the -S flag. For
example, to tell siena_diff to run FAST segmentation with an
increased number of iterations, use -S "-s -I 20"
siena carries out the following steps:
Run bet on the two input images, producing as output, for each input: extracted brain, binary brain mask and skull image. If you need to call BET with a different threshold than the default of 0.5, use -f <threshold>.
Run siena_flirt, a separate script, to register the two brain images. This first calls the FLIRT-based registration script pairreg (which uses the brain and skull images to carry out constrained registration). It then deconstructs the final transform into two half-way transforms which take the two brain images into a space halfway between the two, so that they both suffer the same amount of interpolation-related blurring. Finally the script produces a multi-slice gif picture showing the registration quality, with one transformed image as the background and edges from the other transformed image superimposed in red.
The final step is to carry out change analysis on the registered brain images. This is done using the program siena_diff. (In order to improve slightly the accuracy of the siena_diff program, a self-calibration script siena_cal, described later, is run before this.) siena_diff carries out the following steps:
report.siena
the SIENA log, including the final PBVC
estimate.
report.html
a webpage report including images showing
various stages of the analysis, the final result and a description of
the SIENA method.
A_halfwayto_B_render
a colour-rendered image of edge
motion superimposed on the halfway A image. Red-yellow means brain
volume increase and Blue means brain volume decrease ("atrophy").
A_and_B.gif
a gif image showing the results of the registration,
using one transformed image as the background and the other as the
coloured edges foreground.
A_to_B.mat
the transformation taking A to B, using the brain and
skull images.
B_to_A.mat
the transformation taking B to A, using the brain and
skull images.
A_halfwayto_B.mat
and B_halfwayto_A.mat
the transformations taking
the images to the halfway positions.
A default SIENAX analysis is run by typing:
sienax <input>
The input filename must not contain directory names - i.e. all must be done within the current directory.
Other options are:
-o <output-dir> : set output directory (the default output is <input>_sienax)
-d : debug (don't delete intermediate files)
-B "bet options" : if you want to change the BET defaults,
put BET options inside double-quotes after using the -B flag. For
example, to increase the size of brain estimation, use: -B "-f
0.3"
-2 : two-class segmentation (don't segment grey and white matter separately) - use this if there is poor grey/white contrast
-t2: tell FAST that the input images are T2-weighted and not T1
-t <t>: ignore from t (mm) upwards in MNI152/Talairach space - if you need to ignore the top part of the head (e.g. if some subjects have the top missing and you need consistency across subjects)
-b <b>: ignore from b (mm) downwards in MNI152/Talairach space; b should probably be -ve
-r: tell SIENAX to estimate "regional" volumes as well as global; this produces peripheral cortex GM volume (3-class segmentation only) and ventricular CSF volume
-lm <mask>: use a lesion (or lesion+CSF) mask to remove incorrectly labelled "grey matter" voxels
-S "FAST options" : if you want to change the segmentation
defaults, put FAST options inside double-quotes after using the -S
flag. For example, to increase the number of segmentation
iterations use: -S "-I 20"
sienax carries out the following steps:
report.sienax
the SIENAX log, including the final
volume estimates.
report.html
a webpage report including images showing
various stages of the analysis, the final result and a description of
the SIENAX method.
I_render
a colour-rendered image showing the segmentation
output superimposed on top of the original image.
We have extended SIENA to allow the voxelwise statistical analysis of atrophy across subjects. This takes a SIENA-derived edge "flow image" (edge displacement between the timepoints) for each subject, warps these to align with a standard-space edge image and then carries out voxelwise cross-subject statistical analysis to identify brain edge points which, for example, are signficantly atrophic for the group of subjects as a whole, or where atrophy correlates significantly with age or disease progression.
In order to carry out voxelwise SIENA statistics, do the following:
siena A B
A
and B
).
cd <siena_output_directory>
siena_flow2std A B
siena
, dilates this several times (to "thicken" this edge
flow image), transforms to standard space, and masks with a standard
space edge mask. It then smooths this with a default Gaussian filter of
half-width 5mm before remasking. If you want to change the smoothing
then use the -s option; set the smoothing to zero to turn if off
completely.
A_to_B_flow_to_std
. Merge these into a
single 4D image; for example, if each subject's analysis has so far
been carried out in a subdirectory called subject_*/A_to_B_siena, where the *
could be subject ID or name, use a command such as: fslmerge
-t flow_all_subjects `imglob
subject_*/A_to_B_siena/A_to_B_flow_to_std*`
design.mat
and contrasts file
design.con
. The mask image that you use for randomise
should be ${FSLDIR}/data/standard/MNI152_T1_2mm_edges