海岸線を引く,陸地部分に色付けする †GMT の pscoast みたいなことを GMT や M_MAP に頼らず MATLAB でもやりたい. 海岸線データ取得 †GSHHG のデータを使わせてもらう. $ locate gshhs_c.b.gz /usr/local/src/MATLAB/R2019b/toolbox/map/mapdata/gshhs_c.b.gz さらに細かい解像度が必要な場合は NOAA のここ から拾ってくる. GSHHG と Mapping Toolbox を使って境界を描く,境界面を塗る †GSHHG のファイルから都合の良い領域のデータを切り出す. LatLim = [20,50]; LonLim = [120,155]; S = gshhs('./gshhs_i.b', LatLim, LonLim); S = S([S.Level]==1); figure axis off axm = axesm('eqdconic', 'FLineWidth', 1.0, 'MapLonLimit', LonLim, 'MapLatLimit', LatLim); geoshow(S, 'EdgeColor', 'k','FaceColor',[1.0, 1.0, 0.5]) gridm on framem on mlabel on plabel on axesm(getm(axm),'LabelRotation','on','PLabelLocation',10,'MlabelLocation',10,'PLineLocation',5,'MLineLocation',5) 色は EdgeColor と FaceColor の値を変えれば海岸線・陸地の色が変えられる. 参考 †MathWorks のページ. |