QC_01_11c_TimeVsDiffDepthSpear.m

%% --------
% QC_01_11c_TimeVsDiffDepthSpear(MeasCat)

% Displays scatter plot of time vs depth

% Input: MeasCat - Measured catalogue


function QC_01_11c_TimeVsDiffDepthSpear(MeasCat, varargin)

WindowLength = 20;
FilterLength = 10;

close all

if ~isnumeric(MeasCat)
    error('Load catagoue in the correct format (see readme for description)')
end

PlotLimits = cell2mat(varargin);

MeasCat = MeasCat(MeasCat(:,9)>10,:);

%% definition of variables

% definition of event origin time, latitude, longitude, depth and magnitude
OriginTime = datenum(MeasCat(:,1), MeasCat(:,2), MeasCat(:,3), MeasCat(:,4), MeasCat(:,5), MeasCat(:,6));
EventLat = MeasCat(:, 7);
EventLon = MeasCat(:, 8);
EventDepth = MeasCat(:, 9);
EventMag = MeasCat(:, 10);

% freeing up the memory
clear MeasCat

%% main loop

NumOfEvents = numel(OriginTime);
NumOfSteps = NumOfEvents - WindowLength;

Spearman = zeros(NumOfSteps,1);
RowNum = (1:WindowLength)';

for n = 1:NumOfSteps
    
    DepthForLoop = EventDepth(n:n+WindowLength-1);
    Spearman(n) = corr(RowNum, DepthForLoop, 'Type', 'Spearman');
        
end

%% diferential depths and smoothing

filter = ones(FilterLength,1)/FilterLength;
PercentileFilt = conv(Spearman, filter, 'valid');
PercentileFilt2 = conv(PercentileFilt, filter, 'valid');

%% display plot

% color definition
FirstColor = [0 .47 .95];
SecondColor = [.95 .47 0];
ThirdColor = [.33 .66 0];
Grey = [.7 .7 .7];

% plot
MinTime = min(OriginTime);
MaxTime = max(OriginTime(1:NumOfSteps));

figure('name', 'Time Vs. Percentile diagram', 'Position', [100, 100, 800, 800])
subplot(2,1,1)
plot(OriginTime(1:NumOfSteps), Spearman, 'o', 'Color', FirstColor, 'MarkerSize', 4)
datetick('x',2)

if ~isempty(PlotLimits)
    ylim([PlotLimits(1) PlotLimits(3)])
    xlim([PlotLimits(2) PlotLimits(4)])
else
    xlim([MinTime MaxTime])
    ylim([min(Spearman) max(Spearman)])
end

title('Time Vs. Percentile', 'FontSize', 16, 'FontWeight', 'bold')
xlabel('Time')
ylabel('Depth difference between consecutive events')

subplot(2,1,2)
plot(OriginTime(FilterLength:NumOfSteps-FilterLength+1), PercentileFilt2, '-o','Color', SecondColor, 'MarkerSize', 2, 'LineWidth', 2)
datetick('x',2)

StDev = std(PercentileFilt2);
hold on
plot([MinTime MaxTime], [0 0], 'Color', FirstColor)
plot([MinTime MaxTime], [StDev StDev], ':', 'Color', FirstColor)
plot([MinTime MaxTime], [-StDev -StDev], ':', 'Color', FirstColor)
plot([MinTime MaxTime], [+2*StDev 2*StDev], '--', 'Color', FirstColor)
plot([MinTime MaxTime], [-2*StDev -2*StDev], '--', 'Color', FirstColor)

title('Time Vs. Spearman diagram - Low-pass filtered', 'FontSize', 16, 'FontWeight', 'bold')
xlabel('Time')
ylabel('Depth difference between consecutive events')

if ~isempty(PlotLimits)
    ylim([PlotLimits(1) PlotLimits(3)])
    xlim([PlotLimits(2) PlotLimits(4)])
else
    xlim([MinTime MaxTime])
    ylim([min(PercentileFilt2) max(PercentileFilt2)])
end

print(gcf,'CurrentFigures/QC_01_11c_TimeVsSpearman','-depsc', '-r300')