使用MCD12C1数据绘制全球土地覆盖类型图

使用MCD12C1数据绘制全球土地覆盖类型图

MCD12数据是利用Terra和Aqua两颗卫星遥感获得的土地覆盖类型数据，MCD12C1是其中分辨率约为0.05°的版本，适合做全球尺度的绘图。MCD12Q1为与之相似的但精度更高的数据。
该数据可以在USGS中下载，格式为HDF。
使用matlab(R2017b)。

clear;
clc;

首先读取数据（HDF格式）

Input_way1='F:/MCD12C1.A2011001.051.2012264191639.hdf';    
		%hdf文件在电脑中的保存位置
Input_way2='/MOD12C1/Data Fields/Majority_Land_Cover_Type_1'; 
		%第一种分类数据在hdf文件中的位置                                                                                       %输出路径
Landcover= hdfread(Input_way1,Input_way2, 'Index', {[1  1],[1  1],[3600  7200]});

将数据中的空值（值为255）改为17，方便设置颜色图

Landcover(Landcover==255)=17;
caxis([0 17]);          %设置显示颜色范围
mycolor=[           	%设置颜色图
    0.4 0.4 0.4;            % 0     Water
    0 98/255 65/255;        % 1     Evergreen Needleleaf forest	
    72/255 150/255 32/255;  % 2     Evergreen Broadleaf forest	
    0 160/255 107/255;      % 3     Deciduous Needleleaf forest	
    91/255 189/255 43/255;  % 4     Deciduous Broadleaf forest	
	131/255 199/255 93/255; % 5     Mixed forest	
    0 132/255 137/255;      % 6     Closed shrublands	
    110/255 195/255 201/255;% 7     Open shrublands	
    156/255 153/255 0;      % 8	    Woody savannas	
    252/255 245/255 78/255; % 9	    Savannas	
    243/255 194/255 70/255; % 10	Grasslands	
    160/255 149/255 196/255;% 11	Permanent 
    189/255 107/255 9/255;  % 12	Croplands	
    139/255 0 22/255;       % 13	Urban and built-up	
    236/255 135/255 14/255; % 14	Cropland/Natural vegetation mosaic	 
    115/255 136/255 193/255;% 15	Snow and ice	 
    170/255 135/255 184/255;% 16	Barren or sparsely vegetated	
    1 1 1;                  % 17	Fill Value/Unclassified
    ];     

使用设置好的颜色图绘图

    colormap(mycolor);
    y1=[0,3600];
    x1=[0,7200];
    axis([0 7200 0 3600]);
    image(x1,y1,Landcover);
    daspect([1 1 1]);

设置坐标轴

    xticks([1 900 1800 2700 3600 4500 5400 6300 7200])
    xticklabels({'180°','135°','90°','45°','0°','45°','90°','135°','180°'})
    yticks([1 900 1800 2700 3600])
    yticklabels({'90°','45°','0°','45°','90°'})
    title('全球土地覆盖类型图')
    

设置colorbar

 colorbar('Ticks',[0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5,10.5,11.5,...
		   12.5,13.5,14.5,15.5,16.5,17.5],...
          'TickLabels',{'水体','常绿针叶林','常绿阔叶林','落叶针叶林',...
          '落叶阔叶林','混交林','密集灌木','疏松灌木',...
          '多树荒原','荒原','草原','永久湿地','农田','城市与建成区',...
          '农田与天然植被相交','冰与雪','贫瘠','无数据或未分类'})

完成。