<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style id="owaParaStyle" type="text/css">
<!--
p
        {margin-top:0;
        margin-bottom:0}
p
        {margin-top:0;
        margin-bottom:0}
p
        {margin-top:0;
        margin-bottom:0}
--><!--
P {margin-top:0;margin-bottom:0;}--></style>
</head>
<body ocsi="0" fpstyle="1">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Until you get a better answer ...<br>
<br>
1) whatever subset you get out of Camino, it will not be a unique one (depending on the first direction that you pick; in fact, there may be 60_Choose_30 subset combinations of directions --- this is important if, say, you want to keep the signal in the z-axis
 ...); or,<br>
<br>
2) the first 30 directions in the list may actually be the ones you're after --- I believe they are listed so that if a scan is interrupted, the sampled directions would form a self-contained shell; or,<br>
<br>
3) you can design your very own shell of 30 dirs. by picking the first direction you like from the 60-shell, and eliminating its closest vector, then pick the next closest and eliminate its neighbour, and so on (via arccos of the dot product --- this is after
 all Coulomb's law, which is how the &quot;minimisation of electrostatic repulsion&quot; is based on, no?);
<br>
--- or you could even employ the 30 directions from Camino's PointSet, and find the closest vector of the 60-shell to each of the 30-shell (again, two loops on Matlab?).<br>
<br>
Uran<br>
PS: But it sounds a little dodgy that all that calculation is taking 5h ... even if it churned out all the subset combinations, it should be faster.<br>
<br>
<br>
<br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div style="direction: ltr;" id="divRpF715349"><font color="#000000" face="Tahoma" size="2"><b>From:</b> Grussu, Francesco<br>
<b>Sent:</b> 19 August 2013 7:59 PM<br>
<b>To:</b> Ferizi, Uran; camino-users@www.nitrc.org<br>
<b>Subject:</b> RE: subsetpoints output<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:#000000; font-size:10pt">Thanks for that.<br>
<br>
However my directions do have to be part of the initial 60 direction set, regardless how&nbsp; good or not it is.
<br>
This initial 60 direction set is something I was given with and I have real measures associated to it. Therefore I cannot change those.<br>
<br>
I just wondered how to interpret the outputs of <span style="font-style:italic">subsetpoints</span> .
<br>
<br>
Regards<br>
Francesco<br>
<br>
<div style="font-family:Times New Roman; color:rgb(0,0,0); font-size:16px">
<hr tabindex="-1">
<div id="divRpF253456" style="direction:ltr"><font color="#000000" face="Tahoma" size="2"><b>From:</b> Ferizi, Uran<br>
<b>Sent:</b> 19 August 2013 19:48<br>
<b>To:</b> Grussu, Francesco; camino-users@www.nitrc.org<br>
<b>Subject:</b> RE: subsetpoints output<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
Use Camino's point sets --- you have all the directions (well, from 3 to 256, and 755) ready for you. T<br>
In a terminal window, do:<br>
ls -l camino/PointSets (to obtain the full list Elec003.txt and so on).<br>
<br>
Is this what you wanted?<br>
<br>
<div style="font-family:Times New Roman; color:rgb(0,0,0); font-size:16px">
<hr tabindex="-1">
<div id="divRpF202196" style="direction:ltr"><font color="#000000" face="Tahoma" size="2"><b>From:</b> camino-users-bounces@www.nitrc.org [camino-users-bounces@www.nitrc.org] on behalf of Grussu, Francesco [francesco.grussu.12@ucl.ac.uk]<br>
<b>Sent:</b> 19 August 2013 7:34 PM<br>
<b>To:</b> camino-users@www.nitrc.org<br>
<b>Subject:</b> [Camino-users] subsetpoints output<br>
</font><br>
</div>
<div></div>
<div>
<div style="direction:ltr; font-family:Tahoma; color:rgb(0,0,0); font-size:10pt">
Dear Camino team,<br>
I am currently trying to select a subset of 30 directions from a set of 60 directions. I would like to get the combination of 30 out of 60 directions which is the most uniform possible. Therefore,&nbsp; I adopted the command<span style="font-style:italic"> subsetpoints</span>.
<br>
<br>
My code, which I run on the input file inputDirs.bvec (a text file with 60 rows, 3 colums, one row per direction)&nbsp; is something like<br>
<br>
<br style="font-style:italic">
<span style="font-style:italic">cat inputDirs.bvec | subsetpoints -singlesubset 30 -savestate myCurrentState &gt; outputdirs.bvec</span><br>
<br>
<br>
and after&nbsp; ~ 5hours of processing, it returns the output <span style="font-style:italic">
outputdirs.bvec</span> and a second file called <span style="font-style:italic">orderedPoints.dat</span> .
<br>
<br>
For what I understood from the guide, I think that both first 30 rows of <span style="font-style:italic">
orderedPoints.dat</span>&nbsp; and rows from the 2nd to 31st of<span style="font-style:italic"> outputdirs.bvec</span> are the subset of directions I am looking for. Can you confirm this?
<br>
<br>
Overall, first 30 rows of <span style="font-style:italic">orderedPoints.dat</span>&nbsp; and rows from the 2nd to 31st of<span style="font-style:italic"> outputdirs.bvec</span> are the same points but ordered differently. What's the difference between the two orders?
<br>
<br>
Thanks a lot for helping me out,<br>
all the best<br>
<br>
Francesco<br>
<br>
&nbsp;<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>