地形图以及非线性画法的总结

s1 = load ('D:\task\projects\nnet\last_result\test12_output_nosmoth\1\8834_down.py');
s1 = reshape(s1,256,512);
s2=load ('D:\task\projects\nnet\last_result\test12_output_nosmoth\1\8834_up.py');
s2 = reshape(s2,256,512);
TD1 = abs((s1-s2)*2.62);
TD1 = imresize(TD1, [384 512]);
TD1(TD1<0)=0;
z_min = min(TD1(:))
z_max = max(TD1(:))


s11 = load ('D:\task\projects\nnet\last_result\test12_output_nosmoth\1\8715_down.py');
s11 = reshape(s11,256,512);
s12=load ('D:\task\projects\nnet\last_result\test12_output_nosmoth\1\8715_up.py');
s12 = reshape(s12,256,512);
TD = abs((s11-s12)*2.62);
TD = imresize(TD, [384 512]);
TD(TD<0)=0; 
z_min2 = min(TD(:))
z_max2 = max(TD(:))
TD=TD/z_max;
k = 1/log(11);
TD=k*log(10*TD+1)*80;
figure
image(TD),colormap(mycmap),caxis manual   
%image(TD),colormap(jet(80)),caxis manual   
colorbar
% colorbar('YTickLabel',[(z_min+z_max/8),(z_min+2*z_max/8),(z_min+3*z_max/8),(z_min+4*z_max/8),(z_min+5*z_max/8),(z_min+6*z_max/8),(z_min+7*z_max/8),(z_min+8*z_max/8)]) 
a = [(z_min+z_max/8),(z_min+2*z_max/8),(z_min+3*z_max/8),(z_min+4*z_max/8),(z_min+5*z_max/8),(z_min+6*z_max/8),(z_min+7*z_max/8),(z_min+8*z_max/8)];
b=a/z_max;
% 
z = (exp(b/k)-1)/10;
z = z*z_max;
colorbar('YTickLabel',z)

x1 = 173;
y1 = 175;
r1 = 512/12;
r2 = 512/12*3;
r3 = 512/12*6;
x2 = 365;
y2 = 219;
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',[x1+r1,y1,r3-r1,1],'Curvature',[0,0],'LineWidth',2),axis equal
% rectangle('Position',[x1-r3,y1,r3-r1,1],'Curvature',[0,0],'LineWidth',2),axis equal
% rectangle('Position',[x1,y1-r3,1,r3-r1],'Curvature',[0,0],'LineWidth',2),axis equal
% rectangle('Position',[x1,y1+r1,1,r3-r1],'Curvature',[0,0],'LineWidth',2),axis equal
% hold on
% plot([x1 x2],[y1 y2],'b','LineWidth',2),axis([0,512,0,384]);
%%%%%%%%%第二分区
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);
x9=x(1);
y9=y(1);
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);
x10=x(1);
y10=y(1);
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')
axis('Position',[0,12,0,9]);
axis equal tight