fprintf('***************************************\n');
fprintf('2D Phase Unwrapping Demo\n');
fprintf('Please select the demo:\n');
fprintf('(1) No noise , no ignored region\n');
fprintf(' 2. With noise, no ignored region\n');
fprintf(' 3. No noise , with ignored region\n');
fprintf(' 4. With noise, with ignored region\n');
while (1)
user_input = input('Your selection (1-4): ', 's');
user_input = strip(user_input);
% if the user does not supply anything, select the default
if strcmp(user_input, '')
fprintf('Demo 1 is selected\n');
user_input = '1';
end
if length(user_input) == 1 && sum(user_input == '1234') == 1
break;
else
fprintf('Invalid input\n');
end
end
[X, Y] = meshgrid(linspace(-1, 1, 512) * 5);
img = -(X.*X + Y.*Y);
fprintf('Image size: %dx%d pixels\n', size(img,1), size(img,2));
% add noise
if any(user_input == '24')
img = img + randn(size(X)) * 0.5;
end
% add an ignored region
if any(user_input == '34')
img(end/4:3*end/4,end/4:3*end/4) = nan;
end
% wrap the image
wimg = wrapTo2Pi(img);
tic;
unwrap_img = unwrap_phase(wimg);
toc;
subplot(221);
pcolor(img);
shading flat;
set(gca, 'ydir', 'reverse');
title('Original phase');
subplot(222);
pcolor(wimg);
shading flat;
set(gca, 'ydir', 'reverse');
title('Wrapped phase');
subplot(223);
pcolor(unwrap_img);
shading flat;
set(gca, 'ydir', 'reverse');
title('Unwrapped phase');
subplot(224);
pcolor(wrapTo2Pi(unwrap_img));
shading flat;
set(gca, 'ydir', 'reverse');
title('Rewrap of unwrapped phase');
- 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.
- 34.
- 35.
- 36.
- 37.
- 38.
- 39.
- 40.
- 41.
- 42.
- 43.
- 44.
- 45.
- 46.
- 47.
- 48.
- 49.
- 50.
- 51.
- 52.
- 53.
- 54.
- 55.
- 56.
- 57.
- 58.
- 59.
- 60.
- 61.
- 62.
- 63.
- 64.
- 65.
- 66.
- 67.
- 68.
- 69.
- 70.
- 71.
- 72.
- 73.
- 74.
- 75.
- 76.
- 77.
- 78.
- 79.
- 80.
- 81.
- 82.
- 83.
- 84.
- 85.
- 86.
- 87.
- 88.
- 89.
- 90.
- 91.
- 92.
- 93.
- 94.
- 95.
- 96.
- 97.
- 98.
- 99.
- 100.
- 101.
- 102.
- 103.
- 104.
- 105.
- 106.
- 107.
- 108.
- 109.
- 110.
- 111.
- 112.
- 113.
- 114.
- 115.
- 116.
- 117.
- 118.
- 119.
- 120.
- 121.
- 122.
- 123.
- 124.
- 125.
- 126.
- 127.
- 128.
- 129.