Hi
> Hi Saad, thank you for your response and the script. I have another question:
>
> 1. You said each row of matrix_seeds_to_all_targets Is the number of streamlines from that Voxel that reached the target (each row). In which order are the voxels in that array (the row)? I mean, does the first value of that row corresponds to the first Voxel connected to target that is closer to [1,1,1] (matlab coordinates)? So, if I go over my mask from i=1,j=1,k=1 until i=xsize,j=ysize,k=zsize the first Voxel in that buckle that is connected to the target is going to correspond to the first value in that row?
the rows of this matrix are in the same order as voxels in your mask if you read them in matlab like this: (note that matlab voxels start at 1 but and fslview voxels start at 0).
mask=read_avw('seedmask');
[x,y,z]=sub2ind(size(mask),find(mask));
> 2. If I do that loop, the values are different from the ones in fdt_paths. Why? Although the same # of voxels are different than zero.
I don't understand why you expect fdt_paths and the seed_to_all_targets matrix to give you the same values. The matrix stores the number of samples that reach your targets, while fdt_paths stores number of samples throughout the brain.
> 3. If I add the values of the row in matrix_seeds_..etc it gives the same number as waytotal. But if I add all the numbers next to simbol # in saved_paths it gives a very bigger number. Why? Isn't that the number of streamlines that were not reflected and waytotal the # of generated tracks... Isn't that the same? (sorry if is dumb question)
the "#" symbol in the saved_paths file indicates the number of xyz locations for each streamlines (related to the length of the streamline), not the number of samples.
Cheers
Saad
>
--
Saad Jbabdi
University of Oxford, FMRIB Centre
JR Hospital, Headington, OX3 9DU, UK
(+44)1865-222466 (fax 717)
www.fmrib.ox.ac.uk/~saad
|