% положение различных объектов абсол. пав-не
wall_corners_X=[1.05 4.55 4.55 1.05 1.05];
wall_corners_Y=[0.25 0.25 3.75 3.75 0.25];
% центр абсол. столба (#1)
pillar1_center=[wall_corners_X(1)+0.78 wall_corners_Y(1)+0.84];
% центр столба #2
pillar2_center=[wall_corners_X(3)-1.07 wall_corners_Y(3)-0.75];
% место датчика POS-1 на полке
place_pos1=[wall_corners_X(3)-0.15 wall_corners_Y(3)-0.80];

% ноябрь
tmp=load('pickets_time_XY_dFmean_dFsd_2022Nov.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(1)
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(wall_corners_X,wall_corners_Y,'r','LineWidth',1)
plot(pillar1_center(1),pillar1_center(2),'kd','MarkerFaceColor','w','MarkerSize',10)
plot(place_pos1(1),place_pos1(2),'kd','MarkerFaceColor','w','MarkerSize',10)
hold off
ho=text(pillar1_center(1)+0.00,pillar1_center(2)+0.00,'1','Color','k','FontWeight','Bold');
set(ho,'HorizontalAlignment','center')
ho=text(place_pos1(1)+0.00,place_pos1(2)+0.00,'2','Color','k','FontWeight','Bold');
set(ho,'HorizontalAlignment','center')
% маркировка панели
h0=text(5.0,-0.3,'(a)','FontWeight','Bold','FontSize',10);

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

% апрель

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

x1=0; x2=4.8; y1=0; y2=4.0;

[xi,yi]=meshgrid(x1:0.01:x2,y1:0.01:y2);
dfi=griddata(x,y,dfm,xi,yi,'cubic');
subplot(122)
% V=[[-50:5:-10] [-12:2:12] [15:5:25]];
% [cs,h]=contourf(xi,yi,dfi,V); grid on, colorbar
% V=linspace(-62,+32,51);
% [cs,h]=contourf(xi,yi,dfi,V); grid on, colorbar
V=linspace(-65,+20,51);
[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('03.04.2023, Paratunka, Absolute hut, H = 1 m') 
axis equal

% стены пав-на
hold on
plot(wall_corners_X,wall_corners_Y,'r','LineWidth',1)
plot(pillar1_center(1),pillar1_center(2),'kd','MarkerFaceColor','w','MarkerSize',10)
plot(place_pos1(1),place_pos1(2),'kd','MarkerFaceColor','w','MarkerSize',10)
hold off
ho=text(pillar1_center(1)+0.00,pillar1_center(2)+0.00,'1','Color','k','FontWeight','Bold');
set(ho,'HorizontalAlignment','center')
ho=text(place_pos1(1)+0.00,place_pos1(2)+0.00,'2','Color','k','FontWeight','Bold');
set(ho,'HorizontalAlignment','center')
% маркировка панели
h0=text(5.0,-0.3,'(b)','FontWeight','Bold','FontSize',10);

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

% return

% рисуем распределение разностей двух съёмок
% ноябрь 2022
tmp=load('pickets_time_XY_dFmean_dFsd_2022Nov.dat');
Np1=tmp(:,1); % номер пикета
% координаты (в м) относительно юго-западного угла (пикет 01)
x1=tmp(:,8); y1=tmp(:,9); dfm1=tmp(:,10); % ср. df
[xi1,yi1]=meshgrid(min(x1):0.01:max(x1),min(y1):0.01:max(y1));
dfi1=griddata(x1,y1,dfm1,xi1,yi1,'cubic');

% ноябрь 2022
tmp=load('pickets_time_XY_dFmean_dFsd_2023Apr.dat');
Np2=tmp(:,1); % номер пикета
% координаты (в м) относительно юго-западного угла (пикет 01)
x2=tmp(:,8); y2=tmp(:,9); dfm2=tmp(:,10); % ср. df
[xi2,yi2]=meshgrid(min(x2):0.01:max(x2),min(y2):0.01:max(y2));
dfi2=griddata(x2,y2,dfm2,xi2,yi2,'cubic');

% смотрим разности
% интерполир. измерения в апреле на все точки (ноября)
dfi2=griddata(x2,y2,dfm2,xi1,yi1,'cubic');
dfi_Apr_Nov=dfi2-dfi1; % разность Apr-Nov по интерполир. сетке
figure(2)
[cs,h]=contourf(xi1,yi1,dfi_Apr_Nov,51); grid on; colorbar
hold on, plot(x,y,'ko','MarkerFaceColor','g','MarkerSize',3), hold off
xlabel('W -> E, m'), ylabel('S -> N, m')
% title('03.04.2023, Paratunka, Absolute hut, H = 1 m') 
axis equal
% стены пав-на
% стены пав-на
hold on
plot(wall_corners_X,wall_corners_Y,'r','LineWidth',1)
plot(pillar1_center(1),pillar1_center(2),'kd','MarkerFaceColor','w','MarkerSize',10)
plot(place_pos1(1),place_pos1(2),'kd','MarkerFaceColor','w','MarkerSize',10)
hold off
ho=text(pillar1_center(1)+0.00,pillar1_center(2)+0.00,'1','Color','k','FontWeight','Bold');
set(ho,'HorizontalAlignment','center')
ho=text(place_pos1(1)+0.00,place_pos1(2)+0.00,'2','Color','k','FontWeight','Bold');
set(ho,'HorizontalAlignment','center')
