diskImageView.py

This code demonstrates how to create an image view on disk and write/read data to/from it.

This code demonstrates how to create an image view on disk and write/read data to/from it. The code creates a 3D image of size 100x100x100 without allocating any CPU memory, and writes a line of data in each slice. Finally, it then reads and prints the value of the last written sample.

import iolink
 
import numpy as np
 
 
def diskImageView():
    shape = iolink.VectorXu64(100, 100, 100)
    dt = iolink.DataTypeId.DOUBLE
 
    image = iolink.ImageViewFactory.create_on_disk(shape, dt)
 
    print("Writing some data slice by slice, line by line")
 
    lineSize = iolink.VectorXu64(shape[0], 1, 1)
    valueInit = 0
    for i in range(0, shape[2]):
        for j in range(0, shape[1]):
            lineRegion = iolink.RegionXu64((0, j, i), lineSize)
            # create a buffer corresponding to the line size, and
            # initialize it to the same value
            lineBuffer = np.full((shape[0], 1), valueInit, dtype=float)
            image.write_region(lineRegion, lineBuffer)
            valueInit+=1
 
    print("Writing completed")
 
    return image
 
 
image = diskImageView()
 
print("Reading last sample")
lastSampleValue = image.read((99, 99, 99))
print("Last sample value:", lastSampleValue)
 
print("SUCCESS")