三维图像旋转（Matlab）

Question

我想在Matlab xyz三维空间中旋转2D图像。我想旋转图像我围绕x轴角度θ= i。我用旋转矩阵（定义here）乘法的方式这样做：

这是我的代码：

x = linspace(-1,1,size(I,1)); 
[x0,y0,z0] = meshgrid(x,x,1); 
R = [1 0 0; 0 cosd(-i) -sind(-i); 0 sind(-i) cosd(-i)]; 
xy = [x0(:),y0(:), z0(:)]*R'; 
X = reshape(xy(:,1),size(x0)); 
Y = reshape(xy(:,2),size(y0)); 
Z = reshape(xy(:,3),size(z0)); 
rotatedI = interp3(x0,y0,z0,I,X,Y,Z, 'nearest'); 
rotatedI(isnan(rotatedI)) = 0; 

，我得到的错误是在插线：

Error using griddedInterpolant 
Interpolation requires at least two sample points in each 
dimension. 

那么究竟是什么问题，我该如何解决？

我很困惑，因为这个代码工作很好，当专门到两个维度。任何解释都非常感谢。

Answer 1

问题是您的输入z网格z0包含单个值（1），而您期望的z网格Z包含多个值。因此无法将此单一值插入多个值中。

您可以为输入z网格指定多个值，以便可以进行插值。

I = im2double(imread('cameraman.tif')); 
I = cat(3,I,I,I); 
i = 10; 
x = linspace(-1,1,size(I,1)); 
[x0,y0,z0] = meshgrid(x,x,-1:1); % 3 z values 
R = [1 0 0; 0 cosd(-i) -sind(-i); 0 sind(-i) cosd(-i)]; 
xy = [x0(:),y0(:), z0(:)]*R; 
X = reshape(xy(:,1),size(x0)); 
Y = reshape(xy(:,2),size(y0)); 
Z = reshape(xy(:,3),size(z0)); 
rotatedI = interp3(x0,y0,z0,I,X,Y,Z, 'nearest'); 
rotatedI(isnan(rotatedI)) = 0; 
imshow(rotatedI); 

此外，你可以使用rotx，roty，并rotz，容易产生旋转矩阵。

UPDATE

，如果你想在3D空间中旋转图像二维平面上，你并不需要使用插值，可以这样做：

I = im2double(imread('cameraman.tif')); 
i = 10; 
x = linspace(-1,1,size(I,1)); 
[x0,y0,z0] = meshgrid(x,x,-1:1); % 3 z values 
[x0,y0,z0] = meshgrid(x,x,1); % 3 z values 
R = [1 0 0; 0 cosd(-i) -sind(-i); 0 sind(-i) cosd(-i)]; 
xy = [x0(:),y0(:), z0(:)]*R; 
X = reshape(xy(:,1),size(x0)); 
Y = reshape(xy(:,2),size(y0)); 
Z = reshape(xy(:,3),size(z0)); 
warp(X,Y,Z,I); 
grid on