Yes this part has not been finished. But the robust tesnor fit function
works. It just needs to be modularized

On Thu, Jan 31, 2013 at 12:07 PM, Ariel Rokem <notifications@github.com>wrote:

> In mrDiffusion/utils/dtiComputeImageNoise.m:
>
> > +if strcmp('corner', noiseCalcMethod)
> > +
> > +    % According to Henkelman (1985), the expected signal variance (sigma) can
> > +    % be computed as 1.5267 * SD of the background (thermal) noise.
> > +    sz = size(dwRaw.data);
> > +    x  = 10;
> > +    y  = 10;
> > +    z  = round(sz(3)/2);
> > +
> > +    [x,y,z,s] = ndgrid(x-5:x+5, y-5:y:5, z-5:z+5, 1:sz(4));
> > +    noiseInds = sub2ind(sz, x(:), y(:), z(:), s(:));
> > +    sigma     = 1.5267 * std(double(dwRaw.data(noiseInds)));
> > +
> > +elseif strcmp('b0', noiseCalcMethod)
> > +    % Pull out the data from the brain mask
> > +    data =
>
> This is still broken, no?
>
> —
> Reply to this email directly or view it on GitHub<https://github.com/vistalab/vistasoft/pull/4/files#r2850004>.
>
>