豆粑粑，matlab 找两个散点的相交集合，intersection of two data sets

meatball1982

别人的问题，
两个点集都有不少点。
要找这两个点集相交的集合。

By the way, I will be dead soon. So if u have problem about this code, contact me please.
meatball1982@163.com.

1. clear all
   
2. clc
   
3. clf
   
4. 
   
5. load('../data/data1.mat');
   
6. load('../data/data2.mat');
   
7. a=M6ofR2;
   
8. b=M6ofR3;
   
9. 
   
10. clear M6ofR2 M6ofR3
    
11. 
    
12. % choose part of the points
    
13. ind1 = a(:,1) > -6 & a(:,1) < -2 & a(:,2) > 2 & a(:,2) < 5 ;
    
14. ind2 = b(:,1) > -6 & b(:,1) < -2 & b(:,2) > 2 & b(:,2) < 5 ;
    
15. p1= a(ind1,:);
    
16. p2= b(ind2,:);
    
17. 
    
18. % the multi distance of two point sets
    
19. D = pdist2(p1,p2);
    
20. 
    
21. % chose intersect part (with radius 0.1)
    
22. indsmall1 = min(D)<0.1;
    
23. p2(~indsmall1,:)=[];
    
24. 
    
25. indsmall2 = min(D')<0.1;
    
26. p1(~indsmall2,:)=[];
    
27. 
    
28. %
    
29. points1=p1;
    
30. points2=p2;
    
31. 
    
32. %
    
33. [chull1,cf1,df1]=convhull_nd(points1);
    
34. inconvhull1=~any(cf1*points2'+df1(:,ones(1,size(points2,1)))>0,1);
    
35. inter_points1=points2(inconvhull1,:);
    
36. 
    
37. [chull2,cf2,df2]=convhull_nd(points2);
    
38. inconvhull2=~any(cf2*points1'+df2(:,ones(1,size(points1,1)))>0,1);
    
39. inter_points2=points1(inconvhull2,:);
    
40. 
    
41. % inter points
    
42. inter_points=[inter_points1;inter_points2];
    
43. 
    
44. % the boundary of inter points
    
45. k_all=boundary(inter_points,0.8);
    
46. x_bou = inter_points(k_all,1);
    
47. y_bou = inter_points(k_all,2);
    
48. 
    
49. h=figure(1)
    
50. set(h, 'Position', [1000, 100, 900, 800]);
    
51. hold on
    
52. plot(a(:,1),a(:,2),'c.')
    
53. plot(b(:,1),b(:,2),'m.')
    
54. 
    
55. plot(points1(:,1),points1(:,2),'o','MarkerEdgeColor','b','LineWidth',2)
    
56. plot(points2(:,1),points2(:,2),'o','MarkerEdgeColor','r','LineWidth',2)
    
57. 
    
58. plot(inter_points(:,1),inter_points(:,2),'s','MarkerEdgeColor','g','LineWidth',1,'markersize',10);
    
59. plot(x_bou,y_bou,'k-','linewidth',3)
    
60. 
    
61. leg_str={'dat 1','dat 2','int part in dat 1','int part in dat 2',' int points','boundary points'};
    
62. legend(leg_str,'location','northwest')

复制代码

需要几个函数。