[burgewk]

Thanks for the quick reply!

[burgewk]

It worked, thanks!!

[burgewk]

Ok I got through the svreg pipeline without error. I have also run my test subject through the thicknessPVC.sh pipeline.

I tried to run through the generate_stats_xls_manual file and it works, but I’m not really sure what I should be feeding it.

subbasename=’/Storage/session_1/anat_1/mprage’;
vol_label_file=’/Storage/1_LH_V1.label.nii’; % binary volume file registered to .bfc
ROI_IDs= [1];
surf_ext=’svreg.dfs’; %[subbasename,’.left.mid.cortex.’,surf_ext]
output_ext=’testout’; %[subbasename,’.roiwise.’,output_ext,’.stats.txt’]
generate_stats_xls_manual(subbasename,vol_label_file,ROI_IDs,surf_ext,output_ext);

I’m not really sure what I should be putting for ROI_IDS. I’ve tried several different combinations of numbers. I thought 1 would be correct as the mask is binary. For ROI_IDs = 1 I get:

ROI_ID Mean_Thickness(mm) GM_Volume(mm^3) CSF_Volume(mm^3) WM_Volume(mm^3) Total_Volume(GM+WM)(mm^3) Cortical_Area_mid(mm^2) Cortical_Area_inner(mm^2) Cortical_Area_pial(mm^2)
1 NaN 219.089869 7.097567 9.160284 228.250153 0.000000 0.000000 0.000000

So it has a GM volume, but no thickness value?

Should I be setting my binary mask not equal to 1? I had a look around the xml file with ROI Ids but I don’t really know what I should be setting it to.

Thanks for the help,

Wes

• This reply was modified 8 years, 6 months ago by burgewk.

[burgewk]

I think I figured it out. I tried to run it from the command line and found it to be a X11 issue. I had to download xquartz and it appears to be running the svreg now. I am using OSX v 10.10.5 for anyone else who might encounter this.

Wes

[burgewk]

Here is the output when trying to run the GUI based svreg

Executing: /Applications/BrainSuite15b/svreg/bin/svreg_label_surf_hemi.sh /Volumes/Storage_backup/session_1/anat_1/mprage /Applications/BrainSuite15b/svreg/BCI-DNI_brain_atlas/BCI-DNI_brain left -v1
Executing: /Applications/BrainSuite15b/svreg/bin/svreg_label_surf_hemi.sh /Volumes/Storage_backup/session_1/anat_1/mprage /Applications/BrainSuite15b/svreg/BCI-DNI_brain_atlas/BCI-DNI_brain right -v1
Finished surface labeling 0:00
error reading /Volumes/Storage_backup/session_1/anat_1/mprage.left.mid.cortex.svreg.init.dfs
error reading /Volumes/Storage_backup/session_1/anat_1/mprage.right.mid.cortex.svreg.init.dfs

[burgewk]

Hi Brain Suite community,

Things are going along well with this. I am now at the stage where I would like to extract the information from my ROIs that have now been registered and put into BS space. This is where I run into an error.

Anand suggested I use the utilities file: main.generate_stats_xls_manual.m which has this at line 28:

l=readdfs([subbasename,’.left.mid.cortex.’,surf_ext]);

However, after my cortical extraction I did not have a .mid in my folder. I tried to run the svreg but it also was looking for the .mid files. Is there some other step I should do to generate the .mid files?

Here is an ls of my output folder:

mprage.air mprage.left.pial.cortex.dfs
mprage.bfc.nii.gz mprage.mask.nii.gz
mprage.brain.dfs mprage.nii.gz
mprage.bse.nii.gz mprage.outer_skull.dfs
mprage.cerebrum.mask.nii.gz mprage.pial.cortex.dfs
mprage.cortex.dewisp.mask.nii.gz mprage.pvc.frac.nii.gz
mprage.cortex.scrubbed.mask.nii.gz mprage.pvc.label.nii.gz
mprage.cortex.tca.mask.nii.gz mprage.right.inner.cortex.dfs
mprage.hemi.label.nii.gz mprage.right.pial.cortex.dfs
mprage.init.cortex.mask.nii.gz mprage.scalp.dfs
mprage.inner.cortex.dfs mprage.skull.label.nii.gz
mprage.inner_skull.dfs mprage.warp

I appreciate the help!

Wes

[burgewk]

Hi Anand,

Thanks for getting back to me. I have a few more questions:

What sort of registration should I be doing to get my regions into BS space?

We are trying to confirm our findings from Freesurfer in a separate analysis stream. So we will have both a Freesurfer based reconstruction and a Brain Suite based reconstruction for each participant. The end goal is to extract the thickness values for each region that we are interested in. All of our regions are hand drawn in Freesurfer as labels which we can convert to volumes (.nii) easily enough. Should I be registering their Freesurfer based anatomy to one of the outputs from the ‘Cortical surface extraction sequence’? If so, which output should I aim for.

I hope the above was clear. Please let me know if you need additional info.

Best,

Wes

[Author]

Posts