自动对焦求傅里叶谱的高精度matlab代码

以下是一个自动对焦求傅里叶谱的高精度Matlab代码示例：

% 将图像读入Matlab
image = imread('image.jpg');

% 将图像转换为灰度图像
gray_image = rgb2gray(image);

% 定义图像的大小
[m, n] = size(gray_image);

% 定义高斯滤波器的参数
sigma = 10;
hsize = [5, 5];

% 创建高斯滤波器
h = fspecial('gaussian', hsize, sigma);

% 对图像进行高斯滤波
blurred_image = imfilter(gray_image, h);

% 定义搜索范围
search_range = 20;

% 定义搜索步长
step = 1;

% 初始化最大相关系数和最佳偏移量
max_corr = -Inf;
best_offset = [0, 0];

% 遍历搜索范围
for i = -search_range:step:search_range
    for j = -search_range:step:search_range
        % 对模板进行平移
        shifted_template = circshift(blurred_image, [i, j]);
        
        % 计算模板和平移后的图像的互相关系数
        corr = sum(sum(blurred_image .* shifted_template));
        
        % 更新最大相关系数和最佳偏移量
        if corr > max_corr
            max_corr = corr;
            best_offset = [i, j];
        end
    end
end

% 对原始图像进行平移以获得最佳对焦
aligned_image = circshift(image, best_offset);

% 显示原始图像和对焦后的图像
figure;
subplot(1, 2, 1);
imshow(image);
title('Original Image');
subplot(1, 2, 2);
imshow(aligned_image);
title('Aligned Image');

这个代码使用了自动对焦技术以获得最佳对焦。它首先将图像转换为灰度图像，并对其进行高斯滤波以减少噪声。然后，它在给定的搜索范围内搜索最佳偏移量，以最大化模板和平移后的图像之间的互相关系数。最后，它对原始图像进行平移以获得最佳对焦。