Алгоритм преобразования полутонового изображения в бинарное методом Отса

Шаг 1. Загрузить изображение Ris_bmp в объект Tbitmap.
Шаг 2. Задать высоту H и ширину W изображения.
Шаг 3. Задать одномерный массив p[i],i=1...256 размерности 256 элементов для определения частот каждого уровня яркости;
Шаг 4. Обнулить одномерный массив p[i] ;
Шаг 5. Задать i=0;
Шаг 6. Задать j=0;
Шаг 7. Определить значение яркости t пикселя изображения с координатами (i,j);
Шаг 8. Увеличить значение n[t] на единицу, т.е. n[t]=n[t]+1 ;
Шаг 9. Если значение j<Height то j=j+1 и перейти к шагу 7;
Шаг 10. Если значение i<Weidth то i=i+1 и перейти к шагу 6;
Шаг 11. Определить общее количество пикселей N на изображении N=HxW ;
Шаг 12. Задать t=1 ;
Шаг 13. Вычислить относительную частоту появления уровня яркости по формуле p[t]=n[t]/N;
Шаг 14. Если t<256 ,то t=t+1 и перейти к шагу 13;
Шаг 15. Задать максимальное значение max=-1000000 ;
Шаг 16. Задать начальное значение переменной t=3;
Шаг 17. Задать начальное значение суммы произведений частот на яркости второго класса S=0 ;
Шаг 18. Вычислить сумму произведений частот на яркости второго класса S=S+t*p[t] ;
Шаг 19. Если t<256, то t=t+1 и перейти к шагу 18;
Шаг 20. Задать начальное значение, суммы частот появления уровней яркости для первого класса ω1=p[1];
Шаг 21. Задать начальное значение суммы произведений частот на яркости первого класса S1=p[1] ;
Шаг 22. Задать начальное значение порога t=2 ;
Шаг 23. Вычислить сумму частот появления уровней яркости для первого класса ω1=ω1+p[t];
Шаг 24. Вычислить сумму частот появления уровней яркости для второго класса ω2=1-ω1 ;
Шаг 25. Вычислить сумму произведений частот на яркости для первого класса S1=S1+t*p[t] ;
Шаг 26. Рассчитать дисперсию между первым и вторым классами по формуле σ=(S1*ω2-s2*ω1)/ω1*ω2;
Шаг 28. Если max<σ, то max=σ и запомнить значение оптимального порога: tp=t ;
Шаг 29.Увеличить значение t=t+1 ;
Шаг 26. Вычислить сумму произведений частот на яркости для второго класса S2=S2-t*p[t] ;
Шаг 29. Если t<255 , перейти к шагу 23;
Шаг 30. Задать i=0 ;
Шаг 31. Задать j=0 ;
Шаг 32. Если Y[i,j]>tp ,Y[i,j]=255 то , иначе Y[i,j]=0 ;
Шаг 33. Если j<H, то j=j+1 и перейти к шагу 32
Шаг 34. Если перейти к шагу 31