|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
利用matlab画GLCC（GLOBAL LAND COVER CHARACTERIZATION）数据。
$ v: m/ Y9 r7 A9 G( B& j% \有两种version、两种投影方式。分别需要不同的matlab函数读取。
% v5 k: x) I# G8 x2 F% y5 I1、选择version和地区，下载数据。以2.0、欧亚地区、Lambert投影数据格式为例，下载USGS Land Use/Land Cover Scheme数据。
/ `! I+ L* O2 p& Q" Y( ?6 k: E6 c2、解压数据，得到eausgs2_0la.img% l4 B8 W; k! ^, O! L' M( U1 }
3、matlab程序：4 k6 a% d) }; ^5 g* c
clear;close all;clc
! h7 |9 e6 h) `& d! O, Bfigure
% t# K9 N6 m K/ hworldmap china
& a; k% H9 q) ?6 Y! e8 T8 w' Smstruct = gcm;6 L3 V `9 w5 F2 s. R8 ?
latlim = mstruct.maplatlimit; C9 U5 \# z2 q0 z( G6 T9 K
lonlim = mstruct.maplonlimit;: D$ _# ^7 o9 D' i- s
scalefactor = 2;
+ e9 w4 p2 U4 `2 R1 Y[latgrat, longrat, Z] = avhrrlambert('a', 'eausgs2_0la.img', ...% y# M: k* r) d4 A' G8 t
scalefactor, latlim, lonlim);1 j. B# L- F% A
geoshow(latgrat, longrat, Z, 'DisplayType', 'texturemap');
8 ?, J1 `$ D# H! Y. ^) mgeoshow('landareas.shp','FaceColor','none','EdgeColor','black')
$ a4 R# }! w2 P. l+ yload usgslulegend; colormap(cmap)
! j9 L0 _/ |: ~# \set(handlem('suRFace'),'CDataMapping','Direct'); tightmap
2 k- Y6 K! R* b, W6 |2 c* \8 Ccaxis([.5 24.5]); hcb = colorbar;' I8 g# A/ e @
set(hcb,'YTick',1:24,'YTickLabel',USGSLandUse,...
- G6 p0 d1 k) _' t 'Position',[.7 .1 .02 .8]); K+ R8 e4 Q! L% B% o O
set(gca,'Position',[-.05 .1 .8 .8])5 v D: ~; R! x: F$ U+ T |; Y
! E$ J% g- @. Y4 }+ {
! E; u2 Q* q" N& }6 N: d
: V2 s& `$ i# w4 }! E K
5 I9 y" W. p4 X( F/ ^9 K" I& j3 ]
. Z1 a0 ^- j. m4 s2 Y$ x2 s" h. R' s+ [$ {7 Z
4、主要用到的matlab函数为avhrrlambert（如果是Goode投影则用avhrrgoode函数）+ V; }$ E0 t9 \% C" k
\+ V* v; T% e, w. c2 e/ J' d
|
|
/1 
发表于 2020-4-16 10:31
收藏
淘帖
支持！
反对！