Задайте упакованный (x, y, z) координаты пяти точек ввода. Упакованные координаты хранятся как 5 3 матрица, где первые, вторые, и третьи столбцы содержат x-, y-, и координаты z-, соответственно.

XYZ = [5 25 20;10 5 25;15 10 5;20 15 10;25 20 15];

Задайте функцию обратного отображения, которая принимает и возвращает точки в упакованном (x, y, z) формат.

inverseFcn = @(c) [c(:,1)+c(:,2),c(:,1)-c(:,2),c(:,3).^2];

Создайте 3-D объект геометрического преобразования, tform, который хранит эту функцию обратного отображения.

tform = geometricTransform3d(inverseFcn)

tform = 
  geometricTransform3d with properties:

        InverseFcn: @(c)[c(:,1)+c(:,2),c(:,1)-c(:,2),c(:,3).^2]
        ForwardFcn: []
    Dimensionality: 3

Примените обратное преобразование этого 3-D геометрического преобразования к точкам ввода.

UVW = transformPointsInverse(tform,XYZ)

UVW = 5×3

    30   -20   400
    15     5   625
    25     5    25
    35     5   100
    45     5   225

Задайте x-, y-и векторы z-координаты пяти точек, чтобы преобразовать.

x = [3 5 7 9 11];
y = [2 4 6 8 10];
z = [5 9 13 17 21];

Задайте инверсию и вперед сопоставляющие функции, которые принимают и возвращают точки в упакованном (x, y, z) формат.

inverseFcn = @(c)[c(:,1).^2,c(:,2).^2,c(:,3).^2];
forwardFcn = @(c)[sqrt(c(:,1)),sqrt(c(:,2)),sqrt(c(:,3))];

Создайте 3-D объект геометрического преобразования, tform, который хранит эти инверсия и передает функции отображения.

tform = geometricTransform3d(inverseFcn,forwardFcn)

tform = 
  geometricTransform3d with properties:

        InverseFcn: @(c)[c(:,1).^2,c(:,2).^2,c(:,3).^2]
        ForwardFcn: @(c)[sqrt(c(:,1)),sqrt(c(:,2)),sqrt(c(:,3))]
    Dimensionality: 3

Примените обратное преобразование этого 3-D геометрического преобразования к точкам ввода.

[u,v,w] = transformPointsInverse(tform,x,y,z)

u = 1×5

     9    25    49    81   121

v = 1×5

     4    16    36    64   100

w = 1×5

    25    81   169   289   441

Применяйте прямое геометрическое преобразование к преобразованным точкам u, v и w.

[x,y,z] = transformPointsForward(tform,u,v,w)

x = 1×5

     3     5     7     9    11

y = 1×5

     2     4     6     8    10

z = 1×5

     5     9    13    17    21

Задайте функцию обратного отображения, которая выполняет отражение о горизонтальной оси. Функция должна принять и возвратиться упакованный (x, y, z) координаты, где первые, вторые, и третьи столбцы содержат x-, y-, и z-координаты, соответственно.

inverseFcn = @(xyz)[xyz(:,1),-xyz(:,2),xyz(:,3)];

Создайте 3-D объект геометрического преобразования, tform, который хранит эту функцию обратного отображения.

tform = geometricTransform3d(inverseFcn)

tform = 
  geometricTransform3d with properties:

        InverseFcn: @(xyz)[xyz(:,1),-xyz(:,2),xyz(:,3)]
        ForwardFcn: []
    Dimensionality: 3

Загрузите и отобразите объем MRI, который будет преобразован.

s = load('mri');
mriVolume = squeeze(s.D);

Используйте imwarp, чтобы применять обратное геометрическое преобразование к объему входа MRI.

[mriVolumeTransformed] = imwarp(mriVolume,tform,'nearest','SmoothEdges',true);

Отобразите фрагменты изображения от объема входа MRI как монтаж.

montage(mriVolume,'Size',[4 8],'BackgroundColor','w')
title('Image Slices from 3-D MRI','FontSize',14)

Отобразите фрагменты изображения от преобразованного объема MRI как монтаж. Срезы преобразованного изображения являются отражением входных фрагментов изображения через ось X.

montage(mriVolumeTransformed,'Size',[4 8],'BackgroundColor','w')
title('Image Slices from Inverse Geometric Transformation of 3-D MRI','FontSize',14)