|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。/ ?1 y7 X* N7 a# x' r7 n0 m, f. A' K
有两种version、两种投影方式。分别需要不同的matlab函数读取。 U6 w9 i& ?/ R/ @- N2 l5 V- l; }" ~
1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
7 E, l$ o" z- S- O* m. B7 f2、解压数据，得到eausgs2_0la.img
8 `2 n( @; s$ e- x' B3、matlab程序：
6 ?& ~" B2 Y; c- g* L* t% _clear;close all;clc
1 ]+ q5 y* z1 _: f$ p2 [7 cfigure
! v- H8 F/ g( {9 ~4 G# C% j* R6 rworldmap china
H- @& r+ D# U( ~2 |mstruct = gcm;
I8 N: ?6 v% Z$ Llatlim = mstruct.maplatlimit;/ g) q7 ^7 \% E6 ]" U: S) D. O
lonlim = mstruct.maplonlimit;
" w$ J2 d* s( o& m: rscalefactor = 2;1 p* D) |2 u% e
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...: r! t5 H% M. i0 X; f, i ?) c$ L9 n
scalefactor, latlim, lonlim);
& a! Q0 T! i0 Ngeoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap'); % \ o( X9 n' Y" f* T; r# ?
geoshow('landareas.shp','FaceColor','none','EdgeColor','black')
E, n! T. o6 l1 _* uload usgslulegend; colormap(cmap)
; a6 x1 g3 k9 C0 Eset(handlem('suRFace'),'CDataMapping','Direct'); tightmap
# B6 i0 N+ W o( D Y. N! T" wcaxis([.5 24.5]); hcb = colorbar;" Q7 M: e& y# R
set(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...6 ? a1 J9 D7 ]- n
'Position',[.7 .1 .02 .8]). z& y+ N- f/ c. D
set(gca,'Position',[-.05 .1 .8 .8]), [* J9 f7 a' }# t* d' |
- L3 |; z0 l* Y, j( `( P7 f6 c5 y! Z6 Z4 _# k: [* R
; b+ K- r4 K% ?: l9 V' u
1 g" t- Q3 h6 D: ]! t3 _( g/ G- R! ]' H$ o+ }0 M
& X# B8 j+ @% ^8 v3 n1 P4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）. h5 v) v: z) R7 }- }
, H- X4 Y6 l. n# U. Y8 _* }
|
|
/1 
发表于 2020-4-16 10:31
收藏
淘帖
支持！
反对！