import numpy as np
from skimage import io

# read the image stack
fname = '100x120x150.tif'
data = io.imread(fname)
data.shape