Affine Transformationen mit scikit-image#

Dieses Notebook demonstriert, wie man affine Transformationen auf 3D-Bilder anwendet.

import numpy as np
import matplotlib
import matplotlib.pyplot as plt
from skimage.io import imread

# Laod example data
np_array = imread('../../data/Haase_MRT_tfl3d1.tif')
np_array.shape

(192, 256, 256)

Um eine affine Transformation einzurichten, können Sie dies mit einer 4x4-Transformationsmatrix tun:

transform_matrix = np.asarray([
    [1, 0, 0, 50],
    [0, 2, 0, 0],
    [0, 0, 0.5, 0],
    [0, 0, 0, 1]
])

Scikit-image unterstützt nur 2D-Transformationen, daher wählen wir eine Scheibe aus, um sie zu transformieren:

# pull image stack from GPU and pick a slice
image = np_array[100]

from skimage.io import imshow
imshow(image)

<matplotlib.image.AxesImage at 0x25a6cff0d30>
../_images/b9ec53895c9942461dd35e1b88b95a9557441ceb09c074bd865d8bb6ddfdb74b.png

Wir definieren jetzt eine affine Transformation mit scikit-image und wenden sie auf das Bild an.

from skimage import transform as tf

# define transform with #scikit image
transform = tf.AffineTransform(scale=0.5, translation=[10,0])

transformed_image = tf.warp(image, transform.inverse)
imshow(transformed_image)

<matplotlib.image.AxesImage at 0x25a6edb8d60>
../_images/ee46efb712e6a7f00d1c03698e27693eb7310b59b6c954469d8c2960bfca3103.png

Interoperabilität mit clesperanto#

Als Nächstes übertragen wir dieses einzelne Ebenenbild auf die GPU und transformieren es mit pyclesperanto

import pyclesperanto_prototype as cle

cle.select_device('RTX')

<gfx90c on Platform: AMD Accelerated Parallel Processing (2 refs)>

image_gpu = cle.push(image)

# define transform with #scikit image
from skimage import transform as tf
transform = tf.AffineTransform(scale=0.5, translation=[10,0])

transformed_image = cle.affine_transform(image_gpu, transform=transform)
cle.imshow(transformed_image, color_map="Greys_r")

c:\structure\code\pyclesperanto_prototype\pyclesperanto_prototype\_tier9\_imshow.py:14: UserWarning: The imshow parameter color_map is deprecated. Use colormap instead.
  warnings.warn("The imshow parameter color_map is deprecated. Use colormap instead.")
../_images/de3fcfb06177c0049426480e12bd738ddf30a227e64865662437057637224755.png