我如何在粘性中摆脱球体表面上的锯齿ZAG伪像？

Question

我正在研究太阳系中行星的3D可视化。当我要应用纹理时，我已经手动计算了纹理坐标（texcoords），但是，我得到图像中的锯齿形伪像。

我相信我的计算可能是错误的。我已经附上了texcoords计算：

# Compute Texture Coordinates  
def get_texcoords(vertices):
    texcoords = []
    for v in vertices:
        #thresholding
        for i in range(3):
            if np.abs(v[i]) > 1e-6:
                v[i] = v[i]
            elif np.abs(v[i]) < 1e-6:
                v[i] = 0.0

        # Compute position in uv-space
        radius = np.sqrt(v[0]**2 + v[1]**2 + v[2]**2)
        
        latitude = np.arcsin(v[2]/radius)
        longitude = np.arctan2(v[1],v[0])
            
        # Convert to texture coordinates
        u = round(0.5 + longitude/(2*np.pi),5)
        v = round(0.5 + latitude/np.pi,5)

        texcoords.append([u,v])
            
    return np.array(texcoords)

有什么方法可以删除工件？还是在vispy中获得纹理坐标的智能方法？

Answer 1

也许使用“在顶点中的V for V：”与使用纹理脐带变量名称为“ u”和“ v”？我不确定它会产生该工件，但是绝对不建议在不谨慎的情况下修改迭代器内部的迭代器。实际上不建议以任何方式推荐。那是我的最初看法。我正在使用同样的算法，所以一旦我开始工作，我会知道这是否是您的问题的原因。

我刚刚找到了此链接，该链接解释了紫外线映射算法中潜伏的错误： https：https：// mft：// mft -dev.dk/uv-mapping-sphere/

简而言之，一些网状三角形越过纹理边界，其中一些最终被“翻转”（即他们的脸正常点向内而不是向外）） 。链接中的文章显示了如何找到这些流氓面孔并修复它们。

我一直在挠头，一天左右。

Answer 2

这是一个代码段，它将在没有Zig-Zag伪像的情况下生成适当的网格：

代码来自： https://www.songho.ca/opengl/index.html

colstep = 2 * pi / cols
    rowstep = pi / rows
    num_v = -1  # these counters are for logging/debuging
    num_f = -1
    num_e = -1
    num_eh = -1
    num_ev = -1

    for row in range(0, rows + 1):

        phi = pi / 2 - row * rowstep
        xy = radius[1] * cos(phi)
        z = radius[2] * sin(phi)

        for col in range(0, cols + 1):
            theta = col * colstep
            x = xy * cos(theta)
            y = xy * sin(theta)
            vert = np.array([x, y, z])
            # print(vert)
            self._verts.append(vert)
            self._norms.append(vert / np.sqrt(vert.dot(vert)))
            self._txcds.append(np.array([(col / cols), (row / rows)]))
            num_v += 1
    logging.info('----->>> Generated %r vertices...', num_v)

    for i in range(0, rows):

        k1 = i * (cols + 1)
        k2 = k1 + cols + 1

        for j in range(0, cols):

            if i != 0:
                self._faces.append(np.array([k1, k2, k1 + 1]))
                self._edges.append(np.array([k1, k2]))
                self._edges.append(np.array([k2, k1 + 1]))
                self._edges.append(np.array([k1 + 1, k1]))
                self._edge_colors.append((0, 0, 0, 1))
                self._edge_colors.append((0, 0, 0, 0))
                self._edge_colors.append((0, 0, 0, 1))
                num_f += 1
                num_e += 3
            if i != (rows - 1):
                self._faces.append(np.array([k1 + 1, k2, k2 + 1]))
                self._edges.append(np.array([k1 + 1, k2]))
                self._edges.append(np.array([k2, k2 + 1]))
                self._edges.append(np.array([k2 + 1, k1 + 1]))
                self._edge_colors.append((0, 0, 0, 0))
                self._edge_colors.append((0, 0, 0, 0))
                self._edge_colors.append((0, 0, 0, 0))
                num_f += 1
                num_e += 3

            k1 += 1
            k2 += 1

            self._v_edges.append(np.array([k1, k2]))
            num_ev += 1
            if i != 0:
                self._h_edges.append(np.array([k1, k1 + 1]))
                num_eh += 1