|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
利用matlab画GLCC(GLOBAL LAND COVER CHARACTERIZATION)数据。
4 ]1 R- o, p. x) B有两种version、两种投影方式。分别需要不同的matlab函数读取。3 O0 X% w9 s& |, u: S! Q* O2 R
1、选择version和地区,下载数据。以2.0、欧亚地区、Lambert投影数据格式为例,下载USGS Land Use/Land Cover Scheme数据。1 I& a) u% W) t2 U8 @
2、解压数据,得到eausgs2_0la.img
: @+ r$ j u& ?% M7 j( t3、matlab程序:( G0 y; U1 m, G0 b5 p
clear;close all;clc i/ C4 S7 y: |, y
figure/ s$ P) Q+ n( W9 D, {2 A
worldmap china; s& @9 _; F/ P h4 q+ O1 U
mstruct = gcm;% l: `) N" s( }
latlim = mstruct.maplatlimit;6 l. X; S5 u/ U& c0 c9 j+ A6 l
lonlim = mstruct.maplonlimit;
7 p, |5 J- x9 b% n9 iscalefactor = 2;
! M) d" X! G' b- h* m[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
& n: O5 l+ l. K scalefactor, latlim, lonlim);
7 v0 s: z) U1 H/ Cgeoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap'); " v9 U/ M A$ r
geoshow('landareas.shp','FaceColor','none','EdgeColor','black')/ Q4 a( A9 m; [, J, B6 X6 o7 |
load usgslulegend; colormap(cmap)' x5 u& A. M0 |/ V7 I
set(handlem('suRFace'),'CDataMapping','Direct'); tightmap
9 U, F( g; N0 Vcaxis([.5 24.5]); hcb = colorbar;" r. B- ^: W- C a5 L
set(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
) ~+ a7 a5 ]# c; @& R) f 'Position',[.7 .1 .02 .8])/ D% \, `! n) h- r
set(gca,'Position',[-.05 .1 .8 .8])1 a& y' ~ n. P3 ~6 A
: n, \; E$ g$ ]
6 O& L! M5 ?6 ?0 e4 n1 \" Y) E$ ?2 ~9 W
! m5 A, E* P- N
% \3 a' l2 n# q) n) U4 h9 v+ \+ D# @* w" j% s# f
4、主要用到的matlab函数为avhrrlambert(如果是Goode投影则用avhrrgoode函数)
* \0 q4 N8 i8 o* N3 e% ?. h; @5 r: j4 n9 e3 `
|
|