地形图坐标的非线性变化

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


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=TD/z_max;
z_max2 = max(TD(:))
z_min2 = min(TD(:))
k = 1/log(26)
TD=k*log(25*TD+1)*80;
figure
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)/25
z = z*z_max
colorbar('YTickLabel',z)
%colorbar('YTickLabel',[z_min+z_max/8:z_max/8:z_max]) 
set(gca, 'XTick', 0:512/12:512); 
set(gca, 'XTickLabel', [0:12]);
set(gca, 'YTick', 0:384/9:384); 
set(gca, 'YTickLabel', [0:9]);
axis equal tight
a = z_min+8*z_max/8
b = a/z_max
c = (exp(b/k)-1)/25
d= log(10)

e = exp(2.3026)