Skip to content

Two colors for two conditions and ss threshold #209

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
duwadisudan wants to merge 1 commit into
base: development
Choose a base branch
from
Open

Two colors for two conditions and ss threshold #209

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
Binary file modified PlotProbe2GUI/PlotProbe2.fig
View file
Open in desktop
Binary file not shown.
117 changes: 100 additions & 17 deletions PlotProbe2GUI/PlotProbe2.m
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -22,7 +22,7 @@

% Edit the above text to modify the response to help PlotProbe2

% Last Modified by GUIDE v2.5 22-Mar-2022 16:55:43
% Last Modified by GUIDE v2.5 09-May-2024 07:29:21

% Return if snirf object was not passed
if isempty(varargin)
Expand Down Expand Up @@ -81,6 +81,10 @@ function PlotProbe2_OpeningFcn(hObject, eventdata, handles, varargin)
else
snirfObj = varargin{1};
end

% Update handles structure
guidata(hObject, handles);

dataTypeLabels = {};
dataTypeOrder = {};
measList = [];
Expand Down Expand Up @@ -122,7 +126,7 @@ function PlotProbe2_OpeningFcn(hObject, eventdata, handles, varargin)
snirfObj.data(v).dataTimeSeries(:,u) = snirfObj.data(v).dataTimeSeries(:,u) -(minAmp+ maxAmp)/2;
minAmp = minAmp-(minAmp+ maxAmp)/2;
maxAmp = maxAmp-(minAmp+ maxAmp)/2;
if strcmp(dataTypeLabel,'HRF HbO') || strcmp(dataTypeLabel,'HRF HbO') || strcmp(dataTypeLabel,'HRF HbO')
if strcmp(dataTypeLabel,'HRF HbO')
dataTypeOrder{end+1} = 'HRF HbX';
elseif contains(dataTypeLabel,'dOD')
dataTypeOrder{end+1} = 'dOD';
Expand Down Expand Up @@ -172,6 +176,8 @@ function PlotProbe2_OpeningFcn(hObject, eventdata, handles, varargin)
end
set(handles.edit_minDistForDisplay, 'String', min_dist);
set(handles.edit_maxDistForDisplay, 'String', max_dist);
defualtssFadeThres = 0;
set(handles.edit_ssFadeThres, 'String', num2str(defualtssFadeThres));

condition_names = {' All'};
for u = 1:length(snirfObj.stim)
Expand Down Expand Up @@ -211,11 +217,10 @@ function PlotProbe2_OpeningFcn(hObject, eventdata, handles, varargin)
display(handles);
guidata(hObject,handles)


function display(handles)

if isfield(handles,'data') & isfield(handles.data, 'snirfObj')
if isfield(handles, 'data') && isfield(handles.data, 'snirfObj')

snirfObj = handles.data.snirfObj;
measList = handles.data.measList;
sPos = snirfObj.probe.sourcePos2D;
Expand Down Expand Up @@ -332,6 +337,8 @@ function display(handles)

channel_min_dist = str2double(get(handles.edit_minDistForDisplay, 'String'));
channel_max_dist = str2double(get(handles.edit_maxDistForDisplay, 'String'));
ssFadeThres = str2double(get(handles.edit_ssFadeThres, 'String'));

% contents = cellstr(get(handles.listbox_selectConditions,'String'));
selected_conditions_index = get(handles.listbox_selectConditions,'Value');
data_index = 0;
Expand All @@ -353,14 +360,61 @@ function display(handles)
end
end
if data_index~= 0
SigCh = [1; 36; 48; 50];
for u = 1:length(snirfObj.data(data_index).measurementList)
activityConditionIndex = GetCondition(snirfObj.data(data_index).measurementList(u));
if any(selected_conditions_index == 1) || any(selected_conditions_index == activityConditionIndex+1)

% for HbO
color_condition_1 = [0.862, 0.078, 0.235]; % Red for condition 1
color_condition_2 = [1, 0, 1]; % Magenta for condition 2

% for HbR
color_condition_3 = [0, 0, 0.8]; % Blue for condition 1
color_condition_4 = [0, 1, 1]; % Cyan for condition 2

% for HbT
color_condition_5 = [0.133, 0.545, 0.133]; % Dark green for condition 1
color_condition_6 = [0.5, 1, 0]; % Light pastel green for condition 2

% before selecting conditions, initialize the colors as:

current_color_HbO = color_condition_1;
current_color_HbR = color_condition_3;
current_color_HbT = color_condition_5;

if ~~any(selected_conditions_index ==2, "all")
if any(selected_conditions_index == 2) && activityConditionIndex + 1 == 2
current_color_HbO = color_condition_1; % Red for condition 1
current_color_HbR = color_condition_3; % Blue for condition 2
current_color_HbT = color_condition_5; % for condition 3
elseif any(selected_conditions_index == 3) && activityConditionIndex + 1 == 3
current_color_HbO = color_condition_2; % Magenta for condition 2
current_color_HbR = color_condition_4; % Cyan for condition 2
current_color_HbT = color_condition_6; % for condition 2
end
end


srcIdx = GetSourceIndex(snirfObj.data(data_index).measurementList(u));
detIdx = GetDetectorIndex(snirfObj.data(data_index).measurementList(u));

detIdx = GetDetectorIndex(snirfObj.data(data_index).measurementList(u));
channel_dist = sqrt(sum((sourcePos3D(srcIdx,:) - detectorPos3D(detIdx ,:)).^2));

lineWidth = 0.5;

fade_factor = 0.7;

if channel_dist >= 0 && channel_dist <= ssFadeThres
% Use a faded color or different line style for short separation
% Light gray color for faded effect
current_color_HbO = current_color_HbO + (1-current_color_HbO)*fade_factor;
current_color_HbR = current_color_HbR + (1-current_color_HbR)*fade_factor;
current_color_HbT = current_color_HbT + (1-current_color_HbT)*fade_factor;

else

end

if channel_dist >= channel_min_dist & channel_dist <= channel_max_dist
xa = (sPos(srcIdx,1) + dPos(detIdx ,1))/2 - axXoff;
ya = (sPos(srcIdx,2) + dPos(detIdx ,2))/2 - axYoff;
Expand Down Expand Up @@ -395,6 +449,9 @@ function display(handles)
end
xT = xa-axWid/4 + axWid*((t-minT)/(maxT-minT))/2;
xyas = [xyas; [xa, ya]];



Avg = snirfObj.data(data_index).dataTimeSeries(:,u);
% minAmp=squeeze(min(min(Avg)));
% maxAmp=squeeze(max(max(Avg)));
Expand All @@ -405,25 +462,26 @@ function display(handles)
% cmin = min(AvgT);
% cmax = max(AvgT);
% AvgT = AvgT-(cmin+cmax)/2;



if any(contains(selected_display_activities,dataTypeLabel)) & contains(dataTypeLabel,'HbO')
plot( xT, AvgT,'color',color(1,:));
plot(xT, AvgT, 'color', current_color_HbO, 'LineWidth', lineWidth);
elseif any(contains(selected_display_activities,dataTypeLabel)) & contains(dataTypeLabel,'HbR')
plot( xT, AvgT,'color',color(2,:));
plot( xT, AvgT,'color', current_color_HbR,'LineWidth', lineWidth);
elseif any(contains(selected_display_activities,dataTypeLabel)) & contains(dataTypeLabel,'HbT')
plot( xT, AvgT,'color',color(3,:));
plot( xT, AvgT,'color',current_color_HbT,'LineWidth', lineWidth);
elseif any(contains(selected_display_activities,dataTypeLabel))
plot( xT, AvgT,'color',color(4,:));
plot( xT, AvgT,'color',color(4,:),'LineWidth', lineWidth);
end
end
end
end
end
hold off
if get(handles.checkbox_axisImage,'Value')
%if get(handles.checkbox_axisImage,'Value')
axis image
end
%end
end


Expand Down Expand Up @@ -528,7 +586,7 @@ function edit_Xscale_Callback(hObject, eventdata, handles)

% Hints: get(hObject,'String') returns contents of edit_Xscale as text
% str2double(get(hObject,'String')) returns contents of edit_Xscale as a double

axis image
display(handles)


Expand All @@ -553,7 +611,7 @@ function edit_Yscale_Callback(hObject, eventdata, handles)

% Hints: get(hObject,'String') returns contents of edit_Yscale as text
% str2double(get(hObject,'String')) returns contents of edit_Yscale as a double

axis image
display(handles)


Expand Down Expand Up @@ -768,7 +826,6 @@ function checkbox_axisImage_Callback(hObject, eventdata, handles)
% handles structure with handles and user data (see GUIDATA)

% Hint: get(hObject,'Value') returns toggle state of checkbox_axisImage

display(handles)


Expand Down Expand Up @@ -824,3 +881,29 @@ function listbox_selectActivity_CreateFcn(hObject, eventdata, handles)
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end



function edit_ssFadeThres_Callback(hObject, eventdata, handles)
% hObject handle to edit_ssFadeThres (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of edit_ssFadeThres as text
% str2double(get(hObject,'String')) returns contents of edit_ssFadeThres as a double
display(handles)


% --- Executes during object creation, after setting all properties.
function edit_ssFadeThres_CreateFcn(hObject, eventdata, handles)
% hObject handle to edit_ssFadeThres (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows.
% See ISPC and COMPUTER.
if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor'))
set(hObject,'BackgroundColor','white');
end