# error reading pvti in paraview
**URL:** https://discourse.paraview.org/t/error-reading-pvti-in-paraview/16229
**Category:** ParaView Support
**Tags:** vtk, python
**Created:** [February 24, 2025, 9:26am UTC](https://discourse.paraview.org/t/error-reading-pvti-in-paraview/16229 "2025-02-24T09:26:20Z")
**Posts on this page:** 1
**Page:** 1
### Author:  [@aliz2](https://discourse.paraview.org/u/aliz2)
#### Post date: [February 24, 2025, 9:26am UTC](https://discourse.paraview.org/t/error-reading-pvti-in-paraview/16229/1 "2025-02-24T09:26:20Z")
Hi there,
I’m trying to generate vtk image data in parallel (pvti) to read from series of vti files written by `vtk.vtkXMLImageDataWriter()` from 2d numpy arrays but getting below error related to extent:
```auto
ERROR: In vtkXMLPStructuredDataReader.cxx, line 429
vtkXMLPImageDataReader (0x9b7b0e0): No available piece provides data for the following extents:
0 511 0 511 0 10710
The UpdateExtent cannot be filled.
```
Tried the code with random generated numpy array to simplify but still giving same error:
```auto
ERROR: In vtkXMLPStructuredDataReader.cxx, line 429
vtkXMLPImageDataReader (0x13ba61e0): No available piece provides data for the following extents:
0 99 0 99 0 49
The UpdateExtent cannot be filled.
```
Not sure where I’m missing anything!
The code:
```auto
import vtk
import numpy as np
def numpy_to_vtk_image_data(numpy_array):
# Create a vtkImageData object
image_data = vtk.vtkImageData()
# Set the dimensions of the image data
height, width = numpy_array.shape
image_data.SetDimensions(width, height, 1)
# Allocate scalars for the image data
image_data.AllocateScalars(vtk.VTK_FLOAT, 1)
# Get a pointer to the image data's scalar values
scalars = image_data.GetPointData().GetScalars()
# Flatten the numpy array and set the scalar values
flat_numpy_array = numpy_array.flatten()
for i in range(len(flat_numpy_array)):
scalars.SetTuple1(i, flat_numpy_array[i])
return image_data
def save_vtk_image_data(image_data, filename):
# Create a vtkXMLImageDataWriter object
writer = vtk.vtkXMLImageDataWriter()
# Set the input data for the writer
writer.SetInputData(image_data)
# Set the filename for the output file
writer.SetFileName(filename)
# Write the data to the file
writer.Write()
def create_pvti_file(vti_filenames, pvti_filename, width, height):
with open(pvti_filename, 'w') as f:
f.write('\n')
f.write('\n')
f.write(' \n'.format(
width-1, height-1, len(vti_filenames)-1))
f.write(' \n')
f.write(' \n')
f.write(' \n')
f.write(' \n')
f.write(' \n')
for i, vti_filename in enumerate(vti_filenames):
f.write(' \n'.format(
width-1, height-1, i, i, vti_filename))
f.write(' \n')
f.write('\n')
# Example usage
if __name__ == " __main__":
# Create a series of 2D numpy arrays
series_of_2d_arrays = [np.random.rand(100, 100) for _ in range(50)]
# Save each 2D numpy array to a VTI file
vti_filenames = []
for i, array in enumerate(series_of_2d_arrays):
vtk_image_data = numpy_to_vtk_image_data(array)
vti_filename = f"slice_{i}.vti"
save_vtk_image_data(vtk_image_data, vti_filename)
vti_filenames.append(vti_filename)
# Create a PVTi file that references all the VTI files
create_pvti_file(vti_filenames, "output.pvti", series_of_2d_arrays[0].shape[1], series_of_2d_arrays[0].shape[0])
```
sample generated `vti`:
```auto
_AgAAAACAAABAHAAA *****
```
and the `pvti`:
```auto
```
Any clue what’s wrong?