External Email - Use Caution        

Hi Lilla,

No worries - any help would be appreciated!

We're collecting low-field (64mT) infant volumes and also processing some
NIHPD2 (1.5T) infant volumes.  We've found that iBEAT currently works a bit
better than Infant Freesurfer at the younger ages (<6 mo), at least for
extracting the pial surface.

From another project, we have .annot files specifying our custom atlas
parcellation in fsaverage space, and our end goal is the extraction of pial
surface area values from atlas ROIs.

To get there from iBEAT pial output (one .vtk triangular mesh per
hemisphere), I think I will need to

   - register the iBEAT pial surfaces to fsaverage (possibly minimizing
   curvature difference as in mris_register)
   - resample to fsaverage resolution(?)
   - use the annot files to identify vertices corresponding to our .annot
   ROIs in the iBEAT surfaces
   - extract surface area from the ROIs in native space

I'm currently working on some of this using the PyVista library in Python,
but I'm better with R than Python and would love to avoid reinventing the
wheel if you know ways to streamline this process within Freesurfer's tools.

Thanks!
Jim H.

On Wed, Apr 2, 2025 at 3:55 PM Zollei, Lilla,Ph.D. <lzol...@mgh.harvard.edu>
wrote:

> Hi Jim,
> I am sorry I missed your message. Are your questions still relevant?
> Apologies, Lilla
> ------------------------------
> *From:* freesurfer-boun...@nmr.mgh.harvard.edu <
> freesurfer-boun...@nmr.mgh.harvard.edu> on behalf of James Hengenius <
> james.hengen...@gmail.com>
> *Sent:* Friday, December 13, 2024 2:49 PM
> *To:* Freesurfer support list <freesurfer@nmr.mgh.harvard.edu>
> *Subject:* [Freesurfer] Importing surface meshes from 3rd party cortical
> surface reconstruction tools for analysis in fsaverage template space
>
>
>         External Email - Use Caution
>
> Hi,
>
> I'm working with infant (0-30 month) MR data, and iBEAT (*MailScanner has
> detected a possible fraud attempt from "secure-web.cisco.com" claiming to
> be* 
> https://secure-web.cisco.com/1wgrR6avzhrKloKIk-lK7SYUO89Nk1mZvx65aN_i4tucrr3ozcqPV2lpqx0RMrrvVLewONwH1Ub4ZNaCSPAfdOZS_EbZphRrMmGyT360MBthjZ_x68P-QlQNZ_vSBeFDOswJ9X1oxnr8tqkPbKhWG-ymp2stCJ4RkMLi3AgiQpdWy_RXy_q4W-2Ce4IMjaM1DFiSEaPqfg4aJF_vFoxLISP33EYPqFCCrX62iMTPmHQzNgDIctnG1osifmDJ9GpjMoO26OLESEl7vVMuwzXQ_jNCnKAOBwEhKUgyu8M7vWjkajC1AliPzhGgdVIrLXc5FUIAyuwCpMk0qtUCTmpHEig/https%3A%2F%2Fwww.nature.com%2Farticles%2Fs41596-023-00806-x
> <https://secure-web.cisco.com/1T2WKIelhUu8n_wuNR_5gb37RydVg1xrWEVvlUOdcz2UFWoSKKUO8S38cKClNmhjA13keedLju0Te5qHI70cMu686M8Nz33LUMDSdSc5UQuK9Vnl0b2bAEYY3X2pBwXbzaEnuzeLsX13HY7M0H8yFuJubt2wgqrSefkFgQ-XpZoEYeznT5PRAbS3RHGC8kb90A3ZS7a47Zkqz0eN6_e9ebZcBdzOsLlJYdHJwvolLCKV6tVC7wAXDgXrrODdQ0njXnHygdM9DOoJI9dttTpF1R-vsrPxiFk_XgPafiu9Hwi4KZHgiPrc6A_O77CgY-ZWPbJQ3x70CJAnfzYZH0pVTIg/https%3A%2F%2Fwww.nature.com%2Farticles%2Fs41596-023-00806-x>)
> is currently doing a better job than Infant FreeSurfer (*MailScanner has
> detected a possible fraud attempt from "secure-web.cisco.com" claiming to
> be* 
> https://secure-web.cisco.com/1HPPrCL1b3QKCLpwo16YY3JNX4_vaaGT3zbdy9V_gqn0jnJnQ2Iuz8tMhzRjb-OZpCYqa0HkFWnLi3Mm9KBpOJfV9HXA0fPheWIofIQ7YHZjURR-IGUWqerf2cwNkeunRYhigIZVideDrJmmwRqeq_LYvok5jDf77DzZXW_r-q0_vSLjm4bSnII3GxoCAE1yiJMbERYjIQbVMH4NpZDC_Pw5TDXfVvfR5H9Mufsq6zxKx2uipQMto4yyG5vyVFn9DyVlmgiBAfkyzGvTiypYpsPi1NQwtow7WHcHbm7vnOoyQ6wBDBoS-Zyd7TSBtuskF3LZDMT8YcbR5UL1jmV2etA/https%3A%2F%2Fwww.sciencedirect.com%2Fscience%2Farticle%2Fpii%2FS1053811920304328
> <https://secure-web.cisco.com/1M3VIitCZu3azJ0yMsnC2UvSxSyhcRAkbfTg_v1IUkOJObAa_9G8w89RsGzk4Xf4BqDqrThHe1m2ImfRWqSWSe4wlOc42k4ums78aO1mPCguwvDX3eo1ThINBtkLveTLfasbBXqMxqrNUHA3TRlxkBjgCG5QzeEm7x4lvFWhohrGB58jYWKeL3Hz7nABbqZe8ScicauUZNyGrLAEyNnuBsbcf70LrV4-qfWZ5YqKzqzUHzNE405Uh9rPMnIXfxNbNjc0Ybx41EhAM_VRu82VibTKLb1lMTckg3EU5UqhgXtYnTTiplYcAnHlQ1oqG0qxd_viSreko60OHp4aeerH74A/https%3A%2F%2Fwww.sciencedirect.com%2Fscience%2Farticle%2Fpii%2FS1053811920304328>)
> at cortical surface recon in this age range. However, I have custom surface
> ROI atlases on the fsaverage template and would very much like to leverage
> FreeSurfer's stats tools to map iBEAT subject surfaces to the template and
> extract stats from my custom ROIs.
>
> To do so, I need to
>
>    - Bring my iBEAT VTK surfaces (~47K vertices per hemisphere) into
>    FreeSurfer (mris_convert?)
>    - Resample the meshes to FS surface resolution (~12.5K vertices per
>    hemi)
>    - Compute curvature stats on the iBEAT surfaces
>    - Use mris_surf2surf (or equivalent) to map the resampled iBEAT
>    surface to the fsaverage template using curvature
>    - Compute/extract anatomical stats from the regions equivalent to the
>    ROIs in my template space atlas.
>
> All of the above are presenting challenges.
>
> So far, I tried using mris_convert to bring the VTK into FS format, but
> when I tried using mris_anatomical_stats, I get errors. Note that I also
> ran InfantFS on the data, and that I copied the converted (but not
> downsampled) VTK pial surface from iBEAT to
> $SUBJECTS_DIR/$subjid/surf/lh.ibeat.pial
>
> > mris_anatomical_stats -log test.ibeat.log sub-0001 lh ibeat.pial
>
> outputting results to test.ibeat.log...
> reading volume
> /scratch/t/tpaus/jhengeni/iBEAT/surface_remapping/sub-0001/sub-0001/mri/wm.mgz...
> reading input surface
> /scratch/t/tpaus/jhengeni/iBEAT/surface_remapping/sub-0001/sub-0001/surf/lh.ibeat.pial...
> Using TH3 vertex volume calc
> Total face volume 58098.7
> Total vertex volume 57390.1 (mask=0)
> reading input pial surface
> /scratch/t/tpaus/jhengeni/iBEAT/surface_remapping/sub-0001/sub-0001/surf/lh.pial...
> error: mrisReadTriangleFile opened
> /scratch/t/tpaus/jhengeni/iBEAT/surface_remapping/sub-0001/sub-0001/surf/lh.pial
> okay but surface doesn't match
> /scratch/t/tpaus/jhengeni/iBEAT/surface_remapping/sub-0001/sub-0001/surf/lh.ibeat.pial.
>  nvertices:46854 != mris->nvertices:12240 || nfaces:93704 !=
> mris->nfaces:24476
>
> error: mrisReadTriangleFile failed.
>
> error: mris_anatomical_stats: could not read surface file
> /scratch/t/tpaus/jhengeni/iBEAT/surface_remapping/sub-0001/sub-0001/surf/lh.pial
>
> It appears this approach is failing due to a difference in meshes (the
> iBEAT mesh is at a different resolution, with more elements, than the
> lh.pial from InfantFS).
>
> Any tips for this and the other steps would be a big help!
>
> Thanks,
> Jim H.
> --
> Jim Hengenius
> Research Associate
> Paus Research Group
> CHU Sainte-Justine Research Centre
> My Google Calendar
>
> <http://secure-web.cisco.com/1NyYdK1uBklqcDSsFMhBiedpJnPC-H6m83iYBtfLt-c-oKzBTwSVYkZEtHRJm20jit2kvoHh5jt8-_bpLVWG7tgMxbJ4Su8MGpCufWWEWDvHAcX9dhQFLeDfjGz9D7RqngZSEUkDY2_NJH_IRlQUylA-G1YdI1digS1CPmzl-aHQ6SiJqQJvSSdOCipdYr7g2ij2udg6xI_MLsXGPrJWj7HwQ_HUf0eVaAwHijfxy4oFp4s-6RbZ7r71HBbYQ4lJie5bkP-ZTBgBcbYgEvyAIAMOMIuJtr8nlS5kmo1o7-nwLg0ejnUKHsJ53Z_7FfyNC1OqSyUEoqzBUIUoZwgeCVg/http%3A%2F%2Fwww.google.com%2Fcalendar%2Fembed%3Fsrc%3Djames.hengenius%40gmail.com>
> _______________________________________________
> Freesurfer mailing list
> Freesurfer@nmr.mgh.harvard.edu
> https://secure-web.cisco.com/1wISCCgFsU7PzYbi6Cr7A5XjAj107lDep952lc-Cj9C0lCWzfeFepaNW1wt_Omx4XjGHIDBpteEXekvmMYk4a7bH4-Qmj3P5Bm7y2wbXruiPnkiz1y0fSD90Cc1AtRD4VOw1eo6EDMEHYARZmXYnqlVpy_dp0-Drn6w2sWKeTQGKrfLL4iWN0ujWDtaYNJFDklhjstNhtR6a85OL4WNkJWGN6G9Crkz4xlqW-xT0H1fwFNnSNBolodpXPbftlKn83cliW2EkOV7cpqFmSQ8hjw1J6p9t50hqJaYqX4cWHsDjPUFOGq2sezIRMDOOjjqzwiLmYBX9BJ8xr70uIP7V9gg/https%3A%2F%2Fmail.nmr.mgh.harvard.edu%2Fmailman%2Flistinfo%2Ffreesurfer
> 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 Mass General
> Brigham Compliance HelpLine at
> https://secure-web.cisco.com/1wcJY3sP01uyyfd3Us2oMk9lItq4WIJ5a-x7jh-EUG15zxyfksLnYO3vv24CCmasX-w_jyCiekj7CSItHfQCTvG0LxXwICPerAjgnzrolON3B5qYp1uVWmvCIAXysgMkZ4vjRqDmJAHtN12q3oiZmj2Db3ezAzDa2MT7_f8UKi-x4Mm1zfhe8DT8BDdQ54Nv0G64TuNS-lGwv5DYSTU_jvzvdwJlYhSZxdEzJp0p84d-vzSO6JWIMFtNsp7ag1fnB-NKXPo1VFRDSD8WDDOvIM_k0VAJ5EldVXg1mlJig8ljx7xhXmL2kZ-6YdmUunS1gzmOeTNSy8mD7iTtjmpue5A/https%3A%2F%2Fwww.massgeneralbrigham.org%2Fcomplianceline
>  <
> https://secure-web.cisco.com/1wcJY3sP01uyyfd3Us2oMk9lItq4WIJ5a-x7jh-EUG15zxyfksLnYO3vv24CCmasX-w_jyCiekj7CSItHfQCTvG0LxXwICPerAjgnzrolON3B5qYp1uVWmvCIAXysgMkZ4vjRqDmJAHtN12q3oiZmj2Db3ezAzDa2MT7_f8UKi-x4Mm1zfhe8DT8BDdQ54Nv0G64TuNS-lGwv5DYSTU_jvzvdwJlYhSZxdEzJp0p84d-vzSO6JWIMFtNsp7ag1fnB-NKXPo1VFRDSD8WDDOvIM_k0VAJ5EldVXg1mlJig8ljx7xhXmL2kZ-6YdmUunS1gzmOeTNSy8mD7iTtjmpue5A/https%3A%2F%2Fwww.massgeneralbrigham.org%2Fcomplianceline>
>  .
>


-- 
Jim Hengenius
Research Associate
Paus Research Group
CHU Sainte-Justine Research Centre
My Google Calendar
<http://secure-web.cisco.com/17NlTZvemZhoiUNDZxRb-sfRuh_CMF78ghWQ7Cq5D_b7ZdzW0yKM56Kt4MO0g3uqjlfdwo9AN-Dh7ts-6I_zk5BSQsUq1YA_Y8qfeCRTRg6ozeatWm-Zu4izpxfig2J22K9kE_7chFmbB90-YfzvT35t4nLe8y2seGc5IVAyKufTnfNY2NpXU29sI_VMW3pgC2aFFsLI50DFFWQ-6L8KLcnOFCUvGb9MbUBIuiLKCOm31tWnQSKV-h8mO6wEnEu3TQksPTRr5Mxwa6DnUbec6wFw6Rjv3RvjZ0eU_N1CjYO_ewi_aja_x0Q3YbMkFnHJTPfJkJv1R3ACBfLfaADSWcw/http%3A%2F%2Fwww.google.com%2Fcalendar%2Fembed%3Fsrc%3Djames.hengenius%40gmail.com>
_______________________________________________
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 Mass General Brigham 
Compliance HelpLine at https://www.massgeneralbrigham.org/complianceline 
<https://www.massgeneralbrigham.org/complianceline> .
Please note that this e-mail is not secure (encrypted).  If you do not wish to 
continue communication over unencrypted e-mail, please notify the sender of 
this message immediately.  Continuing to send or respond to e-mail after 
receiving this message means you understand and accept this risk and wish to 
continue to communicate over unencrypted e-mail. 

Reply via email to