用matlab绘制中国地图
reference:https://jingyan.baidu.com/article/870c6fc36fdacfb03ee4be58.html
shp: http://muchong.com/html/201304/5748467.html
MATLAB是著名的科学软件，具有绘图、编程、仿真等强大的功能。现在介绍MATLAB绘制中国地图的三种方式，分别是（1）使用m_map工具箱命令绘制中国地图；（2）使用mapshow和geoshow命令直接绘制中国地图；（3）使用worldmap+geoshow命令绘制中国地图。
由于没有找到合适的shp数据包，P/L用同一个shp；老旧的shp中缺少南海部分；
%E:\wml\matlab_workspace\map\chain_sharp.shp close all, clear all, clc, dbstop if error infoL=shapeinfo(\'E:\matlab_workspace\map\chain_sharp.shp\') infoP=shapeinfo(\'E:\matlab_workspace\map\chain_sharp.shp\') ChinaL=shaperead(\'E:\matlab_workspace\map\chain_sharp.shp\'); ChinaP=shaperead(\'E:\matlab_workspace\map\chain_sharp.shp\'); bou2_4lx=[ChinaL(:).X]; bou2_4ly=[ChinaL(:).Y]; bou2_4px=[ChinaP(:).X]; bou2_4py=[ChinaP(:).Y]; figure(1) subplot(1,2,1); m_proj(\'lambert\',\'lon\',[70,140],\'lat\',[0,60]);m_plot(bou2_4lx,bou2_4ly,\'k\'); m_grid; subplot(1,2,2); m_proj(\'robinson\',\'lon\',[70,140],\'lat\',[0,60]);m_plot(bou2_4px,bou2_4py,\'k\'); m_grid; m_coast(\'color\',\'r\',\'linewidth\',1);
左图使用的兰伯特（lambert）投影，右图使用的罗宾逊（robinson）投影和增加了海岸线（coast）。
fnshp_L=\'E:\matlab_workspace\map\chain_sharp.shp\';
fnshp_P=\'E:\matlab_workspace\map\chain_sharp.shp\';
infoL=shapeinfo(fnshp_L);
infoP=shapeinfo(fnshp_P);
readL=shaperead(fnshp_L);
readP=shaperead(fnshp_P);
subplot(2,2,1);mapshow(readL,...
\'DefaultFaceColor\', \'green\', ...
\'DefaultEdgeColor\', \'black\');
title(\'China mapshow PolyLine_1\')
subplot(2,2,2);mapshow(readP,...
\'DefaultFaceColor\', \'green\', ...
\'DefaultEdgeColor\', \'black\');
title(\'China mapshow Polygon_2\')
subplot(2,2,3);geoshow(readL, ...
\'DefaultFaceColor\', \'red\', ...
\'DefaultEdgeColor\', \'black\');
title(\'China geoshow PolyLine-3\')
subplot(2,2,4);geoshow(readP, ...
\'DefaultFaceColor\', \'red\', ...
\'DefaultEdgeColor\', \'black\');
title(\'China geoshow Polygon-4\')
直接使用MATLAB本身自带的mapshow和geoshow命令绘制中国地图。这两个命令可以直接读取.shp格式的中国国界省界数据。
其中图1和图2使用的mapshow命令画得，并且图1使用的.shp数据为PloyLine属性的，图2使用的.shp数据为Ploygon属性的；图3和图4使用的geoshow命令画得，并且图3使用的.shp数据为PloyLine属性的，图4使用的.shp数据为Ploygon属性的。值得注意的是，title(‘ ’)命令中的下划线可以起到下角标的作用，而中划线就无此作用（对比Figure 2的前两幅图和后两幅图的标题就会发现）。
fnshp_L=\'E:\matlab_workspace\map\chain_sharp.shp\';
fnshp_P=\'E:\matlab_workspace\map\chain_sharp.shp\';
% infoL=shapeinfo(fnshp_L);
% infoP=shapeinfo(fnshp_P);
% readL=shaperead(\'E:\matlab_workspace\map\chain_sharp.shp\');
% readP=shaperead(\'E:\matlab_workspace\map\chain_sharp.shp\');
%figure(3)
subplot(2,2,1);
worldmap(\'China\');
mapshow(fnshp_L,...
\'DefaultFaceColor\', \'green\', ...
\'DefaultEdgeColor\', \'black\');
subplot(2,2,2);
worldmap([0,55],[70,140]);
mapshow(fnshp_P,\'FaceColor\',[0.5,1.0,0.5]);
setm(gca,\'MLineLocation\',10);
setm(gca,\'PLineLocation\',10);
setm(gca,\'MLabelLocation\',20);
setm(gca,\'PLabelLocation\',10);
subplot(2,2,3);
worldmap(\'China\');
geoshow(fnshp_L, ...
\'DefaultFaceColor\', \'white\', ...
\'DefaultEdgeColor\', \'black\');
subplot(2,2,4);
worldmap([0,55],[70,140]);
geoshow(fnshp_P,\'FaceColor\',[0.5,1.0,0.5]);
setm(gca,\'MLineLocation\',10);
setm(gca,\'PLineLocation\',10);
setm(gca,\'MLabelLocation\',20);
setm(gca,\'PLabelLocation\',10);
title(\'中国地图\',\'FontSize\',14,\'FontWeight\',\'Bold\')
其中worldmap([0,55],[70,140])为设置显示纬度经度范围；setm(gca,\'MLineLocation\',10)为设置经度间隔；setm(gca,\'PLineLocation\',10)为设置纬度间隔；setm(gca,\'MLabelLocation\',20)为设置经度标签每隔几度；setm(gca,\'PLabelLocation\',10)为设置纬度标签每隔几度。
可以看出worldmap+mapshow命令不起作用（前两张图不显示）；然后worldmap+geoshow命令起作用（后两张图效果不错）。
