Subversion Repositories gelsvn

Rev

Go to most recent revision | Details | Last modification | View Log | RSS feed

Rev Author Line No. Line
81 jab 1
#include "CGLA/Mat4x4f.h"
2
#include "Util/ArgExtracter.h"
3
#include "Geometry/RGrid.h"
4
#include "Geometry/load_raw.h"
5
#include "Geometry/save_raw.h"
178 bj 6
#include "HMesh/obj_save.h"
7
#include "HMesh/x3d_save.h"
8
#include "HMesh/caps_and_needles.h"
185 jab 9
#include "HMesh/volume_polygonize.h"
10
#include "HMesh/triangulate.h"
81 jab 11
 
12
using namespace CGLA;
13
using namespace Util;
14
using namespace HMesh;
15
using namespace std;
16
using namespace Geometry;
17
 
18
typedef RGrid<unsigned char> RGridb;
19
 
20
int main(int argc, char**argv)
21
{
192 jab 22
 
81 jab 23
		ArgExtracter ae(argc, argv);
97 bj 24
		float iso = 40;
81 jab 25
		ae.extract("-i", iso);
97 bj 26
		int X  = 128;
81 jab 27
		ae.extract("-x", X);
97 bj 28
		int Y  = 128;
81 jab 29
		ae.extract("-y", Y);
97 bj 30
		int Z  = 62;
81 jab 31
		ae.extract("-z", Z);
192 jab 32
 
33
		string ofile = "../../data/isosurf.x3d";
34
		ae.extract("-o", ofile);
35
 
105 jab 36
		string file;
37
		if(ae.no_remaining_args()<2) 
38
				file ="../../data/teddybear.raw";
39
		else
40
				file = ae.get_last_arg();
81 jab 41
 
42
		RGridb grid(Vec3i(X,Y,Z));
178 bj 43
 
44
//		RGridb grid(Vec3i(32));
45
// 		srand(0);
46
// 		for(int i=0;i<32;++i)
47
// 			for(int j=0;j<32;++j)
48
// 				for(int k=0;k<32;++k)		
49
// 					grid[Vec3i(i,j,k)] = 80.0 * rand()/double(RAND_MAX);
50
 
51
 		cout << "loading " << file << endl;
52
 		load_raw(file,grid);
81 jab 53
 
54
		cout << "Polygonizing, iso = " << iso << endl;
55
		Manifold mani;
178 bj 56
		cuberille_polygonize(grid, mani, iso, true);
185 jab 57
		shortest_edge_triangulate(mani);
58
 
81 jab 59
		cout << "Saving" << endl;
192 jab 60
		x3d_save(ofile, mani);
81 jab 61
}