[ create a new paste ] login | about

Link: http://codepad.org/VPOa13St    [ raw code | output | fork ]

Python, pasted on Jul 15: 
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
102
103
104
105

#!/usr/bin/python 

import vtk
import math

def make_rwi():
    """
    makes renderer, window, interactor
    """
    renderer = vtk.vtkRenderer()
    renderer.SetBackground(0.1, 0.2, 0.4)
    
    window = vtk.vtkRenderWindow()
    window.SetSize(600, 600)
    
    window.AddRenderer(renderer)
    interactor = vtk.vtkRenderWindowInteractor()
    interactor.SetRenderWindow( window )    

    return (renderer, window, interactor)    
    
def gen_triangles(cd):
    clsdPoly = vtk.vtkPolyData()
    clsdPoly.SetPoints(cd.GetPoints())
    clsdPoly.SetPolys(cd.GetLines())

    triangles = vtk.vtkTriangleFilter()
    triangles.SetInput(clsdPoly)
    
    return triangles.GetOutput()

def make_contour_data():
    points = vtk.vtkPoints()
    lines  = vtk.vtkCellArray()
    
    n = 20
    
    for i in range(0, n):
        angle = 2.0 * math.pi * float(i) / float(n)
        points.InsertPoint(i, 0.1 * math.cos(angle), 0.1 * math.sin(angle), 0.0 )
        lines.InsertNextCell(i)

    lines.InsertNextCell(0)

    pd = vtk.vtkPolyData()
 
    pd.SetPoints(points)
    pd.SetLines(lines)
    
    return pd
    
def render_things(renderer, window, interactor):
    renderer.ResetCamera()
    window.Render()

    interactor.Initialize()
    interactor.Start()    
    
def write_poly(fname, pd):
    writer = vtk.vtkXMLPolyDataWriter()
    writer.SetInput(pd)
    writer.SetFileName(fname)
    writer.SetDataModeToAscii()
    writer.SetCompressorTypeToNone()
    writer.Write()

def All():
    renderer, window, interactor = make_rwi()

    cd = make_contour_data()
    #write_poly("before.vtp", cd)

    contourRep = vtk.vtkOrientedGlyphContourRepresentation()
    contourRep.GetLinesProperty().SetColor(1, 0, 0) # set color to red
    contourRep.GetProperty().SetColor(0,1,0)
    contourRep.GetActiveProperty().SetColor(0,0,0)
 
    contourWidget = vtk.vtkContourWidget()
    contourWidget.SetInteractor(interactor)
    contourWidget.SetRepresentation(contourRep)
    contourWidget.On()

    contourWidget.Initialize(cd, 1)
    contourWidget.CloseLoop()
    cd = contourWidget.GetRepresentation().GetContourRepresentationAsPolyData()
    write_poly("after.vtp", cd)
    
    contourWidget.Render()

    td = gen_triangles(cd)
    
    mapper = vtk.vtkDataSetMapper()
    mapper.SetInput(td)

    actor = vtk.vtkActor()
    actor.SetMapper(mapper)
    actor.GetProperty().SetColor(0,1,0)
    actor.GetProperty().SetOpacity(0.2)  #actor.GetProperty().SetRepresentationToWireframe()

    renderer.AddActor(actor)

    render_things(renderer, window, interactor)
    
if __name__ == "__main__":
    All()


Output: 
1
2
3
4

Traceback (most recent call last):
  Line 3, in <module>
    import vtk
ImportError: No module named vtk


Create a new paste based on this one


Comments: