Merge pull request #2017 from MRtrix3/fix_matlab_write_mrtrix

write_mrtrix.m: fix handling of .mih outputs

Author: jdtournier

matlab/read_mrtrix.m

@@ -100,7 +100,7 @@
 image.data = ipermute (image.data, order);
 for i=1:size(order,2)
   if layout{i}(1) == '-'
-    image.data = flipdim(image.data, i);
+    image.data = flip(image.data, i);
   end
 end
 

matlab/write_mrtrix.m

@@ -122,26 +122,24 @@ function write_mrtrix (image, filename)
     end
   end
 end
- 
+
+fprintf (fid, '\nfile: ')
 
 if strcmp(filename(end-3:end), '.mif')
-  datafile = filename;
-  dataoffset = ftell (fid) + 24;
-  fprintf (fid, '\nfile: . %d\nEND\n                         ', dataoffset);
+  dataoffset = ftell (fid) + 18;
+  dataoffset += mod((4 - mod(dataoffset, 4)), 4);
+  fprintf (fid, '. %d\nEND\n              ', dataoffset);
+  fseek (fid, dataoffset);
 elseif strcmp(filename(end-3:end), '.mih')
-  datafile = [ filename(end-3:end) '.dat' ];
-  dataoffset = 0;
-  fprintf (fid, '\nfile: %s %d\nEND\n', datafile, dataoffset);
+  datafile = [ filename(1:end-4) '.dat' ];
+  fprintf (fid, '%s 0\nEND\n', datafile);
+  fclose(fid);
+  fid = fopen (datafile, 'w', byteorder);
 else
   fclose(fid);
   error('unknown file suffix - aborting');
 end
 
-fclose(fid);
-
-fid = fopen (datafile, 'r+', byteorder);
-fseek (fid, dataoffset, -1);
-
 if isstruct(image)
   fwrite (fid, image.data, precision);
 else