Sélectionnez des régions spécifiques du fichier .shp en utilisant Geopandas

Question

J'ai un fichier .shp qui contient les limites des océans et des mers. Mais, au lieu de tout les tracer, je ne m'intéresse qu'à 6. Geopandas crée quelque chose de similaire à un dataframe (appelons-le "df"), tout comme Pandas. Est-il possible de créer une nouvelle base de données ("df1") qui aurait juste ces zones sélectionnées de "df"?

from mpl_toolkits.basemap import Basemap 
import numpy as np 
import matplotlib.pyplot as plt 
import geopandas as gp 

tes = gp.read_file(r'your\path\World_Seas_IHO_v1\World_Seas.shp') 

tes1 = tes[(tes.NAME == "North Pacific Ocean"), 
      (tes.NAME == "South Pacific Ocean")] 

tes1.plot() 

plt.show() 
plt.ion() 

Quand je lance ce "TES1" une erreur:

"objets de la série sont mutable, ainsi ils ne peuvent pas être hachés."

Une idée?

Merci!

Answer 1

(tes.NAME == "North Pacific Ocean"), (tes.NAME == "South Pacific Ocean") est une tuple de série booléenne. Vous ne pouvez pas passer cela comme un indexeur. Vous souhaitez utiliser bitwise ou | pour combiner les séries booléennes, puis utiliser le résultat pour découper la trame de données.

from mpl_toolkits.basemap import Basemap 
import numpy as np 
import matplotlib.pyplot as plt 
import geopandas as gp 

tes = gp.read_file(r'your\path\World_Seas_IHO_v1\World_Seas.shp') 

tes1 = tes[(tes.NAME == "North Pacific Ocean") | 
      (tes.NAME == "South Pacific Ocean")] 

tes1.plot() 

plt.show() 
plt.ion() 

Ou vous pouvez utiliser isin

tes = gp.read_file(r'your\path\World_Seas_IHO_v1\World_Seas.shp') 

tes1 = tes[tes.NAME.isin(['North Pacific Ocean', 'South Pacific Ocean'])] 

tes1.plot() 

plt.show() 
plt.ion()