K-means的缺点

昨天记录了使用K-means对网格模型进行分割的步骤和一些简单的结果，从昨天的实验结果来看，使用顶点坐标和顶点法向作为K-means聚类的特征得到的分割效果总体上还是不错的。分割结束后，每个顶点会被赋予一个分割的标号。
但是，只使用顶点距离作为聚类得到的结果并没有语义，因此可能会得到一些不太好的聚类结果。

如上图所示，“马”模型的右后腿和两条左腿被分割到了同一个类，这使得分割的结果不具有连通性，因此在后续的学习中，需要对该算法的分割结果进行调整，或者需找更加合适的分割方法。

分割边界的获取

边界获取的思想挺简单，首先遍历网格顶点，筛选出某一特定区域的顶点，对于这一区域中某一顶点，判断其一邻域中的顶点的标号与其自身标号是否相同，如果存在不同标号与自身标号不同的顶点，那么这样的顶点就是边界顶点。找出边界顶点后，为其设置边界标记。同时，记录区域中非边界顶点的个数，以便于计算区域局部的拉普拉斯矩阵。

寻找边界代码如下

void find_Boundary(PolygonMesh::Mesh * _mesh, int area)
{/*_mesh: 输入的网格area: 区域编号*/num_vtx = _mesh->n_vertices();is_boundray.clear();for ( int i = 0; i < num_vtx; i++ ){//默认所有顶点都不是边界当某个顶点的1领域中有一个顶点标号和它不同时，这个顶点为边界顶点is_boundray.push_back(0);}int v_idx = 0;cnt_non_bdy = 0;//非边界顶点个数for (auto v_it = _mesh->vertices_begin();v_it != _mesh->vertices_end();++v_it,v_idx++){if (all_lables[v_idx] == area){for ( auto vv_it = _mesh->vv_begin(v_it); vv_it != _mesh->vv_end(v_it); ++vv_it){OpenMesh::VertexHandle vertex2  = vv_it.handle();int v2_idx = vertex2.idx();//当某个顶点的1领域中有一个顶点标号和它不同时，这个顶点为边界顶点if ( all_lables[v2_idx] !=area ){is_boundray[v_idx]=1;break;}}   }//区域area中的非边界顶点if (all_lables[v_idx] == area && is_boundray[v_idx] !=area ){cnt_non_bdy++;}}
}

需要注意的是，在计算区域的拉普拉斯矩阵时，如果使用L = D - A这种形式的拉普拉斯矩阵，在计算顶点度矩阵和顶点邻接矩阵A时，都要将其邻域中的边界迪昂点排除在外

计算度矩阵和邻接矩阵代码

    int v_idx = 0;int n_vidx = 0;//对同一区域中的顶点重新编号VectorX byd_neibor;//记录每个顶点的边界邻域byd_neibor.setZero(cnt_non_bdy);for(auto v_it = _mesh->vertices_begin();v_it != _mesh->vertices_end();++v_it,++v_idx){if (all_lables[v_idx] == 1 && is_boundray[v_idx] !=1 ){//计算其边界点邻居个数for ( auto vv_it = _mesh->vv_begin(v_it); vv_it != _mesh->vv_end(v_it); ++vv_it){OpenMesh::VertexHandle vertex2  = vv_it.handle();int v2_idx = vertex2.idx();//当某个顶点的1领域中有一个顶点标号和它不同时，这个顶点为边界顶点if ( all_lables[v2_idx] !=1 ){byd_neibor[n_vidx]++;}}   //获得顶点i的度,然后减去边界点个数int val = _mesh->valence(v_it.handle());D_matrix(n_vidx ,n_vidx ) = val - byd_neibor[n_vidx];//获得顶点的1邻域邻接矩阵OpenMesh::VertexHandle vertex1 = v_it.handle();int v1_idx = vertex1.idx();for ( auto vv_it = _mesh->vv_begin(v_it); vv_it != _mesh->vv_end(v_it); ++vv_it){OpenMesh::VertexHandle vertex2  = vv_it.handle();int v2_idx = vertex2.idx();if (is_boundray[v2_idx] != 1)//筛选非边界的邻居顶点{A_matrix(n_vidx,n_vidx) = 1;}}n_vidx ++;//对同一区域中的顶点进行新的编号}}