input-tables.sh

Go to the documentation of this file.

#!/bin/zsh
script=${0##*/}

# --------------------------------------------------------------------------------------------
#
# Auxiliary script to convert csv files from input-tables to input files.
#
# --------------------------------------------------------------------------------------------

if [ -z $JPP_DIR ]; then
    echo "Variable JPP_DIR undefined."
    exit
fi

source $JPP_DIR/setenv.sh $JPP_DIR >& /dev/null

zmodload zsh/mathfunc

set_variable  DEBUG          2
set_variable  WORKDIR        ${TMPDIR:-/tmp}

if do_usage $*; then
    usage "$script <directory> <detector identifier>"\
          "\nThe directory corresponds to the GIT repositiry \"calibration/input-tables\""\
          "\nand the detector identifier to the DETOID."
fi

if (( $# != 2 )); then
    fatal "Wrong number of arguments."
fi

set_variable  DIR           $1
set_variable  DETOID        $2

if [[ ! -d $DIR/data/${DETOID} ]]; then
    fatal "No directory $DIR/data/${DETOID}"
fi

#
# Create tripod file.
#
# \param  1    input file
# \param  2    output file
#
function tripod()
{
    TRIPOD_CSV=$1
    TRIPOD_TXT=$2

    if [[ -f $TRIPOD_CSV ]]; then

        printf "# Data from %s\n" ${TRIPOD_CSV:t}                  >  $TRIPOD_TXT
        grep -E -v "^#" $TRIPOD_CSV  |  tail -n +2  |  tr ',' ' '  >> $TRIPOD_TXT
    fi
}

#
# Create hydrophone file.
#
# \param  1    input file
# \param  2    output file
#
function hydrophone()
{
    HYDROPHONE_CSV=$1
    HYDROPHONE_TXT=$2

    if [[ -f $HYDROPHONE_CSV ]]; then

        grep -E -v "^#" $HYDROPHONE_CSV  |  tail -n +2  |  tr ',' ' '  >  $WORKDIR/hydrophones.txt

        printf "# Data from %s\n" ${HYDROPHONE_CSV:t}                  >  $HYDROPHONE_TXT

        while read -A BUFFER; do
            
            set_variable STRING     $BUFFER[1]
            set_variable FLOOR      $BUFFER[2]
            set_variable X          $BUFFER[3]
            set_variable Y          $BUFFER[4]
            set_variable ANGLE_DEG  $BUFFER[5]
            set_variable Z          $BUFFER[6]
            
            let "X0 = -1.0 * $Y"
            let "Y0 = -1.0 * $X"

            let "CT = cos($ANGLE_DEG * $PI / 180.0)"
            let "ST = sin($ANGLE_DEG * $PI / 180.0)"

            let "X1 = $ST * $X0  +  $CT * $Y0"
            let "Y1 = $CT * $X0  -  $ST * $Y0"
            
            printf "%4d %2d %6.2f %6.2f %6.2f\n" $STRING $FLOOR $X1 $Y1 $Z  >>  $HYDROPHONE_TXT

        done < $WORKDIR/hydrophones.txt

        rm -f  $WORKDIR/hydrophones.txt
    fi
}

#
# Create transmitter file.
#
# \param  1    input file
# \param  2    output file
#
function transmitter()
{
    TRANSMITTER_CSV=$1
    TRANSMITTER_TXT=$2

    if [[ -f $TRANSMITTER_CSV ]]; then

        grep -E -v "^#" $TRANSMITTER_CSV  |  tail -n +2  |  tr ',' ' '  >  $WORKDIR/transmitters.txt

        printf "# Data from %s\n" ${TRANSMITTER_CSV:t}                  >  $TRANSMITTER_TXT

        while read -A BUFFER; do

            set_variable ID         $BUFFER[1]
            set_variable STRING     $BUFFER[2]
            set_variable FLOOR      $BUFFER[3]
            set_variable X          $BUFFER[4]
            set_variable Y          $BUFFER[5]
            set_variable ANGLE_DEG  $BUFFER[6]
            set_variable Z          $BUFFER[7]

            let "X0 = -1.0 * $Y"
            let "Y0 = -1.0 * $X"

            let "CT = cos($ANGLE_DEG * $PI / 180.0)"
            let "ST = sin($ANGLE_DEG * $PI / 180.0)"

            let "X1 = $ST * $X0  +  $CT * $Y0"
            let "Y1 = $CT * $X0  -  $ST * $Y0"
            
            printf "%2d %4d %2d %6.2f %6.2f %6.2f\n" $ID $STRING $FLOOR $X1 $Y1 $Z  >>  $TRANSMITTER_TXT

        done < $WORKDIR/transmitters.txt

        rm -f  $WORKDIR/transmitters.txt
    fi
}


tripod       $DIR/data/${DETOID}/${DETOID}_tripods.csv             tripod.txt
hydrophone   $DIR/data/${DETOID}/${DETOID}_hydrophones.csv         hydrophone.txt
transmitter  $DIR/data/${DETOID}/${DETOID}_baseAcousticBeacon.csv  transmitter.txt