From 4439b9de1a1922cd3b58255bbf2bfe4ee304d6b9 Mon Sep 17 00:00:00 2001 From: Oliver Lemke <olemke@gmail.com> Date: Tue, 14 Jan 2025 11:42:30 +0100 Subject: [PATCH] Take median around i --- atovs_read_find.m | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/atovs_read_find.m b/atovs_read_find.m index 76062c4..a0eb5e9 100644 --- a/atovs_read_find.m +++ b/atovs_read_find.m @@ -144,10 +144,10 @@ y2 = double(squeeze(data.dsvcounts21)); y3 = double(squeeze(data.dsvcounts31)); y4 = double(squeeze(data.dsvcounts41)); for i = intr_start:intr_end - A(1,i) = median(y1(i:i+10))-(median(y1(i-70:i-60))+median(y1(i+60:i+70)))/2; - A(2,i) = median(y2(i:i+10))-(median(y2(i-70:i-60))+median(y2(i+60:i+70)))/2; - A(3,i) = median(y3(i:i+10))-(median(y3(i-70:i-60))+median(y3(i+60:i+70)))/2; - A(4,i) = median(y4(i:i+10))-(median(y4(i-70:i-60))+median(y4(i+60:i+70)))/2; + A(1,i) = median(y1(i-5:i+5))-(median(y1(i-70:i-60))+median(y1(i+60:i+70)))/2; + A(2,i) = median(y2(i-5:i+5))-(median(y2(i-70:i-60))+median(y2(i+60:i+70)))/2; + A(3,i) = median(y3(i-5:i+5))-(median(y3(i-70:i-60))+median(y3(i+60:i+70)))/2; + A(4,i) = median(y4(i-5:i+5))-(median(y4(i-70:i-60))+median(y4(i+60:i+70)))/2; end M(1) = max(A(1,:)); M(2) = max(A(2,:)); -- GitLab