tmp=load('pickets_time_XY_dFmean_dFsd.dat');
Np=tmp(:,1); % номер пикета
% координаты (в м) относительно юго-западного угла (пикет 01)
x=tmp(:,8); y=tmp(:,9); 
dfm=tmp(:,10); % ср. df
dfs=tmp(:,11); % std. df

[xi,yi]=meshgrid(min(x):0.01:max(x),min(y):0.01:max(y));
dfi=griddata(x,y,dfm,xi,yi,'cubic');
figure(2)
subplot(1,2,1)
% V=[[-50:5:-10] [-12:2:12] [15:5:25]];
% V=linspace(-12,+25,51);
% [cs,h]=contourf(xi,yi,dfi,V); grid on, colorbar
V=linspace(-65,+36,51); % полный диапазон
V=linspace(-65,+20,51); % для сравнения с 03.04.2023
[cs,h]=contourf(xi,yi,dfi,V); grid on, colorbar
% clabel(cs, h, 'labelspacing',100,'FontSize',8);
hold on, plot(x,y,'ko','MarkerFaceColor','g','MarkerSize',3), hold off
xlabel('W -> E, m'), ylabel('S -> N, m')
% title('07-08.11.2022, Paratunka, Absolute hut, H = 1 m') 
axis equal

% стены пав-на
hold on
plot([1.05 4.55 4.55 1.05 1.05],[0.25 0.25 3.75 3.75 0.25],'r','LineWidth',2)
hold off

% text(x-0.1,y+0.1,num2str(Np,'%02d'),'Color','r','FontWeight','Bold')
