|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
5 b" T+ |" k4 j# J& i% J有两种version、两种投影方式。分别需要不同的matlab函数读取。- o$ v& D9 l2 ~2 ~- a" ^
1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。5 t% Q0 g/ g# |; E2 y
2、解压数据，得到eausgs2_0la.img2 i/ `$ c! R. j
3、matlab程序：# `" I9 n. @6 M5 Z: o) d
clear;close all;clc
6 `1 _/ \1 L! t8 p' Ufigure" h9 I* i( G" C( w! {" l; X
worldmap china
. r! f0 ~" c5 \' S8 o* M5 N2 U. Xmstruct = gcm;
, k+ Z5 l: r9 Rlatlim = mstruct.maplatlimit;
0 y; b" i3 d1 W6 p* Zlonlim = mstruct.maplonlimit;
6 `3 W$ I! Y. Y, N6 X% Gscalefactor = 2;
5 b2 Z5 t' n- E. a/ |/ s% p[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...
4 ~& }# g6 {7 b% [4 ~ scalefactor, latlim, lonlim);
' G/ j# A& l. g) b! O9 mgeoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
6 t4 Z5 }: r1 @, H) v( j' cgeoshow('landareas.shp','FaceColor','none','EdgeColor','black')% r- }* {5 `1 P" O, p
load usgslulegend; colormap(cmap)
6 w* B# {* k! a2 hset(handlem('suRFace'),'CDataMapping','Direct'); tightmap
+ i+ S0 y- r f; s2 Wcaxis([.5 24.5]); hcb = colorbar;
5 M2 s: K- h4 [& Jset(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...7 W, e. ?. F5 U1 Q- l" D' E
'Position',[.7 .1 .02 .8])
( A _0 C1 z4 _; |$ b* {9 I9 B! S( uset(gca,'Position',[-.05 .1 .8 .8])
' k7 Y! ~* r' d! V" d+ t4 ~
$ D: K! Q) b& o; l$ b; r" W! r' ?' V. _3 e. O
1 s, t, a$ @& \9 h
' {& Z; H6 ] g2 r1 y
k" X% H9 i, ^+ u3 Y1 f+ S% G" f" |3 c4 `* u5 m7 a/ g
4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）
V; ~9 p; i. s8 w8 J
( k7 U. K/ q& E2 |5 I |
|
/1 
发表于 2020-4-16 10:31
收藏
淘帖
支持！
反对！