|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
, F5 Z3 [( Y/ y有两种version、两种投影方式。分别需要不同的matlab函数读取。
( i" ^4 p) s& ^! t& v9 O, g1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
8 ^: u3 ]9 R; t# V f2、解压数据，得到eausgs2_0la.img1 q, @) y% ~! y# X* e
3、matlab程序：
' i: j: @; R- m% m: aclear;close all;clc
$ p1 s5 I- q- ~5 u" [. ^' W/ b' [2 Xfigure+ z5 u8 l5 j, a1 B# s3 `# I8 N& L8 Y% z$ t
worldmap china9 u& T4 \1 ^4 s+ A3 w
mstruct = gcm;% ~7 e2 R/ N/ w$ k) i5 d
latlim = mstruct.maplatlimit;
) n' ?; x. q! a" Xlonlim = mstruct.maplonlimit;! }$ e/ a% N8 |/ o0 \5 D( G: f
scalefactor = 2;1 g* y, n; @& u, Y6 Q" O2 |' W- ?( W
[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...8 q* |3 M2 f6 [1 U3 ` l" A
scalefactor, latlim, lonlim);6 j$ u+ W* ?7 x' ?% V- `
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
( [; z" _& W& `7 M$ F" Rgeoshow('landareas.shp','FaceColor','none','EdgeColor','black')
- r7 U" z' T/ W/ r+ ` L; Dload usgslulegend; colormap(cmap)
2 C# \% w0 u g1 h" lset(handlem('suRFace'),'CDataMapping','Direct'); tightmap1 J' \) E% A2 H# _% }, [
caxis([.5 24.5]); hcb = colorbar;
; q+ l4 a8 X- n( `$ d: qset(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...4 I3 s; C" C$ c. X0 N
'Position',[.7 .1 .02 .8])
: T& j: Q0 y1 K5 }; m7 M) Lset(gca,'Position',[-.05 .1 .8 .8])5 Q3 P4 I" N. }* y$ ]- ^
- W# E0 t% A$ ~: S
' A3 ?: l) [9 B! i" P3 D: A
! w/ q( n- S& G- Q
9 k( B. @! n& h# |3 U% q# C
* e' L" E8 j4 @% J7 h4 P3 M1 U& L
, `: s8 ]6 r/ [& U3 @& l3 ^! P4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）0 T* F3 O2 c/ [6 Y- V: I
' y/ T3 X) F9 t1 y7 V" c2 y |
|
/1 
发表于 2020-4-16 10:31
收藏
淘帖
支持！
反对！