[Paraview] Getting plate boundary data into Paraview

[David E DeMarle]

Change

pdo.SetStrips(cells)
to
pdo.SetLines(cells)

David E DeMarle
Kitware, Inc.
R&D Engineer
28 Corporate Drive
Clifton Park, NY 12065-8662
Phone: 518-371-3971 x109



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