%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Project     : AiPT
% Module      : WinCamD_DisplayTiff
% Description : Display .tiff file generated by WinCamD.
% References  :
%
%Version|Date    |author|Description
%-------|--------|------|-------------------------------------------------------
%01.01  |20230209|L,YM  |Initial version
%-------|--------|------|-------------------------------------------------------
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clc; clear; %close all;
tic;
%% Coefficient set
FilePath = 'D:\Nick_Documents\v3_aston_save_and_load\MyPapers\v2023_StrongTurbulence\MatlabCode\profile_test\';
FileName = 'X_pol_FixExpose';
FigTiffX = imread([FilePath,FileName,'.tiff']);
FigTiffX = fliplr(flipud(FigTiffX));
% FigTiffX = imgaussfilt(FigTiffX,3);

FileName = 'Y_pol_FixExpose';
FigTiffY = imread([FilePath,FileName,'.tiff']);
FigTiffY = fliplr(flipud(FigTiffY));
% FigTiffY = imgaussfilt(FigTiffY,3);


FileName = 'A_pol_HalfExpose';
FigTiffA = imread([FilePath,FileName,'.tiff']);
FigTiffA = fliplr(flipud(FigTiffA));

%% choose here for x-pol/y-pol/all-pol
% FigTiff = imgaussfilt(FigTiffX,3);%Gaussian filter
% FigTiff = imgaussfilt(FigTiffY,3);%Gaussian filter
FigTiff = imgaussfilt(FigTiffA,3)*0.8;%Gaussian filter


centre = [750 ,482 ];%x,y
span   = [300 ,300 ];%x,y
FigTiff = FigTiff(centre(2)-span(2) : centre(2)+span(2)-1, ...
                  centre(1)-span(1) : centre(1)+span(1)-1);

%% Fig.position
Fig1 = figure;
Fig1.Units = 'inches';
Fig1.InnerPosition = [1, 1, 2, 2];

ax1 = axes;
ax1.InnerPosition = [0 0 1 1];%position with labels

%% plot
pcolor(FigTiff)%2D figure.

AxisMax = mean(max(FigTiffX,[],'all'),max(FigTiffY,[],'all'))*0.8;
caxis([0,AxisMax])
axis equal
axis off
shading interp%flat% faceted% 
set(gca,'FontSize',16);

toc;