利用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数据。
2、解压数据，得到eausgs2_0la.img
5 j- B) o% A& s7 ~3、matlab程序：
* k1 N: k: w% C( h1 @clear;close all;clc
figure
worldmap china
  h! g5 A) T2 t! ?; B1 V, k2 Smstruct = gcm;
9 X% }$ _' `; Tlatlim = mstruct.maplatlimit;
5 p- X8 e1 T+ T" y9 Alonlim = mstruct.maplonlimit;
scalefactor = 2;
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
9 q3 k! y7 u+ e+ v( H( D3 k9 \   scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
geoshow('landareas.shp','FaceColor','none','EdgeColor','black')
# |3 q. W$ u4 A7 w; _3 G% wload usgslulegend; colormap(cmap)
set(handlem('suRFace'),'CDataMapping','Direct'); tightmap
caxis([.5 24.5]); hcb = colorbar;
% q: G( b5 ?# Lset(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
    'Position',[.7 .1 .02 .8])
$ @. X- o* k* e9 j' @, Y- Eset(gca,'Position',[-.05 .1 .8 .8])
0 u, m; @2 d7 s6 r

% y  X3 R9 a5 x  r5 n- D1 z
. P" V) c3 k' U$ z. Q5 v+ f


4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）
$ j5 |& C9 K: y7 F5 c; Z+ t
3 G! y! V1 U9 T2 m: [( ~

sharkN

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