38 lines
1008 B
Matlab
38 lines
1008 B
Matlab
%%%%%%%%%%%%%%%%%%%%%
|
|
% Task: Creating FFT of average R,G,B over a period of 30 seconds
|
|
%
|
|
% Input:
|
|
% - CSV file containing the RGB data for each pixels of the region of interest
|
|
%
|
|
% Output:
|
|
% - Fast Fourier Transform of the colors over time
|
|
%
|
|
% author: Delattre Loïc, Dey Maryne (loic.delattre@ecam.fr, maryne.dey@ecam.fr)
|
|
% date: 06/02/2023
|
|
%%%%%%%%%%%%%%%%%%%%%
|
|
|
|
pkg load io %% to be able to extract from external format (excel)
|
|
|
|
data = csvread('RGB_data_monfront.xlsx')
|
|
standard_deviation = std(data)
|
|
mean_value = mean(data)
|
|
|
|
for i = 1:size(data,1)
|
|
normalized_data_G(i,1) = (data(i,2)-mean_value(2))/standard_deviation(2); %%2 and not 1 because green is the 2nd color
|
|
endfor
|
|
|
|
fft_green = fft(normalized_data_G)
|
|
|
|
plot([1:size(fft_green,1)],fft_green)
|
|
|
|
|
|
|
|
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|
|
%%For all the colors
|
|
%%for j = 1:3
|
|
%% for i = 1:size(data,1)
|
|
%% normalized_data(i,j) = (data(i,j)-mean_value(j))/standard_deviation(j);
|
|
%% endfor
|
|
%%endfor
|
|
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% |