<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'>
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 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> </body>
</html>