#!/bin/bash

# Wrapper for modelfit that provides a simple interface for robustly
# fitting the diffusion tensor using the RESTORE algorithm
# 
# $Id: restore,v 1.1 2008/12/08 17:48:43 bennett Exp $


# Find the directory containing this script.
SCRIPTDIR=${0%/*}

if [ $# -lt 3 ]; then
    echo "restore <data file> <scheme file> <noise std> [<outlier map file>] [options]"
    exit
fi

DATAFILE=$1
SCHEMEFILE=$2
NOISESTD=$3
OUTLIERMAPFILE=$4

OUTLIERMAPARG=

if [ $# -gt 3 ] && [ ${OUTLIERMAPFILE:0:1} != "-" ]; then

    OUTLIERMAPARG="-outliermap $OUTLIERMAPFILE"
    shift 4

else

    shift 3

fi

DATAOPT="-inputfile $DATAFILE"
if [ $DATAFILE == - ]; then
    DATAOPT=
fi


exec ${SCRIPTDIR}/modelfit $* -inversion -2 -schemefile $SCHEMEFILE $DATAOPT $OUTLIERMAPARG -sigma $NOISESTD 


