1a8468988ae35b6d3814288018fc2f9056e74974,dipy/reconst/gqi.py,GeneralizedQSampling,__init__,#GeneralizedQSampling#Any#Any#Any#Any#Any#Any#Any#,54

Before Change

                    QA[i][:l] = peaks[:l]
                    IN[i][:l] = inds[:l]

        if mask==None:
            for (i,s) in enumerate(S):                            
                //Q to ODF
                odf=np.dot(s,q2odf_params)            
                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
       
        if len(datashape) == 4:
            self.QA=QA.reshape(x,y,z,5)    

After Change

            QA = np.zeros((x,5))
            IN = np.zeros((x,5))
            if mask != None:
                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

Instances