利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据

uperrua

利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
有两种version、两种投影方式。分别需要不同的matlab函数读取。
1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
: [# M+ C' F: N! D1 x# i2、解压数据，得到eausgs2_0la.img
3、matlab程序：
clear;close all;clc
8 _' e6 z7 [* }# Mfigure
9 p- B8 j8 ~% I" @, D- pworldmap china
5 n& z' ?- R- B9 Y( f& p5 lmstruct = gcm;
# `, q0 ~# q  ^7 glatlim = mstruct.maplatlimit;
7 }; f- x7 |% z1 o5 M1 _3 Blonlim = mstruct.maplonlimit;
scalefactor = 2;
# y$ m1 K8 r2 J# b# O+ w: l9 D+ [[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
) @/ |; {# g- [   scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
6 y# W( ]* l6 k* L, `geoshow('landareas.shp','FaceColor','none','EdgeColor','black')
5 k2 d% g, s7 h. t; _1 G4 k" \load usgslulegend; colormap(cmap)
set(handlem('suRFace'),'CDataMapping','Direct'); tightmap
  H* @7 u  C  D; I# Mcaxis([.5 24.5]); hcb = colorbar;
1 W$ n0 X$ Z% _8 iset(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
    'Position',[.7 .1 .02 .8])
set(gca,'Position',[-.05 .1 .8 .8])


  J* m% Z& z% N. k$ D6 d$ f



4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）
1 }5 L! ]# W* N0 f% C4 W1 C

sharkN

利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。