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

uperrua

利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
8 I2 g2 R5 M% A9 t/ s6 R有两种version、两种投影方式。分别需要不同的matlab函数读取。
1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
2、解压数据，得到eausgs2_0la.img
3、matlab程序：
) X# H# `0 g$ Y% n7 A" x& Dclear;close all;clc
figure
: v" r) E/ X, q# w: @worldmap china
" a9 x. h, Z' g, L- M, smstruct = gcm;
latlim = mstruct.maplatlimit;
lonlim = mstruct.maplonlimit;
. a+ s. x2 D1 F4 E) T5 t$ Cscalefactor = 2;
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
1 x; s+ z, x) p+ D1 i6 s% |2 k2 w: {& D$ V   scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
# Q! E, B* D3 w5 n' m; B1 tgeoshow('landareas.shp','FaceColor','none','EdgeColor','black')
load usgslulegend; colormap(cmap)
. j) ^, O% l# Cset(handlem('suRFace'),'CDataMapping','Direct'); tightmap
caxis([.5 24.5]); hcb = colorbar;
- r2 Z! g/ n) g3 g( v0 nset(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
- d3 n4 F  l* F5 X8 w    'Position',[.7 .1 .02 .8])
set(gca,'Position',[-.05 .1 .8 .8])

- Q# K8 L. G0 }/ J

& @; ^+ `, T- P8 t# _
# J6 X$ O1 q# |0 i- D

" h( T+ D$ h3 a6 t2 m9 ^; j3 F% r4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）

  L6 a! z! C; d' _6 v+ z$ P

sharkN

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