I know. The translation is not so important to me, because this is going to just a little function that is working in the last step of the whole proces. The whole implementation wont be published and has just the meaning to implement a first alpha version of a face recognition system. It is nearly finished, and this is the only function I am missing. For this purpose I need it to work "EXACTLY".
I know that it doesnt look pythonic. After I know that it works correctly, I am going to write a total Rewrite for the Beta Version. On Sep 19, 1:55 pm, Peter Otten <[EMAIL PROTECTED]> wrote: > azrael wrote: > > Can somone look at this > > def otsu(hi): > > fmax=-1.0 > > border=len(hi) > > for i in range(border): > > if hi[i]!=0:break > > for j in range(border-1,0-1,-1): > > if hi[j] != 0:break > > s = sum([k*hi[k] for k in range(border)]) n = sum(hi) # > > product(im.size) > > n1=n2=csum=0.0 > > for k in range(i,j): > > n1 += hi[k] > > n2 = n - n1 > > csum+= k * hi[k] > > m1 = csum/ n1 > > m2 = (s - csum)/n2 > > As I said in my previous post, try replacing this line > > > sb = n1 * n2 * (m2 - m1) > > with the original > > sb = n1 * n2 * (m1 - m2) * (m1 - m2) > > that has been commented out. > > > if sb > fmax: > > fmax = sb > > V=k+1 > > print V > > I try to implement it from C from this location. > > Personally, I would start with a literal translation and not try to > improve it until it works. > > Peter -- http://mail.python.org/mailman/listinfo/python-list
