Mapping Toolbox メモ

地理に関するを描く時に.

南極を中心点として,2種類の projection 方法で描いてみる.

clear
close all
​
% data 1: land areas
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 -15]);
% % data 1
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
% % data 2
geoshow(lat, lon, some_data,'DisplayType', 'texturemap','FaceAlpha',0.5);
tissot;
axis off
% ------------------------------------------------
​
​
% --- Orthographic Projection --------------------
figure
axesm ('ortho', 'Frame', 'on', 'Grid', 'on', 'Origin',[-90 0], 'MapLatLimit',[-90 -15]);
% % data 1
geoshow(landareas,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);
% % data 2
geoshow(lat, lon, some_data,'DisplayType', 'texturemap','FaceAlpha',0.5);
tissot;
axis off
% ------------------------------------------------

projection_eqaazim.png projection_ortho.png

概要

axesm で projection と表示範囲等を決定し, geoshow で表示.
面的なデータを可視化するときの geoshow は,geoshow(lat, lon, Z) のように lat, lon の順になることに注意.
これ以外にもいくつか描き方がある.詳細は MathWorks のページで確認.

参考

全て MathWorks 公式のページ.


Front page   Edit Diff Attach Copy Rename Reload   New List of pages Search Recent changes   Help   RSS of recent changes
Last-modified: 2019-11-05 (Tue) 04:32:07 (6d)