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 で表示. 軸のラベルを表示するには,MeridianLabel と ParallelLabel に関するプロパティをいじる. 参考 †全て MathWorks 公式のページ. |