| Thread Previous • Date Previous • Date Next • Thread Next |
Question #147820 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/147820
Patrick Riesen proposed the following answer:
in python, I use this to read in a vtu file to an array:
this doesn't need the GetArray("name")
# -------------------------------------------------
import vtk as v
from enthought.tvtk import array_handler as ah
datareader = v.vtkXMLUnstructuredGridReader()
datareader.SetFileName(filename)
datareader.Update()
data = datareader.GetOutput()
# if you have a vector (e.g. velocity):
vtk_values = data.GetPointData().GetVectors()
# if you have a scalar:
# vtk_values = data.GetPointData().GetScalars()
# this is convenient to convert vtkarrays to numpy arrays
solution = ah.vtk2array(data)
# ------------------------------------------------------
regards,
patrick
On 03/04/2011 09:30 PM, JJ Jiang wrote:
> Question #147820 on DOLFIN changed:
> https://answers.launchpad.net/dolfin/+question/147820
>
> Status: Answered => Open
>
> JJ Jiang is still having a problem:
> Hi Johan and Kevin,
> Thanks for the notes. I think you can dump simulated temperature or mass to raw files. That can be manipulated when the mesh data are also given either in a vtu format or in a raw format.
> By the way, I used Paraview to see the simulated "temperature" scalar has been attached to a vtu. Surprisingly, I couldn't read it out by its name. I don't know why.
> I guess Johan's suggestions will give users more flexibilities to handle the data.
> Thanks again to both of you.
> Best regards,
> --JJ
>
> On 03/04/11, Johan Hake wrote:
>> Your question #147820 on DOLFIN changed:
>> https://answers.launchpad.net/dolfin/+question/147820
>>
>> Johan Hake proposed the following answer:
>> Maybe Viper could provide a high level interface to read a vtu file into a
>> NumPy array, and a mesh into a coordinate and cell arrays?
>>
>> Then we eventually can use this to pass data to Viper for plotting in the C++
>> interface. This might fix the parallel ploting issue we have.
>>
>> It would also be cool to implement a helper function in the Python interface
>> of DOLFIN that uses the potential added functionality in Viper to generate a
>> Function, with or without a given mesh. A FunctionSpace could be automatically
>> generated based on the dimensions of the data (is this true?) If a mesh is not
>> given we generate it based on the stored mesh data from the file.
>>
>> Johan
>>
>>
>> On Friday March 4 2011 11:30:52 Kent-Andre Mardal wrote:
>>> Question #147820 on DOLFIN changed:
>>> https://answers.launchpad.net/dolfin/+question/147820
>>>
>>> Status: Open => Answered
>>>
>>> Kent-Andre Mardal proposed the following answer:
>>> You should be able to use GetPointData().GetArray(). Have you looked at the
>>> name of the variable in the vtu file? It is in ASCII format.
>>>
>>> I attach something I did earlier.
>>>
>>> Kent
>>>
>>> On 4 March 2011 19:34, JJ Jiang
>>>
>>> <question147820@xxxxxxxxxxxxxxxxxxxxx>wrote:
>>>> New question #147820 on DOLFIN:
>>>> https://answers.launchpad.net/dolfin/+question/147820
>>>>
>>>> Hi all,
>>>>
>>>> I am running into a problem. I run the advection-diffusion demo. All the
>>>>
>>>> simulated mass transport quantities u were save to a vtk PVD file with an
>>>> unstructured mesh.
>>>>
>>>> My question is how to read out the u from the associated individual vtu
>>>> files. I thought the u is attached to either PointData or CellData so
>>>> that I can read out the scalar u using one of the following:
>>>> grid->GetPointData()->GetArray('name');
>>>> or
>>>> grid->GetCellData()->GetArray('name');
>>>>
>>>> But my attempts have failed. Please provide your advice. Thanks so much
>>>> for your help!
>>>>
>>>> Best regards,
>>>> --JJ
>>>>
>>>>
>>>> --
>>>> You received this question notification because you are a member of
>>>> DOLFIN Team, which is an answer contact for DOLFIN.
>>>>
>>>> _______________________________________________
>>>> Mailing list: https://launchpad.net/~dolfin
>>>> Post to : dolfin@xxxxxxxxxxxxxxxxxxx
>>>> Unsubscribe : https://launchpad.net/~dolfin
>>>> More help : https://help.launchpad.net/ListHelp
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>> https://answers.launchpad.net/dolfin/+question/147820/+confirm?answer_id=1
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/dolfin/+question/147820
>>
>> You received this question notification because you are a direct
>> subscriber of the question.
>
> --
> +++++++++++++++++++++++++++++++++++
> Jingfeng Jiang (JJ)
> University of Wisconsin
> Medical Physics Department
> 1530 Medical Sciences Center
> 1300 University Avenue
> Madison, WI 53706
> Phone: (608) 262-4197
> FAX: (608) 262-2413
> +++++++++++++++++++++++++++++++++++
>
> You received this question notification because you are a member of
> DOLFIN Team, which is an answer contact for DOLFIN.
>
> _______________________________________________
> Mailing list: https://launchpad.net/~dolfin
> Post to : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help : https://help.launchpad.net/ListHelp
--
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.
| Thread Previous • Date Previous • Date Next • Thread Next |
