BUG: adapt to numpy > 1.8

np.diag is not writeable in recent numpy versions
nipy · Mar 11, 2014 · 548f089 · 548f089
1 parent 83668cb
commit 548f089
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/nipy/labs/datasets/volumes/volume_img.py b/nipy/labs/datasets/volumes/volume_img.py
@@ -277,7 +277,7 @@ def xyz_ordered(self, resample=False, copy=True):
             axis_numbers = np.argmax(np.abs(A), axis=0)
 
         # Now make sure the affine is positive
-        pixdim = np.diag(A)
+        pixdim = np.diag(A).copy()
         data = img.get_data()
         if pixdim[0] < 0:
             b[0] = b[0] + pixdim[0]*(data.shape[0] - 1)