Cannot create IJKCoordinateTransformer for 2D dataset (Python)
Description
I am currently playing around with the creation of VDS files. I am especially interested in working via the Python interface and the different coordinate systems. I set up the a small Python script that creates a simple 2D dataset from a NumPy array with random content. Parts of the script is based on the npz_to_vds.py
script from the examples. I would like to convert between inline/crossline coordinates, voxel coordinates and world coordinates.
In my script, the creation of the VDS file is successful. I also see that the file is recognized as 2D file by OpenVDS during writing since the chunks written to the page buffer are 4*brick_size. However, when I want to obtain the IJKCoordinateTransformer
for this file, I run into the following exception
Exception:
Dimension -1 is not a valid dimension. Dimensionality_Max is 6.
When I create a 3D file with only one coordinate in z direction, obtaining the transformer seems to be successful.
Expectation
I obtain the coordinate transformer which allows me to transform between different coordinate systems (ijk, inline/crossline etc.).
Questions
- Is this behavior expected? I assumed that I could still work with the IJK transformer.
- Do I create the file as a 2D file in a wrong way?
- If the behavior is expected, would it be possible to extend the transformer to work with 2D data and/or as a quick fix to make the error message more expressive?
System
- Arm64 MacOS 12.6
- VDS 3.0.3 with Python interface