Hi, Clement-

Glad to know the purpose of this---de-sphinxifying.

The reason for 3dresample is that I find it easier to think about using the change in dataset orientation for the purpose of desphinxifying when starting from a "simple" orientation where the L-R axis is in the spot of the x-coord, the A-P axis in that of the y-coord and the I-S axis in that of the z-coord. 3dresample doesn't do any altering of the data's appearance, it just sets the header to be more easily visualizable. I still find that useful, from a human point of view. Also, the reason I could recommend a possible way to desphinxify is because I tend to 'reset' the orientation to RAI before trying to use orientation to change the way things appear.

So, taking things in steps, and integrating the orientation you said was useful:

# copy the data and remove obliquity
3dcopy DSET_00 DSET_01
3drefit -deoblique DSET_01

# reset the orientation to something simple
3dresample -orient RAI -input DSET_01 -prefix DSET_02

# change the header orientation information to rearrange how the dataset appears
3dcopy DSET_02 DSET_03
3drefit -orient LIP DSET_03

# get a reasonable center of coordinate grid in the center-of-mass of the dset
3dCM -set 0 0 0 -automask DSET_03

--pt