(name suggestions welcome)

As raised in multiple locations:

A convenience script for concatenating multiple DWI series where the intensity scaling of image reconstruction may have varied due to recalibration. This would involve:

• Extracting the b=0 volumes from each series.
  (Easiest would be to require .mifs at input so that gradient tables are embedded)

• Using mrhistmatch scale or mrstats -output median to determine the global intensity fluctuations between acquisitions.

• Concatenate all data into a single 4D image, applying appropriate multiplicative factors.

It should also check for variations in TE and TR between the input series, and warn where appropriate.