#! /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`; }