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,
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