我正在使用oceansdb模块和Plotly在Python中生成地形图。但是,这张地图显示了对地形的粗略解释,我想知道如何在两点之间平滑地图?
这是我尝试过的代码:
import plotly.graph_objects as go
import oceansdb
import numpy as np
from scipy.interpolate import griddata
Xa = np.linspace(29.005,29.405,200)
Ya = np.linspace(-93.6683,-93.2683,200)
db = oceansdb.ETOPO()
dcont = db['topography'].extract(lat=Xa, lon=Ya)
depth = dcont['height']
fig = go.Figure(data=[go.Surface(z=depth, x=Xa, y=Ya)])
fig.show()
我已经尝试过scipy.interpolate.griddata
来平滑情节,但是情节中没有任何变化。
答案 0 :(得分:1)
我认为插值不是您想要的,因为如果您在嘈杂的点内进行插值,仍然会产生噪点。 scipy.ndimage
库具有一些不错的平滑算法(我展示的是一种简单的高斯平滑)。
import plotly.graph_objects as go
import oceansdb
import numpy as np
import scipy.ndimage
Xa = np.linspace(29.005,29.405,200)
Ya = np.linspace(-93.6683,-93.2683,200)
db = oceansdb.ETOPO()
dcont = db['topography'].extract(lat=Xa, lon=Ya)
depth = dcont['height']
sigma = [4, 4]
print(sigma)
depthSmooth = scipy.ndimage.filters.gaussian_filter(depth, sigma)
fig = go.Figure(data=[go.Surface(z=depthSmooth, x=Xa, y=Ya)])
fig.show()
您可以更改sigma
值来更改进行多少平滑处理(更大的sigma
表示进行平滑处理,列表为[sigmax, sigmay]
)。