<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>
Cheers Dave that worked perfectly :)<br><br>Lester<br><br>> From: dave.demarle@kitware.com<br>> Date: Tue, 31 Aug 2010 09:13:18 -0400<br>> Subject: Re: [Paraview] Getting plate boundary data into Paraview<br>> To: lester_anderson1963@hotmail.com<br>> CC: paraview@paraview.org<br>> <br>> Change<br>> <br>> pdo.SetStrips(cells)<br>> to<br>> pdo.SetLines(cells)<br>> <br>> David E DeMarle<br>> Kitware, Inc.<br>> R&D Engineer<br>> 28 Corporate Drive<br>> Clifton Park, NY 12065-8662<br>> Phone: 518-371-3971 x109<br>> <br>> <br>> <br>> On Tue, Aug 31, 2010 at 6:19 AM, Lester Anderson<br>> <lester_anderson1963@hotmail.com> wrote:<br>> > Hi Dave<br>> ><br>> > Thanks for the info.<br>> ><br>> > I added a programmable source and it worked fine to import the data, but it<br>> > only gave me points and not joined to give line data?<br>> ><br>> > Cheers<br>> ><br>> > Lester<br>> ><br>> >> From: dave.demarle@kitware.com<br>> >> Date: Mon, 30 Aug 2010 18:55:45 -0400<br>> >> Subject: Re: [Paraview] Getting plate boundary data into Paraview<br>> >> To: lester_anderson1963@hotmail.com<br>> >> CC: paraview@paraview.org<br>> >><br>> >> Pop this into a python programmable source and it should import your data.<br>> >><br>> >> import os<br>> >> import string<br>> >><br>> >> def insertNextCell():<br>> >> nextpoly = vtk.vtkPolyLine()<br>> >> nextpoly.GetPointIds().DeepCopy(segments)<br>> >> cells.InsertNextCell(nextpoly)<br>> >><br>> >> filename = os.path.normcase("/location/of/your/data.txt")<br>> >> f = open(filename)<br>> >><br>> >> pdo = self.GetPolyDataOutput()<br>> >><br>> >> pts = vtk.vtkPoints()<br>> >><br>> >> #an array for the fault information<br>> >> segmentNames = vtk.vtkStringArray()<br>> >> segmentNames.SetName("Fault Information")<br>> >> segmentNames.SetNumberOfComponents(1)<br>> >><br>> >> cells = vtk.vtkCellArray()<br>> >> segments = vtk.vtkIdList()<br>> >><br>> >> #traverse file<br>> >> for line in f:<br>> >> if line.startswith(">"):<br>> >> #when you find a new fault<br>> >> #save the one you just finished<br>> >> if segments.GetNumberOfIds() != 0:<br>> >> insertNextCell()<br>> >> segments = vtk.vtkIdList()<br>> >> segmentNames.InsertNextValue(line)<br>> >> else:<br>> >> #otherwise just keep remembering coordinates<br>> >> x,y = [float(n) for n in line.strip().split()[:2]]<br>> >> id = pts.InsertNextPoint(x,y,0)<br>> >> segments.InsertNextId(id)<br>> >> f.close()<br>> >><br>> >> if segments.GetNumberOfIds()>0:<br>> >> #don't forget the last fault in progress<br>> >> insertNextCell()<br>> >><br>> >> pdo.SetPoints(pts)<br>> >> pdo.SetStrips(cells)<br>> >> pdo.GetCellData().AddArray(segmentNames)<br>> >><br>> >><br>> >> David E DeMarle<br>> >> Kitware, Inc.<br>> >> R&D Engineer<br>> >> 28 Corporate Drive<br>> >> Clifton Park, NY 12065-8662<br>> >> Phone: 518-371-3971 x109<br>> >><br>> >><br>> >><br>> >> On Wed, Aug 25, 2010 at 5:02 AM, Lester Anderson<br>> >> <lester_anderson1963@hotmail.com> wrote:<br>> >> >> 2014 BOCONO TRANSFORM FAULT, N. SOUTH AMERICA<br>> >> > -58.9415 10.5593<br>> >> > -59.1817 10.5598<br>> >> > -59.4918 10.5899<br>> >> > -59.8086 10.6593<br>> >> > -60.1954 10.7256<br>> >> > -60.4489 10.7065<br>> >> > -60.8159 10.6678<br>> >> > -61.0594 10.6486<br>> >> > -61.4631 10.6068<br>> >> > -61.8535 10.5682<br>> >> > -62.2972 10.5395<br>> >> > -62.8876 10.4915<br>> >> > -63.7494 10.4103<br>> >> >> 2014 BOCONO TRANSFORM FAULT, N. SOUTH AMERICA<br>> >> > -77.5807 0.4294<br>> >> > -77.8448 0.1430<br>> >> > -77.9752 -0.0068<br>> >> > -78.2026 -0.3332<br>> >> > -78.4367 -0.6463<br>> >> > -78.7342 -0.9592<br>> >> > -79.0453 -1.3955<br>> >> > -79.2994 -1.6951<br>> >> > -79.4901 -1.9814<br>> >> > -79.6407 -2.2444<br>> >> > -79.7878 -2.4208<br>> >> > -79.8400 -2.5000<br>> >> > -79.9182 -2.5805<br>> >> > -80.0386 -2.7369<br>> >> > -80.1589 -2.8732<br>> >> > -81.9900 -4.8900<br>> >> >> 2014 SE CAYMAN TROUGH TRANSFORM FAULT<br>> >> > -78.4356 18.4315<br>> >> > -78.3248 18.4569<br>> >> > -78.1413 18.4710<br>> >> > -78.0147 18.4664<br>> >> > -77.9386 18.4565<br>> ><br> </body>
</html>