这里写目录标题

openmesh的文件操作读文件读文件修改文件1-read_data（）读文件修改文件2- mesh_scaling() 写文件

openmesh的文件操作

//读文件 OpenMesh::IO::read_mesh(mesh, filename); //写文件 OpenMesh::IO::write_mesh(mesh, filename);

读文件

Mesh src_mesh; Mesh tar_mesh; std::string src_file; std::string tar_file; std::string out_file; RegParas paras; //读文件 Registration* reg;//我新建了一个类，要读取参数并进行初始化 bool use_demean = false; read_data(src_file, src_mesh, use_demean); read_data(tar_file, tar_mesh, use_demean); if(src_mesh.n_vertices()==0 || tar_mesh.n_vertices()==0) { exit(0); } double scale; scale = mesh_scaling(src_mesh, tar_mesh); scale = 1; //写文件 write_data(out_file.c_str(), src_mesh, use_demean, scale); std::cout << "write result to " << out_file << std::endl; delete reg;

read_data(src_file, src_mesh, use_demean);

bool read_data(const std::string filename, Mesh& mesh, bool use_demean) { //将文件读到mesh中 bool read_OK = OpenMesh::IO::read_mesh(mesh, filename); std::cout << "filename = " << filename << std::endl; if (read_OK) { mesh.request_vertex_status(); mesh.request_edge_status(); mesh.request_face_status(); mesh.request_face_normals(); mesh.request_vertex_normals(); printBasicMeshInfo(mesh, use_demean); mesh.update_face_normals(); mesh.update_vertex_normals(); return true; } std::cout << "#vertices = " << mesh.n_vertices() << std::endl; return false; }

读文件修改文件1-read_data（）

printBasicMeshInfo()中更新点的坐标： $$MeshScales=[\frac{1}{n}\sum _{i=0}^N{x}_i,\frac{1}{n}\sum _{i=0}^N{y}_i,\frac{1}{n}\sum _{i=0}^N{z}_i]$$ $$V_i^{\prime }=V_i-MeshScales$$ 即： $$V_i^{\prime }=[x_i,y_i,z_i]-[\frac{1}{n}\sum _{i=0}^N{x}_i,\frac{1}{n}\sum _{i=0}^N{y}_i,\frac{1}{n}\sum _{i=0}^N{z}_i]$$

void printBasicMeshInfo(Mesh& mesh, bool use_demean) { if (mesh.n_vertices() == 0) printf("No Mesh\n"); switch(checkMeshMode(mesh)) { case TRIANGLE: printf("Triangle Mesh.\n"); break; case QUAD: printf("Quadrilateral Mesh.\n"); break; default: printf("General Mesh.\n"); break; } printf("Information of the input mesh:\nVertex : %d;\nFace : %d;\nEdge : %d, HalfEdge : %d\n\n", mesh.n_vertices(),mesh.n_faces(),mesh.n_edges(),mesh.n_halfedges()); OpenMesh::Vec3d MeshScales; MeshScales[0] = 0.0; MeshScales[1] = 0.0; MeshScales[2] = 0.0; for (Mesh::VertexIter v_it = mesh.vertices_begin(); v_it != mesh.vertices_end(); ++v_it) { MeshScales += mesh.point(*v_it); } MeshScales /= mesh.n_vertices(); mesh.data(mesh.vertex_handle(0)).mesh_trans = MeshScales; if(use_demean) { for (Mesh::VertexIter v_it = mesh.vertices_begin(); v_it != mesh.vertices_end(); ++v_it) { mesh.point(*v_it) = mesh.point(*v_it)- MeshScales; } } }

读文件修改文件2- mesh_scaling()

mesh_scaling()更新点的坐标： 找到src中最大的点$v_{max}$和最小的点$v_{min}$ $$scale=|v_{max}-v_{min}|$$ $$[x_i^{\prime },y_i^{\prime },z_i^{\prime }]=\frac{1}{scale}[x_i,y_i,z_i]$$

double mesh_scaling(Mesh& src_mesh, Mesh& tar_mesh) { OpenMesh::Vec3d max(-1e12, -1e12, -1e12); OpenMesh::Vec3d min(1e12, 1e12, 1e12); for(auto it = src_mesh.vertices_begin(); it != src_mesh.vertices_end(); it++) { for(int j = 0; j < 3; j++) { if(src_mesh.point(*it)[j] > max[j]) { max[j] = src_mesh.point(*it)[j]; } if(src_mesh.point(*it)[j] < min[j]) { min[j] = src_mesh.point(*it)[j]; } } } for(auto it = tar_mesh.vertices_begin(); it != tar_mesh.vertices_end(); it++) { for(int j = 0; j < 3; j++) { if(tar_mesh.point(*it)[j] > max[j]) { max[j] = tar_mesh.point(*it)[j]; } if(tar_mesh.point(*it)[j] < min[j]) { min[j] = tar_mesh.point(*it)[j]; } } } double scale = (max-min).norm(); for(auto it = src_mesh.vertices_begin(); it != src_mesh.vertices_end(); it++) { OpenMesh::Vec3d p = src_mesh.point(*it); p = p/scale; src_mesh.set_point(*it, p); } for(auto it = tar_mesh.vertices_begin(); it != tar_mesh.vertices_end(); it++) { OpenMesh::Vec3d p = tar_mesh.point(*it); p = p/scale; tar_mesh.set_point(*it, p); } return scale; }

写文件

将src网格写到文件中去

if (use_demean) src_mesh.data(src_mesh.vertex_handle(0)).mesh_trans = tar_mesh.data(tar_mesh.vertex_handle(0)).mesh_trans; write_data(out_file.c_str(), src_mesh, use_demean, scale); std::cout << "write result to " << out_file << std::endl;

bool write_data(const char* filename, Mesh& mesh, bool use_demean, double scale)

bool write_data(const char* filename, Mesh& mesh, bool use_demean, double scale) { OpenMesh::Vec3d MeshScales; MeshScales = mesh.data(mesh.vertex_handle(0)).mesh_trans; for (Mesh::VertexIter v_it = mesh.vertices_begin(); v_it != mesh.vertices_end(); ++v_it) { mesh.point(*v_it) = mesh.point(*v_it)*scale; } if(use_demean) { for (Mesh::VertexIter v_it = mesh.vertices_begin(); v_it != mesh.vertices_end(); ++v_it) { mesh.point(*v_it) = mesh.point(*v_it) + MeshScales; } } bool ok = OpenMesh::IO::write_mesh(mesh, filename); return ok; }