EDA365电子论坛网
标题:
利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据
[打印本页]
作者:
uperrua
时间:
2020-4-16 10:31
标题:
利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据
利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
, g* |2 H% ?- ]) ~! E
有两种version、两种投影方式。分别需要不同的matlab函数读取。
/ B- M E* K0 o8 y' W% C
1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
# h4 e6 x4 K1 T( R
2、解压数据，得到eausgs2_0la.img
! R v: `7 ~+ C$ z' D Z
3、matlab程序：
. ^6 ]' q2 G# y# H/ U2 A0 Q. P
clear;close all;clc
/ R+ b8 e6 q E2 N
figure
1 `; q, U; ?* o* r% q/ D5 R
worldmap china
- k C' i. l7 Z
mstruct = gcm;
/ g6 W: C u. a" x# D" e
latlim = mstruct.maplatlimit;
5 }) j! [3 @/ {; l
lonlim = mstruct.maplonlimit;
4 ?0 q& A, Z8 M
scalefactor = 2;
8 r2 g* V- j& H" @3 \5 a( d( ^! C
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
* c( I' N' \( |, |$ p
scalefactor, latlim, lonlim);
& G$ }5 v- B; O9 r8 ^5 }+ ?# J
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
0 E% D% K! F& ?
geoshow('landareas.shp','FaceColor','none','EdgeColor','black')
t8 v. B. H' Y: S/ P! T0 M/ J
load usgslulegend; colormap(cmap)
9 z: L& k* o0 D D2 P) d+ `; U, T
set(handlem('surface'),'CDataMapping','Direct'); tightmap
( p7 X; M; f' \, ?0 t
caxis([.5 24.5]); hcb = colorbar;
- m, C6 y6 y+ o0 k" p- m
set(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
8 u2 Z6 u+ W% m1 O! ?
'Position',[.7 .1 .02 .8])
( V+ ]) J3 U8 q; [! x) C& X
set(gca,'Position',[-.05 .1 .8 .8])
: G8 C6 I/ ?& m0 b$ S8 d7 P
}5 j7 y/ [8 c5 ~# z0 j( P) w" D
. V4 k4 X/ O7 W8 Z# W
m3 P, ~! N% A. M7 @0 ~
7.jpg
(24.13 KB, 下载次数: 9)
下载附件
保存到相册
2020-4-16 10:31 上传
4 |" f% K, _" {+ N3 e
7 b) P w* {8 h/ M5 S: _
; `! ]/ n: f: X$ @
4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）
" g# X2 Y% D6 Y: V* w, u6 y
: L# I& I1 a j& Y
作者:
sharkN
时间:
2020-4-16 18:23
利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2
