clc;clear all;close all;

county_name={'达孜县','曲水县','尼木县','堆龙德庆县','林周县','墨竹工卡县','城关区','日喀则市','拉孜县','江孜县',...

   '谢通门县','南木林县','白朗县','乃东县','扎囊县','桑日县','贡嘎县','琼结县'};

max_h_county=xlsread('耕地分布上限.xlsx');

hf1= figure('Color',[1 1 1],'PaperPositionMode','auto','Units','Centimeter','Position',[1.5 1.5 15 18]);

t_year=1990:5:2010;

for i=1:6

   for j=1:3

        k=(i-1)*3+j;

        subplot(6,3,k);  

        plot(t_year,max_h_county(k,:),'b-o','linewidth',1.5);

        axis([1990 2010 min(max_h_county(k,:)) max(max_h_county(k,:))]);

        h=title(county_name(k),'FontName','Times New Roman','Fontsize',9);

        n=18-2.8*i;

        m=(j-1)*4.6+1.5;

        if k==1 || k==4 || k==7 || k==10 || k==13 || k==16

            ylabel('耕地上限(km)','Fontsize',8);

        end

        if k>=16

            xlabel('Year','Fontsize',8);

        end

        set(gca,'Units','Centimeter','Position',[m n 3.5 1.9],'Fontname','Times New Roman','TickLength',[0.02 0.02],'linewidth',0.8,'Fontsize',8)

   end

end

print(gcf,'-dtiff','-r600','耕地分布上限.tif');