clear
[xielv,R] = geotiffread('D:\data.tif');%读取线性趋势
info=geotiffinfo('D:\data.tif');
[m,n] = size(xielv);
xielv = reshape(xielv,m*n,1);
xielv_valid_zheng = xielv(xielv >= 0); %正变化趋势
[Num_zheng,xaxis_zheng] = hist(xielv_valid_zheng,50);
b2=bar(xaxis_zheng,Num_zheng,0.5,'FaceColor',[71/255,145/255,186/255],'EdgeColor',[71/255,145/255,186/255])
hold on
xielv_valid_fu = xielv(xielv < 0); %负变化趋势
[Num_fu,xaxis_fu] = hist(xielv_valid_fu,80);
b3=bar(xaxis_fu,Num_fu,0.5,'FaceColor',[218/255,103/255,81/255],'EdgeColor',[218/255,103/255,81/255])
% set(gca,'XLim', [-0.05 0.03])
% set(gca,'XTick', [-0.05 -0.04 -0.03 -0.02 -0.01 -0.005 0 0.005 0.01 0.015 0.02 0.03])
% set(gca,'Xticklabel',{'-0.05','-0.04','-0.03','-0.02','-0.01','-0.005','0','0.005','0.01','0.015','0.02','0.03'})
set(gca,'XLim', [-0.02 0.02])
set(gca,'XTick', [ -0.02 -0.015 -0.01 -0.005 0 0.005 0.01 0.015 0.02 ])
set(gca,'Xticklabel',{'-0.02','-0.015','-0.01','-0.005','0','0.005','0.01','0.015','0.02'})
set(gca,'YLim', [0 400000])
set(gca,'YTick',0:100000:400000)
set(gca,'Yticklabel',{'0','1','2','3','4'})
oldLabels = str2num(char(get(gca,'YTickLabel')));
Yscale = 10^5;
newLabels = num2str(oldLabels*Yscale);
set(gca,'YTickLabel',newLabels,'YTickLabelMode','auto');
set(gca,'FontSize',10)
xlabel('NDVI变化趋势 yr^-^1');
ylabel('频数');
legend([b2,b3],'正趋势','负趋势')
print(gcf,['D:\趋势直方图.tif'],'-dtiff','-r600');
close
