AFNI Message Board

Dear AFNI users-

We are very pleased to announce that the new AFNI Message Board framework is up! Please join us at:

https://discuss.afni.nimh.nih.gov

Existing user accounts have been migrated, so returning users can login by requesting a password reset. New users can create accounts, as well, through a standard account creation process. Please note that these setup emails might initially go to spam folders (esp. for NIH users!), so please check those locations in the beginning.

The current Message Board discussion threads have been migrated to the new framework. The current Message Board will remain visible, but read-only, for a little while.

Sincerely, AFNI HQ

History of AFNI updates  

Previous Message Next Message
Message List New Topic Print View
|
Damien Mannion
October 04, 2012 05:24AM
Hi,

I'm trying to run 3dDFT on a surface file but I'm running into some problems. Running 3dDFT prints a lot of '** GBMMP, bad dtype' errors to the screen, but outputs a file OK. If I load the file into SUMA, it looks like the ABS frequency spectrum I would expect. If I try and extract the phase using 3dcalc, though, the resulting file looks identical to the ABS. I'm wondering if 3dDFT is supported on surface files?

Thanks.

More details-
The input file: 
++ 3dinfo: AFNI version=AFNI_2011_12_21_1014 (Sep 25 2012) [64-bit]

Dataset File:    lh_func.niml.dset
Identifier Code: XYZ_Dqxwiex5g9hn6HSl-GwsHg  Creation Date: Thu Oct  4 18:02:18 2012
Template Space:  
Dataset Type:    Echo Planar (-epan)
Byte Order:      LSB_FIRST {assumed} [this CPU native = LSB_FIRST]
Storage Mode:    NI_SURF_DSET
Storage Space:   26,029,200 (26 million) bytes
Data Axes Tilt:  Unspecified, assumed plumb
Data Axes Orientation:
  first  (x) = Right-to-Left
  second (y) = Anterior-to-Posterior
  third  (z) = Inferior-to-Superior   [-orient RAI]
R-to-L extent:     0.000     -to- 65072.000 [L] -step-     1.000 mm [65073 voxels]
A-to-P extent:     0.000     -to-     0.000     -step-     1.000 mm [  1 voxels]
I-to-S extent:     0.000     -to-     0.000     -step-     1.000 mm [  1 voxels]
Number of time steps = 100  Time step = 3.00000s  Origin = 0.00000s
  -- At sub-brick #0 '#7' datum type is float
  -- At sub-brick #1 '#8' datum type is float
  -- At sub-brick #2 '#9' datum type is float
** For info on all 100 sub-bricks, use '3dinfo -verb' **

----- HISTORY -----
[fmri@mannion-fmri: Thu Oct  4 18:02:15 2012] 3dVol2Surf -spec /home/fmri/stemp/k1002/SUMA/k1002_lh.spec -surf_A smoothwm -surf_B pial -map_func ave -f_steps 15 -f_index nodes -sv k1002_SurfVol_Alnd_Exp+orig -grid_parent func+orig -out_niml lh_func.niml.dset

The 3dDFT: 
3dDFT -prefix lh_dft.niml.dset lh_func.niml.dset 
++ Authored by: Kevin Murphy & Zhark the Transformer
++ Data length = 100 ; FFT length = 100
++ Output complex-valued dataset: ./lh_dft.niml.dset
** GBMMP, bad dtype
#  (... * 100 ) trimmed by DJM

The 3dDFT output: 
++ 3dinfo: AFNI version=AFNI_2011_12_21_1014 (Sep 25 2012) [64-bit]

Dataset File:    lh_dft.niml.dset
Identifier Code: XYZ_E00GTjyzbF5c9DQxiTQ54A  Creation Date: Thu Oct  4 18:04:43 2012
Template Space:  
Dataset Type:    Echo Planar (-epan)
Byte Order:      LSB_FIRST {assumed} [this CPU native = LSB_FIRST]
Storage Mode:    NI_SURF_DSET
Storage Space:   26,029,200 (26 million) bytes
Data Axes Tilt:  Unspecified, assumed plumb
Data Axes Orientation:
  first  (x) = Right-to-Left
  second (y) = Anterior-to-Posterior
  third  (z) = Inferior-to-Superior   [-orient RAI]
R-to-L extent:     0.000     -to- 65072.000 [L] -step-     1.000 mm [65073 voxels]
A-to-P extent:     0.000     -to-     0.000     -step-     1.000 mm [  1 voxels]
I-to-S extent:     0.000     -to-     0.000     -step-     1.000 mm [  1 voxels]
Number of time steps = 100  Time step = 0.00333s  Origin = 0.00000s
  -- At sub-brick #0 '#7' datum type is float
  -- At sub-brick #1 '#8' datum type is float
  -- At sub-brick #2 '#9' datum type is float
** For info on all 100 sub-bricks, use '3dinfo -verb' **

----- HISTORY -----
[fmri@mannion-fmri: Thu Oct  4 18:03:10 2012] 3dDFT -prefix lh_dft.niml.dset lh_func.niml.dset;

The 3dcalc command: 
3dcalc -cx2r PHASE -a lh_dft.niml.dset -expr a -prefix lh_phase.niml.dset
Reply Quote
RSS
Subject Author Posted

3dDFT on surfaces

Damien Mannion October 04, 2012 05:24AM

Re: 3dDFT on surfaces

ziad October 05, 2012 09:13AM

Re: 3dDFT on surfaces

Damien Mannion October 05, 2012 07:02PM

Re: 3dDFT on surfaces

ziad October 09, 2012 10:58AM

Re: 3dDFT on surfaces

Damien Mannion October 09, 2012 07:19PM