<div dir="ltr"><div><div>Hi Arjen,<br><br>Thanks for your message. I did align the mri to Talairach; as you can see from <a href="http://i.imgur.com/26nyHYZ.png">http://i.imgur.com/26nyHYZ.png</a>, the volume conduction model and sourcespace are both expressed in the same coordinate system (i.e., everything's pointing in the same direction) but they're just not sitting on top of one another. If anyone has any ideas on where that problem was introduced (or how to re-align them now), I would greatly appreciate it. Below I have some more details about how I processed that data, if it helps.<br>
<br></div>I'm trying to go through the data one step at a time and track where the problem might have happened. When I compare the sourcespace before having applied any transformation (i.e., the headshape from <subject>-oct-6-src.fif) to the original mri (orig-nomask.mgz), they look ok (I don't know how to plot them together, but see <a href="http://i.imgur.com/LGW7YnJ.png">http://i.imgur.com/LGW7YnJ.png</a> and <a href="http://i.imgur.com/JUoTxc9.png">http://i.imgur.com/JUoTxc9.png</a> -- things at least look like they're on more or less the same plane). Then I re-register the mri to CTF (<a href="http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate#source_modelco-registration_of_the_source_space_to_the_sensor-based_head_coordinate_system">http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate#source_modelco-registration_of_the_source_space_to_the_sensor-based_head_coordinate_system</a>); after that, in mri_nom_ctf, the axes are all going in the right direction but the whole head is tilted forward and the origin of the axes is no longer at the anterior commisure (see <a href="http://i.imgur.com/CTZNOTk.png">http://i.imgur.com/CTZNOTk.png</a> for the realigned MRI).  Applying the transformation matrix T to the sourcespace also seems to tilt it like that (<a href="http://i.imgur.com/bcNIpf3.png">http://i.imgur.com/bcNIpf3.png</a>), although as can be seen from the first image in this message it doesn't quite line up with the volume conduction model in the end. As for the volume conduction model, here is what it looks like at first (<a href="http://i.imgur.com/ZX3m38b.png">http://i.imgur.com/ZX3m38b.png</a>) and here is what it looks like after applying the transformation matrix (<a href="http://i.imgur.com/vXa3Cnc.png">http://i.imgur.com/vXa3Cnc.png</a>). Obviously the transformation matrix is doing something, but it's not getting the sourcespace and volume conduction model lined up; since it's the same transformation matrix, all I can guess is that there was some pre-existing difference between the source mesh (the .fif file) and the anatomical mri (orig-nomask.mgz), but I'm not sure when that came in.<br>
<br></div>Another minor issue: when I first compared the volume conduction model and the sourcespace, they were expressed in different units even though I followed the code in the tutorial. See <a href="http://i.imgur.com/orwgcTJ.png">http://i.imgur.com/orwgcTJ.png</a>: the sourcespace looks 10x smaller than the volume conduction model, which I assume is because it is expressed in cm whereas the volume conduction model is expressed in mm. To get the figure linked at the very beginning of this message, I had to convert the units of the volume conduction model to cm, even though that's not in the tutorial.<br>
<br>I notice that the tutorial on the wiki hasn't been edited since October 2012 (other than a few edits I made this month which were just correcting typos in the prose). Is it possible that what's on the wiki is out of date? (Also cc'ing Lilla on this.)<br>
<br>Thanks,<br>Steve<br><div><div><br><div><br><br>> Message: 1<br>> Date: Fri, 31 May 2013 08:11:23 +0200 (CEST)<br>> From: "Stolk, A." <<a href="mailto:a.stolk@fcdonders.ru.nl">a.stolk@fcdonders.ru.nl</a>><br>
> To: FieldTrip discussion list <<a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a>><br>> Subject: Re: [FieldTrip] Question about minimum norm estimate pipeline<br>> Message-ID:<br>>         <<a href="mailto:1914237354.1292588.1369980683984.JavaMail.root@sculptor.zimbra.ru.nl">1914237354.1292588.1369980683984.JavaMail.root@sculptor.zimbra.ru.nl</a>><br>
> Content-Type: text/plain; charset="utf-8"<br>><br>> Hi Steve, A quick guess; did you correctly align your resliced mri to Talairach space by indicating the commissures ( <a href="http://imaging.mrc-cbu.cam.ac.uk/imaging/FindingCommissures">http://imaging.mrc-cbu.cam.ac.uk/imaging/FindingCommissures</a> ) and, if I'm correct, a point in the same place, e.g. between the hemispheres? This should update the transformation matrix. Best regards, Arjen ----- Oorspronkelijk bericht -----<br>
> > Van: "Stephen Politzer-Ahles" <<a href="mailto:politzerahless@gmail.com">politzerahless@gmail.com</a>><br>> > Aan: <a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a><br>> > Verzonden: Vrijdag 31 mei 2013 05:53:45<br>
> > Onderwerp: [FieldTrip] Question about minimum norm estimate pipeline<br>> > Hello all,<br>> > I have not yet gotten a response to my question below, but in the<br>> > meantime I have another question about the minimum norm estimate<br>
> > workflow--specifically, about the coordinate system for the<br>> > skull-stripped anatomy in the step described at<br>> > <a href="http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate#preprocessing_of_the_anatomical_mrisave_to_disk">http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate#preprocessing_of_the_anatomical_mrisave_to_disk</a><br>
> > . I'm confused by the following bit of code:<br>> > % ensure that the skull-stripped anatomy is expressed in the same<br>> > coordinate system as the anatomy<br>> > seg.transform = mri_tal.transform;<br>
> > In my data, mri_tal.coordsys is 'spm' (I presume this is the result of<br>> > re-aligning to Talairach in the previous step?) whereas seg.coordsys<br>> > is 'ctf' (as a result of re-aligning to CTF several steps earlier).<br>
> > (But mri_tal also has a field mri_tal.transformorig, which seg does<br>> > not have.) So should I really be using the same transform for both, as<br>> > shown in the tutorial?<br>> > Apologies if this question is pretty basic; I'm just trying to<br>
> > pinpoint where the mis-alignment described in my message below<br>> > occurred, so I want to make sure I understand each step of the<br>> > workflow correctly<br>> > Best,<br>> > Steve<br>
> > > Message: 1<br>> > > Date: Sat, 25 May 2013 08:11:18 -0500<br>> > > From: Stephen Politzer-Ahles < <a href="mailto:politzerahless@gmail.com">politzerahless@gmail.com</a> ><br>> > > To: <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>
> > > Subject: [FieldTrip] Sourcespace and volume conductor misaligned<br>> > > Message-ID:<br>> > > <CAJT2k_9-hd_sM=hp4P-CUu+=aduSOMyZV7XPG0= <a href="mailto:fFk0ouR0wzA@mail.gmail.com">fFk0ouR0wzA@mail.gmail.com</a><br>
> > > ><br>> > > Content-Type: text/plain; charset="utf-8"<br>> > ><br>> > > Hello all,<br>> > ><br>> > > I am going through the workflow at<br>> > > <a href="http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate">http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate</a> . After<br>
> > > making<br>> > > the volume conduction model using ft_prepare_headmodel(), I noticed<br>> > > that<br>> > > although the volume conduction model and sourcespace have the same<br>> > > orientation and overall size/shape (after I converted the volume<br>
> > > conduction<br>> > > model to cm, which wasn't in the tutorial but my original model came<br>> > > out in<br>> > > mm), they don't quite line up, as you can see in this figure:<br>
> > ><br>> > > <a href="http://i.imgur.com/mGEtLOa.png">http://i.imgur.com/mGEtLOa.png</a><br>> > ><br>> > > I did interactively re-align the data to CTF (twice--in step 2 of<br>> > > "Preprocessing of the anatomical MRI" and in step 4 of "Source<br>
> > > model")<br>> > > using fiducials, and to Talairach (step 5 of "Preprocessing of the<br>> > > anatomical data"), so I'm not sure how it ended up this way. The<br>> > > code I've<br>
> > > used at each step is basically the same as that in the tutorial.<br>> > ><br>> > > Is there any way to line up my volume conduction model and<br>> > > sourcespace now,<br>> > > without going back and re-running most of the workflow?<br>
> > ><br>> > > Best,<br>> > > Steve<br>> > ><br>> > > --<br>> > > Stephen Politzer-Ahles<br>> > > University of Kansas<br>> > > Linguistics Department<br>
> > > <a href="http://people.ku.edu/~sjpa/">http://people.ku.edu/~sjpa/</a><br>> > On Sat, May 25, 2013 at 1:56 PM, < <a href="mailto:fieldtrip-request@science.ru.nl">fieldtrip-request@science.ru.nl</a> ><br>
> > wrote:<br>> > ><br>> > > Send fieldtrip mailing list submissions to<br>> > > <a href="mailto:fieldtrip@science.ru.nl">fieldtrip@science.ru.nl</a><br>> > ><br>> > > To subscribe or unsubscribe via the World Wide Web, visit<br>
> > > <a href="http://mailman.science.ru.nl/mailman/listinfo/fieldtrip">http://mailman.science.ru.nl/mailman/listinfo/fieldtrip</a><br>> > > or, via email, send a message with subject or body 'help' to<br>
> > > <a href="mailto:fieldtrip-request@science.ru.nl">fieldtrip-request@science.ru.nl</a><br>> > ><br>> > > You can reach the person managing the list at<br>> > > <a href="mailto:fieldtrip-owner@science.ru.nl">fieldtrip-owner@science.ru.nl</a><br>
> > ><br>> > > When replying, please edit your Subject line so it is more specific<br>> > > than "Re: Contents of fieldtrip digest..."<br>> > ><br>> > ><br>> > > Today's Topics:<br>
> > ><br>> > > 1. Sourcespace and volume conductor misaligned<br>> > > (Stephen Politzer-Ahles)<br>> > > 2. Re: fieldtrip Digest, Vol 30, Issue 31 (Johanna Zumer)<br>> > ><br>> > ><br>
> > > ----------------------------------------------------------------------<br>> > ><br>> > > Message: 1<br>> > > Date: Sat, 25 May 2013 08:11:18 -0500<br>> > > From: Stephen Politzer-Ahles < <a href="mailto:politzerahless@gmail.com">politzerahless@gmail.com</a> ><br>
> > > To: <a href="mailto:fieldtrip@donders.ru.nl">fieldtrip@donders.ru.nl</a><br>> > > Subject: [FieldTrip] Sourcespace and volume conductor misaligned<br>> > > Message-ID:<br>> > > <CAJT2k_9-hd_sM=hp4P-CUu+=aduSOMyZV7XPG0= <a href="mailto:fFk0ouR0wzA@mail.gmail.com">fFk0ouR0wzA@mail.gmail.com</a><br>
> > > ><br>> > > Content-Type: text/plain; charset="utf-8"<br>> > ><br>> > > Hello all,<br>> > ><br>> > > I am going through the workflow at<br>> > > <a href="http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate">http://fieldtrip.fcdonders.nl/tutorial/minimumnormestimate</a> . After<br>
> > > making<br>> > > the volume conduction model using ft_prepare_headmodel(), I noticed<br>> > > that<br>> > > although the volume conduction model and sourcespace have the same<br>> > > orientation and overall size/shape (after I converted the volume<br>
> > > conduction<br>> > > model to cm, which wasn't in the tutorial but my original model came<br>> > > out in<br>> > > mm), they don't quite line up, as you can see in this figure:<br>
> > ><br>> > > <a href="http://i.imgur.com/mGEtLOa.png">http://i.imgur.com/mGEtLOa.png</a><br>> > ><br>> > > I did interactively re-align the data to CTF (twice--in step 2 of<br>> > > "Preprocessing of the anatomical MRI" and in step 4 of "Source<br>
> > > model")<br>> > > using fiducials, and to Talairach (step 5 of "Preprocessing of the<br>> > > anatomical data"), so I'm not sure how it ended up this way. The<br>> > > code I've<br>
> > > used at each step is basically the same as that in the tutorial.<br>> > ><br>> > > Is there any way to line up my volume conduction model and<br>> > > sourcespace now,<br>> > > without going back and re-running most of the workflow?<br>
> > ><br>> > > Best,<br>> > > Steve<br>> > ><br>> > > --<br>> > > Stephen Politzer-Ahles<br>> > > University of Kansas<br>> > > Linguistics Department<br>
> > > <a href="http://people.ku.edu/~sjpa/">http://people.ku.edu/~sjpa/</a></div></div></div></div>