导入地理空间数据

1import geodatasets
2import pandas as pd
3
4datasets = pd.DataFrame(geodatasets.data.flatten()).T
5
6print(datasets[['name', 'geometry_type', 'description']].to_html())

1import geodatasets
2
3natural_earth = geodatasets.get_path('naturalearth.land')
4print(f'Natural Earth location: {natural_earth}\n')
5
6nybb = geodatasets.get_path('ny.bb')
7print(f'New York Borough Boundaries location: {nybb}')

代码示例

1zipfile = zip:///local_path/zippedfile.zip!folder/filename

1import geopandas as gpd
2
3# 读取形状文件
4gdf = gpd.read_file('path/to/shapefile.shp')

1import geopandas as gpd
2import geodatasets
3
4# 打开 Natural Earth 数据集
5n_earth = geodatasets.get_path('naturalearth.land')
6gdf = gpd.read_file(n_earth)
7
8# 应用函数裁剪几何（用于显示目的）
9gdf['geometry_str'] = gdf.geometry.map(lambda x: str(x)[:50])
10
11# 预览数据框
12print(gdf.head().to_html())

代码示例

1gdf.set_geometry('列名')

1import fiona
2
3# 获取支持的驱动程序
4print(fiona.supported_drivers)

1import geopandas as gpd
2
3gdf = gpd.read_file('https://d2ad6b4ur7yvpq.cloudfront.net/naturalearth-3.3.0/ne_110m_admin_1_states_provinces_shp.geojson')
4
5ax = gdf.plot(column='iso_3166_2', figsize=(7, 5))
6ax.set_ylabel('纬度（度）')
7ax.set_xlabel('经度（度）')
8
9ax.figure.savefig('output/states.png')