如何找一个点附近的四个点？

hugegraywolf

假设在一个点周围随机分布着一些点且坐标都是已知的，如何快速定位它附近的四个点？这里，“附近”的涵义是以该点为原点，将平面划分为四个象限，分别从这四个象限中找一个与该点距离最近的点出来。

xiezhh

我编了一个函数，代码如下：

1. function  Point = FindPoint(xy0, xy)
   
2. % xy0为多行2列的矩阵，xy为包含2个元素的向量
   
3. 
   
4. Point = cell(1,4);
   
5. x = xy0(:,1) - xy(1);
   
6. y = xy0(:,2) - xy(2);
   
7. % 第一象限
   
8. id1 = (x>0 & y>0);
   
9. Point{1} = subfun(id1,x,y,xy0);
   
10. % 第二象限
    
11. id2 = (x<0 & y>0);
    
12. Point{2} = subfun(id2,x,y,xy0);
    
13. % 第三象限
    
14. id3 = (x<0 & y<0);
    
15. Point{3} = subfun(id3,x,y,xy0);
    
16. % 第四象限
    
17. id4 = (x>0 & y<0);
    
18. Point{4} = subfun(id4,x,y,xy0);
    
19. 
    
20. % 子函数
    
21. function xyPoint = subfun(id,x,y,xy0)
    
22. xyPoint = [];
    
23. x1 = x(id);
    
24. if ~isempty(x1)
    
25.     y1 = y(id);
    
26.     distance = x1.^2 + y1.^2;
    
27.     xy0_1 = xy0(id,:);
    
28.     xyPoint = xy0_1(find(distance == min(distance)),:);
    
29. end

复制代码

用法：

1. >> xy0 = normrnd(0,1,100,2);
   
2. >> xy = [0 0];
   
3. >> p = FindPoint(xy0,xy)
   
4. 
   
5. p =
   
6. 
   
7.     [1x2 double]    [1x2 double]    [1x2 double]    [1x2 double]
   
8. 
   
9. >> p{:}
   
10. 
    
11. ans =
    
12. 
    
13.     0.0689    0.0327
    
14. 
    
15. 
    
16. ans =
    
17. 
    
18.    -0.3368    0.2209
    
19. 
    
20. 
    
21. ans =
    
22. 
    
23.    -0.0896   -0.2915
    
24. 
    
25. 
    
26. ans =
    
27. 
    
28.     0.0618   -0.0210

复制代码

[ 本帖最后由 xiezhh 于 2010-1-29 14:40 编辑 ]