% These commands are optimized for DICOM images saved from Philips EPIQ 7C
% We extract the B-mode images, ECG, Time Stamp;
% Resolution of B model images in the X and Y direction are also calculated 

% Local name the Data file
Data_file                   = 'input';

% Exact filename
file_name                   = varargin{1};

% if file does not exist, return
if ~isfile(file_name)
    fprintf('File does NOT exist in this location: %s\n', file_name)
    handles.file_exist = false;
    return; 
else
    handles.file_exist = true;
end

% copy the data from source to the this folder and saveas Data_file
copyfile(file_name, Data_file);

fprintf('Loading the DICOM IMAGES. It might take time!\n');

% Load the data into MATLAB via helper files from Philips
% This function was written by Jonathan (Philips)
[~, dataout, ~]             = DICOMReaderGeneral(Data_file);

if (size(dataout{1}.time{1}, 2) ~= 1)
    
    % This variable is set to true if the input os bmode.
    handles.Bmode           = true;
    
    % Title of the reults plots
    handles.Waveform_Title.String = "Diameter [mm] Vs Time [seconds]";
    
    % Update handles about the current position
    handles.Message_Bar.String = "Loaded BMode images";
    
    % Extract b mode images 
    script_extract_B_mode_images;
    
    handles.Play_Button.String = 'Process B-Mode Data';
    
else
    
    % This variable is set to true if the input os bmode.
    handles.Bmode           = false;
    
    % Title of the reults plots
    handles.Waveform_Title.String = "Velocity [m/s] Vs Time [seconds]";
    
    % Extract pulsed Doppler data
    script_extract_Pulse_Doppler_data;
    
    % Update handles about the current position
    handles.Message_Bar.String = "Loaded Pulsed Doppler Data";
    
    handles.Play_Button.String = 'Process Pulsed Doppler Data';
   
end

% Location of the path to save results
handles.save_path           = script_create_save_path_name(file_name, handles.Bmode);

% Load ECG from the file. ECG is normalized using maximum value
physio                   	= DICOMPhysioReader(Data_file);

if ~(size(physio.time, 1) == 0)
    time_ecg            	= (physio.time-physio.time(1))/1e6;
    ECG                 	= double(physio.channels(:,1));
    % Subsample ECG to # frames of B mode image
    time_frame            	= linspace(0, max(time_ecg), Num_Time_Stamps)';
    interpolate_ECG       	= interp1(time_ecg, ECG, time_frame); 
    handles.time         	= time_frame;
    handles.ECG           	= interpolate_ECG/max(interpolate_ECG);
else
    fprintf('NO ECG Data\n');
    time_ecg              	= [];
    handles.ECG          	= [];
    Fs                      = 42;
    time_frame             	= (0:1/Fs:(Num_Time_Stamps-1)/Fs)'; 
    handles.time         	= time_frame;
end