#! /usr/bin/env perl

use strict;

my @files = glob("*.mnc");
my $backupdir = "old_mncs";

print "mkdir $backupdir\n";
`mkdir $backupdir`;

foreach my $file (@files){

 print "cp -p $file $backupdir/$file\n";
 `cp -p $file $backupdir/$file`;

  my $dx = `mincinfo -attvalue xspace:step $file`; chomp( $dx );
  my $dy = `mincinfo -attvalue yspace:step $file`; chomp( $dy );
  my $dz = `mincinfo -attvalue zspace:step $file`; chomp( $dz );
  my $nx = `mincinfo -dimlength xspace $file`; chomp( $nx );
  my $ny = `mincinfo -dimlength yspace $file`; chomp( $ny );
  my $nz = `mincinfo -dimlength zspace $file`; chomp( $nz );
  my $sx = `mincinfo -attvalue xspace:start $file`; chomp( $sx );
  my $sy = `mincinfo -attvalue yspace:start $file`; chomp( $sy );
  my $sz = `mincinfo -attvalue zspace:start $file`; chomp( $sz );

  my $xstart = -0.5 * $nx * $dx;
  my $ystart = -0.5 * $ny * $dy;
  my $zstart = -0.5 * $nz * $dz;

  print "minc_modify_header -dinsert xspace:start=$xstart -dinsert yspace:start=$ystart -dinsert zspace:start=$zstart $file\n";
  `minc_modify_header -dinsert xspace:start=$xstart -dinsert yspace:start=$ystart -dinsert zspace:start=$zstart $file`;
}