用多条直线去描述一条曲线

1. #_*_coding:utf-8
   
2. """
   
3. 制作人：盖房子的猫熊
   
4. 用多条直线去描述一条曲线,采用的算法为二分法，和递归：
   
5. accuracy 为直线的精确度。
   
6. """
   
7. 
   
8. import rhinoscriptsyntax as rs
   
9. 
   
10. 
    
11. def curve_line(curve,accuracy):
    
12.     #定义了一个方法去判定曲线是否需要优化，利用的数学原理为两点之间直线最短
    
13.     curve_length=rs.CurveLength(curve)
    
14.     end_point=rs.CurveEndPoint(curve)
    
15.     star_point=rs.CurveStartPoint(curve)
    
16.     line=rs.AddLine(end_point,star_point)
    
17.     line_length=rs.CurveLength(line)
    
18.     if curve_length == 0:     #为防止曲线长度为零报错，提前捕获，也可采用try 语法。
    
19.         line_curve = 1
    
20.     else:
    
21.         line_curve=line_length/curve_length
    
22.     if line_curve>accuracy:
    
23.         return False          #不需要优化返回False
    
24.     else:
    
25.         return True          #需要优化返回True
    
26.    
    
27. 
    
28. def divide_curve(L,accuracy,line_list):
    
29.     #定义了一个递归方法。L为递归的曲线，accuracy为精确度，line_list为返回的直线
    
30. 
    
31.     L_curve=[]  #定义一个空的列表，用于存放等分的曲线
    
32.     bool=True  #设定一个布尔值，用于判定是否返回直线，还是进入下一次递归
    
33.     for i in range(len(L)):
    
34.         if curve_line(L[i],accuracy):
    
35.             bool = False
    
36.             if rs.IsCurve(L[i]):   
    
37.                 domain = rs.CurveDomain(L[i])
    
38.                 parameter = (domain[1]-domain[0]) / 2.0+domain[0]
    
39.                 #求出曲线的中点，此处有陷阱，不可直接使用domain[1]/2
    
40.                 split_curve=rs.SplitCurve( L[i],parameter)
    
41.                 L_curve.extend(split_curve) #列表方法将一个列表的内容追加给前一个列表
    
42.                
    
43.         else:
    
44.             end_point=rs.CurveEndPoint(L[i])
    
45.             star_point=rs.CurveStartPoint(L[i])
    
46.             line=rs.AddLine(end_point,star_point)
    
47.             line_list.append(line)  #求出精确描述的每一段直线
    
48.             
    
49.     if bool == True:
    
50.        return line_list
    
51.     else:
    
52.         return divide_curve(L_curve,accuracy,line_list)
    
53.         
    
54.         
    
55. if __name__ == "__main__":
    
56.     if accuracy>0.999:
    
57.         accuracy=0.999
    
58.     elif accuracy is None:
    
59.         accuracy=0.9
    
60.     line_list=[]
    
61.     lines=divide_curve(curve,accuracy,line_list)

复制代码

猫熊 发表于 2016-3-10 16:10
你看，是吧.

找到原因了，问题已修正，感谢反馈。