使用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',{'水体','常绿针叶林','常绿阔叶林','落叶针叶林',...
'落叶阔叶林','混交林','密集灌木','疏松灌木',...
'多树荒原','荒原','草原','永久湿地','农田','城市与建成区',...
'农田与天然植被相交','冰与雪','贫瘠','无数据或未分类'})
完成。