使用MCD12C1資料繪製全球土地覆蓋型別圖

MCD12資料是利用Terra和Aqua兩顆衛星遙感獲得的土地覆蓋型別資料，MCD12C1是其中解析度約為0.05°的版本，適合做全球尺度的繪圖。MCD12Q1為與之相似的但精度更高的資料。
該資料可以在USGS或NASA中下載，格式為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',{'水體','常綠針葉林','常綠闊葉林','落葉針葉林',...
          '落葉闊葉林','混交林','密集灌木','疏鬆灌木',...
          '多樹荒原','荒原','草原','永久溼地','農田','城市與建成區',...
          '農田與天然植被相交','冰與雪','貧瘠','無資料或未分類'})

完成。