Fixed concat_images so it can work on large sequences of images.

Author: Yannick Schwartz

nibabel/funcs.py

@@ -11,6 +11,7 @@
 
 from .orientations import (io_orientation, orientation_affine, flip_axis,
                            apply_orientation, OrientationError)
+from .loadsave import load
 
 
 def squeeze_image(img):
@@ -92,7 +93,7 @@ def concat_images(images, check_affines=True):
     Parameters
     ----------
     images : sequence
-       sequence of ``SpatialImage``\s
+       sequence of paths to be loaded as ``SpatialImage``\s
     check_affines : {True, False}, optional
        If True, then check that all the affines for `images` are nearly
        the same, raising a ``ValueError`` otherwise.  Default is True
@@ -104,13 +105,14 @@ def concat_images(images, check_affines=True):
        dimension
     '''
     n_imgs = len(images)
-    img0 = images[0]
+    img0 = load(images[0])
     i0shape = img0.shape
     affine = img0.get_affine()
     header = img0.get_header()
     out_shape = (n_imgs, ) + i0shape
     out_data = np.empty(out_shape)
     for i, img in enumerate(images):
+        img = load(img)
         if check_affines:
             if not np.all(img.get_affine() == affine):
                 raise ValueError('Affines do not match')