Coding report@20231103 --> Add spatial visualization notebook. Still exist bugs when visualizing the spatial frequency of different samples.

[KeranLi]

counts = df.groupby(['Latitud', 'Logitud'])['Latitud'].count()

counts = counts.reset_index(name='counts')

data = df.merge(counts, on=['Latitud', 'Logitud'])

# 载入地形数据
grid = pygmt.datasets.load_earth_relief(resolution='01m', region=[70, 140, 0, 55])

# 创建图形
fig = pygmt.Figure()

# 绘制地形图  
fig.grdimage(grid=grid, projection='M15c', frame='a', cmap='geo')

# 在地图上绘制点 
fig.plot(x=data['Logitud'], y=data['Latitud'], c=data['counts'], cmap='hot', style='c0.5c')

# 添加颜色条标注  
fig.colorbar(frame=["a2000", "x+lElevation", "y+lm"])

fig.colorbar(frame='af+lFrequency')

# 显示结果
fig.show()

Canceled future for execute_request message before replies were done
在当前单元格或上一个单元格中执行代码时 Kernel 崩溃。请查看单元格中的代码，以确定故障的可能原因。有关详细信息，请单击 此处。有关更多详细信息，请查看 Jupyter log。