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

uperrua

利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
: ?4 S6 }0 r5 J* ~' c有两种version、两种投影方式。分别需要不同的matlab函数读取。
1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
) w7 V' \7 p* h; k! ?# ~! B2、解压数据，得到eausgs2_0la.img
/ l- y* K& s+ I2 T3、matlab程序：
1 ?( ^+ k! q1 G* d5 Xclear;close all;clc
" |& @/ d4 u9 V! Q5 o  Yfigure
worldmap china
mstruct = gcm;
/ j- |2 ?& f$ \5 platlim = mstruct.maplatlimit;
8 b" N* k; W  ]) Q9 f/ Plonlim = mstruct.maplonlimit;
6 @/ b8 O  X' H0 @scalefactor = 2;
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
( E2 Y3 A& }: C; d. d5 P. r8 h8 z   scalefactor, latlim, lonlim);
8 g8 e. o/ n" Z' c9 J2 Kgeoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
; Y# J8 p) J" f2 ~. n* L6 t; ugeoshow('landareas.shp','FaceColor','none','EdgeColor','black')
' X0 y; n1 a2 d$ {2 N! gload usgslulegend; colormap(cmap)
$ b/ f7 d6 ?4 }2 b, wset(handlem('suRFace'),'CDataMapping','Direct'); tightmap
caxis([.5 24.5]); hcb = colorbar;
set(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
4 B7 t9 e6 @! ~5 _% G# s3 d6 r    'Position',[.7 .1 .02 .8])
6 x0 U) o6 E" ]set(gca,'Position',[-.05 .1 .8 .8])
1 p% _5 U3 u1 c! o
. C+ p1 l3 O1 P' I- R3 [& V
6 D% Y3 C" d& Q6 H+ C+ v" b5 F
9 [! b" r) ~2 K- W0 x5 p4 `
+ G! b3 H/ R. n' [. n3 [: [
+ K8 f% e6 j/ C- z* J  F
4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）
8 _3 [2 K- X3 O0 z& q4 O8 w

sharkN

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