Почему мой шейп-файл содержит другой многоугольник, чем предполагалось? ⇐ Python

[Anonymous]

У меня есть данные о рисках, которые я хотел бы раскрасить на карте в соответствии с уровнем риска. Я читаю шейп-файл и файл данных CSV, которые затем объединяю. Это очень хорошо работает с шейп-файлом adm1.
Когда я запускаю тот же сценарий с шейп-файлом adm2, результаты совершенно странные: цветные многоугольники находятся далеко от многоугольника с данными. У меня есть прикрепленная карта с тремя примерами цветных многоугольников и черными точками, где указано расположение данных.
Я буду признателен, если кто-нибудь сможет дать представление о том, что происходит, и если это возможно. как решить эту проблему.
Скрипт Python, использованный для создания графика:

Код: Выделить всё

#!/home/zmumba/anaconda3/bin/python

import pandas as pd
import geopandas as gpd
import matplotlib.pyplot as plt
from matplotlib.colors import ListedColormap
from mpl_toolkits.basemap import Basemap  # optional ?

map_df = gpd.read_file("/home/zmumba/DA/Dekad_Progs/Shapefiles/Lesotho/geoBoundaries-LSO-ADM2-all/geoBoundaries-LSO-ADM2.shx")
risks_df=pd.read_csv("/home/zmumba/DA/Dekad_Progs/Output/H1dRrisks.csv")

map_df["risk"] = map_df.merge(risks_df, left_on="shapeName", right_on="District")["risk"]

colors = {1: "green", 2: "yellow", 3: "orange", 4: "red"}  # or a list
labels = {1: "no risk", 2: "low risk", 3: "medium risk", 4: "high risk"}
catego = map_df["risk"].astype(str).str.cat(map_df["risk"].map(labels), sep="- ")

fig, ax = plt.subplots(figsize=(5, 5))
plt.title(f'Risk of Heavy 24hr Rain: 20-24Nov', y=1.04)

map_df.plot(
column=catego,
categorical=True,
edgecolor="k",
linewidths=0.8,
alpha=0.7,
cmap=ListedColormap([c for r,c in colors.items() if r in map_df["risk"].unique()]),
legend=True,
legend_kwds={
"title": "Risk Level",
"shadow": True,
"loc": "lower right",
"fontsize": 10,
},
ax=ax,
)
ax.set_axis_off()
plt.savefig('H1dRriskmap.png', dpi=300)

plt.show()

Подробнее здесь: https://stackoverflow.com/questions/792 ... n-intended