front_main_multiple_pmrch[-]

close all; clear all; clc;
addpath('C:\Users\iptea\Documents\MATLAB\Nonreplica September_ROUGHNESS\Roughness_Amore_2\General');

sbj = {'01_KB','02_YJ','03_PU','04_HG','05_HR',...
    '06_HG','07_EJ','08_MA','09_NY','10_MR'};
part = {'forearm','AF_forearm','forehead','AF_forehead','backhand','AF_backhand'};

%% Attention:
% diff_idl should ne greater than diff_rl
% size of median filter neighborhood should be greater than desired radius

%%
data_pth = 'C:\Users\iptea\Documents\MATLAB\AMORE_ROUGHNESS';

sht_name = 'rgh_frnt.xlsx';
ttl = {'part','sbj','avg_rgh','rms_rgh','max_rgh'};
ttl_rng = 'B1:F1';
xlswrite(sht_name,ttl,ttl_rng);

max_no = 1;
avg_rgh = zeros(1,max_no);
rms_rgh = zeros(1,max_no);
max_rgh = zeros(1,max_no);
counter = 0;

fprintf('Roughness calculation started:\n');

for part_no = 1:length(part)
    for sbj_no = 1:length(sbj)
        %fprintf('%i/%i\n', part_no*sbj_no);
        fld_pth = strcat(data_pth,'\',sbj{sbj_no},'\',part{part_no},'\');
        for img_no = 1:1 %max_no
            img_name = strcat('Image',num2str(img_no-1,'%04d'),'.jpg');
            img = imread(strcat(fld_pth,img_name));
            % figure()
            % imshow(img);
            [H,W] = size(img);

            %% Region of interest
            l_crop = 130;
            t_crop = 20;
            r_crop = 0;
            b_crop = 300;

            [img_crp,H,W] = crop(img,H,W,l_crop,r_crop,t_crop,b_crop);
            % figure()
            % imshow(img_crop);
            %% Prefiltering
            gsigma = 6;
            img_gau = imgaussfilt(img_crp,gsigma);
            %% Finding ideal surface
            med = 60;
            diff_idl = 15;
            img_med = medfilt2(img_gau,[med*2 med*2],'symmetric');
            img_med = filt_diff(img_med,H,W,diff_idl);
            [~,ind_idl,~] = find_srf(img_med,H,W); %img_idl&avg_idl not used
            % figure()
            % imshow(img_idl);

            %% Flattening surface according to ideal line
            % Note: We put gaussian image directly
            mrg = 50;
            [img_flt,H,W] = flatten(img_crp,H,W,ind_idl,mrg);
            % figure()
            % imshow(img_flt,[]);

            %% Prefiltering flattened image
            gsigma = 6;
            img_gau = imgaussfilt(img_flt,gsigma);

            %% Finding real surface
            diff_rl = 5;
            img_diff = filt_diff(img_gau,H,W,diff_rl);
            [~,ind_rl,avg_rl] = find_srf(img_diff,H,W); %img_rl
            % figure()
            % imshow(img_rl,[])
            % hold on
            % plot(avg_rl)

            %% Calculating roughness
            [avg_rgh(img_no),rms_rgh(img_no),max_rgh(img_no)] = rgh(avg_rl,ind_rl);
        end
        mean_avg_rgh = mean(avg_rgh);
        mean_rms_rgh = mean(rms_rgh);
        mean_max_rgh = mean(max_rgh);

        A = [part(1,part_no),sbj(1,sbj_no),num2str(mean_avg_rgh),num2str(mean_rms_rgh),...
        num2str(mean_max_rgh)];

        %xlRangeVariables = strcat('B',num2str((nm2L-1)*10+nm1L+1),':','G',num2str((nm2L-1)*10+nm1L+1));
        if mod(part_no,2) == 1
            rw_bef = part_no - 1;
            vrbl_rng = strcat('B',num2str(rw_bef*10+sbj_no+1),':','F',num2str(rw_bef*10+sbj_no+1));
        else
            rw_aft = part_no - 2;
            vrbl_rng = strcat('H',num2str(rw_aft*10+sbj_no+1),':','L',num2str(rw_aft*10+sbj_no+1));
        end
        xlswrite(sht_name,A,vrbl_rng);
        counter = counter + 1;
        fprintf('%i out of %i\n', counter, length(part)*length(sbj));
    end
end

fprintf('Process finished.\n');
%%
% sht_name = 'rgh_frnt.xlsx';
% T = readtable(sht_name,'ReadRowNames',false);
% A = T{1:10,4}
% B = T{1:10,9}
% C = A > B
% xlswrite(sht_name,C,vrbl_rng);
% sht_name = 'rgh_frnt - Copy.xlsx';
%cmpr_rng = strcat('N',num2str(rw_aft*10+sbj_no+1));


T = readtable(sht_name,'ReadRowNames',false);
A = T{1:10,4}
B = T{1:10,9}
C = A > B
xlswrite(sht_name,C,'N2:N11');
A = T{21:30,4}
B = T{21:30,9}
C = A > B
xlswrite(sht_name,C,'N22:N31');
A = T{41:50,4}
B = T{41:50,9}
C = A > B
xlswrite(sht_name,C,'N42:N51');

% bef = T{rw_aft*10+sbj_no+1:rw_aft*10+sbj_no+1,4};
% aft = T{rw_aft*10+sbj_no+1:rw_aft*10+sbj_no+1,9};
% C = bef > aft;

%%
% for part_no = 1:length(part)
%     for sbj_no = 1:length(sbj)
%
%     end
% end
%%