peaks,inds=rp.peak_finding(odf,odf_faces)
glob_norm_param=max(np.max(odf),glob_norm_param)
//remove the isotropic part
peaks = peaks - np.min(odf)
l=min(len(peaks),5)
QA[i][:l] = peaks[:l]IN[i][:l] = inds[:l]
//normalize
QA/=glob_norm_param
After Change
if mask.shape[0]==datashape[0]:
msk=mask.ravel().copy()
if mask == None:
mask=np.ones(datashape[:1])
msk=mask.ravel().copy()
glob_norm_param = 0
self.q2odf_params=q2odf_params
//Calculate Quantitative Anisotropy and