//?DouglasPeuker.cpp?:?Defines?the?entry?point?for?the?console?application.
//

#include?“stdafx.h“
#include?
#include?
#include?
using?namespace?std;

typedef?struct?point
{
	double?x;
	double?y;
}POINT;

vector?initialPoints;
vector?FinalPoints;
bool?*bTag=NULL;
//////////////////////////////////////////////////////////////////////////
void?Readdata（char?*）;
void?DouglasPeuker（int?int?double）;
double?cacuDistance（int?int?int）;
void?Writedata（char?*）;

//////////////////////////////////////////////////////////////////////////
int?main（int?argc?char*?argv[]）
{
	int?nums=0;
	char?path1[]=“.\\jpl_vegetation_grass2.txt“;//打開(kāi)文件的路徑
	char?path2[]=“.\\2.txt“;//輸出文件的路徑
	double?torelance;
	int?i;
	int?nums2=0;

	Readdata（path1）;

	nums=initialPoints.size（）;
	bTag=（bool?*）malloc（sizeof（bool）*nums）;
	memset（bTag0nums）;
	torelance=0.00002;

????DouglasPeuker（0nums-1torelance）;

	for?（i=0;i	{
		if（bTag[i]）
		{
			FinalPoints.push_back（initialPoints[i]）;
		}
	}

	Writedata（path2）;
/*	nums2=FinalPoints.size（）;
	printf（“final?data:%d?\n“nums2）;
	for?（i=0;i	{
	
		printf（“%lf?%lf\n“FinalPoints[i].xFinalPoints[i].y）;
	}*/
	free（bTag）;
	return?0;
}
void?Readdata（char?*filename）
{

	char?predata1[256];
	char?predata2[256];
	char?predata3[256];

	FILE?*filein;
//	int?i;
	POINT?point0;
	//float?xiyi;

	memset（predata10256）;
	memset（predata20256）;
	memset（predata30256）;
	filein=fopen（filename“r“）;
//	for?（i=0;i<3;i++）
	{
		fscanf（filein“%[^\n]c\n“predata1）;
	//	fscanf（filein“%[^\n]c\n“predata2）;
	//	fscanf（filein“%[^\n]c\n“predata3）;
	}


	point0.x=0.00000;//啟用浮點(diǎn)庫(kù)，不加報(bào)錯(cuò)
	point0.y=0.00000;

	while?（fscanf（filein“%lf%lf“&point0.x&point0.y）!=EOF）
	//while?（fscanf（filein“%f%f“&xi&yi）!=EOF）
	{
		initialPoints.push_back（point0）;
		printf（“%f?%f\n“point0.xpoint0.y）;
	}
	fclose（filein）;
}
void?DouglasPeuker（int?leftpointint?rightpointdouble?tolerance）
{
	int?imaxindex;
	double?dismaxdis;
	maxdis=0;
	maxindex=0;
	for?（i=leftpoint;i	{
		dis=cacuDistance（leftpointrightpointi）;
		if?（dis>maxdis）
		{
			maxdis=dis;
			maxindex=i;
		}
	}

	if?（maxdis>tolerance）
	{
		bTag[maxindex]=1;
	//	p=initialPoints（maxindex）;
???????//?FinalPoints.push_back（initialPoints[maxindex]）;
		DouglasPeuker（leftpointmaxindextolerance）;
		DouglasPeuker（maxindexrightpointtolerance）;//回調(diào)
	}
	else
	{
	//	FinalPoints.push_back（initialPoints[leftpoint]）;
	//	FinalPoints.push_back（initialPoints[rightpoint]）;
	}

}

double?cacuDistance（int?lint?rint?n）
{
	double?abc;
	//POINT?p1p2p;
	double?x1x2y1y2xy;
	double?dis=0;

	
	y1=initialPoints[l].y*1000;
	x1=initialPoints[l].x*1000;
	x2=initialPoints[r].x*1000;
	y2=initialPoints[r].y*1000;

	x=initialPoints[n].x*1000;
	y=initialPoints[n].y*1000;

	if?（x1==x2）
	{
		dis=abs（x-x1）;//兩點(diǎn)式失效的情況
		r

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----

?????文件????????239??2011-08-25?21:56??DouglasPeuker\1.txt

?????文件???????3859??2011-08-26?09:52??DouglasPeuker\2.txt

?????文件?????200770??2011-08-26?09:52??DouglasPeuker\Debug\DouglasPeuker.exe

?????文件?????276840??2011-08-26?09:52??DouglasPeuker\Debug\DouglasPeuker.ilk

?????文件??????31825??2011-08-26?09:52??DouglasPeuker\Debug\DouglasPeuker.obj

?????文件?????187000??2011-08-25?21:23??DouglasPeuker\Debug\DouglasPeuker.pch

?????文件?????508928??2011-08-26?09:52??DouglasPeuker\Debug\DouglasPeuker.pdb

?????文件???????1656??2011-08-25?21:23??DouglasPeuker\Debug\StdAfx.obj

?????文件??????66560??2011-08-26?09:52??DouglasPeuker\Debug\vc60.idb

?????文件??????94208??2011-08-26?09:52??DouglasPeuker\Debug\vc60.pdb

?????文件???????3484??2011-08-26?09:52??DouglasPeuker\DouglasPeuker.cpp

?????文件???????4620??2011-08-25?20:57??DouglasPeuker\DouglasPeuker.dsp

?????文件????????549??2011-08-25?20:57??DouglasPeuker\DouglasPeuker.dsw

?????文件??????41984??2011-11-09?11:44??DouglasPeuker\DouglasPeuker.ncb

?????文件??????53760??2011-11-09?11:44??DouglasPeuker\DouglasPeuker.opt

?????文件???????1352??2011-08-26?09:52??DouglasPeuker\DouglasPeuker.plg

?????文件???????9360??2011-08-25?12:22??DouglasPeuker\jpl_vegetation_grass.txt

?????文件???????9289??2011-08-25?21:32??DouglasPeuker\jpl_vegetation_grass2.txt

?????文件???????1250??2011-08-25?20:57??DouglasPeuker\ReadMe.txt

?????文件????????300??2011-08-25?20:57??DouglasPeuker\StdAfx.cpp

?????文件????????667??2011-08-25?20:57??DouglasPeuker\StdAfx.h

?????目錄??????????0??2011-08-26?09:52??DouglasPeuker\Debug

?????目錄??????????0??2011-11-09?11:44??DouglasPeuker

-----------?---------??----------?-----??----

??????????????1498500????????????????????23