Mapping Toolbox メモ †
地理に関するを描く時に．
例 †
南極を中心点として，2種類の projection 方法で描いてみる．
clear
close all
% data 1
landareas = shaperead('landareas.shp','UseGeoCoords',true);
% data 2
[lon,lat] = meshgrid(-150:5:-90, -75:2.5:-45);
some_data = peaks(25);
% --- Lambert Azimuthal Equal-Area Projection ----
figure
axesm('eqaazim', 'Frame', 'on', 'Grid', 'on', 'Origin',[-90 0], 'MapLatLimit',[-90 -30]);
% % data 1
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
% % data 2
geoshow(lat, lon, some_data, 'DisplayType','texturemap', 'FaceAlpha',0.5);
axis off
% --- Orthographic Projection --------------------
figure
axm = axesm('ortho', 'Frame', 'on', 'Grid', 'on', 'Origin',[-90 0], 'MapLatLimit',[-90 -30]);
% % data 1
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
% % data 2
geoshow(lat, lon, some_data, 'DisplayType','texturemap', 'FaceAlpha',0.5);
axis off
% add labels
axesm(getm(axm),'MeridianLabel','on','ParallelLabel','on','LabelRotation','on','PLabelMeridian','prime')
概要 †
axesm で projection と表示範囲等を決定し， geoshow で表示．
面的なデータを可視化するときの geoshow は，geoshow(lat, lon, Z) のように lat, lon の順になることに注意．
これ以外にもいくつか描き方がある．詳細は MathWorks のページで確認．
軸のラベルを表示するには，MeridianLabel と ParallelLabel に関するプロパティをいじる．
どのようにいじれば良いかは Map Axes Properties をよく読む．
参考 †
全て MathWorks 公式のページ．
