CalFrechetDist

// FileRead.cpp : Defines the entry point for the console application.
//

#include "stdafx.h"
#include <vector>
#include <iterator>
#include <fstream>
#include <iostream>
using namespace std;

void LoadLine(char* fname,vector<double>& generatedPoints);
double CalFrechetDist(char* fname1, char* fname2, long& pnum);

struct D_DOT {
double x;
double y;
};

typedef double* DBLPtr;

//void LoadLine(char* fname,vector<double>& generatedPoints)
//{
// char buf[1024];
// fstream infile;
// char sep[2]={0};
// double d;
// strcpy(sep,",");
//
// infile.open(fname,ios::in);
// while(1)
// {
// buf[0]=0;
// infile.getline(buf,1024);
// if(strlen(buf)<1)
// break;
//
// d=atof(strtok(buf,sep));
// generatedPoints.push_back(d);
//
// d=atof(strtok(NULL,sep));
// generatedPoints.push_back(d);
// }
//
// infile.close();
//}

void LoadLine(char* fname,vector<double>& generatedPoints)
{
//char buf[1024];
fstream infile;
char sep[2]={0};
double d;
strcpy(sep,",");

infile.open(fname); //,ios::in
printf("d:%s ",fname);

//char str[34];
//fstream fin;
//fin.open(fname);
//if (!fin)
//{
// cout<<"error "<<endl;
// exit(1);
//}
//while(fin.getline(str,34))
//{
// cout<<str<<endl;
// /*d=atof(strtok(str,sep));
// generatedPoints.push_back(d);
// printf("d:%d ",d);

// d=atof(strtok(NULL,sep));
// printf("d:%d ",d);
// generatedPoints.push_back(d); */
//}
//fin.close();

char buf[34];
while(infile.getline(buf,34))//1
{
//cout<<buf<<endl;
//printf("buf:%s ",buf);

if(buf!=" "&&strlen(buf)>=1)
{
d=atof(strtok(buf,sep));
//cout<<"d:"<<d<<endl;
generatedPoints.push_back(d);
d=atof(strtok(NULL,sep));
generatedPoints.push_back(d);
}
}
infile.close(); //很重要别忘了关
}

int _tmain(int argc, _TCHAR* argv[])
{
int i;
double d1;
vector<double> dlist;
vector<long> pnumlist;

char ss[100] = { 0 };
long pnum;
for (i = 1; i<40; i++)
{
sprintf(ss, "C:\Users\Administrator\Desktop\testtxt\TestTxt%d.txt", i + 1); //C:\Users\lenovo\Desktop\wo\YYI_%d.txt
d1 = CalFrechetDist("C:\Users\Administrator\Desktop\testtxt\TestTxt1.txt", ss, pnum); //double CalFrechetDist(char* fname1, char* fname2, long& pnum)
//cout << "距离:" << d1 << endl;
printf("距离：%d ",d1);
dlist.push_back(d1);
pnumlist.push_back(pnum);
}
system("pause");
return 0;
}

double pi_pj(D_DOT& pi, D_DOT& pj)
{
double dx = pi.x - pj.x;
double dy = pi.y - pj.y;

return sqrt(dx*dx + dy*dy);
}

double c(D_DOT *zb1, D_DOT *zb2, int i, int j, double **ca)
{
int curI, curJ;
double **caao;
caao = (double **)new DBLPtr[i + 1];
for (curI = 0; curI <= i; curI++)
{
caao[curI] = new double[j + 1];
}

for (curI = 1; curI <= i; curI++)
{
for (curJ = 1; curJ <= j; curJ++)
{
caao[curI][curJ] = pi_pj(zb1[curI], zb2[curJ]);
}
}

ca[1][1] = caao[1][1];
for (curJ = 1; curJ <= j; curJ++)
{
for (curI = 1; curI <= i; curI++)
{
if (curJ == 1)
{
if (curI>1)
{
ca[curI][curJ] = caao[curI][curJ];
if (ca[curI][curJ]<ca[curI - 1][1])
{
ca[curI][curJ] = ca[curI - 1][1];
}
}
}
else if (curI == 1)
{
if (curJ>1)
{
ca[curI][curJ] = caao[curI][curJ];
if (ca[curI][curJ]<ca[1][curJ - 1])
{
ca[curI][curJ] = ca[1][curJ - 1];
}
}
}
else if (curI>1 && curJ>1)
{
ca[curI][curJ] = caao[curI][curJ];
double tmp1 = ca[curI - 1][curJ];
double tmp2 = ca[curI - 1][curJ - 1];
double tmp3 = ca[curI][curJ - 1];
if (tmp1<tmp2&&tmp1<tmp3)
{
if (ca[curI][curJ]<tmp1)
ca[curI][curJ] = tmp1;
}
else if (tmp2<tmp1&&tmp2<tmp3)
{
if (ca[curI][curJ]<tmp2)
ca[curI][curJ] = tmp2;
}
else
{
if (ca[curI][curJ]<tmp3)
ca[curI][curJ] = tmp3;
}
}
else ca[curI][curJ] = 9999;
}
}

for (curI = 0; curI <= i; curI++)
{
delete[] caao[curI];
}
delete[] caao;

return ca[i][j];
}

double Frechet(D_DOT *zb1,int n1,D_DOT *zb2,int n2,double** fm)
{
int i,j;
double **caa;
double result;

caa = (double **)new DBLPtr[n1+1];
for(i=0;i<=n1;i++)
{
caa[i]=new double[n2+1];
}
for(i=0;i<=n1;i++)
{
for(j=0;j<=n2;j++)
{
caa[i][j] = -1;
}
}
result = c(zb1,zb2,n1,n2,caa);
for(i=0;i<=n1;i++)
{
for(j=0;j<=n2;j++)
{
fm[j][i] = caa[i][j];
}
}

for(i=0;i<=n1;i++)
{
delete [] (caa[i]);
}

delete [] caa;
return result;
}

double CalFrechetDist(char* fname1, char* fname2, long& pnum) //计算Frechet距离
{
int i;
long len1, len2;
double d;
vector<double> generatedPoints1;
vector<double> generatedPoints2;

LoadLine(fname1, generatedPoints1); //读取文件fname1中的内容到generatedPoints1
LoadLine(fname2, generatedPoints2);

len1 = generatedPoints1.size() / 2 + 1;
D_DOT *lxy1 = new D_DOT[len1];
for (i = 0; i<len1 - 1; i++)
{
lxy1[i + 1].x = generatedPoints1[i * 2];
lxy1[i + 1].y = generatedPoints1[i * 2 + 1];
}
len2 = generatedPoints2.size() / 2 + 1;
D_DOT *lxy2 = new D_DOT[len2];
for (i = 0; i<len2 - 1; i++)
{
lxy2[i + 1].x = generatedPoints2[i * 2];
lxy2[i + 1].y = generatedPoints2[i * 2 + 1];
}

DBLPtr* fm = new DBLPtr[len2]; //指向数组的指针?no!指针数组？no!指向指针数组的指针
for (i = 0; i<len2; i++)
{
fm[i] = new double[len1];
}

d = Frechet(lxy1, len1 - 1, lxy2, len2 - 1, (double**)fm);

delete[] lxy1;
delete[] lxy2;
for (i = 0; i<len2; i++)
{
delete[] fm[i];
}
delete[] fm;

generatedPoints1.clear();
generatedPoints2.clear();
pnum = len2;
return d;
}