利用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数据。
) ?4 f' y8 m! s$ _& u+ m& \2、解压数据，得到eausgs2_0la.img
3、matlab程序：
clear;close all;clc
0 F  P1 @. i3 d2 b3 Jfigure
worldmap china
mstruct = gcm;
+ b) S. N# X: Glatlim = mstruct.maplatlimit;
lonlim = mstruct.maplonlimit;
+ Q' u* ?; Y0 T9 U. E+ Q# [scalefactor = 2;
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
8 V# |! l, m4 V0 Y   scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
  @, U8 ~+ a: r/ ngeoshow('landareas.shp','FaceColor','none','EdgeColor','black')
load usgslulegend; colormap(cmap)
2 S. `9 j$ R4 S1 o' Gset(handlem('suRFace'),'CDataMapping','Direct'); tightmap
caxis([.5 24.5]); hcb = colorbar;
9 u. g6 G6 M1 kset(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
* J& G, N, ~3 O    'Position',[.7 .1 .02 .8])
set(gca,'Position',[-.05 .1 .8 .8])
1 Z: ~) E- P6 B' w
4 y2 Z( _4 l6 t: z& M; ^
# v8 v  U3 P- g4 \! x3 L4 E



4 G4 D9 E( E2 j# z- x$ O8 Y8 d- p/ h4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）

1 P  e% W1 _1 G4 ^, ?: q3 u

sharkN

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