1

如何在转换（旋转、平移、缩放和剪切）图像后自动裁剪黑色边框？

我正在使用基于强度的自动图像配准（imregtform、优化器、imregister ....）来旋转一个图像以适应另一个图像。现在旋转后，我需要裁剪两个图像以具有相同的大小，并且仍然排成一行。

旋转图像的最简单表示形式是应用了变换矩阵的白色正方形。

可以在这里看到我正在使用的一组实际图像的示例。

我假设图像边缘可能具有零值像素的事实可能会使操作复杂化，尽管我假设由于噪声，整个图像将是非零的。

也许我可以使用变换矩阵来实际计算需要裁剪的边界。在上面的例子中，矩阵是：

 0,999428374496743      0,00888472048904662   0
-0,00888472048904659    0,999428374496743     0
 3,79626401832983      -0,493066986575474     1

该矩阵列在工作区中的“1x1 仿射 2d”对象中。我无法从那里找到使用它的语法。

4

1 回答 1

1

我已经想出了一种方法，虽然它不是很优雅，因此我仍然对解决这个问题的好方法感兴趣。

可以将变换矩阵应用于图像的四个角点，并将这些新点用作裁剪限制。应确保裁剪角在原始图像内或需要固定到边缘。

% width and length of the input image
width = size(img_fixed,1);
height = size(img_fixed,2);

% transform the four corners of the image to find crop area
[x1,y1] = transformPointsForward(T,0,0);
[x2,y2] = transformPointsForward(T,width,0);
[x3,y3] = transformPointsForward(T,width,height);
[x4,y4] = transformPointsForward(T,0,height);

% find inner most borders for a rectangular crop
if max([x1,x4]) < 0
    x_left = 0;      
else
    x_left = ceil(max([x1,x4]));
end

if min([x2,x3]) > width
    x_right = width;     
else
    x_right = floor(min([x2,x3]));
end

if max([y1,y2]) < 0
    y_top = 0;   
else
    y_top = ceil(max([y1,y2])); 
end

if min([y3,y4]) > height
    y_bottom = height;      
else
    y_bottom = floor(min([y3,y4]));
end

img_fixed_crop = imcrop(img_fixed,[x_left y_top x_right-x_left y_bottom-y_top]);
img_moving_crop = imcrop(img_moving,[x_left y_top x_right-x_left y_bottom-y_top]);
于 2015-09-16T11:18:34.293 回答