igllib 202_gaussian curvature

最新推荐文章于 2025-06-08 09:21:36 发布

seamanj

最新推荐文章于 2025-06-08 09:21:36 发布

阅读量1.3k

点赞数

CC 4.0 BY-SA版权

分类专栏： IGL library 数学积累

本文链接：https://blog.youkuaiyun.com/seamanj/article/details/51550084

数学积累同时被 2 个专栏收录

51 篇文章

订阅专栏

IGL library

21 篇文章

订阅专栏

本文指出了离散高斯曲率公式中存在的错误，并通过解析论文与代码示例进行详细说明。针对《三角网格上的离散微分几何算子》一文中提到的方法进行了讨论，并提供了计算平均高斯曲率的正确实现。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

the original tutor is like:

there are some points I wanna mention:

1. the discrete Gaussian curvature formula it listed is wrong. How? listen to me explain.

According to the paper "Discrete Differential-Geometry Operators for Triangulated 2-Manifolds" section 4.2, it says

In addition, the code also told us about that:

#include <igl/gaussian_curvature.h>
#include <igl/massmatrix.h>
#include <igl/invert_diag.h>
#include <igl/readOFF.h>
#include <igl/viewer/Viewer.h>
#include <igl/jet.h>
#include "tutorial_shared_path.h"

int main(int argc, char *argv[])
{
  using namespace Eigen;
  using namespace std;
  MatrixXd V;
  MatrixXi F;
  igl::readOFF(TUTORIAL_SHARED_PATH "/bumpy.off",V,F);

  VectorXd K;
  // Compute integral of Gaussian curvature
  igl::gaussian_curvature(V,F,K);
  // Compute mass matrix
  SparseMatrix<double> M,Minv;
  igl::massmatrix(V,F,igl::MASSMATRIX_TYPE_DEFAULT,M);
  igl::invert_diag(M,Minv);
  // Divide by area to get integral average, note this line
  K = (Minv*K).eval();

  // Compute pseudocolor
  MatrixXd C;
  igl::jet(K,true,C);

  // Plot the mesh with pseudocolors
  igl::viewer::Viewer viewer;
  viewer.data.set_mesh(V, F);
  viewer.data.set_colors(C);
  viewer.launch();
}

2. I want to offer some material to help readers understand the geometrical meaning