在matlab中裁剪图像

Question

我需要对matlab中的图像执行以下操作：

1. 加载图像
2. 计算图像的FFT（快速傅立叶变换）
3. 将频率分量移到中心
4. 像follow一样裁剪图像（如果图像分辨率为1000x1000，则图像所需的部分类似于以下坐标：100100800800。这是一个较小的图像。（应用过滤器去除高频的想法）
5. 反向移位
6. 傅里叶逆变换

. . .

我的代码如下所示：

I = imread('2.jpg'); %loading

ID = im2double(I);
FID = fft2(ID); %FFT
F = fftshift(FID); %shifting

F = imcrop(F,[100, 100, 800, 800]);

FID = ifftshift(F); %inverse of shifting
IFID = ifft2(FID); %inverse of FFT

I8 = im2uint8(IFID);

问题是，当我想裁剪图像时，我的功能无法裁剪类型为“复杂双”的图像，因为我认为...

错误：

使用im作物时出错

imcrop错误

输入错误（第93行）[x， y， a， cm，spatial_rect，h_image，placement_cancelled]=parseInput（varargin{：}）；

有帮助吗？...还有裁剪的另一个功能吗？

Answer 1

问题确实是张量F中的复数值。一个证明是将im作物（F，[100, 100, 800, 800]）替换为im作物（abs（F），[100, 100, 800, 800]）。

我的建议是将信号F分解为幅值和相位部分，并在这两幅图像上使用imcrop，然后进行重建。就理解“imcrop”而言，这应该会产生相同的结果。

代码将如下所示：

I = imread('PersonalPhoto.JPG'); %loading

ID = im2double(I);
FID = fft2(ID); %FFT
F = fftshift(FID); %shifting

Fabs = imcrop(abs(F),[100, 100, 800, 800]);
Fang = imcrop(angle(F),[100, 100, 800, 800]);

F = Fabs.*exp(1j.*Fang);

FID = ifftshift(F); %inverse of shifting
IFID = ifft2(FID); %inverse of FFT

I8 = im2uint8(IFID);

Answer 2

裁剪图像的实值分解（幅度和相位，或实部和虚部）将避免im作物报告的错误。

但是，由于您表示要过滤高频分量，因此应该注意，裁剪也会产生副作用，即减少生成的时域图像。

为了保持图像大小，您可以考虑将这些频率单元归零（不需要使用更精细的滤波器设计）。您可以通过将频域数据与一个矩阵相乘来实现这一点，该矩阵充当要保留的频率分量上的掩码。

实现这一点的代码如下所示：

% create a mask to zero-out high frequency components
H = zeros(size(I));
halfwidth = 400;
xmin = size(H,2)/2+1 - halfwidth;
xmax = size(H,2)/2+1 + halfwidth;
ymin = size(H,1)/2+1 - halfwidth;
ymax = size(H,1)/2+1 + halfwidth;
H(xmin:xmax, ymin:ymax, :) = 1;

% apply mask
F = F.*H;

这将为您提供以下整体代码：

I = imread('2.jpg'); %loading

ID = im2double(I);
FID = fft2(ID); %FFT
F = fftshift(FID); %shifting

% create a mask to zero-out high frequency components
H = zeros(size(I));
halfwidth = 400;
xmin = size(H,2)/2+1 - halfwidth;
xmax = size(H,2)/2+1 + halfwidth;
ymin = size(H,1)/2+1 - halfwidth;
ymax = size(H,1)/2+1 + halfwidth;
H(xmin:xmax, ymin:ymax, :) = 1;

% apply mask
F = F.*H;

FID = ifftshift(F); %inverse of shifting
IFID = real(ifft2(FID)); %inverse of FFT

I8 = im2uint8(IFID);