If I remember correctly, this routine updates the estimates of the Gaussians, and is part of spm_segment.m. In SPM5, spm_segment.m is not actually used. Take a look at spm_preproc.m, which is described in: J. Ashburner and K.J. Friston. Unified segmentation. NeuroImage, 26:839-851, 2005. Best regards, -John On Tuesday 20 May 2008 14:10, Wang wrote: > Hi, professors, > I am a newer in spm ,now I have some questions to ask. > Now I am analysing the code of segment.I have some question to ask .In > the function CP = update_cp_est(CP,P,dat,msk,p) , CP.mom2(:,:,i,p) in the > line 525,and it related code is CP.cv(:,:,i).I cost many time to analyse it > ,but I failed ,may I don¡¯t know the algorithm of segment.I really need for > your help.Inthe In the function BP = init_bp(VF,co,reg), BP.IC0 is said to > Cost function based on sum of squares of 4th derivatives.I also cost my > many time to analyse,but I think I can¡¯t related to the algorithm. Above > all,I need the algorithm of the segment.May be it is easy to you ,but I > really greet for your help. Thank you > > > --------------------------------- > ÑÅ»¢ÓÊÏ䣬ÄúµÄÖÕÉúÓÊÏ䣡