?? k_averagedoc.cpp
字號(hào):
// K_averageDoc.cpp : implementation of the CK_averageDoc class
//
#include "stdafx.h"
#include "K_average.h"
#include "K_averageDoc.h"
#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif
/////////////////////////////////////////////////////////////////////////////
// CK_averageDoc
IMPLEMENT_DYNCREATE(CK_averageDoc, CDocument)
BEGIN_MESSAGE_MAP(CK_averageDoc, CDocument)
//{{AFX_MSG_MAP(CK_averageDoc)
ON_COMMAND(ID_MENUI_KAVERAGE, OnMenuiKaverage)
ON_COMMAND(ID_MENUI_SAVE_CLASSRESULT, OnMenuiSaveClassresult)
ON_COMMAND(ID_MENUI_RESET, OnMenuiReset)
ON_COMMAND(ID_MENUI_SAVERESULT, OnMenuiSaveresult)
ON_COMMAND(ID_MENUI_DUOLEI, OnMenuiDuolei)
ON_COMMAND(ID_MENUI_MULTIEPYSAVE, OnMenuiMultiepysave)
ON_COMMAND(ID_MENUI_HELP, OnMenuiHelp)
//}}AFX_MSG_MAP
END_MESSAGE_MAP()
/////////////////////////////////////////////////////////////////////////////
// CK_averageDoc construction/destruction
CK_averageDoc::CK_averageDoc()
{
// TODO: add one-time construction code here
}
CK_averageDoc::~CK_averageDoc()
{
}
BOOL CK_averageDoc::OnNewDocument()
{
if (!CDocument::OnNewDocument())
return FALSE;
// TODO: add reinitialization code here
// (SDI documents will reuse this document)
return TRUE;
}
/////////////////////////////////////////////////////////////////////////////
// CK_averageDoc serialization
void CK_averageDoc::Serialize(CArchive& ar)
{
if (ar.IsStoring())
{
// TODO: add storing code here
}
else
{
// TODO: add loading code here
}
}
/////////////////////////////////////////////////////////////////////////////
// CK_averageDoc diagnostics
#ifdef _DEBUG
void CK_averageDoc::AssertValid() const
{
CDocument::AssertValid();
}
void CK_averageDoc::Dump(CDumpContext& dc) const
{
CDocument::Dump(dc);
}
#endif //_DEBUG
/////////////////////////////////////////////////////////////////////////////
// CK_averageDoc commands
BOOL CK_averageDoc::OnOpenDocument(LPCTSTR lpszPathName)
{
if (!CDocument::OnOpenDocument(lpszPathName))
return FALSE;
FILE *file_swatch; //實(shí)驗(yàn)數(shù)據(jù)文件文件
file_swatch=fopen(lpszPathName,"rt"); //以讀的方式打開樣本文件
int num_swatch,num_class; //實(shí)驗(yàn)數(shù)據(jù)數(shù)目和類別數(shù)目
if (file_swatch==NULL)
{
AfxMessageBox("文件不能正常打開！",MB_OKCANCEL,NULL);
return false;
}
fscanf(file_swatch,"%d%d",&num_swatch,&num_class);
for (int i=0;i<num_swatch;i++)
{
fscanf(file_swatch,"%d%d",&feature_point.x,&feature_point.y);
vec_feature_point.push_back(feature_point);
}
UpdateAllViews(NULL);
// TODO: Add your specialized creation code here
fclose(file_swatch);
return TRUE;
}
void CK_averageDoc::OnMenuiKaverage()
{
BOOL flags=FALSE;
CPoint m_ptz1; //第一類初始中心
CPoint m_ptz2;//第二類初始中心
CPoint m_pt_tmp;
m_ptz1=vec_feature_point[0];//賦初始值
m_ptz2=vec_feature_point[1];
vec_point1.push_back(m_ptz1);
vec_point2.push_back(m_ptz2);
float x1,y1,x2,y2; //the center coordinate(x,y) of the two class
float x1_tmp,y1_tmp,x2_tmp,y2_tmp;//中心臨時(shí)變量
float dist1,dist2; //the distance to the center of the two class
HCURSOR m_hcur;
m_hcur=::SetCursor(LoadCursor(NULL,IDC_WAIT));
x1=m_ptz1.x;y1=m_ptz1.y;//中心初始賦值
x2=m_ptz2.x;y2=m_ptz2.y;//中心初始賦值
vec_point_class.clear();
do
{
for (int i=0;i<vec_feature_point.size();i++)
{
dist1=(x1-vec_feature_point[i].x)*(x1-vec_feature_point[i].x)+
(y1-vec_feature_point[i].y)*(y1-vec_feature_point[i].y);
dist2=(x2-vec_feature_point[i].x)*(x2-vec_feature_point[i].x)+
(y2-vec_feature_point[i].y)*(y2-vec_feature_point[i].y);
m_pt_tmp.x=vec_feature_point[i].x;
m_pt_tmp.y=vec_feature_point[i].y;
if (dist1<=dist2)
{
vec_point1.push_back(m_pt_tmp);
}
else
{
vec_point2.push_back(m_pt_tmp);
}
}
float total1=0;float total2=0;
for (int m=0;m<vec_point1.size();m++)
{
total1+=vec_point1[m].x;
total2+=vec_point1[m].y;
}
int num_vec1;
num_vec1=vec_point1.size();
x1_tmp=total1/num_vec1;
y1_tmp=total2/num_vec1;
float total3=0;float total4=0;
for (int n=0;n<vec_point2.size();n++)
{
total3+=vec_point2[n].x;
total4+=vec_point2[n].y;
}
int num_vec2;
num_vec2=vec_point2.size();
x2_tmp=total3/num_vec2;
y2_tmp=total4/num_vec2;
if (x1!=x1_tmp||y1!=y1_tmp||x2!=x2_tmp||y2!=y2_tmp)
{
x1=x1_tmp;y1=y1_tmp;x2=x2_tmp;y2=y2_tmp;
vec_point2.clear();
vec_point1.clear();
flags=true;
}
else
{
flags=false;
}
}
while(flags);
::SetCursor(m_hcur);
UpdateAllViews(NULL);
}
BOOL CK_averageDoc::OnSaveDocument(LPCTSTR lpszPathName)
{
// TODO: Add your specialized code here and/or call the base class
FILE *vec_file;
vec_file=fopen(lpszPathName,"wt");
int num_swatch,num_class=2;
num_swatch=vec_feature_point.size();
fprintf(vec_file,"%d %d\n",num_swatch,num_class);
for (int i=0;i<vec_feature_point.size();i++)
{
fprintf(vec_file,"%d %d\n",vec_feature_point[i].x,vec_feature_point[i].y);
}
fclose(vec_file);
return true;
// return CDocument::OnSaveDocument(lpszPathName);
}
void CK_averageDoc::OnMenuiSaveClassresult()
{
// TODO: Add your command handler code here
CFileDialog save_result_filedlg(FALSE,"txt",NULL,0,"(文本文件*.txt)|*.txt|DATA型文件(*.dat)|*.dat|所有文件(*.*)|*.*||",NULL);
CString sav_filname;
LPCTSTR lpszPathName;
FILE *sav_result_file;
int num_vec_point1,num_vec_point2,num_total;
int num_class=2;
num_total=vec_feature_point.size();
num_vec_point1=vec_point1.size();
num_vec_point2=vec_point2.size();
if (save_result_filedlg.DoModal()==IDOK)
{
sav_filname=save_result_filedlg.GetPathName();
lpszPathName=(LPCTSTR)sav_filname;
sav_result_file=fopen(lpszPathName,"wt");
fprintf(sav_result_file,"/**********************************************************************************/\n\n");
fprintf(sav_result_file," /*版權(quán)所有*/\n /*盜版不究*/\n /\n\n");
fprintf(sav_result_file,"/**********************************************************************************/\n\n");
fprintf(sav_result_file,"樣本總個(gè)數(shù)為:%d\n樣本分成了%d類\n第一類的樣本個(gè)數(shù)為:%d\n特征向量具體如下：\n",num_total,num_class,num_vec_point1);
for (int i=0;i<vec_point1.size();i++)
{
fprintf(sav_result_file,"%d %d\n",vec_point1[i].x,vec_point1[i].y);
}
fprintf(sav_result_file,"第二類樣本個(gè)數(shù)為：%d\n特征向量具體如下：\n",num_vec_point2);
for (int j=0;j<vec_point2.size();j++)
{
fprintf(sav_result_file,"%d %d\n",vec_point2[j].x,vec_point2[j].y);
}
fclose(sav_result_file);
return;
}
else if (save_result_filedlg.DoModal()==IDCANCEL)
{
return;
}
}
void CK_averageDoc::OnMenuiReset()
{
// TODO: Add your command handler code here
vec_point2.clear();
vec_point1.clear();
vec_feature_point.clear();
vec_point_class.clear();
UpdateAllViews(NULL);
}
void CK_averageDoc::OnMenuiSaveresult()
{
// TODO: Add your command handler code here
CFileDialog save_result_filedlg(FALSE,"txt",NULL,0,"(文本文件*.txt)|*.txt|DATA型文件(*.dat)|*.dat|所有文件(*.*)|*.*||",NULL);
CString sav_filname;
LPCTSTR lpszPathName;
FILE *sav_result_file;
int num_vec_point1,num_vec_point2,num_total;
int num_class=2;
num_total=vec_feature_point.size();
num_vec_point1=vec_point1.size();
num_vec_point2=vec_point2.size();
if (save_result_filedlg.DoModal()==IDOK)
{
sav_filname=save_result_filedlg.GetPathName();
lpszPathName=(LPCTSTR)sav_filname;
sav_result_file=fopen(lpszPathName,"wt");
fprintf(sav_result_file,"/**********************************************************************************/\n\n");
fprintf(sav_result_file," /*版權(quán)所有*/\n /*盜版不究*/\n /\n\n");
fprintf(sav_result_file,"/**********************************************************************************/\n\n");
fprintf(sav_result_file,"樣本總個(gè)數(shù)為:%d\n樣本分成了%d類\n第一類的樣本個(gè)數(shù)為:%d\n特征向量具體如下：\n",num_total,num_class,num_vec_point1);
for (int i=0;i<vec_point1.size();i++)
{
fprintf(sav_result_file,"%d %d\n",vec_point1[i].x,vec_point1[i].y);
}
fprintf(sav_result_file,"第二類樣本個(gè)數(shù)為：%d\n特征向量具體如下：\n",num_vec_point2);
for (int j=0;j<vec_point2.size();j++)
{
fprintf(sav_result_file,"%d %d\n",vec_point2[j].x,vec_point2[j].y);
}
}
fclose(sav_result_file);
}
/*
解決分成多類的問題，可以動(dòng)態(tài)更改分成的類別
*/
void CK_averageDoc::OnMenuiDuolei()
{
// TODO: Add your command handler code here
BOOL flags=FALSE;
CPoint m_pt_tmp;
CclassDlg m_classdlg;
HCURSOR m_hcur;
vec_point2.clear();
vec_point1.clear();
if (m_classdlg.DoModal()==IDOK)
{
m_class_num=m_classdlg.m_numclass;
m_first_center=m_classdlg.m_center;
float *x1= new float[m_class_num];
float*y1=new float[m_class_num]; //the center coordinate(x,y) of the two class
float *x1_tmp=new float[m_class_num];
float *y1_tmp=new float[m_class_num];//中心臨時(shí)變量
// Cvec_pt1 *vec_mypt=new Cvec_pt1[m_class_num];
m_hcur=::SetCursor(LoadCursor(NULL,IDC_WAIT));
for (int num=0;num<m_class_num;num++)//獲取初始特征中心
{
x1[num]=vec_feature_point[num+m_first_center-1].x;
y1[num]=vec_feature_point[num+m_first_center-1].y;
}
do
{
vec_point_class.resize(m_class_num);
for (int i=0;i<vec_feature_point.size();i++)
{
float *dist=new float[m_class_num]; //the distance to the center of the two class
for (int j=0;j< m_class_num;j++)
{
dist[j]=(x1[j]-vec_feature_point[i].x)*(x1[j]-vec_feature_point[i].x)+
(y1[j]-vec_feature_point[i].y)*(y1[j]-vec_feature_point[i].y);//計(jì)算到中心距離
m_pt_tmp.x=vec_feature_point[i].x; //獲取樣本特征向量
m_pt_tmp.y=vec_feature_point[i].y;
}
float dist_tmp=dist[0];
int dist_min_hao;
for (int k=0;k<m_class_num;k++) //遍歷所有類別
{
if (dist[k]<=dist_tmp) //比較到各個(gè)類別中心的距離
{
dist_tmp=dist[k];
dist_min_hao=k; //記錄符合要求類別號(hào)
}
}
delete[]dist;
vec_point_class[dist_min_hao].push_back(m_pt_tmp); //將樣本插入到符合要求的類別中
}
float *total1=new float[m_class_num];
float *total2=new float[m_class_num];
memset(total1,0,m_class_num*sizeof(float)); //設(shè)置初始值為零
memset(total2,0,m_class_num*sizeof(float));
flags=FALSE;//標(biāo)志中心是否改變，初始狀態(tài)為未改變
for (int vec1_num=0;vec1_num<m_class_num;vec1_num++) //遍歷所有樣本類別
{
for (int m=0;m<vec_point_class[vec1_num].size();m++)
{ total1[vec1_num]+=vec_point_class[vec1_num][m].x;//樣本特征求和
total2[vec1_num]+=vec_point_class[vec1_num][m].y;
}
int num_vec1;
num_vec1=vec_point_class[vec1_num].size(); //計(jì)算某一類中樣本數(shù)量
x1_tmp[vec1_num]=total1[vec1_num]/num_vec1;//計(jì)算新的中心
y1_tmp[vec1_num]=total2[vec1_num]/num_vec1;
if (x1[vec1_num]!=x1_tmp[vec1_num]||y1[vec1_num]!=y1_tmp[vec1_num]) //檢查中心是否更改
{
x1[vec1_num]=x1_tmp[vec1_num];
y1[vec1_num]=y1_tmp[vec1_num];
vec_point_class[vec1_num].clear(); //清除小容器中的樣本數(shù)據(jù)
flags=true;
}
// else{flags=false;}
}
if (flags==TRUE) {vec_point_class.clear();} //把大容器中的所有小容器清空
delete[]total1;//刪除變量
delete[]total2;
}
while(flags);
::SetCursor(m_hcur); //恢復(fù)光標(biāo)
delete[]x1;delete []y1; //刪除變量
delete[]x1_tmp;delete []y1_tmp;
}
UpdateAllViews(NULL);
}
//////////////////////////////////////////////////////////////////////////
/*
save the result of classfication
*/
void CK_averageDoc::OnMenuiMultiepysave()
{
// TODO: Add your command handler code here
CFileDialog save_result_filedlg(FALSE,"txt",NULL,0,"(文本文件*.txt)|*.txt|DATA型文件(*.dat)|*.dat|所有文件(*.*)|*.*||",NULL);
CString sav_filname;
LPCTSTR lpszPathName;
FILE *sav_result_file;
int num_vec_point,num_vec;
num_vec_point=vec_feature_point.size();
num_vec=vec_point_class.size();
if (save_result_filedlg.DoModal()==IDOK)
{
sav_filname=save_result_filedlg.GetPathName();
lpszPathName=(LPCTSTR)sav_filname;
sav_result_file=fopen(lpszPathName,"wt");
fprintf(sav_result_file,"/**********************************************************************************/\n\n");
fprintf(sav_result_file," /*版權(quán)所有*/\n /*盜版不究*/\n /\n\n");
fprintf(sav_result_file,"/**********************************************************************************/\n\n");
fprintf(sav_result_file,"樣本總個(gè)數(shù)為:%d\n分成的類別數(shù)目為：%d\n\n",num_vec_point,num_vec);
for (int i=0;i<vec_point_class.size();i++)
{
fprintf(sav_result_file,"/////////////////////////////////////////////////////////////////////////////////\n");
fprintf(sav_result_file,"這第%d類樣本具體結(jié)果如下：\n本類樣本數(shù)量為：%d\n特征向量為：\n",i+1,vec_point_class[i].size());
for (int j=0;j<vec_point_class[i].size();j++)
{
fprintf(sav_result_file,"%d %d\n",vec_point_class[i][j].x,vec_point_class[i][j].y);
}
}
}
fclose(sav_result_file);
}
void CK_averageDoc::OnMenuiHelp()
{
// TODO: Add your command handler code here
ShellExecute(NULL,"open","幫助文件.CHM",NULL,NULL,SW_SHOW);
}
?? 快捷鍵說明
復(fù)制代碼
Ctrl + C
搜索代碼
Ctrl + F
全屏模式
F11
切換主題
Ctrl + Shift + D
顯示快捷鍵
?
增大字號(hào)
Ctrl + =
減小字號(hào)
Ctrl + -