File: poissonrecon.py

package info (click to toggle)
vistrails 2.1.1-1
  • links: PTS
  • area: main
  • in suites: jessie, jessie-kfreebsd
  • size: 74,208 kB
  • ctags: 46,250
  • sloc: python: 316,267; xml: 52,512; sql: 3,627; php: 731; sh: 260; makefile: 108
file content (101 lines) | stat: -rw-r--r-- 3,583 bytes parent folder | download
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
101
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.checkInputPort('pts')
        f1 = self.getInputFromPort('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.setResult('error_mesh_1', error_mesh_output[0])
            self.setResult('error_mesh_2', error_mesh_output[1])

        if need_histograms:
            self.setResult('error_hist_1', error_hist_output[0])
            self.setResult('error_hist_2', error_hist_output[1])

        self.need_error_meshes = need_error_meshes
        self.need_histograms = need_histograms

        self.setResult('hausdorff_error', float(results[0]))
        self.setResult('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)