[Paraview] Getting plate boundary data into Paraview

[Lester Anderson]

Cheers Dave that worked perfectly :)

Lester

> From: dave.demarle at kitware.com
> Date: Tue, 31 Aug 2010 09:13:18 -0400
> Subject: Re: [Paraview] Getting plate boundary data into Paraview
> To: lester_anderson1963 at hotmail.com
> CC: paraview at paraview.org
> 
> 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
> >
 		 	   		  
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.paraview.org/pipermail/paraview/attachments/20100831/cff2e31f/attachment.htm>