//?MyFirstCode.cpp?:?定義控制臺應用程序的入口點。
//

#include?“stdafx.h“

#include?
#include?
#include?
#include?
#include???
#include?
#include?
#include?
#include?


using?namespace?std;

class?JointPoint
{
public:
	double?jointPoint_X;
	double?jointPoint_Y;
	double?jointPoint_Z;

	JointPoint（）{}
	~JointPoint（）{}
};

//K均值聚類操作
void?CreateRandomArray（int?dataNum?int?clusterNum?int*?center）;//創建隨機初始聚類中心
void?CopyCenter（int?dimNum?int?clusterNum?double**?cluCenterDataCopy?double**?cluCenterData）;//復制聚類中心
void?AddToCluster（int?index?double*?inputData?int?dimNum?int*?Top?double***?clusterResult）;//加入一個數據到一個Cluster[index]集合
void?UpdateCluster（int?clusterNum?int?frameNum?double**?frameData?double**?Center?int?*Top?int?dim?double***?clusterResult）;
void?UpdateCenter（int?clusterNum?int?dim?double**?Center?int?*Top?double***?clusterResult）;//重新計算質心集合，對每一簇集合中的元素加總求平均即可
bool?IsEqual（double**?center1?double**?center2?int?clusterNum?int?dim）;//判斷兩個聚類中心是否相等，用來判斷是否結束
void?Print（int?clusterNum?int?dim?double**?Center?int?*Top?double***?clusterResult）;//打印聚合結果
//獲得骨架關節角度
double?AljointAngle（JointPoint?point1?JointPoint?point2?JointPoint?point3）;//計算關節角度

//獲取文件名
string?getFileName（int?actI?int?peoJ）;
string?getVideoFileName（int?act?int?peo）;
//讀取文件中的數據
bool?getFileData（string?fileName?int&?frameNum?int&?skeJointNum?double?**?&frameData）;

//獲取關鍵幀
void?getMainframe（int*?&?mainframe?int?frameNum?int?skeJointNum?int?clusterNum?double**?cluCenterData?double**?frameData）;
void?printMainframe（int*?mainframe?int?clusterNum?int?frameNum?int?act?int?peo）;//打印關鍵幀

int?_tmain（int?argc?_TCHAR*?argv[]）
{
	/*輸出文件流*/
	ofstream?outAngleFile;
	ofstream?outPositionFile;
	//打開要輸出的文件
	outAngleFile.open（“getAngleFeature20.csv“?ios::out?|?ios::trunc）;
	outPositionFile.open（“getPositionFeature20.csv“?ios::out?|?ios::trunc）;
	if?（!outAngleFile）return?1;
	if?（!outPositionFile）return?1;
	int?fileNum?=?0;
	int?frameNumTemp?=?0;
	/*****************************讀取所有文件**************************/
	for?（int?actI?=?1;?actI?<=?16;?actI++）
	{
		for?（int?peoJ?=?1;?peoJ?<=?10;?peoJ++）
		{
			//cout?<				/********************************************讀取數據************************************************/
				int?frameNum;
				int?skeJointNum;
				string?fileName;

				fileName?=?getFileName（actI?peoJ）;
				//cout?<
				double?**frameData;
				if?（!getFileData（fileNameframeNumskeJointNumframeData））{
					return?1;
				}
				int?dim?=?skeJointNum?*?3;				
				/**********************************************K均值聚類*******************************************/
				/*K均值聚類可以先選擇?k進行聚類，驗證結果是否正確
				*/
				//聲明變量
				int?clusterNum?=?frameNum/10;

				double**?cluCenterData;?//聚類中心
				cluCenterData?=?new?double*[clusterNum];
				for?（int?i?=?0;?i?