基于geopandas库的地理信息处理
在这里插入代码片
1、导入所需要的库
import geopandas as gpd
from learntools.core import binder
binder.bind(globals())
from learntools.geospatial.ex1 import *
2、数据的加载
使用geopandas库自带的read_file函数加载地理信息数据
loans_filepath = "../input/geospatial-learn-course-data/kiva_loans/kiva_loans/kiva_loans.shp"
world_loans = gpd.read_file(loans_filepath)
world_loans.head()
加载世界地图信息
world_filepath = gpd.datasets.get_path('naturalearth_lowres')
world = gpd.read_file(world_filepath)
world.head()
loans_filepath 为shapefile文件的地址。
read_file为geopandas库的函数,用于加载shapefile文件。
world_loans.head()检查加载数据的前五行数据。
3、绘制地图
使用plot()函数根据所加载的地理信息数据继续绘制地图。
ax = world.plot(figsize=(20,20), color='whitesmoke', linestyle=':', edgecolor='black')
world_loans.plot(ax=ax, markersize=2)
4、加载菲律宾地区信息
加载菲律宾地理信息
PHL_loans=world_loans.loc[world_loans.country=="Philippines"].copy()
加载菲律宾地图信息
gpd.io.file.fiona.drvsupport.supported_drivers['KML'] = 'rw'
PHL = gpd.read_file("../input/geospatial-learn-course-data/Philippines_AL258.kml", driver='KML')
PHL.head()
5、绘制菲律宾地区地图
ax = PHL.plot(figsize=(12,12), color='whitesmoke', linestyle=':', edgecolor='black')
PHL_loans.plot(ax=ax, markersize=2)
7总结
第一次学习了geopandas库的地理信息加载以及创建地图。使用了read_file()函数来继续加载文件,以及根据所加载的地理文件信息使用plot()函数进行绘制地图并是地理信息体现在地图上。