首先确保已经安装了geopand as 和seaborn库。然后按照以下步骤操作：

1.导入必要库

import geopand as as gpd import seaborn as sns %mat pl otlib inline

2.读取地理数据

读取shp文件

gdf = gpd.read_file('file.shp')

3.将坐标系转换为经纬度

gdf = gdf.to_ cr s({'init': 'epsg:4326'})

4.绘制 热力图

sns.kde pl ot(gdf.geometry.x, gdf.geometry.y, cmap='Reds', shade=True, bw=.01)

其中，gdf.geometry.x和gdf.geometry.y分别表示地理数据中所有点的经度和纬度信息，cmap表示颜色映射，shade表示是否阴影化，bw表示核密度估计的带宽。可以根据实际需要调整这些参数。

运行以上代码即可绘制核密度图 热力图 。