Freesurfer May 2007

freesurfer@nmr.mgh.harvard.edu

53 participants
97 discussions

vertex normals
by burak ozkalayci 08 May '07

08 May '07

Hi all, I want to get an ascii file listing the normal vectors of each vertex of the surface (white or pial). In http://www.mail-archive.com/freesurfer@nmr.mgh.harvard.edu/msg02210.html mri_surf2surf is adviced to used but I could not find any example usage for creating vertex normals. And I appreciate any help on converting the output of the mri_surf2surf to ascii file. One more question about the vertex normals. using tksurfer to see the surfaces, the normals of the selected vertices are printed on the console. to use mri_surf2surf the sphere.reg files are needed and they are created at the 3rd step of reconall. If there is any straightway to get those vertex normals from the surface definition like tksurfer do, I prefer it since I do not want to wait for the reconall 3. Thanks in advance. Burak

5 7

mri_surf2surf resampling does not give same number of nodes
by David Perlman 08 May '07

08 May '07

I think I must be misunderstanding something about the way mri_surf2surf works. I thought it was supposed to resample onto the target mesh, which I would think would mean that resampling two different data files onto the same target mesh would give outputs with the same numbers of nodes. But when I do this and then read the resampled data into Matlab, the nodes do not match up. Here are the commands I used: mri_surf2surf --srcsubject sp01t7 --hemi lh --srcsurfval thickness -- src_type curv --trgsubject sp01MULTt1 --trgsurfval ../../sp01MULTt1/ surf/t7tot1-lh.thickness.w --trg_type w mri_surf2surf --srcsubject sp01t6 --hemi lh --srcsurfval thickness -- src_type curv --trgsubject sp01MULTt1 --trgsurfval ../../sp01MULTt1/ surf/t6tot1-lh.thickness.w --trg_type w And here's what I see in Matlab: >> [vt6to1,ft6to1]=read_wfile('t6tot1-lh.thickness.w'); >> [vt7to1,ft7to1]=read_wfile('t7tot1-lh.thickness.w'); >> size(vt6to1) ans = 124463 1 >> size(vt7to1) ans = 124134 1 What am I missing here? -- -dave---------------------------------------------------------------- Science arose from poetry... when times change the two can meet again on a higher level as friends. -Göthe

2 3

(no subject)
by hankee 07 May '07

07 May '07

hi, Do FreeSurfer have Gray matter volume map,just like thickness map? Because I want to overlay the gray matter volume to inflate map,like thickness map? hankee，duhanjian(a)sohu.com 2007-05-05

1 0

(no subject)
by hankee 07 May '07

07 May '07

Hi, In the ?h.aparc.stats file,the Gray Matter Volume's units is mm,not mm^3. Is it right? Can you tell me the reason,because I want to use this option.Thank you very much! hankee，duhanjian(a)sohu.com 2007-05-05

2 1

ascii and mri_surf2surf??
by David Perlman 07 May '07

07 May '07

Forgive me if I'm simply missing the obvious, but as far as I can tell, mri_surf2surf is not capable of writing output in *any* of the surface formats listed in the Freesurfer surface formats document on the wiki, except .w paint format. Also, mris_convert doesn't accept .w as input; so there's no way at all to get output from mri_surf2surf into ascii format? I'm also confused about why mri_surf2surf only writes out in volume formats like analyze and mgh (other than .w). If there are surface versions of these filetypes, where can I find documentation about those other surface formats, since they're not in the existing surface formats document? Actually, this is all extra confusing because the output of --help for mri_surf2surf conflicts with the documentation on the wiki. The wiki says output can be curv format; the --help says not. It doesn't seem to be possible. Also some of the options have different names in the two locations. Just out of curiosity, why does it work like that? I'd think it would make more sense simply to have it write out in the standard surface formats (surface or curv or ascii)... All I'm trying to do is to subtract one set of curv-format data from another. (I have cortical thickness measures for the same subjects at two different times, and I want to be able to see the changes for individual subjects.) I have successfully registered the time 1 data onto the time 2 surface, and I can visualize it overlaid no problem. But I am struggling to find a way to subtract the aligned data... which seems like it should be easy... which makes me think I'm missing something obvious, but I've been searching through documentation for hours. I can see how to do it with Matlab, using mri_surf2surf to generate .w format output and reading it into Matlab, but I was hoping there would be a more straightforward way... Also, as an aside, the comments for read_wfile.m are actually the comments for write_wfile.m. -- -dave---------------------------------------------------------------- After all, it is not *that* inexpressible. -H.H. The Dalai Lama

3 2

Cortical sheet mask?
by David Perlman 07 May '07

07 May '07

I'm new to freesurfer, so I'm hoping this question isn't too obvious. I've searched the mailing list archives and the wiki but haven't found any leads. I'm wondering if there's a way to generate a volume representation of the cortical sheet, based on the data from the white surface and the pial surface, so that the sheet will be guaranteed to have good topology, as well as retaining the proper thickness at each point? It seems like all the information for this is there already since we have the two surfaces, but I don't know how to go about figuring out how to convert it into volume data. Here's a small image that might help give an idea of what I need: http://brainimaging.waisman.wisc.edu/~perlman/corticalsheet.png In case you're curious, the reason I want this is because we're trying to test the novel method for cortical thickness measurement given in "Three-Dimensional Mapping of Cortical Thickness Using Laplace's Equation", Jones, Buchbinder and Aharon, Human Brain Mapping 11:12-32(2000). -- -dave---------------------------------------------------------------- After all, it is not *that* inexpressible. -H.H. The Dalai Lama

3 5

Re: [Freesurfer] Cortical sheet mask?
by David Perlman 05 May '07

05 May '07

On May 5, 2007, at 8:37 AM, Gheorghe Postelnicu wrote: > PS: this might take a couple of days, as I can foresee small technical > problems occurring with the resolution increase. Speaking of small technical problems: I manually increased the resolution of brain.mgz and orig.mgz (I read in two different places conflicting statements about which of those two is used as the resolution template) and then ran mris_volmask. It ran for a while and then gave a "malloc failed" error: dyn-120-159:~/data/freesurfer/sp01MULTt1/mri omgoleus$ time mris_volmask --out_root laplaceribbon --label_left_white 3 -- label_right_white 3 --label_left_ribbon 2 --label_right_ribbon 2 sp01MULTt1 CTABreadASCII(/Applications/freesurfer/FreeSurferColorLUT.txt): Line 178: Duplicate structure index 222, was Ectorhinal 197 150 250 255 mris_volmask.bin(10271,0xa000cfc0) malloc: *** vm_allocate (size=536870912) failed (error code=3) mris_volmask.bin(10271,0xa000cfc0) malloc: *** error: can't allocate region mris_volmask.bin(10271,0xa000cfc0) malloc: *** set a breakpoint in szone_error to debug ERROR: In /usr/pubsw/packages/vtk/CVS/src/VTK/Common/ vtkDataArrayTemplate.txx, line 95 vtkFloatArray (0x8cbcaf80): Unable to allocate 134217728 elements of size 4 bytes. Bus error real 15m58.366s user 13m59.934s sys 0m8.270s dyn-120-159:~/data/freesurfer/sp01MULTt1/mri omgoleus$ On a different computer (both Intel Macs) it's been running for a couple of hours, but that computer only has 1 gig of memory and it has been paging like crazy. So I guess I have to wait for the modified code... thanks for the help with this -- -dave---------------------------------------------------------------- After all, it is not *that* inexpressible. -H.H. The Dalai Lama

1 0

glibc version compatability
by Paul Greenberg 05 May '07

05 May '07

2 1

Re: [Freesurfer] Cortical sheet mask?
by David Perlman 05 May '07

05 May '07

That sounds great, of course! Is mris_volmask a new script? It doesn't seem to be present in my installation, which I downloaded last week. Where can I get that program? When I fed a higher-resolution volume as template to mri_surf2vol with --fillribbon, it actually worked just fine, except for the same graininess bug mentioned. So if there's already a bug fix version of mri_surf2vol available I'd sure like to have it now so I can try it out. It is probably sufficient for my current needs. Where can I get the new mri_surf2vol? Thanks so much for your help... On May 5, 2007, at 8:37 AM, Gheorghe Postelnicu wrote: > It can do that, although this is not supported currently. > > Also, the current binary for the ribbon generation is called > mris_volmask. > > I think the best thing to do to satisfy your request is to introduce a > new parameter for voxel scaling. This would be applied to the template > voxel size (which is now by default orig.mgz). > > How would that sound? > > PS: this might take a couple of days, as I can foresee small technical > problems occurring with the resolution increase. > > On 5/5/07, Bruce Fischl <fischl(a)nmr.mgh.harvard.edu> wrote: >> Hi David, >> >> there are two issues here. One is a known bug in which the ribbon >> isn't >> filled properly. This has been fixed by Gheorghe Postelnicu and >> Nick can >> post you a new version. The resolution is a different issue. >> Gheorghe: will >> your code handle generating higher resolution versions of the ribbon? >> >> thanks, >> Bruce >> >> >> On Sat, 5 May 2007, David Perlman wrote: >> >> > Is there a way to re-do the generation of the ribbon but >> generate it with a >> > smaller voxel size? I've already had success resampling the >> T1.mgz to 0.25mm >> > and using it as a template for mri_surf2vol with --mkmask and -- >> fillribbon, >> > but the problem is that only projects out from the vertices, so >> the resulting >> > ribbon is sparse and thready, and doesn't fill the whole >> cortical ribbon. >> > (Likewise the ?h.ribbon.mgz file seems to be missing a lot of >> voxels in its >> > interior.) >> > >> > Ideally what I want is to be able to overlay the surfaces onto a >> volume of my >> > chosen resolution, and then fill the voxels between the >> surface. I've been >> > looking at matlab code, and found this program: >> > http://tinyurl.com/2bbyap >> > surface2volume.m which allegedly does what I want, but it >> doesn't actually >> > seem to work right when I try it: I get only a few scattered >> voxels filled. >> > >> > Does anyone have any idea of how to do this? Any help is greatly >> > appreciated! >> > >> > >> > On May 3, 2007, at 10:31 AM, Bruce Fischl wrote: >> > >> >> Hi David, >> >> >> >> we do generate such a thing called ?h.ribbon.mgz, but note that >> once you go >> >> back into the volume the topology is no longer preserved. >> >> >> >> cheers, >> >> Bruce >> >> On Thu, 3 May 2007, David Perlman wrote: >> >> >> >>> I'm new to freesurfer, so I'm hoping this question isn't too >> obvious. >> >>> I've searched the mailing list archives and the wiki but >> haven't found any >> >>> leads. >> >>> >> >>> I'm wondering if there's a way to generate a volume >> representation of the >> >>> cortical sheet, based on the data from the white surface and >> the pial >> >>> surface, so that the sheet will be guaranteed to have good >> topology, as >> >>> well as retaining the proper thickness at each point? It >> seems like all >> >>> the information for this is there already since we have the >> two surfaces, >> >>> but I don't know how to go about figuring out how to convert >> it into >> >>> volume data. >> >>> >> >>> Here's a small image that might help give an idea of what I need: >> >>> http://brainimaging.waisman.wisc.edu/~perlman/corticalsheet.png >> >>> >> >>> >> >>> In case you're curious, the reason I want this is because >> we're trying to >> >>> test the novel method for cortical thickness measurement given in >> >>> "Three-Dimensional Mapping of Cortical Thickness Using Laplace's >> >>> Equation", Jones, Buchbinder and Aharon, Human Brain Mapping >> >>> 11:12-32(2000). >> >>> >> >>> -- >> >>> - >> dave---------------------------------------------------------------- >> >>> After all, it is not *that* inexpressible. >> >>> -H.H. The Dalai Lama >> >>> >> >>> >> >>> >> >>> _______________________________________________ >> >>> Freesurfer mailing list >> >>> Freesurfer(a)nmr.mgh.harvard.edu >> >>> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer >> >>> >> >>> >> >>> >> >> _______________________________________________ >> >> Freesurfer mailing list >> >> Freesurfer(a)nmr.mgh.harvard.edu >> >> https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer >> > >> > -- >> > - >> dave---------------------------------------------------------------- >> > Science arose from poetry... when times change the two can meet >> again >> > on a higher level as friends. -Göthe >> > >> > >> > >> > _______________________________________________ >> > Freesurfer mailing list >> > Freesurfer(a)nmr.mgh.harvard.edu >> > https://mail.nmr.mgh.harvard.edu/mailman/listinfo/freesurfer >> > >> > >> > > > > -- > Gheorghe Postelnicu, PhD > MGH, Harvard Medical School -- -dave---------------------------------------------------------------- Science arose from poetry... when times change the two can meet again on a higher level as friends. -Göthe

2 1

(no subject)
by hankee 05 May '07

05 May '07

Hi, How many parts consisted of the occipital lobe in freesurfer's ?h.aparc.stats file and what the name of this parts in this file? Thank you! hankee，duhanjian(a)sohu.com 2007-05-05

2 1

← Newer
1
2
3
4
5
6
7
8
9
10
Older →

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

2006

2005

2004

Freesurfer May 2007