[Paraview] Plot over Line source

Sun Sep 13 09:42:27 EDT 2015

Hi Richard,

Please keep your replies on the list with a "reply-all". Someone else may
have a similar question as you (or better yet, an answer) and this may come
in handy in a future search of the list serve.

I'm not sure how to annotate the position on the graph as you seek to do as
there is no mechanism for this. You could perhaps hack this in by creating
another column in your table where all the values are zero except for the
time that you are showing. Setting this column to the maximum of your
plotted y values in the row at the time you are showing would result in an
almost vertical line showing up at the time you are displaying if your data
samples in x are dense enough.

HTH,
Cory

On Fri, Sep 11, 2015 at 10:58 AM, no name <keinepostnurmuell at gmail.com>
wrote:

> Hey Cory
>
> I have one more question I am trieing to solve since we started mailing
> ... I got everything working so far using a programmable filter and I read
> my files into a vtkTable now, but I want to read more than one information
> and hand it to my LineChartView. At the moment I am doing:
>
> A = np.rec.fromarrays(switching_data_clean, dtype= {'names': ['Time', 'x',
> 'y'], 'formats':['int','f8', 'f8']})
>
>
> for name in A.dtype.names:
>
>     array = A[name]
>
>     output.RowData.append(array, name)
>
>  But I want to add a point data aswell (basically I am plotting a time dependent external field, and I want to point out where I am at the moment, so I have to be able to read one more Point in every timestep). I tried to do a new array with only that one point inside in a second Programmable Filter as vtkDATA, which didn't really work. Do you have another hint for me?
>
> Thanks
>
> Richard
>
>
> On 10 September 2015 at 15:34, Cory Quammen <cory.quammen at kitware.com>
> wrote:
>
>> Richard,
>>
>> I'm glad you got it working.
>>
>> cory
>>
>> On Wed, Sep 9, 2015 at 1:45 PM, no name <keinepostnurmuell at gmail.com>
>> wrote:
>>
>>> I just managed to do this myself and I had to do it the same way again,
>>> redo everything by copy and paste, there seems to be a bug in here too, but
>>> I am not capable of reproducing it :( I just wanted to let you know.
>>>
>>> Thanks for your Help!
>>>
>>> Richard
>>>
>>> On 9 September 2015 at 18:20, no name <keinepostnurmuell at gmail.com>
>>> wrote:
>>>
>>>> this is straight what I was looking for! Thank you SO MUCH for pointing
>>>> out this bug to me. This drove me mad :(
>>>>
>>>> One last question: Line Chart != PlotOverLine, so its hard for me to
>>>> plot x and y axis, or like in my script mx over hx, as nothing happens when
>>>> I choose to use one of my coordinates as "X-Axis Parameter" and just let
>>>> the other one being drawn.
>>>>
>>>> Can't I make a real Plot over Line?
>>>>
>>>> Thanks again for your time!
>>>> Richard
>>>>
>>>> On 9 September 2015 at 15:12, Cory Quammen <cory.quammen at kitware.com>
>>>> wrote:
>>>>
>>>>> Hi Richard,
>>>>>
>>>>> If you don't mind, let's keep the discussion on the list so that
>>>>> others may benefit from the discussion.
>>>>>
>>>>> Outputting a vtkTable is the right thing to do. What I suspect is
>>>>> happening with the first error message is that if you execute a
>>>>> programmable source with the output set to one type of data set (e.g.,
>>>>> vtkPolyData), there is a bug that doesn't let you switch to a different
>>>>> type of output. As a workaround, you can create a new Programmable Source,
>>>>> copy and paste your code from the original Programmable Source into the new
>>>>> one, and set the output type to vtkTableData *before* you click Apply.
>>>>>
>>>>> To view line plots of your data, you'll need to create a Line Chart
>>>>> View by splitting the 3D Render View. To do this, look for the icon on the
>>>>> top right of the 3D Render View that has two vertically-oriented rectangles
>>>>> side by side (tool tip "Split Horizontal"). In the list of views that comes
>>>>> up, click on "Line Chart View". If your line plots are not visible right
>>>>> away, you will need to make them visible by clicking on the eye icon next
>>>>> to the Programmable Source. Voila, you should have a line plot.
>>>>>
>>>>> Hope that helps,
>>>>> Cory
>>>>>
>>>>> On Wed, Sep 9, 2015 at 5:23 AM, no name <keinepostnurmuell at gmail.com>
>>>>> wrote:
>>>>>
>>>>>> Hey Cory
>>>>>>
>>>>>> Thanks for the hint! I did some coding and I got stuck at the
>>>>>> vtkTable as Output DataSet Type.
>>>>>>
>>>>>>
>>>>>> I used -> Programmable Source and then vtkTable which gives me the
>>>>>> following error code:
>>>>>>
>>>>>> ERROR: In
>>>>>> C:\DBD\pvs-x64\paraview\src\paraview\VTK\Common\ExecutionModel\vtkDemandDrivenPipeline.cxx,
>>>>>> line 810
>>>>>>
>>>>>> vtkPVDataRepresentationPipeline (000000009F2FC720): Input for
>>>>>> connection index 0 on input port index 0 for algorithm
>>>>>> vtkGeometryRepresentationWithFaces(00000000946A85F0) is of type vtkTable,
>>>>>> but a vtkDataSet is required.
>>>>>>
>>>>>>
>>>>>>
>>>>>> ERROR: In
>>>>>> C:\DBD\pvs-x64\paraview\src\paraview\VTK\Common\ExecutionModel\vtkDemandDrivenPipeline.cxx,
>>>>>> line 810
>>>>>>
>>>>>> vtkPVDataRepresentationPipeline (000000009F2FAB00): Input for
>>>>>> connection index 0 on input port index 0 for algorithm
>>>>>> vtkGeometryRepresentation(000000004FAC3570) is of type vtkTable, but a
>>>>>> vtkDataSet is required.
>>>>>>
>>>>>> for the code I am using:
>>>>>>
>>>>>> import numpy as np
>>>>>>
>>>>>>
>>>>>> switching_log = 'C:\\Users\\raphael\\Desktop\\logs auswertung\\newsims\\flat_Aex_1_E_-11\\flat_12\\switching.log'
>>>>>>
>>>>>> class DataFile(object):
>>>>>>
>>>>>>     def __init__(self, filename, h_xyz, m_xyz):
>>>>>>
>>>>>>         self._filename = filename
>>>>>>
>>>>>>         self._h_xyz = h_xyz
>>>>>>
>>>>>>         self._m_xyz = m_xyz
>>>>>>
>>>>>>                     def read_dat_file(self):
>>>>>>
>>>>>>         data = np.genfromtxt(self._filename, delimiter='', \
>>>>>>
>>>>>>         usecols=(0, self._h_xyz, self._m_xyz), unpack=True)
>>>>>>
>>>>>>         return data
>>>>>>
>>>>>>
>>>>>> log = DataFile(switching_log, 9, 13)
>>>>>>
>>>>>> data = log.read_dat_file()
>>>>>>
>>>>>>
>>>>>> A = np.rec.fromarrays(data, dtype= {'names': ['Time', 'Hx', 'Mx'], 'formats':['int','f8', 'f8']})
>>>>>>
>>>>>> print A.dtype.names
>>>>>>
>>>>>>
>>>>>> for name in A.dtype.names:
>>>>>>
>>>>>>     array = A[name]
>>>>>>
>>>>>>     output.RowData.append(array, name)
>>>>>>
>>>>>>
>>>>>> Now I could use the "vtkMultiBlockSet" or "vtkPolyData" as output
>>>>>> DataSet, but then I get:
>>>>>>
>>>>>>
>>>>>> Traceback (most recent call last):
>>>>>>
>>>>>> File "<string>", line 19, in <module>
>>>>>>
>>>>>> File "<string>", line 41, in RequestData
>>>>>>
>>>>>> File "C:\Program Files (x86)\ParaView
>>>>>> 4.3.1\lib\paraview-4.3\site-packages\vtk\numpy_interface\dataset_adapter.py",
>>>>>> line 121, in __getattr__
>>>>>>
>>>>>> return getattr(self.VTKObject, name)
>>>>>>
>>>>>> AttributeError: RowData
>>>>>>
>>>>>>
>>>>>> Because output.rowdata needs a vtkTable. I found the vtkTable in the
>>>>>>  Paraview Guide in Chapter 13.2.1 "Reading a CSV" file - which is basically
>>>>>> my thing, but I dont exactly understand where the difference between the
>>>>>> example and my situation lies ...
>>>>>>
>>>>>>
>>>>>> Thanks for your help!
>>>>>>
>>>>>>
>>>>>> Richard
>>>>>>
>>>>>> On 8 September 2015 at 15:00, Cory Quammen <cory.quammen at kitware.com>
>>>>>> wrote:
>>>>>>
>>>>>>> Hi Richard,
>>>>>>>
>>>>>>> You won't be able to set a vtkDataArray subclass as the input to a
>>>>>>> PlotOverLine filter. Instead, you should read the log files in a
>>>>>>> Programmable Source, and connect the Programmable Source to the
>>>>>>> PlotOverLine filter.
>>>>>>>
>>>>>>> See Chapter 13 of the ParaView Guide [1].
>>>>>>>
>>>>>>> Hope that helps,
>>>>>>> Cory
>>>>>>>
>>>>>>> [1]
>>>>>>> http://www.paraview.org/paraview-downloads/download.php?submit=Download&version=v4.3&type=data&os=all&downloadFile=TheParaViewGuide-v4.3-CC-Edition.pdf
>>>>>>>
>>>>>>>
>>>>>>> On Tue, Sep 8, 2015 at 3:57 AM, no name <keinepostnurmuell at gmail.com
>>>>>>> > wrote:
>>>>>>>
>>>>>>>> Hey there
>>>>>>>>
>>>>>>>> I have some hundred VTU files and I want to plot some data in a
>>>>>>>> graph (Plot over Line - POL) next to it. The Problem is, the data from
>>>>>>>> these POLs is coming from a simple log file, basically conaining only x and
>>>>>>>> y values (like two np.arrays).
>>>>>>>>
>>>>>>>> I can transform the arrays via
>>>>>>>>
>>>>>>>> from vtk.util import numpy_support
>>>>>>>> values = np.array((x,y))
>>>>>>>> vtk_values = numpy_support.numpy_to_vtk(values, array_type=vtk.VTK_FLOAT)
>>>>>>>> vtk_values.setName('values_2_show')
>>>>>>>> vtk_values.setComponentName(1, 'x_value')
>>>>>>>> vtk_values.setComponentName(2, 'y_value')
>>>>>>>>
>>>>>>>> But so far I have no success with how to change the input (source)
>>>>>>>>
>>>>>>>> VTUFILE = FindSource('VTUFILE*')
>>>>>>>> plotOverLine1 = PlotOverLine(Input=VTUFILE, Source='High Resolution Line Source')
>>>>>>>>
>>>>>>>> to a simple array.
>>>>>>>>
>>>>>>>> Thanks for any help!
>>>>>>>> Richard
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> Powered by www.kitware.com
>>>>>>>>
>>>>>>>> Visit other Kitware open-source projects at
>>>>>>>> http://www.kitware.com/opensource/opensource.html
>>>>>>>>
>>>>>>>> Please keep messages on-topic and check the ParaView Wiki at:
>>>>>>>> http://paraview.org/Wiki/ParaView
>>>>>>>>
>>>>>>>> Search the list archives at: http://markmail.org/search/?q=ParaView
>>>>>>>>
>>>>>>>> Follow this link to subscribe/unsubscribe:
>>>>>>>> http://public.kitware.com/mailman/listinfo/paraview
>>>>>>>>
>>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> --
>>>>>>> Cory Quammen
>>>>>>> R&D Engineer
>>>>>>> Kitware, Inc.
>>>>>>>
>>>>>>
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> Cory Quammen
>>>>> R&D Engineer
>>>>> Kitware, Inc.
>>>>>
>>>>
>>>>
>>>
>>
>>
>> --
>> Cory Quammen
>> R&D Engineer
>> Kitware, Inc.
>>
>
>

-- 
Cory Quammen
R&D Engineer
Kitware, Inc.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://public.kitware.com/pipermail/paraview/attachments/20150913/2ca049b3/attachment.html>