[vtkusers] Fwd: AMR, ImageData and millions of files

[Peter Delmont]

Hello,

I`m extending the output of a fortran code, using Adaptive Mesh
Refinement, so that the output is .vtm metafile, type
vtkHierarchicalBoxDataSet.  Each of the datasets (which means: the
data of each of the grids) is saved in a different .vti-file.  Now,
that works fine.  Except for one small problem.  I`ve got millions of
grids, which means I also got millions of files.  Trying to cure this
problem, I thought that about putting lots of grids together in one
file, using the <piece>s.  My first idea was to put all data from the
same level in the same file.  The problem with this idea is that the
"Whole Extent" isn`t a rectangle anymore.  When I make the
<WholeExtent> a rectangle, which is not fully covered by the <piece>s,
There seem to happen strange things with the parts wich are not
covered by the <piece>s: They get random colours (while I thought that
the group argument would be used to look at the underlying grids, from
coarser levels.)

So my question is the following.  Is it necessary that the
<wholeExtent> is fully covered by its <piece Extent>s?  Is there a way  
around this problem?

I attached artificial .vtm and .vti files, where you can see the problem:
(The <WholeExtent> for the last file is not covered by all of its
<Piece Extent>s. The missing part seems to be coloured ramdomly.)

Greetings,
and thanks in advance,
(any reaction can be helpful)




Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm

-------------- next part --------------
<?xml version="1.0"?>
<VTKFile type ="vtkHierarchicalBoxDataSet" version="0.1" byte_order="LittleEndian">
<vtkHierarchicalDataSet>
<RefinementRatio level="0" refinement="2"/>
<DataSet group="0" dataset="0" amr_box="0 3 0 3 0 0" file="lasthope_0.vti"/>
<DataSet group="1" dataset="0" amr_box="0 3 4 7 0 0" file="lasthope_1.vti"/>
<DataSet group="1" dataset="1" amr_box="4 7 4 7 0 0" file="lasthope_1.vti"/>
<DataSet group="1" dataset="2" amr_box="4 7 0 3 0 0" file="lasthope_1.vti"/>
</vtkHierarchicalDataSet>
</VTKFile>
-------------- next part --------------
<?xml version="1.0"?>
<VTKFile type="ImageData" version="0.1" byte_order="LittleEndian">
<ImageData WholeExtent="0 4 0 4 0 0" Origin="0 0 0" Spacing="1 1 0">
<Piece Extent="0 4 0 4 0 0">
<PointData>
</PointData>
<CellData>
<DataArray type= "Float32" Name="nonsens" format="ascii">
0 1 2 3 
1 2 3 4 
2 3 4 5 
3 4 5 6 
</DataArray>
</CellData>
</Piece>
</ImageData>
</VTKFile>
-------------- next part --------------
<?xml version="1.0"?>
<VTKFile type="ImageData" version="0.1" byte_order="LittleEndian">
<ImageData WholeExtent="0 8 0 8 0 0" Origin="0 0 0" Spacing="0.5 0.5 0">
<Piece Extent="0 4 4 8 0 0">
<PointData>
</PointData>
<CellData>
<DataArray type= "Float32" Name="nonsens" format="ascii">
2 3 4 3
3 4 5 4 
4 5 4 3
3 4 3 2
</DataArray>
</CellData>
</Piece>
<Piece Extent="4 8 4 8 0 0">
<PointData>
</PointData>
<CellData>
<DataArray type= "Float32" Name="nonsens" format="ascii">
4 5 6 4
5 6 7 6 
6 7 6 5 
5 6 5 4 
</DataArray>
</CellData>
</Piece>
</ImageData>
</VTKFile>