1. 随机确定⼀个初始阈值 k1 = 127(可以为随机值【0，255】） 2. 根据阈值k1把原始图像分为两部分（⼩于 k1 的像素部分，⼤于k1的背景部分），并分别求其均值avgPix, avgBac 3. 计算新的阈值 k2 = (avgPix + avgBac) / 2 4. 判断 k1 与 k2 是否相等，如果不等，则更新k1, 重新执⾏步骤2-3 直到 k1 == k2
5. 根据计算得到的阈值 k1 (或 k2) 对图像进⾏⼆值化
#include "JpegDecoder.h"
#include <stdio.h>
#include <opencv2/highgui.hpp>
#include <math.h>
using namespace JpegCodec;
static cv::Mat ConvertToMat(Matrix &mat)
{
int channel = CV_8UC3;
if (mat.channal == 1) channel = CV_8UC1;
cv::Mat ws, ls, channel); // create a new matrix
for (int i = 0; i < ws * ls * mat.channal; i++)
{
img.data[i] = mat.data[i];
}
return img;
}
void ShowImage(Matrix &mat)
{
cv::Mat img = ConvertToMat(mat);
cv::imshow("Bitmap", img);
}
// ⼆值化
void Binary(Matrix &mat)
{
int k1 = 127, k2 = 0;
int cnt = ws * ls * mat.channal;
// 计算阈值
while (k1 != k2)
{
k2 = k1;
int avgPix = 0;
int avgBac = 0;
int size = 0;
for (int i = 0; i < cnt; i++)
{
if (mat.data[i] < k2)
{
avgPix += mat.data[i]; // 像素数据
size++;
}
else
董鸡
else
avgBac += mat.data[i]; // 背景数据
}
// 计算像素阈值与背景阈值中国选煤论坛
avgPix /= size;
avgBac /= (cnt - size);
// 计算 k1
k1 = (avgPix + avgBac) / 2;
}
// ⼆值化
for (int i = 0; i < cnt; i++)
{
if (mat.data[i] < k1) mat.data[i] = 0;
else mat.data[i] = 255;
}
}
// 灰度化
void Gray(Matrix &dst, Matrix &src)
{
wwer
dst.ws, ls, 1);
for (int i = 0; i < ws; i++)
{
for (int j = 0; j < ls; j++)
{
int idx = (i * ws + j) * 3;
dst.data[idx / 3] = (src.data[idx] + src.data[idx + 1] + src.data[idx + 2]) / 3; }
}
}
int main(int argc, char *arrv[])
{
JpegDecoder decoder("01.jpg");
Matrix mat, dst;
温州民房倒塌decoder.Decoder(mat);
Gray(dst, mat);
Binary(dst);
ShowImage(dst);
cvWaitKey(0);
return 0;
}
JpegDecoder:
ips细胞>屈服极限运⾏⽰例
转载于:my.oschina/tigerBin/blog/1488493
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议