计算两多边形是否香蕉的代码……

yangff

1. module CV
   
2.         module_function
   
3.         class Point
   
4.                 def initialize(x,y)
   
5.                         @x=x.to_f
   
6.                         @y=y.to_f
   
7.                 end
   
8.                 def x
   
9.                         return @x
   
10.                 end
    
11.                 def y
    
12.                         return @y
    
13.                 end
    
14.                 def x=(v)
    
15.                         return (@x=v)
    
16.                 end
    
17.                 def y=(v)
    
18.                         return (@y=v)
    
19.                 end
    
20. 
    
21.                 def in_polygon(p)
    
22.                         p=p.points
    
23.                         nCross=0
    
24.                         for i in 0...p.size
    
25.                                 p1=p[i]
    
26.                                 p2=p[(i+1) % p.size]
    
27.                                 next if p1.y == p2.y
    
28.                                 next if y < [p1.y, p2.y].min
    
29.                                 next if y >= [p1.y, p2.y].max
    
30.                                
    
31.                                 x1=(y - p1.y)*(p2.x - p1.x)/(p2.y - p1.y) + p1.x
    
32.                                 #p  (p2.y - p1.y)
    
33.                                 nCross+=1 if x1>x
    
34.                         end
    
35.                         #p nCross
    
36.                         return (nCross % 2 == 1)
    
37.                 end
    
38.         end
    
39.         class Grap
    
40.        
    
41.         end
    
42.         class Polygon<Grap
    
43.                 def initialize(*args)
    
44.                         @points=args
    
45.                 end
    
46.                 def points
    
47.                         return @points
    
48.                 end
    
49.                 def &(v)
    
50.                         if v.is_a?(Grap)
    
51.                                 if v.is_a?(Polygon)
    
52.                                         for i in @points
    
53.                                                 if i.in_polygon(v)==true
    
54.                                                         return true
    
55.                                                 end
    
56.                                         end
    
57.                                         for i in v.points
    
58.                                                 if i.in_polygon(self)==true
    
59.                                                         return true
    
60.                                                 end
    
61.                                         end
    
62.                                 end
    
63.                         end
    
64.                         return false
    
65.                 end
    
66.         end
    
67.         class Vector
    
68.                
    
69.         end
    
70. 
    
71. end
    
72. poy=CV::Polygon.new(CV::Point.new(0,0),CV::Point.new(1,2),CV::Point.new(5,5),CV::Point.new(5,15))
    
73. poy1=CV::Polygon.new(CV::Point.new(300,500),CV::Point.new(100,200),CV::Point.new(500,500),CV::Point.new(500,1500))
    
74. p poy&poy1

复制代码

效率一般……没做粗判定……