1.首先区分三个通道,画出每个通道的直方图;
2.平滑直方图,像素值等于前后共七种颜色数量的平均值;
3.取波峰,波谷。(条件:颜色数亮s(i-2)<s(i-1)<s(i)>s(i+1)>s(i+2)或s(i-2)>s(i-1)>s(i)<s(i+1)<s(i+2)))
4.计算算有颜色的组合,并存入一个一维矩阵中;
5.颜色合并(利用map<int,int>);
6.把原图所有颜色与合并后的颜色计算距离,得出距离最近的合并颜色,在另一个原图大小的1通道矩阵中,记录相应位置颜色对应合并颜色的序号;
7.把原图序号相同的颜色加和,求平均值,就是类聚一次后的颜色。
8.把聚类后的颜色存入相应的序号的位置。
即为一次聚类后的结果,效果不很好。是由于聚类次数过少。
利用Kmean2算法,效果会很好。
#include "cv.h"
#include"highgui.h"
#include<map>
#include<vector>
#include<iostream>
using namespace std;
template<typename T>T sqr(T x){return x*x;}
CvMat* mat=0;
IplImage* image=0;
IplImage *histimg=cvCreateImage(cvSize(256,300),IPL_DEPTH_8U,3);
IplImage* DrawHist(map<int,int>colornum,int n);
map<int,int> Fine(CvMat* mat1u);
//CvMat* Merge(CvMat* Fine,int FineNum);
int main()
{
image=cvLoadImage("D:/7.jpg",1);
cvNamedWindow("Show",1);
cvShowImage("Show",image);
IplImage *hist0=0;//cvCreateImage(cvSize(300,300),IPL_DEPTH_8U,1);
IplImage *hist1=0;//cvCreateImage(cvSize(300,300),IPL_DEPTH_8U,1);
CvMat* r=cvCreateMat(image->height,image->width,CV_8UC1);
CvMat* g=cvCreateMat(image->height,image->width,CV_8UC1);
CvMat* b=cvCreateMat(image->height,image->width,CV_8UC1);
CvMat* rgb=cvCreateMat(image->height,image->width,CV_32FC3);
CvMat* result=cvCreateMat(image->height,image->width,CV_32FC3);
CvMat* result00=cvCreateMat(image->height,image->width,CV_8UC3);
cvConvertScale(image,rgb,1.0/255);
cvSplit(image,r,g,b,NULL);
map<int,int>fine_r=Fine(r);
map<int,int>fine