画地形图以及几种图像的融合

 

s11 = load ('D:\task\projects\nnet\last_result\test12_output_nosmoth\1\150_down.py');
s11 = reshape(s11,256,512);
s12=load ('D:\task\projects\nnet\last_result\test12_output_nosmoth\1\150_up.py');
s12 = reshape(s12,256,512);
TD = abs((s11-s12)*2.62);

TD = imresize(TD, [384 512]);
TD(TD<0)=0; 
%TD = round(TD/6);
%TD = TD/2;
DIV=3;
figure
%load MyColormaps8.mat
%image(TD), colormap(choroid_map8)
image(TD/DIV), colormap(jet(80))
colorbar
z_min = min(TD(:));
z_max = max(TD(:));
if DIV==1
    colorbar('YTickLabel',[z_min+z_max/8:z_max/8:z_max]) 
else
    colorbar('YTickLabel',[z_max-z_max/DIV:(z_max-(z_max-z_max/DIV))/7:z_max])
end

%%%%%%%%%%%%%黄斑onh分区画三个圆
x1 = 50;
y1 = 190;
r1 = 512/12;
r2 = 512/12*3;
r3 = 512/12*6;
x2 = 400;
y2 = 215;
rectangle('Position',[x1-r1,y1-r1,2*r1,2*r1],'Curvature',[1,1],'LineWidth',2),axis equal
rectangle('Position',[x1-r2,y1-r2,2*r2,2*r2],'Curvature',[1,1],'LineWidth',2),axis equal
rectangle('Position',[x1-r3,y1-r3,2*r3,2*r3],'Curvature',[1,1],'LineWidth',2),axis equal
hold on
%%%%%%%%%%第一种分区
%rectangle('Position',[x+r1,y,r3-r1,2],'Curvature',[0,0],'LineWidth',1),axis equal
%rectangle('Position',[x-r3,y,r3-r1,2],'Curvature',[0,0],'LineWidth',1),axis equal
%rectangle('Position',[x,y-r3,1,r3-r1],'Curvature',[0,0],'LineWidth',2),axis equal
%rectangle('Position',[x,y+r1,1,r3-r1],'Curvature',[0,0],'LineWidth',2),axis equal
%%%%%%%%%第二分区
k = -((x2-x1)/(y2-y1));
% y=k*x+y1-k*x1
% plot(x,y)
syms x y;
eq1=(x-x1)^2+(y-y1)^2-r1^2;
eq1=subs(eq1);
eq2=y-k*x-y1+k*x1;
eq2=subs(eq2);
[x,y]=solve(eq1,eq2);
x=vpa(x,5);
y=vpa(y,5);
x3=x(1);
x4=x(2);
y3=y(1);
y4=y(2);

syms x y;
eq1=(x-x1)^2+(y-y1)^2-r3^2;
eq1=subs(eq1);
eq2=y-k*x-y1+k*x1;
eq2=subs(eq2);
[x,y]=solve(eq1,eq2);
x=vpa(x,5);
y=vpa(y,5);
x5=x(1);
x6=x(2);
y5=y(1);
y6=y(2);

plot([x1 x2],[y1 y2],'b','LineWidth',2),axis([0,512,0,384]);
plot([x4 x6],[y4 y6],'k','LineWidth',2),axis([0,512,0,384]);
plot([x3 x5],[y3 y5],'k','LineWidth',2),axis([0,512,0,384]);

%%%%%%%%%%%%%%%%%%%%
k = (y2-y1)/(x2-x1);
syms x y;
eq1=(x-x1)^2+(y-y1)^2-r1^2;
eq1=subs(eq1);
eq2=y-k*x-y1+k*x1;
eq2=subs(eq2);
[x,y]=solve(eq1,eq2);
x=vpa(x,5);
y=vpa(y,5);
x7=x(2);
y7=y(2);
k = (y2-y1)/(x2-x1);
syms x y;
eq1=(x-x1)^2+(y-y1)^2-r3^2;
eq1=subs(eq1);
eq2=y-k*x-y1+k*x1;
eq2=subs(eq2);
[x,y]=solve(eq1,eq2);
x=vpa(x,5);
y=vpa(y,5);
x8=x(2);
y8=y(2);

plot([x7 x8],[y7 y8],'k','LineWidth',2),axis([0,512,0,384]);
% 
% 
% 
% %%%%%%%%%%%%%%设置参数
set(gca, 'XTick', 0:512/12:512); 
set(gca, 'XTickLabel', [0:12]);
set(gca, 'YTick', 0:384/9:384); 
set(gca, 'YTickLabel', [0:9]);
set(gca,'clipping','on')
axes('Position',[0.00001,12,0.00001,9]);
axis equal tight
% 
% %figure
% %hist(TD(:)*6,300), axis([0 350  0 25000])
% 
%