blob: 7e14c8b63eb3da48fc8d7c032db33b5c84f5922a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
|
//Caption:Scilab code to intergchange phase information between two images
//Example4.6
//page 174-175
clc;
close;
a = imread('E:\DIP_JAYARAMAN\Chapter4\lena.png'); //SIVP toolbox
b = imread('E:\DIP_JAYARAMAN\Chapter4\baboon.png');
a = rgb2gray(a);
b = rgb2gray(b);
a = imresize(a,0.5);
b = imresize(b,0.5);
figure(1)
ShowImage(a,'Original lena Image'); //IPD toolbox
title('Original lena Image');
figure(2)
ShowImage(b,'Original baboon Image');
title('Original baboon Image')
ffta = fft2d(double(a));
fftb = fft2d(double(b));
mag_a = abs(ffta);
mag_b = abs(fftb);
ph_a = atan(imag(ffta),real(ffta));
ph_b = atan(imag(fftb),real(fftb));
newfft_a = mag_a.*(exp(%i*ph_b));
newfft_b = mag_b.*(exp(%i*ph_a));
rec_a = ifft2d(newfft_a);
rec_b = ifft2d(newfft_b);
figure(3)
ShowImage(uint8(rec_a),'lena Image after phase reversal');
title('lena Image after phase reversal')
figure(4)
ShowImage(uint8(rec_b),'baboon Image after phase reversal');
title('baboon Image after phase reversal')
|
