Приближение сферы в OpenGL

Question

Я пытаюсь приблизить сферу, используя инструкции на http://local.wasp.uwa.edu.au/~pbourke/miscellaneous/sphere_cylinder/,, но она выглядит неправильно.Это мой код:

def draw_sphere(facets, radius=100):
    """approximate a sphere using a certain number of facets"""

    dtheta = 180.0 / facets
    dphi = 360.0 / facets

    global sphere_list
    sphere_list = glGenLists(2)
    glNewList(sphere_list, GL_COMPILE)

    glBegin(GL_QUADS)

    for theta in range(-90, 90, int(dtheta)):
        for phi in range(0, 360, int(dphi)):
            print theta, phi
            a1 = theta, phi
            a2 = theta + dtheta, phi
            a3 = theta + dtheta, phi + dphi
            a4 = theta, phi + dphi

            angles = [a1, a2, a3, a4]

            print 'angles: %s' % (angles)


            glColor4f(theta/360.,phi/360.,1,0.5)

            for angle in angles:
                x, y, z = angle_to_coords(angle[0], angle[1], radius)
                print 'coords: %s,%s,%s' % (x, y, z)
                glVertex3f(x, y, z)



    glEnd()

    glEndList()


def angle_to_coords(theta, phi, radius):  
    """return coordinates of point on sphere given angles and radius"""

    x = cos(theta) * cos(phi)
    y = cos(theta) * sin(phi)
    z = sin(theta)

    return x * radius, y * radius, z * radius

Кажется, что некоторые из квадратов не простые, то есть ребра пересекаются, но изменение порядка вершин, кажется, не имеет никакого значения.*

Answer 1

У меня нет системы, которая могла бы запускать Python и OpenGL вместе, но я все равно вижу некоторые проблемы:

Вы округляете dphi и dtheta в операторах range. Это означает, что фасеты всегда будут начинаться с целого градуса, но тогда, когда вы добавляете необоснованные значения дельты, дальний край не гарантирует этого. Это вероятная причина ваших совпадений.

Было бы лучше, чтобы значения вашего диапазона пошли от 0 .. facets-1, а затем умножили эти индексы на 360 / фасет (или 180 для линий широты) точно. Это позволит избежать ошибок округления, например ::10000

dtheta = 180.0 / facets
dphi = 360.0 / facets

for y in range(facets):
    theta = y * dtheta - 90
    for x in range(facets):
        phi = x * dphi
        ...

Кроме того, вы переводите в радианы в другом месте? Стандартные триггерные функции Python принимают радианы, а не градусы.