#!/bin/sh

sizex=`fslval $1 pixdim1`
sizey=`fslval $1 pixdim2`
sizez=`fslval $1 pixdim3`

voxvol=`echo "$sizex * $sizey * $sizez" | bc -l` 
echo "Voxel volume is $voxvol (in mm^3)"

nx=`fslval $1 dim1`
ny=`fslval $1 dim2`
nz=`fslval $1 dim3`

imvol=`echo "$nx * $ny * $nz * $voxvol" | bc -l`
echo "Image volume is $imvol (in mm^3)"

fovx=`echo "$sizex * $nx" | bc -l`
fovy=`echo "$sizey * $ny" | bc -l`
fovz=`echo "$sizez * $nz" | bc -l`

echo "FOV is $fovx by $fovy by $fovz (in mm)"