图像的H和V方向分开来单独计算,mark一下,matlab仿真一个方向的,就是设计一个高通滤波器,把滤波出来的像素相加,
滤波器采用FDAtool设计,FIR或者IIR。
im1 = imread('1.jpg');
im2 = imread('2.jpg');
im3 = imread('3.jpg');
im4 = imread('4.jpg');
im1gray = im2double(rgb2gray(im1));
im2gray = im2double(rgb2gray(im2));
im3gray = im2double(rgb2gray(im3));
im4gray = im2double(rgb2gray(im4));
im1line = im1gray(1500,[1:4032]);
im2line = im2gray(1500,[1:4032]);
im3line = im3gray(1500,[1:4032]);
im4line = im4gray(1500,[1:4032]);
im1fftmag = abs(fftshift(fft(im1line,4032)));
im2fftmag = abs(fftshift(fft(im2line,4032)));
im3fftmag = abs(fftshift(fft(im3line,4032)));
im4fftmag = abs(fftshift(fft(im4line,4032)));
f = [1:4032]/4032;
figure;
%subplot(4,1,1);
%plot(f,im1fftmag);
subplot(5,1,1);
plot(f-0.5,im1fftmag,'b');
title('谱分析');
subplot(5,1,2);
plot