全国各省的难抵点(最纵深处)分别在什么地方

  • 2020-09-15 16:15:47
  • 作者: 地图故事
  • 来源: 麻辣GIS

难抵点(极)这个词听着比较艰涩,但是其英文名称就容易理解的多,即"pole of inaccessibility",指的是对一个多边形而言(可以是凸多边形也可以是凹多边形),存在某个点,在这个点可以画出位于多边形内的半径最大的圆。

比如,欧亚大陆的难抵点在这个位置:

全国各省难抵点分析

那么,我们或许会有兴趣知道,全国各省的难抵点分别在什么地方?所谓难抵点,其实就是多边形最纵深的地方。

方法也并不算难,首先,获取各省省界的经纬度坐标,形成一个polygon。数据来自gadm网站的公开数据。

然后,获取各省省界的最大经纬度、最小经纬度,确定地理范围,在该地理范围内均匀地选取300*300个点,看哪个点和周围的最小距离是最大的(有点拗口也有点费解,其实大概类似于hausdorff distance),前提是,该点必须要位于polygon内部。

用matplotlib判断点是否在多边形内比较容易,这几句代码就可以做到:

然后,就是将难抵点绘制出来,并且把相应的内接圆绘制出来,显得更为直观。

河南省

我们可以分省看看。比如先看看河南:

河南的最纵深的点是在许昌,而不是省会郑州,不过该点和郑州的距离也不算远就是了。

安徽省

安徽的最纵深点是在巢湖,而不是省会合肥,不过距离也不算远。

四川省

四川省最纵深的点是在甘孜藏族自治州,与省会成都的距离也相当不近,不过可以看到,三个民族自治州占据了全省一半以上的面积,所以这样的结果也并不算意外。

以上的省份抛砖引玉,其他省的情况,各位看官可以自行分析。