plotearth.m

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

%HAMMER Hammer Projection
 
  %This projection is equal-area.  The world is displayed as an ellipse, with 
  %curved parallels and meridians.  It is neither conformal nor equal area.  
  %The only point free of distortion is the center point.  Distortion of shape 
  %and area are moderate throughout.
  
  %This projection was presented by H.H. Ernst von Hammer in 1892.  It is a 
  %modification of the Lambert azimuthal equal area projection.  Inspired by 
  %Aitoff projection, it is also known as the Hammer-Aitoff.  It in turn 
  %inspired the Briesemeister, a modified oblique Hammer projection.  John 
  %Bartholomew's Nordic projection is an oblique Hammer centered on 45 degrees 
  %north and the Greenwich meridian.  The Hammer projection is used in 
  %whole-world maps and astronomical maps in galactic coordinates.
  
  

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



% set circle's latitute and longtitue
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];

% set circle's radium
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];

% plot and filled
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];