topview_colormap_ssoctdata

tic;
clc, clear all, close all
addpath('D:\Labview Projects\matlab_script_skin');

%filePath = 'C:\Users\iptea\Documents\MATLAB\MEMS_SD';
filePath = 'C:\Users\iptea\Documents\MATLAB\MEMS_SD';

imgFilePath = dir([filePath '/*.jpg']);
imgFilesNo = size(imgFilePath,1);
imgName = imgFilePath(1).name;
srcimg = imread([filePath,'\',imgName]);

l_crop = 78;
r_crop = 112;
t_crop = 35;
b_crop = 47;

srcimg = rgb2gray(srcimg);
[H,W] = size(srcimg);
[srcimg,H,W] = crop(srcimg,H,W,l_crop,r_crop,t_crop,b_crop);

img3d = zeros(H,W,imgFilesNo);
imgWidth = W;

mrg = 60;

diff_idl = 2;
diff_rl = 2;
gsigma1 = 2;
gsigma2 = 2;
med = 5;

%  hbar = parfor_progressbar(imgFilesNo+imgWidth,'Please wait...'); %create the progress bar
% parfor_progress(imgFilesNo); % Initialize
% two_di_arr = zeros(imgFilesNo,385); %300

for step = 1:imgFilesNo
        %fprintf('.');
%         if mod(step,100) == 0
%             fprintf('%i\n',step);
%         end

        imgName = imgFilePath(step).name;
        img = imread([filePath,'\',imgName]);
        img = rgb2gray(img);
        [H,W] = size(img);

        [img,H,W] = crop(img,H,W,l_crop,r_crop,t_crop,b_crop);

        img = uint8(img);

        img3d(:,:,step) = img;

        img_gau = imgaussfilt(img,gsigma1);
        %imshow(img_gau)
        img_med = medfilt2(img_gau,[med*2 med*2],'symmetric');

        img_med = filt_diff(img_med,H,W,diff_idl);
        %imshow(img_med,[])
        %[~,ind_idl,~] = find_srf(img_med,H,W); %img_idl
        [~,ind_idl,~] = find_srf(img_med,H,W); %img_idl
        %imshow(img_idl)
        [img_flt,H,W] = flatten(img,H,W,ind_idl,mrg);

        img_gau = imgaussfilt(img_flt,gsigma2);

        img_diff = filt_diff(img_gau,H,W,diff_rl);
        %imshow(img_diff,[])
        [~,ind_rl,avg_rl] = find_srf(img_diff,H,W); %img_rl
        %imshow(img_rl,[])
        %[avg_rghn_fr(step),rms_rghn_fr(step),max_rghn_fr(step)] = rgh(avg_rl,ind_rl);
        two_di_arr(step,:) = ind_rl;
        waitbar(step / imgFilesNo)
%         pause(rand); % computation
%         hbar.iterate(1);
end

%two_di_arr_norma = two_di_arr-min(two_di_arr(:));
%two_di_arr_norma = two_di_arr/max(two_di_arr(:));

med_two_di_arr = medfilt2(two_di_arr,[2 2]);
%med_two_di_arr_norma = medfilt2(two_di_arr_norma,[7 7]);


figure()
image(med_two_di_arr,'CDataMapping' , 'scaled')
caxis([59 63])
% xlabel('x axis', 'FontSize', 24);
% ylabel('y axis', 'FontSize', 24);
% set(gca,'DefaultTextFontSize',24)

colormap jet
colorbar

wholetime = toc;
% close(hbar);
%%
max(two_di_arr(:))
min(two_di_arr(:))
%%