利用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
3、matlab程序：
clear;close all;clc
; U7 K7 x0 l, [* u! x  C7 yfigure
; H$ d$ o7 G; @& Tworldmap china
7 T& [8 I, y- z. H6 D( j/ rmstruct = gcm;
, M3 c! H  |% olatlim = mstruct.maplatlimit;
* M& P+ ~" [6 olonlim = mstruct.maplonlimit;
" d$ i9 m6 p4 U6 V$ n: ]% Wscalefactor = 2;
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
   scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
geoshow('landareas.shp','FaceColor','none','EdgeColor','black')
$ B/ d, d% G! |- p/ W: N3 Y# Mload usgslulegend; colormap(cmap)
! |) \# P6 ~# z  Eset(handlem('suRFace'),'CDataMapping','Direct'); tightmap
caxis([.5 24.5]); hcb = colorbar;
# S% Z5 |  [! F6 O% A# Gset(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
% B/ H7 `0 {1 U7 \: M0 s" ~    'Position',[.7 .1 .02 .8])
set(gca,'Position',[-.05 .1 .8 .8])
; l: o. H* g7 l/ g5 M! R3 B


9 m3 k& A* F! c' _7 T# i
/ E8 S, i2 D% w: A: [6 @

/ k) h4 y9 t. D& a7 [! y4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）
% Y8 x' _3 R0 E4 ]- u$ s& i7 k

sharkN

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