<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
mso-margin-top-alt:auto;
margin-right:0in;
mso-margin-bottom-alt:auto;
margin-left:0in;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle18
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri",sans-serif;}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Works well for me in 3D Slicer’s embedded Python 2.7 64 bit, both with latest VTK master and VTK 7.1.0.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Andras<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif"> vtkusers <vtkusers-bounces@public.kitware.com>
<b>On Behalf Of </b>Andrea Gavana<br>
<b>Sent:</b> Wednesday, September 26, 2018 10:23 AM<br>
<b>To:</b> vtkusers@public.kitware.com<br>
<b>Subject:</b> [vtkusers] Simple 2D unstructured grid<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div>
<div>
<p class="MsoNormal">Hello list,<o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"> I m trying to create a very, very simple 2D unstructured grid but every time I run my script - no matter what modifications I do - Python crashes hard with a Windows message "Python has stopped working...".<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">A simplified version of my grid is in the picture below (done in Excel, sorry...):<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal"><img width="481" height="289" style="width:5.0104in;height:3.0104in" id="_x0000_i1025" src="cid:image001.png@01D4558D.56BD6730" alt="image.png"><o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">I attach a copy of my script below. Could anyone offer some suggestions on what I am doing wrong? I am using VTK 7.1.0 on Windows 7 64 bit, Python 2.7 64 bit - VTK comes from Christophe Gohlke very useful Python binary wheels.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Thank you in advance for your help.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Andrea.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># -- BEGIN CODE -- #<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<div>
<p class="MsoNormal">import numpy<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">import vtk<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">import vtk.util.numpy_support as numpy_support<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># X, Y coordinates of a simple grid - Z is zero<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">X = [600.4957421, 600.5, 600.5, 600.4957421, 600.4913048, 600.5,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 600.5, 600.4957421, 600.4913048, 600.4868674, 600.5, 600.5,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 600.4913048, 600.4868674, 600.4824301, 600.5, 600.5, 600.4868674,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 600.4824301, 600.4779928, 600.5, 600.5, 600.4824301, 600.4779928,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 600.4735554, 600.5, 600.5, 600.4779928, 600.4735554]<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Y = [2940.5, 2940.5, 2940.404044, 2940.5, 2940.6, 2940.6, 2940.5,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 2940.5, 2940.6, 2940.7, 2940.7, 2940.6, 2940.6, 2940.7, 2940.8,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 2940.8, 2940.7, 2940.7, 2940.8, 2940.9, 2940.9, 2940.8, 2940.8,<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> 2940.9, 2941, 2941, 2940.9, 2940.9, 2941]<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># First face is triangular, everything else is a 4-sided polygon<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">IDS = [range(4)]<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">IDS += numpy.split(numpy.arange(4, len(X)), 5)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">npoints = len(X)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">matrix = numpy.zeros((npoints, 3), numpy.float32)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">matrix[:, 0] = X<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">matrix[:, 1] = Y<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># Create the grid points<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">vtk_pts = vtk.vtkPoints()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">vtk_pts.SetData(numpy_support.numpy_to_vtk(matrix, deep=1))<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># Create the unstructured grid<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">grid = vtk.vtkUnstructuredGrid()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">grid.SetPoints(vtk_pts)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># Allocate space for the cells in the grid<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">nc = len(IDS)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">grid.Allocate(nc)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># Loop through all cells<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">for i in xrange(nc):<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> cell_ids = IDS[i]<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> ncoords = len(cell_ids)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> grid.InsertNextCell(vtk.VTK_POLYGON, ncoords, cell_ids)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">print grid.GetNumberOfCells(), nc<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">mapper = vtk.vtkDataSetMapper()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">mapper.SetInputData(grid)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">actor = vtk.vtkActor()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">actor.SetMapper(mapper)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">ren = vtk.vtkRenderer()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">renWin = vtk.vtkRenderWindow()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">renWin.AddRenderer(ren)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">iren = vtk.vtkRenderWindowInteractor()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">iren.SetRenderWindow(renWin)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># Add the actors to the renderer, set the background and size<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">ren.AddActor(actor)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">ren.SetBackground(0, 0, 0)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># This allows the interactor to initalize itself. It has to be<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"># called before an event loop.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">iren.Initialize()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">ren.ResetCamera()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">renWin.Render()<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"># Start the event loop.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">iren.Start()<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</div>
</body>
</html>