% uiopen('C:\Users\ThinkBook 13s\Desktop\ISI\isi_LCT.tif',1)%'E:\lab_work\Gd-PAA_ISI\figures\20x_isi-LCT.tif',1)
I_isi=imread('C:\Users\86185\Desktop\ISI\isi_LCT.tif');I_ivi=imread('C:\Users\86185\Desktop\IVI\ivi_LCT.tif');
zoom=8;
II_isi=imresize(double(I_isi),1/zoom);II_ivi=imresize(double(I_ivi),1/zoom);   
r_isi=II_isi(:,:,1);r_ivi=II_ivi(:,:,1);
g_isi=II_isi(:,:,2);g_ivi=II_ivi(:,:,2);
b_isi=II_isi(:,:,3);b_ivi=II_ivi(:,:,3);

% plot3(r_isi(:),g_isi(:),b_isi(:),'r.');plot3(r_ivi(:),g_ivi(:),b_ivi(:),'b.');
grid on

ps_isi=size(II_isi);ps_ivi=size(II_ivi);
rr_isi=reshape(r_isi,1,ps_isi(1)*ps_isi(2));rr_ivi=reshape(r_ivi,1,ps_ivi(1)*ps_ivi(2));
gg_isi=reshape(g_isi,1,ps_isi(1)*ps_isi(2));gg_ivi=reshape(g_ivi,1,ps_ivi(1)*ps_ivi(2));
bb_isi=reshape(b_isi,1,ps_isi(1)*ps_isi(2));bb_ivi=reshape(b_ivi,1,ps_ivi(1)*ps_ivi(2));

scatter3(rr_isi,gg_isi,bb_isi,'.r') 
hold on
scatter3(rr_ivi,gg_ivi,bb_ivi,'.y')

hold on
xlabel('Red')
ylabel('Green')
zlabel('Blue')

% x=y=z
t=linspace(0,300);
plot3(t,t,t)

% z=f(x)+g(y)
% [rr_isi, gg_isi] = deal(linspace(-2, 2, 120));
% zMax = 2*sin(acos(rr_isi/2));
% [X Y] = meshgrid(rr_isi, gg_isi);
% Z = bsxfun(@plus, zMax, zMax');
% surf(X, Y, Z)