全國各省的難抵點（最縱深處）分別在什麼地方？

難抵點這個詞聽著比較艱澀，但是其英文名稱就容易理解的多，即"pole of inaccessibility"，指的是對乙個多邊形而言（可以是凸多邊形也可以是凹多邊形），存在某個點，在這個點可以畫出位於多邊形內的半徑最大的圓。

比如，歐亞大陸的難抵點在這個位置：

上圖來自維基百科。

那麼，我們或許會有興趣知道，全國各省的難抵點分別在什麼地方？所謂難抵點，其實就是多邊形最縱深的地方。

方法也並不算難，首先，獲取各省省界的經緯度座標，形成乙個polygon。資料來自gadm**的公開資料。

然後，獲取各省省界的最大經緯度、最小經緯度，確定地理範圍，在該地理範圍內均勻地選取300*300個點，看哪個點和周圍的最小距離是最大的（有點拗口也有點費解，其實大概類似於hausdorff distance），前提是，該點必須要位於polygon內部。

用matplotlib判斷點是否在多邊形內比較容易，這幾句**就可以做到：

from matplotlib import path
def pointin********6(p, t):
polygon = path.path(t)
point = [p[0], p[1]]
res = polygon.contains_point(point)
return res

然後，就是將難抵點繪製出來，並且把相應的內接圓繪製出來，顯得更為直觀。

然後，我們可以分省看看。

比如先看看河南：

河南的最縱深的點是在許昌，而不是省會鄭州，不過該點和鄭州的距離也不算遠就是了。

安徽的最縱深點是在巢湖，而不是省會合肥，不過距離也不算遠。

四川省最縱深的點是在甘孜藏族自治州，與省會成都的距離也相當不近，不過可以看到，三個民族自治州佔據了全省一半以上的面積，所以這樣的結果也並不算意外。

其他省的情況，各位看官也可以自行分析。