几个缺的图请移驾博客园：cnblogs.com/Gou-Hailong/p/13822460.html

参考：
气象家园： https://mp.weixin.qq.com/s/JWi7VllWr93wDw6BStk8Dw
M-Map：https://www.eoas.ubc.ca/~rich/map.html
找shp的网站：https://gadm.org/download_country_v3.html
中国科学院资源环境数据平台：http://www.resdc.cn/Default.aspx
matlab利用m_map工具包画中国地图及散点云图：https://www.cnblogs.com/righdflf/p/11484189.html
百度经验出图：https://jingyan.baidu.com/article/870c6fc36fdacfb03ee4be58.html

快应用：
shp文件下载：
世界国家+中国大陆部分：https://download.csdn.net/download/Gou_Hailong/12744519
国家基础地理数据：https://download.csdn.net/download/tessykandy/9112343
代码下载：
求内切圆外接圆：https://download.csdn.net/download/Gou_Hailong/12744427
缩放矩阵或图像：https://download.csdn.net/download/Gou_Hailong/12744441
扣shp：https://download.csdn.net/download/Gou_Hailong/12744543

目前搞到的shape文件有：中国科学院资源环境数据、从CSDN上下载的、奇哥给的。

这个数据精度很高，比如，数据2的P2 长度为925，这个数据P2长度是4461；当然精度高同时画起来也超级慢，我的lj 电脑画个图就可以吃一碗泡面了。
它配套一个xls文件：里面存了省市代码。P2数据存在下面的问题：有的省名字缺失，P3数据有的也缺失，但是其对应的代码查询不到，所以不管了。下面补一下P2的省名字。

修正后的shp文件名为：sheng.shp。

下载链接：https://download.csdn.net/download/tessykandy/9112343

结果：

可以看到地级市只有内部线，想要使用，得和省界线配合使用

奇哥给的和上面那个是一样的，多了世界国家。下面将两个整合到一起。

代码用的是 3 中的功能1，然后将 改成 即可。调用：

常用数据路径汇总：

下面对中国行政单位划分做个总结：

https://www.cnblogs.com/Gou-Hailong/p/13545885.html
详细描述请移驾：https://jingyan.baidu.com/article/59a015e398ed3ff794886589.html
行政区划代码查询：https://xingzhengquhua.51240.com/

下载链接：https://www.eoas.ubc.ca/~rich/map.html

下好工具箱解压后，将文件夹放在一个不常动的地方，我放在了中，然后打开matlab，点击，将路径加进去就行了。测试是否成功的话，就在matlab命令空间敲括号里面可以是1-15，这实际上m_map中附带的一个函数，在文件，看看可不可以画出来图。

User Guide：https://www.eoas.ubc.ca/~rich/mapug.html#p2

可用的投影方式有：

参数	中文描述	图
Stereographic	保角，常用于极区
Orthographic	既不是等面积也不是保角的，而是类似于地球的透视图
Azimuthal Equal-area	朗伯方位等面积投影
Azimuthal Equidistant	等距
Gnomonic	点切投影，地图上所有的直线(不仅仅是那些通过中心的)都是大圆路线
Satellite	地球的透视图，由卫星在特定的高度所看到。指定视点高度而不是为地图指定半径
Albers Equal-Area Conic	等积阿波斯圆锥
Lambert Conformal Conic	兰伯特等角圆锥
Mercator	墨卡托，正形切圆柱
Miller Cylindrical	米勒圆柱
Equidistant Cylindrical	等距圆柱
Cylindrical Equal-Area	等积圆柱
Oblique Mercator	斜轴墨卡托
Transverse Mercator	横轴墨卡托
Sinusoidal	伪圆柱投影
Gall-Peters	纬线和子午线的平行线都呈现为直线，但垂直尺度被扭曲
Hammer-Aitoff	具有弯曲子午线和平行线的等面积投影
Mollweide	外德投影，也称为椭圆或等积投影。在这个投影中，平行线是直的。
Robinson	洛滨孙，养眼
UTM	UTM投影

注：本人自己翻译，有的没见过，可能不准，图片来源于上述链接。

通常，全世界的地图都是墨卡托（Mercator），尽管米勒圆柱投影通常看起来更好，因为它没有强调极地区域。 另一个选择是Hammer-Aitoff或Mollweide（使子午线在两极附近弯曲在一起）。 两者都是等面积的。 将这些投影用于在中点附近没有赤道的地图可能不是一个好主意。 Robinson投影不是等面积或等角投影，而是国家地理杂志的选择（无论如何），并且也出现在IPCC报告中。
如果您要绘制的东西的北/南范围较大，但不是很宽（例如北美洲和南美洲，或者北大西洋和南大西洋），那么正弦或Mollweide投影将看起来非常不错。 另一个选择是“横轴墨卡托”，尽管通常只用于非常大比例尺的地图（即，覆盖很小面积的地图）。
对于一个半球或其他半球内较小的区域（例如，澳大利亚，美国，地中海，北大西洋），您可能会选择圆锥投影。 两种可用的圆锥投影之间的差异非常细微，如果您对投影不太了解，那么使用哪种圆锥投影可能不会有太大的区别。
如果您要画的区域变得更小，则使用哪种投影都无关紧要。 我认为在很多情况下有用的一种投影是斜墨卡托，因为您可以将其沿着长（但窄）的沿海地区对齐。 如果沿着经度/纬度线的地图限制正常，请使用横轴墨卡托或圆锥投影。 如果要以米为单位建立网格，则UTM投影也很有用，因为投影坐标（即“东和北”）以米为单位。
极地地区传统上是使用Stereographic投影来绘制地图的，因为每个人似乎都认为具有“靶心”式的纬线模式看起来不错（经度线在极点汇聚在一起的事实也使它们成为极地科学家几乎不可抗拒的目的地 和探险者）。

全代码：https://blog.csdn.net/Gou_Hailong/article/details/108209764

调用实例：

函数全代码请移驾：https://blog.csdn.net/Gou_Hailong/article/details/108209395

先搞个简单的：河南省的
1.第一步，把河南省的shp提取出来，注意看其范围

可以看到，其大致范围是。

2.从 tif 里面将该大致范围搞出来，并提取shp边界，之后用搞出逻辑矩阵

3.把边界外的数据搞成

4.画图

PS：出逻辑矩阵的时候超级慢，这还是一个简单的边界，目前只能这么搞了，还没有其他方法。
跑全中国的跑十几分钟还跑步不来。我丢，中国大陆矩阵大小，本来以为是三百多万，从开始跑，跑到跑了177万，我觉得明天早上应该就跑完了，有点期待夹杂着莫名的兴奋。到了跑到了270万，我想着在有俩小时估计就跑完了，但是10点之后，我发现还没跑完，跑到了310万。我去，不是308万吗，又仔细数了数零，发现是3080万，当场哭晕，那不得10天吗，就没管它，它自己在那跑，到跑到390万停，结果如下图：

算是把雄鸡鸡冠画出来了吧。。

突然想到一个新方法，中国边界内的最大内接圆一定在中国里面，中国边界的最小外接圆一定在中国外边，所以最后只需判断中间夹的那个小圆环就好了。

这是裁北京市的图，小圆太小，大圆太大，效果还是不太理想，求重心求得精确一点就好了。

改进了求重心的算法，效果如上图（右）还可以。
这个功能的算法，目前来说，我是没招了，全代码链接如下：

https://blog.csdn.net/Gou_Hailong/article/details/108209764

但是，对于中国大陆来说，两个圆还是有点少：

上图虽简单，但上面花了我好几个小时，搞出来了这7个圆的信息

我梦寐以求的小公鸡，快点出来吧，这是用mapshow(shp)画出来的。 还在苦逼等待程序运行中…（真累，明天换个活，先歇歇。）
经历了大约两天时间，终于跑完了，加上那七个圆，算法运行时间大大减少了，之前的话少说也得10天，现在大约两天跑好了，电脑烫的厉害。下面看下结果：

这张图片具有重大的“历史纪念”意义，希望CSDN不要和谐掉。。