Matlab 亚像素提取代码示例

以下是一个简单的示例代码，用于在 Matlab 中进行亚像素提取：

% 读取图像
image = imread('image.jpg');
% 将图像转换为灰度图像
grayImage = rgb2gray(image);
% 计算梯度
[gradX, gradY] = gradient(double(grayImage));
% 计算亚像素位移
subPixelX = subpixel(gradX);
subPixelY = subpixel(gradY);
% 显示结果
imshow(image);
hold on;
quiver(subPixelX, subPixelY, 'r');
hold off;

function subPixel = subpixel(grad)
    % 计算梯度的一阶导数
    d1 = diff(grad, 1, 2);
    % 计算梯度的二阶导数
    d2 = diff(grad, 2, 2);
    % 计算亚像素位移
    subPixel = -d1 ./ d2;
end

这段代码首先读取一张图像，然后将其转换为灰度图像。接下来，通过计算灰度图像的梯度，得到图像的一阶导数。然后，通过计算一阶导数的差分和二阶导数的差分，得到亚像素位移。最后，使用quiver函数在原始图像上绘制亚像素位移的箭头。