spline_filter1d#
- scipy.ndimage.spline_filter1d(input, order=3, axis=-1, output=<class 'numpy.float64'>, mode='mirror')[源代码]#
沿给定轴计算一维样条滤波器。
沿给定轴的数组行通过样条滤波器进行滤波。样条的阶数必须 >= 2 且 <= 5。
- 参数:
- inputarray_like
输入数组。
- orderint,可选
样条的阶数,默认为 3。
- axisint,可选
应用样条滤波器的轴。默认为最后一个轴。
- outputndarray 或 dtype,可选
放置输出的数组,或返回数组的 dtype。默认为
numpy.float64
。- mode{‘reflect’,‘grid-mirror’,‘constant’,‘grid-constant’,‘nearest’,‘mirror’,‘grid-wrap’,‘wrap’},可选
mode 参数确定输入数组如何在其边界之外扩展。默认值为 ‘mirror’。每个有效值的行为如下(有关更多绘图和详细信息,请参见边界模式)
- ‘reflect’ (d c b a | a b c d | d c b a)
通过反射最后一个像素的边缘来扩展输入。此模式有时也称为半样本对称。
- ‘grid-mirror’
这是 ‘reflect’ 的同义词。
- ‘constant’ (k k k k | a b c d | k k k k)
通过使用由 cval 参数定义的相同常数值填充超出边缘的所有值来扩展输入。在输入的边缘之外不执行插值。
- ‘grid-constant’ (k k k k | a b c d | k k k k)
通过使用由 cval 参数定义的相同常数值填充超出边缘的所有值来扩展输入。插值也发生在输入范围之外的样本上。
- ‘nearest’ (a a a a | a b c d | d d d d)
通过复制最后一个像素来扩展输入。
- ‘mirror’ (d c b | a b c d | c b a)
通过反射最后一个像素的中心来扩展输入。此模式有时也称为全样本对称。
- ‘grid-wrap’ (a b c d | a b c d | a b c d)
通过环绕到相对边缘来扩展输入。
- ‘wrap’ (d b c d | a b c d | b c a b)
通过环绕到相对边缘来扩展输入,但方式使得最后一个点和初始点完全重叠。在这种情况下,在重叠点处选择哪个样本没有明确的定义。
- 返回:
- spline_filter1dndarray
滤波后的输入。
另请参阅
spline_filter
多维样条滤波器。
注释
ndimage 中的所有插值函数都对输入图像进行样条插值。如果使用 order > 1 的 B 样条,则必须首先将输入图像值转换为 B 样条系数,这是通过沿输入的所有轴顺序应用此一维滤波器来实现的。所有需要 B 样条系数的函数都会自动过滤其输入,这种行为可以通过 prefilter 关键字参数来控制。对于接受 mode 参数的函数,只有当它与过滤时使用的 mode 匹配时,结果才是正确的。
对于复数值 input,此函数独立处理实部和虚部。
1.6.0 版本新增:添加了对复数值的支持。
示例
我们可以使用沿给定轴的一维样条来过滤图像
>>> from scipy.ndimage import spline_filter1d >>> import numpy as np >>> import matplotlib.pyplot as plt >>> orig_img = np.eye(20) # create an image >>> orig_img[10, :] = 1.0 >>> sp_filter_axis_0 = spline_filter1d(orig_img, axis=0) >>> sp_filter_axis_1 = spline_filter1d(orig_img, axis=1) >>> f, ax = plt.subplots(1, 3, sharex=True) >>> for ind, data in enumerate([[orig_img, "original image"], ... [sp_filter_axis_0, "spline filter (axis=0)"], ... [sp_filter_axis_1, "spline filter (axis=1)"]]): ... ax[ind].imshow(data[0], cmap='gray_r') ... ax[ind].set_title(data[1]) >>> plt.tight_layout() >>> plt.show()