ccolin/pgai_tp1
Archived
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

Compare commits

base: ccolin:6310c2c551f9eb41397b588e2891a5bae477c160
Branches Tags
ccolin:master
..
compare: ccolin:1772005c1b72b9ea10b0c73954269ef60ee4d8aa
Branches Tags
ccolin:master

No commits in common. "6310c2c551f9eb41397b588e2891a5bae477c160" and "1772005c1b72b9ea10b0c73954269ef60ee4d8aa" have entirely different histories.
6310c2c551 ... 1772005c1b

1 changed files with 5 additions and 49 deletions
Show Stats Expand all files Collapse all files

54
src/main.cpp
View File

@ -1,11 +1,8 @@
#include "my_mesh.h" #include "my_mesh.h"
#include <ctgmath>
#include <iostream> #include <iostream>
#include <fstream> #include <fstream>
#define PI 3.14159265
bool check_faces_are_triangles(MyMesh &mesh) { bool check_faces_are_triangles(MyMesh &mesh) {
for (auto f_it = mesh.faces_begin(); f_it != mesh.faces_end(); ++f_it) { for (auto f_it = mesh.faces_begin(); f_it != mesh.faces_end(); ++f_it) {
@ -92,39 +89,11 @@ void stats_surface_area(MyMesh &mesh) {
} }
void stats_n_neighbors(MyMesh &mesh) {
for (const VertexHandle &vh : mesh.vertices()) {
unsigned count = 0;
for (auto vv_it = mesh.vv_iter(vh); vv_it.is_valid(); ++vv_it) {
count++;
}
std::cout << count << " ";
}
std::cout << std::endl;
}
void stats_normal_deviation(MyMesh &mesh) {
mesh.update_normals();
for (const VertexHandle &vh : mesh.vertices()) {
MyMesh::Normal normal = mesh.normal(vh);
float max = 0;
for (auto vf_it = mesh.vf_iter(vh); vf_it.is_valid(); ++vf_it) {
float angle = acos(OpenMesh::dot(mesh.normal(*vf_it), normal)) * 180.0 / PI;
if (angle > max)
max = angle;
}
std::cout << max << " ";
}
std::cout << std::endl;
}
int main(int argc, char *argv[]) { int main(int argc, char *argv[]) {
using namespace std; using namespace std;
if (argc != 2 && argc != 3) { if (argc != 2) {
cerr << "Utilisation : " << argv[0] << " fichier.obj [surface|valence|deviation]" << endl; cerr << "Utilisation : " << argv[0] << " [fichier.obj]" << endl;
return 1; return 1;
} }
MyMesh mesh; MyMesh mesh;
@ -141,21 +110,8 @@ int main(int argc, char *argv[]) {
cerr << "Les arêtes sont sur une face : " cerr << "Les arêtes sont sur une face : "
<< (check_edges_arent_lonely(argv[1]) ? "oui" : "non") << (check_edges_arent_lonely(argv[1]) ? "oui" : "non")
<< endl; << endl;
float area = total_area(mesh); cerr << "Aire totale : " << total_area(mesh) << endl;
cerr << "Aire totale : " << area << endl; // stats_to_csv(stats_surface_area(mesh));
cerr << "Aire moyenne : " << area / mesh.n_faces() << endl; stats_surface_area(mesh);
if (argc == 3) {
string cmd = argv[2];
if (cmd == "surface")
stats_surface_area(mesh);
else if (cmd == "valence")
stats_n_neighbors(mesh);
else if (cmd == "deviation")
stats_normal_deviation(mesh);
else {
cerr << "Commande inconnue : " << cmd << endl;
return 1;
}
}
return 0; return 0;
} }