-
Notifications
You must be signed in to change notification settings - Fork 1
/
poissonrecon.py
101 lines (77 loc) · 3.5 KB
/
poissonrecon.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
import core.modules
import core.modules.module_registry
from core.modules.vistrails_module import Module, ModuleError
import core.bundles
import core.requirements
import os
import re
from core.configuration import ConfigurationObject
configuration = ConfigurationObject(executable_path=(None, str))
identifier = 'edu.utah.sci.cscheid.poissonrecon'
version = '0.1'
name = 'Poisson Surface Reconstruction'
##############################################################################
class PoissonRecon(Module):
def __init__(self):
Module.__init__(self)
def compute(self):
self.check_input('pts')
f1 = self.get_input('mesh1')
make_file = self.interpreter.filePool.create_file
metro_output = make_file(suffix='.ply')
values = ['', f1.name, f2.name, '-q']
need_error_meshes = False
if (self.outputPorts.has_key('error_mesh_1') or
self.outputPorts.has_key('error_mesh_2')):
need_error_meshes = True
error_mesh_output = (make_file(suffix='.ply'),
make_file(suffix='.ply'))
values.append('-c')
values.append(error_mesh_output[0].name)
values.append(error_mesh_output[1].name)
need_histograms = False
if (self.outputPorts.has_key('error_hist_1') or
self.outputPorts.has_key('error_hist_2')):
need_histograms = True
error_hist_output = (make_file(suffix='.csv'),
make_file(suffix='.csv'))
values.append('-h')
values.append(error_hist_output[0].name)
values.append(error_hist_output[1].name)
print values
print ('%s ' * len(values))
cmdline = ('%s ' * len(values)) % tuple(values)
cmdline += '> '
cmdline += metro_output.name
print cmdline
result = os.system(cmdline)
if result != 0:
raise ModuleError(self, "Execution failed")
l = [x for x in file(metro_output.name).readlines()
if x.startswith('Hausdorff')][0]
regexp = re.compile(r'Hausdorff distance: ([0123456789.]+) \(([0123456789.]+) wrt bounding box diagonal\)')
results = regexp.match(l).groups()
if need_error_meshes:
self.set_output('error_mesh_1', error_mesh_output[0])
self.set_output('error_mesh_2', error_mesh_output[1])
if need_histograms:
self.set_output('error_hist_1', error_hist_output[0])
self.set_output('error_hist_2', error_hist_output[1])
self.need_error_meshes = need_error_meshes
self.need_histograms = need_histograms
self.set_output('hausdorff_error', float(results[0]))
self.set_output('normalized_hausdorff_error', float(results[1]))
##############################################################################
def initialize():
if not core.requirements.executable_file_exists('PoissonRecon'):
raise core.requirements.MissingRequirement("PoissonRecon")
reg = core.modules.module_registry
reg.add_module(PoissonRecon)
reg.add_input_port(PoissonRecon, 'pts',
core.modules.basic_modules.File)
reg.add_input_port(PoissonRecon, 'binary',
core.modules.basic_modules.Boolean)
reg.add_input_port(PoissonRecon, 'depth',
core.modules.basic_modules.Integer)
reg.add_output_port(PoissonRecon, 'output',
core.modules.basic_modules.File)