大津阈值分割算法

#include <iostream.h>
#include <math.h>
#include <cv.h>
#include <highgui.h>
#include <stdio.h>
#include <windows.h>

//#define Struct_size 11

//int flag_color;
int hei;
int wid;


void otsu(IplImage* A, IplImage* B)
{
long N = hei * wid;
int h[256];
double p[256],u[256],w[256];

for(int i = 0; i < 256; i++)
{
   h[i] = 0;
   p[i] = 0;
   u[i] = 0;
   w[i] = 0;
}

for(i = 0; i < hei; i++)
{