<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=us-ascii"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:#0563C1;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:#954F72;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri",sans-serif;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link="#0563C1" vlink="#954F72"><div class=WordSection1><p class=MsoNormal>Dear FT experts,<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>I am not able to realign a T1 scan using the fiducial points and headshapes of a 4D/BTI system, can anyone help me on this? The “interactive” realignment helped a lot, but the end result was not perfect!<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Here’s what I have tried,:<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>individual_mri = ft_read_mri(mripath); % ..\MEG_4D\C-101_T1.nii</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>ft_sourceplot([], individual_mri);<o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:12.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg = [];</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.method = </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'flip'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>;</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>mri.coordsys = </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'4d'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>;</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>individual_mri = ft_volumereslice(cfg, individual_mri);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>ft_sourceplot([], individual_mri);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>individual_mri = ft_convert_units(individual_mri, </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'m'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>headshape = ft_read_headshape(headpath); % ..\MEG_4D\hs_file</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>mri_fids = headshape.fid.pos;<o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'><o:p> </o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg = [];</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.fiducial.nas = mri_fids(1,:);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.fiducial.lpa = mri_fids(2,:);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.fiducial.rpa = mri_fids(3,:);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.coordsys = </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'4d'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>;</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.method = </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'fiducial'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>;</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>mri_realigned = ft_volumerealign(cfg, individual_mri);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg = [];</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.coordsys = </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'4d'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>;</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.headshape = headshape;</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>cfg.method = </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'headshape'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>;</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>mri_realigned = ft_volumerealign(cfg, individual_mri);<o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'><o:p> </o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:forestgreen'>% check that the MRI is consistent after realignment</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>ft_determine_coordsys(mri_realigned, </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'interactive'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>, </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>'no'</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>hold </span><span style='font-size:10.0pt;font-family:"Courier New";color:#A020F0'>on</span><span style='font-size:10.0pt;font-family:"Courier New";color:black'>; </span><span style='font-size:10.0pt;font-family:"Courier New";color:forestgreen'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>drawnow; </span><span style='font-size:10.0pt;font-family:"Courier New";color:forestgreen'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:10.0pt;font-family:"Courier New";color:black'>ft_plot_headshape(headshape);</span><span style='font-size:12.0pt;font-family:"Courier New"'><o:p></o:p></span></p><p class=MsoNormal style='text-autospace:none'><span style='font-size:12.0pt;font-family:"Courier New"'><o:p> </o:p></span></p><p class=MsoNormal>Data can be found at,<o:p></o:p></p><p class=MsoNormal>https://www.dropbox.com/sh/e74pt2uc08q6z4q/AAB3fkp_9gvJxi-YtO0muRSca?dl=0<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal>Best, <o:p></o:p></p><p class=MsoNormal>Vahab<o:p></o:p></p><p class=MsoNormal><o:p> </o:p></p><p class=MsoNormal><o:p> </o:p></p></div></body></html>