delauney三角形剖分程序文档格式.docx

1、 static long z=17; x=（171*x）%30269; y=（172*y）%30307; z=（170*z）%30323; long x1=x; long y1=y; long z1=z; double T,R; T=double（x1）/double（30269）+double（y1）/double（30307）+double（z1）/double（30323）; R=T-long（T）; return R;angle（T_line line,long p） double c; double a; double b; c=distance（line.p0,line.p1）;

2、a=distance（line.p0,p）; b=distance（line.p1,p）; double z=0; z=acos（b*b+a*a-c*c）/（2*a*b）; return z;long T_TIN:Max（T_line line,vector usableList） long index=-1; double xitaMax=0; double xita=0; for（long i=0;iusableList.size（）;i+） xita=angle（line,usableListi）; if（xitaMaxSelectObject（brush）; / pDC-SetPixe

3、l（pp,RGB（255,0,0）;Ellipse（&pp）;createTin（） vector usalbePlointList; /生成第0个三角形 T_tri tri; long k=0; /第L0边 T_line line0; line0.p0=0; line0.p1=1; line0.useCount=1; lineList.push_back（line0）; tri.L0=lineList.size（）-1; / for（long i=2;i+）/生成可用点的连表 usalbePlointList.push_back（i）; index=Max（line0,usalbePloin

4、tList）; /生成 L1，L2边 T_line line1; T_line line2; line1.p0=0; line1.p1=index; line1.useCount=1; lineList.push_back（line1）; tri.L1=lineList.size（）-1; line2.p0=1; line2.p1=index; line2.useCount=1; lineList.push_back（line2）; tri.L2=lineList.size（）-1; triList.push_back（tri）; /CString s; /s.Format（连接顶点0,1,%

5、d,index）; /AfxMessageBox（s）; /* 对于第k个三角形为拓展三角形时 */ while（1） /*第L0边*/ if（lineListtriListk.L0.useCount2）/是可拓展边 /找拓展顶点 long tuo=-1; if（lineListtriListk.L0.p0=lineListtriListk.L1.p0 |lineListtriListk.L0.p1=lineListtriListk.L1.p0） tuo=lineListtriListk.L1.p1; if（lineListtriListk.L0.p0=lineListtriListk.L1.

6、p1 |lineListtriListk.L0.p1=lineListtriListk.L1.p1） tuo=lineListtriListk.L1.p0; usalbePlointList.clear（）; double Fuhao=0; double Fu=0; Fuhao=F（lineListtriListk.L0,tuo）; for（long i=0;i+）/找到可用的顶点 if（ lineListtriListk.L0.p0!=i /不是三角形的顶点 &lineListtriListk.L0.p1!=ilineListtriListk.L1.p0!lineListtriListk.L

7、1.p1!lineListtriListk.L2.p0!lineListtriListk.L2.p1!=i） Fu=F（lineListtriListk.L0,i）; if（Fuhao*Fu0） int idx=-1; idx=Max（lineListtriListk.L0,usalbePlointList）; long lineIdx1=-1; long lineIdx2=-1; int cout1=0; int cout2=0; T_line line1; T_line line2; line1.p0=lineListtriListk.L0.p0; line1.p1=idx; line2.

8、p0=lineListtriListk.L0.p1; line2.p1=idx; for（long i=0;lineList.size（）; if（line1=lineListi） lineIdx1=i; cout1=lineListlineIdx1.useCount; if（line2=lineListi） lineIdx2=i; cout2=lineListlineIdx2.useCount; if（cout12&cout22） T_tri tri; tri.L0=triListk.L0; if（cout1=0） T_line line; line.p0=lineListtriListk.

9、L0.p0; line.p1=idx; line.useCount=1; lineList.push_back（line）; tri.L1=lineList.size（）-1; if（cout1=1） lineListlineIdx1.useCount=2; tri.L1=lineIdx1; if（cout2=0） line.p0=lineListtriListk.L0.p1; tri.L2=lineList.size（）-1; if（cout2=1） lineListlineIdx2.useCount=2; tri.L2=lineIdx2; triList.push_back（tri）; l

10、ineListtriListk.L0.useCount=2; /*第L1边*/ if（lineListtriListk.L1.useCount if（lineListtriListk.L1.p0=lineListtriListk.L0.p0 |lineListtriListk.L1.p1=lineListtriListk.L0.p0） tuo=lineListtriListk.L0.p1; if（lineListtriListk.L1.p0=lineListtriListk.L0.p1 |lineListtriListk.L1.p1=lineListtriListk.L0.p1） tuo=li

11、neListtriListk.L0.p0; Fuhao=F（lineListtriListk.L1,tuo）; Fu=F（lineListtriListk.L1,i）; idx=Max（lineListtriListk.L1,usalbePlointList）; line1.p0=lineListtriListk.L1.p0; line2.p0=lineListtriListk.L1.p1; tri.L0=triListk.L1; line.p0=lineListtriListk.L1.p0; line.p0=lineListtriListk.L1.p1; lineListtriListk.L

12、1.useCount=2; /*第L2边*/ if（lineListtriListk.L2.useCount if（lineListtriListk.L2.p0=lineListtriListk.L0.p0 |lineListtriListk.L2.p1=lineListtriListk.L0.p0） if（lineListtriListk.L2.p0=lineListtriListk.L0.p1 |lineListtriListk.L2.p1=lineListtriListk.L0.p1） Fuhao=F（lineListtriListk.L2,tuo）; Fu=F（lineListtriListk.L2,i）; idx=Max（lineListtriListk.L2,usalbePlointList）; line1.p0=lineListtriListk.L2.p0; line2.p0=lineListtriListk.L2.p1; tri.L0=triListk.L2; line.p0=lineListtriListk.L2.p0; if（