Mein MATLAB Forum

function earthmap

figure('renderer','opengl')
%axes('units','normalized','position',[0,0,1,1],'color','blue')
%axis off

proj='aitoff';

countries=shaperead('landareas.shp', 'UseGeoCoords', true);
for i=1:numel(countries)
   lon=countries(i).Lon/180*pi;
   lat=countries(i).Lat/180*pi;
   
   subs=[0 find(isnan(lon))]; % indices of nans which seperate independend polygons
   for j=1:numel(subs)-1
      sublat=lat(subs(j)+1:subs(j+1)-1);
      sublon=lon(subs(j)+1:subs(j+1)-1);
     
      lambda=sublon;
      phi=sublat;

      if strcmpi(proj,'equi')
         phi1=acos(2/pi);

         x=lambda.*cos(phi1);
         y=phi;
         z=0*x;

         p=[x(:) y(:) z(:)];
      end

      if strcmpi(proj,'aitoff')

         alpha=acos(cos(phi).*cos(lambda/2));
         
         x=2*cos(phi).*sin(lambda/2)./sinc(alpha);
         y=sin(phi)./sinc(alpha);
         z=0*x;

         p=[x(:) y(:) z(:)];
      end        

      patch('xdata',p(:,1),'ydata',p(:,2),'zdata',p(:,3),...
         'edgecolor','none','facecolor',rand(3,1));
   end
end

end