ChristianGaser / cat12

Computational Anatomy Toolbox for SPM12
https://neuro-jena.github.io/cat
GNU General Public License v2.0
33 stars 5 forks source link

Surface recon works for LH but fails for RH #4

Closed arkhanna closed 1 year ago

arkhanna commented 1 year ago

Hi all -- having an odd issue with running CAT12 for use in LEAD-DBS. Surface reconstruction works just fine for LH but fails for RH, claiming no WM voxel exists. Running on Mac 13.0.1 via MATLAB R2022b. Tried with numerous datasets so doubt it's a scan quality issue. Has anyone encountered something like this? Any ideas?

CAT12.8.2 r2137

Amap using initial SPM12 segmentations (MRF filter strength 0.07)    50s
    AMAP peaks: [CSF,GM,WM] = [0.36±0.09,0.68±0.06,0.97±0.06]
Final cleanup (gcutstr=0.25)                                      
  Level 1 cleanup (ROI estimation)                                    4s
  Level 1 cleanup (brain masking)                                     4s
  Level 2 cleanup (CSF correction)                                    2s
  Level 3 cleanup (CSF/WM PVE)                                        3s
                                                                     13s
Shooting registration is not required.
Write result maps                                                     1s
Surface and thickness estimation                                  
lh:
  Thickness estimation (0.50 mm³)                                 
    WM distance:                                                     35s
    CSF distance:                                                    15s
    PBT2x thickness:                                                 46s
                                                                    109s
  Create initial surface                                            158s
  Topology correction:                                              696s
  Surface refinement:                                               283s
  Reduction of surface collisions with optimization:                280s
  Spherical mapping with areal smoothing                            141s
  Spherical registration                                            229s
rh:
  Thickness estimation (0.50 mm³)                                 
    WM distance:                                                  WARNING:cat_vol_eidist: Found no object for distance estimation!
WARNING:cat_vol_eidist: Found no object for distance estimation!
    0s
    CSF distance:                                                 WARNING:cat_vol_eidist: Found no object for distance estimation!
WARNING:cat_vol_eidist: Found no object for distance estimation!
    0s
    PBT2x thickness:                                              
------------------------------------------------------------------------
CAT Preprocessing error for anat_t1:
------------------------------------------------------------------------
ERROR: no WM voxel

------------------------------------------------------------------------
  326 - cat_vol_pbt
  506 - cat_surf_createCS2
  795 - cat_main1639
 1269 - cat_run_job1639
   40 - cat_run_newcatch
  341 - run_job
ChristianGaser commented 1 year ago

Hi,

this is indeed weird. Please can you send me a PDF of the report file. Furthermore, I would try the newest CAT12 version which is also tested on Matlab 2022b.

Christian

ChristianGaser commented 1 year ago

I have only recently seen this error and it was caused by a wrong TPM.nii, probably replacing the original TPM.nii that comes with SPM12.