%%
clear;clc;close all
load geoid
% Create a figure with an Eckert projection.
figure
axesm eckert4; %注意axesm后面的m了吗?,可以使用maps命令查看所有的地图投影的方式,然后选一个
framem; gridm;%显示框架和网格线,注意后面都多了个m,表示map
axis off %关闭外部坐标轴,外部坐标轴不同于map axes
% Display the geoid as a texture map.
geoshow(geoid, geoidrefvec, 'DisplayType', 'texturemap');
% Create a colorbar and title.
hcb = colorbar('southoutside');
set(get(hcb,'Xlabel'),'String','EGM96 Geoid Heights in Meters.')
% Mask out all the land.
geoshow('landareas.shp', 'FaceColor', 'white');
clear;clc;close all
maps %查看当前可用的地图投影方式
worldmap world %绘制世界地图
load coast %导入数据,全球海岸线
patchm(lat,long,'y'); %创建补片图像对象，是底层的图像函数，lat和long为经纬度，y表示黄色填充
%% 设置属性
setm(gca);%查看当前可以设置的所有图形坐标轴(map axes)的属性
setm(gca,'Frame','on');%使框架可见
getm(gca,'Frame');%使用getm可以获取指定的图形坐标轴的属性
setm(gca,'Grid','on');%打开网格
setm(gca,'MLabelLocation',60);%标上经度刻度标签,每隔60度
setm(gca,'MeridianLabel','on');%设置经度刻度标签可见
setm(gca,'PLabelLocation',[-90:30:90])%标上经度刻度标签,[-90:30:90]
setm(gca,'ParallelLabel','on');%设置经度刻度标签可见
setm(gca,'MLabelParallel','south');%将经度刻度标签放在南方,即下部
setm(gca,'Origin',[0,90,0]);%设置地图的中心位置和绕中心点和地心点的轴旋转角度[latitude longitude orientation]
setm(gca,'PLabelMeridian',90);%将纬度标签放置在经度为90度的地方
worldmap('world')
load coast %导入数据,全球海岸线
plotm(lat,long)
figure%将图画在另外的一张图表上，以免叠加
worldmap('world')
scaleruler
ax=worldmap('world')
load coast
geoshow(ax,lat,long,'Displaytype','Polygon','facecolor',[.45 .60 .30])
states=shaperead('worldstatelo.shp','UseGeoCoords',true);
for k=1:numel(states)
states(k).Number=k;
end
facecolors=makesymbolspec('Polygon',{'Number',[1 numel(states)],'facecolor',polcmap(numel(states))});
geoshow(ax,states,'Displaytype','Polygon','symbolspec',facecolors)
