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

uperrua

利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
有两种version、两种投影方式。分别需要不同的matlab函数读取。
  K( O/ P9 W" k/ v5 r6 I1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
2、解压数据，得到eausgs2_0la.img
! K! ~6 z( y* |& J+ A8 f, v3、matlab程序：
clear;close all;clc
figure
worldmap china
( c: G2 F3 w, m( |3 l' S# Vmstruct = gcm;
latlim = mstruct.maplatlimit;
9 I. o6 E3 B" s( ulonlim = mstruct.maplonlimit;
! L& \/ \$ v' N! P% p, e- c2 o9 Uscalefactor = 2;
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
   scalefactor, latlim, lonlim);
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
% ^9 ~) z8 j. I9 k+ Qgeoshow('landareas.shp','FaceColor','none','EdgeColor','black')
* S) I( [% \6 {4 J- c: r5 lload usgslulegend; colormap(cmap)
& H- A2 D0 G, z' i% g5 {/ h$ bset(handlem('suRFace'),'CDataMapping','Direct'); tightmap
6 W* u  g+ ^0 S$ g) H3 Q/ h+ @caxis([.5 24.5]); hcb = colorbar;
set(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
    'Position',[.7 .1 .02 .8])
, [' D, G. D/ n4 e; oset(gca,'Position',[-.05 .1 .8 .8])


* A. N  T6 L4 _! d


+ m, e) q9 @$ \2 K9 F# J8 c5 ~+ z/ H$ A
' K" G; b) g* Y. T4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）
" a/ D' f) ~; ~( [

sharkN

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