分享
 
 
 

优化了的过关键点的光滑曲线拟合算法

王朝vc·作者佚名  2006-01-08
窄屏简体版  字體: |||超大  

这个是我一个数学老师(教授,数学高手,经常自己做算法)给我的例子,用于多个离散点拟合光滑曲线的,他优化了追赶法,这个例子适用于闭合和不闭合两种情况。当时由于工程情况,写的急,代码不好看,但是很好用。为了方便传递参数,我做了一个链表,用时候根据自己情况可以修改,核心算法不动即可。

class CFoldPoint

{public:

double X; double Y;

};

typedef CTypedPtrList CFoldPointList;

typedef CArray CDoubleArray;

三个函数,SPLine 调用另外两个。用时候直接调用SPLine函数,入口pList是已知离散点链表,pDestList是生成的点的链表。SM是在两个点中间插入点的数目,continue=0是采样点无规律,要求生成闭合曲线。1是采样点x坐标连续 2是y连续

void ZG(CDoubleArray *A,CDoubleArray *B,CDoubleArray *C,CDoubleArray *G,int &LOGI)

{

//追赶法

register long I;

int N;

N=A-GetSize();

if(LOGI==0)

{

(*C)[0]=(*C)[0]/(*B)[0];

for(I=1;I

{

(*B)[I]=(*B)[I]-(*A)[I]*(*C)[I-1];

(*C)[I]=(*C)[I]/(*B)[I];

}

(*A)[0]=0.;

(*C)[N-1]=0.;

LOGI=1;

}

(*G)[0]=(*G)[0]/(*B)[0];

for(I=1;I

{

(*G)[I]=((*G)[I]-(*A)[I]*(*G)[I-1])/(*B)[I];

}

for(I=N-2;I-1;I--)//DO 30 I=N-1,1,-1

{

(*G)[I]=(*G)[I]-(*C)[I]*(*G)[I+1];

}

return;

}

void SPLine4(CDoubleArray *X,CDoubleArray *Y,double &XI,double&YI,CDoubleArray *A,CDoubleArray *B,CDoubleArray *C,CDoubleArray *G,int &LOGI,int MD)

{

register long I;

double W1,W2,H;

int N=X-GetSize();

if(LOGI==0)

{

for(I=1;I

{

(*B)[I]=(*X)[I]-(*X)[I-1];

(*C)[I]=((*Y)[I]-(*Y)[I-1])/(*B)[I];

}

for(I=1;I

{

(*A)[I]=(*B)[I]+(*B)[I+1];

(*G)[I]=6.*((*C)[I+1]-(*C)[I])/(*A)[I];

(*A)[I]=(*B)[I]/(*A)[I];

}

for(I=1;I

{

(*C)[I]=1.-(*A)[I];

(*B)[I]=2.;

}

(*B)[0]=2.;

(*B)[N-1]=2.;

if(MD==3)

{

(*C)[0]=-1.;

(*A)[N-1]=-1.;

(*A)[0]=0.;

(*C)[N-1]=0.;

}

ZG(A,B,C,G,LOGI);

}

for(I=1;I

{

if(XI=(*X)[I-1] && XI

{

H=(*X)[I]-(*X)[I-1];

W1=(*X)[I]-XI;

W2=XI-(*X)[I-1];

YI=W1*W1*W1*(*G)[I-1]/6./H;

YI=YI+W2*W2*W2*(*G)[I]/6./H;

YI=YI+W1*((*Y)[I-1]-(*G)[I-1]*H*H/6.)/H;

YI=YI+W2*((*Y)[I]-(*G)[I]*H*H/6.)/H;

}

}

}

void SPLine(CFoldPointList *pList,CFoldPointList *pDestList,int SM,int Continue=0)

{

CFoldPoint *pFoldHead,*pFoldTail;

POSITION pos;

CDoubleArray A,B,C,G,X,Y,T;

double XI,YI,XX,YY;

register long i;

long N;

int LOGI;

long RealSM;

long Bei,Yu;

CFoldPoint *pFold;

file://赋初值

N=pList-GetCount();

A.SetSize(N);

B.SetSize(N);

C.SetSize(N);

G.SetSize(N);

X.SetSize(N);

Y.SetSize(N);

T.SetSize(N);

RealSM=(N-1)*SM+N;

pos=pList-GetHeadPosition();

for(i=0;i

{

pFold=pList-GetNext(pos);

X[i]=pFold-X;

Y[i]=pFold-Y;

}

pFoldHead=pList-GetHead();

pFoldTail=pList-GetTail();

if(Continue==0)//pFoldHead-X==pFoldTail-X && pFoldHead-Y==pFoldTail-Y)

{ file://闭

T[0]=0;

for(i=0;i

{

T[i+1]=T[i]+CalculateDistance(X[i],Y[i],X[i+1],Y[i+1])+0.000000001;

}

LOGI=0;

YI=0;

for(i=0;i

{

Bei=i/(SM+1);

Yu=i%(SM+1);

if(Yu!=0)

{

XI=T[Bei]+(T[Bei+1]-T[Bei])/(SM+1)*Yu;

SPLine4(&T,&Y,XI,YI,&A,&B,&C,&G,LOGI,3);

YY=YI;//+Y[Bei];

}

else

{

YY=Y[Bei];

}

pFold=new CFoldPoint;

pFold-Y=YY;

pDestList-AddTail(pFold);

}

LOGI=0;

YI=0;

pos=pDestList-GetHeadPosition();

for(i=0;i

{

Bei=i/(SM+1);

Yu=i%(SM+1);

if(Yu!=0)

{

XI=T[Bei]+(T[Bei+1]-T[Bei])/(SM+1)*Yu;

SPLine4(&T,&X,XI,YI,&A,&B,&C,&G,LOGI,3);

YY=YI;//+X[Bei];

}

else

{

YY=X[Bei];

}

pFold=pDestList-GetNext(pos);

pFold-X=YY;

}

}

else if(Continue==1)

{

file://x连续

LOGI=0;

YI=0;

for(i=0;i

{

Bei=i/(SM+1);

Yu=i%(SM+1);

if(Yu!=0)

{

XI=X[Bei]+(X[Bei+1]-X[Bei])/(SM+1)*Yu;

SPLine4(&X,&Y,XI,YI,&A,&B,&C,&G,LOGI,3);

XX=XI;

YY=YI;

}

else

{

XX=X[Bei];

YY=Y[Bei];

}

pFold=new CFoldPoint;

pFold-X=XX;

pFold-Y=YY;

pDestList-AddTail(pFold);

}

}

else

{

file://y连续

LOGI=0;

YI=0;

for(i=0;i

{

Bei=i/(SM+1);

Yu=i%(SM+1);

if(Yu!=0)

{

XI=Y[Bei]+(Y[Bei+1]-Y[Bei])/(SM+1)*Yu;

SPLine4(&Y,&X,XI,YI,&A,&B,&C,&G,LOGI,3);

XX=YI;

YY=XI;

}

else

{

XX=X[Bei];

YY=Y[Bei];

}

pFold=new CFoldPoint;

pFold-X=XX;

pFold-Y=YY;

pDestList-AddTail(pFold);

}

}

return;

}

 
 
 
免责声明:本文为网络用户发布,其观点仅代表作者个人观点,与本站无关,本站仅提供信息存储服务。文中陈述内容未经本站证实,其真实性、完整性、及时性本站不作任何保证或承诺,请读者仅作参考,并请自行核实相关内容。
2023年上半年GDP全球前十五强
 百态   2023-10-24
美众议院议长启动对拜登的弹劾调查
 百态   2023-09-13
上海、济南、武汉等多地出现不明坠落物
 探索   2023-09-06
印度或要将国名改为“巴拉特”
 百态   2023-09-06
男子为女友送行,买票不登机被捕
 百态   2023-08-20
手机地震预警功能怎么开?
 干货   2023-08-06
女子4年卖2套房花700多万做美容:不但没变美脸,面部还出现变形
 百态   2023-08-04
住户一楼被水淹 还冲来8头猪
 百态   2023-07-31
女子体内爬出大量瓜子状活虫
 百态   2023-07-25
地球连续35年收到神秘规律性信号,网友:不要回答!
 探索   2023-07-21
全球镓价格本周大涨27%
 探索   2023-07-09
钱都流向了那些不缺钱的人,苦都留给了能吃苦的人
 探索   2023-07-02
倩女手游刀客魅者强控制(强混乱强眩晕强睡眠)和对应控制抗性的关系
 百态   2020-08-20
美国5月9日最新疫情:美国确诊人数突破131万
 百态   2020-05-09
荷兰政府宣布将集体辞职
 干货   2020-04-30
倩女幽魂手游师徒任务情义春秋猜成语答案逍遥观:鹏程万里
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案神机营:射石饮羽
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案昆仑山:拔刀相助
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案天工阁:鬼斧神工
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案丝路古道:单枪匹马
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案镇郊荒野:与虎谋皮
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案镇郊荒野:李代桃僵
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案镇郊荒野:指鹿为马
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案金陵:小鸟依人
 干货   2019-11-12
倩女幽魂手游师徒任务情义春秋猜成语答案金陵:千金买邻
 干货   2019-11-12
 
推荐阅读
 
 
 
>>返回首頁<<
 
靜靜地坐在廢墟上,四周的荒凉一望無際,忽然覺得,淒涼也很美
© 2005- 王朝網路 版權所有