工控编程吧

标题: 上位机MFC凸壳生成源代码 [打印本页]

---

作者: qq263946146    时间: 2019-9-16 12:02
标题: 上位机MFC凸壳生成源代码
凸壳可以看作是点集合的边界。
其精确定义如下：
设集合S是n维空间的k个点组成的集合，即S={x1,x2,...xk},xi是n维向量。
定义S的凸壳Conv(S)为：
Conv(S)={x=λ1*x1+λ2*x2+...+λk*xk | λ1+λ2+ . . .+λk=1}
VC6.0编译运行，实例运行后，打开文件就可以生成凸壳。可以下载学习 (, 下载次数: 0)

上传

点击文件名下载附件
这个例程实现凸壳生成功能。

1. /*         for(j=0;j<pts.GetSize();j++)
   
2.          {
   
3.                  for(k=0;k<final_point.GetSize();k++)
   
4.                  {
   
5.                          change_pts.Add( (CPoint)final_point.GetAt(k) );
   
6.                  }
   
7.             
   
8.                  for(p=0,k=0;k<final_point.GetSize();k++)
   
9.          {
   
10.                          if( pts.GetAt(j)==final_point.GetAt(k) )   p=1;
    
11.                  }
    
12.                  
    
13.                  if(p==0)   
    
14.                  {
    
15.                          change_pts.Add( (CPoint)pts.GetAt(j) );
    
16.                  }
    
17.          }
    
18.    
    
19. */
    
20. 
    
21. 
    
22. 
    
23.            pti=(CPoint)change_pts.GetAt(i);
    
24.      ptii=(CPoint)change_pts.GetAt(i+1);
    
25.      pt_ii=(CPoint)change_pts.GetAt(i-1);
    
26. 
    
27.        w=( ptii.x-pti.x ) * ( pti.x-pt_ii.x )+( ptii.y-pti.y ) * ( pti.y-pt_ii.y );
    
28.        q1x=(ptii.x-pti.x)*(ptii.x-pti.x);          q1y=(ptii.y-pti.y)*(ptii.y-pti.y);
    
29.        q2x=(pti.x-pt_ii.x)*(pti.x-pt_ii.x);        q2y=(pti.y-pt_ii.y)*(pti.y-pt_ii.y);
    
30.        q1cos=w/ ( sqrt(q1x+q1y)*sqrt(q2x+q2y) );
    
31.       
    
32.       
    
33.            for(j=i+1,p=i+1; j<change_pts.GetSize();j++)
    
34.               {   pt_j_i=(CPoint)change_pts.GetAt(j);
    
35.                pti=(CPoint)change_pts.GetAt(i);
    
36.              
    
37.               w=(pt_j_i.x-pti.x) * (pti.x-pt_ii.x)+(pt_j_i.y-pti.y) * (pti.y-pt_ii.y);
    
38.           q1x=(pt_j_i.x-pti.x)*(pt_j_i.x-pti.x);     
    
39.                   q1y=(pt_j_i.y-pti.y)*(pt_j_i.y-pti.y);
    
40.         
    
41.                   q1cos=w/ ( sqrt(q1x+q1y)*sqrt(q2x+q2y) );
    
42.       
    
43.               if(  (q1cos!=1)&&((q1cos-q0cos)>0)  )
    
44.                   {
    
45.             q0cos=q1cos;   p=j;
    
46.                   }
    
47. 
    
48.            }
    
49.    final_point.Add( (CPoint)change_pts.GetAt(p) );
    
50.   point_change=(CPoint)final_point.GetAt( final_point.GetSize()-1 );
    
51.   i++;

复制代码

---

欢迎光临 工控编程吧 (https://www.gkbc8.com/)

Powered by Discuz! X3.4