Я пишу алгоритм поиска фазы в Matlab, у меня есть несколько разных изображений, которые представляют разные миллилитры воды на поверхности листа (влажность поверхности листа), которые необходимо было вставить в алгоритм, чтобы увидеть фазовое изменение/фазовый сдвиг, который произошел в каждом миллиметр воды на поверхности листа, но я не совсем уверен, как увидеть фазовое изменение/фазовый сдвиг каждого изображения и сравнить его. Представляю ли я фазовое изменение/фазовый сдвиг на графике или есть другие способы увидеть фазовое изменение/фазовый сдвиг? Я рисую фазу или получаю значение? Могу ли я получить помощь в этом? Спасибо.
РЕДАКТИРОВАТЬ: тема, которой я сейчас занимаюсь, - это интерферометрия, и она создает полосы на каждом из изображений. Каждая из этих полос будет содержать значения/информацию о фазе, мой код должен получить информацию/значения о фазе каждого изображения, а затем мне нужно увидеть изменение фазы/сдвиг фазы, которое произошло на каждом изображении, но я Я не уверен, как я могу это сделать, и я не могу этого сделать. На основе алгоритма мне нужно доказать, что при увеличении/уменьшении уровня влажности значение фазы будет увеличиваться/уменьшаться/значение фазы будет смещаться в какую-то сторону. Теперь я не знаю, должен ли я видеть изменение фазы/сдвиг фазы каждого изображения с помощью графика или других методов. Я прикрепил изображение, чтобы показать пример входных изображений. /a> Я надеюсь, что смогу получить помощь в этом. Спасибо.
РЕДАКТИРОВАТЬ 2: изображение образца графика
Вот мой код:
close all; clc; workspace;
Iterations = 1100; %Number of iterations
p = 0.1; %Time to pause between display of images
tic; %Timer Initialization
error = [];
%This whole part is to generate a mesh plot, i still dont know what info
%can get from here but just leave it first.
x = linspace(-10,10,256);
y = linspace(-10,10,256);
[X,Y] = meshgrid(x,y);
x0 = 0; % Center
y0 = 0; % Center
sigma = 2; % Beam Waistfringes
A = 1; % Peak of the Beam
%These 2 lines, res and Source, is basically input intensity. U need
%this in order to get the exponential image to process ltr.
res = ((X-x0).^2 + (Y-y0).^2)./(2*sigma^2);
Source = A .* exp(-res);
%These 2 lines is for mesh plot also, ignore first
surf(Source);
shading interp
%Reading Fringe Image (Target)
Target = imread('0ml(1).jpeg');
Target = rgb2gray(Target);
Target = im2double(Target);
%Get the size of Target
[rowsimgA, colsimgA, numberOfColorChannelsimgA] = size(Source);
[rowsimgB, colsimgB, numberOfColorChannelsimgB] = size(Target);
%imresize function is basically to resize any image to a particular size.
%In this code u can see that i resize the Source beam input to the same
%size as the target image in the format of (Target, [rows, columns])
if rowsimgB ~= rowsimgA || colsimgA ~= colsimgB
Target = imresize(Target, [rowsimgA colsimgA]);
end
%FFT process
A = fftshift(ifft2(fftshift(Target)));
%Initiate loop
% The loop im not too sure yet, but basically when u multiplay the source
% and the phase, u produce exponential image, and then u do fft and ifft in
% order to reconstruct the image.
for i = 1:Iterations
B = abs(Source) .* exp(1i*angle(A));
C = fftshift(fft2(fftshift(B)));
D = abs(Target) .* exp(1i*angle(C));
A = fftshift(ifft2(fftshift(D)));
error = [error; (1/sqrt(rowsimgA*colsimgA)*sqrt(sum(sum(abs(C)-abs(Target)).^2)))];
end
%Display Outputs
figure(2), imagesc(Target), colorbar, title('Original Fringe')
figure(3), i = 1:1:i; plot(i,(error')); title('Error');
%Phase Mask
figure(4), imagesc(angle(A)), title('Phase Mask');
%Last Pattern
figure(5), imagesc(abs(C)), colorbar, title('Reconstructed Image');
%Reconstructed Image (Phase Only)
figure(6), imagesc(angle(D)), colorbar, title('Reconstructed Image (Phase Only)');
%Reconstructed Image (Amplitude Only)
figure(7), imagesc(abs(D)), colorbar, title('Reconstructed Image (Amplitude Only)');
%Display Total Time Taken for Phase Retrieval
toc;
`
