Subversion Repositories gelsvn

Rev

Rev 192 | Go to most recent revision | Details | Compare with Previous | 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();
205 jab 41
 
42
		//RGrid<unsigned short> grid(Vec3i(X,Y,Z));
81 jab 43
		RGridb grid(Vec3i(X,Y,Z));
178 bj 44
 
45
//		RGridb grid(Vec3i(32));
46
// 		srand(0);
47
// 		for(int i=0;i<32;++i)
48
// 			for(int j=0;j<32;++j)
49
// 				for(int k=0;k<32;++k)		
50
// 					grid[Vec3i(i,j,k)] = 80.0 * rand()/double(RAND_MAX);
51
 
52
 		cout << "loading " << file << endl;
53
 		load_raw(file,grid);
81 jab 54
 
55
		cout << "Polygonizing, iso = " << iso << endl;
56
		Manifold mani;
205 jab 57
		//cuberille_polygonize(grid, mani, iso, false);
58
		mc_polygonize(grid, mani, iso);
59
 
60
				//shortest_edge_triangulate(mani);
185 jab 61
 
81 jab 62
		cout << "Saving" << endl;
192 jab 63
		x3d_save(ofile, mani);
81 jab 64
}