Subversion Repositories gelsvn

Rev

Rev 90 | Only display areas with differences | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 90 Rev 94
1
#include "CGLA/Mat4x4f.h"
1
#include "CGLA/Mat4x4f.h"
2
#include "Util/ArgExtracter.h"
2
#include "Util/ArgExtracter.h"
3
#include "Geometry/RGrid.h"
3
#include "Geometry/RGrid.h"
4
#include "Geometry/load_raw.h"
4
#include "Geometry/load_raw.h"
5
#include "Geometry/save_raw.h"
5
#include "Geometry/save_raw.h"
6
#include "HMeshUtil/obj_save.h"
6
#include "HMeshUtil/obj_save.h"
7
#include "HMeshUtil/x3d_save.h"
7
#include "HMeshUtil/x3d_save.h"
8
#include "HMeshUtil/caps_and_needles.h"
8
#include "HMeshUtil/caps_and_needles.h"
9
#include "HMeshUtil/fair_polygonize.h"
9
#include "HMeshUtil/fair_polygonize.h"
10
 
10
 
11
using namespace CGLA;
11
using namespace CGLA;
12
using namespace HMeshUtil;
12
using namespace HMeshUtil;
13
using namespace Util;
13
using namespace Util;
14
using namespace HMesh;
14
using namespace HMesh;
15
using namespace std;
15
using namespace std;
16
using namespace Geometry;
16
using namespace Geometry;
17
 
17
 
18
typedef RGrid<unsigned char> RGridb;
18
typedef RGrid<unsigned char> RGridb;
19
 
19
 
20
int main(int argc, char**argv)
20
int main(int argc, char**argv)
21
{
21
{
22
		ArgExtracter ae(argc, argv);
22
		ArgExtracter ae(argc, argv);
23
		float iso = 30;
23
		float iso = 30;
24
		ae.extract("-i", iso);
24
		ae.extract("-i", iso);
25
		int X  = 256;
25
		int X  = 256;
26
		ae.extract("-x", X);
26
		ae.extract("-x", X);
27
		int Y  = 256;
27
		int Y  = 256;
28
		ae.extract("-y", Y);
28
		ae.extract("-y", Y);
29
		int Z  = 256;
29
		int Z  = 256;
30
		ae.extract("-z", Z);
30
		ae.extract("-z", Z);
31
		string file = ae.get_last_arg();
31
		string file = ae.get_last_arg();
32
		
32
		
33
		RGridb grid(Vec3i(X,Y,Z));
33
		RGridb grid(Vec3i(X,Y,Z));
34
		load_raw(file,grid);
34
		load_raw(file,grid);
35
		
35
		
36
		
36
		
37
		cout << "Polygonizing, iso = " << iso << endl;
37
		cout << "Polygonizing, iso = " << iso << endl;
38
		Manifold mani;
38
		Manifold mani;
39
		fair_polygonize(grid, mani, iso);
39
		fair_polygonize(grid, mani, iso);
40
 
40
 
41
		cout << "Removing caps and needles ... " << endl;
41
		cout << "Removing caps and needles ... " << endl;
42
		remove_caps_from_trimesh(mani, M_PI * 0.85);
42
		remove_caps_from_trimesh(mani, M_PI * 0.85);
43
		remove_needles_from_trimesh(mani, 1e-2);
43
		remove_needles_from_trimesh(mani, 1e-2);
44
		
44
		
45
		cout << "Saving" << endl;
45
		cout << "Saving" << endl;
46
		obj_save("isosurf.obj", mani);
46
		obj_save("isosurf.obj", mani);
47
}
47
}
48
 
48