【Opencv】⼤律法求阈值+分割图像计算像素⽐版权声明：本⽂为博主原创⽂章，转载请注明出处。 blog.csdn/u013617144/article/details/44863123 #include "opencv2/highgui/highgui.hpp"
#include "opencv2/imgproc/imgproc.hpp"
#include "opencv2/core/core.hpp"
#include <opencv\ml.h>
#include <iostream>
#include <vector>
#include <math.h>
#include <string.h>
using namespace std;
using namespace cv;;
//⼤津法求阈值
int otsu1(Mat frame)
{
int ws;
int ls;
int pixelCount[256]={0};
rspfloat pixelPro[256]={0};
int i, j, pixelSum = width * height, threshold = 0;
uchar* data = frame.data;
float w0, w1, u0tmp, u1tmp, u0, u1, deltaTmp, deltaMax = 0;
/
/统计每个灰度级中像素的个数
for(i = 0; i < height; i++)
{
for(j = 0;j < width;j++)
{
pixelCount[(int)data[i * width + j]]++;
}
}
自由交易//计算每个灰度级的像素数⽬占整幅图像的⽐例
for(i = 0; i < 256; i++)
{
pixelPro[i] = (float)pixelCount[i] / pixelSum;
}
//遍历灰度级[0,255]，寻合适的threshold
for(i = 0; i < 256; i++)
{
w0 = w1 = u0tmp = u1tmp = u0 = u1 = deltaTmp = 0;
for(j = 0; j < 256; j++)
{
if(j <= i) //背景部分
{
w0 += pixelPro[j];
u0tmp += j * pixelPro[j];
}
else //前景部分
{
w1 += pixelPro[j];
u1tmp += j * pixelPro[j];
}
}
u0 = u0tmp / w0;
日本海啸预警
u1 = u1tmp / w1;
deltaTmp = (float)(w0 *w1* pow((u0 - u1), 2)) ;
if(deltaTmp > deltaMax)
{
deltaMax = deltaTmp;
threshold = i;
}
}
return threshold;
return threshold;
}
//计算每个⼩块所占⽐例，存⼊向量
float Bili1(Mat src)
{
float BL,counter = 0;
int Sqaure = ls * ws;
Mat_<uchar>::iterator it = src.begin<uchar>();
Mat_<uchar>::iterator itend = d<uchar>(); for (; it!=itend; ++it)
{
if((*it)>0) counter+=1;
}
BL = counter/Sqaure;
return BL;
}
vector<float> Devide441(Mat img)
{
//归⼀化，灰度化，⼆值化
resize(img,img,Size(16,32));
cvtColor(img,img,COLOR_BGR2GRAY);
int th = otsu1(img);
cv::threshold(img,img,th,255,THRESH_BINARY);
//切割为4*4
int a = ls/4;
int b = ws/4;
vector<float> reback;
for(int i = 0; i<4 ;++i)
{
for(int j = 0; j<4 ;++j)
{
Rect roi_rect = Rect(Point(j*a,i*b),Size(a,b));
Mat roi = img(roi_rect);
float BL = Bili1(roi);
曾康霖
reback.push_back(BL);
}
}
return reback;
}
int main()
{
Mat img1 = imread("train/7.bmp",1);
imshow("img1",img1);
vector<float> a = Devide441(img1);
Mat A = Mat(a);
cout<<A<<endl;
Mat img2 = imread("LP8.jpg",1);
imshow("img2",img2);
vector<float> b = Devide441(img2);
Mat B = Mat(b);
cout<<B<<endl;
waitKey();
return 0;
公众安全感}
-——————————————————————————————————————————————————————————————
