Right click on a graph can now be used for manual QC. It differs from…

… the left click in that instead of highlighting what is included in the box, it highlights what is outside the box.

Graph/DepthProfile/highlightDepthProfileGeneric.m

@@ -19,7 +19,8 @@
 %               data. If no data points lie within the highlight region, 
 %               an empty matrix is returned.
 %
-% Author: Paul McCarthy <[email protected]>
+% Author:       Paul McCarthy <[email protected]>
+% Contributor:  Guillaume Galibert <[email protected]>
 %
 
 %
@@ -39,4 +40,48 @@
 % along with this program.
 % If not, see <https://www.gnu.org/licenses/gpl-3.0.en.html>.
 %
-highlight = highlightTimeSeriesGeneric(region, data, variable, type);
+narginchk(4, 4);
+
+if ~isnumeric(region) || ~isvector(region) || length(region) ~= 4
+  error('region must be a numeric vector of length 4');
+end
+
+if ~ishandle(data),     error('data must be a graphics handle'); end
+if ~isstruct(variable), error('variable must be a struct');      end
+if ~ischar(type),       error('type must be a string');          end
+
+xdata = get(data(1), 'XData'); % data(1) retrieves the current graphic handles only, in case extra sample is selected
+ydata = get(data(1), 'YData');
+
+if iscell(xdata)
+   xdata = cell2mat(xdata)'; 
+   ydata = cell2mat(ydata)';
+end
+
+if strcmp(type, 'alt') % right click
+    % figure out indices of all data points outside the X range but within
+    % the Y range
+    xidx  = (xdata <  region(1)) | (xdata >  region(3));
+    yidx  = (ydata >= region(2)) & (ydata <= region(4));
+else
+    % figure out indices of all data points within the X and Y range
+    xidx  = (xdata >= region(1)) & (xdata <= region(3));
+    yidx  = (ydata >= region(2)) & (ydata <= region(4));
+end
+
+% figure out indices of all the points to be highlighted
+idx = xidx & yidx;
+
+if ~any(idx)
+    % return nothing if no points to plot
+    highlight = [];
+else
+    % create the highlight
+    highlight = line(xdata(idx),ydata(idx), ...
+        'UserData',        idx, ...
+        'Parent',          gca, ...
+        'LineStyle',       'none', ...
+        'Marker',          'o', ...
+        'MarkerEdgeColor', 'white', ...
+        'MarkerFaceColor', 'white');
+end

Graph/TimeSeries/highlightTimeSeriesGeneric.m

@@ -55,20 +55,16 @@
    ydata = cell2mat(ydata)';
 end
 
-% on right click highlight, only highlight 
-% unflagged data points in the region
-if strcmp(type, 'alt')    
-    f = variable.flags;
-    f = f == 0;
-
-    xdata = xdata(f);
-    ydata = ydata(f);
+if strcmp(type, 'alt') % right click
+    % figure out indices of all data points within the X range but outside the Y range
+    xidx  = (xdata >= region(1)) & (xdata <= region(3));
+    yidx  = (ydata <  region(2)) | (ydata >  region(4));
+else
+    % figure out indices of all data points within the X and Y range
+    xidx  = (xdata >= region(1)) & (xdata <= region(3));
+    yidx  = (ydata >= region(2)) & (ydata <= region(4));
 end
 
-% figure out indices of all data points within the range
-xidx  = (xdata >= region(1) & xdata <= region(3));
-yidx  = (ydata >= region(2) & ydata <= region(4));
-
 % figure out indices of all the points to be highlighted
 idx = xidx & yidx;
 

Graph/TimeSeries/highlightTimeSeriesTimeDepth.m

@@ -61,21 +61,18 @@
 Ydata = repmat(ydata', [size(xdata), 1]);
 clear xdata ydata;
 
-% figure out indices of all data points within the range
-X = ((Xdata >= region(1)) & (Xdata <= region(3)));
-Y = ((Ydata >= region(2)) & (Ydata <= region(4)));
+if strcmp(type, 'alt') % right click
+    % figure out indices of all data points within the X range but outside the Y range
+    X = (Xdata >= region(1)) & (Xdata <= region(3));
+    Y = (Ydata <  region(2)) | (Ydata >  region(4));
+else
+    % figure out indices of all data points within the X and Y range
+    X = (Xdata >= region(1)) & (Xdata <= region(3));
+    Y = (Ydata >= region(2)) & (Ydata <= region(4));
+end
 idx = X & Y;
 clear X Y;
 
-% on right click, only highlight unflagged points
-if strcmp(type, 'alt')
-
-  % see if any points in the region haven't been flagged
-  idx = (variable.flags == 0) & idx;
-  if ~any(any(variable.flags)), return; end
-
-end
-
 if ~any(any(idx))
     highlight = [];
 else

Graph/TimeSeries/highlightTimeSeriesTimeFrequencyDirection.m

@@ -50,21 +50,18 @@
 Xdata = get(data, 'XData');
 Ydata = get(data, 'YData');
 
-% figure out indices of all data points within the range
-X = ((Xdata >= region(1)) & (Xdata <= region(3)));
-Y = ((Ydata >= region(2)) & (Ydata <= region(4)));
+if strcmp(type, 'alt') % right click
+    % figure out indices of all data points outside the X and Y range
+    X = ((Xdata < region(1)) | (Xdata > region(3)));
+    Y = ((Ydata < region(2)) | (Ydata > region(4)));
+else
+    % figure out indices of all data points within the X and Y range
+    X = (Xdata >= region(1)) & (Xdata <= region(3));
+    Y = (Ydata >= region(2)) & (Ydata <= region(4));
+end
 idx = X & Y;
 clear X Y;
 
-% on right click, only highlight unflagged points
-if strcmp(type, 'alt')
-
-  % see if any points in the region haven't been flagged
-  idx = (variable.flags == 0) & idx;
-  if ~any(any(variable.flags)), return; end
-
-end
-
 if ~any(any(idx))
     highlight = [];
 else