EllipseFit.py 文件源码

def elltest(scale=0.8,off=0.2):
    #generate an example, random, non-self-intersecting polygon.
    #This is done by first generating
    #it in polar coordinates and than translating it
    #to cartesian.
    Theta1,R1=linspace(0,2*pi,30),rand(30)*scale+off
    X1,Y1=R1*cos(Theta1),R1*sin(Theta1)
    X1=append(X1,X1[0])
    Y1=append(Y1,Y1[0])

    p.plot(X1,Y1,".-",ms=10)


    a2,b2,ecc2,alpha2=ellfit(X1,Y1,showFig=False)

    Xe,Ye=ellipse(b2,a2,-alpha2,X1.mean(),Y1.mean(),Nb=40)

    p.plot(Xe,Ye,"r.-")


    p.grid(True)
    p.show()
    pass