plotearth.asv

为画地球椭球二维图的程序

figure
axesm hammer
framem; gridm; mlabel; plabel
setm(gca,'GLineWidth',.01,'GLineStyle','-',...
'MLineLocation',15,'MLabelParallel','equator',...
'MLabelLocation',-120:60:180)
% N +
% S -
% W +
% E -

glat=[47 39 42 55 13 43 -29 43 -82 0 -35 34 -49 -37 -9 21 30 64 -1 -18 -44 5 -7 47 34];
glong=[33 -24 -60 97 43 42 -74 28 -113 -68 -27 -15 -103 38 45.51 86 -96 -51 -48 -23 -47 -55 60 30 41];

c=[54 40 38 38 20 28 24 16 64 44 22 14 32 18 12 54 64 34 22 28 26 28 36 36 22];

scatterm(glat,glong,c,'filled')

glat=[47 39 42 55 13 43 -29 43 -82 0 -35 34 -49 -37 -9 21 30 64 -1 -18 -44 5 -7 47 34];
glong=[33 -24 -60 97 43 42 -74 28 -113 -68 -27 -15 -103 38 45.51 86 -96 -51 -48 -23 -47 -55 60 30 41];

c=[54 40 38 38 20 28 24 16 64 44 22 14 32 18 12 54 64 34 22 28 26 28 36 36 22];



%运行结果为：