Найти все точки кривой, принадлежащие диапазону / matlab

figure(2)
[x1,y1,z1]=ginput(1)
[x2,y2,z2]=ginput(1)
 
 
figure(3)
m=round(y1):round(y2);% выбираем фрагмент изображения,в котором будем искать кривые изменения интенсивности сигнала от времени
n=round(x1):round(x2);
count =1;
z=1
for j=min(n):max(n) ;
    for i=min(m):max(m);%пробегаем по всем значениям внутри фрагмента
        subplot(length(n),length(m),count);
        curve1=squeeze(dataAy.dataAy(i,j,z,:));% выводим кривые изменения интенсивности
        plot(curve1)
        count =count+1;
    end
end
 
figure(4)
m=round(y1):round(y2);
n=round(x1):round(x2);
count =1;
z=1
for j=min(n):max(n) ;
    for i=min(m):max(m);
        subplot(length(n),length(m),count);
        curve1=squeeze(dataAy.dataAy(i,j,z,3:10));% выбираем часть кривой до провала сигнала
        plot(curve1) 
     
        precontrastMean(count)=mean(curve1); %средняя линия от каждой кривой
        precontrastSigma(count)=std(precontrastMean); % среднеквадратичное отклонение от среднего
 
hold on
     plot(precontrastSigma(count),'r+')
hold off
 
 threshold=3*precontrastSigma;%порог для преконтраста 3дельта
 maxLim=precontrastMean+threshold; %верхний предел
 minLim=precontrastMean-threshold; %нижний предел
 
 
 
     count =count+1;
    end
end
precontrastMean,
precontrastSigma

threshold=3*precontrastSigma;%порог для преконтраста 3дельта
 maxLim=precontrastMean+threshold; %верхний предел
 minLim=precontrastMean-threshold; %нижний предел
 
t = 1:length(sigma1);
L = (minLim<=sigma1)&(sigma1<=maxLim); % логическое условие попадания внутрь диапазона
plot( t(L), sigma1(L),'r',... % точки которые попали в диапазон
       t(~L), sigma1(~L),'b' ) % точки которые не попали в диапазон