function lemi=lemi_gps_correct(lemi);
% Корректируем пропуски в данных LEMI во времЯ GPS-измерений
%     каждый час в 4 мин. 22 сек. - данные отсутствуют
% ЗаполнЯем линейной интерполЯцией между соседними 
%     измерениЯми в 21 и 23 сек.
%
% Вызов:   lemi=lemi_gps_correct(lemi);

% изменен анализ при наличии измерениЯ в HH:04:22

hh=lemi(:,1); mm=lemi(:,2); ss=lemi(:,3);
dx=lemi(:,4); dy=lemi(:,5); dz=lemi(:,6);
t1=lemi(:,7); t2=lemi(:,8);

n=length(dx);
% ищем индексы длЯ строк с mm=4 и ss=21
ki=find(mm==4 & ss==21); 
% проверЯем, всегда ли были пропуски hh:04:22
k22=find(mm==4 & ss==22);
% исключаем случаи с измерением в hh:04:22
if ~isempty(k22),
   for k=1:length(k22), 
      ki=ki(ki~=(k22(k)-1)); 
   end
%    ki=ki(ki~=(k22-1));
end
nc=length(ki);

dxc=zeros(nc,1); dyc=zeros(nc,1); dzc=zeros(nc,1);
hhc=zeros(nc,1); mmc=ones(nc,1)*4; ssc=ones(nc,1)*22;
t1c=zeros(nc,1); t2c=zeros(nc,1);
for k=1:nc,
    dxc(k)=mean(dx(ki(k):ki(k)+1));
    dyc(k)=mean(dy(ki(k):ki(k)+1));
    dzc(k)=mean(dz(ki(k):ki(k)+1));
    hhc(k)=hh(ki(k));
    t1c(k)=mean(t1(ki(k):ki(k)+1));
    t2c(k)=mean(t2(ki(k):ki(k)+1));
end
hh=[hh;hhc]; mm=[mm;mmc]; ss=[ss;ssc]; tl=dms_deg(hh,mm,ss);
[tl,ii]=sort(tl);
dx=[dx;dxc]; dy=[dy;dyc]; dz=[dz;dzc]; t1=[t1;t1c]; t2=[t2;t2c];
dx=dx(ii); dy=dy(ii); dz=dz(ii); t1=t1(ii); t2=t2(ii);
hh=hh(ii); mm=mm(ii); ss=ss(ii);

lemi=[hh mm ss dx dy dz t1 t2];
