% itype = 2 for Pulsed Wave Doppler data

itype                       = 2;

% Velocity Distribution [Units = cm/s]

Velocity_axis               = linspace(0, dataout{itype}.meta{1}.Vrange(2), size(dataout{itype}.img{1}, 1))';

% Spectogram [Rows represents the Velocity and columns are the time stamp]

Pulsed_Doppler_Data         = 10*log10(dataout{itype}.img{1}(:,:) + 1);

% Total number of time stamps

Num_Time_Stamps             = size(Pulsed_Doppler_Data, 2);

%%Plot the spectogram

%imagesc((1:Num_Time_Stamps)', Velocity_axis, Pulsed_Doppler_Data); set(gca,'YDir','normal');

% Find the maximum of the Pulsed Doppler data in each column .

% Here, Column represents the time stamps.

B = find(Pulsed_Doppler_Data == max(Pulsed_Doppler_Data), size(Pulsed_Doppler_Data, 2), 'first');

% Repeat the velocity vector for each column

new_Velocity_axis           = repmat(Velocity_axis, [1, size(Pulsed_Doppler_Data, 2)]);

% Calculate the velcotiy data [Units = m/s]

handles.Velocity_Data       = new_Velocity_axis(B)/100;

bframe                      = Pulsed_Doppler_Data;

% This saves the single frame 

%bframe                      = permute(dataout{1}.img{1}, [2, 1, 3]);