Skip to content

discrete wavelet display - uses patches instead of rectangle to display the discrete wavelet #43

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
Edouard2laire wants to merge 9 commits into
base: master
Choose a base branch
from
Open

discrete wavelet display - uses patches instead of rectangle to display the discrete wavelet #43

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
c8514e2
uses patches instead of rectangle to display the discrete wavelet
Edouard2laire Sep 5, 2025
5b0db60
Organize the display of boxes, by scales
Edouard2laire Sep 5, 2025
4c7fb9d
Small reorganization - still display the same thing
Edouard2laire Sep 5, 2025
ff43b82
WIL - use one patch per scale
Edouard2laire Sep 5, 2025
4321c78
fix off-by-one error
Edouard2laire Sep 5, 2025
39e5969
Continue the optimization
Edouard2laire Oct 16, 2025
3e42c8a
keep the spacing between scale
Edouard2laire Oct 16, 2025
e8a6252
keep only one window active at a time
Edouard2laire Oct 16, 2025
ac4b857
change marker for used boxes
Edouard2laire Oct 16, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
166 changes: 96 additions & 70 deletions best/wmem/display/be_display_time_scale_boxes.m
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
function [hp,hptab] = be_display_time_scale_boxes(obj,OPTIONS)
%DISPLAY_ENTROPY_DROPS displays a graphical window with MEM entropy drops at
% each iteration
function [hp, hptab] = be_display_time_scale_boxes(obj, OPTIONS)
%be_display_time_scale_boxes displays the discrete wavelet
%representation
%
% INPUTS:
% - obj
Expand Down Expand Up @@ -39,92 +39,118 @@
% N : nomber of samples
% TFboxes is a structure .k ans .j are (j,k) coordinates

hp = obj.hfig;
hptab = obj.hfigtab;

if ~strcmp(obj.data_type,'discrete_wavelet')
disp('NO WAVELET DISPLAY');
return
end

for ii = 1:length(OPTIONS.mandatory.DataTypes)

hhh = uitab(hptab,'title', OPTIONS.mandatory.DataTypes{ii});
hpc = uipanel('Parent', hhh, ...
'Units', 'Normalized', ...
'Position', [0.01 0.01 0.98 0.98], ...
'FontWeight','demi');
set(hpc,'Title', 'Time-frequency boxes', 'FontSize', 8);

Tmin = obj.t0 - (obj.info_extension.start-1)/OPTIONS.automatic.sampling_rate;
Tmax = Tmin + (size(obj.data{ii},2)-1)/OPTIONS.automatic.sampling_rate;

scales = unique(OPTIONS.automatic.Modality(ii).selected_jk(2,:));
J = max(scales);

ax = axes( 'parent', hpc, ...
'outerPosition', [0.01 0.01 0.98 0.98], ...
'YTick', 0.5:J-0.5, ...
'YTickLabel', num2cell(1:J), ...
'Xlim', [Tmin, Tmax], ...
'Ylim', [0, J], ...
'Box','on');

sBox = create_rectangles(obj, OPTIONS, ii);
patch(ax, sBox);
xlabel(ax, 'time (s)'); ylabel(ax, 'scale j');

for scl = 1:length(OPTIONS.wavelet.selected_scales)

sj = OPTIONS.wavelet.selected_scales(scl);
bj = find(OPTIONS.automatic.Modality(ii).selected_jk(2,:)==sj);
tt = OPTIONS.automatic.Modality(ii).selected_jk(6,bj);
vv = OPTIONS.automatic.selected_values{ii}(2,bj);
tv = OPTIONS.automatic.selected_values{ii}(3,bj);

if isempty(tt)
% if there is no selected box for that scale, we skip.
continue;
end

hhh = uitab(hptab,'title', [' scale ' num2str(OPTIONS.wavelet.selected_scales(scl)) ' ']);
hpc = uipanel( 'Parent', hhh, ...
'Units', 'Normalized', ...
'Position', [0.01 0.01 0.98 0.98], ...
'FontWeight', 'demi');

set(hpc,'Title', ' Wavelet coefficients (% of power) ', 'FontSize', 8);
ax = axes('parent',hpc, 'outerPosition',[0.01 0.01 0.98 0.98]);

hold(ax,'on')
stem(ax, tt(tv==0), vv(tv==0), 'x', 'filled', 'markersize', 8, 'MarkerFaceColor', 'black', 'Color','black');
stem(ax, tt(tv==1), vv(tv==1), 'o', 'filled', 'markersize', 8, 'MarkerFaceColor', 'red', 'Color','red');
hold(ax,'off')
end

drawnow
end

warning('off','all');
hp = obj.hfig;
hptab = obj.hfigtab;

if ~strcmp(obj.data_type,'discrete_wavelet')
disp('NO WAVELET DISPLAY');
return
end

for ii = 1:length(OPTIONS.mandatory.DataTypes)

onglet{ii} = uitab(hptab,'title',OPTIONS.mandatory.DataTypes{ii});
function sBox = create_rectangles(obj, OPTIONS, iMod)

Tmin = obj.t0 - (obj.info_extension.start-1)/OPTIONS.automatic.sampling_rate;
Tmax = Tmin + (size(obj.data{ii},2)-1)/OPTIONS.automatic.sampling_rate;
J = size(OPTIONS.automatic.scales,2);
N = size(obj.data{ii},2);
T = Tmax-Tmin;
Tmax = Tmin + (size(obj.data{iMod},2)-1)/OPTIONS.automatic.sampling_rate;
N = size(obj.data{iMod},2);
T = Tmax - Tmin;
e = 0.05;
MMM = colormap(jet(size(OPTIONS.automatic.selected_values{iMod},2)));
MMM = MMM(end:-1:1, :);
selection = OPTIONS.automatic.Modality(iMod).selected_jk;


hpc = uipanel('Parent', onglet{ii}, ...
'Units', 'Normalized', ...
'Position', [0.01 0.01 0.98 0.98], ...
'FontWeight','demi');

set(hpc,'Title', 'Time-frequency boxes', 'FontSize', 8);

ax = axes( 'parent', hpc, ...
'outerPosition', [0.01 0.01 0.98 0.98], ...
'YTick', 0.5:J-0.5, ...
'YTickLabel', num2cell(1:J));

axis(ax, 'fill');
box(ax,'on');
xlim(ax,[Tmin Tmax]);
ylim(ax,[0 J]);
xlabel(ax,'time (s)');
ylabel(ax,'scale j');

hold(ax, 'on');
sBox = struct('Vertices', [],'Faces',[], 'FaceVertexCData', [],'FaceColor','flat', 'EdgeColor', 'none');
iBox = 1;

MMM = colormap(jet(size(OPTIONS.automatic.selected_values{ii},2)));
selection = OPTIONS.automatic.Modality(ii).selected_jk;
for sj = 1:max(OPTIONS.automatic.Modality(iMod).selected_jk(2,:))

for b=1:length(OPTIONS.automatic.selected_values{ii})
l = T/N*2^selection(2,b);
rectangle('parent', ax, ...
'Position', [Tmin+(selection(3,b)-1)*l,selection(2,b)-1+e,l,1-2*e],...
'FaceColor', MMM(b,end:-1:1), ...
'LineStyle','none');
end
hold off

for scl = 1:length(OPTIONS.wavelet.selected_scales)
bj = find(OPTIONS.automatic.Modality(iMod).selected_jk(2,:)==sj);
tt = OPTIONS.automatic.Modality(iMod).selected_jk(6,bj);

sj = OPTIONS.wavelet.selected_scales(scl);
bj = find(OPTIONS.automatic.Modality(ii).selected_jk(2,:)==sj);
tt = OPTIONS.automatic.Modality(ii).selected_jk(6,bj);
vv = OPTIONS.automatic.selected_values{ii}(2,bj);
tv = OPTIONS.automatic.selected_values{ii}(3,bj);

if isempty(tt)
% if there is no selected box for that scale, we skip.
continue;
end

box_length = T/N*2^sj;
box_width = 1-2*e;

hhh = uitab(hptab,'title', [' scale ' num2str(OPTIONS.wavelet.selected_scales(scl)) ' ']);
hpc = uipanel( 'Parent', hhh, ...
'Units', 'Normalized', ...
'Position', [0.01 0.01 0.98 0.98], ...
'FontWeight','demi');
[val, I] = sort(selection(3,bj));

set(hpc,'Title', ' Wavelet coefficients (% of power) ', 'FontSize', 8);
ax = axes('parent',hpc, 'outerPosition',[0.01 0.01 0.98 0.98]);
color_scale = MMM(bj, :);
color_scale = color_scale(I, :);

for b=1:length(I)

hold(ax,'on')
stem(ax, tt, vv,'xk','filled','markersize',8);
plot(ax, tt(tv==1), vv(tv==1),'xr');
hold(ax,'off')
end
sBox.Vertices(end+1, :) = [Tmin + ((val(b)-1) * box_length), sj - box_width] + [ 0, 0];
sBox.Vertices(end+1, :) = [Tmin + ((val(b)-1) * box_length), sj - box_width] + [ box_length, 0];
sBox.Vertices(end+1, :) = [Tmin + ((val(b)-1) * box_length), sj - box_width] + [ box_length, box_width];
sBox.Vertices(end+1, :) = [Tmin + ((val(b)-1) * box_length), sj - box_width] + [ 0, box_width];

drawnow
pause(1)
end
sBox.Faces(end+1, :) = [iBox, iBox + 1, iBox + 2 , iBox + 3, iBox];
iBox = iBox + 4;

sBox.FaceVertexCData(end+1, :) = color_scale(b, :);
end
end
end
35 changes: 29 additions & 6 deletions gui/panel_brainentropy.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1224,11 +1224,22 @@ function setOptions(OPTIONS)

%% ===== CANCEL BUTTON =====
function ButtonCancel_Callback()

hfig = findall(0, 'Type', 'Figure', 'Tag', 'wMEM-TF');
if ~isempty(hfig)
close(hfig)
end
gui_hide(panelName);
end

%% ===== OK BUTTON =====
function ButtonOk_Callback(varargin)
function ButtonOk_Callback(varargin)

hfig = findall(0, 'Type', 'Figure', 'Tag', 'wMEM-TF');
if ~isempty(hfig)
close(hfig)
end

% Release mutex and keep the panel opened
bst_mutex('release', panelName);
be_print_best(OPTIONS);
Expand Down Expand Up @@ -1821,15 +1832,29 @@ function clean_bsl_info()
end

function VisualizeOscillation()
hfig = [];

if length(MEMglobal.DataToProcess) == 1
hfig = findall(0, 'Type', 'Figure', 'Tag', 'wMEM-TF');
if isempty(hfig) || length(hfig) > 1
hfig = uifigure("Name", "Time-Frequency Representation",'Tag', 'wMEM-TF', "Position", [725 198 560 420]);
else
clf(hfig)
end
end

for iData = 1:length(MEMglobal.DataToProcess)

if isempty(hfig)
hfig = uifigure("Name", "Time-Frequency Representation",'Tag', 'wMEM-TF', "Position", [725 198 560 420]);
end

for iData = 1:length(MEMglobal.DataToProcess)
hfig = uifigure("Name", "Time-Frequency Representation", "Position", [725 198 560 420]);
hfigtab = uitabgroup(hfig, "Position",[0 0 hfig.Position(3) hfig.Position(4)]);

try
focus(hfig)
catch
figure(hfig);

end

sInput = in_bst_data(MEMglobal.DataToProcess{iData});
Expand Down Expand Up @@ -1916,8 +1941,6 @@ function VisualizeOscillation()
[OPTIONS_wav, obj] = be_wdata_preprocessing(obj, OPTIONS_wav);
be_display_time_scale_boxes(obj, OPTIONS_wav);
end
drawnow; % force event queue to process
pause(1)
end

end
Expand Down