Re: [Freesurfer] recon-all bug
External Email - Use Caution Hi freesurfer team, I've tried running recon-all as mentioned before with the -xf flag under different conditions. As far as I can tell, any line using the UpdateNeeded function. e.g.: "set ud = `UpdateNeeded $thickness $pial $white`" fails with the following if statement because the variable $ud includes paths as strings in it instead of a 0/1 integer. Setting it to a 1 or 0 beforehand also seems to solve the problem. The function UpdateNeeded only seems to return an integer as an exit status, so I'm not sure why there are path strings in it causing it to fail. If I print out the $ud variable before it fails and the result is: "/subjectDirectory/subjectID_FREESURFER/mri mri mri_FREESURFER /subjectDirectory/subjectID_FREESURFER/mri/mri_FREESURFER/stats 1" It returns the exit status of 1 preceded by the path of the subject's mri directory, just "mri", "mri_FREESURFER" and the subject's stats directory (all subjects in this particular study have the suffix "_FREESURFER"). It prints this out these strings no matter the variables as .mgz file names, or number of variables passed to the UpdateNeeded function. I still don't have a clue what would cause this behavior. The UpdateNeeded function seems to work fine and as intended as far as I can tell, so I still don't know why these paths/strings are being added and causing it to crash. Best, Kody Zalewski Message: 3 Date: Fri, 02 Apr 2021 03:49:02 +0200 From: fsbuild Subject: Re: [Freesurfer] recon-all bug To: freesurfer@nmr.mgh.harvard.edu Cc: cm...@neurosurgery.washington.edu Message-ID: <1617328142.6066780ef2...@trashmail.com> Content-Type: text/plain; charset="utf-8" External Email - Use Caution Hello Kody, One thing you could try to get more information printed about what is going on is to edit the very first line in the recon-all script to be: #! /bin/tcsh -xf ? instead of #! /bin/tcsh -f Adding the x option will echo each line in the script as it?s run with the variable values expanded, and seeing that output could help with debugging. - R. On Apr 1, 2021, at 20:29, zalewk zal...@uw.edu wrote:External Email - Use CautionHi Freesurfer team,My lab recently switched over from FS6 to FS7 and I think I've encountered a bug with recon-all..I'm running subjects stepwise and editing them, running recon-all with the -autorecon2 and -autorecon3 flags initially are no problem, they run fine. Trying to re-run them or use the -autorecon2-wm, -autorecon-pial etc. flags lead to crashing at certain parts of the program. The error message I get is "if: expression syntax". Digging into the code, it seems to crash when it hits a line of code "if($ud || $ForceUpdate)". The syntax superficially looks right (I don't know tcsh), checking the variables it seems like $ForceUpdate is never set, but removing it from recon-all doesn't fix the problem.It seems the $ud variable is for checking if updates are ! needed fo r a set of volumes passed to the program. I don't know why that would be a problematic argument, removing that part seems to fix it, but this also seems like an important variable and I don't know enough about tcsh to understand why it works after removing it. Can one of the folks who wrote the recon-all script help me out? I'm running the commands step-wise as a work-around, but that is much less convenient.Best,Kody Zalewski___Freesurfer mailing listfreesur...@nmr.mgh.harvard.eduhttps://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
[Freesurfer] recon-all bug
External Email - Use Caution Hi Freesurfer team, My lab recently switched over from FS6 to FS7 and I think I've encountered a bug with recon-all. I'm running subjects stepwise and editing them, running recon-all with the -autorecon2 and -autorecon3 flags initially are no problem, they run fine. Trying to re-run them or use the -autorecon2-wm, -autorecon-pial etc. flags lead to crashing at certain parts of the program. The error message I get is "if: expression syntax". Digging into the code, it seems to crash when it hits a line of code "if($ud || $ForceUpdate)". The syntax superficially looks right (I don't know tcsh), checking the variables it seems like $ForceUpdate is never set, but removing it from recon-all doesn't fix the problem. It seems the $ud variable is for checking if updates are needed for a set of volumes passed to the program. I don't know why that would be a problematic argument, removing that part seems to fix it, but this also seems like an important variable and I don't know enough about tcsh to understand why it works after removing it. Can one of the folks who wrote the recon-all script help me out? I'm running the commands step-wise as a work-around, but that is much less convenient. Best, Kody Zalewski ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
[no subject]
...@nmr.mgh.harvard.edu X-MXTHUNDER-Identifier: X-MXTHUNDER-IP-Rating: 0, 155.52.251.87, Ugly c=0 p=0 Source New X-MXTHUNDER-Scan-Result: 100 X-MXTHUNDER-Rules: 100-75201-11157-11176-m 100-75201-0-26847-f X-MXTHUNDER-Group: Bulk Mail --===1012566416093711800== Content-Language: en-US Content-Type: multipart/alternative; boundary="_000_CY4PR08MB3639C3F96BCDD4D07831A789AC400CY4PR08MB3639namp_" --_000_CY4PR08MB3639C3F96BCDD4D07831A789AC400CY4PR08MB3639namp_ Content-Type: text/plain; charset="iso-8859-1" MIME-Version: 1.0 Content-Transfer-Encoding: quoted-printable External Email - Use Caution= Thanks for getting back to me Doug/Bruce. Sounds good, yes I definitely would want to wait until the bugs have been w= orked out before utilizing latest release. If your dev team needs contributors to document/write the ex-vivo processin= g stream I'd be happy to step in and help, I've been running various ex-viv= o scan acquisitions through freesurfer for a year at this point. I also hav= e experience in C and shell scripting if that helps, as it seems that's wha= t much of freesurfer is written in. Not sure who I would talk to about that. Best, Kody Zalewski Message: 11 Date: Tue, 19 Mar 2019 14:53:30 + From: "Greve, Douglas N.,Ph.D." Subject: Re: [Freesurfer] Omitting aseg regions from segmentationTo: "frees= ur...@nmr.mgh.harvard.edu" Message-ID: <235d9f13-1eda-1d76-920f-d30b9c467...@mgh.harvard.edu> Content-Type: text/plain; charset=3D"windows-1252" ... but know that the dev version is not fully tested -- there may be signi= ficant bugs that will not be uncovered until it is tested. As such, we may = not be able to support any problems you run into. On 3/18/19 8:12 PM, Bruce Fischl wrote: Hi Kody there definitely are flags like this in the dev version. e.g.: case "-no-cerebellum": case "-nocerebellum": case "-NO-CEREBELLUM": case "-NOCEREBELLUM": set nocerebellum=3D( -nocerebellum) breaksw ; Maybe you should just grab the dev version? And sorry, we haven't documented the ex vivo processing stream yet as it is= too much of a work-in-progress. cheers Bruce On Thu, 14 Mar 2019, zalewk wrote: External Email - Use Caution Hi Bruce, In response to your previous email, in v6.0 it doesn't look like recon-all = has a "-no-cerebellum" flag (or any cerebellum-related flags), I did verify that the recon-all scr= ipt in the bin folder had hemisphere-specific flags available however. Is there any documentation on this? It still seems as if I'm missing someth= ing here, was it maybe included in previous versions of freesurfer and removed? Best, Kody Zalewski Date: Mon, 4 Mar 2019 15:07:26 -0500 (EST) From: Bruce Fischl <mailto:fis...@nmr.mgh.harva= rd.edu> Subject: Re: [Freesurfer] Omitting aseg regions from segmentation To: Freesurfer support list <mailto:freesur= f...@nmr.mgh.harvard.edu> Message-ID: Content-Type: text/plain; charset=3D"iso-8859-1" Hi Kody yes, there are some hidden options to recon-all for this. Take a look at the script for things like -no-cerebellum, -rh-only, and see if they are sufficient for what you are trying to do cheers Bruce --_000_CY4PR08MB3639C3F96BCDD4D07831A789AC400CY4PR08MB3639namp_ Content-Type: text/html; charset="iso-8859-1" MIME-Version: 1.0 Content-Transfer-Encoding: quoted-printable <!-- P {margin-top:0;margi= n-bottom:0;} --> External Email -= Use Caution<= /p> Thanks for getting back to me Doug/Bruce. Sounds good, yes I definitely would want to wait until the bugs have been w= orked out before utilizing latest release. If your dev team needs contributors to document/write the ex-vivo processin= g stream I'd be happy to step in and help, I've been running various ex-viv= o scan acquisitions through freesurfer for a year at this point. I also hav= e experience in C and shell scripting if that helps, as it seems that's what much of freesurfer is written in. <= br> Not sure who I would talk to about that. Best, Kody Zalewski Message: 11 Date: Tue, 19 Mar 2019 14:53:30 From: Greve, Douglas N.,Ph.D. dgr...@mgh.harvard.edu Subject: Re: [Freesurfer] Omitting aseg regions from segmentationTo: = freesurfer@nmr.mgh.harvard.edu freesurfer@nmr.mgh.harvard.edu= Message-ID: 235d9f13-1eda-1d76-920f-d30b9c467...@mgh.harvard.edu Content-Type: text/plain; charset=3Dwindows-1252 ... but know that the dev version is not fully tested -- there may be signi= ficant bugs that will not be uncovered until it is tested. As such, we may = not be able to support any problems you run into. On 3/18/19 8:12 PM, Br= uce Fischl wrote: Hi Kody there definitely are flags like this in the dev version. e.g.: case -no-cerebellum: case -nocerebel
[no subject]
, __HAS_HTML 0, __HAS_MSGID 0, __HAS_XOIP 0, __HTML_TAG_DIV 0, __MIME_HTML 0, __MIME_TEXT_H 0, __MIME_TEXT_H1 0, __MIME_TEXT_H2 0, __MIME_TEXT_P 0, __MIME_TEXT_P1 0, __MIME_TEXT_P2 0, __MIME_VERSION 0, __RDNS_WEBMAIL 0, __SANE_MSGID 0, __STOCK_PHRASE_24 0, __STYLE_RATWARE_NEG 0, __STYLE_TAG 0, __TAG_EXISTS_HTML 0, __TO_MALFORMED_2 0, __TO_NAME 0, __TO_NO_NAME 0, __URI_NO_WWW 0, __URI_NS ' X-NMR-MGH-MailScanner: Not scanned: please contact your Internet E-Mail Service Provider for details, Not scanned: please contact your Internet E-Mail Service Provider for details X-NMR-MGH-MailScanner-SpamCheck: not spam, SpamAssassin (not cached, score=3.066, required 5, autolearn=disabled, HTML_MESSAGE 0.00, RP_8BIT 3.06), not spam, SpamAssassin (cached, score=2.067, required 5, autolearn=disabled, BAD_ENC_HEADER 0.00, HTML_MESSAGE 0.00, MAILING_LIST_MULTI -1.00, RP_8BIT 3.06) X-NMR-MGH-MailScanner-SpamScore: sss, ss X-Spam-Status: No, No Subject: [Freesurfer] Omitting aseg regions from segmentation X-BeenThere: freesurfer@nmr.mgh.harvard.edu X-Mailman-Version: 2.1.9 Precedence: list List-Id: Freesurfer support list List-Unsubscribe: <https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer>, <mailto:freesurfer-requ...@nmr.mgh.harvard.edu?subject=unsubscribe> List-Archive: <http://mail.nmr.mgh.harvard.edu/pipermail/freesurfer> List-Post: <mailto:freesurfer@nmr.mgh.harvard.edu> List-Help: <mailto:freesurfer-requ...@nmr.mgh.harvard.edu?subject=help> List-Subscribe: <https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer>, <mailto:freesurfer-requ...@nmr.mgh.harvard.edu?subject=subscribe> Reply-To: Freesurfer support list Content-Type: multipart/mixed; boundary="===9188205466225368560==" Sender: freesurfer-boun...@nmr.mgh.harvard.edu Errors-To: freesurfer-boun...@nmr.mgh.harvard.edu X-NMR-MGH-MailScanner-Information: Please contact the ISP for more information X-NMR-MGH-MailScanner-ID: x2EI71sN018720 X-NMR-MGH-MailScanner-From: freesurfer-boun...@nmr.mgh.harvard.edu X-MXTHUNDER-Identifier: X-MXTHUNDER-IP-Rating: 0, 155.52.251.20, Ugly c=0 p=0 Source New X-MXTHUNDER-Scan-Result: 100 X-MXTHUNDER-Rules: 100-75201-12926-12945-m 100-75201-0-25630-f X-MXTHUNDER-Group: Bulk Mail --===9188205466225368560== Content-Language: en-US Content-Type: multipart/alternative; boundary="_000_BN6PR08MB3635363BE6C81EC56F06F8B5AC4B0BN6PR08MB3635namp_" --_000_BN6PR08MB3635363BE6C81EC56F06F8B5AC4B0BN6PR08MB3635namp_ Content-Type: text/plain; charset="iso-8859-1" MIME-Version: 1.0 Content-Transfer-Encoding: quoted-printable External Email - Use Caution= Hi Bruce, In response to your previous email, in v6.0 it doesn't look like recon-all = has a "-no-cerebellum" flag (or any cerebellum-related flags), I did verify= that the recon-all script in the bin folder had hemisphere-specific flags = available however. Is there any documentation on this? It still seems as if I'm missing someth= ing here, was it maybe included in previous versions of freesurfer and remo= ved? Best, Kody Zalewski Date: Mon, 4 Mar 2019 15:07:26 -0500 (EST) From: Bruce Fischl Subject: Re: [Freesurfer] Omitting aseg regions from segmentation To: Freesurfer support list Message-ID: Content-Type: text/plain; charset=3D"iso-8859-1" Hi Kody yes, there are some hidden options to recon-all for this. Take a look at the script for things like -no-cerebellum, -rh-only, and see if they are sufficient for what you are trying to do cheers Bruce On Mon, 4 Mar 2019, zalewk wrote: > > External Email - Use Caution > > Hello Freesurfer team, > > I am currently working processing the volumetrics of exvivo brains, many = of which do not have a > brainstem/cerebellum. > It seems this often results in a subcortical alignment that isn't optimal= , and fails to create the > pial/wm surface where freesurfer attempts to label more ventral non-zero = voxels as > brainstem/cerebellum. Currently, I have to manually edit the aseg.auto_no= CCseg.label_intensities.mgz > and aseg.presurf.mgz to fix this. > > Is there a way to omit aseg labels (like the brainstem/cerebellum) from t= he processing stream that > I've overlooked? If not, would there be a potential alternate solution su= ch as creating a second > atlas from the original sans these ROIs and have the cross-sectional stre= am still run normally? > > Best, > Kody Zalewski --_000_BN6PR08MB3635363BE6C81EC56F06F8B5AC4B0BN6PR08MB3635namp_ Content-Type: text/html; charset="iso-8859-1" MIME-Version: 1.0 Content-Transfer-Encoding: quoted-printable <!-- P {margin-top:0;margi= n-bottom:0;} --> External
[no subject]
x-uw-ex-currentlocation-365-recipient: Outside x-microsoft-exchange-diagnostics: =?iso-8859-1?Q?1; CY4PR08MB3573; 23:4qVRtf+4Frk8+HkoSkWpxiBI3Aktx4BluHS3p5r?= =?iso-8859-1?Q?Y9pwKqXqpNNthdH0OWUp6OBtvK1l/YV1ElGbh7XgLoXOL3hSC0itzucGip?= =?iso-8859-1?Q?p+wDHeAm0zUT5UWpPv2w4rFLyg8NfIEjU1A1cygRV2xQ2Ol7hLoOXvy6rA?= =?iso-8859-1?Q?sEA4Im0UhI8Csa6NEdNTAP22KGtPL4hETtiiLnSaW5DGKCTGSREjkNCtJ0?= =?iso-8859-1?Q?18HvvK7POiaz+vnUIGaOs4C/w5UbtGOE4dQWIw9OYP3fRdrPqtfumZq1b+?= =?iso-8859-1?Q?Opp8jlxZic1wbdhkb0txLLINMtuYn4Y/Sn0e6LXV2TwD7RLK8g0TxUV1UW?= =?iso-8859-1?Q?qSCEhsIQ3odu7DeSUhTdHKnPOxrMJeIadLRS95PGAi70j8BoWFXzo/Yl35?= =?iso-8859-1?Q?fCN0UmqG4acCz+xJ20q/pYJY8iErX1bBICC9ot5KxfW5rgrFsi+uQa3HLM?= =?iso-8859-1?Q?8thwBBIJa1xzarXaTZKYHfQXwAtUS9kx6wR9lfEtg7odAUJAy0yh5Os1Ge?= =?iso-8859-1?Q?YOahMdFZmPsdsThbxs2NZ1S3OzPDu7yjUmH3exPG8+ga4kaH1gssZBHp9Z?= =?iso-8859-1?Q?Z6H91mdMF/snNJb0INHdo5nbjMWZ3ntJK6F5VOJYbipOfSGsi39kHP8LuV?= =?iso-8859-1?Q?1GRw7l2c9k+Llh3nwsjYo7CxDNjCo2VsF6DW7Wh3KN/Ve1IF8i3mIrMtuh?= =?iso-8859-1?Q?uj9hNEpuASB8YQibnAd3mftdR9SNE2eYMI1kZe3bfPUyvqpkYJ/nlxt3Fx?= =?iso-8859-1?Q?nMbJ8fkFkV34J5dnB3GVmNn0OzphlV8San71aDYcuzoHR2clETw9iDWU0/?= =?iso-8859-1?Q?0LlM4hPv6PGQgVxpF2ubbbXq97mKY8DuuRMmWaUi0T9NdrzB7/gfaBr69K?= =?iso-8859-1?Q?vdMV7t5mtQlglKZb/ft/qXbcRLeWUXy0oPqGPMVIFWyMVfWL2eHWQnvw+1?= =?iso-8859-1?Q?o+d3KTkrdczyNSJNSao18xPxSW0fJRbgJEsgqfkg/WeMu5AA/D3oXC3P8W?= =?iso-8859-1?Q?3yYTs2+FxP9tkQ3v7d7K6OyAKfDrPWlek/9IiQT8o1r4Xc6Ip/I9zgakTg?= =?iso-8859-1?Q?FjfLxSevYwTaHU8peuSQf17Wfjyeh8xAhIWQPn2LQAueob3hTOT1uz5K2C?= =?iso-8859-1?Q?G5/1Op7C6A04QS0hS/f42Zby8j4T74im2ZVcoxm5Q6fs0JlQTQZDE1wQCy?= =?iso-8859-1?Q?a1Xu1nJ8vdEH3Ome1KCWpgVvu4WEpDNydY8OwLJVsgBZvLiJ7H+XpzAtAP?= =?iso-8859-1?Q?ThcoZLUxN/cQMM1QM+eKqd/kUn0POyu9qP9rMplYhHP2NAgW+pfwZdfVBq?= =?iso-8859-1?Q?IVRBK2xT7WX1bGOVTh3P2VaA9BmPySMzb2eqpVujxvRsIwTgfsMvvi7dLo?= =?iso-8859-1?Q?spS0Q3bPaaxX1bkoAgCL/6Vgfjp4jIh3HuKoB7JpGP/LBoYBn8w=3D=3D?= x-microsoft-antispam-prvs: x-forefront-prvs: 09669DB681 x-forefront-antispam-report: SFV:NSPM; SFS:(10009020)(396003)(3986042)(136003)(366004)(376002)(346002)(199004)(189003)(564073)(97736004)(99286004)(6916009)(6606003)(2906002)(7120041)(7119041)(19627405001)(106356001)(25786009)(105586002)(6436002)(75432002)(2501003)(52536013)(476003)(33656002)(8936002)(47861)(486006)(86362001)(14454004)(4744005)(566032)(68736007)(186003)(26005)(102836004)(53936002)(6116002)(3846002)(6506007)(7736002)(8676002)(74316002)(81156014)(2171002)(55016002)(66066001)(173073)(81166006)(9686003)(54896002)(86152003)(316002)(786003)(88552002)(2351001)(256004)(7696005)(53336002); DIR:OUT; SFP:1101; SCL:1; SRVR:CY4PR08MB3573; H:CY4PR08MB3639.namprd08.prod.outlook.com; FPR:; SPF:None; LANG:en; PTR:InfoNoRecords; A:1; MX:1; received-spf: None (protection.outlook.com: uw.edu does not designate permitted sender hosts) x-ms-exchange-senderadcheck: 1 x-microsoft-antispam-message-info: jCnFTSMixpgSZQAG3fC7hiywhgemVoPP2h+IKQ0jmXpLtKqZPjh1GVz2c7wA7JKAqb8f+uEMt5YeFDfbHKEDGyfN028BRmG2boXXrfSZLcLuTyB7HuCNZJSOtSKEJJqc8y+oEw8M9SbTKT1YsZ8tNgb1E1saA8Ge42lLed4rVoC3WUgo6u3BzApL88loE2eWPXoj8T9ECO/+YU73toZdFboWCOF8MiUaHWjsxaVMFSF+t9j0corZ4IwcwC1hBFjAtMShmmFZi/EekTixmOVRSEO4ZEsPM3182eLJr6I7dlZQyfskJ67pTvf1lIwndvleUkQUlNZ4urMVQkNp5F1P5Zdr7QsYk2hR7+G+ggSOwkPi7q8wKlwRqnZX8j0T4ysUP8lK7dtwMx9/5uMs/wvc/RIFN4Cepr4H4vNITS9A5xM= MIME-Version: 1.0 X-MS-Exchange-CrossTenant-Network-Message-Id: 6f79dd33-9f1d-4afd-afa1-08d6a0ce85d5 X-MS-Exchange-CrossTenant-originalarrivaltime: 04 Mar 2019 18:23:38.7744 (UTC) X-MS-Exchange-CrossTenant-fromentityheader: Hosted X-MS-Exchange-CrossTenant-id: f6b6dd5b-f02f-441a-99a0-162ac5060bd2 X-MS-Exchange-CrossTenant-mailboxtype: HOSTED X-MS-Exchange-Transport-CrossTenantHeadersStamped: CY4PR08MB3573 X-OriginatorOrg: uw.edu X-PMX-Version: 6.4.6.2792898, Antispam-Engine: 2.7.2.2107409, Antispam-Data: 2019.3.4.181216, AntiVirus-Engine: 5.58.0, AntiVirus-Data: 2019.3.4.5580001 X-PMX-Server: mxout25.s.uw.edu X-Uwash-Spam: Gauge=, Probability=8%, Report=' HTML_NO_HTTP 0.1, FROM_NAME_ONE_WORD 0.05, SUPERLONG_LINE 0.05, BODYTEXTH_SIZE_1_LESS 0, BODYTEXTP_SIZE_3000_LESS 0, BODY_SIZE_2000_2999 0, BODY_SIZE_5000_LESS 0, BODY_SIZE_7000_LESS 0, DKIM_SIGNATURE 0, FROM_EDU_TLD 0, NO_URI_HTTPS 0, SPF_NEUTRAL 0, WEBMAIL_SOURCE 0, WEBMAIL_XOIP 0, WEBMAIL_X_IP_HDR 0, __ANY_URI 0, __CT 0, __CTYPE_HAS_BOUNDARY 0, __CTYPE_MULTIPART 0, __CTYPE_MULTIPART_ALT 0, __DQ_NEG_HEUR 0, __DQ_NEG_IP 0, __FUR_RDNS_OUTLOOK 0, __HAS_FROM 0, __HAS_HTML 0, __HAS_MSGID 0, __HAS_XOIP 0, __HTML_TAG_DIV 0, __MIME_HTML 0, __MIME_TEXT_H 0, __MIME_TEXT_H1 0, __MIME_TEXT_H2 0, __MIME_TEXT_P 0, __MIME_TEXT_P1 0,
[no subject]
n-Result: 100 X-MXTHUNDER-Rules: 100-75201-11621-11640-m 100-75201-0-19892-f X-MXTHUNDER-Group: Bulk Mail --===3208088125263462719== Content-Language: en-US Content-Type: multipart/alternative; boundary="_000_CY4PR08MB36392D74B065DC6A604AC7A6AC710CY4PR08MB3639namp_" --_000_CY4PR08MB36392D74B065DC6A604AC7A6AC710CY4PR08MB3639namp_ Content-Type: text/plain; charset="iso-8859-1" MIME-Version: 1.0 Content-Transfer-Encoding: quoted-printable External Email - Use Caution= Hello Freesurfer Team, Sorry to bother you guys again, any word on what you think is the best appr= oach going forward here? Best, Kody Zalewski From: zalewk Sent: Thursday, June 7, 2018 10:39 AM To: freesurfer@nmr.mgh.harvard.edu Subject: Longitudinal analysis of multiple scanner acquisitions Good afternoon Freesurfer team, I want to run the longitudinal pipeline on multiple groups of subjects, the= problem is we are at the mercy of the scanners available to us at the vari= ous time-points. We have two different Siemens 1.5 scanners (same type of s= canner, there was an effort made to match them) at time point 1 and 2 with = a Phillips 3T at time point 3. Is it worth the time to attempt a longitudin= al analysis on this cohort? Best, Kody Zalewski --_000_CY4PR08MB36392D74B065DC6A604AC7A6AC710CY4PR08MB3639namp_ Content-Type: text/html; charset="iso-8859-1" MIME-Version: 1.0 Content-Transfer-Encoding: quoted-printable <!-- P {margin-top:0;margi= n-bottom:0;} --> External Email -= Use Caution<= /p> Hello Freesurfer Team, Sorry to bother you guys again, any word on what you think is the best appr= oach going forward here? Best, Kody Zalewski From: zalewk Sent: Thursday, June 7, 2018 10:39 AM To: freesurfer@nmr.mgh.harvard.edu Subject: Longitudinal analysis of multiple scanner acquisitions Good afternoon Freesurfer team,<= br> I want to run the longitudinal pipeline on multiple groups of subjects, the= problem is we are at the mercy of the scanners available to us at the vari= ous time-points. We have two different Siemens 1.5 scanners (same type of s= canner, there was an effort made to match them) at time point 1 and 2 with a Phillips 3T at time point 3. I= s it worth the time to attempt a longitudinal analysis on this cohort? Best, Kody Zalewski --_000_CY4PR08MB36392D74B065DC6A604AC7A6AC710CY4PR08MB3639namp_-- --===3208088125263462719== Content-Type: text/plain; charset="us-ascii" MIME-Version: 1.0 Content-Disposition: inline ___ 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. --===3208088125263462719==--
[Freesurfer] Longitudinal analysis of multiple scanner acquisitions
External Email - Use Caution Good afternoon Freesurfer team, I want to run the longitudinal pipeline on multiple groups of subjects, the problem is we are at the mercy of the scanners available to us at the various time-points. We have two different Siemens 1.5 scanners (same type of scanner, there was an effort made to match them) at time point 1 and 2 with a Phillips 3T at time point 3. Is it worth the time to attempt a longitudinal analysis on this cohort? Best, Kody Zalewski ___ 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.
Re: [Freesurfer] clarification
My apologies as I misstated some of the data in my first email, The voxel sizes are relatively similar, (two of the acquisitions are 1mm vs one with 0.9mm). Whereas the FOV is different between the scans being 256x256x176 on one acquisition, 224x256x176 on the second and 240x240x180 on the third. KZ On Tue, 5 Jul 2016, zal...@u.washington.edu wrote: > Hello Freesurfer team, > > I have the unfortunate task of performing longitudinal volumetric analysis on > subjects collected on three different scanners with different voxel > numbers/sizes. I didn't find this out until attempting to run the > longitudinal base and receiving this message: > > WARNING: Image geometries differ across time, maybe due to aquisition changes? > This can potentially bias a longitudinal study! Will continue in 10s. > > By re-registering the original nifti and recreating the norm.mgz/talairach I > can correct the geometry, but FSL's Flirt will only work with niftis and not > mgz volumes, ergo I can't re-register the manually edited brainmask/wm > volumes properly. > I was wondering if there were any freesurfer tools I could use to simply > re-register the .mgz volumes to the new norm.mgz? If so, how would I go about > doing that? Seeing as the brainmask/wm are both registered to 256^3/1mm > anyhow, is this something I should even be concerned about? > > Thanks, > KZ > > > ___ 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.
[Freesurfer] Registering scans with different voxel sizes/numbers
Hello Freesurfer team, I have the unfortunate task of performing longitudinal volumetric analysis on subjects collected on three different scanners with different voxel numbers/sizes. I didn't find this out until attempting to run the longitudinal base and receiving this message: WARNING: Image geometries differ across time, maybe due to aquisition changes? This can potentially bias a longitudinal study! Will continue in 10s. By re-registering the original nifti and recreating the norm.mgz/talairach I can correct the geometry, but FSL's Flirt will only work with niftis and not mgz volumes, ergo I can't re-register the manually edited brainmask/wm volumes properly. I was wondering if there were any freesurfer tools I could use to simply re-register the .mgz volumes to the new norm.mgz? If so, how would I go about doing that? Seeing as the brainmask/wm are both registered to 256^3/1mm anyhow, is this something I should even be concerned about? Thanks, KZ ___ 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.
[Freesurfer] QA measures mgz
Excellent, thank you for sharing! I have one more question. When running the mri_cnr function I tried using a number of different .mgz files with varying results, which would you recommend for the most accurate gm/wm and gm/csf CNR? Best, KZ -- Message: 1 Date: Fri, 17 Jul 2015 18:24:28 +0200 From: Jacek Manko j...@wp.pl Subject: Re: [Freesurfer] QA measures To: Freesurfer support list freesurfer@nmr.mgh.harvard.edu Cc: freesurfer@nmr.mgh.harvard.edu Message-ID: 55a92c3ce71501.52368...@wp.pl Content-Type: text/plain; charset=iso-8859-2 You may find this paper useful: http://www.biomedcentral.com/1471-2342/12/27 Authors describe there some measures such as Euler number and #8220;Contrast to Noise ratio (CNR). Actually, they dont belong to QA package, but they are definitely worth looking at. Best, Jacek Manko Dnia 17-07-2015 o godz. 0:03 zal...@u.washington.edu napisa?(a): Dear Freesurfer Team, We would like to have a Q/A measure (or two) that would give us indication of how well your software worked on each of our data sets. Obviously some image sets will be noisier than others for a variety of reasons. Could you point us towards what options we may have for extracting such a measure from your processing pipeline? We were playing around with QAtools but wanted to understand if there is a systematic computation of say a handful of variables that would give us a global representation of how well things worked for each case we run? Something that we could compare across our patients and even across data sets? Best, KZ ___ 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 -- ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer
[Freesurfer] QA measures
Dear Freesurfer Team, We would like to have a Q/A measure (or two) that would give us indication of how well your software worked on each of our data sets. Obviously some image sets will be noisier than others for a variety of reasons. Could you point us towards what options we may have for extracting such a measure from your processing pipeline? We were playing around with QAtools but wanted to understand if there is a systematic computation of say a handful of variables that would give us a global representation of how well things worked for each case we run? Something that we could compare across our patients and even across data sets? Best, KZ ___ 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.
[Freesurfer] QAtools error
Dear Freesurfer experts, I am attempting to use QAtools to retrieve QA measures for my subjects. I can obtain WM SNR/avg/sd, but it appears the structural segmentation table is from a previous version of freesurfer (as the file order is different) and doesn't appear to work with my subjects processed recently in v5.3. To fix this I would imagine one would need to create a new template off of a subject with the -gen-asegLUT flag or use the -asegLUT flag. Attempting to use -asegLUT flag with the ASegStatsLUT.txt results in: Processing each subcortical label ... mysubject has 35 outliers ... (With all structures as outliers, and without reporting the SNR/avg/sd for each segmentation) When attempting to write out a new segmentation table with the -gen-asegLUT flag the error I receive is: gawk: /usr/local/freesurfer/QAtools/QAtools_v1.1/gstd.awk:22: (FILENAME=- FNR=1) fatal: division by zero attempted (standard_in) 1: syntax error (standard_in) 1: syntax error The command is here: $QA_tools/recon_checker -s mysubject -asegLUT $FREESURFER_HOME/ASegStatsLUT.txt Am I going about this correctly? Best, KZ ___ 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.
[Freesurfer] Editing ROI labels
Dear Freesurfer experts, My PI would like to know, when editing the ROI labels from the freesurfer color look-up table/aparc+aseg.mgz file with freeview voxel editor (at the end of -autorecon3), what commands need to be re-run to apply the changes to reflect the accurate volumetric data of each label? As, if I'm not mistaken, running -autorecon3 would overwrite all of the edits made at this point. Best, Kody J. Zalewski ___ 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.
Re: [Freesurfer] Stats file question
Hi Doug, I sincerely apologize for being unclear! I'm still relatively new to freesurfer and learning the ropes. I've attached the three aforementioned .stats files. The lh.aparc.stats and DKTatlas40.stats were generated from running recon-all. When comparing them it appears the grey matter volume data between the same ROIs can vary significantly. I have also been running the commands individually from the recon-all v5.3 dev table (listed here https://surfer.nmr.mgh.harvard.edu/fswiki/ReconAllDevTable). When I ran the mri_segstats command to generate the wmparc.stats file: mri_segstats --seg mri/wmparc.mgz --sum stats/wmparc.stats --pv mri/norm.mgz --excludeid 0 --brain-vol-from-seg --brainmask mri/brainmask.mgz --in mri/norm.mgz --in-intensity-name norm --in-intensity-units MR --etiv --subject TB0204U_FREESURFER_THIRDRUN --surf-wm-vol --ctab /usr/local/freesurfer/FreeSurferColorLUT.txt The cortical ROIs were listed alongside the white matter segmentations in wmparc.stats, and now have different grey matter volumes than either of the previous two .stats files, which has left me even more confused. Since these cortical volumes all differ from one another, it makes me wonder which .stats file has the most accurate cortical volumetric data. Best, Kody J. Zalewski Message: 4 Date: Wed, 27 May 2015 14:58:44 -0400 From: Douglas N Greve gr...@nmr.mgh.harvard.edu Subject: Re: [Freesurfer] Stats file question To: freesurfer@nmr.mgh.harvard.edu Message-ID: 556613e4.3060...@nmr.mgh.harvard.edu Content-Type: text/plain; charset=windows-1252; format=flowed Can you list the 3 stats files you are talking about? I can't quite follow from below doug On 05/26/2015 02:21 PM, zal...@u.washington.edu wrote: Hello Freesurfer team, When generating the .stats files for grey matter volumetric data, this information appears in multiple places and I'm not sure which data I should be using for my final results. There appears to be cortical information generated in both the .stats file based off the DKT atlas and the l/rh.aparc.stats file. The information in these files is quite different from one another. Using the FreeSurferColorLUT.txt (as listed in the recon-all dev table) instead of WMParcStatsLUT.txt (which is used when running recon-all) in generating the mri_segstats wmparc.stats file also produces cortical volumetric data. This third set of data differs significantly from the DKT Atlas, and differs less-so from the cortical volumes produced in the aparc.stats file. Which of the .stats files containing grey matter volume data is most accurate? Best, Kody J. Zalewski ___ Freesurfer mailing list Freesurfer@nmr.mgh.harvard.edu https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer -- Douglas N. Greve, Ph.D. MGH-NMR Center gr...@nmr.mgh.harvard.edu Phone Number: 617-724-2358 Fax: 617-726-7422 lh.aparc.DKTatlas40.stats Description: Binary data lh.aparc.stats Description: Binary data wmparc.stats Description: Binary data ___ 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.
[Freesurfer] Stats file question
Hello Freesurfer team, When generating the .stats files for grey matter volumetric data, this information appears in multiple places and I'm not sure which data I should be using for my final results. There appears to be cortical information generated in both the .stats file based off the DKT atlas and the l/rh.aparc.stats file. The information in these files is quite different from one another. Using the FreeSurferColorLUT.txt (as listed in the recon-all dev table) instead of WMParcStatsLUT.txt (which is used when running recon-all) in generating the mri_segstats wmparc.stats file also produces cortical volumetric data. This third set of data differs significantly from the DKT Atlas, and differs less-so from the cortical volumes produced in the aparc.stats file. Which of the .stats files containing grey matter volume data is most accurate? Best, Kody J. Zalewski ___ 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.
[Freesurfer] Motion Correction with T2
Hello again Freesurfer team, Is it possible to perform a motion correction with only a T1 and a T2 scan, as my lab only has one available T1 scan per participant. Since the intensities of T1 and T2 are opposite of one another I've had no luck so far, but perhaps your team knows of a work-around. Best, Kody J. Zalewski ___ 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.
[Freesurfer] Cortical volumes/ROI labelling question
Hello Freesurfer team, I've been previously learning freesurfer by going through the individual recon-all commands listed on the freesurfer website v5.3 table here: http://freesurfer.net/fswiki/ReconAllTableStableV5.3 and recently just ran through the recon-all autorecon steps 1, 2 and 3. Running the command individually as listed on the website: mri_segstats --seg mri/wmparc.mgz --sum stats/wmparc.stats --pv mri/norm.mgz --excludeid 0 --brain-vol-from-seg --brainmask mri/brainmask.mgz --in mri/norm.mgz --in-intensity-name norm --in-intensity-units MR --etiv --subject TB0204U_FREESURFER_THIRDRUN --surf-wm-vol --ctab /usr/local/freesurfer/FreeSurferColorLUT.txt This produces both cortical and subcortical volumetric data in the same wmparc.stats file. (attached to this email as wmparc.stats) Whereas, the default command from recon-all is: mri_segstats --seg mri/wmparc.mgz --sum stats/wmparc.stats --pv mri/norm.mgz --excludeid 0 --brainmask mri/brainmask.mgz --in mri/norm.mgz --in-intensity-name norm --in-intensity-units MR --subject TB0204U_FREESURFER_THIRDRUN --surf-wm-vol --ctab /usr/local/freesurfer/WMParcStatsLUT.txt --etiv The resulting file is attached to this email as wmparc2.stats along with the DK40 stats. The only differences are the --brain-vol-from-seg flag (which doesn't seem relevant, based on the description of the flag), and the lookup table used in the generation of the wmparc.stats file. During the recon-all run, the wmparc.stats file only listed subcortical regions, whereas the DK40 atlas is the only source listing cortical volumes. While I would just use the DK40 atlas to retrieve this information, the problem here is that there appears to be a significant difference between the mm^3 volume of hemispherical cortical areas listed in the DK40 atlas compared to those produced when running the command with the flags. Unfortunately, I am not familiar enough with cortical volumes to recognize which estimate is more reasonable one. What .stats file or look-up table should I be using to obtain the most accurate cortical data? Also, when fixing errant ROI label voxels with freeview voxel editor, what commands need to be run post-editing to apply the corrections? I was under the impression that the aparc+aseg.stats file was edited when changing the ROIs. Although, I was rather confused when I went back and regenerated the aparc+aseg.mgz file and my edits were still present. Best, -Kody J. Zalewski wmparc.stats Description: Binary data wmparc2.stats Description: Binary data lh.aparc.DKTatlas40.stats Description: Binary data rh.aparc.DKTatlas40.stats Description: Binary data ___ 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.
[Freesurfer] recon-all compared to dev table
Dear freesurfer team, I'm very new to freesurfer, and I am currently working my way through the recon-all dev table functions (listed here: http://freesurfer.net/fswiki/ReconAllDevTable) step-by-step in preparation for processing a large amount of participant volumetric data. While doing so I've noticed some inconsistencies between the dev table and the series of functions called by the recon-all command. For example: -The apas2aseg function listed in the dev table does not exist in my bin file (nor does it appear to be called by recon-all). -The cortical ribbon output occurs after the three cortical parcellation steps/creation of the anatomical.stats files (instead of prior to these steps like the dev table suggests) -The MNI intensity correction (the mri_nu_correct.mni function) is not called following the skull-stripping. I downloaded the latest version of freesurfer (5.3.0) last week so it should be up-to-date. Thank you for taking the time to answer my questions, -Kody J. Zalewski ___ 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.