clear,clc
grs80 = almanac('earth','grs80','km');
domeRadius = 3000;
domeLat = 39;
domeLon = -77;
domeAlt = 0;
ax = axesm('globe','Geoid',grs80);
set(ax,'Position',[0 0 1 1]);
axis equal off
view([-167 28])
set(gcf,'Renderer','opengl')
load topo
geoshow(topo,topolegend,'DisplayType','texturemap')
demcmap(topo)
land = shaperead('landareas','UseGeoCoords',true);
plotm([land.Lat],[land.Lon],'Color','black')
rivers = shaperead('worldrivers','UseGeoCoords',true);
plotm([rivers.Lat],[rivers.Lon],'Color','blue'); |