我正在尝试从here中读取一些.gdb文件(文件夹?)。
我使用GeoPandas并执行以下操作:
# file from local path
mallard = gpd.read_file('./bird-species/E00039600_mallard.gdb/')
# geopandas included map, filtered to just this hemisphere
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
westhem = world[(world['continent'] == 'North America') |
(world['continent'] == 'South America')]
# making sure the coordinates line up:
mallard = mallard.to_crs(world.crs)
#establishing figure axes
base = westhem.plot(color='white', edgecolor='black',figsize=(11,11))
# cmap because I'd LIKE the multiple layers to exist
bbw_duck.plot(ax=base, cmap = 'Reds');
输出看起来像这样:
通常,在GeoPandas或Python(Jupyter Notebook)中是否可以查看所有图层?
答案 0 :(得分:4)
是的,GeoPandas支持图层。由于各层的名称非常长,因此我建议使用各层的顺序。
mallard_0 = gpd.read_file('./bird-species/E00039600_mallard.gdb/', layer=0)
mallard_1 = gpd.read_file('./bird-species/E00039600_mallard.gdb/', layer=1)
# geopandas included map, filtered to just this hemisphere
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
westhem = world[(world['continent'] == 'North America') |
(world['continent'] == 'South America')]
# making sure the coordinates line up:
mallard_0 = mallard_0.to_crs(world.crs)
mallard_1 = mallard_1.to_crs(world.crs)
# establishing figure axes
base = westhem.plot(color='white', edgecolor='black', figsize=(11, 11))
# cmap because I'd LIKE the multiple layers to exist
mallard_0.plot(ax=base, color='red', alpha=.5)
mallard_1.plot(ax=base, color='blue', alpha=.5)
编辑:Geopandas正在使用引擎盖下的Fiona进行文件处理,因此,如果要查看图层列表,请使用
import fiona
fiona.listlayers('./bird-species/E00039600_mallard.gdb')
Edit2:遍历所有图层将如下所示:
import fiona
# geopandas included map, filtered to just this hemisphere
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
westhem = world[(world['continent'] == 'North America') |
(world['continent'] == 'South America')]
base = westhem.plot(color='white', edgecolor='black', figsize=(11, 11))
layers = fiona.listlayers('./bird-species/E00039600_mallard.gdb')
for l in layers:
mallard = gpd.read_file('./bird-species/E00039600_mallard.gdb', layer=l)
mallard = mallard.to_crs(world.crs)
mallard.plot(ax=base, color='red', alpha=.1)