我已经从网页https://gadm.org/download_country_v3.html下载了一个国家(厄瓜多尔)的形状文件“ .shp”。 首先,我使用功能将.shp文件转换为熊猫数据框,其中“ sf”是读取的.shp文件。
def read_shapefile(sf):
fields = [x[0] for x in sf.fields][1:]
records = sf.records()
shps = [s.points for s in sf.shapes()]
df = pd.DataFrame(columns=fields, data=records)
df = df.assign(coords=shps)
return df
我已经使用下一个函数绘制了地图:
def plot_map(sf, x_lim = None, y_lim = None, figsize = (20,9)):
plt.figure(figsize = figsize)
id=0
for shape in sf.shapeRecords():
x = [i[0] for i in shape.shape.points[:]]
y = [i[1] for i in shape.shape.points[:]]
plt.plot(x, y, 'k')
if (x_lim == None) & (y_lim == None):
x0 = np.mean(x)
y0 = np.mean(y)
plt.text(x0, y0, id, fontsize=10)
id = id+1
if (x_lim != None) & (y_lim != None):
plt.xlim(x_lim)
plt.ylim(y_lim)
问题在于,该图显示了奇怪的线条,如您在图像中看到的那样。
如何删除这些奇怪的行?
预先感谢您的帮助。