Riesz变换MATLAB代码

function [orientation, phase, amplitude] = riesz(image)

[imHeight, imWidth] = size(image);

%approx riesz, from Riesz Pyramids for Fast Phase-Based VideoMagnification

dxkernel = zeros(size(image));

dxkernel(1, 2)=-0.5;

dxkernel(1,imWidth) = 0.5;

dykernel = zeros(size(image));

dykernel(2, 1) = -0.5;

dykernel(imHeight, 1) = 0.5;

R1 = ifft2(fft2(image) .* fft2(dxkernel));

R2 = ifft2(fft2(image) .* fft2(dykernel));

orientation = zeros(im