clear 
clc
close all
warning('OFF');

tic

% load mat files
[file_name,path_name] = uigetfile('*','Select file for Kors transformation');
path_name_saving=path_name;

%% file import
file_ID          = strsplit(file_name,'.');
file_path        = fullfile(path_name,file_name);


matfile          = matfile(file_path)


%% ===================== load variables from .mat file ========================
ECG12Lead      = matfile.ECG12Lead;

XYZ_O=kors(ECG12Lead);

save(file_path,'-append','XYZ_O');


%% =============================== 12 Lead Plot ====================================

    ECG12L = figure('visible','on');
    ax(1)=subplot(3,4,1);plot(ECG12Lead(:,1));ylabel('Lead I');
    ax(2)=subplot(3,4,2);plot(ECG12Lead(:,2));ylabel('Lead II');
    ax(3)=subplot(3,4,3);plot(ECG12Lead(:,3));ylabel('Lead III');
    ax(4)=subplot(3,4,4);plot(ECG12Lead(:,4));ylabel('Lead aVR');
    ax(5)=subplot(3,4,5);plot(ECG12Lead(:,5));ylabel('Lead aVL');
    ax(6)=subplot(3,4,6);plot(ECG12Lead(:,6));ylabel('Lead aVF');
    ax(7)=subplot(3,4,7);plot(ECG12Lead(:,7));ylabel('Lead V1');
    ax(8)=subplot(3,4,8);plot(ECG12Lead(:,8));ylabel('Lead V2');
    ax(9)=subplot(3,4,9);plot(ECG12Lead(:,9));ylabel('Lead V3');
    ax(10)=subplot(3,4,10);plot(ECG12Lead(:,10));ylabel('Lead V4');
    ax(11)=subplot(3,4,11);plot(ECG12Lead(:,11));ylabel('Lead V5');
    ax(12)=subplot(3,4,12);plot(ECG12Lead(:,12));ylabel('Lead V6');
    linkaxes(ax,'x');
%    saveas(ECG12L,strcat(images_folder,name0,'_12Lead'),'fig');


%% =============================== XYZ Leads Plot ====================================
    
    ECG3L = figure('visible','on');
    subplot(3,1,1);plot(XYZ_O(:,1));ylabel('Lead X');
    subplot(3,1,2);plot(XYZ_O(:,2));ylabel('Lead Y');
    subplot(3,1,3);plot(XYZ_O(:,3));ylabel('Lead Z');
%    saveas(ECG3L,strcat(images_folder,name0,'_3Lead'),'fig');


%% ============================ Function - Kors Matrix ============================

function transform_k = kors(leads12) %data in columns 12 leads
korsMatrix = [0.38, -0.07, 0, 0, 0, 0, -0.13, 0.05, -0.01, 0.14, 0.06,0.54;
              -0.07, 0.93, 0, 0, 0, 0, 0.06, -0.02, -0.05, 0.06, -0.17, 0.13;
              0.11, -0.23, 0, 0, 0, 0, -0.43, -0.06,-0.14,-0.20,-0.11,0.31];
transform_k = leads12 * korsMatrix';
end