diff --git a/Seance1.m b/Seance1.m new file mode 100644 index 0000000000000000000000000000000000000000..dac05180b9a32a5fd8412681c753ab8a9ed09d7b --- /dev/null +++ b/Seance1.m @@ -0,0 +1,126 @@ +%% +clear all; +I = imread('Amosite5.jpg'); +figure(1) +imagesc(I) +colormap(gray) +colorbar +img_fft = fftshift(fft2(I)); +figure(2) +imagesc(20*log10(abs(img_fft))) +%% +clear all; +load('angiocor.mat'); +figure(1) +imagesc(angiocor) +colormap(gray) +colorbar +img_fft = fftshift(fft2(angiocor)); +figure(2) +imagesc(20*log10(abs(img_fft))) +%% +clear all; + +load('barbara.mat'); +figure(1) +imagesc(X) +colormap(gray) +colorbar +figure(2) +img_fft = fftshift(fft2(X)); +imagesc(20*log10(abs(img_fft))) + +%% +clear all; + +load('cuir.mat'); +figure(1) +imagesc(X) +colormap(gray) +colorbar +figure(2) +img_fft = fftshift(fft2(X)); +imagesc(20*log10(abs(img_fft))) + +%% +clear all; + +load('empreinte_roi.mat'); +figure(1) +imagesc(image) +colormap(gray) +colorbar +img_fft = fftshift(fft2(image)); +figure(2) +imagesc(20*log10(abs(img_fft))) + +Z = 20*log10(abs(fft2(image))); +[M,I] = max(Z(:)); +[I_row, I_col] = ind2sub(size(Z),I); +N = size(Z); +vu1 = (I_col - 1)/N(1); +vu2 = (I_row - 1)/N(1); +rho = sqrt(vu1*vu1 + vu2*vu2); +phi = atan(vu2/vu1); + +% +x = [-N/2:1:N/2]; +y = transpose([-N/2:1:N/2]); + +H = exp(-x*x+y*y) + +%% +clear all; +load('head.mat'); +figure(1) +imagesc(X) +colormap(gray) +colorbar + +img_fft = fftshift(fft2(X)); +figure(2) +imagesc(20*log10(abs(img_fft))) + +%% +clear all; +load('lena512.mat'); +figure(1) +imagesc(lena512) +colormap(gray) +colorbar +img_fft = fftshift(fft2(lena512)); +figure(2) +imagesc(20*log10(abs(img_fft))) + +%% +clear all; +load('poupee.mat'); +figure(1) +imagesc(X) +colormap(gray) +colorbar +figure(2) +img_fft = fftshift(fft2(X)); +imagesc(20*log10(abs(img_fft))) + +%% +clear all; +load('tartan.mat'); +figure(1) +imagesc(X) +colormap(gray) +colorbar +figure(2) +img_fft = fftshift(fft2(X)); +imagesc(20*log10(abs(img_fft))) + +%% +clear all; +load('tyre.mat'); +figure(1) +imagesc(X) +colormap(gray) +colorbar +figure(2) +img_fft = fftshift(fft2(X)); +imagesc(20*log10(abs(img_fft))) \ No newline at end of file