Re: [Paraview] Writing XML VTK Binary files from C++

[Mark Olesen]

When you say "writing binary", you need to distinguish between three 
possibilities.

1) writing binary "inline" (actually base64 encoded)
2) writing binary "append" (actually base64 encoded)
3) writing raw binary "append" (really raw binary)

Since you already have ASCII writing working and its content is 
"inline", it won't take much more to get binary inline working.
By "inline", I mean when the output is placed between the   
markers. Eg,


content


For the binary case, the content is written as base64-encoded data, 
which means that your output writer for these sections needs to pass the 
content through a base64 layer to do the encoding for you.


If it helps as reference, we have the same thing in OpenFOAM, except 
that we only write vtu and vtp files (we don't have rectilinear meshes).


In the repo https://develop.openfoam.com/Development/OpenFOAM-plus
we have a foamVtkBase64Formatter and a foamVtkBase64Layer (both under 
src/fileFormats/vtk/format/) that add a base64Layer to encode and output 
as base64 (src/OpenFOAM/db/IOstreams/hashes/base64Layer.[CH]).


You'll see that the foamVtkBase64Layer and base64Layer are quite low 
level means of adding an tiny encoding buffer (3 chars size) to 
intercept output prior to sending through to a std::ostream. It take 
very little effort to adopt for your output and thus quite easy to drop 
in instead of your current ASCII outputter.  For it too work easily, 
however, you should make sure that you need to generate your output 
content with a write() method instead of using '<<'. This allows 
somewhat easy switching between something like a foamVtkAsciiFormatter 
and the binary version, but more importantly it makes it easier to track 
the output state.


When browsing through the code, you may also notice that we have support 
for writing in appended format (raw and base64). However, I would not 
advise you to tackle that immediately. There are a few more things to 
watch out for here, but more importantly it will change many more things 
on the calling side.


I hope this information helps you.
/mark

--
Dr Mark OLESEN
Principal Engineer, ESI-OpenCFD
ESI GmbH | Einsteinring 24 | 85609 Munich | GERMANY
www.openfoam.com | www.esi-group.com | mark.ole...@esi-group.com


On 11/30/17 18:00, Stegmeier, Nicholas wrote:

Hello,

I am new to Paraview and C++ coming from a mostly mathematics 
background. I am emailing to get resources or help on writing binary XML 
VTK files from C++.


I have finally succeeded in using the ASCII XML VTK format for a 2D 
rectilinear CFD application. My ".pvtr" file is shown below.


How can I write this file and my other XML VTK files in binary from C++? 
Do I need a special C++ library?

___
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

Re: [Paraview] unable to RGB voxel from tiff/jpeg slice to volume

[Cory Quammen]

Joe,

This is more of a question for the ITK discussion forum [1], but I'll
give you a potential solution here.

The problem is that

typedef itk::Image  ImageType;

declares a 3-dimensional scalar-valued image, not a 3-component (RGB)
2-dimensional image. To declare an RGB ITK image, try this:

typedef itk::RGBPixel< unsigned char > RGBPixelType;
typedef itk::Image< RGBPixelType, 3 > ImageType;

I hope that helps.

Cory

[1] https://discourse.itk.org/

On Thu, Nov 30, 2017 at 8:14 PM, joe kozak  wrote:
> grabbed frames from a standard youtube color video via ffmpeg
> used something like $./imageSeriesToVolume -sz 1.0 -o test.vtk -i
> input*.jpeg
> loaded test.vtk in paraview 5.0.1 and also paraview 5.4.1
> on information tab i see:   scalars  unsigned char [33,239]
>
> "Map Scalars"  on or off I get no RGB colored voxels, of course, because
> there is only one datum of type unsigned char.
>
> Where are the other two scalars?!?
>
> The important part you will ask me is here:
>
>  91   io->SetFileName (filenames[0].c_str());
>  92   io->ReadImageInformation();
>  93
>  94   if (io->GetNumberOfDimensions()==2) {
>  95
>  96 typedef itk::Image  ImageType;
>  97 typedef itk::ImageSeriesReader SeriesReaderType;
>  98
>  99
> 100 SeriesReaderType::Pointer reader = SeriesReaderType::New();
> 101 std::cout << "Adding:\n";
> 102 for (unsigned int i=0; i 103 {
> 104   std::cout << filenames[i] << std::endl;
> 105   reader->AddFileName ( filenames[i].c_str() );
> 106 }
> 107
>
> Looks pretty standard, but no 3 scalar datums for an RGB voxel for
> non-false-colored volume.
> It is perported to work, and all other examples i see are none different
> from the code i used to generate a .vtk file.
>
>
> I cannot use IsoVolume If I load the jpegs directly in paraview,  but I can
> get RGB color voxels.
> The GOAL is to use IsoVolume on a volume of RGB points from frames from a
> video sequence.
> I have tried on and off for MANY YEARS to figure this out...(its for an
> art project, metalcasting &3d Printing)
>
> ___
> 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
Staff R Engineer
Kitware, Inc.
___
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

[Paraview] unable to RGB voxel from tiff/jpeg slice to volume

[joe kozak]

grabbed frames from a standard youtube color video via ffmpeg
used something like $./imageSeriesToVolume -sz 1.0 -o test.vtk -i 
input*.jpeg
loaded test.vtk in paraview 5.0.1 and also paraview 5.4.1
on information tab i see:   scalars  unsigned char [33,239]

"Map Scalars"  on or off I get no RGB colored voxels, of course, because there 
is only one datum of type unsigned char.

Where are the other two scalars?!?

The important part you will ask me is here:

 91   io->SetFileName (filenames[0].c_str());
 92   io->ReadImageInformation();
 93
 94   if (io->GetNumberOfDimensions()==2) {
 95
 96 typedef itk::Image  ImageType;
 97 typedef itk::ImageSeriesReader SeriesReaderType;
 98
 99
100 SeriesReaderType::Pointer reader = SeriesReaderType::New();
101 std::cout << "Adding:\n";
102 for (unsigned int i=0; iAddFileName ( filenames[i].c_str() );
106 }
107

Looks pretty standard, but no 3 scalar datums for an RGB voxel for 
non-false-colored volume.
It is perported to work, and all other examples i see are none different from 
the code i used to generate a .vtk file.


I cannot use IsoVolume If I load the jpegs directly in paraview,  but I can get 
RGB color voxels.
The GOAL is to use IsoVolume on a volume of RGB points from frames from a video 
sequence.
I have tried on and off for MANY YEARS to figure this out...(its for an art 
project, metalcasting &3d Printing)
___
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

Re: [Paraview] post-process of CFL3D result with surfaces in paraview?

[kenichiro yoshimi]

Hi Cheng,

You can specify the i-j-k min/max indices in a structured grid dataset
to extract curviliner grid as wall using the vtkExtractGrid filter
through Programmable. Specifically, after adding Programmable Filter
and Changing the output type to vtkPolyDatat, there is something like
the following as Script:
-
import vtk

mdi = self.GetInput()
sdo = self.GetOutput()
append = vtk.vtkAppendPolyData()

blockNo = 0
sg = mdi.GetBlock(blockNo)
ext = list(sg.GetExtent())
print ext

ext[0] = 0
ext[1] = 0
extractor = vtk.vtkExtractGrid()
extractor.SetInputData(sg)
extractor.SetVOI(ext)

surface0 = vtk.vtkGeometryFilter()
surface0.SetInputConnection(extractor.GetOutputPort())
surface0.Update()
append.AddInputData(surface0.GetOutput())

ext = list(sg.GetExtent())
ext[2] = 0
ext[3] = 0
extractor.SetVOI(ext)
extractor.Modified()
surface1 = vtk.vtkGeometryFilter()
surface1.SetInputConnection(extractor.GetOutputPort())
surface1.Update()
append.AddInputData(surface1.GetOutput())

append.Update()

sdo.ShallowCopy(append.GetOutput())
-

Thanks

2017-11-30 5:32 GMT+09:00 程迪 :
> Hi paraviewers,
>
> I am using paraview to render results in CFL3D, which is in Plot3D format,
> including `abc.xyz` grid file and `abc.q` solution file. Besides, I have a
> boundary condition file `abc.xyz.fvbnd` in FieldView 1.3 format which was
> generated by grid pre-processor.
>
> My problem is: How to extract the surface of the aircraft?
>
> Because the Plot3D format is a multi-block curvilinear grid. When I use
> `Extract Surface` filter, it will extract the block interfaces. My current
> solution to combine blocks by `Merge Blocks` filter. Then I can extract the
> surfaces, including the outer boundary and the inner aircraft walls. After
> that, I need to clip the farfield and keep the aircraft walls. Is there a
> better solution?
>
> However, what I need to do is to evaluate the force on specific parts. So I
> need to use the information from the boundary condition file `abc.xyz.fvbnd`
> which indicates ranges of indexes in curvilinear grid as wall of certain
> parts. My problem is how to use it to generate several surfaces in paraview
> for rendering and integration.
>
>
> p.s. I am using paraview 5.4.1. The grid consists of farfield outer boundary
> and several aircraft parts' wall inner boundaries.
>
>
> Di Cheng
> Engineer of Research Center
> China Academy of Aerospace Aerodynamics
> Phone: +86-l58Ol5949ll
>
> Address: No.17, YunGang West Road, Fengtai District, Beijing, China
>
> Zip Code：100074
>
>
> ___
> 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
>
___
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

[Paraview] Writing XML VTK Binary files from C++

[Stegmeier, Nicholas]

Hello,


I am new to Paraview and C++ coming from a mostly mathematics background. I am 
emailing to get resources or help on writing binary XML VTK files from C++.


I have finally succeeded in using the ASCII XML VTK format for a 2D rectilinear 
CFD application. My ".pvtr" file is shown below.


How can I write this file and my other XML VTK files in binary from C++? Do I 
need a special C++ library?


filename: PVTK.pvtr


contents:





















Thank you,
Nicholas Stegmeier


___
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

Re: [Paraview] Assemble Variables into Tensor

[Utkarsh Ayachit]

You're most welcome :).

Utkarsh

On Thu, Nov 30, 2017 at 10:20 AM, Nicholas Richmond
 wrote:
> Hi Utkarsh,
>
> Using the script you wrote with the Programmable Filter does the job nicely!
> Thanks so much for going above and beyond to help me assemble the arrays
> into a tensor. Programmable Filter never occurred to me and I am now looking
> forward to exploring Chapter 13 of the Paraview Guide to see what else I can
> do with it.
>
> Many thanks,
> Nick Richmond
>
> On Wed, Nov 29, 2017 at 9:56 PM, Utkarsh Ayachit
>  wrote:
>>
>> Nick,
>>
>> You can use the Programmable Filter instead with something like the
>> following script as the Script:
>>
>>
>> from paraview.vtk.numpy_interface import dataset_adapter as dsa
>>
>> import numpy
>>
>> def make_tensor(xx,yy,zz, xy, yz, xz):
>>
>>   t = numpy.vstack([xx,yy,zz,xy, yz, xz]).transpose().view(dsa.VTKArray)
>>
>>   t.DataSet = xx.DataSet
>>
>>   t.Association = xx.Association
>>
>>   return t
>>
>>
>> xx = inputs[0].PointData["sigma_xx"]
>>
>> yy = inputs[0].PointData["sigma_yy"]
>>
>> zz = inputs[0].PointData["sigma_zz"]
>>
>> xy = inputs[0].PointData["sigma_xy"]
>>
>> yz = inputs[0].PointData["sigma_yz"]
>>
>> xz = inputs[0].PointData["sigma_xz"]
>>
>> output.PointData.append(make_tensor(xx,yy,zz,xy,yz,xz), "tensor")
>>
>>
>> Utkarsh
>>
>> On Wed, Nov 29, 2017 at 5:01 PM, Nicholas Richmond
>>  wrote:
>>>
>>> Greetings,
>>>
>>> I have a vtk Polygonal Mesh with six separate arrays (sigma_xx, sigma_yy,
>>> sigma_zz, sigma_xy, sigma_yz, sigma_xz). I'd like to assemble these into a
>>> tensor so that I may apply the Tensor Glyph filter.
>>>
>>> After reading that the paraview.simple.TensorGlyph expects a symmetric
>>> tensor with the order XX, YY, ZZ, XY, YZ, XZ, I tried to assemble the six
>>> separate arrays into a single array with the Python Calculator, but have
>>> been unsuccessful.
>>>
>>> I've tried a variety of approaches, including:
>>>
>>> Expression: np.array([sigma_xx, sigma_yy, sigma_zz, sigma_xy, sigma_yz,
>>> sigma_xz])
>>> Array Association: Point Data
>>>
>>> I can make a vector of the normal stresses using the "make_vector"
>>> function from paraview.vtk.numpy_interface.algorithms, but there is no
>>> "make_tensor" function to assemble both normal and shear stresses into a
>>> single symmetric tensor.
>>>
>>> Any guidance would be most appreciated.
>>> Many thanks,
>>> Nick Richmond
>>>
>>> ___
>>> 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
>>>
>>
>
___
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

Re: [Paraview] Assemble Variables into Tensor

[Nicholas Richmond]

Hi Utkarsh,

Using the script you wrote with the Programmable Filter does the job
nicely! Thanks so much for going above and beyond to help me assemble the
arrays into a tensor. Programmable Filter never occurred to me and I am now
looking forward to exploring Chapter 13 of the Paraview Guide to see what
else I can do with it.

Many thanks,
Nick Richmond

On Wed, Nov 29, 2017 at 9:56 PM, Utkarsh Ayachit <
utkarsh.ayac...@kitware.com> wrote:

> Nick,
>
> You can use the *Programmable Filter* instead with something like the
> following script as the *Script*:
>
>
> from paraview.vtk.numpy_interface import dataset_adapter as dsa
>
> import numpy
>
> def make_tensor(xx,yy,zz, xy, yz, xz):
>
>   t = numpy.vstack([xx,yy,zz,xy, yz, xz]).transpose().view(dsa.VTKArray)
>
>   t.DataSet = xx.DataSet
>
>   t.Association = xx.Association
>
>   return t
>
>
> xx = inputs[0].PointData["sigma_xx"]
>
> yy = inputs[0].PointData["sigma_yy"]
>
> zz = inputs[0].PointData["sigma_zz"]
>
> xy = inputs[0].PointData["sigma_xy"]
>
> yz = inputs[0].PointData["sigma_yz"]
>
> xz = inputs[0].PointData["sigma_xz"]
>
> output.PointData.append(make_tensor(xx,yy,zz,xy,yz,xz), "tensor")
>
>
> Utkarsh
>
> On Wed, Nov 29, 2017 at 5:01 PM, Nicholas Richmond <
> nicholas.richm...@maine.edu> wrote:
>
>> Greetings,
>>
>> I have a vtk Polygonal Mesh with six separate arrays (sigma_xx, sigma_yy,
>> sigma_zz, sigma_xy, sigma_yz, sigma_xz). I'd like to assemble these into a
>> tensor so that I may apply the Tensor Glyph filter.
>>
>> After reading that the paraview.simple.TensorGlyph
>> 
>>  expects a symmetric tensor with the order XX, YY, ZZ, XY, YZ, XZ, I
>> tried to assemble the six separate arrays into a single array with the
>> Python Calculator, but have been unsuccessful.
>>
>> I've tried a variety of approaches, including:
>>
>> *Expression:* np.array([sigma_xx, sigma_yy, sigma_zz, sigma_xy,
>> sigma_yz, sigma_xz])
>> *Array Association:* Point Data
>>
>> I can make a vector of the normal stresses using the "make_vector"
>> function from paraview.vtk.numpy_interface.algorithms, but there is no
>> "make_tensor" function to assemble both normal and shear stresses into a
>> single symmetric tensor.
>>
>> Any guidance would be most appreciated.
>> Many thanks,
>> Nick Richmond
>>
>> ___
>> 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
>>
>>
>
___
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

Re: [Paraview] particle traces from point data

[Alex Barrie]

Yes, I have one file per timestep and each file is a list of points, X, Y,
Z, C with one point per line. I am loading it with the 'Point3D' input.
Thanks!
Alex


On Thu, Nov 30, 2017 at 9:18 AM, Joachim Pouderoux <
joachim.pouder...@kitware.com> wrote:

> Alex,
>
> Could you specify how is your data split over time? Do you have one file
> per timestep? How do you load your data with ParaView?
>
> Joachim
>
> *Joachim Pouderoux*, PhD
>
> *Technical Expert - Scientific Computing Team*
> *Kitware SAS *
>
>
> 2017-11-30 9:55 GMT-04:00 Alex Barrie :
>
>> Thanks, I see it now - I missed it because it wasn't in the "Temporal"
>> tab.
>>
>> When I apply this filter, however, I get an error about a missing
>> DATA_TIME_STEPS information key. Is this something I can specify or does it
>> have to be in the raw data somehow?
>>
>> Thanks,
>> Alex
>>
>>
>> On Wed, Nov 29, 2017 at 3:13 PM, Joachim Pouderoux <
>> joachim.pouder...@kitware.com> wrote:
>>
>>> Alex,
>>>
>>> I think here the filter "Temporal Particles to Pathlines" should do the
>>> trick.
>>>
>>> Best,
>>> Joachim
>>>
>>> *Joachim Pouderoux*, PhD
>>>
>>> *Technical Expert - Scientific Computing Team*
>>> *Kitware SAS *
>>>
>>>
>>> 2017-11-29 14:33 GMT-04:00 Alex Barrie :
>>>
 Hi,

 I have animated particles which are imported from point data (x,y,z). I
 want to add a particle trace, but I am not sure what to put for the 'seed'.
 I am not using a vector field, but rather just a list of particle positions
 that updates over time. How can I just add a trace to each particle in the
 list?

 Thanks,
 Alex

 ___
 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


>>>
>>
>
___
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

Re: [Paraview] particle traces from point data

[Joachim Pouderoux]

Alex,

Could you specify how is your data split over time? Do you have one file
per timestep? How do you load your data with ParaView?

Joachim

*Joachim Pouderoux*, PhD

*Technical Expert - Scientific Computing Team*
*Kitware SAS *


2017-11-30 9:55 GMT-04:00 Alex Barrie :

> Thanks, I see it now - I missed it because it wasn't in the "Temporal"
> tab.
>
> When I apply this filter, however, I get an error about a missing
> DATA_TIME_STEPS information key. Is this something I can specify or does it
> have to be in the raw data somehow?
>
> Thanks,
> Alex
>
>
> On Wed, Nov 29, 2017 at 3:13 PM, Joachim Pouderoux <
> joachim.pouder...@kitware.com> wrote:
>
>> Alex,
>>
>> I think here the filter "Temporal Particles to Pathlines" should do the
>> trick.
>>
>> Best,
>> Joachim
>>
>> *Joachim Pouderoux*, PhD
>>
>> *Technical Expert - Scientific Computing Team*
>> *Kitware SAS *
>>
>>
>> 2017-11-29 14:33 GMT-04:00 Alex Barrie :
>>
>>> Hi,
>>>
>>> I have animated particles which are imported from point data (x,y,z). I
>>> want to add a particle trace, but I am not sure what to put for the 'seed'.
>>> I am not using a vector field, but rather just a list of particle positions
>>> that updates over time. How can I just add a trace to each particle in the
>>> list?
>>>
>>> Thanks,
>>> Alex
>>>
>>> ___
>>> 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
>>>
>>>
>>
>
___
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

Re: [Paraview] particle traces from point data

[Alex Barrie]

Thanks, I see it now - I missed it because it wasn't in the "Temporal" tab.

When I apply this filter, however, I get an error about a missing
DATA_TIME_STEPS information key. Is this something I can specify or does it
have to be in the raw data somehow?

Thanks,
Alex


On Wed, Nov 29, 2017 at 3:13 PM, Joachim Pouderoux <
joachim.pouder...@kitware.com> wrote:

> Alex,
>
> I think here the filter "Temporal Particles to Pathlines" should do the
> trick.
>
> Best,
> Joachim
>
> *Joachim Pouderoux*, PhD
>
> *Technical Expert - Scientific Computing Team*
> *Kitware SAS *
>
>
> 2017-11-29 14:33 GMT-04:00 Alex Barrie :
>
>> Hi,
>>
>> I have animated particles which are imported from point data (x,y,z). I
>> want to add a particle trace, but I am not sure what to put for the 'seed'.
>> I am not using a vector field, but rather just a list of particle positions
>> that updates over time. How can I just add a trace to each particle in the
>> list?
>>
>> Thanks,
>> Alex
>>
>> ___
>> 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
>>
>>
>
___
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

Re: [Paraview] Best file fomat for big data

[Chris Richardson]

I would recommend using xdmf (www.xdmf.org). It works well for large parallel datasets, and supports the same structures as vtk.ChrisOn 30 Nov 2017 10:13, a...@yandex.ru wrote:Hi! I just finish my own CFD simulation on 256 processors. The program made graphics output in 2000 time points. So large number of outputs I need to prepare smooth animation. I use XML VTK ascii format. So I have now 2000 pvtu files and 2000*256=512000 vtu files. I have problems with such big number of files. Now I understand that may be it was not the best idea to use VTK format. I have the one fixed unstructured grid with ~10^7 cells, about ten CFD variables on the cells and several thousands output points. My program is MPI parallel and output time points are irregular over time. So I'd like to keep value of time in the files (as I know it is impossible in the VTK format). And I'd like to use ParaView for the visualization. What is the better format for storing big amount of CFD data that I be able to write in parallel and load in the ParaView? Thank in advance! 
___
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

Re: [Paraview] Best file fomat for big data

[Mark Olesen]

A few months back there was a suggestion on the mailing list that you 
can apparently write a single vtk geometry and then use something like 
"map data" to load the new, updated fields (perhaps Ken posted it? - 
don't know any other details).


I agree that these vtk format limitations can be quite frustrating and 
cause quite a bit of data bloat on the disk. For these situations, I 
would normally resort to using EnSight format. This gives you full 
separation of geometry and data, with moving or non-moving geometries 
etc. EnSight does support a server-of-servers data format, which makes 
it possible to a sub-case for each process and then load them together, 
but I've never used it. Instead we have the master process write a 
single file and stream the geometry and fields information across.
The additional overhead of moving data across the network is offset by 
the fact that the filer only has a single process writing to it.


If it helps, you can take a look at the foamToEnsight utility for a 
quick overview: https://develop.openfoam.com
The harder bits are located under src/conversion/ensight/mesh/, 
src/conversion/ensight/output/, src/conversion/ensight/part/


If you find another possibility, please let us know since it will be of 
interest for other people too.


Cheers,
/mark

--
Dr Mark OLESEN
Principal Engineer, ESI-OpenCFD
ESI GmbH | Einsteinring 24 | 85609 Munich | GERMANY
Mob. +49 171 9710 149
www.openfoam.com | www.esi-group.com | mark.ole...@esi-group.com


On 11/30/17 11:13, a...@yandex.ru wrote:
Hi! I just finish my own CFD simulation on 256 processors. The program 
made graphics output in 2000 time points. So large number of outputs I 
need to prepare smooth animation. I use XML VTK ascii format. So I have 
now 2000 pvtu files and 2000*256=512000 vtu files. I have problems with 
such big number of files. Now I understand that may be it was not the 
best idea to use VTK format.
I have the one fixed unstructured grid with ~10^7 cells, about ten CFD 
variables on the cells and several thousands output points. My program 
is MPI parallel and output time points are irregular over time. So I'd 
like to keep value of time in the files (as I know it is impossible in 
the VTK format). And I'd like to use ParaView for the visualization.
What is the better format for storing big amount of CFD data that I be 
able to write in parallel and load in the ParaView?

Thank in advance!

___
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

[Paraview] Best file fomat for big data

[as92]

Hi! I just finish my own CFD simulation on 256 processors. The program made graphics output in 2000 time points. So large number of outputs I need to prepare smooth animation. I use XML VTK ascii format. So I have now 2000 pvtu files and 2000*256=512000 vtu files. I have problems with such big number of files. Now I understand that may be it was not the best idea to use VTK format. I have the one fixed unstructured grid with ~10^7 cells, about ten CFD variables on the cells and several thousands output points. My program is MPI parallel and output time points are irregular over time. So I'd like to keep value of time in the files (as I know it is impossible in the VTK format). And I'd like to use ParaView for the visualization. What is the better format for storing big amount of CFD data that I be able to write in parallel and load in the ParaView? Thank in advance! 
___
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

[Paraview] Registered texture images

[postgurke]

Hello,

 

I am using python scripts to import data into a ParaView (5.4.1) session. In order to keep the texture images in a saved state, I am registering them:

 


pathToTextureImage='path\\imagefile'
Display = Show(testvtp, renderView1)
Display.SetRepresentationType('Surface')
texProxy = servermanager.CreateProxy('textures','ImageTexture')
texProxy.GetProperty('FileName').SetElement(0, pathToTextureImage)
texProxy.UpdateVTKObjects()
servermanager.Register(texProxy, registrationName='ImageName')

 

Now, when I call this script multiple times (e.g. because I made changes to the vtp defining the texture plane), I will have multiple instances of my image in the pull-down menu for the texture. Is there a way to avoid this, either by checking if the resitrationName already exists or by deleting the multiple texture image from the list afterwards?

 

I have found this https://public.kitware.com/pipermail/paraview/2007-October/006126.html , but I cannot figure out how to use the code correctly in my case.

 

Also, is it possible to show only the images associated with the respective object in the pull-down menu? When I have many profiles with images loaded, the menu gets quite crowded and there is a risk of choosing an image not belonging to the current profile. 

 

Any suggestions on how to solve this are very much apprechiated. Thanks in advance!

Cheers,

Venke

 

 

 

 

___
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