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

uperrua

利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
/ g* ^: D3 i8 s' e8 @# l# o2 ~2 ^& ^0 G有两种version、两种投影方式。分别需要不同的matlab函数读取。
" x0 P. M% m6 _# H# P1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
2、解压数据，得到eausgs2_0la.img
% A, B3 a" K* P! j) E# z5 D) a) n8 ]3、matlab程序：
clear;close all;clc
figure
' ^1 y4 ]) s' j( i) u& e3 k0 F0 jworldmap china
mstruct = gcm;
; w9 l$ h' x. K, }9 Xlatlim = mstruct.maplatlimit;
1 Y- E, i+ z) {4 t' I% _lonlim = mstruct.maplonlimit;
7 ]3 U* G% H7 r/ b' iscalefactor = 2;
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
$ Q, Y& r# w/ Y; ~( @   scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
" |6 C: T" D" h/ Ygeoshow('landareas.shp','FaceColor','none','EdgeColor','black')
2 Z5 e, W3 y1 `1 o) e" eload usgslulegend; colormap(cmap)
3 `/ u7 E! l1 y# A! V0 y- |set(handlem('suRFace'),'CDataMapping','Direct'); tightmap
caxis([.5 24.5]); hcb = colorbar;
& H% Z( u" y4 Y9 h4 S6 @8 Iset(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
    'Position',[.7 .1 .02 .8])
set(gca,'Position',[-.05 .1 .8 .8])


: |! ?7 ]3 ^) C

- L# d" `) w9 J
+ d! ^7 F+ Q* }' N# _2 C
! q  u: C0 V: Y, I2 K4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）

$ \! Y% K0 Z' v! h

sharkN

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