学习OpenCV——一次颜色聚类

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