GQI_GFA.nii.gz
R-to-L extent:  -254.000 [R] -to-    -0.000     -step-     2.000 mm [128 voxels]
ROI_ICMAP_GMI+orig
R-to-L extent:  -129.041 [R] -to-   124.959 [L] -step-     2.000 mm [128 voxels]

#!/bin/tcsh

set refset    = [whatever is a reference data set for you with correct grid/location/etc., such as your *GMI* file]
set all_dsis  = `ls GQI*nii.gz`  # or whatever your DSI-Studio-named files are-- take all of them and process in loop below

# ref set's orient
set refo      = `3dinfo -orient $refset`

foreach ff ( $all_dsis )

    echo "\n\n+++ Proc'ing dset: $ff\n\n"

    set gg   = $ff:gas/.nii.gz//
    set fout = ${gg}_${refo}.nii.gz

    # new orient for dset (this information is initially *correct* in
    # the input dset as it is, so we change both BRIK+HEAD info);
    # create a new dset whose name reflects the new orient
    3dresample             \
        -orient $refo      \
        -input  $ff        \
        -prefix $fout

    # now change origin of the dset (we are saying this information is
    # *incorrect* in the dset, so we just change the header; the
    # distiction may be a bit trivial, since the origin is juuust info
    # in the header, but there you have it); *overwrites* input dset.
    3drefit \
        -duporigin $refset \
        $fout
end