如何平滑Plotly 3d表面图?

时间:2019-10-18 15:06:18

标签: python scipy plotly interpolation

我正在使用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来平滑情节,但是情节中没有任何变化。

1 个答案:

答案 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])。