{"topic":{"id":13851,"title":"几何部分 loop 细分部分不正常","slug":"topic","categoryId":56,"createdAt":"2024-11-11 11:55:54.922617","bumpedAt":"2024-11-12 01:04:28.756401","views":183,"postsCount":2,"highestPostNumber":2,"likeCount":0,"closed":false,"pinned":false,"pinnedGlobally":false,"pinnedAt":null,"pinnedUntil":null,"isBanner":false,"excerpt":"我在几何部分第二章 loop 面细分的时候，在 cube 细分达不到文档的要求，面是平滑了，但是八个顶点突出。这是一次细分，可以看到八个顶点明显这是四次细分，八个顶点更为明显我觉得应该是新顶点计算有问题，我的实现如下： if (v->is_new) { Eigen::Vector3f V1AddV2 = Eigen::Vector3f::Zero(); Eigen::Vector3f V3AddV4 = Eigen::Vec…","deletedAt":null,"isAnonymous":false,"category":{"id":56,"name":"计算机图形学","slug":"cg","color":"92278F"},"author":{"id":6145,"username":"SunXueHao","name":"rss","avatarUrl":null},"tags":[]},"posts":[{"id":43171,"postNumber":1,"raw":"我在几何部分第二章 loop 面细分的时候，在 cube 细分达不到文档的要求，面是平滑了，但是八个顶点突出。\n![QQ 图片 20241111193432|635x500, 50%](/uploads/image/1b/0f/1b0f0e8a326d2292775650eb89cd06ffe899425200d546a7a187fa3b857cd830.avif)\n这是一次细分，可以看到八个顶点明显\n![QQ 图片 20241111193438|608x488, 50%](/uploads/image/3a/9d/3a9d78bc2aee3f14cbe06b10a4fba2fc1f5b8c740c12a210382d642e8f9d1381.avif)\n这是四次细分，八个顶点更为明显\n我觉得应该是新顶点计算有问题，我的实现如下：\n```\nif (v->is_new)\n {\n Eigen::Vector3f V1AddV2 = Eigen::Vector3f::Zero();\n Eigen::Vector3f V3AddV4 = Eigen::Vector3f::Zero();\n Eigen::Vector3f V1V2V3V4 = Eigen::Vector3f::Zero();\n Halfedge* h = v->halfedge;\n do\n {\n if(h->edge->is_new)\n {\n V1V2V3V4 += h->inv->from->pos;\n }\n h = h->inv->next;\n } while (h != v->halfedge);\n V1AddV2 = 2.0f * v->pos;\n V3AddV4 = V1V2V3V4 - V1AddV2;//不知道为啥这里要减去 V3AddV4 才第一次 cube 正常，之后就不正常，别的物体不正常。\n v->new_pos = 3.0f/8.0f * V1AddV2 + 1.0f/8.0f * V3AddV4;\n }\n```\n![QQ 图片 20241111195046|552x500](/uploads/image/ec/a8/eca8d848a3574f91118bc4cbc5115ea87d1e059d3441bd1844c9d73f0ab15738.avif)\n\n方法是在 split 生成的新点坐标正好是 1/2（V1+V2），就可以作为一个整体来计算 V1+V2 和 V3+V4，通过对一个新顶点周围顶点对 split 产生的新边遍历求坐标和（图中红色顶点的另外四个虚线的点的坐标和），正好是 V1+V2+V3+V4，V1+V2 是它中点新顶点坐标的两倍，然后这俩求差就可以得到 V3+V4;我觉得这么计算没有问题。但是遇到了怪异的现象：\n1.V3AddV4 = V1V2V3V4 - V1AddV2;我把后面减去的 V1+V2 换成 V3+V4（这个显然是错误的），cube 第一次细分和文档一摸一样，但是后续细分错误，其他的 cow 和兔子错误，尖锐和凹陷明显。\n2.V3AddV4 = V1V2V3V4 - V1AddV2;仍然用这个，cube 细分错误，但是兔子和 cow 细分正确，他们的表面仔细检查十分平滑，没有突起和凹陷。\n![QQ 图片 20241111193427|690x494, 50%](/uploads/image/3c/29/3c2997b72c6af458f2a2936799456eee9e7f1dfc1ca63cc230950ee2725f8113.avif)\n想知道是什么原因导致的部分正确，按理说要正确就全部正确吗。\n```\nif (n == 3)\n {\n v->new_pos = (1.0f - (9.0f / 16.0f)) * v->pos + (3.0f / 16.0f) * Eneighbors;\n }\n else\n {\n // For other vertices, apply the general rule based on degree\n v->new_pos = (1.0f - (3.0f / 8.0f)) * v->pos + (0.375f / n) * Eneighbors;\n }\n```\nEneighbor 是通过 do {}while 遍历求出来的临近所有顶点的坐标和\n这是我处理旧顶点的方式，和文档计算一摸一样，我觉得不是这个的问题导致细分错误。\n另外在进行一次 loop 后，我会把所有的顶点和边全部置为旧的，这个操作也没问题。","cooked":"

我在几何部分第二章 loop 面细分的时候，在 cube 细分达不到文档的要求，面是平滑了，但是八个顶点突出。
\n $\"QQ$
\n这是一次细分，可以看到八个顶点明显
\n $\"QQ$
\n这是四次细分，八个顶点更为明显
\n我觉得应该是新顶点计算有问题，我的实现如下：

if (v->is_new)\n        {\n            Eigen::Vector3f V1AddV2 = Eigen::Vector3f::Zero();\n            Eigen::Vector3f V3AddV4 = Eigen::Vector3f::Zero();\n            Eigen::Vector3f V1V2V3V4 = Eigen::Vector3f::Zero();\n            Halfedge* h = v->halfedge;\n            do\n            {\n                if(h->edge->is_new)\n                {\n                    V1V2V3V4 += h->inv->from->pos;\n                }\n                h = h->inv->next;\n            } while (h != v->halfedge);\n            V1AddV2 = 2.0f * v->pos;\n            V3AddV4 = V1V2V3V4 - V1AddV2;//不知道为啥这里要减去 V3AddV4 才第一次 cube 正常，之后就不正常，别的物体不正常。\n            v->new_pos = 3.0f/8.0f * V1AddV2 + 1.0f/8.0f * V3AddV4;\n        }\n

$\"QQ$

方法是在 split 生成的新点坐标正好是 1/2（V1+V2），就可以作为一个整体来计算 V1+V2 和 V3+V4，通过对一个新顶点周围顶点对 split 产生的新边遍历求坐标和（图中红色顶点的另外四个虚线的点的坐标和），正好是 V1+V2+V3+V4，V1+V2 是它中点新顶点坐标的两倍，然后这俩求差就可以得到 V3+V4;我觉得这么计算没有问题。但是遇到了怪异的现象：
\n1.V3AddV4 = V1V2V3V4 - V1AddV2;我把后面减去的 V1+V2 换成 V3+V4（这个显然是错误的），cube 第一次细分和文档一摸一样，但是后续细分错误，其他的 cow 和兔子错误，尖锐和凹陷明显。
\n2.V3AddV4 = V1V2V3V4 - V1AddV2;仍然用这个，cube 细分错误，但是兔子和 cow 细分正确，他们的表面仔细检查十分平滑，没有突起和凹陷。
\n $\"QQ$
\n想知道是什么原因导致的部分正确，按理说要正确就全部正确吗。

if (n == 3)\n                {\n                    v->new_pos = (1.0f - (9.0f / 16.0f)) * v->pos + (3.0f / 16.0f) * Eneighbors;\n                }\n                else\n                {\n                    // For other vertices, apply the general rule based on degree\n                    v->new_pos = (1.0f - (3.0f / 8.0f)) * v->pos + (0.375f / n) * Eneighbors;\n                }\n

Eneighbor 是通过 do {}while 遍历求出来的临近所有顶点的坐标和
\n这是我处理旧顶点的方式，和文档计算一摸一样，我觉得不是这个的问题导致细分错误。
\n另外在进行一次 loop 后，我会把所有的顶点和边全部置为旧的，这个操作也没问题。

\n","createdAt":"2024-11-11 11:55:55.104771","updatedAt":"2024-11-11 11:55:55.104771","replyToPostNumber":null,"replyCount":0,"likeCount":0,"wordCount":630,"liked":false,"bookmarked":false,"polls":[],"reactions":[],"wiki":false,"version":1,"publicVersion":1,"deletedAt":null,"isAnonymous":false,"user":{"id":6145,"username":"SunXueHao","name":"rss","avatarUrl":null,"admin":false,"moderator":false,"signatureRaw":null,"signatureCooked":null}},{"id":43184,"postNumber":2,"raw":"曲面细分对于 cube 有个整体往回缩的趋势，八个顶点突出考虑是因为旧顶点坐标计算错了往回缩的不够多。可能要看下更新顺序问题？比如要在对每条边分裂之前先计算旧顶点新坐标。\n对于新顶点的计算，如果你是将所有边分裂完之后再用你现在的算法，你应该注意到由于旁边的边已经分裂，在半边网络数据里产生了新顶点，公式里的 v3 和 v4 两个顶点与新顶点已经**不邻接**，你的 V3addV4 可能变成了相邻的两个新顶点相加，你应该在分裂前就计算。","cooked":"

曲面细分对于 cube 有个整体往回缩的趋势，八个顶点突出考虑是因为旧顶点坐标计算错了往回缩的不够多。可能要看下更新顺序问题？比如要在对每条边分裂之前先计算旧顶点新坐标。
\n对于新顶点的计算，如果你是将所有边分裂完之后再用你现在的算法，你应该注意到由于旁边的边已经分裂，在半边网络数据里产生了新顶点，公式里的 v3 和 v4 两个顶点与新顶点已经不邻接，你的 V3addV4 可能变成了相邻的两个新顶点相加，你应该在分裂前就计算。

\n","createdAt":"2024-11-12 00:24:35.799802","updatedAt":"2024-11-12 01:04:28.711828","replyToPostNumber":null,"replyCount":0,"likeCount":0,"wordCount":195,"liked":false,"bookmarked":false,"polls":[],"reactions":[],"wiki":false,"version":6,"publicVersion":6,"deletedAt":null,"isAnonymous":false,"user":{"id":3740,"username":"Lskyi","name":null,"avatarUrl":null,"admin":false,"moderator":false,"signatureRaw":null,"signatureCooked":null}}],"hasMore":false,"lastReadPostNumber":1,"hasUrlPostNumber":false}