Hi,
I wish to split the putamen in an anterior part and a posterior part in >60 subjects which have been analyzed by the recon-all stream (FS version 5.3.0). Ideally, I would split the putamen in the atlas (e.g., at the coronal plane through the anterior commissure) and apply this "new" atlas to all subjects, one of the advantages being the bias that would be introduced by splitting the putamena manually in all subjects would be avoided. After some time trying to make this happen, I now think this could be achieved by adjusting RB_all_2008-03-26.gca and applying something like
"mri_ca_label -align subject_norm.mgz subject_talairach.m3z RB_all_2008-03-26.gca subject_aseg.auto_noCCseg.PUTAMENA_SEPARATED.mgz".
The problem is that till now I wasn't able to figure out how to correctly read the binary file RB_all_2008-03-26.gca. I tried to use the MATLAB function load_gca.m, but that seems to only work up to version 4 of the .gca file format (as GCA_VERSION is set to 4 in the function), and RB_all_2008-03-26.gca seems to be version 5. Also, when trying to read the file anyway in MATLAB using the code written for version 4 and below (by setting GCA_VERSION to 5.0), the code fails (the first gibbs_nlabels is already way too large, although I suspect the error occurs earlier).
My questions: 1. Do you think I am following the right track to achieve a sensible splitting of the putamena this way? Or should I try something else? 2. Is there a newer load_gca.m file around which is able to read a version 5 .gca file? Or could you point me to the difference in the design of the version 5 .gca file compared to that of the version 4 file so that I can rewrite load_gca.m?
Cheers,
Jelmer ________________________________ De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
Hi Jelmer
if you have manually split asegs you can use them as input to mri_ca_train to create a new gca that has the separated putamen. You'll need to add a couple of entries to FreeSurferColorLUT.txt to specify how to display and name the new entries as well
cheers Bruce
On Thu, 25 Feb 2016, Kok, JG (neuro) wrote:
Hi,
I wish to split the putamen in an anterior part and a posterior part in >60 subjects which have been analyzed by the recon-all stream (FS version 5.3.0). Ideally, I would split the putamen in the atlas (e.g., at the coronal plane through the anterior commissure) and apply this "new" atlas to all subjects, one of the advantages being the bias that would be introduced by splitting the putamena manually in all subjects would be avoided. After some time trying to make this happen, I now think this could be achieved by adjusting RB_all_2008-03-26.gca and applying something like
"mri_ca_label -align subject_norm.mgz subject_talairach.m3z RB_all_2008-03-26.gca subject_aseg.auto_noCCseg.PUTAMENA_SEPARATED.mgz".
The problem is that till now I wasn't able to figure out how to correctly read the binary file RB_all_2008-03-26.gca. I tried to use the MATLAB function load_gca.m, but that seems to only work up to version 4 of the .gca file format (as GCA_VERSION is set to 4 in the function), and RB_all_2008-03-26.gca seems to be version 5. Also, when trying to read the file anyway in MATLAB using the code written for version 4 and below (by setting GCA_VERSION to 5.0), the code fails (the first gibbs_nlabels is already way too large, although I suspect the error occurs earlier).
My questions:
- Do you think I am following the right track to achieve a sensible
splitting of the putamena this way? Or should I try something else? 2. Is there a newer load_gca.m file around which is able to read a version 5 .gca file? Or could you point me to the difference in the design of the version 5 .gca file compared to that of the version 4 file so that I can rewrite load_gca.m?
Cheers,
Jelmer
De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
Hi Bruce,
Thanks for your quick reply. I do not have manually split asegs, I wish to create them all automatically. I'm hoping it is possible to read RB_all_2008-03-26.gca (preferably in MATLAB, using fread), and adjust all putamen values posterior from the coronal plane through the anterior commissure from (FreeSurferColorLUT)values 12 and 51 to for example values 10012 and 10051 (and subsequently save the new file with fwrite). Do you think this is possible?
Cheers,
Jelmer
________________________________________ Van: freesurfer-bounces@nmr.mgh.harvard.edu [freesurfer-bounces@nmr.mgh.harvard.edu] namens Bruce Fischl [fischl@nmr.mgh.harvard.edu] Verzonden: donderdag 25 februari 2016 20:32 Aan: Freesurfer support list Onderwerp: Re: [Freesurfer] Minor adjustment to RB_all_2008-03-26.gca
Hi Jelmer
if you have manually split asegs you can use them as input to mri_ca_train to create a new gca that has the separated putamen. You'll need to add a couple of entries to FreeSurferColorLUT.txt to specify how to display and name the new entries as well
cheers Bruce
On Thu, 25 Feb 2016, Kok, JG (neuro) wrote:
Hi,
I wish to split the putamen in an anterior part and a posterior part in >60 subjects which have been analyzed by the recon-all stream (FS version 5.3.0). Ideally, I would split the putamen in the atlas (e.g., at the coronal plane through the anterior commissure) and apply this "new" atlas to all subjects, one of the advantages being the bias that would be introduced by splitting the putamena manually in all subjects would be avoided. After some time trying to make this happen, I now think this could be achieved by adjusting RB_all_2008-03-26.gca and applying something like
"mri_ca_label -align subject_norm.mgz subject_talairach.m3z RB_all_2008-03-26.gca subject_aseg.auto_noCCseg.PUTAMENA_SEPARATED.mgz".
The problem is that till now I wasn't able to figure out how to correctly read the binary file RB_all_2008-03-26.gca. I tried to use the MATLAB function load_gca.m, but that seems to only work up to version 4 of the .gca file format (as GCA_VERSION is set to 4 in the function), and RB_all_2008-03-26.gca seems to be version 5. Also, when trying to read the file anyway in MATLAB using the code written for version 4 and below (by setting GCA_VERSION to 5.0), the code fails (the first gibbs_nlabels is already way too large, although I suspect the error occurs earlier).
My questions:
- Do you think I am following the right track to achieve a sensible
splitting of the putamena this way? Or should I try something else? 2. Is there a newer load_gca.m file around which is able to read a version 5 .gca file? Or could you point me to the difference in the design of the version 5 .gca file compared to that of the version 4 file so that I can rewrite load_gca.m?
Cheers,
Jelmer
De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
_______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.
________________________________ De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
Hi Jelmer
no, that won't work as the gca contains more information than that. YOu will need to use mri_ca_train to create a new gca, which should be relatively straightforward
cheers Bruce On Thu, 25 Feb 2016, Kok, JG (neuro) wrote:
Hi Bruce,
Thanks for your quick reply. I do not have manually split asegs, I wish to create them all automatically. I'm hoping it is possible to read RB_all_2008-03-26.gca (preferably in MATLAB, using fread), and adjust all putamen values posterior from the coronal plane through the anterior commissure from (FreeSurferColorLUT)values 12 and 51 to for example values 10012 and 10051 (and subsequently save the new file with fwrite). Do you think this is possible?
Cheers,
Jelmer
Van: freesurfer-bounces@nmr.mgh.harvard.edu [freesurfer-bounces@nmr.mgh.harvard.edu] namens Bruce Fischl [fischl@nmr.mgh.harvard.edu] Verzonden: donderdag 25 februari 2016 20:32 Aan: Freesurfer support list Onderwerp: Re: [Freesurfer] Minor adjustment to RB_all_2008-03-26.gca
Hi Jelmer
if you have manually split asegs you can use them as input to mri_ca_train to create a new gca that has the separated putamen. You'll need to add a couple of entries to FreeSurferColorLUT.txt to specify how to display and name the new entries as well
cheers Bruce
On Thu, 25 Feb 2016, Kok, JG (neuro) wrote:
Hi,
I wish to split the putamen in an anterior part and a posterior part in >60 subjects which have been analyzed by the recon-all stream (FS version 5.3.0). Ideally, I would split the putamen in the atlas (e.g., at the coronal plane through the anterior commissure) and apply this "new" atlas to all subjects, one of the advantages being the bias that would be introduced by splitting the putamena manually in all subjects would be avoided. After some time trying to make this happen, I now think this could be achieved by adjusting RB_all_2008-03-26.gca and applying something like
"mri_ca_label -align subject_norm.mgz subject_talairach.m3z RB_all_2008-03-26.gca subject_aseg.auto_noCCseg.PUTAMENA_SEPARATED.mgz".
The problem is that till now I wasn't able to figure out how to correctly read the binary file RB_all_2008-03-26.gca. I tried to use the MATLAB function load_gca.m, but that seems to only work up to version 4 of the .gca file format (as GCA_VERSION is set to 4 in the function), and RB_all_2008-03-26.gca seems to be version 5. Also, when trying to read the file anyway in MATLAB using the code written for version 4 and below (by setting GCA_VERSION to 5.0), the code fails (the first gibbs_nlabels is already way too large, although I suspect the error occurs earlier).
My questions:
- Do you think I am following the right track to achieve a sensible
splitting of the putamena this way? Or should I try something else? 2. Is there a newer load_gca.m file around which is able to read a version 5 .gca file? Or could you point me to the difference in the design of the version 5 .gca file compared to that of the version 4 file so that I can rewrite load_gca.m?
Cheers,
Jelmer
De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.
De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
Hi Bruce,
I'll have a look at mri_ca_train. Thanks again!
Cheers,
Jelmer ________________________________________ Van: freesurfer-bounces@nmr.mgh.harvard.edu [freesurfer-bounces@nmr.mgh.harvard.edu] namens Bruce Fischl [fischl@nmr.mgh.harvard.edu] Verzonden: donderdag 25 februari 2016 20:58 Aan: Freesurfer support list Onderwerp: Re: [Freesurfer] Minor adjustment to RB_all_2008-03-26.gca
Hi Jelmer
no, that won't work as the gca contains more information than that. YOu will need to use mri_ca_train to create a new gca, which should be relatively straightforward
cheers Bruce On Thu, 25 Feb 2016, Kok, JG (neuro) wrote:
Hi Bruce,
Thanks for your quick reply. I do not have manually split asegs, I wish to create them all automatically. I'm hoping it is possible to read RB_all_2008-03-26.gca (preferably in MATLAB, using fread), and adjust all putamen values posterior from the coronal plane through the anterior commissure from (FreeSurferColorLUT)values 12 and 51 to for example values 10012 and 10051 (and subsequently save the new file with fwrite). Do you think this is possible?
Cheers,
Jelmer
Van: freesurfer-bounces@nmr.mgh.harvard.edu [freesurfer-bounces@nmr.mgh.harvard.edu] namens Bruce Fischl [fischl@nmr.mgh.harvard.edu] Verzonden: donderdag 25 februari 2016 20:32 Aan: Freesurfer support list Onderwerp: Re: [Freesurfer] Minor adjustment to RB_all_2008-03-26.gca
Hi Jelmer
if you have manually split asegs you can use them as input to mri_ca_train to create a new gca that has the separated putamen. You'll need to add a couple of entries to FreeSurferColorLUT.txt to specify how to display and name the new entries as well
cheers Bruce
On Thu, 25 Feb 2016, Kok, JG (neuro) wrote:
Hi,
I wish to split the putamen in an anterior part and a posterior part in >60 subjects which have been analyzed by the recon-all stream (FS version 5.3.0). Ideally, I would split the putamen in the atlas (e.g., at the coronal plane through the anterior commissure) and apply this "new" atlas to all subjects, one of the advantages being the bias that would be introduced by splitting the putamena manually in all subjects would be avoided. After some time trying to make this happen, I now think this could be achieved by adjusting RB_all_2008-03-26.gca and applying something like
"mri_ca_label -align subject_norm.mgz subject_talairach.m3z RB_all_2008-03-26.gca subject_aseg.auto_noCCseg.PUTAMENA_SEPARATED.mgz".
The problem is that till now I wasn't able to figure out how to correctly read the binary file RB_all_2008-03-26.gca. I tried to use the MATLAB function load_gca.m, but that seems to only work up to version 4 of the .gca file format (as GCA_VERSION is set to 4 in the function), and RB_all_2008-03-26.gca seems to be version 5. Also, when trying to read the file anyway in MATLAB using the code written for version 4 and below (by setting GCA_VERSION to 5.0), the code fails (the first gibbs_nlabels is already way too large, although I suspect the error occurs earlier).
My questions:
- Do you think I am following the right track to achieve a sensible
splitting of the putamena this way? Or should I try something else? 2. Is there a newer load_gca.m file around which is able to read a version 5 .gca file? Or could you point me to the difference in the design of the version 5 .gca file compared to that of the version 4 file so that I can rewrite load_gca.m?
Cheers,
Jelmer
De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
The information in this e-mail is intended only for the person to whom it is addressed. If you believe this e-mail was sent to you in error and the e-mail contains patient information, please contact the Partners Compliance HelpLine at http://www.partners.org/complianceline . If the e-mail was sent to you in error but does not contain patient information, please contact the sender and properly dispose of the e-mail.
De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
_______________________________________________ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer ________________________________ De inhoud van dit bericht is vertrouwelijk en alleen bestemd voor de geadresseerde(n). Anderen dan de geadresseerde(n) mogen geen gebruik maken van dit bericht, het niet openbaar maken of op enige wijze verspreiden of vermenigvuldigen. Het UMCG kan niet aansprakelijk gesteld worden voor een incomplete aankomst of vertraging van dit verzonden bericht.
The contents of this message are confidential and only intended for the eyes of the addressee(s). Others than the addressee(s) are not allowed to use this message, to make it public or to distribute or multiply this message in any way. The UMCG cannot be held responsible for incomplete reception or delay of this transferred message.
freesurfer@nmr.mgh.harvard.edu
-
Bruce Fischl -
Kok, JG (neuro)